Site-specific hydrogen-bonding interaction between N-acetylproline amide and protic solvent molecules: comparisons of IR and VCD measurements with MD simulations

The effects of solute-solvent interactions on solution structures of small peptides have been paid a great deal of attention. To study the effect of hydrogen-bonding interactions on peptide solution structures, we measured the amide I IR and VCD spectra of N-acetylproline amide (AP) in various protic solvents, i.e., D2O, MeOD, EtOD, and PrOD, and directly compared them with theoretically simulated ones. The numbers of protic solvent molecules hydrogen-bonded to the two peptide bonds in the AP were quantitatively determined by carrying out the molecular dynamics (MD) simulations and then compared with the spectral analyses of the experimentally measured amide I bands. The two peptides in the AP have different propensities of forming H-bonds with protic solvent molecules, and the H-bond population distribution is found to be strongly site-specific and solvent-dependent. However, it is found that adoption of the polyproline II (PII) conformation by AP in protic solvents does not strongly depend on the hydrogen bond network-forming ability of protic solvents nor on the solvent polarity. We present a brief discussion on the validity as well as limitation of the currently available force field parameters used for the present MD simulation study.     

Thermal denaturation of polyalanine peptide in water by molecular dynamics simulations and theoretical prediction of infrared spectra: helix-coil transition kinetics

Perspectives in the helix-coil transition kinetics of secondary structures are examined by temperature-jump molecular dynamics (T-jump MD) simulations and theoretically calculated infrared (IR) spectra. Homopolymeric polyalanine, Ac-(A)(21)-NHMe (A21), is unfolded in water by T-jumps from 273 to 300 K approximately 450 K using AMBER ff99 and ff03 force fields. MD simulation results provide in silico evidence that 3(10)-helix and type I beta-turn motifs are highly probable in both ff99 and ff03 results. Temperature-dependent difference IR spectra of A21 do not possess an isosbestic point in both results, and isotope-labeled difference IR spectra in ff03 results predict characteristic profiles observed in experiments. Unfolding rates obtained from simulated time-resoled IR spectra are in good agreement with those estimated by helical contents, but they are still an order of magnitude smaller than experimental values. We demonstrate that the conventional criteria such as single-exponential fit of transient amide I absorbance, sigmoidal fit of temperature-dependent amide I absorbance, and Arrhenius plot of relaxation rates cannot guarantee the validity of assuming a two-state mechanism. We suggest a way of determining T(m) by the temperature dependence of center frequency and full width at half-maximum of amide I band. Overall, both ff99 and ff03 force fields give consistent results in reproducing key aspects concerned experimentally, but are not predominantly satisfactory in quantitative aspects.     

Strike a balance: optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides

Based on the AMBER polarizable model (ff02), we have re-optimized the parameters related to the main-chain (Phi, Psi) torsion angles by fitting to the Boltzmann-weighted average quantum mechanical (QM) energies of the important regions (i.e., beta, P(II), alpha(R), and alpha(L) regions). Following the naming convention of the AMBER force field series, this release will be called ff02pol.rl The force field has been assessed both by energetic comparison against the QM data and by the replica exchange molecular dynamics simulations of short alanine peptides in water. For Ace-Ala-Nme, the simulated populations in the beta, P(II) and alpha(R) regions were approximately 30, 43, and 26%, respectively. For Ace-(Ala)(7)-Nme, the populations in these three regions were approximately 24, 49, and 26%. Both were in qualitative agreement with the NMR and CD experimental conclusions. In comparison with the previous force field, ff02pol.rl demonstrated good balance among these three important regions. The optimized torsion parameters, together with those in ff02, allow us to carry out simulations on proteins and peptides with the consideration of polarization.     

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.     

A test on peptide stability of AMBER force fields with implicit solvation

We used replica exchange molecular dynamics (REMD) simulations to evaluate four different AMBER force fields and three different implicit solvent models. Our aim was to determine if these physics-based models captured the correct secondary structures of two alpha-helical and two beta-peptides: the 14-mer EK helix of Baldwin and co-workers, the C-terminal helix of ribonuclease, the 16-mer C-terminal hairpin of protein G, and the trpzip2 miniprotein. The different models gave different results, but generally we found that AMBER ff96 plus the implicit solvent model of Onufriev, Bashford, and Case gave reasonable structures, and is fairly well-balanced between helix and sheet. We also observed differences in the strength of ion pairing in the solvent models, we but found that the native secondary structures were retained even when salt bridges were prevented in the conformational sampling. Overall, this work indicates that some of these all-atom physics-based force fields may be good starting points for protein folding and protein structure prediction.     

