The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure

Summary A new simple empirical function has been developed that estimates the free energy of binding for a given protein-ligand complex of known 3D structure. The function takes into account hydrogen bonds, ionic interactions, the lipophilic protein-ligand contact surface and the number of rotatable bonds in the ligand. The dataset for the calibration of the function consists of 45 protein-ligand complexes. The new energy function reproduces the binding constants (ranging from 2.5·10^-2 to 4·10^-14 M, corresponding to binding energies between -9 and -76 kJ/mol) of the dataset with a standard deviation of 7.9 kJ/mol, corresponding to 1.4 orders of magnitude in binding affinity. The individual contributions to protein-ligand binding obtained from the scoring function are: ideal neutral hydrogen bond: -4.7 kJ/mol; ideal ionic interaction: -8.3 kJ/mol; lipophilic contact: -0.17 kJ/mol Ų; one rotatable bond in the ligand: +1.4 kJ/mol. The function also contains a constant contribution (+5.4 kJ/mol) which may be rationalized as loss of translational and rotational entropy. The function can be evaluated very fast and is therefore also suitable for application in a 3D database search or de novo ligand design program such as LUDI.     

Binding of acetylcholine and tetramethylammonium to a cyclophane receptor: anion's contribution to the cation-pi interaction

The interaction of the lipophilic cyclophane 1 with several acetylcholine (ACh) and tetramethylammonium (TMA) salts has been investigated in deuteriochloroform to ascertain the influence of the counterion on the cation-pi interaction. Reliable association constants have been measured for 17 salts of commonly used anions; corresponding binding free energies -DeltaG degrees ranged from over 8 kJ mol(-1) down to the limit of detection. The dramatic dependence of the binding energy on the anion showed that the latter takes part in the process with a passive and adverse contribution, which inhibits cation binding even to complete suppression in unfavorable cases. Thermodynamic parameters for the association of 1 with TMA picrate demonstrate that binding is enthalpic in origin, showing a substantial enthalpy gain (DeltaH degrees = -16.7 kJ mol(-1)) and an adverse entropic contribution (DeltaS degrees = -27.9 J mol(-1) K(-1)). A correlation has been found between the "goodness" of anions as binding partners and the solubility of their salts. Conversion of the anion into a more charge-dispersed species, for example, conversion of chloride into dialkyltrichlorostannate, improves cation binding substantially, indicating that charge dispersion is a main factor determining the influence of the anion on the cation-pi interaction. DFT computational studies show that the variation of the binding free energy of TMA with the counterion is closely accounted for by the electrostatic potential (EP) of the ion pair: guest binding appears to respond to the cation's charge density exposed to the receptor, which is determined by the anion's charge density through a polarization mechanism. A value of -DeltaG degrees = 38.6 kJ mol(-1) has been extrapolated for the free energy of binding of TMA to 1 in chloroform but in the absence of a counterion. The transmission of electrostatic effects from the ion pair to the cation-pi interaction demonstrates that host-guest association is governed by Coulombic attraction, as long as factors (steric, entropic, solvation, etc.) other than pure electrostatics are not prevalent.     

How Can Rotaxanes Be Modified by Varying Functional Groups at the Axle?—A Combined Theoretical and Experimental Analysis of Thermochemistry and Electronic Effects

We present theoretically as well as experimentally determined thermochemical data of the non-covalent interactions in different axle-substituted pseudorotaxanes. The overall interaction energy lies in the region of 35 kJ mol(-1), independent of the substitution pattern at the axle. Because rearrangement energies of 7 and 3 kJ mol(-1) are required for wheel and axle, respectively, the sum of the net interactions of individual non-covalent bonds must exceed 10 kJ mol(-1) to achieve a successful host-guest interaction. The geometrical analysis shows three hydrogen bonds, and the close inspection of the individual dipole moments as well as the individual hydrogen bonds reveals trends according to the different functional groups at the axle. The individual trends for the different hydrogen bonds almost lead to a cancellation of the substitution effects. From solvent-effect considerations it can be predicted that the pseudorotaxane is stable in CHCl(3) and CH(2)Cl(2), whereas it would dethread in water. Comparing experimentally and theoretically calculated Gibbs free enthalpies, we find reasonable agreement if an exchange reaction of one solvent molecule instead of the direct formation reaction is considered.     

Prediction of binding constants of protein ligands: A fast method for the prioritization of hits obtained from de novo design or 3D database search programs

A dataset of 82 protein–ligand complexes of known 3D structure and binding constant K_i was analysed to elucidate the important factors that determine the strength of protein–ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein–ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein–ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 Å in the length and up to 30 °Cs in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 × 10^-2 M to 1 × 10^-14 M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.     

Computational studies of Cu(II)/Met and Cu(I)/Met binding motifs relevant for the chemistry of Alzheimer's disease

A systematic study of the binding motifs of Cu(II) and Cu(I) to a methionine model peptide, namely, N-formylmethioninamide 1, has been carried out by quantum chemical computations. Geometries of the coordination modes obtained at the B3LYP/6-31G(d) level of theory are discussed in the context of copper coordination by the peptide backbone and the S atom of a methionine residue in peptides with special emphasis on Met35 of the amyloid-beta peptide (Abeta) of Alzheimer's disease. The relative binding free energies in the gas phase, DeltaG(g), are calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level of theory, and the solvation affects are included by means of the COSMO model to obtain the relative binding energies in solution, DeltaG(aq). A free energy of binding, DeltaG(aq) = -19.4 kJ mol(-1), relative to aqueous Cu(II) and the free peptide is found for the most stable Cu(II)/Met complex, 12. The most stable Cu(I)/Met complex, 23, is bound by -15.6 kJ mol(-1) relative to the separated species. The reduction potential relative to the standard hydrogen electrode is estimated to be E degrees (12/23) = 0.41 V. On the basis of these results, the participation of Met35 as a low affinity binding site of Cu(II) in Abeta, and its role in the redox chemistry underlying Alzheimer's disease is discussed.     

Computational evidence for methyl-donated hydrogen bonds and hydrogen-bond networking in 1,2-ethanediol-dimethyl sulfoxide

The 1:1 complex of 1,2-ethanediol with dimethyl sulfoxide was studied using density functional theory. A network of three hydrogen bonds holds the complex together, including two in which each methyl group donates to the same hydroxyl oxygen. Four lines of evidence support the existence of methyl-donated hydrogen bonds. The interaction energy is 36 +/- 5 kJ/mol using Becke's three parameter hybrid theory with the 1991 nonlocal correlation functional of Perdew and Wang, and a moderately large basis set (B3PW91/6-311++G**//B3PW91/6-31+G**). To determine the energy of each hydrogen bond, a relaxed potential energy scan was performed in a smaller basis set to break the weaker hydrogen bonds by forced systematic rotation of the methyl groups. Two cross-checking analyses show cooperative effects that cause individual hydrogen bond energies in the network to be nonadditive. When one methyl hydrogen bond is broken, the remaining interactions stabilize the complex by storing an additional 2-3 kJ/mol. With all hydrogen bonds intact, the O[bond]H...O[bond]S hydrogen bond contributes 26 +/- 2 kJ/mol stability, and each weak methyl bond stores 5 +/- 2 kJ/mol.     

Influence of water-protein hydrogen bonding on the stability of Trp-cage miniprotein. A comparison between the TIP3P and TIP4P-Ew water models

We report extensive replica exchange molecular dynamics (REMD) simulations on the folding/unfolding equilibrium of Trp-cage miniprotein using the Amber ff99SB all atom forcefield and TIP3P and TIP4P-Ew explicit water solvent models. REMD simulation-lengths in the 500 ns to the microsecond regime per replica are required to adequately sample the folding/unfolding equilibrium. We observe that this equilibrium is significantly affected by the choice of the water model. Compared with experimental data, simulations using the TIP3P solvent describe the stability of the Trp-cage quite realistically, providing a melting point which is just a few Kelvins above the experimental transition temperature of 317 K. The TIP4P-Ew model shifts the equilibrium towards the unfolded state and lowers the free energy of unfolding by about 3 kJ mol(-1) at 280 K, demonstrating the need to fine-tune the protein-forcefield depending on the chosen water model. We report evidence that the main difference between the two water models is mostly due to the different solvation of polar groups of the peptide. The unfolded state of the Trp-cage is stabilized by an increasing number of hydrogen bonds, destabilizing the α-helical part of the molecule and opening the R-D salt bridge. By reweighting the strength of solvent-peptide hydrogen bonds by adding a hydrogen bond square well potential, we can fully recover the effect of the different water models and estimate the shift in population as due to a difference in hydrogen bond-strength of about 0.4 kJ mol(-1) per hydrogen bond.     

