Modeling Human Cardiac Thin Filament Structures

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.

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

"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).

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.

Hybrid structure fitting methods combine data from cryo-electron microscopy and X-ray crystallography with molecular dynamics simulations for the determination of all-atom structures of large biomolecular complexes. Evaluating the quality-of-fit obtained from hybrid fitting is computationally demanding, particularly in the context of a multiplicity of structural conformations that must be evaluated. Existing tools for quality-of-fit analysis and visualization have previously targeted small structures and are too slow to be used interactively for large biomolecular complexes of particular interest today such as viruses or for long molecular dynamics trajectories as they arise in protein folding. We present new data-parallel and GPU-accelerated algorithms for rapid interactive computation of quality-of-fit metrics linking all-atom structures and molecular dynamics trajectories to experimentally-determined density maps obtained from cryo-electron microscopy or X-ray crystallography. We evaluate the performance and accuracy of the new quality-of-fit analysis algorithms vis-à-vis existing tools, examine algorithm performance on GPU-accelerated desktop workstations and supercomputers, and describe new visualization techniques for results of hybrid structure fitting methods.

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

Molecular Dynamics Simulations of the AAA Protein P97

The sarcomere is the fundamental unit of cardiac and skeletal muscle contraction. During the last ten years, there has been growing awareness of the etiology of skeletal and cardiac muscle diseases originating in the sarcomere, an important evolving field. Many sarcomeric diseases affect newborn children, i. e. are congenital myopathies. The discovery and characterization of several myopathies caused by mutations in myosin heavy chain genes, coding for the major component of skeletal muscle thick filaments, has led to the introduction of a new entity in the field of neuromuscular disorders: myosin myopathies. Recently, mutations in genes coding for skeletal muscle thin filaments, associated with various clinical features, have been identified. These mutations evoke distinct structural changes within the sarcomeric thin filament. Current knowledge regarding contractile protein dysfunction as it relates to disease pathogenesis has failed to decipher the mechanistic links between mutations identified in sarcomeric proteins and skeletal myopathies, which will no doubt require an integrated physiological approach. The discovery of additional genes associated with myopathies and the elucidation of the molecular mechanisms of pathogenesis will lead to improved and more accurate diagnosis, including prenatally, and to enhanced potential for prognosis, genetic counseling and developing possible treatments for these diseases. The goal of this review is to present recent progress in the identification of gene mutations from each of the major structural components of the sarcomere, the thick and thin filaments, related to skeletal muscle disease. The genetics and clinical manifestations of these disorders will be discussed.

Human α-defensins are proteins of the innate immune system that suppress viral and bacterial infections by multiple mechanisms including membrane disruption. For viruses that lack envelopes, such as human adenovirus (HAdV), other, less well defined, mechanisms must be involved. A previous structural study on the interaction of an α-defensin, human α-defensin 5 (HD5), with HAdV led to a proposed mechanism in which HD5 stabilizes the vertex region of the capsid and blocks uncoating steps required for infectivity. Studies with virus chimeras comprised of capsid proteins from sensitive and resistant serotypes supported this model. To further characterize the critical binding site, we determined subnanometer resolution cryo-electron microscopy (cryoEM) structures of HD5 complexed with both neutralization-sensitive and -resistant HAdV chimeras. Models were built for the vertex regions of these chimeras with monomeric and dimeric forms of HD5 in various initial orientations. CryoEM guided molecular dynamics flexible fitting (MDFF) was used to restrain the majority of the vertex model in well-defined cryoEM density. The RGD-containing penton base loops of both the sensitive and resistant virus chimeras are predicted to be intrinsically disordered, and little cryoEM density is observed for them. In simulations these loops from the sensitive virus chimera, interact with HD5, bridge the penton base and fiber proteins, and provides significant stabilization with a three-fold increase in the intermolecular nonbonded interactions of the vertex complex. In the case of the resistant virus chimera, simulations revealed fewer bridging interactions and reduced stabilization by HD5. This study implicates a key dynamic region in mediating a stabilizing interaction between a viral capsid and a protein of the innate immune system with potent anti-viral activity.

Skeletal muscle actin (ACTA1) mutations are a prevalent cause of skeletal myopathies consistent with ACTA1’s high expression in skeletal muscle. Rare de novo mutations in ACTA1 associated with combined cardiac and skeletal myopathies have been reported, but ACTA1 represents only ~20% of the total actin pool in cardiomyocytes, making its role in cardiomyopathy controversial. Here we demonstrate how a mutation in an actin isoform expressed at low levels in cardiomyocytes can cause cardiomyopathy by focusing on a unique ACTA1 mutation, R256H. We previously identified this mutation in multiple family members with dilated cardiomyopathy (DCM), who had reduced systolic function without clinical skeletal myopathy. Using a battery of multiscale biophysical tools, we show that R256H has potent functional effects on ACTA1 function at the molecular scale and in human cardiomyocytes. Importantly, we demonstrate that R256H acts in a dominant manner, where the incorporation of small amounts of mutant protein into thin filaments is sufficient to disrupt molecular contractility, and that this effect is dependent on the presence of troponin and tropomyosin. To understand the structural basis of this change in regulation, we resolved a structure of R256H filaments using Cryo-EM, and we see alterations in actin’s structure that have the potential to disrupt interactions with tropomyosin. Finally, we show that ACTA1R256H/+ human induced pluripotent stem cell cardiomyocytes demonstrate reduced contractility and sarcomeric disorganization. Taken together, we demonstrate that R256H has multiple effects on ACTA1 function that are sufficient to cause reduced contractility and establish a likely causative relationship between ACTA1 R256H and clinical cardiomyopathy.

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π⋯π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN⋯Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50 ns MD simulation, the π⋯π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.Communicated by Ramaswamy H. Sarma

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Computational modeling of cooperative regulation of both strands in thin filaments.

Cryo-electron microscopy (EM) and X-ray studies proposed different mechanisms for annexin-induced membrane aggregation. In this work, molecular dynamics (MD) simulation technique was utilized to gain an insight into the calcium-induced conformational changes on annexin I and their implication in membrane aggregation mechanism. MD simulations were performed on the Ca2+ -free annexin I with the N-terminal domain buried inside the core (System 1), the Ca2+ -bound annexin I without N-terminal domain (System 2) and the Ca2+ -bound annexin I with the N-terminal domain exposed (System 3). Our results indicated that calcium binding increases the flexibility of annexin I core domain residues including the calcium coordinating residues. As a result, annexin I was activated to interact with the negatively charged membrane. The exposed N-terminal domain was very flexible and gradually lost the secondary structure during MD simulation, suggesting that the N-terminal may adopt a favorable conformation to bind a second membrane and also explaining the failure of attempts to crystallize the full-length annexin I in the presence of calcium ions. The measured dimensions of the averaged simulation structure of the Ca2+ -bound annexin I with the N-terminal exposed (System 3) support the proposed membrane aggregation mechanism based on X-ray studies.

Cryo_fit: Democratization of flexible fitting for cryo-EM