Using PC clusters to evaluate the transferability of molecular mechanics force fields for proteins

The transferability of molecular mechanics parameters derived for small model systems to larger biopolymers such as proteins can be difficult to assess. Even for small peptides, molecular dynamics simulations are typically too short to sample structures significantly different than initial conformations, making comparison to experimental data questionable. We employed a PC cluster to generate large numbers of native and non-native conformations for peptides with experimentally measured structural data, one predominantly helical and the other forming a beta-hairpin. These atomic-detail sets do not suffer from slow convergence, and can be used to rapidly evaluate important force field properties. In this case a suspected bias toward alpha-helical conformations in the ff94 and ff99 force fields distributed with the AMBER package was verified. The sets provide critical feedback not only on force field transferability, but may also predict modifications for improvement. Such predictions were used to modify the ff99 parameter set, and the resulting force field was used to test stability and folding of model peptides. Structural behavior during molecular dynamics with the modified force field is found to be very similar to expectations, suggesting that these basis sets of conformations may themselves have significant transferability among force fields. We continue to improve and expand this data set and plan to make it publicly accessible. The calculations involved in this process are trivially parallel and can be performed using inexpensive personal computers with commodity components.     

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.     

Choice of Force Fields and Water Models for Sampling Solution Conformations of Bacteriophage T4 Lysozyme

A protein may exist as an ensemble of different conformations in solution, which cannot be represented by a single static structure. Molecular dynamics (MD) simulation has become a useful tool for sampling protein conformations in solution, but force fields and water models are important issues. This work presents a case study of the bacteriophage T4 lysozyme (T4L). We have found that MD simulations using a classic AMBER99SB force field and TIP4P water model cannot well describe hinge-bending domain motion of the wild-type T4L at the timescale of one microsecond. Other combinations, such as a residue-specific force field called RSFF2+ and a dispersion-corrected water model TIP4P-D, are able to sample reasonable solution conformations of T4L, which are in good agreement with experimental data. This primary study may provide candidates of force fields and water models for further investigating conformational transition of T4L.     

A force field for monosaccharides and (1 → 4) linked polysaccharides

A force field for monosaccharides that can be extended to (1 → 4) linked polysaccharides has been developed for the AMBER potential function. The resulting force field is consistent with the existing AMBER force field for proteins and nucleic acids. Modifications to the standard AMBER OH force constant and to the Lennard-Jones parameters were made. Furthermore, a 10–12 nonbonded term was included between the hydroxyl hydrogen of the saccharide and the water oxygen (TIP3P, SPC/E, etc.) to reproduce better the water–saccharide intermolecular distances. STO-3G electrostatic potential (ESP) charges were used to represent the electrostatic interactions between the saccharide and its surrounding environment. To obtain charges for polysaccharides, a scheme was developed to piece together saccharide residues through 1 → 4 connections while still retaining a net neutral charge on the molecule as a whole. Free energy perturbation (FEP) simulations of D -glucose and D -mannose in water were performed to test the resulting force field. The FEP simulations demonstrate that AMBER overestimates intramolecular interaction energies, suggesting that further improvements are needed in this part of the force field. To test further the reliability of the parameters, a molecular dynamics (MD) simulation of α-D -glucose in water was also performed. The MD simulation was able to produce structural and conformational results that are in accord with experimental evidence and previous theoretical results. Finally, a relaxed conformational map of β-maltose was assembled and it was found that the present force field is consistent with available theoretical and experimental results. © 1994 by John Wiley & Sons, Inc.     

Simulation of the thermodynamics of folding and unfolding of the Trp-cage mini-protein TC5b using different combinations of force fields and solvation models

