Abstract:Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free-energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography-based restraints into the MD protocol. Because these restraints are aimed at the ensemble-average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native-like fashion as dictated by the original force field. To validate this approach, we have used the data from solid-state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well-established model protein, ubiquitin. The ensemble-restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid-state chemical shifts, and backbone order parameters. The predictions for N-15 R1 relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high-resolution crystallographic data can be viewed as protein-specific empirical corrections to the standard force fields.

What Is the Protonation State of Proteins in Crystals? Insights from Constant pH Molecular Dynamics Simulations

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Functional Protein Dynamics in a Crystal

Joint neutron crystallographic and NMR solution studies of Tyr residue ionization and hydrogen bonding: Implications for enzyme-mediated proton transfer

Protonation effects in protein-ligand complexes - a case study of endothiapepsin and pepstatin A with computational and experimental methods

Constant-pH Molecular Dynamics Simulations of Closed and Open States of a Proton-gated Ion Channel

Biomolecular simulations at constant pH

Proton first: rationalizing a proton transfer in a protein-fragment complex

Protonation State-Induced Unfolding of Protein Secondary Structure in the Gas Phase

Perspective: Structure determination of protein-ligand complexes at room temperature using X-ray diffraction approaches

Charge Interactions in a Highly Charge-depleted Protein

Ensemble MD Simulations Restrained Via Crystallographic Data: Accurate Structure Leads to Accurate Dynamics

A New Route for the Determination of Protein Structure in Physiological Environment

Nonlinear impact of electrolyte solutions on protein dynamics

Introducing titratable water to all-atom molecular dynamics at constant pH.

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications

pH-dependent structural dynamics of the neuropeptide Y in aqueous solution

Slow Conformational Exchange and Overall Rocking Motion in Ubiquitin Protein Crystals

High-Resolution Electron Diffraction of Hydrated Protein Crystals at Room Temperature

Direct determination of protonation states and visualization of hydrogen bonding in a glycoside hydrolase with neutron crystallography

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins