Car-Parrinello MD simulations for the Na+ -phenylalanine complex in aqueous solution

Car-Parrinello molecular dynamics (CPMD) calculations are presented for a Na (+)(Phe) complex in aqueous solution and for various stable Na (+)(Phe) complexes and Na (+)(H 2O) n clusters in the gas phase (with up to six water molecules). The CPMD results are compared to available experimental and ab initio reference data, to DFT results obtained with various combinations of density functionals and basis sets, and to previous classical mechanics MD simulations. The agreement with the reference data in the gas phase validates the CPMD method, showing that it is a valid approach for studying these systems and that it describes correctly the competing Na (+)-Phe and Na (+)-H 2O interactions. Analysis of MD trajectories reveals that the Na (+)(Phe) complex in aqueous solution maintains a stable configuration in which the Na (+) cation hovers above the phenyl ring, at an average distance of 3.85 A from the ring center, while remaining strongly bound to one of the carboxylic oxygens of Phe. Constrained MD simulations indicate that the free energy barrier opposing dissociation of the complex exceeds 5.5 kcal/mol. We thus confirm that "cation- pi" interactions between alcali cations and the pi ring, combined with other kinds of interactions, may allow aromatic amino acids to overcome the competition with water in binding a cation.     

Cation-pi interactions with a model for the side chain of tryptophan: structures and absolute binding energies of alkali metal cation-indole complexes

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.     

Development and validation of an integrated computational approach for the study of ionic species in solution by means of effective two-body potentials. The case of Zn2+, Ni2+, and Co2+ in aqueous solutions

In this paper we have developed an effective computational procedure for the structural and dynamical investigation of ions in aqueous solutions. Quantum mechanical potential energy surfaces for the interaction of a transition metal ion with a water molecule have been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective ion-water interactions have been fitted by suitable analytical potentials, and have been utilized in molecular dynamics (MD) simulations to obtain structural and dynamical properties of the ionic aqueous solutions. This procedure has been successfully applied to the Co2+-H2O open-shell system and, for the first time, Co-oxygen and Co-hydrogen pair potential functions have been determined and employed in MD simulations. The reliability of the whole procedure has been assessed by applying it also to the Zn2+ and Ni2+ aqueous solutions, and the structural and dynamical properties of the three systems have been calculated by means of MD simulations and have been found to be in very good agreement with experimental results. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to extended X-ray absorption fine structure (EXAFS) experiments.     

A valence bond model for aqueous Cu(II) and Zn(II) ions in the AMOEBA polarizable force field

A general molecular mechanics (MM) model for treating aqueous Cu²⁺ and Zn²⁺ ions was developed based on valence bond (VB) theory and incorporated into the atomic multipole optimized energetics for biomolecular applications (AMOEBA) polarizable force field. Parameters were obtained by fitting MM energies to that computed by ab initio methods for gas‐phase tetra‐ and hexa‐aqua metal complexes. Molecular dynamics (MD) simulations using the proposed AMOEBA‐VB model were performed for each transition metal ion in aqueous solution, and solvent coordination was evaluated. Results show that the AMOEBA‐VB model generates the correct square‐planar geometry for gas‐phase tetra‐aqua Cu²⁺ complex and improves the accuracy of MM model energetics for a number of ligation geometries when compared to quantum mechanical (QM) computations. On the other hand, both AMOEBA and AMOEBA‐VB generate results for Zn²⁺–water complexes in good agreement with QM calculations. Analyses of the MD trajectories revealed a six‐coordination first solvation shell for both Cu²⁺ and Zn²⁺ ions in aqueous solution, with ligation geometries falling in the range reported by previous studies. © 2012 Wiley Periodicals, Inc.     

The role of cation-pi interactions in biomolecular association. Design of peptides favoring interactions between cationic and aromatic amino acid side chains.

