Abstract:Due to hierarchic nature of biomolecular systems, their computational modeling calls for multiscale approaches, in which coarse-grained (CG) simulations are used to address long-time dynamics of large systems. Here, we review recent developments and applications of CG modeling methods, focusing on our methods primarily for proteins, DNA, and their complexes. These methods have been implemented in the CG biomolecular simulator, CafeMol. Our CG model has resolution such that ∼10 non-hydrogen atoms are grouped into one CG particle on average. For proteins, each amino acid is represented by one CG particle. For DNA, one nucleotide is simplified by three CG particles, representing sugar, phosphate, and base. The protein modeling is based on the idea that proteins have a globally funnel-like energy landscape, which is encoded in the structure-based potential energy function. We first describe two representative minimal models of proteins, called the elastic network model and the classic Go̅ model. We then present a more elaborate protein model, which extends the minimal model to incorporate sequence and context dependent local flexibility and nonlocal contacts. For DNA, we describe a model developed by de Pablo's group that was tuned to well reproduce sequence-dependent structural and thermodynamic experimental data for single- and double-stranded DNAs. Protein-DNA interactions are modeled either by the structure-based term for specific cases or by electrostatic and excluded volume terms for nonspecific cases. We also discuss the time scale mapping in CG molecular dynamics simulations. While the apparent single time step of our CGMD is about 10 times larger than that in the fully atomistic molecular dynamics for small-scale dynamics, large-scale motions can be further accelerated by two-orders of magnitude with the use of CG model and a low friction constant in Langevin dynamics. Next, we present four examples of applications. First, the classic Go̅ model was used to emulate one ATP cycle of a molecular motor, kinesin. Second, nonspecific protein-DNA binding was studied by a combination of elaborate protein and DNA models. Third, a transcription factor, p53, that contains highly fluctuating regions was simulated on two perpendicularly arranged DNA segments, addressing intersegmental transfer of p53. Fourth, we simulated structural dynamics of dinucleosomes connected by a linker DNA finding distinct types of internucleosome docking and salt-concentration-dependent compaction. Finally, we discuss many of limitations in the current approaches and future directions. Especially, more accurate electrostatic treatment and a phospholipid model that matches our CG resolutions are of immediate importance.

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

Learning Effective Molecular Models from Experimental Observables.

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

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Consistent View of Protein Fluctuations from All-Atom Molecular Dynamics and Coarse-Grained Dynamics with Knowledge-Based Force-Field.

Modeling Structural Dynamics of Biomolecular Complexes by Coarse-Grained Molecular Simulations.

Accurate computational design of three-dimensional protein crystals

Functional Protein Dynamics in a Crystal

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics.

Anisotropic coarse-grain Monte Carlo simulations of lysozyme, lactoferrin, and NISTmAb by precomputing atomistic models

Resolution-Adapted All-Atomic and Coarse-Grained Model for Biomolecular Simulations

Molecular Dynamics Simulations of Peptides and Proteins with Amplified Collective Motions.

Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding

Collective Variables for Conformational Polymorphism in Molecular Crystals

Coarse graining atomistic simulations of plastically deforming amorphous solids

Protein Simulation Using Coarse-Grained Two-Bead Multipole Force Field with Polarizable Water Models.

Molecular Dynamics Simulation of Protein Cages.

Molecular dynamics simulation of rigid molecules

Molecular-dynamics simulations of macromolecular diffraction, part I: Preparation of protein crystal simulations

Mesoscopic model for mechanical characterization of biological protein materials

Crystallizing protein assemblies via free and grafted linkers