Standard All-atom Force Fields Research Articles

Sampling of the conformational landscape of small proteins with Monte Carlo methods

Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.

Extended Charge-On-Particle Optimized Potentials for Liquid Simulation Acetone Model: The Case of Acetone-Water Mixtures.

It is well-known that classical molecular dynamics simulations of acetone-water mixtures lead to a strong phase separation when using most of the standard all-atom force fields, despite the well-known experimental fact that acetone is miscible with water in any proportion at room temperature. We describe here the use of a charge-on-particle model for accounting for the induced polarization effect in acetone-water mixtures which can solve the demixing problem at all acetone molar fractions. The polarizability effect is introduced by means of a virtual site (VS) on the carbonyl group of the acetone molecule, which increases its dipole moment and leads to a better affinity with water molecules. The VS parameter is set by fitting the density of the mixture at different acetone molar fractions. The main novelty of the VS approach lies on the transferability and universality of the model because the polarizability can be controlled without modifying the force field adopted, like previous efforts did. The results are satisfactory also in terms of the transport properties, that is, diffusivity and viscosity coefficients of the mixture.

Overcoming potential energy distortions in constrained internal coordinate molecular dynamics simulations.

The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

Direct Mixing of Atomistic Solutes and Coarse-Grained Water.

We present a new dual-resolution approach for coupling atomistic and coarse-grained models in molecular dynamics simulations of hydrated systems. In particular, a coarse-grained point dipolar water model is used to solvate molecules represented with standard all-atom force fields. A unique characteristic of our methodology is that the mixing of resolutions is direct, meaning that no additional or ad hoc scaling factors, intermediate regions, or extra sites are required. To validate the methodology, we compute the hydration free energy of 14 atomistic small molecules (analogs of amino acid side chains) solvated by the coarse-grained water. Remarkably, our predictions reproduce the experimental data as accurately as the predictions from state-of-the-art fully atomistic simulations. We also show that the hydration free energy of the coarse-grained water itself is in comparable or better agreement with the experimental value than the predictions from all but one of the most common multisite atomistic models. The coarse-grained water is then applied to solvate a typical atomistic protein containing both α-helix and β-strand elements. Moreover, parallel tempering simulations are performed to investigate the folding free energy landscape of a representative α helical and a β hairpin structure. For the simulations considered in this work, our dual-resolution method is found to be 3 to 6 times more computationally efficient than corresponding fully atomistic approaches.

Consistent free energy landscapes and thermodynamic properties of small proteins based on a single all-atom force field employing an implicit solvation

We have attempted to improve the PARAM99 force field in conjunction with the generalized Born (GB) solvation model with a surface area correction for more consistent protein folding simulations. For this purpose, using an extended alphabeta training set of five well-studied molecules with various folds (alpha, beta, and betabetaalpha), a previously modified version of PARAM99/GBSA is further refined, such that all native states of the five training species correspond to their lowest free energy minimum states. The resulting modified force field (PARAM99MOD5/GBSA) clearly produces reasonably acceptable conformational free energy surfaces of the training set with correct identifications of their native states in the free energy minimum states. Moreover, due to its well-balanced nature, this new force field is expected to describe secondary structure propensities of diverse folds in a more consistent manner. Remarkably, temperature dependent behaviors simulated with the current force field are in good agreement with the experiment. This agreement is a significant improvement over the existing standard all-atom force fields. In addition, fundamentally important thermodynamic quantities, such as folding enthalpy (DeltaH) and entropy (DeltaS), agree reasonably well with the experimental data.

Free energy surfaces of miniproteins with a ββα motif: Replica exchange molecular dynamics simulation with an implicit solvation model

Designed miniproteins with a betabetaalpha motif, such as BBA5, 1FSD, and 1PSV can serve as a benchmark set to test the validity of all-atom force fields with computer simulation, because they contain all the basic structural elements in protein folding. Unfortunately, it was found that the standard all-atom force fields with the generalized Born (GB) implicit solvation model tend to produce distorted free energy surfaces for the betabetaalpha proteins, not only because energetically those proteins need to be described by more balanced weights of the alpha- and beta-strands, but also because the GB implicit solvation model suffers from overestimated salt bridge effects. In an attempt to resolve these problems, we have modified one of the standard all-atom force fields in conjunction with the GB model, such that each native state of the betabetaalpha proteins is in its free energy minimum state with reasonable energy barriers separating local minima. With this modified energy model, the free energy contour map in each protein was constructed from the replica exchange molecular dynamics REMD simulation. The resulting free energy surfaces are significantly improved in comparison with previous simulation results and consistent with general views on small protein folding behaviors with realistic topology and energetics of all three proteins.