Cation-pi interactions between amino acid side chains are increasingly being recognized as important structural and functional features of proteins and other biomolecules. Although these interactions have been found in static protein structures, they have not yet been detected in dynamic biomolecular systems. We determined, by (1)H NMR spectroscopic titrations, the energies of cation-pi interactions of the amino acid derivative AcLysOMe (1) with AcPheOEt (2) and with AcTyrOEt (3) in aqueous and three organic solvents. The interaction energy is substantial; it ranges from -2.1 to -3.4 kcal/mol and depends only slightly on the dielectric constant of the solvent. To assess the effects of auxiliary interactions and structural preorganization on formation of cation-pi interactions, we studied these interactions in the association of pentapeptides. Upon binding of the positively-charged peptide AcLysLysLysLysLysNH(2) (5) to the negatively-charged partner AcAspAspXAspAspNH(2) (6), in which X is Leu (6a), Tyr (6b), and Phe (6c), multiple interactions occur. Association of the two pentapeptides is dynamic. Free peptides and their complex are in fast exchange on the NMR time-scale, and 2D (1)H ROESY spectra of the complex of the two pentapeptides do not show intermolecular ROESY peaks. Perturbations of the chemical shifts indicated that the aromatic groups in peptides 6b and 6c were affected by the association with 5. The association constants K(A) for 5 with 6a and with 6b are nearly equal, (4.0 +/- 0.7) x 10(3) and (5.0 +/- 1.0) x 10(3) M(-)(1), respectively, while K(A) for 5 with 6c is larger, (8.3 +/- 1.3) x 10(3) M(-)(1). Molecular-dynamics (MD) simulations of the pentapeptide pairs confirmed that their association is dynamic and showed that cation-pi contacts between the two peptides are stereochemically possible. A transient complex between 5 and 6 with a prominent cation-pi interaction, obtained from MD simulations, was used as a template to design cyclic peptides C(X) featuring persistent cation-pi interactions. The cyclic peptide C(X) had a sequence in which X is Tyr, Phe, and Leu. The first two peptides do, but the third does not, contain the aromatic residue capable of interacting with a cationic Lys residue. This covalent construct offered conformational stability over the noncovalent complexes and allowed thorough studies by 2D NMR spectroscopy. Multiple conformations of the cyclic peptides C(Tyr) and C(Phe) are in slow exchange on the NMR time-scale. In one of these conformations, cation-pi interaction between Lys3 and Tyr9/Phe9 is clearly evident. Multiple NOEs between the side chains of residues 3 and 9 are observed; chemical-shift changes are consistent with the placement of the side chain of Lys3 over the aromatic ring. In contrast, the cyclic peptide C(Leu) showed no evidence for close approach of the side chains of Lys3 and Leu9. The cation-pi interaction persists in both DMSO and aqueous solvents. When the disulfide bond in the cyclic peptide C(Phe) was removed, the cation-pi interaction in the acyclic peptide AC(Phe) remained. To test the reliability of the pK(a) criterion for the existence of cation-pi interactions, we determined residue-specific pK(a) values of all four Lys side chains in all three cyclic peptides C(X). While NOE cross-peaks and perturbations of the chemical shifts clearly show the existence of the cation-pi interaction, pK(a) values of Lys3 in C(Tyr) and in C(Phe) differ only marginally from those values of other lysines in these dynamic peptides. Our experimental results with dynamic peptide systems highlight the role of cation-pi interactions in both intermolecular recognition at the protein-protein interface and intramolecular processes such as protein folding.     

Molecular dynamics simulations of helix-forming, amine-functionalized m-poly(phenyleneethynylene)s

We present here the results of all-atom and united-atom molecular dynamics (MD) simulations that were used to examine the folding behavior of an amine-functionalized m-poly(phenyleneethynylene) (m-PPE) oligomer in aqueous environment. The parallelized GROMACS MD simulation code and OPLS force field were used for multiple MD simulations of m-PPE oligomers containing 24 phenyl rings in extended, coiled and helix conformations separately in water to determine the minimum energy conformation of the oligomer in aqueous solvent and what interactions are most important in determining this structure. Simulation results showed that the helix is the preferred minimum energy conformation of a single oligomer in water and that Lennard-Jones interactions are the dominant forces for the stabilization of the helix. In addition, these solvophobic interactions are strong enough to maintain the helix conformation at temperatures up to 523 K.     

Electronic Structure Insights into the Solvation of Magnesium Ions with Cyclic and Acyclic Carbonates

