Energetics of protein backbone hydrogen bonds and their local electrostatic environment

MD simulation study of several peptides including a polyalanine,a helix(pdb:2I9M),and a leucine zipper were carried out to investigate hydrogen bond energetics using dynamic polarized protein-specific charge(DPPC)to account for the polarization effect in protein dynamics.Results show that the backbone hydrogen-bond strength is generally correlated with its specific local electrostatic environment,measured by the number of water molecules near the hydrogen bond in the first solvation shell.The correlation coefficient is found to be 0.89,0.78,and 0.80,respectively,for polyalanine,2I9M protein,and leucine zipper.In the polyalanine,the energies of the backbone hydrogen bonds are very similar to each other due to their similar local electrostatic environment.The current study helps demonstrate and support the understanding that hydrogen bonds are stronger in a hydrophobic surrounding than in a hydrophilic one.For comparison,the result from simulation using standard force field shows a much weaker correlation between hydrogen bond energy and local electrostatic environment due to the lack of polarization effect in the force field.     

Simulant molecules with trivalent or pentavalent phosphorus atoms: bond dissociation energies and other thermodynamic and structural properties from quantum chemical models

The CBS-QB3 and G4 thermochemical models have been used to generate energetic, structural, and spectroscopic data on a set of molecules with trivalent or pentavalent phosphorus atoms that can serve as simulants of chemical warfare agents. Based on structural data, the conformational stabilities of these molecules are explained in terms of the anomeric interaction within the OPOC and OPSC fragments. For those cases where experimental data are available, comparisons have been made between calculated and previously reported vibrational frequencies. All varieties of bond dissociation energies have been examined except those for C-H and P═O bonds. In trivalent phosphorus molecules, the O-C and S-C bonds have the lowest dissociation energies. In the pentavalent phosphorus set, the S-C bonds, followed by P-S bonds, have the lowest dissociation energies. In the fluorinated simulant molecules, the P-F bond is strongest, and the P-C or O-C bonds are weakest.     

Structure and hydrogen bonding in neat N-methylacetamide: classical molecular dynamics and Raman spectroscopy studies of a liquid of peptidic fragments

The results of classical molecular dynamics (MD) simulations and Raman spectroscopy studies of neat liquid N-methylacetamide (NMA), the simplest model system relevant to the peptides, are reported as a function of temperature and pressure. The MD simulations predict that near ambient conditions, the molecules form a hydrogen bond network consisting primarily of linear chains. Both the links between molecules within the hydrogen-bonded chains and the associations between chains are stabilized by weak methyl-donated "improper" hydrogen bonds. The three-dimensional structural motifs observed in the liquid show some similarity to protein beta-sheets. The temperature and pressure dependence of the hydrogen bond network, as probed by the mode frequency of the experimentally determined amide-I Raman band, blue shifts on heating and red shifts under compression, respectively, suggesting weakened and enhanced hydrogen bonding in response to temperature and pressure increases. Disruption of the hydrogen-bonding network is clearly observed in the simulation data as temperature is increased, whereas the improper hydrogen bonding is enhanced under compression to reduce the energetic cost of increasing the packing fraction. Because of the neglect of polarizability in the molecular model, the computed dielectric constant is underestimated compared to the experimental value, indicating that the simulation may underestimate dipolar coupling in the liquid.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Theoretical Prediction of Adsorption Properties of Carmustine Drug on Various Sites of the Outer Surface of the Single-Walled Boron Nitride Nanotube and Investigation of Urea Effect on Drug Delivery by DFT and MD

