Comparison of solvation-effect methods for the simulation of peptide interactions with a hydrophobic surface

In this study we investigated the interaction behavior between thirteen different small peptides and a hydrophobic surface using three progressively more complex methods of representing solvation effects: a united-atom implicit solvation method [CHARMM 19 force field (C19) with Analytical Continuum Electrostatics (ACE)], an all-atom implicit solvation method (C22 with GBMV), and an all-atom explicit solvation method (C22 with TIP3P). The adsorption behavior of each peptide was characterized by the calculation of the potential of mean force as a function of peptide-surface separation distance. The results from the C22/TIP3P model suggest that hydrophobic peptides exhibit relatively strong adsorption behavior, polar and positively-charged peptides exhibit negligible to relatively weak favorable interactions with the surface, and negatively-charged peptides strongly resist adsorption. Compared to the TIP3P model, the ACE and GBMV implicit solvent models predict much stronger attractions for the hydrophobic peptides as well as stronger repulsions for the negatively-charged peptides on the CH(3)-SAM surface. These comparisons provide a basis from which each of these implicit solvation methods may be reparameterized to provide closer agreement with explicitly represented solvation in simulations of peptide and protein adsorption to functionalized surfaces.     

Evaluation of DNA Force Fields in Implicit Solvation

DNA structural deformations and dynamics are crucial to its interactions in the cell. Theoretical simulations are essential tools to explore the structure, dynamics, and thermodynamics of biomolecules in a systematic way. Molecular mechanics force fields for DNA have benefited from constant improvements during the last decades. Several studies have evaluated and compared available force fields when the solvent is modeled by explicit molecules. On the other hand, few systematic studies have assessed the quality of duplex DNA models when implicit solvation is employed. The interest of an implicit modeling of the solvent consists in the important gain in the simulation performance and conformational sampling speed. In this study, respective influences of the force field and the implicit solvation model choice on DNA simulation quality are evaluated. To this end, extensive implicit solvent duplex DNA simulations are performed, attempting to reach both conformational and sequence diversity convergence. Structural parameters are extracted from simulations and statistically compared to available experimental and explicit solvation simulation data. Our results quantitatively expose the respective strengths and weaknesses of the different DNA force fields and implicit solvation models studied. This work can lead to the suggestion of improvements to current DNA theoretical models.     

Vibrational spectra of alpha-amino acids in the zwitterionic state in aqueous solution and the solid state: DFT calculations and the influence of hydrogen bonding

The zwitterionic forms of the two simplest alpha-amino acids, glycine and l-alanine, in aqueous solution and the solid state have been modeled by DFT calculations. Calculations of the structures in the solid state, using PW91 or PBE functionals, are in good agreement with the reported crystal structures, and the vibrational spectra computed at the optimized geometries provide a good fit to the observed IR and Raman spectra in the solid state. DFT calculations of the structures and vibrational spectra of the zwitterions in aqueous solution at the B3-LYP/cc-pVDZ level were found to require both explicit and implicit solvation models. Explicit solvation was modeled by inclusion of five hydrogen-bonded water molecules attached to each of the five possible hydrogen-bonding sites in the zwitterion and the integration equation formalism polarizable continuum model (IEF-PCM) was employed, providing a satisfactory fit to observed IR and Raman spectra. Band assignments are reported in terms of potential-energy distributions, which differ in some respects to those previously reported for glycine and l-alanine.     

2,6-二巯基嘌呤互变异构体热力学稳定性的密度泛函理论研究

利用密度泛函(DFT)B3LYP/6-311G(d,p)方法,水相计算采用自洽反应场(SCRF)中的Onsager模型,对气相和水相中可能存在的13种2,6-二巯基嘌呤互变异构体进行了全优化,并计算了各异构体的热力学参数、偶极矩及原子净电荷。计算结果表明,不论是气相还是水相,二硫酮DTP(1,3,7)是最稳定的异构体。溶剂化效应使各异构体的稳定性均增强,偶极矩大者其稳定性显著增大。溶剂化吉布斯自由能与异构体在两相中偶极矩之差存在相关性。二硫酮DTP(1,3,7)在水相中与致癌物BPDE进行亲核取代反应时,二硫酮DTP(1,3,7)中的S10原子优先进攻亲电试剂BPDE.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

The combined simulation approach of atomistic and continuum models for the thermodynamics of lysozyme crystals

We have studied the thermodynamic properties of hen egg white lysozyme crystals using a novel simulation method combining atomistic Monte Carlo simulation to calculate van der Waals interactions and the boundary element method to solve the Poisson-Boltzmann equation for the electrostatic interactions. For computational simplicity, we treat the protein as a rigid body, using the crystallographic coordinates of all non-hydrogen atoms of the protein to describe the detailed shape. NVT Monte Carlo simulations are carried out for tetragonal and orthorhombic crystals to obtain the van der Waals energy, incorporating an implicit solvation effect. For crystal phases, an optimally linearized Poisson-Boltzmann equation is used to include the effect of the Donnan equilibrium of the salt ions. The Helmholtz energy is obtained from expanded ensemble Monte Carlo simulations. By using the force field parameters that had previously been tuned for the solution properties, reasonable agreement with experiment is found for the crystallization energy of the tetragonal form. The prediction of the entropy is also reasonable with a slight underestimation suggesting the release of a few water molecules per protein during the crystallization. However, the predictions of the properties of the orthorhombic crystal are poor, probably due to differences in the solvation structure as indicated by experiments, and also as a result of the approximate force field used.     