Free energy calculations using flexible-constrained, hard-constrained and non-constrained molecular dynamics simulations.

A comparison of different treatments of bond-stretching interactions in molecular dynamics simulation is presented. Relative free energies from simulations using rigid bonds maintained with the SHAKE algorithm, using partially rigid bonds maintained with a recently introduced flexible constraints algorithm, and using fully flexible bonds are compared in a multi-configurational thermodynamic integration calculation of changing liquid water into liquid methanol. The formula for the free energy change due to a changing flexible constraint in a flexible constraint simulation is derived. To allow for a more direct comparison between these three methods, three different pairs of models for water and methanol were used: a flexible model (simulated without constraints and with flexible constraints), a rigid model (simulated with standard hard constraints), and an alternative flexible model (simulated with flexible constraints and standard hard constraints) in which the ideal or constrained bond lengths correspond to the average bond lengths obtained from a short simulation of the unconstrained flexible model. The particular treatment of the bonds induces differences of up to 2 % in the liquid densities, whereas (excess) free energy differences of up to 5.7 (4.3) kJ mol(-1) are observed. These values are smaller than the differences observed between the three different pairs of methanol/water models: up to 5 % in density and up to 8.5 kJ mol(-1) in (excess) free energy.     

Accurate quantum chemical energies for the interaction of hydrocarbons with oxide surfaces: CH(4)/MgO(001)

We examine the adsorption of CH(4) on the MgO(001) surface by a hybrid approach. It combines MP2 calculations with extrapolation to the complete basis set limit for the adsorption site and the CH(4)-CH(4) pair interactions in the adsorbate layer, with DFT+dispersion calculations under periodic boundary conditions for the whole system. To the total binding energy of 10.7 kJ mol(-1), the DFT+D(ispersion) correction contributes 0.7 kJ mol(-1) only, showing that the Mg(9)O(9) two-layer surface model is an excellent choice and that the interaction between the CH(4) molecules in the adsorbate layer is dominated by pair interactions. Contributions due to relaxation of the atom positions of 0.6 kJ mol(-1) (evaluated at DFT+dispersion) and of higher order correlation effects of 2.0 kJ mol(-1) (evaluated by CCSD(T)) yield a final estimate of 13.3 kJ mol(-1). To this total adsorption energy, the lateral interactions between the CH(4) molecules in the adsorbate layer contribute substantially, 4.1 kJ mol(-1)."Observed" desorption energies of 15.3 and 16.0 kJ mol(-1) have been derived from the observed Arrhenius desorption barriers (12.6 and 13.1 kJ mol(-1)) using thermal enthalpy contributions and a substantial zero-point energy (4.2 kJ mol(-1)) calculated from DFT+D vibrational frequencies. The comparison shows that our final hybrid MP2?:?PBE+D+ΔCCSD(T) estimate has reached chemical accuracy. It misses 2-3 kJ mol(-1) of binding only, which is most likely due to missing higher order correlation effects.PBE+D(ispersion) itself yields an adsorption energy that agrees within 1 kJ mol(-1) with our final hybrid MP2?:?PBE+D+ΔCCSD(T) estimate.     

Prediction of Henry's law constant of benzene derivatives using quantum chemical continuum‐solvation models

