Kei Moritsugu

Minimum free energy path of ligand-induced transition in adenylate kinase.

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.

REACH Coarse-Grained Biomolecular Simulation: Transferability between Different Protein Structural Classes

Coarse graining of protein interactions provides a means of simulating large biological systems. The REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-graining method, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation, involves direct mapping between scales without the need for iterative optimization. Here, the transferability of the REACH force field is examined between protein molecules of different structural classes. As test cases, myoglobin (all alpha), plastocyanin (all beta), and dihydrofolate reductase (alpha/beta) are taken. The force constants derived are found to be closely similar in all three proteins. An MD version of REACH is presented, and low-temperature coarse-grained (CG) REACH MD simulations of the three proteins are compared with atomistic MD results. The mean-square fluctuations of the atomistic MD are well reproduced by the CGMD. Model functions for the CG interactions, derived by averaging over the three proteins, are also shown to produce fluctuations in good agreement with the atomistic MD. The results indicate that, similarly to the use of atomistic force fields, it is now possible to use a single, generic REACH force field for all protein studies, without having first to derive parameters from atomistic MD simulation for each individual system studied. The REACH method is thus likely to be a reliable way of determining spatiotemporal motion of a variety of proteins without the need for expensive computation of long atomistic MD simulations.

REACH coarse-grained simulation of a cellulose fiber.

A molecular level understanding of the structure, dynamics and mechanics of cellulose fibers can aid in understanding the recalcitrance of biomass to hydrolysis in cellulosic biofuel production. Here, a residue-scale REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-grained force field was derived from all-atom molecular dynamics (MD) simulations of the crystalline Iβ cellulose fibril. REACH maps the atomistic covariance matrix onto coarse-grained elastic force constants. The REACH force field was found to reproduce the positional fluctuations and low-frequency vibrational spectra from the all-atom model, allowing elastic properties of the cellulose fibril to be characterized using the coarse-grained force field with a speedup of >20 relative to atomistic MD on systems of the same size. The calculated longitudinal/transversal Youngs modulus and the velocity of sound are in agreement with experiment. The persistence length of a 36-chain cellulose microcrystal was estimated to be ~380 μm. Finally, the normal-mode analysis with the REACH force field suggests that intrinsic dynamics might facilitate the deconstruction of the cellulose fibril from the hydrophobic surface.

REACH Coarse-Grained Normal Mode Analysis of Protein Dimer Interaction Dynamics

The REACH (realistic extension algorithm via covariance Hessian) coarse-grained biomolecular simulation method is a self-consistent multiscale approach directly mapping atomistic molecular dynamics simulation results onto a residue-scale model. Here, REACH is applied to calculate the dynamics of protein-protein interactions. The intra- and intermolecular fluctuations and the intermolecular vibrational densities of states derived from atomistic molecular dynamics are well reproduced by the REACH normal modes. The phonon dispersion relations derived from the REACH lattice dynamics model of crystalline ribonuclease A are also in satisfactory agreement with the corresponding all-atom results. The REACH model demonstrates that increasing dimer interaction strength decreases the translational and rotational intermolecular vibrational amplitudes, while their vibrational frequencies are relatively unaffected. A comparative study of functionally interacting biological dimers with crystal dimers, which are formed artificially via crystallization, reveals a relation between their static structures and the interprotein dynamics: i.e., the consequence of the extensive interfaces of biological dimers is reduction of the intermonomer translational and rotational amplitudes, but not the frequencies.

Extended Phase-Space Methods for Enhanced Sampling in Molecular Simulations: A Review

Molecular Dynamics simulations are a powerful approach to study biomolecular conformational changes or protein–ligand, protein–protein, and protein–DNA/RNA interactions. Straightforward applications, however, are often hampered by incomplete sampling, since in a typical simulated trajectory the system will spend most of its time trapped by high energy barriers in restricted regions of the configuration space. Over the years, several techniques have been designed to overcome this problem and enhance space sampling. Here, we review a class of methods that rely on the idea of extending the set of dynamical variables of the system by adding extra ones associated to functions describing the process under study. In particular, we illustrate the Temperature Accelerated Molecular Dynamics (TAMD), Logarithmic Mean Force Dynamics (LogMFD), and Multiscale Enhanced Sampling (MSES) algorithms. We also discuss combinations with techniques for searching reaction paths. We show the advantages presented by this approach and how it allows to quickly sample important regions of the free-energy landscape via automatic exploration.

Extended string-like binding of the phosphorylated HP1α N-terminal tail to the lysine 9-methylated histone H3 tail

The chromodomain of HP1α binds directly to lysine 9-methylated histone H3 (H3K9me). This interaction is enhanced by phosphorylation of serine residues in the N-terminal tail of HP1α by unknown mechanism. Here we show that phosphorylation modulates flexibility of HP1α’s N-terminal tail, which strengthens the interaction with H3. NMR analysis of HP1α’s chromodomain with N-terminal tail reveals that phosphorylation does not change the overall tertiary structure, but apparently reduces the tail dynamics. Small angle X-ray scattering confirms that phosphorylation contributes to extending HP1α’s N-terminal tail. Systematic analysis using deletion mutants and replica exchange molecular dynamics simulations indicate that the phosphorylated serines and following acidic segment behave like an extended string and dynamically bind to H3 basic residues; without phosphorylation, the most N-terminal basic segment of HP1α inhibits interaction of the acidic segment with H3. Thus, the dynamic string-like behavior of HP1α’s N-terminal tail underlies the enhancement in H3 binding due to phosphorylation.