Implicit solvent models: Combining an analytical formulation of continuum electrostatics with simple models of the hydrophobic effect

The recent development of approximate analytical formulations of continuum electrostatics opens the possibility of efficient and accurate implicit solvent models for biomolecular simulations. One such formulation (ACE, Schaefer & Karplus, J. Phys. Chem., 1996, 100:1578) is used to compute the electrostatic contribution to solvation and conformational free energies of a set of small solutes and three proteins. Results are compared to finite-difference solutions of the Poisson equation (FDPB) and explicit solvent simulations and experimental data where available. Small molecule solvation free energies agree with FDPB within 1–1.5 kcal/mol, which is comparable to differences in FDPB due to different surface treatments or different force field parameterizations. Side chain conformation free energies of aspartate and asparagine are in qualitative agreement with explicit solvent simulations, while 74 conformations of a surface loop in the protein Ras are accurately ranked compared to FDPB. Preliminary results for solvation free energies of small alkane and polar solutes suggest that a recent Gaussian model could be used in combination with analytical continuum electrostatics to treat nonpolar interactions. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 322–335, 1999     

Simultaneous detection of two triplets: a time-resolved resonance Raman study

Solvents are known to affect the triplet state structure and reactivity. In this paper, we have employed time-resolved resonance Raman (TR3) spectroscopy to understand solvent-induced subtle structural changes in the lowest excited triplet state of thioxanthone. Density functional theory (DFT) combined with the self-consistent reaction field (SCRF) implicit solvation model has been used to calculate the vibrational frequencies in the solvents. Here, we report a unique observation of the coexistence of two triplets, which has been substantiated by the probe wavelength-dependent Raman experiments. The coexistence of two triplets has been further supported by photoreduction experiments carried out at various temperatures.     

Partition coefficients of methylated DNA bases obtained from free energy calculations with molecular electron density derived atomic charges

Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have also been used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study, atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in a vacuum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases, the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform–water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER‐99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account, and when the energy needed to polarize the electron density of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogs. Comparison of the two partitioning methods, Hirshfeld‐I and Minimal Basis Iterative Stockholder (MBIS), revealed some deficiencies in the Hirshfeld‐I method related to the unstable isolated anionic nitrogen pro‐atom used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. © 2018 Wiley Periodicals, Inc.     

Hydration of l-tyrosine in aqueous medium. An experimental and theoretical study by mixed quantum mechanical/molecular mechanics methods

Quantum mechanical/molecular mechanics (QM/MM) calculations were carried out in order to study the theoretical structures of l-tyrosine in both gas phase and in aqueous solution and observe the changes that occur on the structural and vibrational properties in two phases. Therefore, the molecule was characterized by infrared and Raman spectroscopy in solid phase and aqueous solution. Optimized geometries and relative stabilities for the zwitterion l-tyrosine derivatives have been calculated taking into account the solvent effects by using the self-consistent reaction field (SCRF) theory. For a complete assignment of the IR and Raman spectra of l-tyrosine in solid and aqueous solution phases, density functional theory (DFT) calculations were combined with Pulay's scaled quantum mechanical force field (SQMFF) methodology in order to fit the theoretical wavenumber values to the experimental ones. A good agreement between theoretical and available experimental results is found.     

Contribution of Solvation Energy in Protein‐Peptide Recognition Systems

The contribution of solvation energy to guiding molecular recognition for six rigid protein‐peptide systems bad been evaluated by the variation in the number of the identified native‐like configurations and in the driving force of specific interaction resulting from the addition of the explicit solvation term in the force field function. The AMBER force field energy and the total energy including the force field energy and the WZS solvation energy were calculated for sampled configurations. The results obtained by the calculations of both force field and total energies were compared with each other. It suggests that the contribution of solvation energy is important to guiding the specific recognition of the systems in which the ligands possess larger hydrophobic or aromatic residues while the protein receptors provide the active surfaces with hydrophobic property.     

Molecular dynamics study of peptides in implicit water: ab initio folding of beta-hairpin, beta-sheet, and beta beta alpha-motif

In this communication, we have demonstrated that molecular dynamics simulations using a GB implicit solvation model with the all-atom based force field (CHARMM19) can describe the spontaneous folding of small peptides in aqueous solution. The native structures of peptides with various structural motifs (beta-hairpin, beta-sheet, and betabetaalpha-moiety) were successfully predicted within reasonable time scales by MD simulations at moderately elevated temperatures. It is expected that the present simulations provide further insight into mechanism/pathways of the peptide folding.     

Relationship between solvation energy, chemical potential and hardness variations