A computational framework to rank the solvation behavior of Mg²⁺ in carbonates by using molecular dynamics simulations and density functional theory is reported. Based on the binding energies and enthalpies of solvation calculated at the M06‐2X/6‐311++G(d,p) level of theory and the free energies of solvation from ABF‐MD simulations, we find that ethylene carbonate (EC) and the ethylene carbonate:propylene carbonate (EC:PC) binary mixture are the best carbonate solvents for interacting with Mg²⁺. Natural bond orbital and quantum theory of atoms in molecules analyses support the thermochemistry calculations with the highest values of charge transfer, perturbative stabilization energies, electron densities, and Wiberg bond indices being observed in the Mg²⁺(EC) and Mg²⁺(EC:PC) complexes. The plots of the noncovalent interactions indicate that those responsible for the formation of Mg²⁺ carbonate complexes are strong‐to‐weak attractive interactions, depending on the regions that are interacting. Finally, density of state calculations indicate that the interactions between Mg²⁺ and the carbonate solvents affects the HOMO and LUMO states of all carbonate solvents and moves them to more negative energy values.     

Halide anion recognition in water by a hexaprotonated octaaza-cryptand: a molecular dynamics investigation

Based on molecular dynamics simulations, we describe the F- versus Cl- complexation by an hexaprotonated cryptand L6+ in aqueous solution, in order to elucidate their structures, solvation properties and the status of external halide counterions. In water, F- and Cl- simulated inclusive complexes adopt a structure somewhat different from the solid state structure of the F- complex: The anion binding involves two diammonium bridges only, and the accompanying counterions are dissociated from the +5 charged complex. A remarkable result is obtained for the dissociated L6+,3F- ,3Cl- system, where spontaneous complexation of F- (the anion which forms the most stable complex with L6+) takes place during the dynamics. The resulting complex is of facial type; this suggests that the equilibrium involves multiple binding modes and structures in aqueous solution. The question of F-/Cl- binding selectivity is investigated by free energy perturbations simulations which nicely reproduce the spectacular preference for F- over Cl-. Two different methodologies used for the treatment of electrostatics (standard versus Ewald calculations) yield similar conclusions.     

Calculations of the free energy of interaction of the c-Fos-c-Jun coiled coil: effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein-protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson-Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos-c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos-c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.     

Molecular dynamics study of the inclusion of cholesterol into cyclodextrins

The interaction of three cyclodextrins (CDs), viz. beta-CD, heptakis (2,6-di-O-methyl)-beta-CD (DM-beta-CD), and 2-hydroxypropyl-beta-CD (HP-beta-CD), with cholesterol was investigated using molecular dynamics (MD) simulations. The free energy along the reaction pathway delineating the inclusion of cholesterol into each CD was computed using the adaptive biasing force method. The association constant and the corresponding association free energy were derived by integrating the potential of mean force (PMF) over a representative ordering parameter. The results show that the free energy profiles possess two local minima corresponding to roughly equally probable binding modes. Among the three CDs, DM-beta-CD exhibits the highest propensity to associate with cholesterol. Ranking for binding cholesterol, viz. DM-beta-CD > HP-beta-CD > beta-CD, agrees nicely with experiment. Partitioning of the PMF into free energy components illuminates that entering of cholesterol into the CD cavity is driven mainly by electrostatic interactions, whereas deeper inclusion results from van der Waals forces and solvation effects. Additional MD simulations were performed to investigate the structural stability of the host-guest complexes near the free energy minima. The present results demonstrate that association of cholesterol and CDs follows two possible binding modes. Although the latter are thermodynamically favorable for all CDs, one of the two inclusion complexes appears to be preferred kinetically in the case of DM-beta-CD.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附. 设置不同的水层厚度, 观察固液界面和气液界面吸附的差异. 模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面, 受碳链和固体表面之间相互作用的影响形成表面活性剂分子层, 并依据吸附量的大小形成不同的聚集结构; 在水层足够厚的情况下, 由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构; 计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关, 解离能垒远大于结合能垒, 引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中; 无论气液还是固液界面, 极性头均伸向水相, 与水分子形成不同类型的氢键. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观结构信息.     

The effects of charge transfer on the aqueous solvation of ions

