GPU-accelerated analysis and visualization of large structures solved by molecular dynamics flexible fitting.

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.DOI: http://dx.doi.org/10.7554/eLife.16105.001

Advances in cryoelectron microscopy (cryo-EM) have revolutionized the structural investigation of large macromolecular assemblies. In this review, we first provide a broad overview of modeling methods used for flexible fitting of molecular models into cryo-EM density maps. We give special attention to approaches rooted in molecular simulations-atomistic molecular dynamics and Monte Carlo. Concise descriptions of the methods are given along with discussion of their advantages, limitations, and most popular alternatives. We also describe recent extensions of the widely used molecular dynamics flexible fitting (MDFF) method and discuss how different model-building techniques could be incorporated into new hybrid modeling schemes and simulation workflows. Finally, we provide two illustrative examples of model-building and refinement strategies employing MDFF, cascade MDFF, and RosettaCM. These examples come from recent cryo-EM studies that elucidated transcription preinitiation complexes and shed light on the functional roles of these assemblies in gene expression and gene regulation.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Striated muscle contraction is regulated in a calcium-dependent manner through dynamic motions of the tropomyosin/troponin polymer, a multicomponent complex wrapped around actin-containing thin filaments. Tropomyosin/troponin sterically blocks myosin-binding at low-calcium concentrations but moves to expose myosin-binding sites at high-calcium concentrations leading to force development. Understanding the key intermolecular interactions that define these dynamic motions will promote our understanding of mutation-induced contractile dysfunction that eventually leads to hypertrophic cardiomyopathy, dilated cardiomyopathy, and skeletal myopathies. Advancements in cryoelectron microscopy (cryoEM) have resulted in a partial elucidation of structures of the thin filament, revealing many atomic-level interactions between the component proteins and critical calcium-dependent conformational alterations. However, building models at the resolutions achieved can be challenging since landmarks in the maps are often missing or ambiguous. Therefore, current computational analyses including de novo structure prediction, protein-protein docking, molecular dynamics flexible fitting, and molecular dynamics simulations are needed to ensure good quality models. We review here our efforts to model the troponin T domain spanning the head-to-tail overlap domain of tropomyosin, improving previous models. Next, we refined the published cryoEM modeled structures, which had mistakenly compressed alpha helices, with a model that has expected helical parameters while matching densities in the cryoEM volume. Lastly, we used this model to reinterpret the interactions between tropomyosin and troponin I showing key features that hold the tropomyosin cable in its low-calcium, sterically blocking position. These revised thin filament models show improved intermolecular interactions in the key low- and high-calcium regulatory states, providing novel insights into function.

An Integrative Approach to Single-Molecule FRET Spectroscopy and Molecular Dynamics Simulations for the Study of Intrinsically Disordered Proteins

"Everything that living things do can be understood in terms of jigglings and wigglings of atoms." Richard Feynman's remarks in the early 1960's summarize what is today widely accepted, namely, that biological processes can be described by the dynamics of biomolecules. Molecular dynamics (MD) simulation, in this regard, is the main methodology employed in structural biology to explore the dynamical behavior of macromolecules at a microscopic level. Aided by MD, researchers have been able, for instance, to resolve atomic structures of multi-protein complexes from cryo-EM densities, thus unveiling the atomistic details of enzymatic mechanisms and characterize the binding of small molecules to proteins. To achieve all this, the capabilities of MD packages are constantly evolving, providing a multitude of complex simulation and analysis techniques, e.g., enhanced sampling and free energy calculations. Although applicable to a wide variety of research problems, a broader usage of MD is hindered by a steep initial learning curve imposed by nearly every MD software. To reduce this initial barrier and make the methodology more accessible to the general community of biomolecular researchers, we developed an intuitive tool named QwikMD (1), which assists the users in the preparation, execution, and analysis of biomolecular MD simulations. Among many other features, QwikMD automatically checks the initial structure for structural inconsistencies, facilitates structure manipulations such as point mutations and partial deletions, simplifies the protein insertion in lipid membranes and enables the visualization and analysis of MD simulations on the fly. The user-friendly graphical interface of QwiKMD allows the preparation of MD simulations in a point-and-click fashion, offering the user multiple MD protocols, such as unbiased MD simulations, Steered MD, MD Flexible Fitting (MDFF), and, most recently, hybrid QM/MM simulations. The latter exploits the recently developed VMD and NAMD interface to common quantum mechanics software packages. QwikMD facilitates performing MD simulations for nearly any user, novice or expert. While assisting the user, QwikMD ensures reproducibility of the results by recording all parameters and steps into two log files, one in a script-like format and another in a "methods section" format. QwikMD also serves as a learning tool, providing the theoretical background of the different MD protocols and options in many "info buttons". J. V Ribeiro et al., QwikMD — Integrative Molecular Dynamics Toolkit for Novices and Experts. Sci. Rep. 6, 26536 (2016).

