Molecular Dynamics Simulations In Explicit Water Research Articles

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties.

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

Dimerization of the p53 Oligomerization Domain: Identification of a Folding Nucleus by Molecular Dynamics Simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.

In silico protein design by combinatorial assembly of protein building blocks.

Utilizing concepts of protein building blocks, we propose a de novo computational algorithm that is similar to combinatorial shuffling experiments. Our goal is to engineer new naturally occurring folds with low homology to existing proteins. A selected protein is first partitioned into its building blocks based on their compactness, degree of isolation from the rest of the structure, and hydrophobicity. Next, the protein building blocks are substituted by fragments taken from other proteins with overall low sequence identity, but with a similar hydrophobic/hydrophilic pattern and a high structural similarity. These criteria ensure that the designed protein has a similar fold, low sequence identity, and a good hydrophobic core compared with its native counterpart. Here, we have selected two proteins for engineering, protein G B1 domain and ubiquitin. The two engineered proteins share approximately 20% and approximately 25% amino acid sequence identities with their native counterparts, respectively. The stabilities of the engineered proteins are tested by explicit water molecular dynamics simulations. The algorithm implements a strategy of designing a protein using relatively stable fragments, with a high population time. Here, we have selected the fragments by searching for local minima along the polypeptide chain using the protein building block model. Such an approach provides a new method for engineering new proteins with similar folds and low homology.

All-Atom Folding Simulations of the Villin Headpiece from Stochastically Selected Coarse-Grained Structures

We discuss the application of a novel efficient protocol for the numerical simulation of the folding dynamics of single domain proteins from the only knowledge of primary sequence. Our approach is based on the combination of a Monte Carlo (MC) coarse-grained evolution followed by all-atom molecular dynamics (MD) simulations in explicit solvent. The coarse-grained model simplifies the protein's energy landscape and allows it to evolve rapidly toward viable starting conformations for MD. A general fine-graining algorithm is then used to reconstruct the full atomic detail of the protein. All atom MD simulations in explicit water are then employed to investigate the protein's conformational evolution toward the native state. We discuss the application of this novel approach to the Villin headpiece, a widely studied test system for folding studies, for which we obtain and maintain an RMS deviation from the NMR structure of 2.4 Å for the core region and 3.7 Å for the whole protein. Finally, the analysis of the MC−MD trajectories provides valuable insight into important aspects of the folding process with regards to the appearance and docking of locally secondary structure elements.

Molecular Dynamics Simulations of Surface Tensions of Aqueous Electrolytic Solutions

We present molecular-dynamics (MD) simulations of surface tension of aqueous ionic solutions of different concentrations for two different salts with different anion size: NaCl and NaF at 300 K in explicit water molecules. Our MD results agree quite closely with the experimental results. In particular, our MD results predict that the salt with a larger anion (NaCl) shows a smaller increase in surface tension than the salt with a smaller anion (NaF). To check the limits of applicability of the solvent primitive model (SPM), we solve exactly the system with the approximation that the solvent (water) is a background dielectric continuum using MD simulations. Compared with explicit water MD simulations, SPM gives qualitatively incorrect results for the dependence of surface tension on the type of salt. More precisely, with the solvent as a dielectric continuum, an NaCl solution shows a larger increase in surface tension than does an NaF solution. Strong water structuring around the solvated ions is the propo...

Dynamics and entropy of a calmodulin-peptide complex studied by NMR and molecular dynamics.

All-atom, explicit water molecular dynamics simulations of calcium-loaded calmodulin complexed with a peptide corresponding to the smooth muscle myosin light chain kinase target were carried out at 295 and 346 K. Amide and side chain methyl angular generalized order parameters were calculated and analyzed in the context of the protein's structure and dynamics. The agreement between amide order parameters measured by NMR and those from the simulations was found to be good, especially at the higher temperature, indicating both better convergence for the latter and excellent transferrability of the CHARMM parameters to the higher temperature. Subtle dynamical features such as helix fraying were reproduced. A large range of order parameters for the nine calmodulin methionines was observed in the NMR, and reproduced quite well in the simulations. The major determinant of the methionine order parameter was found to be the proximity to side chains of aromatic residues. An upper bound estimate of the difference in backbone entropy between loop and helical regions was extracted from the order parameters using a model of motion in an effective potential. Although loop regions are more flexible than helical regions, it was found that the entropy loss per residue upon folding was only approximately 20% less for loops than for helices. Pairwise correlated motions, which could significantly lower entropy estimates obtained from order parameter analysis alone, were found to be largely absent.

Comparison of the potentials of mean force for alanine tetrapeptide between integral equation theory and simulation

The dielectrically consistent reference interaction site model (DRISM) integral equation theory is applied to determine the potential of mean force (PMF) for an alanine tetramer. A stochastic dynamics simulation of the alanine tetramer using this PMF is then compared with an explicit water molecular dynamics simulation. In addition, comparison is also done with simulations using other solvent models like the extended reference interaction site model (XRISM) theory, constant dielectric and linear distance-dependent dielectric models. The results show that the DRISM method offers a fairly accurate and computationally inexpensive alternative to explicit water simulations for studies on small peptides.