Molecular dynamics simulations based on AMBER force fields(ff96 and ff03) and generalized Born models(igb1 and igb5) have been carried out in order to study folding/unfolding of the Trp-cage mini-protein TC5b.The thermodynamic properties of TC5b were found to be sensitive to the specific version of the solvation model and force field employed.When the ff96/igb5 combination was used,the predicted melting temperature from unfolding simulations was in good agreement with the experimental value of 315 K,but the...     

Computational spectroscopy of ubiquitin: comparison between theory and experiments

Using the constrained molecular dynamics simulation method in combination with quantum chemistry calculation, Hessian matrix reconstruction, and fragmentation approximation methods, the authors have established computational schemes for numerical simulations of amide I IR absorption, vibrational circular dichroism (VCD), and two-dimensional (2D) IR photon echo spectra of the protein ubiquitin in water. Vibrational characteristic features of these spectra in the amide I vibration region are discussed. From the semiempirical quantum chemistry calculation results on an isolated ubiquitin, amide I local mode frequencies and vibrational coupling constants were fully determined. It turns out that the amide I local mode frequencies of ubiquitin in both gas phase and aqueous solution are highly heterogeneous and site dependent. To directly test the quantitative validity of thus obtained spectroscopic properties, they compared the experimentally measured amide I IR, 2D IR, and electronic circular dichroism spectra with experiments, and found good agreements between theory and experiments. However, the simulated VCD spectrum is just qualitatively similar to the experimentally measured one. This indicates that, due to delicate cancellations between the positive and negative VCD contributions, the prediction of protein VCD spectrum is critically relied on quantitative accuracy of the theoretical model for predicting amide I local mode frequencies. On the basis of the present comparative investigations, they found that the site dependency of amide I local mode frequency, i.e., diagonal heterogeneity of the vibrational Hamiltonian matrix in the amide I local mode basis, is important. It is believed that the present computational methods for simulating various vibrational and electronic spectra of proteins will be of use in further refining classical force fields and in addressing the structure-spectra relationships of proteins in solution.     

Assessment of biomolecular force fields for molecular dynamics simulations in a protein crystal

Different biomolecular force fields (OPLS‐AA, AMBER03, and GROMOS96) in conjunction with SPC, SPC/E and TIP3P water models are assessed for molecular dynamics simulations in a tetragonal lysozyme crystal. The root mean square deviations for the C_a atoms of lysozymes are about 0.1 to 0.2 nm from OPLS‐AA and AMBER03, smaller than 0.4 nm from GROMOS96. All force fields exhibit similar pattern in B‐factors, whereas OPLS‐AA and AMBER03 accurately reproduce experimental measurements. Despite slight variations, the primary secondary structures are well conserved using different force fields. Water diffusion in the crystal is approximately ten‐fold slower than in bulk phase. The directional and average water diffusivities from OPLS‐AA and AMBER03 along with SPC/E model match fairly well with experimental data. Compared to GROMOS96, OPLS‐AA and AMBER03 predict larger hydrophilic solvent‐accessible surface area of lysozyme, more hydrogen bonds between lysozyme and water, and higher percentage of water in hydration shell. SPC, SPC/E and TIP3P water models have similar performance in most energetic and structural properties, but SPC/E outperforms in water diffusion. While all force fields overestimate the mobility and electrical conductivity of NaCl, a combination of OPLS‐AA for lysozyme and the Kirkwood‐Buff model for ions is superior to others. As attributed to the steric restraints and surface interactions, the mobility and conductivity in the crystal are reduced by one to two orders of magnitude from aqueous solution. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A conformational study of the trisaccharide beta-D-Glcp-(1-->2)[beta-D-Glcp-(1-->3)]alpha-D-Glcp-OMe by NMR NOESY and TROESY experiments, computer simulations, and X-ray crystal structure analysis