The idea of protein “sectors” posits that sparse subsets of amino acid residues form cooperative networks that are key elements of protein stability, ligand binding, and allosterism. To date, protein sectors have been calculated by the statistical coupling analysis (SCA) method of Ranganathan and co-workers via the spectral analysis of conservation-weighted evolutionary covariance matrices obtained from a multiple sequence alignments of homologous families of proteins. SCA sectors, a knowledge-based protocol, have been indentified with functional properties and allosterism for a number of systems. In this study, we investigate the utility of the sector idea for the analysis of physics-based molecular dynamics (MD) trajectories of proteins. Our test case for this procedure is PSD95- PDZ3, one of the smallest proteins for which allosterism has been observed. It has served previously as a model system for a number of prediction algorithms, and is well characterized by X-ray crystallography, NMR spectroscopy and site specific mutagenisis. All-atom MD simulations were performed for a total of 500 nanoseconds using AMBER, and MD-calculated covariance matrices for the fluctuations of residue displacements and non-bonded interaction energies were subjected to spectral analysis in a manner analogous to that of SCA. The composition of MD sectors was compared with results from SCA, site specific mutagenesis, and allosterism. The concordance indicates that MD sectors are a viable protocol for analyzing MD trajectories and provide insight into the physical origin of the phenomenon.Communicated by Ramaswamy H. Sarma

Cryo-electron microscopy (cryo-EM) has for a long time been a technique of choice for determining structure of large and flexible macromolecular complexes that were difficult to study by other experimental techniques such as X-ray crystallography or nuclear magnetic resonance. However, a fast development of instruments and software for cryo-EM in the last decade has allowed that a large range of complexes can be studied by cryo-EM, and that their structures can be obtained at near-atomic resolution, including the structures of small complexes (e.g., membrane proteins) whose size was earlier an obstacle to cryo-EM. Image analysis to identify multiple coexisting structures in the same specimen (multiconformation reconstruction) is now routinely done both to solve structures at near-atomic resolution and to study conformational dynamics. Methods for multiconformation reconstruction and latest examples of their applications are the focus of this review.

Scorpion venom is a rich source of toxins which have great potential to develop new therapeutic agents. Scorpion chloride channel toxins (ClTxs), such as Chlorotoxin selectively inhibit human Matrix Methaloproteinase-2 (hMMP-2). The inhibitors of hMMP-2 have potential use in cancer therapy. Three new ClTxs, meuCl14, meuCl15 and meuCl16, derived from the venom transcriptome of Iranian scorpion, M. eupeus (Buthidea family), show high sequence identity (71.4%) with Chlorotoxin. Here, 3-D homology model of new ClTxs were constructed. The models were optimized by Molecular Dynamics simulation based on MDFF (molecular dynamics flexible fitting) method. New ClTxs indicate the presence of CSαβ folding of other scorpion toxins. A docking followed by steered molecular dynamics (SMD) simulations to investigate the interactions of meuCl14, meuCl15, and meuCl16 with hMMP-2 was applied. The current study creates a correlation between the unbinding force and the inhibition activities of meuCl14, meuCl15 and meuCl16 to shed some insights as to which toxin may be used as a drug deliverer. To this aim, SMD simulations using Constant Force Pulling method were carried out. The SMD provided useful details related to the changes of electrostatic, van de Waals (vdW), and hydrogen-bonding (H-bonding) interactions between ligands and receptor during the pathway of unbinding. According to SMD results, the interaction of hMMP-2 with meuCl14 is more stable. In addition, Arginine residue was found to contribute significantly in interaction of ClTxs with hMMP-2. All in all, the present study is a dynamical approach whose results are capable of being implemented in structure-based drug design.

Molecular Dynamics Simulations of the AAA Protein P97