In this work, the adsorption behavior of Carmustine (BCNU) drug over the (6,0) zigzag single-wall boron nitride nanotube (SWBNNT) is studied by means of density functional theory calculations and molecular dynamics simulations (MD). The calculated adsorption energies proved that the adsorption of BCNU molecule on SWBNNT is a physisorption process. The natural bond orbital calculations demonstrated that existence of a charge transfer from the SWBNNT to the BCNU molecule. Moreover, quantum theory of atoms in molecules showed that the hydrogen bonds and electrostatic interactions are two major factors contributed to the overall stabilities of the complexes. Furthermore, interaction of BCNU with the surface of single wall BNNT at 310 K and 1 bar in the present of water and different concentration of Urea molecules has been studied by MD simulation. The MD results confirm that the highest number of hydrogen bond and the lowest value of Lennard-Jones (L-J) energy between nanotube and drug exist in the simulation system with concentration of 1 mol L^?1 Urea.     

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

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

Transition state dynamics of N2O4 <==> 2NO2 in liquid state

Transition state dynamics of dissociation and association reactions N2O4 <==> 2NO2 in liquid state are studied by classical molecular dynamics simulations of reactive liquid NO(2) at 298 K. An OSPP+LJ potential between NO(2) molecules proposed in paper I [J. Chem. Phys. 115, 10852 (2001)], which takes into account the orientational sensitivity of the chemical bond, has been used in the simulation. The trajectory and energy evolution of various reactions are studied in the OSPP+LJ liquid, which reproduces both the observed liquid phase equilibrium constant and Raman band shape of the dissociation mode. It is found that a NO(2) pair in reactive liquid NO(2) is bound when E(T)<0 and dissociates when E(T)>0, and the dissociation of a reactant pair occurs when the transition state (TS) surface of E(T)=0 is crossed from negative to positive, where E(T) is the sum of the potential and kinetic energies of intermolecular motion of the pair. Two types of dissociation are found depending on the source of energy for dissociation; the first type D is the dissociation via collisional activation of the reactive mode by solvent molecules, and the second type T is the dissociation via bond transfer from a dimer to a monomer NO(2) through the TS of NO(2) trimer. It is concluded that the type T dissociation is found to be much more probable than the type D dissociation because of easy energy conservation. The reactant experiences the TS of NO(2) trimer for a long time (1-10 ps) in NO(2) mediated bond transfer reactions, and crossing and recrossing trajectories and dynamics in the TS neighborhood are studied.     

热力学抑制剂作用下甲烷水合物分解过程的分子动力学模拟

采用正则系综(NVT)分子动力学方法模拟研究277.0 K、11.45 mol·L-1的热力学抑制剂乙二醇(EG)溶液作用下甲烷水合物分解微观过程. 模拟显示甲烷水合物的分解从甲烷水合物固体表面开始, 逐渐向内部推移, 固态水合物在分解过程中逐渐缩小, 直至消失. 固态水合物的分解从晶格扭曲变形开始, 之后笼形框架结构破裂, 最后形成笼形结构碎片. 同时已经分解的甲烷水合物在外层形成水膜, 包裹里层正在分解的甲烷水合物, 增大里层甲烷水合物分解传质阻力.     

On a definition of bond energies based on localized molecular orbitals

A theoretical model is presented for defining bond energies based on localized molecular Orbitals. These bond energies are obtained by rearranging the total SCF energy including the nuclear repulsion term to a sum over orbital and orbital interaction terms and then to total orbital terms, which can be interpreted as the energies of localized orbitals in a molecule. A scaling procedure is used to obtain a direct connection with experimental bond dissociation energies. Two scale parameters are employed, the C-C and the C-H bond dissociation energy in C₂H₆ for A-B and C-H type bonds, respectively. The implications of this scaling procedure are discussed. Numerical applications to a number of organic molecules containing no conjugated bonds gives in general a very satisfactory agreement between experimental and theoretical bond energies.     

Water at a hydrophilic solid surface probed by ab initio molecular dynamics: inhomogeneous thin layers of dense fluid

