Molecular dynamics-based model refinement and validation for sub-5 Å cryo-electron microscopy maps
- University of Illinois at Urbana Champaign, United States
- University of Illinois at Urbana-Champaign, United States
- University of California, San Francisco, United States
Abstract
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 convergence radius of ~25Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2Å and 3.4Å resolution. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the respective macromolecules studied, captured employing local root mean square fluctuations. The MDFF tools are made available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.
Data availability
-
NIZN[FE4S4] AND NINI[FE4S4] CLUSTERS IN CLOSED AND OPEN ALPHA SUBUNITS OF ACETYL-COA SYNTHASE/CARBON MONOXIDE DEHYDROGENASEPublicly available at the Protien Data Bank (accession no. 1OAO).
-
Structure of TRPV1 ion channel determined by single particle electron cryo-microscopyPublicly available at the Protien Data Bank (accession no. 3J5P).
-
Structure of the capsaicin receptor, TRPV1, determined by single particle electron cryo-microscopyPublicly available at the EMDataBank (accesion no. EMD-5778).
-
2.2 A resolution cryo-EM structure of beta-galactosidase in complex with a cell-permeant inhibitorPublicly available at the Protien Data Bank (accession no. 5A1A).
-
2.2 A resolution cryo-EM structure of beta-galactosidase in complex with a cell-permeant inhibitorPublicly available at the EMDataBank (accesion no. EMD-2984).
-
Cryo-EM structure of the human gamma-secretase complex at 3.4 angstrom resolution.Publicly available at the Protien Data Bank (accession no. 5A63).
-
Cryo-EM structure of the human gamma-secretase complex at 3.4 angstrom resolutionPublicly available at the EMDataBank (accesion no. EMD-3061).
-
Structure of a extracellular domainPublicly available at the Protien Data Bank (accession no. 4UPC).
-
Three-dimensional structure of human gamma-secretase at 4.5 angstrom resolutionPublicly available at the EMDataBank (accesion no. EMD-2677).
-
Thermoplasma acidophilum 20S proteasomePublicly available at the Protien Data Bank (accession no. 3J9I).
-
3D reconstruction of archaeal 20S proteasomePublicly available at the EMDataBank (accesion no. EMD-5623).
-
Structure of beta-galactosidase at 3.2-A resolution obtained by cryo-electron microscopyPublicly available at the Protien Data Bank (accession no. 3J7H).
-
Structure of beta-galactosidase at 3.2-A resolution obtained by cryo-electron microscopyPublicly available at the EMDataBank (accesion no. EMD-5995).
Article and author information
Author details
Copyright
© 2016, singharoy et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,806
- views
-
- 1,096
- downloads
-
- 144
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Molecular dynamics-based model refinement and validation for sub-5 Å cryo-electron microscopy maps
eLife 5:e16105.
https://doi.org/10.7554/eLife.16105
Further reading
-
- Neuroscience
- Structural Biology and Molecular Biophysics
We present near-atomic-resolution cryoEM structures of the mammalian voltage-gated potassium channel Kv1.2 in open, C-type inactivated, toxin-blocked and sodium-bound states at 3.2 Å, 2.5 Å, 3.2 Å, and 2.9 Å. These structures, all obtained at nominally zero membrane potential in detergent micelles, reveal distinct ion-occupancy patterns in the selectivity filter. The first two structures are very similar to those reported in the related Shaker channel and the much-studied Kv1.2–2.1 chimeric channel. On the other hand, two new structures show unexpected patterns of ion occupancy. First, the toxin α-Dendrotoxin, like Charybdotoxin, is seen to attach to the negatively-charged channel outer mouth, and a lysine residue penetrates into the selectivity filter, with the terminal amine coordinated by carbonyls, partially disrupting the outermost ion-binding site. In the remainder of the filter two densities of bound ions are observed, rather than three as observed with other toxin-blocked Kv channels. Second, a structure of Kv1.2 in Na+ solution does not show collapse or destabilization of the selectivity filter, but instead shows an intact selectivity filter with ion density in each binding site. We also attempted to image the C-type inactivated Kv1.2 W366F channel in Na+ solution, but the protein conformation was seen to be highly variable and only a low-resolution structure could be obtained. These findings present new insights into the stability of the selectivity filter and the mechanism of toxin block of this intensively studied, voltage-gated potassium channel.
-
- Structural Biology and Molecular Biophysics
Aberrant signaling of BRAFV600E is a major cancer driver. Current FDA-approved RAF inhibitors selectively inhibit the monomeric BRAFV600E and suffer from tumor resistance. Recently, dimer-selective and equipotent RAF inhibitors have been developed; however, the mechanism of dimer selectivity is poorly understood. Here, we report extensive molecular dynamics (MD) simulations of the monomeric and dimeric BRAFV600E in the apo form or in complex with one or two dimer-selective (PHI1) or equipotent (LY3009120) inhibitor(s). The simulations uncovered the unprecedented details of the remarkable allostery in BRAFV600E dimerization and inhibitor binding. Specifically, dimerization retrains and shifts the αC helix inward and increases the flexibility of the DFG motif; dimer compatibility is due to the promotion of the αC-in conformation, which is stabilized by a hydrogen bond formation between the inhibitor and the αC Glu501. A more stable hydrogen bond further restrains and shifts the αC helix inward, which incurs a larger entropic penalty that disfavors monomer binding. This mechanism led us to propose an empirical way based on the co-crystal structure to assess the dimer selectivity of a BRAFV600E inhibitor. Simulations also revealed that the positive cooperativity of PHI1 is due to its ability to preorganize the αC and DFG conformation in the opposite protomer, priming it for binding the second inhibitor. The atomically detailed view of the interplay between BRAF dimerization and inhibitor allostery as well as cooperativity has implications for understanding kinase signaling and contributes to the design of protomer selective RAF inhibitors.