Proton-proton cross-relaxation rates have been measured for the trisaccharide beta-D-Glcp-(l --> 2)[beta-D-Glcp-(1 --> 3)]alpha-D-Glcp-OMe in D2O as well as in D2O/[D6]DMSO 7:3 solution at 30 degrees C by means of one-dimensional NMR pulsed field gradient 1H,1H NOESY and TROESY experiments. Interatomic distances for the trisaccharide in D2O were calculated from the cross-relaxation rates for two intraresidue and three interglycosidic proton pairs, using the isolated spin-pair approximation. In the solvent mixture one intraresidue and three interglycosidic distances were derived without the use of a specific molecular model. In this case the distances were calculated from the cross-relaxation rates in combination with "model-free" motional parameters previously derived from 13C relaxation measurements. The proton-proton distances for interglycosidic pairs were compared with those averaged from Metropolis Monte Carlo and Langevin Dynamics simulations with the HSEA, PARM22, and CHEAT95 force fields. The crystal structure of the trisaccharide was solved by analysis of X-ray data. Interresidue proton pairs from the crystal structure and those observed by NMR experiments were similar. However, the corresponding proton-proton distances generated by computer simulations were longer. For the (1 --> 2) linkage the glycosidic torsion angles of the crystal structure were found in a region of conformational space populated by all three force fields, whereas for the (1 --> 3) linkage they occupied a region of low population density, as seen from the simulations.     

Raman,infrared and force field studies of K2(12)C2O4 x H2O and K2(13)C2O4 x H2O in the solid state and in aqueous solution,and of (NH4)2(12)C2O4 x H2O and (NH4)2(13)C2O4 x H2O in the solid state

The Raman and infrared spectra of solid K2(12)C2O4 x H2O are reported together with, for the first time, the corresponding Raman and infrared spectra of solid K2(13)C2O4 x H2O. Raman spectra of aqueous solutions of both isotopomers are also reported. In the solid state the oxalate anion is planar with D2h symmetry in this salt, whereas in aqueous solution the Raman spectra of the anion are best interpreted on the basis of D2d symmetry. The Raman spectra of solid (NH4)2(12)C2O4 x H2O and (NH4)2(13)C2O4 x H2O, in which the oxalate anion is twisted from planarity by 28 degrees about the CC bond, have also been recorded. Several reassignments have been made. The harmonic force field for the oxalate anion in the D2h, D2 and D2d geometries has been determined in part, and approximate values of key valence force constants determined. All the observed band wavenumbers and 12C/13C isotopic shifts are well reproduced by the force fields. The potential energy distribution of the totally symmetric normal modes of planar oxalate indicates that each mode consists of extensively mixed symmetry corrdinates and that the labels previously used for the bands seen here at 475 and 879 cm(-1) would better be described as v(CC) and deltaS(CO2), respectively, putting them in the same wavenumber order as v(NN) and deltaS(NO2) for the isoelectronic and isostructural molecule N2O4. The stretching force constants of N2O4 and planar C2O4(2-) are established to be in the order f(NN) < f(CC) and f(NO) > f(CO), consistent with the known relative bond lengths.     

Solution structure of a peptide nucleic acid duplex from NMR data: features and limitations

This paper describes the results of a 1D and 2D NMR spectroscopy study of a palindromic 8-base pair PNA duplex GGCATGCC in H2O and H2O-D2O solutions. The (1)H NMR peaks have been assigned for most of the protons of the six central base pairs, as well as for several amide protons of the backbone. The resulting 36 interbase and base-backbone distance restraints were used together with Watson-Crick restraints to generate the PNA duplex structure in the course of 10 independent simulated annealing runs followed by restrained molecular dynamics (MD) simulations in explicit water. The resulting PNA structures correspond to a P-type helix with helical parameters close to those observed in the crystal structures of PNA. Based on the current limited number of restraints obtained from NMR spectra, alternative structures obtained by MD from starting PNA models based on DNA cannot be ruled out and are also discussed.     

栗钙土胡敏酸分子三维结构模型构建及其优化

本研究的目的在于确定胡敏酸的分子结构,探究其结构与性质的关系。综合应用了元素分析、红外光谱、13C核磁共振波谱等仪器分析方法对内蒙古草原样地栗钙土土壤的胡敏酸分子结构特征进行比较分析,并基于实验结构数据构建了土壤胡敏酸结构单元模型,通过模拟光谱与实际光谱比较验证了构建的二维结构有很大的可适性。借助Hyperchem软件,在MM+、AMBER、OPLS三种力场下应用分子模拟方法对三维单体结构进行分子力学、分子动力学结构优化。结果表明土壤胡敏酸的三维结构单体在AMBER力场下生成热最小,等于-1638.87 kJ/mol,所以其构象最为稳定。同时模拟其处于水环境下的分子结构,总势能及各种能量在优化前后有变化,但总体来说各种状态都较稳定。     