Ab initio-based charge partitioning of ionic systems results in ions with non-integer charges. This charge-transfer (CT) effect alters both short- and long-range interactions. Until recently, the effects of CT have been mostly neglected in molecular dynamics (MD) simulations. The method presented in this paper for including charge transfer between ions and water is consistent with ab initio charge partitioning and does not add significant time to the simulation. The ions of sodium, potassium, and chloride are parameterized to reproduce dimer properties and aqueous structures. The average charges of the ions from MD simulations (0.900, 0.919, and -0.775 for Na(+), K(+), and Cl(-), respectively) are consistent with quantum calculations. The hydration free energies calculated for these ions are in agreement with experimental estimates, which shows that the interactions are described accurately. The ions also have diffusion constants in good agreement with experiment. Inclusion of CT results in interesting properties for the waters in the first solvation shell of the ions. For all ions studied, the first shell waters acquire a partial negative charge, due to the difference between water-water and water-ion charge-transfer amounts. CT also reduces asymmetry in the solvation shell of the chloride anion, which could have important consequences for the behavior of chloride near the air-water interface.     

Cation-pi interactions and the gas-phase thermochemistry of the Na(+)/phenylalanine complex

The complex of Na(+) with phenylalanine (Phe) is a prototype for the participation of cation-pi interactions in metal-ion binding to biological molecules. A recent comparison of this complex with the Na(+)/alanine (Na(+)/Ala) counterpart suggested only a small contribution of the phenyl ring interaction to binding, casting doubt on the extent of the cation-pi effect. The present work reexamines this thermochemistry using ligand-exchange equilibrium measurements in the Fourier transform ion cyclotron resonance (FT-ICR) ion trapping mass spectrometer. An increment of 7 +/- 2 kcal mol(-1) was found in the Ala/Phe comparison of binding enthalpies, confirming the importance of cation-pi binding enhancement in the Phe case. Absolute Na(+) binding enthalpies of 38 +/- 2 and 45 +/- 2 kcal mol(-1) were assigned for Ala and Phe, respectively, using pyridine as the thermochemical reference ligand. All of these results were supported by quantum calculations using both density functional and Hartree-Fock/MP2 methods, improved in several respects over previous calculations. Alanine methyl ester (AlaMe) was also observed, and found to have an Na(+) ion affinity larger by 2.3 kcal mol(-1) than Ala. New, lower energy conformations of neutral Phe were discovered in the computations.     

The performance of ANI-ML potentials for ligand-n(H2O) interaction energies and estimation of hydration free energies from end-point MD simulations

Here, we investigate the performance of “Accurate NeurAl networK engINe for Molecular Energies” (ANI), trained on small organic compounds, on bulk systems including non-covalent interactions and applicability to estimate solvation (hydration) free energies using the interaction between the ligand and explicit solvent (water) from single-step MD simulations. The method is adopted from ANI using the Atomic Simulation Environment (ASE) and predicts the non-covalent interaction energies at the accuracy of wb97x/6-31G(d) level by a simple linear scaling for the conformations sampled by molecular dynamics (MD) simulations of ligand-n(H₂O) systems. For the first time, we test ANI potentials' abilities to reproduce solvation free energies using linear interaction energy (LIE) formulism by modifying the original LIE equation. Our results on ~250 different complexes show that the method can be accurate and have a correlation of R² = 0.88–0.89 (MAE <1.0 kcal/mol) to the experimental solvation free energies, outperforming current end-state methods. Moreover, it is competitive to other conventional free energy methods such as FEP and BAR with 15-20 × fold reduced computational cost.     

水化Na-蒙脱石和Na/Mg-蒙脱石的分子动力学模拟

利用分子动力学方法(MD)研究了Na-蒙脱石和Na/Mg-蒙脱石层间的补偿阳离子和水分子的结构及扩散性质. 模拟结果表明, 在一定水含量范围内Na-蒙脱石和Na/Mg-蒙脱石表现出不同的膨胀形式, 特别是层间水分子数目在48～72之间时, Na/Mg-蒙脱石的层间距比Na-蒙脱石有较为明显的增大. Na/Mg-蒙脱石两层水化物的层间水分子与Mg2+形成了明显的两层水合壳; 而与Na+只形成了一层平面的水合壳. 在Na/Mg-蒙脱石中, Na+和 Mg2+的扩散方式不同, Na+的扩散范围相对更广, 自扩散系数更大. Na/Mg-蒙脱石比相同水含量下的Na-蒙脱石层间水的自扩散系数小. 由于Mg2+和Na+对层间结构的强烈影响, 从而使有少量Mg2+取代Na+的Na/Mg-蒙脱石与Na-蒙脱石表现出不同的膨胀性质和层间物质的扩散性质.     