Semiempirical (SM2, SM5.4A, MST‐AM1, COSMO‐AM1) and ab inito (HF/PCM‐vdW, MP2//PCM‐vdW, COSMO‐DFT) dielectric continuum‐solvation models as well as the surface‐tension model SM5.0R are analyzed with respect to predicting Henry's law constant at 25°C using a compound set of benzene and 39 benzene derivatives. Both hydrophilic and hydrophobic compounds are covered with a total variation in Henry's law constant of almost eight orders of magnitude corresponding to 44 kJ/mol, and the data set is selected such that there are cases where subtle changes in the molecular structure result in substantial changes of the free energy of solvation. The calculations with SM2, COSMO‐AM1, and COSMO‐DFT include solution‐phase geometry optimization, and the ab initio results refer to polarized basis sets of double‐zeta quality, with two gradient‐corrected functionals (BPW and BLYP) being used for the DFT‐based models. The results show considerable differences in performance between the different continuum‐solvation models, and among the methods yielding solvation free energies the systematic error ranges from −0.9 kJ/mol (SM5.0R) to 12.1 kJ/mol (MP2//PCM‐vdW). In particular, the nonelectrostatic solvation energy contributions of SM2, SM5.4A, MST‐AM1, and PCM‐vdW do not correlate with each other, and with PCM‐vdW omission of the nonelectrostatic component significantly improves the relative trend. The best statistics after scaling through linear regression are achieved with the electrostatic component of MP2//PCM‐vdW (r=0.94) and with COSMO‐DFT (r=0.93). The discussion includes detailed analyses of pecularities associated with certain functional groups, deviations from the expected relationship between dipole moment and solvation energy, and a simple approach to model dispersion interaction and cavitation energy by surface area terms that differentiate between individual atom types. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 17–34, 2000     

最小冰聚体结构单元模型及其对水性质影响的量化分析

水的一些奇特性质主要源于水分子之间存在的氢键, 但在分子尺度上的氢键结构和数据仍是目前研究和争论的焦点. 统计分析了目前文献中普遍采用的水分子和氢键结构数据, 并在此基础上应用AutoCAD图形软件模拟出(H2O)10结构的最小冰聚体结构单元(Minimum Ice Structural Unit, MISU)模型, 以及由MISU聚合而成的冰晶体三维模型. 根据MISU模型, 可以计算得到冰在0 ℃融化为水、水由0 ℃加热至100 ℃、水在100 ℃汽化为水蒸气的三相转化过程中分别需要吸收5.86, 4.40和24.94 kJ•mol－1的能量以断裂16.7%, 12.5%和70.8%的氢键. 如若不考虑氢键的影响, 那么计算得到水的融化热和汽化热分别为0.15和15.73 kJ•mol－1, 与VIA族氢化物H2S, H2Se, H2Te的融化热和汽化热基本呈线性关系. 另外, 由MISU模型计算得到冰在0 ℃融化为水时, 密度由923.17 kg•m－3增至999.89 kg•m－3, 亦与实际测量数据基本一致.     

Effects of hydrogen-bonding and stacking interactions with amino acids on the acidity of uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Conformational behaviour, hydrogen bond competition and intramolecular dynamics in vanillin derivatives: acetovanillone and 6-hydroxy-3-methoxyacetophenone

The conformational landscape of the structural isomers acetovanillone (apocynin, AV) and 6-hydroxy-3-methoxyacetophenone (HMAP) has been investigated in a supersonic jet using Fourier transform microwave spectroscopy. Two conformers have been detected in the jet-cooled expansion for each molecule (s-cis and s-trans in AV; s-trans and a-trans for HMAP), differing in the relative orientation of the acetyl and methoxy groups. Both molecules are stabilized by O-H···O or O-H···O=C hydroxyl intramolecular hydrogen bonds, either constraining the local conformations of the methoxy group in AV, or that of the acetyl group in HMAP. Internal rotation splittings have been observed in both conformers of each molecule, originated by the acetyl group, that yield information on the influence of the intramolecular hydrogen bonds on the methyl torsion. The similar internal rotation barriers in both molecules (6.6 and 7.4 kJ mol(-1) in AV; 7.3 and 7.0 kJ mol(-1) in HMAP) suggest that the acetyl torsion is only slightly affected by intramolecular hydrogen bonding. The absence of torsional tunnellings due to the methoxy group indicates torsional barriers above 10.2 and 8.9 kJ mol(-1) for AV conformers, 10.1 and 10.4 kJ mol(-1) for HMAP. Conformational ratios and relative free energies have been estimated from relative intensity measurements of the spectral lines. Ab initio (MP2) and density functional calculations using the recent M05-2X empirical functional have been used to aid the experimental work in describing the structures, internal rotation barriers and isomerization potentials.     

Structural properties of water: comparison of the SPC, SPCE, TIP4P, and TIP5P models of water

Molecular-dynamics simulations were carried out for the SPC, SPCE, TIP4P, and TIP5P models of water at 298 K. From these results we determine the following quantities: the absolute entropy using the two-particle approximation, the mean lifetime of the hydrogen bond, the mean number of hydrogen bonds per molecule, and the mean energy of the hydrogen bond. From the entropy calculations we find that nearly all contributions to the total entropy originates from the orientation effects. Moreover, we determine the contributions to the total entropy which originate from the first, second, and higher solvation shells. It is interesting that the limits between solvation shells are clearly visible. The first solvation shell (0.22 < r < 0.36 nm) contributes approximately 43 J mol K to the total entropy; the second solvation shell (0.36 < r < 0.60 nm) contributes approximately 12 J mol K, while contributions from the third and other solvation shells are very small, approximately 2 J mol K in summary. This indicates that water molecules are strongly ordered up to 0.55-0.6 nm around the central water molecule, and beyond this limit the ordering diminishes. The results of calculations (entropy and hydrogen bonds) are compared with the experimental data for the choosing of the best water model. We find that the SPC and TIP4P models reproduce the best experimental values, and we recommend these models for computer simulations of the aqueous solution of biomolecules.     

Infrared multiple photon dissociation spectroscopy of potassiated proline

The structure of proline in [proline + K]+ has been investigated in the gas phase using high level DFT and MP2 calculations and infrared photo dissociation spectroscopy with a free electron laser (FELIX). The respective FELIX spectrum of [proline + K]+ matches convincingly the calculated spectra of two structurally closely related and nearly iso-energetic zwitterionic salt bridge (SB) structures. An additional unresolved band at approximately 1725 cm(-1) matching with the characteristic CO stretching mode of charge solvation (CS) structures points toward the presence of a minor population of these conformers of proline in [proline + K]+. However, theory predicts a significant energy gap of 18.9 kJ mol(-1) (B3LYP/6-311++G(2d,2p)) or 15.6 kJ mol(-1) (MP2) between the lowest CS conformer of proline and the clearly favored SB structure.     

Solvation properties of N-substituted cis and trans amides are not identical: significant enthalpy and entropy changes are revealed by the use of variable temperature 1H NMR in aqueous and chloroform solutions and ab initio calculations

The cis/trans conformational equilibrium of N-methyl formamide (NMF) and the sterically hindered tert-butylformamide (TBF) was investigated by the use of variable temperature gradient 1H NMR in aqueous solution and in the low dielectric constant and solvation ability solvent CDCl3 and various levels of first principles calculations. The trans isomer of NMF in aqueous solution is enthalpically favored relative to the cis (deltaH(o) = -5.79 +/- 0.18 kJ mol(-1)) with entropy differences at 298 K (298 x deltaS(o) = -0.23 +/- 0.17 kJ mol(-1)) playing a minor role. The experimental value of the enthalpy difference strongly decreases (deltaH(o) = -1.72 +/- 0.06 kJ mol(-1)), and the contribution of entropy at 298 K (298 x deltaS(o) = -1.87 +/- 0.06 kJ mol(-1)) increases in the case of the sterically hindered tert-butylformamide. The trans isomer of NMF in CDCl3 solution is enthalpically favored relative to the cis (deltaH(o) = -3.71 +/- 0.17 kJ mol(-1)) with entropy differences at 298 K (298 x deltaS(o) = 1.02 +/- 0.19 kJ mol(-1)) playing a minor role. In the sterically hindered tert-butylformamide, the trans isomer is enthalpically disfavored (deltaH(o) = 1.60 +/- 0.09 kJ mol(-1)) but is entropically favored (298 x deltaS(o) = 1.71 +/- 0.10 kJ mol(-1)). The results are compared with literature data of model peptides. It is concluded that, in amide bonds at 298 K and in the absence of strongly stabilizing sequence-specific inter-residue interactions involving side chains, the free energy difference of the cis/trans isomers and both the enthalpy and entropy contributions are strongly dependent on the N-alkyl substitution and the solvent. The significant decreasing enthalpic benefit of the trans isomer in CDCl3 compared to that in H2O, in the case of NMF and TBF, is partially offset by an adverse entropy contribution. This is in agreement with the general phenomenon of enthalpy versus entropy compensation. B3LY/6-311++G** and MP2/6-311++G** quantum chemical calculations confirm the stability orders of isomers and the deltaG decrease in going from water to CHCl3 as solvent. However, the absolute calculated values, especially for TBF, deviate significantly from the experimental values. Consideration of the solvent effects via the PCM approach on NMF x H2O and TBF x H2O supermolecules improves the agreement with the experimental results for TBF isomers, but not for NMF.     

Isomer- and species-selective infrared spectroscopy of jet-cooled 7H- and 9H-2-aminopurine and 2-aminopurine·H2O clusters

The infrared (IR) spectra of the supersonic-jet cooled 9H- and 7H-tautomers of 2-aminopurine (2AP) and of the 9H-2-aminopurine·H(2)O monohydrate clusters have been measured by mass- and species-selective IR-UV double resonance spectroscopy in the 3200-3900 cm(-1) region, covering the N-H and O-H stretching vibrations. The spectra are complemented by density functional (B3LYP and PW91) and by second-order M?ller-Plesset (MP2) calculations of the electronic energies and vibrational frequenciesof the respective 2AP tautomers and clusters. The 9H- and 7H-2-aminopurine tautomers were definitively identified by the shifts of their NH and NH(2) symmetric and asymmetric stretching frequencies and by comparison to the B3LYP/TZVP calculated IR spectra. The H-bond topologies of the two previously observed 9H-2-aminopurine·H(2)O isomers (Sinha. R. K.; et al. J. Phys. Chem. A2011, 115, 6208) are definitively identified as the "sugar-edge" isomer A and the "trans-amino-bound" isomer B by comparing their IR spectra to the calculated frequencies and IR intensities of the cluster isomers A, B, C, and D, as well as to the IR spectrum of 9H-2AP. The sugar-edge isomer A involves N9-H···OH(2) and HOH···N3 hydrogen bonds and is predicted to be the most stable form. The amino-bound isomer B involves NH(2)···OH(2) and HOH···N1 hydrogen bonds and is calculated to lie 2.5 kJ/mol above isomer A. The H-bond topology of the "cis-amino-bound" isomer C is symmetrically related to isomer B, with a hydrogen bond to the N3 of the pyrimidine group. However, it is calculated to lie 7 kJ/mol above isomer A and indeed is not observed in the supersonic jet. Isomer D involves a single H-bond to the N7 position, is predicted to be 14 kJ/mol above A and is therefore not observed.     

Binding affinities for models of biologically available potential Cu(II) ligands relevant to Alzheimer's disease: an ab initio study

A systematic study of the binding affinities of the model biological ligands X: = (CH3)2S, CH3S-, CH3NH2, 4-CH3-imidazole (MeImid), C6H5O-, and CH3CO2- to (NH3)i(H2O)3-iCu(II)-H2O (i = 3, 2, 1, 0) complexes has been carried out using quantum chemical calculations. Geometries have been obtained at the B3LYP/ 6-31G(d) level of theory, and binding energies, Delta, relative to H2O as a ligand, have been calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level. Solvation effects have been included using the COSMO model, and the relative binding free energies in aqueous solution (Delta) have been determined at pH 7 for processes that are pH dependent. CH3S- (Delta = -16.0 to -53.5 kJ mol(-1)) and MeImid (Delta = -18.5 to -35.2 kJ mol(-1)) give the largest binding affinities for Cu(II). PhO- and (CH3)2S are poor ligands for Cu(II), Delta = 20.6 to -9.7 and 19.8 to -3.7 kJ mol(-1), respectively. The binding affinities for CH3NH2 range from -0.8 to -15.0 kJ mol(-1). CH3CO2- has Cu(II) binding affinities in the ranges Delta = -13.5 to -32.4 kJ mol(-1) if an adjacent OH bond is available for hydrogen bonding and Delta = 10.1 to -4.6 kJ mol(-1) if this interaction is not present. In the context of copper coordination by the Abeta peptide of Alzheimer's disease, the binding affinities suggest preferential binding of Cu(II) to the three histidine residues plus a lysine or the N-terminus. For a 3N1O Cu(II) ligand arrangement, it is more probable that the oxygen ligand comes from an aspartate/glutamate residue side chain than from the tyrosine at position 10. Methionine appears unlikely to be a Cu(II) ligand in Abeta.