草酸桥联双核Ni(Ⅱ)配合物的合成、光谱和结构

合成并表征了两个含有不同阴离子的双核镍(Ⅱ)配合物{[(tacn)Ni(H~2O)]~2(μ-C~2O~4)}I~2·2H~2O(1)和{[(tacn)Ni(H~2O)]~2(μ-C~2O~4)}(ClO~4)~2·2H~2O(2)(tacn=1,4,7-三氮杂环壬烷)。晶体结构分析表明这两个配合物中,两个Ni离子通过草酸根桥联,每个Ni离子还与一个大环配体tacn上的三个氮原子和一个水分子配位形成变形八面体结构。结晶水和配位水之间通过氢键相连。在紫外-可见区测定了配合物的固体反射谱和溶液吸收谱。     

Density functional theory studies on tautomeric stability and infrared and Raman spectra of some purine derivatives

The molecular vibrations of 6-hydroxy-purine (6HP) and 6-amino-purine (6AP) were investigated in polycrystalline sample, at room temperature by Fourier transform infrared (FTIR) and FT-Raman spectroscopy. The spectra of the above compounds have been recorded in the region 4000-50, 3500-100 cm(-1), respectively. They were interpreted with the aid of normal coordinate analysis following full structure optimization and force field calculations based on density functional theory (DFT) using HF/6-31G* and B3LYP/6-311+G** methods and basis set combinations. The results of the calculations were applied to simulated infrared and Raman spectra of the title compounds, which showed excellent agreement with the observed spectra. The dipole moment and the tautometric stability of purine derivatives were also studied.     

Water-ketones hydrogen bonding: the rotational spectrum of cyclobutanone-water

The hydrogen-bonded complex cyclobutanone-water has been studied by Fourier-transform molecular-beam microwave spectroscopy in the frequency range of 6-18.5 GHz. The rotational spectra of ten isotopomers have been assigned and measured. Five of them have been obtained from different isotopic species (or configurations) of water (H2O, D2O, DOH, HOD, and H2 18O). The remaining five correspond to the four singly substituted 13C and to the 18O species of cyclobutanone, observed in natural abundance. For all species the inertial defect is in the range from -10.44 to -10.50 uA2, showing that the cyclobutanone frame is effectively planar and that the water molecule is coplanar to this frame. The hydrogen bond, almost linear, is formed between a water proton and one of the lone pairs of the cyclobutanone oxygen.     

Comparison of Secondary Structure Formation Using 10 Different Force Fields in Microsecond Molecular Dynamics Simulations

We have compared molecular dynamics (MD) simulations of a β-hairpin forming peptide derived from the protein Nrf2 with 10 biomolecular force fields using trajectories of at least 1 μs. The total simulation time was 37.2 μs. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. The MD simulations were done in explicit solvent with a 16-mer Nrf2 β-hairpin forming peptide using Amber ff99SB-ILDN, Amber ff99SB*-ILDN, Amber ff99SB, Amber ff99SB*, Amber ff03, Amber ff03*, GROMOS96 43a1p, GROMOS96 53a6, CHARMM27, and OPLS-AA/L force fields. The effects of charge-groups, terminal capping, and phosphorylation on the peptide folding were also examined. Despite using identical starting structures and simulation parameters, we observed clear differences among the various force fields and even between replicates using the same force field. Our simulations show that the uncapped peptide folds into a native-like β-hairpin structure at 310 K when Amber ff99SB-ILDN, Amber ff99SB*-ILDN, Amber ff99SB, Amber ff99SB*, Amber ff03, Amber ff03*, GROMOS96 43a1p, or GROMOS96 53a6 were used. The CHARMM27 simulations were able to form native hairpins in some of the elevated temperature simulations, while the OPLS-AA/L simulations did not yield native hairpin structures at any temperatures tested. Simulations that used charge-groups or peptide capping groups were not largely different from their uncapped counterparts with single atom charge-groups. On the other hand, phosphorylation of the threonine residue located at the β-turn significantly affected the hairpin formation. To our knowledge, this is the first study comparing such a large set of force fields with respect to β-hairpin folding. Such a comprehensive comparison will offer useful guidance to others conducting similar types of simulations.