DFT:B3LYP ab initio molecular dynamics study of the Zundel and Eigen proton complexes, H5O2+ and H9O4+, in the triplet state in gas phase and solution

DFT:B3LYP ab initio molecular dynamics (MD) approach is used to elucidate the properties of the Zundel and Eigen, H5O2+ and H9O4+, proton complexes in the triplet state. The simulation considers the complexes in the gas phase (isolated complexes) and inside the clusters composed of 32, 64, and 128 water molecules, mimicking the behavior of aqueous solutions. MD simulations reveal three distinct periods. For the complex in solutions, the periods are smoothed out. The H5O2+ and H9O4+ complexes in the triplet state undergo structural rearrangements, which eventually result in hydrogen elimination. For the H5O2+, the hydrogen is eliminated from the center of the water cluster, whereas for the H9O4+ it is removed from a near-surface water molecule. The rate of hydrogen elimination decreases with increasing the number of water molecules surrounding the complex.     

Perturbation of local solvent structure by a small dication: a theoretical study on structural, vibrational, and reactive properties of beryllium ion in water

A molecular based understanding of beryllium chemistry in the context of biomolecules is necessary for gaining progress in prevention and treatment of chronic beryllium disease. One aspect that has hindered the theoretical progress has been the lack of a simple classical two-body potential for the aqueous beryllium ion (Be2+) to be used with biomolecular simulations. We provide new parameters for Be2+ that capture the structural and reactive properties of this small dication. Using classical molecular dynamics simulations, we show that these parameters reproduce the correct radial distribution function and coordination numbers for this cation in explicit aqueous solution when compared to published diffraction and NMR measurements. The geometrical parameters obtained using classical simulations are also in agreement with ab initio calculations. We successfully predict the vibrational modes of the tetra aqua Be2+ dication from ab initio calculations on solvated structures obtained from the simulations. The calculated vibrational modes show better agreement with experiments compared to any published work. This new potential also produces a well-established hydrogen bonding between the first and second solvation shells. More importantly, when the molecular dynamics (MD) and ab initio results are interpreted in concert, the dynamics and nature of interactions between the first and second shells capture the pivotal role they play on the reactivity of aqua-Be complexes.     

Insights into uranyl chemistry from molecular dynamics simulations

First-principles and purely classical molecular dynamics (MD) simulations for complexes of the uranyl ion (UO(2)(2+)) are reviewed. Validation of Car-Parrinello MD simulations for small uranyl complexes in aqueous solution is discussed. Special attention is called to the mechanism of ligand-exchange reactions at the uranyl centre and to effects of solvation and hydration on coordination and structural properties. Large-scale classical MD simulations are surveyed in the context of liquid-liquid extraction, with uranyl complexes ranging from simple hydrates to calixarenes, and nonaqueous phases from simple organic solvents and supercritical CO(2) to ionic liquids.     

Decomposition of the solvation free energies of deoxyribonucleoside triphosphates using the free energy perturbation method

Free energy perturbation (FEP) calculations using the Amber 95 force field and the TIP3P water model were carried out to evaluate the solvation free energy of deoxyribonucleoside triphosphates in aqueous solution. Solvation free energies of -307.5, -311.5, -314.1, and -317.0 kcal/mol were calculated for the (Mg x dTTP)2-, (Mg x dATP)2-, (Mg x dCTP)2-, and (Mg x dGTP)2- complexes, respectively. Structural origins of the relative solvation free energies of deoxyribonucleoside phosphates were examined by calculating the contribution of the interaction of the base moiety with its surroundings. We showed that for each nucleobase the magnitude of this contribution is unaffected by substituting the 5'-OH group of the corresponding nucleoside with the charged mono- or triphosphate groups. This free energy contribution was further decomposed into the sum of free energies originating from the interactions of the base with itself, its substituent, water, and Na+ ions. Although the sum of these components was nearly constant over a wide range of solutes the individual free energy constituents varied significantly. Furthermore, this decomposition showed a high degree of additivity. Computational conditions necessary for obtaining additive free energy decomposition for the systems studied here within the framework of the FEP method included the use of a single mutation pathway and a subdivision of the FEP protocol into 51 or more windows.