The structural flexibility of DNA plays a key role in many biological processes of DNA, such as protein-DNA interactions, DNA packaging in viruses and nucleosome positioning on genomic DNA. Some experimental techniques have been employed to investigate the structural flexibility of DNA with the combination of elastic models, but these experiments could only provide the macroscopic properties of DNA, and thus, it is still difficult to understand the corresponding microscopic mechanisms. Recently, all-atom molecular dynamics (MD) simulation has emerged as a useful tool to investigate not only the macroscopic properties of DNA, but also the microscopic description of the flexibility of DNA at an atomic level. The most important issue in all-atom MD simulations of DNA is to choose an appropriate force field for simulating DNA. Very recently, a new force field for DNA has been developed based on the last generation force field of Amber bsc0, which was named Amber bsc1. In this work, all-atom MD simulations are employed to study the flexibility of a 30-bp DNA with the force fields of Amber bsc1 and Amber bsc0 in a comparative way. Our aim of the research is to examine the improvement of the new development of force field (Amber bsc1) in the macroscopic and microscopic properties of DNA, in comparison with the corresponding experimental measurements. All the MD simulations are performed with Gromacs 4.6 and lasted with a simulation time of 600 ns. The MD trajectories are analyzed with Curves+ for the last 500 ns, since the system reaches equilibrium approximately after ~100 ns. Our results show that the new force field (Amber bsc1) can lead to the improvements in the macroscopic parameters of DNA flexibility, i.e., stretch modulus S and twist-stretch coupling D become closer to experimental measurements, while bending persistence lengths lp and torsional persistence lengths C from the two force fields (bsc1 and bsc0) are both in good agreement with experimental data. Our microscopic analyses show that the microscopic structure parameters of DNA from the MD simulation with the Amber bsc1 force field are closer to the experimental values than those with the Amber bsc0 force field, except for slide, and the obvious improvements are observed in some microscopic parameters such as twist and inclination. Our further analyses show that the improvements in macroscopic flexibility from the Amber bsc1 force field are tightly related to the microscopic parameters and their fluctuations. This study would be helpful in understanding the performances of Amber bsc1 and bsc0 force fields in the description of DNA flexibility at both macroscopic and microscopic level.

How motions in enzymes might be linked to catalytic function is of considerable general interest. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating functionally relevant motions in proteins that are not easily amenable to study using NMR methods. We use 3D variability analysis (3DVA) of a cryo-EM structure for human asparagine synthetase (ASNS) to identify conformational changes in a single side chain that mediates interconversion of the open and closed forms of a catalytically relevant intramolecular tunnel. Our 3DVA results are consistent with those obtained independently from molecular dynamics (MD) simulations. These MD trajectories also suggest that formation of a key reaction intermediate acts to stabilize the open form of the tunnel in ASNS, thereby permitting ammonia translocation and asparagine formation. These findings illustrate the power of cryo-EM to identify localized motions and hence dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics might regulate function in multi-domain enzymes possessing multiple active sites.

Bacterial chaperonin GroEL is a complex ring-shaped protein oligomer that promotes the folding of other proteins by encapsulating them in the cavity. There is very little structural information about the disordered C-terminal fragment of the GroEL subunits, which is involved in the folding of the substrate protein. A 3D reconstruction of the GroEL apo-form was obtained by cryo-electron microscopy (cryo-EM) with a resolution of 3.02 Å and supplemented by molecular dynamics (MD) calculations. The results of cryo-EM and MD are in good agreement and demonstrate a diverse mobility of the protein subunit domains. The MD results predict the dynamics and the network of intramolecular contacts of the C-terminal sections of the protein. These results are of great importance for the subsequent study of the mechanism of protein folding in the GroEL cavity.

Multiscale Computational Modeling of Tubulin-Tubulin Lateral Interaction

The Chikungunya virus is carried by mosquitos and can cause a number of diseases in humans including encephalitis, which can be fatal in some cases, and severe arthritis. A recent mutation in the E1 protein of the virus has allowed it to efficiently reproduce in a different species of mosquitos, leading to a Chikungunya epidemic in Réunion Island in 2005 and the subsequent infection of millions of individuals in Africa and Asia. The virus also has the potential to spread to many areas of Europe and the Americas. Chikungunya virus has a single-stranded RNA genome that codes for four non-structural proteins and five structural proteins. Based on this knowledge it has been possible to develop virus-like particles that can be used to immunise non-human primates against Chikungunya infection by inducing antibody production. However, the development of vaccines for Chikungunya in humans will require a deeper understanding of how these antibodies produced by the vaccine interact with the virus and more detailed information about the structures of the virus and antibodies. Sun et al. have used two techniques – X-ray crystallography and electron cryo-microscopy – to determine the structure of Chikungunya virus-like particles, and to obtain new insights into the interactions of these particles with four related antibodies. Electron cryo-microscopy was used to figure out the structure of the particles at near atomic resolution, and X-ray crystallography was used to determine the atomic resolution structures of two of the four Fab antibodies that neutralize the Chikungunya virus. Electron cryo-microscopy was also used to probe the complex formed by the interactions between the virus-like particles and the antibodies. Sun et al. were able to identify the likely viral receptor site that is blocked by three of the antibodies when they neutralize the virus; the fourth antibody is thought to act by immobilizing one of the domains of protein E2, thereby hiding the “fusion loop” that allows the virus to enter and infect human tissue. It is hoped that these findings will contribute to efforts to combat the spread of the Chikungunya virus worldwide.