We present a microscopic model of the interface between liquid water and a hydrophilic, solid surface, as obtained from ab initio molecular dynamics simulations. In particular, we focused on the (100) surface of cubic SiC, a leading semiconductor candidate for biocompatible devices. Our results show that in the liquid in contact with the clean substrate, molecular dissociation occurs in a manner unexpectedly similar to that observed in the gas phase. After full hydroxylation takes place, the formation of a thin (approximately 3 A) interfacial layer is observed, which has higher density than bulk water and forms stable hydrogen bonds with the substrate. The presence of this thin layer points at rather weak effects on the structural properties of water induced by a one-dimensional confinement between approximately 1.3 nm hydrophilic substrates. In addition, our results show that the liquid does not uniformly wet the surface, but molecules preferably bind along directions parallel to the Si dimer rows.     

Homopolar- and Heteropolar Bond Dissociation Energies and Heats of Formation of Radicals and Ions in the Gas Phase. I. Data on organic molecules

The literature data on heteropolar and homopolar 2-center bond dissociation energies in organic molecules in the gas phase and the corresponding heats of formation of radicals and ions have been critically evaluated. Data for more than 500 bonds are represented in tabular form together with the pertinent literature references. Selected electron affinities and π-bond dissociation energies have also been incorporated. The follow-up paper will discuss some empirical general aspects of these data particularly regarding the effect of structure on the bond dissociation energies.     

On the structure and dynamics of the hydrated sulfite ion in aqueous solution--an ab initio QMCF MD simulation and large angle X-ray scattering study

Theoretical ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism has been applied in conjunction to experimental large angle X-ray scattering to study the structure and dynamics of the hydrated sulfite ion in aqueous solution. The results show that there is a considerable effect of the lone electron-pair on sulfur concerning structure and dynamics in comparison with the sulfate ion with higher oxidation number and symmetry of the hydration shell. The S-O bond distance in the hydrated sulfite ion has been determined to 1.53(1) ? by both methods. The hydrogen bonds between the three water molecules bound to each sulfite oxygen are only slightly stronger than those in bulk water. The sulfite ion can therefore be regarded as a weak structure maker. The water exchange rate is somewhat slower for the sulfite ion than for the sulfate ion, τ(0.5) = 3.2 and 2.6 ps, respectively. An even more striking observation in the angular radial distribution (ARD) functions is that the for sulfite ion the water exchange takes place in close vicinity of the lone electron-pair directed at its sides, while in principle no water exchange did take place of the water molecules hydrogen bound to sulfite oxygens during the simulation time. This is also confirmed when detailed pathway analysis is conducted. The simulation showed that the water molecules hydrogen bound to the sulfite oxygens can move inside the hydration shell to the area outside the lone electron-pair and there be exchanged. On the other hand, for the hydrated sulfate ion in aqueous solution one can clearly see from the ARD that the distribution of exchange events is symmetrical around the entire hydration sphere.     

An efficient method for computing excess free energy of liquid

We present a new theoretical method for efficient calculation of free energy of liquid. This interaction entropy method allows one to compute entropy and free energy of liquid from standard single step MD (molecular dynamics) simulation directly in liquid state without the need to perform MD simulations at many intermediate states as required in thermodynamic integration or free energy perturbation methods. In this new approach, one only needs to evaluate the interaction energy of a single (fixed) liquid molecule with the rest of liquid molecules as a function of time from a standard MD simulation of liquid and the fluctuation of distribution of this interaction energy is then used to calculate the interaction entropy of the liquid. Explicit theoretical derivation of this interaction entropy approach is provided and numerical calculations for the benchmark liquid water system were carried out using three different water models. Numerical analysis of the result was performed and comparison of the computational result with experimental data and other theoretical results were provided. Excellent agreement of calculated free energies with the experimental data using TIP4P model is obtained for liquid water.     

Polarization charge densities provide a predictive quantification of hydrogen bond energies