Starting from a density functional theory (DFT) formalism describing the energy change from one ground state representing an isolated solute, to another one representing the same solute in the field of the solvent, it is possible to obtain a simple and useful expression for the solvation energy in terms of the variation of the electronic chemical potential and global hardness associated to the change from gas to solution phase. Since both properties may be obtained from an orbital theory within the approximate Self Consistent Reaction Field (SCRF) methodology, the proposed model is expected to be useful for the analysis of chemical reactivity in solution.     

Solvent effects on tautomerism equilibria in heterocycles

. High-level ab initio quantum mechanical methods have been used to analyze the tautomeric preferences in the gas phase and in aqueous solution of three important five-member heterocycles: 4-(5-)methylimidazole, 5-hydroxyisoxazole, and 3-hydroxypyrazole. Solvent effects have been introduced by means of self-consistent reaction field (SCRF) calculations at the ab initio level using our parametrized version of the polarizable continuum model developed by Miertus, Scrocco and Tomasi (MST), including geometry relaxation upon solvation. The extent to which the MST model, and SCRF methods in general, are suitable for the study of processes of this type is discussed.     

Optimized parameters for continuum solvation calculations with carbohydrates

Continuum solvent models have been shown to be an efficient method for the calculation of the energetics of biomolecules in solution. However, for these methods to produce accurate results, an appropriate set of atomic radii or volumes is needed. While these have been developed for proteins and nucleic acids, the same is not true of carbohydrates. Here, a set of optimized parameters for continuum solvation calculations of carbohydrates in conjunction with the Carbohydrate Solution Force Field are presented. Explicit solvent free-energy perturbation simulations were performed on a set of hexapyranose sugars and used to fit atomic radii for Poisson-Boltzmann and generalized-Born calculations, and to fit atomic volumes for use with the analytical continuum electrostatics model. The solvation energetics computed with the optimized radii and a Poisson-Boltzmann model show remarkable agreement with explicit solvent simulation, with a root-mean-square error of 1.19 kcal/mol over a large test set of sugars in many conformations. The generalized-Born model gives slightly poorer agreement, but still correlates very strongly, with an error of 1.69 kcal/mol. The analytical continuum electrostatics model correlates well with the explicit solvent results, but gives a larger error of 4.71 kcal/mol. The remarkable agreement between the solvation free energies computed in explicit and implicit solvent provides strong motivation for the use of implicit solvent models in the simulation of carbohydrate-containing systems.     

Theoretical determination of the standard reduction potentials of pheophytin-a in N,N-dimethyl formamide and membrane

Quantum mechanical/molecular mechanics (QM/MM) calculations were performed on the neutral, anionic, and dianionic forms of Pheophytin-a (Pheo-a) in N,N-dimethyl formamide (DMF) in order to calculate the absolute free energy of reduction of Pheo-a in solution. The geometry of the solvated species was optimized by restricted open-shell density functional treatment (ROB3LYP) using the 6-31G(d) basis set for the molecular species while the primary solvent shell consisting of 45 DMF molecules was treated by the MM method using the universal force field (UFF). Electronic energies of the neutral, anionic, and dianionic species were obtained by carrying out single point density functional theory (DFT) calculations using the 6-311+G(2d,2p) basis set on the respective ONIOM optimized geometries. The CHARMM27 force field was used to account for the dynamical nature of the primary solvation shell of 45 DMF molecules. In the calculations using solvent shells, the required atomic charges for each solvent molecule were obtained from ROB3LYP/6-31G(d) calculation on a single isolated DMF molecule. Randomly sampled configurations generated by the Monte Carlo (MC) technique were used to determine the contribution of the primary shell to the free energy of solvation of the three species. The dynamical nature of the primary shell significantly corrects the free energy of solvation. Frequency calculations at the ROB3LYP/6-31G(d) level were carried out on the optimized geometries of truncated 47-atom models of the neutral and ionic species in vacuum so as to determine the differences in thermal energy and molecular entropy. The Born energy of ion-dielectric interaction, the Onsager energy of dipole-dielectric interaction, and the Debye-Hückel energy of ion-ionic cloud interaction for the pheophytin-solvent aggregate were added as perturbative corrections. The Born interaction also makes a large contribution to the absolute free energy of reduction. An implicit solvation model (DPCM) was also employed for the calculation of standard reduction potentials in DMF. Both the models were successful in reproducing the standard reduction potentials. An explicit solvent treatment(QM/MM/MC + Born + Onsager + Debye corrections) yielded the one electron reduction potential of Pheo-a as -0.92 +/- 0.27 V and the two electron reduction potential as -1.34 +/- 0.25 V at 298.15 K, while the implicit solvent treatment yielded the corresponding values as -1.03 +/- 0.17 and -1.30 +/- 0.17 V, respectively. The calculated values more or less agree with the experimental midpoint potentials of -0.90 and -1.25 V, respectively. Moreover, a numerical finite difference Poisson-Boltzmann solver (FDPB) along with the DPCM methodology was employed to obtain the reduction potential of pheophytin in the thylakoid membrane. The calculated reduction potential value of -0.58 V is in excellent agreement with the reported value -0.61 V.