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.

Competition Between Supramolecular Interaction and Protein–Protein Interaction in Protein Crystallization: Effects of Crystallization Method and Small Molecular Bridge

Protein crystalline frameworks with controllable interpenetration directed by dual supramolecular interactions

Multi-Scaled Explorations of Binding-Induced Folding of Intrinsically Disordered Protein Inhibitor IA3 to Its Target Enzyme

The Inducing Ligand Strategy for Supramolecular Protein Assembly

Mechanism and kinetic modeling study on the crystallization of concomitant polymorphs

Supramolecular Crystals and Crystallization with Nanosized Motifs of Giant Molecules

Crystallizing protein assemblies via free and grafted linkers

Influence of chemical modifications of the crystallophore on protein nucleating properties and supramolecular interactions network.

Crystallization of Precise Side-Chain Giant Molecules with Tunable Sequences and Functionalities

Crystallization of Small Molecules in Lyotropic Liquid Crystals

Controlling the crystal growth of DNA molecules via strategic chemical modifications

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

Solvent Effect on Molecular Conformational Evolution and Polymorphic Manipulation of Cimetidine

Competition of Hydrogen Bonds, Halogen Bonds, and π–π Interactions in Crystal Structures. Exploring the Effect of One Atom Substitution

Noncovalent Self-Assembly of Protein Crystals with Tunable Structures

Theoretical Study on Conformation Dynamics of Three-Station Molecular Shuttle in Different Environments and Its Influence on NMR Chemical Shifts and Binding Interactions.

Computational method for highly-constrained molecular dynamics of rigid bodies: coarse-grained simulation of auxetic two-dimensional protein crystals

Conformational Changes Accompanying with the Binding of Protein and Ligand

Conformational Selectivity and Evolution Affected by the Desolvation Process

Influence of Divalent Cations in the Protein Crystallization Process Assisted by Lanthanide-Based Additives