A systematic density functional theory based study of hydrogen bond energies of 2465 single hydrogen bonds has been performed. In order to be closer to liquid phase conditions, different from the usual reference state of individual donor and acceptor molecules in vacuum, the reference state of donors and acceptors embedded in a perfect conductor as simulated by the COSMO solvation model has been used for the calculation of the hydrogen bond energies. The relationship between vacuum and conductor reference hydrogen bond energies is investigated and interpreted in the light of different physical contributions, such as electrostatic energy and dispersion. A very good correlation of the DFT/COSMO hydrogen bond energies with conductor polarization charge densities of separated donor and acceptor atoms was found. This provides a method to predict hydrogen bond strength in solution with a root mean square error of 0.36 kcal mol(-1) relative to the quantum chemical dimer calculations. The observed correlation is broadly applicable and allows for a predictive quantification of hydrogen bonding, which can be of great value in many areas of computational, medicinal and physical chemistry.     

Optimization and application of lithium parameters for the reactive force field, ReaxFF

To make a practical molecular dynamics (MD) simulation of the large-scale reactive chemical systems of Li-H and Li-C, we have optimized parameters of the reactive force field (ReaxFF) for these systems. The parameters for this force field were obtained from fitting to the results of density functional theory (DFT) calculations on the structures and energy barriers for a number of Li-H and Li-C molecules, including Li(2), LiH, Li(2)H(2), H(3)C-Li, H(3)C-H(2)C-Li, H(2)C=C-LiH, HCCLi, H(6)C(5)-Li, and Li(2)C(2), and to the equations of state and lattice parameters for condensed phases of Li. The accuracy of the developed ReaxFF was also tested by comparison to the dissociation energies of lithium-benzene sandwich compounds and the collision behavior of lithium atoms with a C(60) buckyball.     

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.     

Perturbation of phospholipid bilayer properties by ethanol at a high concentration

Atomistic molecular dynamics (MD) simulations have been carried out at 30 degrees C on a fully hydrated liquid crystalline lamellar phase of dimyrystoylphosphatidylcholine (DMPC) lipid bilayer with embedded ethanol molecules at 1:1 composition, as well as on the pure bilayer phase. The ethanol molecules are found to exhibit a preference to occupy regions near the upper part of the lipid acyl chains and the phosphocholine headgroups. The calculations revealed that the phosphocholine headgroup dipoles (P- --> N+) of the lipids prefer to orient more toward the aqueous layer in the presence of ethanol. It is noticed that the ethanol molecules modify the dynamic properties of both lipids as well as the water molecules in the hydration layer of the lipid headgroups. Both the in-plane "rattling" and out-of-plane "protrusion" motions of the lipids have been found to increase in the presence of ethanol. Most importantly, it is observed that the water molecules within the hydration layer of the lipid headgroups exhibit faster translational and rotational motions in the presence of ethanol. This arises due to faster dynamics of hydrogen bonds between lipid headgroups and water in the presence of ethanol.     

MO theoretical studies of electron pair decorrelation processes. I. Pair correlation energy and single bond dissociation

The homolytic dissociation of a single bond involves the decorrelation of one electron pair. Thus, the contribution of electron correlation to dissociation energies is large. In the present paper a new procedure is presented which allows the computation of the (within the given basis) complete correlation energy of one optimized electron pair. The method which requires only modest computational effort has been applied to the calculation of dissociation energies of a number of bonds of different types. The results show that the correlation of the electron pair of the bond which is broken contributes about 50–80% to the change of the total correlation energy occuring during the dissociation process which amounts to 20–70 kcal/mol. The fraction of correlation contributed by the bond electron pair as well as the relative importance of the left-right correlation within the bond depend very much on the type of the bond. In the case of CC and CH single bonds our method yields dissociation energies which are low by only about 5 kcal/mol. Thus, the method seems to be well suited for the calculation of potential surfaces of non-concerted organic chemical reactions which involve diradicals as intermediates.     

Direct assessment of quantum nuclear effects on hydrogen bond strength by constrained-centroid ab initio path integral molecular dynamics

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.