Energy Landscape of All-Atom Protein-Protein Interactions Revealed by Multiscale Enhanced Sampling

Protein-protein interactions are regulated by a subtle balance of complicated atomic interactions and solvation at the interface. To understand such an elusive phenomenon, it is necessary to thoroughly survey the large configurational space from the stable complex structure to the dissociated states using the all-atom model in explicit solvent and to delineate the energy landscape of protein-protein interactions. In this study, we carried out a multiscale enhanced sampling (MSES) simulation of the formation of a barnase-barstar complex, which is a protein complex characterized by an extraordinary tight and fast binding, to determine the energy landscape of atomistic protein-protein interactions. The MSES adopts a multicopy and multiscale scheme to enable for the enhanced sampling of the all-atom model of large proteins including explicit solvent. During the 100-ns MSES simulation of the barnase-barstar system, we observed the association-dissociation processes of the atomistic protein complex in solution several times, which contained not only the native complex structure but also fully non-native configurations. The sampled distributions suggest that a large variety of non-native states went downhill to the stable complex structure, like a fast folding on a funnel-like potential. This funnel landscape is attributed to dominant configurations in the early stage of the association process characterized by near-native orientations, which will accelerate the native inter-molecular interactions. These configurations are guided mostly by the shape complementarity between barnase and barstar, and lead to the fast formation of the final complex structure along the downhill energy landscape.

Motion Tree Delineates Hierarchical Structure of Protein Dynamics Observed in Molecular Dynamics Simulation

Molecular dynamics (MD) simulations of proteins provide important information to understand their functional mechanisms, which are, however, likely to be hidden behind their complicated motions with a wide range of spatial and temporal scales. A straightforward and intuitive analysis of protein dynamics observed in MD simulation trajectories is therefore of growing significance with the large increase in both the simulation time and system size. In this study, we propose a novel description of protein motions based on the hierarchical clustering of fluctuations in the inter-atomic distances calculated from an MD trajectory, which constructs a single tree diagram, named a “Motion Tree”, to determine a set of rigid-domain pairs hierarchically along with associated inter-domain fluctuations. The method was first applied to the MD trajectory of substrate-free adenylate kinase to clarify the usefulness of the Motion Tree, which illustrated a clear-cut dynamics picture of the inter-domain motions involving the ATP/AMP lid and the core domain together with the associated amplitudes and correlations. The comparison of two Motion Trees calculated from MD simulations of ligand-free and -bound glutamine binding proteins clarified changes in inherent dynamics upon ligand binding appeared in both large domains and a small loop that stabilized ligand molecule. Another application to a huge protein, a multidrug ATP binding cassette (ABC) transporter, captured significant increases of fluctuations upon binding a drug molecule observed in both large scale inter-subunit motions and a motion localized at a transmembrane helix, which may be a trigger to the subsequent structural change from inward-open to outward-open states to transport the drug molecule. These applications demonstrated the capabilities of Motion Trees to provide an at-a-glance view of various sizes of functional motions inherent in the complicated MD trajectory.

Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB

The multidrug transporter AcrB transports a broad range of drugs out of the cell by means of the proton-motive force. The asymmetric crystal structure of trimeric AcrB suggests a functionally rotating mechanism for drug transport. Despite various supportive forms of evidence from biochemical and simulation studies for this mechanism, the link between the functional rotation and proton translocation across the membrane remains elusive. Here, calculating the minimum free energy pathway of the functional rotation for the complete AcrB trimer, we describe the structural and energetic basis behind the coupling between the functional rotation and the proton translocation at atomic resolution. Free energy calculations show that protonation of Asp408 in the transmembrane portion of the drug-bound protomer drives the functional rotation. The conformational pathway identifies vertical shear motions among several transmembrane helices, which regulate alternate access of water in the transmembrane as well as peristaltic motions that pump drugs in the periplasm.

REACH: A program for coarse-grained biomolecular simulation☆

Abstract REACH ( R ealistic E xtension A lgorithm via C ovariance H essian) is a program package for residue-scale coarse-grained biomolecular simulation. The program calculates the force constants of a residue-scale elastic network model in single-domain proteins using the variance–covariance matrix obtained from atomistic molecular dynamics simulation. Secondary-structure dependence of the force constants is integrated. The method involves self-consistent, direct mapping of atomistic simulation results onto a coarse-grained force field in an efficient automated procedure without requiring iterative fits and avoiding system dependence. Program summary Program title: REACH Catalogue identifier: AEDA_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEDA_v1_0.html Program obtainable from: CPC Program Library, Queens University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 42 244 No. of bytes in distributed program, including test data, etc.: 3 682 118 Distribution format: tar.gz Programming language: FORTRAN 77 Computer: x86 PC Operating system: GNU/Linux, SUSE and Red Hat RAM: Depends on the system size to be calculated Word size: 32 or 64 bits Classification: 3 External routines: LAPACK, BLAS Nature of problem: A direct calculation of force field for residue-scale coarse-grained biomolecular simulation derived from atomistic molecular dynamics trajectory. Solution method: A variance–covariance matrix and the associated Hessian (second-derivative) matrix are calculated from an atomistic molecular dynamics trajectory of single-domain protein internal motion and the off-diagonal Hessian matrix is fitted to that of a residue-scale elastic network model. The resulting force constants for the residue pair interactions are expressed as model functions as a function of pairwise distance. Running time: Depends on the system size and the number of MD trajectory frames used. The test run provided with the distribution takes only a few seconds to execute.