Implicit Solvation Parameters Derived from Explicit Water Forces in Large-Scale Molecular Dynamics Simulations (original) (raw)
Related papers
Evaluation of a fast implicit solvent model for molecular dynamics simulations
Proteins: Structure, Function, and Genetics, 2001
A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail. Proteins 2002;46:24-33.
Methods in molecular biology (Clifton, N.J.)
The effects of solvation on molecular recognition are investigated from different perspectives, ranging from methods to analyse explicit solvent dynamical behaviour at the protein surface to methods for the implicit treatment of solvent effects associated with the conformational behaviour of biomolecules. The here presented implicit solvation method is based on an analytical approximation of the Solvent Accessible Surface Area (SASA) of solute molecules, which is computationally efficient and easy to parametrise. The parametrised SASA solvation method is discussed in the light of protein design and ligand binding studies. The POPS program for the SASA computation on single molecules and complex interfaces is described in detail. Explicit solvent behaviour is described here in the form of solvent density maps at the protein surface. We highlight the usefulness of that approach in defining the organisation of specific water molecules at functional sites and in determining hydrophobicity scores for the identification of potential interaction patches.
Design and application of implicit solvent models in biomolecular simulations
Current Opinion in Structural Biology, 2014
We review implicit solvent models and their parametrisation by introducing the concepts and recent devlopments of the most popular models with a focus on parametrisation via force matching. An overview of recent applications of the solvation energy term in protein dynamics, modelling, design and prediction is given to illustrate the usability and versatility of implicit solvation in reproducing the physical behaviour of biomolecular systems. Limitations of implicit modes are discussed through the example of more challenging systems like nucleic acids and membranes.
A semi-implicit solvent model for the simulation of peptides and proteins
Journal of Computational Chemistry, 2004
We present a new model of biomolecules hydration based on macroscopic electrostatic theory, that can both describe the microscopic details of solvent-solute interactions and allow for an efficient evaluation of the electrostatic hydration free energy. This semi-implicit model considers the solvent as an ensemble of polarizable pseudoparticles whose induced dipole describe both the electronic and orientational solvent polarization. In the presented version of the model, there is no mutual dipolar interaction between the particles, and they only interact through short-ranged Lennard-Jones interactions. The model has been integrated into a molecular dynamics code, and offers the possibility to simulate efficiently the conformational evolution of biomolecules. It is able to provide estimations of the electrostatic solvation free energy within short time windows during the simulation. It has been applied to the study of two small peptides, the octaalanine and the N-terminal helix of ribonuclease A, and two proteins, the bovine pancreatic trypsin inhibitor and the B1 immunoglobin-binding domain of streptococcal protein G. Molecular dynamics simulations of these biomolecules, using a slightly modified Amber force field, provide stable and meaningful trajectories in overall agreement with experiments and all-atom simulations. Correlations with respect to Poisson-Boltzmann electrostatic solvation free energies are also presented to discuss the parameterization of the model and its consequences.
Implicit solvent models for flexible protein-protein docking by molecular dynamics simulation
Proteins-structure Function and Bioinformatics, 2003
The suitability of three implicit solvent models for flexible protein–protein docking by procedures using molecular dynamics simulation is investigated. The three models are (i) the generalized Born (GB) model implemented in the program AMBER6.0; (ii) a distance-dependent dielectric (DDD) model; and (iii) a surface area-dependent model that we have parameterized and call the NPSA model. This is a distance-dependent dielectric model modified by neutralizing the ionizable side-chains and adding a surface area-dependent solvation term. These solvent models were first tested in molecular dynamics simulations at 300 K of the native structures of barnase, barstar, segment B1 of protein G, and three WW domains. These protein structures display a range of secondary structure contents and stabilities. Then, to investigate the performance of the implicit solvent models in protein docking, molecular dynamics simulations of barnase/barstar complexation, as well as PIN1 WW domain/peptide complexation, were conducted, starting from separated unbound structures. The simulations show that the NPSA model has significant advantages over the DDD and GB models in maintaining the native structures of the proteins and providing more accurate docked complexes. Proteins 2003;50:158–169. © 2002 Wiley-Liss, Inc.
Biophysical Journal, 2003
Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming ''pH memory,'' i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, ;10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.
Journal of Physical Chemistry B, 2008
The potentials of mean force of 21 heterodimers of the molecules modeling hydrophobic amino acid side chains: ethane (for alanine), propane (for proline), isobutane (for valine), isopentane (for leucine and isoleucine), ethylbenzene (for phenylalanine), methyl propyl sulfide (for methionine), and indole (for tryptophane) were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation. Analytical expressions consisting of the Gay-Berne term to represent effective van der Waals interactions and the cavity term proposed in our earlier work were fitted to the potentials of mean force. The positions and depths of the contact minima and the positions and heights of the desolvation maxima, including their dependence on the orientation of the molecules, are well represented by the analytical expressions for all systems; large deviations between the MD-determined PMF and the analytical approximations are observed for pairs involving the least spheroidal solutes: ethylbenzene, indole, and methyl propyl sulfide at short distances at which the PMF is high and, consequently, these regions are rarely visited. When data from the PMF within only 10 kcal/mol above the global minimum are considered, the standard deviation between the MD-determined and the fitted PMF is from 0.25 to 0.55 kcal/mol (the relative standard deviation being from 4% to 8%); it is larger for pairs involving nonspherical solute molecules. The free energies of contact computed from the PMF surfaces are well correlated with those determined from proteincrystal data with a slope close to that relating the free energies of transfer of amino acids (from water to n-octanol) to the average contact free energies determined from protein-crystal data. These observations justify future use of the determined potentials in coarse-grained protein-folding simulations.
Journal of Chemical Theory and Computation, 2015
Modeling of macromolecular structures and interactions represents an important challenge for computational biology, involving different time and length scales. However, this task can be facilitated through the use of coarse-grained (CG) models, which reduce the number of degrees of freedom and allow efficient exploration of complex conformational spaces. This article presents a new CG protein model named SIRAH, developed to work with explicit solvent and to capture sequence, temperature, and ionic strength effects in a topologically unbiased manner. SIRAH is implemented in GROMACS, and interactions are calculated using a standard pairwise Hamiltonian for classical molecular dynamics simulations. We present a set of simulations that test the capability of SIRAH to produce a qualitatively correct solvation on different amino acids, hydrophilic/ hydrophobic interactions, and long-range electrostatic recognition leading to spontaneous association of unstructured peptides and stable structures of single polypeptides and protein−protein complexes. a N indicates the number of residues. Values reported correspond to the average calculated over the last 100 ns of each μs. Parentheses and square brackets indicate standard deviations and values calculated from the experimental structure, respectively.
Molecular mechanics and dynamics of biomolecules using a solvent continuum model
Journal of Computational Chemistry, 2001
An easy implementation of molecular mechanics and molecular dynamics simulation using a continuum solvent model is presented that is particularly suitable for biomolecular simulations. The computation of solvation forces is made using the linear Poisson-Boltzmann equation (polar contribution) and the solvent-accessible surface area approach (nonpolar contribution). The feasibility of the methodology is demonstrated on a small protein and a small DNA hairpin. Although the parameters employed in this model must be refined to gain reliability, the performance of the method, with a standard choice of parameters, is comparable with results obtained by explicit water simulations.