1
|
Walton T, Doran MH, Brown A. Structural determination and modeling of ciliary microtubules. Acta Crystallogr D Struct Biol 2024; 80:220-231. [PMID: 38451206 PMCID: PMC10994176 DOI: 10.1107/s2059798324001815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Accepted: 02/24/2024] [Indexed: 03/08/2024] Open
Abstract
The axoneme, a microtubule-based array at the center of every cilium, has been the subject of structural investigations for decades, but only recent advances in cryo-EM and cryo-ET have allowed a molecular-level interpretation of the entire complex to be achieved. The unique properties of the nine doublet microtubules and central pair of singlet microtubules that form the axoneme, including the highly decorated tubulin lattice and the docking of massive axonemal complexes, provide opportunities and challenges for sample preparation, 3D reconstruction and atomic modeling. Here, the approaches used for cryo-EM and cryo-ET of axonemes are reviewed, while highlighting the unique opportunities provided by the latest generation of AI-guided tools that are transforming structural biology.
Collapse
Affiliation(s)
- Travis Walton
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Matthew H. Doran
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Alan Brown
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
2
|
Chen Z, Shiozaki M, Haas KM, Skinner WM, Zhao S, Guo C, Polacco BJ, Yu Z, Krogan NJ, Lishko PV, Kaake RM, Vale RD, Agard DA. De novo protein identification in mammalian sperm using in situ cryoelectron tomography and AlphaFold2 docking. Cell 2023; 186:5041-5053.e19. [PMID: 37865089 PMCID: PMC10842264 DOI: 10.1016/j.cell.2023.09.017] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Revised: 08/02/2023] [Accepted: 09/16/2023] [Indexed: 10/23/2023]
Abstract
To understand the molecular mechanisms of cellular pathways, contemporary workflows typically require multiple techniques to identify proteins, track their localization, and determine their structures in vitro. Here, we combined cellular cryoelectron tomography (cryo-ET) and AlphaFold2 modeling to address these questions and understand how mammalian sperm are built in situ. Our cellular cryo-ET and subtomogram averaging provided 6.0-Å reconstructions of axonemal microtubule structures. The well-resolved tertiary structures allowed us to unbiasedly match sperm-specific densities with 21,615 AlphaFold2-predicted protein models of the mouse proteome. We identified Tektin 5, CCDC105, and SPACA9 as novel microtubule-associated proteins. These proteins form an extensive interaction network crosslinking the lumen of axonemal doublet microtubules, suggesting their roles in modulating the mechanical properties of the filaments. Indeed, Tekt5 -/- sperm possess more deformed flagella with 180° bends. Together, our studies presented a cellular visual proteomics workflow and shed light on the in vivo functions of Tektin 5.
Collapse
Affiliation(s)
- Zhen Chen
- Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, CA, USA; Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA.
| | - Momoko Shiozaki
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Kelsey M Haas
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA; J. David Gladstone Institutes, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, USA
| | - Will M Skinner
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Shumei Zhao
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Caiying Guo
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Benjamin J Polacco
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, USA
| | - Zhiheng Yu
- Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Nevan J Krogan
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA; J. David Gladstone Institutes, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, USA
| | - Polina V Lishko
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Robyn M Kaake
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA; J. David Gladstone Institutes, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, USA
| | - Ronald D Vale
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, CA, USA; Janelia Research Campus, Howard Hughes Medical Institute, Ashburn, VA, USA.
| | - David A Agard
- Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, CA, USA; Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, USA.
| |
Collapse
|
3
|
Tworek JW, Elcock AH. Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides a Fast but Still Accurate Treatment of Hydrodynamic Interactions in Brownian Dynamics Simulations of Biological Macromolecules. J Chem Theory Comput 2023; 19:5099-5111. [PMID: 37409946 PMCID: PMC10413861 DOI: 10.1021/acs.jctc.3c00476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Indexed: 07/07/2023]
Abstract
The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HIs). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here, we explore the use of an alternative way to accelerate the calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA-RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA-RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast, approximate way of including HIs in BD simulations of intermediate-scale systems.
Collapse
Affiliation(s)
- John W. Tworek
- Department of Biochemistry
& Molecular Biology, University of Iowa, Iowa City, Iowa 52242, United States
| | - Adrian H. Elcock
- Department of Biochemistry
& Molecular Biology, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
4
|
Tworek JW, Elcock AH. An Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides A Fast But Still Accurate Treatment Of Hydrodynamic Interactions In Brownian Dynamics Simulations Of Biological Macromolecules. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.21.537865. [PMID: 37162930 PMCID: PMC10168278 DOI: 10.1101/2023.04.21.537865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HI). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here we explore the use of an alternative way to accelerate calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast approximate way of including HIs in BD simulations of intermediate-scale systems.
Collapse
|
5
|
Mei K, Guo W. Modeling the Cryo-EM Structure of the Exocyst Complex. Methods Mol Biol 2023; 2557:247-262. [PMID: 36512220 DOI: 10.1007/978-1-0716-2639-9_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Multi-subunit tethering complexes (MTCs) are a family of evolutionarily conserved large protein complexes that function to tether intracellular vesicles from the donor compartments to the membrane of receptor compartments. The exocyst complex is an octameric MTC that tethers the post-Golgi secretory vesicles to the plasma membrane. To learn the function and regulation of the exocyst complex, it is crucial to understand the structure of the complex. We have solved the cryo-EM structure of the exocyst complex at 4.4 Angstrom (Å) resolution and detected the spatial relationship between the eight subunits using chemical cross-linking mass spectrometry. Here, we describe the method of modeling and validating the cryo-EM structure of the exocyst complex. This method could provide a guide for modeling of other protein complexes of which the structures are solved at medium to near-atomic resolution.
Collapse
Affiliation(s)
- Kunrong Mei
- School of Pharmaceutical Science and Technology, Tianjin University, Tianjin, China.
| | - Wei Guo
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
6
|
Beton JG, Cragnolini T, Kaleel M, Mulvaney T, Sweeney A, Topf M. Integrating model simulation tools and
cryo‐electron
microscopy. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Joseph George Beton
- Centre for Structural Systems Biology (CSSB) Leibniz‐Institut für Virologie (LIV) Hamburg Germany
| | - Tristan Cragnolini
- Institute of Structural and Molecular Biology, Birkbeck and University College London London UK
| | - Manaz Kaleel
- Centre for Structural Systems Biology (CSSB) Leibniz‐Institut für Virologie (LIV) Hamburg Germany
| | - Thomas Mulvaney
- Centre for Structural Systems Biology (CSSB) Leibniz‐Institut für Virologie (LIV) Hamburg Germany
| | - Aaron Sweeney
- Centre for Structural Systems Biology (CSSB) Leibniz‐Institut für Virologie (LIV) Hamburg Germany
| | - Maya Topf
- Centre for Structural Systems Biology (CSSB) Leibniz‐Institut für Virologie (LIV) Hamburg Germany
| |
Collapse
|
7
|
Abstract
Adeno-associated virus (AAV) has a single-stranded DNA genome encapsidated in a small icosahedrally symmetric protein shell with 60 subunits. AAV is the leading delivery vector in emerging gene therapy treatments for inherited disorders, so its structure and molecular interactions with human hosts are of intense interest. A wide array of electron microscopic approaches have been used to visualize the virus and its complexes, depending on the scientific question, technology available, and amenability of the sample. Approaches range from subvolume tomographic analyses of complexes with large and flexible host proteins to detailed analysis of atomic interactions within the virus and with small ligands at resolutions as high as 1.6 Å. Analyses have led to the reclassification of glycan receptors as attachment factors, to structures with a new-found receptor protein, to identification of the epitopes of antibodies, and a new understanding of possible neutralization mechanisms. AAV is now well-enough characterized that it has also become a model system for EM methods development. Heralding a new era, cryo-EM is now also being deployed as an analytic tool in the process development and production quality control of high value pharmaceutical biologics, namely AAV vectors.
Collapse
Affiliation(s)
- Scott
M. Stagg
- Department
of Biological Sciences, Florida State University, Tallahassee, Florida 32306, United States
- Institute
of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, United States
| | - Craig Yoshioka
- Department
of Biomedical Engineering, Oregon Health
& Science University, Portland Oregon 97239, United States
| | - Omar Davulcu
- Environmental
Molecular Sciences Laboratory, Pacific Northwest
National Laboratory, 3335 Innovation Boulevard, Richland, Washington 99354, United States
| | - Michael S. Chapman
- Department
of Biochemistry, University of Missouri, Columbia, Missouri 65211, United States
| |
Collapse
|
8
|
Krieger JM, Sorzano COS, Carazo JM, Bahar I. Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0. ACTA CRYSTALLOGRAPHICA SECTION D STRUCTURAL BIOLOGY 2022; 78:399-409. [PMID: 35362464 PMCID: PMC8972803 DOI: 10.1107/s2059798322001966] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 02/18/2022] [Indexed: 11/24/2022]
Abstract
New computational biophysics pipelines for analysing the global dynamics of structural ensembles and large, dynamic complexes resolved by cryoEM are reviewed. Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.
Collapse
|
9
|
Bandyopadhyay H, Deng Z, Ding L, Liu S, Uddin MR, Zeng X, Behpour S, Xu M. Cryo-shift: reducing domain shift in cryo-electron subtomograms with unsupervised domain adaptation and randomization. Bioinformatics 2022; 38:977-984. [PMID: 34897387 DOI: 10.1093/bioinformatics/btab794] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Revised: 10/18/2021] [Accepted: 11/17/2021] [Indexed: 02/05/2023] Open
Abstract
MOTIVATION Cryo-Electron Tomography (cryo-ET) is a 3D imaging technology that enables the visualization of subcellular structures in situ at near-atomic resolution. Cellular cryo-ET images help in resolving the structures of macromolecules and determining their spatial relationship in a single cell, which has broad significance in cell and structural biology. Subtomogram classification and recognition constitute a primary step in the systematic recovery of these macromolecular structures. Supervised deep learning methods have been proven to be highly accurate and efficient for subtomogram classification, but suffer from limited applicability due to scarcity of annotated data. While generating simulated data for training supervised models is a potential solution, a sizeable difference in the image intensity distribution in generated data as compared with real experimental data will cause the trained models to perform poorly in predicting classes on real subtomograms. RESULTS In this work, we present Cryo-Shift, a fully unsupervised domain adaptation and randomization framework for deep learning-based cross-domain subtomogram classification. We use unsupervised multi-adversarial domain adaption to reduce the domain shift between features of simulated and experimental data. We develop a network-driven domain randomization procedure with 'warp' modules to alter the simulated data and help the classifier generalize better on experimental data. We do not use any labeled experimental data to train our model, whereas some of the existing alternative approaches require labeled experimental samples for cross-domain classification. Nevertheless, Cryo-Shift outperforms the existing alternative approaches in cross-domain subtomogram classification in extensive evaluation studies demonstrated herein using both simulated and experimental data. AVAILABILITYAND IMPLEMENTATION https://github.com/xulabs/aitom. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Hmrishav Bandyopadhyay
- Department of Electronics and Telecommunication Engineering, Jadavpur University, Kolkata 700032, India
| | - Zihao Deng
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Leiting Ding
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Sinuo Liu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Mostofa Rafid Uddin
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Xiangrui Zeng
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Sima Behpour
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| |
Collapse
|
10
|
Krakow J, Hammel M, Zhu Y, Hillier BJ, Paolella B, Desmarais A, Wall R, Chen THT, Pei R, Karunatilake C, DuBridge R, Vinogradova M. Structural arrangement of the VH and VL domains in the COBRA™ T-cell engaging single-chain diabody. Antib Ther 2022; 5:1-10. [PMID: 35005430 PMCID: PMC8719580 DOI: 10.1093/abt/tbab028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Revised: 11/19/2021] [Accepted: 11/22/2021] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND COBRA™ (COnditional Bispecific Redirected Activation) T-cell engagers are designed to target solid tumors as a single polypeptide chain prodrug that becomes activated by proteolysis in the tumor microenvironment. One COBRA molecule comprises seven Ig domains: three single-domain antibodies (sdAbs) recognizing a tumor target or human serum albumin (HSA), and CD3ε-binding variable fragment heavy chain (VH) and variable fragment light chain (VL) and their inactivated counterparts, VHi and VLi. Pairing of VH and VL, and VLi and VHi into single-chain variable fragments (Fv) is prevented by shortened inter-domain linkers. Instead, VH and VL are expected to interact with VLi and VHi, respectively, thus making a diabody whose binding to CD3ε on the T-cells is impaired. METHODS We analyzed the structure of an epidermal growth factor receptor (EGFR) COBRA in solution using negative stain electron microscopy (EM) and small-angle X-ray scattering (SAXS). RESULTS We found that this EGFR COBRA forms stable monomers with a very dynamic interdomain arrangement. At most, only five domains at a time appeared ordered, and only one VH-VL pair was found in the Fv orientation. Nonenzymatic posttranslational modifications suggest that the CDR3 loops in the VL-VHi pair are exposed but are buried in the VH-VLi pair. The MMP9 cleavage rate of the prodrug when bound to recombinant EGFR or HSA is not affected, indicating positioning of the MMP9-cleavable linker away from the EGFR and HSA binding sites. CONCLUSION Here, we propose a model for EGFR COBRA where VH and VLi form an Fv, and VL and VHi do not, possibly interacting with other Ig domains. SAXS and MMP9 cleavage analyses suggest that all COBRA molecules tested have a similar structural architecture.
Collapse
Affiliation(s)
- Jessica Krakow
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Michal Hammel
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Ying Zhu
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Brian J Hillier
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Bryce Paolella
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Austin Desmarais
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Rusty Wall
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Tseng-Hui T Chen
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Rex Pei
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Chulani Karunatilake
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Robert DuBridge
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| | - Maia Vinogradova
- Maverick Therapeutics, Inc., a wholly owned subsidiary of Takeda Pharmaceutical Company Limited, Brisbane, CA, USA
| |
Collapse
|
11
|
Miner JC, Fenimore PW, Fischer WM, McMahon BH, Sanbonmatsu KY, Tung CS. Integrative structural studies of the SARS-CoV-2 spike protein during the fusion process (2022). Curr Res Struct Biol 2022; 4:220-230. [PMID: 35765663 PMCID: PMC9221923 DOI: 10.1016/j.crstbi.2022.06.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Revised: 06/15/2022] [Accepted: 06/17/2022] [Indexed: 11/26/2022] Open
|
12
|
Zeng Y, Howe G, Yi K, Zeng X, Zhang J, Chang YW, Xu M. UNSUPERVISED DOMAIN ALIGNMENT BASED OPEN SET STRUCTURAL RECOGNITION OF MACROMOLECULES CAPTURED BY CRYO-ELECTRON TOMOGRAPHY. PROCEEDINGS. INTERNATIONAL CONFERENCE ON IMAGE PROCESSING 2021; 2021:106-110. [PMID: 35350462 PMCID: PMC8959888 DOI: 10.1109/icip42928.2021.9506205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Cellular cryo-Electron Tomography (cryo-ET) provides three-dimensional views of structural and spatial information of various macromolecules in cells in a near-native state. Subtomogram classification is a key step for recognizing and differentiating these macromolecular structures. In recent years, deep learning methods have been developed for high-throughput subtomogram classification tasks; however, conventional supervised deep learning methods cannot recognize macromolecular structural classes that do not exist in the training data. This imposes a major weakness since most native macromolecular structures in cells are unknown and consequently, cannot be included in the training data. Therefore, open set learning which can recognize unknown macromolecular structures is necessary for boosting the power of automatic subtomogram classification. In this paper, we propose a method called Margin-based Loss for Unsupervised Domain Alignment (MLUDA) for open set recognition problems where only a few categories of interest are shared between cross-domain data. Through extensive experiments, we demonstrate that MLUDA performs well at cross-domain open-set classification on both public datasets and medical imaging datasets. So our method is of practical importance.
Collapse
Affiliation(s)
- Yuchen Zeng
- Computational Biology Department, Carnegie Mellon University, United States
| | - Gregory Howe
- Computational Biology Department, Carnegie Mellon University, United States
| | - Kai Yi
- King Abdullah University of Science and Technology, Saudi Arabia
| | - Xiangrui Zeng
- Computational Biology Department, Carnegie Mellon University, United States
| | - Jing Zhang
- Department of Computer Science, University of California Irvine, United States
| | - Yi-Wei Chang
- Perelman School of Medicine, University of Pennsylvania, United States
| | - Min Xu
- Computational Biology Department, Carnegie Mellon University, United States
| |
Collapse
|
13
|
Woll KA, Van Petegem F. Calcium Release Channels: Structure and Function of IP3 Receptors and Ryanodine Receptors. Physiol Rev 2021; 102:209-268. [PMID: 34280054 DOI: 10.1152/physrev.00033.2020] [Citation(s) in RCA: 79] [Impact Index Per Article: 26.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Ca2+-release channels are giant membrane proteins that control the release of Ca2+ from the endoplasmic and sarcoplasmic reticulum. The two members, ryanodine receptors (RyRs) and inositol-1,4,5-trisphosphate Receptors (IP3Rs), are evolutionarily related and are both activated by cytosolic Ca2+. They share a common architecture, but RyRs have evolved additional modules in the cytosolic region. Their massive size allows for the regulation by tens of proteins and small molecules, which can affect the opening and closing of the channels. In addition to Ca2+, other major triggers include IP3 for the IP3Rs, and depolarization of the plasma membrane for a particular RyR subtype. Their size has made them popular targets for study via electron microscopic methods, with current structures culminating near 3Å. The available structures have provided many new mechanistic insights int the binding of auxiliary proteins and small molecules, how these can regulate channel opening, and the mechanisms of disease-associated mutations. They also help scrutinize previously proposed binding sites, as some of these are now incompatible with the structures. Many questions remain around the structural effects of post-translational modifications, additional binding partners, and the higher-order complexes these channels can make in situ. This review summarizes our current knowledge about the structures of Ca2+-release channels and how this informs on their function.
Collapse
Affiliation(s)
- Kellie A Woll
- Department of Biochemistry and Molecular Biology, University of British Columbia, Vancouver, BC, Canada
| | - Filip Van Petegem
- Department of Biochemistry and Molecular Biology, University of British Columbia, Vancouver, BC, Canada
| |
Collapse
|
14
|
Wang X, Alnabati E, Aderinwale TW, Maddhuri Venkata Subramaniya SR, Terashi G, Kihara D. Detecting protein and DNA/RNA structures in cryo-EM maps of intermediate resolution using deep learning. Nat Commun 2021; 12:2302. [PMID: 33863902 PMCID: PMC8052361 DOI: 10.1038/s41467-021-22577-3] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2020] [Accepted: 03/19/2021] [Indexed: 12/21/2022] Open
Abstract
An increasing number of density maps of macromolecular structures, including proteins and DNA/RNA complexes, have been determined by cryo-electron microscopy (cryo-EM). Although lately maps at a near-atomic resolution are routinely reported, there are still substantial fractions of maps determined at intermediate or low resolutions, where extracting structure information is not trivial. Here, we report a new computational method, Emap2sec+, which identifies DNA or RNA as well as the secondary structures of proteins in cryo-EM maps of 5 to 10 Å resolution. Emap2sec+ employs the deep Residual convolutional neural network. Emap2sec+ assigns structural labels with associated probabilities at each voxel in a cryo-EM map, which will help structure modeling in an EM map. Emap2sec+ showed stable and high assignment accuracy for nucleotides in low resolution maps and improved performance for protein secondary structure assignments than its earlier version when tested on simulated and experimental maps.
Collapse
Affiliation(s)
- Xiao Wang
- Department of Computer Science, Purdue University, West Lafayette, IN, USA
| | - Eman Alnabati
- Department of Computer Science, Purdue University, West Lafayette, IN, USA
| | - Tunde W Aderinwale
- Department of Computer Science, Purdue University, West Lafayette, IN, USA
| | | | - Genki Terashi
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA
| | - Daisuke Kihara
- Department of Computer Science, Purdue University, West Lafayette, IN, USA.
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA.
| |
Collapse
|
15
|
Chiu W, Schmid MF, Pintilie GD, Lawson CL. Evolution of standardization and dissemination of cryo-EM structures and data jointly by the community, PDB, and EMDB. J Biol Chem 2021; 296:100560. [PMID: 33744287 PMCID: PMC8050867 DOI: 10.1016/j.jbc.2021.100560] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 02/08/2021] [Accepted: 03/16/2021] [Indexed: 01/04/2023] Open
Abstract
Cryogenic electron microscopy (cryo-EM) methods began to be used in the mid-1970s to study thin and periodic arrays of proteins. Following a half-century of development in cryo-specimen preparation, instrumentation, data collection, data processing, and modeling software, cryo-EM has become a routine method for solving structures from large biological assemblies to small biomolecules at near to true atomic resolution. This review explores the critical roles played by the Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in partnership with the community to develop the necessary infrastructure to archive cryo-EM maps and associated models. Public access to cryo-EM structure data has in turn facilitated better understanding of structure–function relationships and advancement of image processing and modeling tool development. The partnership between the global cryo-EM community and PDB and EMDB leadership has synergistically shaped the standards for metadata, one-stop deposition of maps and models, and validation metrics to assess the quality of cryo-EM structures. The advent of cryo-electron tomography (cryo-ET) for in situ molecular cell structures at a broad resolution range and their correlations with other imaging data introduce new data archival challenges in terms of data size and complexity in the years to come.
Collapse
Affiliation(s)
- Wah Chiu
- Department of Bioengineering, Stanford University, Stanford, California, USA; Division of CryoEM and Bioimaging, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California, USA.
| | - Michael F Schmid
- Division of CryoEM and Bioimaging, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California, USA
| | - Grigore D Pintilie
- Department of Bioengineering, Stanford University, Stanford, California, USA
| | - Catherine L Lawson
- Institute for Quantitative Biomedicine and Research Collaboratory for Structural Bioinformatics, Rutgers, The State University of New Jersey, Piscataway, New Jersey, USA
| |
Collapse
|
16
|
Singla J, White KL, Stevens RC, Alber F. Assessment of scoring functions to rank the quality of 3D subtomogram clusters from cryo-electron tomography. J Struct Biol 2021; 213:107727. [PMID: 33753204 DOI: 10.1016/j.jsb.2021.107727] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Revised: 03/12/2021] [Accepted: 03/17/2021] [Indexed: 11/17/2022]
Abstract
Cryo-electron tomography provides the opportunity for unsupervised discovery of endogenous complexes in situ. This process usually requires particle picking, clustering and alignment of subtomograms to produce an average structure of the complex. When applied to heterogeneous samples, template-free clustering and alignment of subtomograms can potentially lead to the discovery of structures for unknown endogenous complexes. However, such methods require scoring functions to measure and accurately rank the quality of aligned subtomogram clusters, which can be compromised by contaminations from misclassified complexes and alignment errors. Here, we provide the first study to assess the effectiveness of more than 15 scoring functions for evaluating the quality of subtomogram clusters, which differ in the amount of structural misalignments and contaminations due to misclassified complexes. We assessed both experimental and simulated subtomograms as ground truth data sets. Our analysis showed that the robustness of scoring functions varies largely. Most scores were sensitive to the signal-to-noise ratio of subtomograms and often required Gaussian filtering as preprocessing for improved performance. Two scoring functions, Spectral SNR-based Fourier Shell Correlation and Pearson Correlation in the Fourier domain with missing wedge correction, showed a robust ranking of subtomogram clusters without any preprocessing and irrespective of SNR levels of subtomograms. Of these two scoring functions, Spectral SNR-based Fourier Shell Correlation was fastest to compute and is a better choice for handling large numbers of subtomograms. Our results provide a guidance for choosing an accurate scoring function for template-free approaches to detect complexes from heterogeneous samples.
Collapse
Affiliation(s)
- Jitin Singla
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, 520 Boyer Hall, Los Angeles, CA 90095, USA; Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, 1050 Childs Way, Los Angeles, CA 90089, USA; Department of Biological Sciences, Bridge Institute, Michelson Center for Convergent Bioscience, University of Southern California, Los Angeles, CA 90089, USA
| | - Kate L White
- Department of Biological Sciences, Bridge Institute, Michelson Center for Convergent Bioscience, University of Southern California, Los Angeles, CA 90089, USA
| | - Raymond C Stevens
- Department of Biological Sciences, Bridge Institute, Michelson Center for Convergent Bioscience, University of Southern California, Los Angeles, CA 90089, USA
| | - Frank Alber
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, 520 Boyer Hall, Los Angeles, CA 90095, USA; Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, 1050 Childs Way, Los Angeles, CA 90089, USA.
| |
Collapse
|
17
|
Kulik M, Mori T, Sugita Y. Multi-Scale Flexible Fitting of Proteins to Cryo-EM Density Maps at Medium Resolution. Front Mol Biosci 2021; 8:631854. [PMID: 33842541 PMCID: PMC8025875 DOI: 10.3389/fmolb.2021.631854] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Accepted: 01/26/2021] [Indexed: 11/13/2022] Open
Abstract
Structure determination using cryo-electron microscopy (cryo-EM) medium-resolution density maps is often facilitated by flexible fitting. Avoiding overfitting, adjusting force constants driving the structure to the density map, and emulating complex conformational transitions are major concerns in the fitting. To address them, we develop a new method based on a three-step multi-scale protocol. First, flexible fitting molecular dynamics (MD) simulations with coarse-grained structure-based force field and replica-exchange scheme between different force constants replicas are performed. Second, fitted Cα atom positions guide the all-atom structure in targeted MD. Finally, the all-atom flexible fitting refinement in implicit solvent adjusts the positions of the side chains in the density map. Final models obtained via the multi-scale protocol are significantly better resolved and more reliable in comparison with long all-atom flexible fitting simulations. The protocol is useful for multi-domain systems with intricate structural transitions as it preserves the secondary structure of single domains.
Collapse
Affiliation(s)
- Marta Kulik
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Japan.,RIKEN Center for Computational Science, Kobe, Japan.,RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| |
Collapse
|
18
|
Woods DC, Rodríguez-Ropero F, Wereszczynski J. The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome Array. J Mol Biol 2021; 433:166902. [PMID: 33667509 DOI: 10.1016/j.jmb.2021.166902] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Revised: 02/15/2021] [Accepted: 02/20/2021] [Indexed: 02/08/2023]
Abstract
Linker histones bind to nucleosomes and modify chromatin structure and dynamics as a means of epigenetic regulation. Biophysical studies have shown that chromatin fibers can adopt a plethora of conformations with varying levels of compaction. Linker histone condensation, and its specific binding disposition, has been associated with directly tuning this ensemble of states. However, the atomistic dynamics and quantification of this mechanism remains poorly understood. Here, we present molecular dynamics simulations of octa-nucleosome arrays, based on a cryo-EM structure of the 30-nm chromatin fiber, with and without the globular domains of the H1 linker histone to determine how they influence fiber structures and dynamics. Results show that when bound, linker histones inhibit DNA flexibility and stabilize repeating tetra-nucleosomal units, giving rise to increased chromatin compaction. Furthermore, upon the removal of H1, there is a significant destabilization of this compact structure as the fiber adopts less strained and untwisted states. Interestingly, linker DNA sampling in the octa-nucleosome is exaggerated compared to its mono-nucleosome counterparts, suggesting that chromatin architecture plays a significant role in DNA strain even in the absence of linker histones. Moreover, H1-bound states are shown to have increased stiffness within tetra-nucleosomes, but not between them. This increased stiffness leads to stronger long-range correlations within the fiber, which may result in the propagation of epigenetic signals over longer spatial ranges. These simulations highlight the effects of linker histone binding on the internal dynamics and global structure of poly-nucleosome arrays, while providing physical insight into a mechanism of chromatin compaction.
Collapse
Affiliation(s)
- Dustin C Woods
- Department of Chemistry and the Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL 60616, United States
| | - Francisco Rodríguez-Ropero
- Department of Physics and the Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL 60616, United States
| | - Jeff Wereszczynski
- Department of Physics and the Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL 60616, United States.
| |
Collapse
|
19
|
Du X, Wang H, Zhu Z, Zeng X, Chang YW, Zhang J, Xing E, Xu M. Active Learning to Classify Macromolecular Structures in situ for Less Supervision in Cryo-Electron Tomography. Bioinformatics 2021; 37:2340-2346. [PMID: 33620460 DOI: 10.1093/bioinformatics/btab123] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2020] [Revised: 01/14/2021] [Accepted: 02/22/2021] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Cryo-Electron Tomography (cryo-ET) is a 3D bioimaging tool that visualizes the structural and spatial organization of macromolecules at a near-native state in single cells, which has broad applications in life science. However, the systematic structural recognition and recovery of macromolecules captured by cryo-ET are difficult due to high structural complexity and imaging limits. Deep learning based subtomogram classification have played critical roles for such tasks. As supervised approaches, however, their performance relies on sufficient and laborious annotation on a large training dataset. RESULTS To alleviate this major labeling burden, we proposed a Hybrid Active Learning (HAL) framework for querying subtomograms for labelling from a large unlabeled subtomogram pool. Firstly, HAL adopts uncertainty sampling to select the subtomograms that have the most uncertain predictions. This strategy enforces the model to be aware of the inductive bias during classification and subtomogram selection, which satisfies the discriminativeness principle in AL literature. Moreover, to mitigate the sampling bias caused by such strategy, a discriminator is introduced to judge if a certain subtomogram is labeled or unlabeled and subsequently the model queries the subtomogram that have higher probabilities to be unlabeled. Such query strategy encourages to match the data distribution between the labeled and unlabeled subtomogram samples, which essentially encodes the representativeness criterion into the subtomogram selection process. Additionally, HAL introduces a subset sampling strategy to improve the diversity of the query set, so that the information overlap is decreased between the queried batches and the algorithmic efficiency is improved. Our experiments on subtomogram classification tasks using both simulated and real data demonstrate that we can achieve comparable testing performance (on average only 3% accuracy drop) by using less than 30% of the labeled subtomograms, which shows a very promising result for subtomogram classification task with limited labeling resources. AVAILABILITY https://github.com/xulabs/aitom.
Collapse
Affiliation(s)
- Xuefeng Du
- Department of Computer Science, University of Wisconsin-Madison, Madison, 53706, USA
| | - Haohan Wang
- Language Technologies Institute, Carnegie Mellon University, Pittsburgh, 15213, USA
| | - Zhenxi Zhu
- Department of Computer Science, Beijing University of Posts and Telecommunications, 100876, China
| | - Xiangrui Zeng
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, 15213, USA
| | - Yi-Wei Chang
- Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, 19104, USA
| | - Jing Zhang
- Department of Computer Science, University of California - Irvine, Irvine, 92697, USA
| | - Eric Xing
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, 15213, USA
| | - Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, 15213, USA
| |
Collapse
|
20
|
Hodge CD, Rosenberg DJ, Wilamowski M, Joachimiak A, Hura GL, Hammel M. Rigid monoclonal antibodies improve detection of SARS-CoV-2 nucleocapsid protein. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2021.01.13.426597. [PMID: 33469584 PMCID: PMC7814821 DOI: 10.1101/2021.01.13.426597] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Monoclonal antibodies (mAbs) are the basis of treatments and diagnostics for pathogens and other biological phenomena. We conducted a structural characterization of mAbs against the N-terminal domain of nucleocapsid protein (NP NTD ) from SARS-CoV-2 using small angle X-ray scattering (SAXS). Our solution-based results distinguished the mAbs' flexibility and how this flexibility impacts the assembly of multiple mAbs on an antigen. By pairing two mAbs that bind different epitopes on the NP NTD , we show that flexible mAbs form a closed sandwich-like complex. With rigid mAbs, a juxtaposition of the Fabs is prevented, enforcing a linear arrangement of the mAb pair, which facilitates further mAb polymerization. In a modified sandwich ELISA, we show the rigid mAb-pairings with linear polymerization led to increased NP NTD detection sensitivity. These enhancements can expedite the development of more sensitive and selective antigen-detecting point-of-care lateral flow devices (LFA), key for early diagnosis and epidemiological studies of SARS-CoV-2 and other pathogens.
Collapse
Affiliation(s)
- Curtis D Hodge
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Daniel J Rosenberg
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Mateusz Wilamowski
- Center for Structural Genomics of Infectious Diseases, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, USA
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA
| | - Andrzej Joachimiak
- Center for Structural Genomics of Infectious Diseases, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, USA
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA
- Structural Biology Center, X-ray Science Division, Argonne National Laboratory, Argonne, IL, USA
| | - Greg L Hura
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Michal Hammel
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
21
|
Hodge CD, Rosenberg DJ, Grob P, Wilamowski M, Joachimiak A, Hura GL, Hammel M. Rigid monoclonal antibodies improve detection of SARS-CoV-2 nucleocapsid protein. MAbs 2021; 13:1905978. [PMID: 33843452 PMCID: PMC8043170 DOI: 10.1080/19420862.2021.1905978] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Revised: 03/10/2021] [Accepted: 03/16/2021] [Indexed: 12/13/2022] Open
Abstract
Monoclonal antibodies (mAbs) are the basis of treatments and diagnostics for pathogens and other biological phenomena. We conducted a structural characterization of mAbs against the N-terminal domain of nucleocapsid protein (NPNTD) from SARS-CoV-2 using small-angle X-ray scattering and transmission electron microscopy. Our solution-based results distinguished the mAbs' flexibility and how this flexibility affects the assembly of multiple mAbs on an antigen. By pairing two mAbs that bind different epitopes on the NPNTD, we show that flexible mAbs form a closed sandwich-like complex. With rigid mAbs, a juxtaposition of the antigen-binding fragments is prevented, enforcing a linear arrangement of the mAb pair, which facilitates further mAb polymerization. In a modified sandwich enzyme-linked immunosorbent assay, we show that rigid mAb-pairings with linear polymerization led to increased NPNTD detection sensitivity. These enhancements can expedite the development of more sensitive and selective antigen-detecting point-of-care lateral flow devices, which are critical for early diagnosis and epidemiological studies of SARS-CoV-2 and other pathogens.
Collapse
Affiliation(s)
- Curtis D. Hodge
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Daniel. J. Rosenberg
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Graduate Group in Biophysics, University of California, Berkeley, CA, USA
| | - Patricia Grob
- Howard Hughes Medical Institute, UC Berkeley, Berkeley, CA, USA
- California Institute for Quantitative Biosciences, University of California, Berkeley, CA, USA
| | - Mateusz Wilamowski
- Center for Structural Genomics of Infectious Diseases, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, USA
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA
| | - Andrzej Joachimiak
- Center for Structural Genomics of Infectious Diseases, Consortium for Advanced Science and Engineering, University of Chicago, Chicago, IL, USA
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, IL, USA
- Structural Biology Center, X-ray Science Division, Argonne National Laboratory, Argonne, IL, USA
| | - Greg L. Hura
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Chemistry and Biochemistry Department, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Michal Hammel
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
22
|
Seffernick JT, Lindert S. Hybrid methods for combined experimental and computational determination of protein structure. J Chem Phys 2020; 153:240901. [PMID: 33380110 PMCID: PMC7773420 DOI: 10.1063/5.0026025] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Accepted: 11/10/2020] [Indexed: 02/04/2023] Open
Abstract
Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.
Collapse
Affiliation(s)
- Justin T. Seffernick
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, USA
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
23
|
Jung T, Shin B, Tamo G, Kim H, Vijayvargia R, Leitner A, Marcaida MJ, Astorga-Wells J, Jung R, Aebersold R, Peraro MD, Hebert H, Seong IS, Song JJ. The Polyglutamine Expansion at the N-Terminal of Huntingtin Protein Modulates the Dynamic Configuration and Phosphorylation of the C-Terminal HEAT Domain. Structure 2020; 28:1035-1050.e8. [PMID: 32668197 PMCID: PMC11059206 DOI: 10.1016/j.str.2020.06.008] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Revised: 04/01/2020] [Accepted: 06/23/2020] [Indexed: 11/15/2022]
Abstract
The polyQ expansion in huntingtin protein (HTT) is the prime cause of Huntington's disease (HD). The recent cryoelectron microscopy (cryo-EM) structure of HTT-HAP40 complex provided the structural information on its HEAT-repeat domains. Here, we present analyses of the impact of polyQ length on the structure and function of HTT via an integrative structural and biochemical approach. The cryo-EM analysis of normal (Q23) and disease (Q78) type HTTs shows that the structures of apo HTTs significantly differ from the structure of HTT in a HAP40 complex and that the polyQ expansion induces global structural changes in the relative movements among the HTT domains. In addition, we show that the polyQ expansion alters the phosphorylation pattern across HTT and that Ser2116 phosphorylation in turn affects the global structure and function of HTT. These results provide a molecular basis for the effect of the polyQ segment on HTT structure and activity, which may be important for HTT pathology.
Collapse
Affiliation(s)
- Taeyang Jung
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), KI for the BioCentury, Daejeon 34141, Korea; School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, 141 52 Huddinge, Sweden; Department of Biosciences and Nutrition, Karolinska Institutet, 141 83 Huddinge, Sweden
| | - Baehyun Shin
- Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02114, USA
| | - Giorgio Tamo
- Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Hyeongju Kim
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), KI for the BioCentury, Daejeon 34141, Korea
| | - Ravi Vijayvargia
- Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02114, USA
| | - Alexander Leitner
- Department of Biology, Institute of Molecular Systems Biology, ETH Zürich, 8093 Zürich, Switzerland
| | - Maria J Marcaida
- Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Juan Astorga-Wells
- Department of Medical Biochemistry & Biophysics, Karolinska Institutet, 171 65 Solna, Sweden; HDxperts AB, 183 48 Täby, Sweden
| | - Roy Jung
- Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02114, USA
| | - Ruedi Aebersold
- Department of Biology, Institute of Molecular Systems Biology, ETH Zürich, 8093 Zürich, Switzerland; Faculty of Science, University of Zürich, Zürich, Switzerland
| | - Matteo Dal Peraro
- Institute of Bioengineering, School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Hans Hebert
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, 141 52 Huddinge, Sweden; Department of Biosciences and Nutrition, Karolinska Institutet, 141 83 Huddinge, Sweden.
| | - Ihn Sik Seong
- Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA 02114, USA; Department of Neurology, Harvard Medical School, Boston, MA 02114, USA.
| | - Ji-Joon Song
- Department of Biological Sciences, Korea Advanced Institute of Science and Technology (KAIST), KI for the BioCentury, Daejeon 34141, Korea.
| |
Collapse
|
24
|
Woods DC, Wereszczynski J. Elucidating the influence of linker histone variants on chromatosome dynamics and energetics. Nucleic Acids Res 2020; 48:3591-3604. [PMID: 32128577 PMCID: PMC7144933 DOI: 10.1093/nar/gkaa121] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Revised: 02/12/2020] [Accepted: 02/14/2020] [Indexed: 12/23/2022] Open
Abstract
Linker histones are epigenetic regulators that bind to nucleosomes and alter chromatin structures and dynamics. Biophysical studies have revealed two binding modes in the linker histone/nucleosome complex, the chromatosome, where the linker histone is either centered on or askew from the dyad axis. Each has been posited to have distinct effects on chromatin, however the molecular and thermodynamic mechanisms that drive them and their dependence on linker histone compositions remain poorly understood. We present molecular dynamics simulations of chromatosomes with the globular domain of two linker histone variants, generic H1 (genGH1) and H1.0 (GH1.0), to determine how their differences influence chromatosome structures, energetics and dynamics. Results show that both unbound linker histones adopt a single compact conformation. Upon binding, DNA flexibility is reduced, resulting in increased chromatosome compaction. While both variants enthalpically favor on-dyad binding, energetic benefits are significantly higher for GH1.0, suggesting that GH1.0 is more capable than genGH1 of overcoming the large entropic reduction required for on-dyad binding which helps rationalize experiments that have consistently demonstrated GH1.0 in on-dyad states but that show genGH1 in both locations. These simulations highlight the thermodynamic basis for different linker histone binding motifs, and details their physical and chemical effects on chromatosomes.
Collapse
Affiliation(s)
- Dustin C Woods
- Department of Chemistry and the Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL 60616, USA
| | - Jeff Wereszczynski
- Department of Physics and the Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, IL 60616, USA
| |
Collapse
|
25
|
Shape-preserving elastic solid models of macromolecules. PLoS Comput Biol 2020; 16:e1007855. [PMID: 32407309 PMCID: PMC7297265 DOI: 10.1371/journal.pcbi.1007855] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Revised: 06/09/2020] [Accepted: 04/07/2020] [Indexed: 02/07/2023] Open
Abstract
Mass-spring models have been a standard approach in molecular modeling for the last few decades, such as elastic network models (ENMs) that are widely used for normal mode analysis. In this work, we present a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as ENMs in producing the equilibrium dynamics and moreover, offers some significant new features that may greatly benefit the research community. ESM is different from ENM in that it treats macromolecules as elastic solids. Our particular version of ESM presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape and thus makes normal mode computations and visualization of extremely large complexes more manageable. Secondly, as a solid model, ESM’s close link to finite element analysis renders it ideally suited for studying mechanical responses of macromolecules under external force. Lastly, we show that ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy. The complete MATLAB code of αESM is provided. Mass-spring models have been a standard approach in classical molecular modeling where atoms are modeled as spheres with a mass and their interactions modeled as springs. The models have been extremely successful. Thinking ahead, however, as molecular systems of our interest grow more quickly in size or dimension than what our computation resources can keep up with, some adjustments in methodology are timely. This work presents a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as mass-spring models in producing the equilibrium dynamics and moreover, offers some unique features that make it suitable for much larger systems. ESM is different from ENMs in that it treats macromolecules as elastic solids. Our particular version of ESM model presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape. ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy.
Collapse
|
26
|
Dodd T, Yan C, Ivanov I. Simulation-Based Methods for Model Building and Refinement in Cryoelectron Microscopy. J Chem Inf Model 2020; 60:2470-2483. [PMID: 32202798 DOI: 10.1021/acs.jcim.0c00087] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Thomas Dodd
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States.,Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302, United States
| | - Chunli Yan
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States.,Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302, United States
| | - Ivaylo Ivanov
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States.,Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30302, United States
| |
Collapse
|
27
|
Akbar S, Mozumder S, Sengupta J. Retrospect and Prospect of Single Particle Cryo-Electron Microscopy: The Class of Integral Membrane Proteins as an Example. J Chem Inf Model 2020; 60:2448-2457. [DOI: 10.1021/acs.jcim.9b01015] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shirin Akbar
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
| | - Sukanya Mozumder
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Jayati Sengupta
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
28
|
Martínez M, Jiménez-Moreno A, Maluenda D, Ramírez-Aportela E, Melero R, Cuervo A, Conesa P, Del Caño L, Fonseca YC, Sánchez-García R, Strelak D, Conesa JJ, Fernández-Giménez E, de Isidro F, Sorzano COS, Carazo JM, Marabini R. Integration of Cryo-EM Model Building Software in Scipion. J Chem Inf Model 2020; 60:2533-2540. [PMID: 31994878 DOI: 10.1021/acs.jcim.9b01032] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Advances in cryo-electron microscopy (cryo-EM) have made it possible to obtain structures of large biological macromolecules at near-atomic resolution. This "resolution revolution" has encouraged the use and development of modeling tools able to produce high-quality atomic models from cryo-EM density maps. Unfortunately, many practical problems appear when combining different packages in the same processing workflow, which make difficult the use of these tools by non-experts and, therefore, reduce their utility. We present here a major extension of the image processing framework Scipion that provides inter-package integration in the model building area and full tracking of the complete workflow, from image processing to structure validation.
Collapse
Affiliation(s)
- M Martínez
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | | | - D Maluenda
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | | | - R Melero
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | - A Cuervo
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | - P Conesa
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | - L Del Caño
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | | | | | - D Strelak
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain.,Institute of Computer Science, Masaryk University, Botanická 68a, 60200 Brno, Czech Republic
| | - J J Conesa
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | | | | | | | - J M Carazo
- CNB-CSIC, C/Darwin 3, 28049 Madrid, Spain
| | - R Marabini
- Escuela Politécnica, Universidad Autónoma de Madrid, C/Francisco Tomás y Valiente 11, 28049 Madrid, Spain
| |
Collapse
|
29
|
Alnabati E, Kihara D. Advances in Structure Modeling Methods for Cryo-Electron Microscopy Maps. Molecules 2019; 25:molecules25010082. [PMID: 31878333 PMCID: PMC6982917 DOI: 10.3390/molecules25010082] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Revised: 12/20/2019] [Accepted: 12/20/2019] [Indexed: 01/16/2023] Open
Abstract
Cryo-electron microscopy (cryo-EM) has now become a widely used technique for structure determination of macromolecular complexes. For modeling molecular structures from density maps of different resolutions, many algorithms have been developed. These algorithms can be categorized into rigid fitting, flexible fitting, and de novo modeling methods. It is also observed that machine learning (ML) techniques have been increasingly applied following the rapid progress of the ML field. Here, we review these different categories of macromolecule structure modeling methods and discuss their advances over time.
Collapse
Affiliation(s)
- Eman Alnabati
- Department of Computer Science, Purdue University, West Lafayette, IN 47907, USA
| | - Daisuke Kihara
- Department of Computer Science, Purdue University, West Lafayette, IN 47907, USA
- Department of Biological Sciences, Purdue University, West Lafayette, IN 47907, USA
- Correspondence:
| |
Collapse
|
30
|
Zhang K, Li S, Kappel K, Pintilie G, Su Z, Mou TC, Schmid MF, Das R, Chiu W. Cryo-EM structure of a 40 kDa SAM-IV riboswitch RNA at 3.7 Å resolution. Nat Commun 2019; 10:5511. [PMID: 31796736 PMCID: PMC6890682 DOI: 10.1038/s41467-019-13494-7] [Citation(s) in RCA: 76] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 11/03/2019] [Indexed: 01/17/2023] Open
Abstract
Specimens below 50 kDa have generally been considered too small to be analyzed by single-particle cryo-electron microscopy (cryo-EM). The high flexibility of pure RNAs makes it difficult to obtain high-resolution structures by cryo-EM. In bacteria, riboswitches regulate sulfur metabolism through binding to the S-adenosylmethionine (SAM) ligand and offer compelling targets for new antibiotics. SAM-I, SAM-I/IV, and SAM-IV are the three most commonly found SAM riboswitches, but the structure of SAM-IV is still unknown. Here, we report the structures of apo and SAM-bound SAM-IV riboswitches (119-nt, ~40 kDa) to 3.7 Å and 4.1 Å resolution, respectively, using cryo-EM. The structures illustrate homologies in the ligand-binding core but distinct peripheral tertiary contacts in SAM-IV compared to SAM-I and SAM-I/IV. Our results demonstrate the feasibility of resolving small RNAs with enough detail to enable detection of their ligand-binding pockets and suggest that cryo-EM could play a role in structure-assisted drug design for RNA.
Collapse
Affiliation(s)
- Kaiming Zhang
- Department of Bioengineering, and James H. Clark Center, Stanford University, Stanford, CA, 94305, USA
| | - Shanshan Li
- Department of Bioengineering, and James H. Clark Center, Stanford University, Stanford, CA, 94305, USA
| | - Kalli Kappel
- Biophysics Program, Stanford University, Stanford, CA, 94305, USA
| | - Grigore Pintilie
- Department of Bioengineering, and James H. Clark Center, Stanford University, Stanford, CA, 94305, USA
| | - Zhaoming Su
- Department of Bioengineering, and James H. Clark Center, Stanford University, Stanford, CA, 94305, USA
| | - Tung-Chung Mou
- Center for Biomolecular Structure and Dynamics, University of Montana, Missoula, MT, 59812, USA
| | - Michael F Schmid
- Division of CryoEM and Bioimaging, SSRL, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, CA, 94025, USA
| | - Rhiju Das
- Department of Biochemistry, Stanford University, Stanford, CA, 94305, USA.
- Department of Physics, Stanford University, Stanford, CA, 94305, USA.
| | - Wah Chiu
- Department of Bioengineering, and James H. Clark Center, Stanford University, Stanford, CA, 94305, USA.
- Division of CryoEM and Bioimaging, SSRL, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, CA, 94025, USA.
| |
Collapse
|
31
|
Lü Y, Zeng X, Zhao X, Li S, Li H, Gao X, Xu M. Fine-grained alignment of cryo-electron subtomograms based on MPI parallel optimization. BMC Bioinformatics 2019; 20:443. [PMID: 31455212 PMCID: PMC6712796 DOI: 10.1186/s12859-019-3003-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2018] [Accepted: 07/19/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Cryo-electron tomography (Cryo-ET) is an imaging technique used to generate three-dimensional structures of cellular macromolecule complexes in their native environment. Due to developing cryo-electron microscopy technology, the image quality of three-dimensional reconstruction of cryo-electron tomography has greatly improved. However, cryo-ET images are characterized by low resolution, partial data loss and low signal-to-noise ratio (SNR). In order to tackle these challenges and improve resolution, a large number of subtomograms containing the same structure needs to be aligned and averaged. Existing methods for refining and aligning subtomograms are still highly time-consuming, requiring many computationally intensive processing steps (i.e. the rotations and translations of subtomograms in three-dimensional space). RESULTS In this article, we propose a Stochastic Average Gradient (SAG) fine-grained alignment method for optimizing the sum of dissimilarity measure in real space. We introduce a Message Passing Interface (MPI) parallel programming model in order to explore further speedup. CONCLUSIONS We compare our stochastic average gradient fine-grained alignment algorithm with two baseline methods, high-precision alignment and fast alignment. Our SAG fine-grained alignment algorithm is much faster than the two baseline methods. Results on simulated data of GroEL from the Protein Data Bank (PDB ID:1KP8) showed that our parallel SAG-based fine-grained alignment method could achieve close-to-optimal rigid transformations with higher precision than both high-precision alignment and fast alignment at a low SNR (SNR=0.003) with tilt angle range ±60∘ or ±40∘. For the experimental subtomograms data structures of GroEL and GroEL/GroES complexes, our parallel SAG-based fine-grained alignment can achieve higher precision and fewer iterations to converge than the two baseline methods.
Collapse
Affiliation(s)
- Yongchun Lü
- University of Chinese Academy of Sciences, Beijing, China
- Institute of Computing Technology of the Chinese Academy of Sciences, Beijing, China
- Key Laboratory of Intelligent Information Processing, CAS, Beijing, China
| | - Xiangrui Zeng
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, USA
| | - Xiaofang Zhao
- University of Chinese Academy of Sciences, Beijing, China
- Institute of Computing Technology of the Chinese Academy of Sciences, Beijing, China
| | - Shirui Li
- Institute of Computing Technology of the Chinese Academy of Sciences, Beijing, China
- Key Laboratory of Intelligent Information Processing, CAS, Beijing, China
| | - Hua Li
- University of Chinese Academy of Sciences, Beijing, China
- Institute of Computing Technology of the Chinese Academy of Sciences, Beijing, China
- Key Laboratory of Intelligent Information Processing, CAS, Beijing, China
| | - Xin Gao
- King Abdullah University of Science and Technology (KAUST), Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, Thuwal, Saudi Arabia
| | - Min Xu
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, USA
| |
Collapse
|
32
|
Zhao Y, Zeng X, Guo Q, Xu M. An integration of fast alignment and maximum-likelihood methods for electron subtomogram averaging and classification. Bioinformatics 2019; 34:i227-i236. [PMID: 29949977 PMCID: PMC6022576 DOI: 10.1093/bioinformatics/bty267] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Motivation Cellular Electron CryoTomography (CECT) is an emerging 3D imaging technique that visualizes subcellular organization of single cells at sub-molecular resolution and in near-native state. CECT captures large numbers of macromolecular complexes of highly diverse structures and abundances. However, the structural complexity and imaging limits complicate the systematic de novo structural recovery and recognition of these macromolecular complexes. Efficient and accurate reference-free subtomogram averaging and classification represent the most critical tasks for such analysis. Existing subtomogram alignment based methods are prone to the missing wedge effects and low signal-to-noise ratio (SNR). Moreover, existing maximum-likelihood based methods rely on integration operations, which are in principle computationally infeasible for accurate calculation. Results Built on existing works, we propose an integrated method, Fast Alignment Maximum Likelihood method (FAML), which uses fast subtomogram alignment to sample sub-optimal rigid transformations. The transformations are then used to approximate integrals for maximum-likelihood update of subtomogram averages through expectation–maximization algorithm. Our tests on simulated and experimental subtomograms showed that, compared to our previously developed fast alignment method (FA), FAML is significantly more robust to noise and missing wedge effects with moderate increases of computation cost. Besides, FAML performs well with significantly fewer input subtomograms when the FA method fails. Therefore, FAML can serve as a key component for improved construction of initial structural models from macromolecules captured by CECT. Availability and implementation http://www.cs.cmu.edu/mxu1
Collapse
Affiliation(s)
- Yixiu Zhao
- Computational Biology and Computer Science Departments, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xiangrui Zeng
- Computational Biology and Computer Science Departments, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Qiang Guo
- Department of Molecular Structural Biology, Max Planck Institute of Biochemistry, Martinsried, Germany
| | - Min Xu
- Computational Biology and Computer Science Departments, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| |
Collapse
|
33
|
Lin R, Zeng X, Kitani K, Xu M. Adversarial domain adaptation for cross data source macromolecule in situ structural classification in cellular electron cryo-tomograms. Bioinformatics 2019; 35:i260-i268. [PMID: 31510673 PMCID: PMC6612867 DOI: 10.1093/bioinformatics/btz364] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
MOTIVATION Since 2017, an increasing amount of attention has been paid to the supervised deep learning-based macromolecule in situ structural classification (i.e. subtomogram classification) in cellular electron cryo-tomography (CECT) due to the substantially higher scalability of deep learning. However, the success of such supervised approach relies heavily on the availability of large amounts of labeled training data. For CECT, creating valid training data from the same data source as prediction data is usually laborious and computationally intensive. It would be beneficial to have training data from a separate data source where the annotation is readily available or can be performed in a high-throughput fashion. However, the cross data source prediction is often biased due to the different image intensity distributions (a.k.a. domain shift). RESULTS We adapt a deep learning-based adversarial domain adaptation (3D-ADA) method to timely address the domain shift problem in CECT data analysis. 3D-ADA first uses a source domain feature extractor to extract discriminative features from the training data as the input to a classifier. Then it adversarially trains a target domain feature extractor to reduce the distribution differences of the extracted features between training and prediction data. As a result, the same classifier can be directly applied to the prediction data. We tested 3D-ADA on both experimental and realistically simulated subtomogram datasets under different imaging conditions. 3D-ADA stably improved the cross data source prediction, as well as outperformed two popular domain adaptation methods. Furthermore, we demonstrate that 3D-ADA can improve cross data source recovery of novel macromolecular structures. AVAILABILITY AND IMPLEMENTATION https://github.com/xulabs/projects. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Ruogu Lin
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xiangrui Zeng
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Kris Kitani
- Robotics Institute, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA, USA
| |
Collapse
|
34
|
Zhang Y, Xia K, Cao Z, Gräter F, Xia F. A new method for the construction of coarse-grained models of large biomolecules from low-resolution cryo-electron microscopy data. Phys Chem Chem Phys 2019; 21:9720-9727. [PMID: 31025999 DOI: 10.1039/c9cp01370a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The rapid development of cryo-electron microscopy (cryo-EM) has led to the generation of significant low-resolution electron density data of biomolecules. However, the atomistic details of huge biomolecules usually cannot be obtained because it is very difficult to construct all-atom models for MD simulations. Thus, it is still a challenge to make use of the rich low-resolution cryo-EM data for computer simulation and functional study. In this study, we proposed a new method called Convolutional and K-means Coarse-Graining (CK-CG) for the efficient coarse-graining of large biological systems. Using the CK-CG method, we could directly map the cryo-EM data into coarse-grained (CG) beads. Furthermore, the CG beads were parameterized with an empirical harmonic potential to construct a new CG model. We subjected the CK-CG models of the fibrillar protein assemblies F-actin and collagen to external forces in pulling dynamic simulations to assess their mechanical response. The agreement between the estimated tensile stiffness between CG models and experiments demonstrates the validity of the CK-CG method. Thus, our method provides a practical strategy for the direct construction of a structural model from low-resolution data for biological function studies.
Collapse
Affiliation(s)
- Yuwei Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.
| | | | | | | | | |
Collapse
|
35
|
Xu M, Singla J, Tocheva EI, Chang YW, Stevens RC, Jensen GJ, Alber F. De Novo Structural Pattern Mining in Cellular Electron Cryotomograms. Structure 2019; 27:679-691.e14. [PMID: 30744995 DOI: 10.1016/j.str.2019.01.005] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2018] [Revised: 07/27/2018] [Accepted: 01/14/2019] [Indexed: 11/16/2022]
Abstract
Electron cryotomography enables 3D visualization of cells in a near-native state at molecular resolution. The produced cellular tomograms contain detailed information about a plethora of macromolecular complexes, their structures, abundances, and specific spatial locations in the cell. However, extracting this information in a systematic way is very challenging, and current methods usually rely on individual templates of known structures. Here, we propose a framework called "Multi-Pattern Pursuit" for de novo discovery of different complexes from highly heterogeneous sets of particles extracted from entire cellular tomograms without using information of known structures. These initially detected structures can then serve as input for more targeted refinement efforts. Our tests on simulated and experimental tomograms show that our automated method is a promising tool for supporting large-scale template-free visual proteomics analysis.
Collapse
Affiliation(s)
- Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA.
| | - Jitin Singla
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology, and Molecular Genetics, University of California, Los Angeles, CA 90095, USA; Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, CA 90089, USA
| | - Elitza I Tocheva
- Department of Microbiology and Immunology, Life Sciences Institute, The University of British Columbia, Vancouver, BC V6T 1Z3, Canada
| | - Yi-Wei Chang
- Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Raymond C Stevens
- Department of Biological Sciences and Department of Chemistry, Bridge Institute, University of Southern California, Los Angeles, CA 90089, USA
| | - Grant J Jensen
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA 91125, USA; Howard Hughes Medical Institute, Pasadena, CA 91125, USA
| | - Frank Alber
- Institute for Quantitative and Computational Biosciences, Department of Microbiology, Immunology, and Molecular Genetics, University of California, Los Angeles, CA 90095, USA; Quantitative and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, CA 90089, USA.
| |
Collapse
|
36
|
Amado E, Muth G, Arechaga I, Cabezón E. The FtsK-like motor TraB is a DNA-dependent ATPase that forms higher-order assemblies. J Biol Chem 2019; 294:5050-5059. [PMID: 30723158 DOI: 10.1074/jbc.ra119.007459] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Revised: 02/01/2019] [Indexed: 11/06/2022] Open
Abstract
TraB is an FtsK-like DNA translocase responsible for conjugative plasmid transfer in mycelial Streptomyces Unlike other conjugative systems, which depend on a type IV secretion system, Streptomyces requires only TraB protein to transfer the plasmid as dsDNA. The γ-domain of this protein specifically binds to repeated 8-bp motifs on the plasmid sequence, following a mechanism that is reminiscent of the FtsK/SpoIIIE chromosome segregation system. In this work, we purified and characterized the enzymatic activity of TraB, revealing that it is a DNA-dependent ATPase that is highly stimulated by dsDNA substrates. Interestingly, we found that unlike the SpoIIIE protein, the γ-domain of TraB does not confer sequence-specific ATPase stimulation. We also found that TraB binds G-quadruplex DNA structures with higher affinity than TraB-recognition sequences (TRSs). An EM-based structural analysis revealed that TraB tends to assemble as large complexes comprising four TraB hexamers, which might be a prerequisite for DNA translocation across cell membranes. In summary, our findings shed light on the molecular mechanism used by the DNA-translocating motor TraB, which may be shared by other membrane-associated machineries involved in DNA binding and translocation.
Collapse
Affiliation(s)
- Eric Amado
- From the Departamento de Biología Molecular and Instituto de Biomedicina y Biotecnología de Cantabria (IBBTEC), Universidad de Cantabria-Consejo Superior de Investigaciones Científicas, Santander, Spain and
| | - Günther Muth
- Interfakultaeres Institut für Mikrobiologie und Infektionsmedizin Tuebingen IMIT, Mikrobiologie/Biotechnologie, Eberhard Karls Universitaet Tuebingen, 72074 Tuebingen, Germany
| | - Ignacio Arechaga
- From the Departamento de Biología Molecular and Instituto de Biomedicina y Biotecnología de Cantabria (IBBTEC), Universidad de Cantabria-Consejo Superior de Investigaciones Científicas, Santander, Spain and
| | - Elena Cabezón
- From the Departamento de Biología Molecular and Instituto de Biomedicina y Biotecnología de Cantabria (IBBTEC), Universidad de Cantabria-Consejo Superior de Investigaciones Científicas, Santander, Spain and
| |
Collapse
|
37
|
Steele TWE, Samsó M. The FKBP12 subunit modifies the long-range allosterism of the ryanodine receptor. J Struct Biol 2019; 205:180-188. [PMID: 30641143 DOI: 10.1016/j.jsb.2018.12.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2018] [Revised: 12/20/2018] [Accepted: 12/21/2018] [Indexed: 02/06/2023]
Abstract
Ryanodine receptors (RyRs) are large conductance intracellular channels controlling intracellular calcium homeostasis in myocytes, neurons, and other cell types. Loss of RyR's constitutive cytoplasmic partner FKBP results in channel sensitization, dominant subconductance states, and increased cytoplasmic Ca2+. FKBP12 binds to RyR1's cytoplasmic assembly 130 Å away from the ion gate at four equivalent sites in the RyR1 tetramer. To understand how FKBP12 binding alters RyR1's channel properties, we studied the 3D structure of RyR1 alone in the closed conformation in the context of the open and closed conformations of FKBP12-bound RyR1. We analyzed the metrics of conformational changes of existing structures, the structure of the ion gate, and carried out multivariate statistical analysis of thousands of individual cryoEM RyR1 particles. We find that under closed state conditions, in the presence of FKBP12, the cytoplasmic domain of RyR1 adopts an upward conformation, whereas absence of FKBP12 results in a relaxed conformation, while the ion gate remains closed. The relaxed conformation is intermediate between the RyR1-FKBP12 complex closed (upward) and open (downward) conformations. The closed-relaxed conformation of RyR1 appears to be consistent with a lower energy barrier separating the closed and open states of RyR1-FKBP12, and suggests that FKBP12 plays an important role by restricting conformations within RyR1's conformational landscape.
Collapse
Affiliation(s)
- Tyler W E Steele
- Department of Physiology and Biophysics, Virginia Commonwealth University, Richmond, VA 23298, United States
| | - Montserrat Samsó
- Department of Physiology and Biophysics, Virginia Commonwealth University, Richmond, VA 23298, United States.
| |
Collapse
|
38
|
Anand DV, Meng Z, Xia K. A complex multiscale virtual particle model based elastic network model (CMVP-ENM) for the normal mode analysis of biomolecular complexes. Phys Chem Chem Phys 2019; 21:4359-4366. [DOI: 10.1039/c8cp07442a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The CMVP-ENM for virus normal mode analysis. With a special ratio parameter, CMVP-ENM can characterize the multi-material properties of biomolecular complexes and systematically enhance or suppress the modes for different components.
Collapse
Affiliation(s)
- D. Vijay Anand
- Division of Mathematical Sciences
- School of Physical and Mathematical Sciences
- Nanyang Technological University
- Singapore
| | - Zhenyu Meng
- School of Biological Sciences
- Nanyang Technological University
- Singapore
| | - Kelin Xia
- Division of Mathematical Sciences
- School of Physical and Mathematical Sciences
- Nanyang Technological University
- Singapore
- School of Biological Sciences
| |
Collapse
|
39
|
Ramesh R, Lim XX, Raghuvamsi PV, Wu C, Wong SM, Anand GS. Uncovering metastability and disassembly hotspots in whole viral particles. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2018; 143:5-12. [PMID: 30553754 DOI: 10.1016/j.pbiomolbio.2018.12.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Revised: 10/26/2018] [Accepted: 12/12/2018] [Indexed: 01/18/2023]
Abstract
Viruses are metastable macromolecular assemblies that toggle between multiple conformational states through molecular rearrangements that are critical for mediating viral host entry. Viruses respond to different host specific environmental cues to form disassembly intermediates for the eventual release of genomic material required for replication. Although static snapshots of these intermediates have been captured through structural techniques such as X-ray crystallography and cryo-EM, the mechanistic details of these conformational rearrangements underpinning viral metastability have been poorly understood. Amide hydrogen deuterium exchange mass spectrometry (HDXMS) is a powerful tool that measures hydrogen bonding propensities to probe changes in the dynamics of different macromolecular interactions. Chaotropic agents such as urea can be used to disrupt hydrogen bonds between different subunits, thereby ranking regions of the virus that are critical in maintaining viral stability. By controlled urea denaturation with HDXMS, we have identified specific loci in a Turnip Crinkle Virus (TCV) model showing increased deuterium exchange with even minimally disruptive concentrations of urea. These loci represent dynamic disassembly hotspots. These hotspots are predominantly present at the quaternary contacts at the 3-fold and 5-fold axes. This approach can be applied to detect vulnerabilities in virus icosahedral structures to uncover the molecular mechanism of viral disassembly.
Collapse
Affiliation(s)
- Ranita Ramesh
- Department of Biological Sciences, National University of Singapore, 117543, Singapore
| | - Xin Xiang Lim
- Department of Biological Sciences, National University of Singapore, 117543, Singapore
| | | | - Chao Wu
- Department of Biological Sciences, National University of Singapore, 117543, Singapore
| | - Sek Man Wong
- Department of Biological Sciences, National University of Singapore, 117543, Singapore; Temasek Life Sciences Laboratory, Singapore, 117604, Singapore
| | | |
Collapse
|
40
|
Che C, Lin R, Zeng X, Elmaaroufi K, Galeotti J, Xu M. Improved deep learning-based macromolecules structure classification from electron cryo-tomograms. MACHINE VISION AND APPLICATIONS 2018; 29:1227-1236. [PMID: 31511756 PMCID: PMC6738941 DOI: 10.1007/s00138-018-0949-4] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Revised: 01/16/2018] [Accepted: 05/18/2018] [Indexed: 05/30/2023]
Abstract
Cellular processes are governed by macromolecular complexes inside the cell. Study of the native structures of macromolecular complexes has been extremely difficult due to lack of data. With recent breakthroughs in Cellular Electron Cryo-Tomography (CECT) 3D imaging technology, it is now possible for researchers to gain accesses to fully study and understand the macro-molecular structures single cells. However, systematic recovery of macromolecular structures from CECT is very difficult due to high degree of structural complexity and practical imaging limitations. Specifically, we proposed a deep learning-based image classification approach for large-scale systematic macromolecular structure separation from CECT data. However, our previous work was only a very initial step toward exploration of the full potential of deep learning-based macromolecule separation. In this paper, we focus on improving classification performance by proposing three newly designed individual CNN models: an extended version of (Deep Small Receptive Field) DSRF3D, donated as DSRF3D-v2, a 3D residual block-based neural network, named as RB3D, and a convolutional 3D (C3D)-based model, CB3D. We compare them with our previously developed model (DSRF3D) on 12 datasets with different SNRs and tilt angle ranges. The experiments show that our new models achieved significantly higher classification accuracies. The accuracies are not only higher than 0.9 on normal datasets, but also demonstrate potentials to operate on datasets with high levels of noises and missing wedge effects presented.
Collapse
Affiliation(s)
- Chengqian Che
- The Robotics Institute, Carnegie Mellon University,Pittsburgh, USA
| | - Ruogu Lin
- Department of Automation, Tsinghua University, Beijing, China
| | - Xiangrui Zeng
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, USA
| | - Karim Elmaaroufi
- Department of Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, USA
| | - John Galeotti
- The Robotics Institute, Carnegie Mellon University,Pittsburgh, USA
| | - Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, USA
| |
Collapse
|
41
|
Kappel K, Liu S, Larsen KP, Skiniotis G, Puglisi EV, Puglisi JD, Zhou ZH, Zhao R, Das R. De novo computational RNA modeling into cryo-EM maps of large ribonucleoprotein complexes. Nat Methods 2018; 15:947-954. [PMID: 30377372 PMCID: PMC6636682 DOI: 10.1038/s41592-018-0172-2] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2018] [Accepted: 07/31/2018] [Indexed: 12/19/2022]
Abstract
Increasingly, cryo-electron microscopy (cryo-EM) is used to determine the structures of RNA-protein assemblies, but nearly all maps determined with this method have biologically important regions where the local resolution does not permit RNA coordinate tracing. To address these omissions, we present de novo ribonucleoprotein modeling in real space through assembly of fragments together with experimental density in Rosetta (DRRAFTER). We show that DRRAFTER recovers near-native models for a diverse benchmark set of RNA-protein complexes including the spliceosome, mitochondrial ribosome, and CRISPR-Cas9-sgRNA complexes; rigorous blind tests include yeast U1 snRNP and spliceosomal P complex maps. Additionally, to aid in model interpretation, we present a method for reliable in situ estimation of DRRAFTER model accuracy. Finally, we apply DRRAFTER to recently determined maps of telomerase, the HIV-1 reverse transcriptase initiation complex, and the packaged MS2 genome, demonstrating the acceleration of accurate model building in challenging cases.
Collapse
Affiliation(s)
- Kalli Kappel
- Biophysics Program, Stanford University, Stanford, CA, USA
| | - Shiheng Liu
- Electron Imaging Center for Nanomachines, California NanoSystems Institute, University of California, Los Angeles (UCLA), Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, UCLA, Los Angeles, CA, USA
| | - Kevin P Larsen
- Biophysics Program, Stanford University, Stanford, CA, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Georgios Skiniotis
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
- Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
| | | | - Joseph D Puglisi
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Z Hong Zhou
- Electron Imaging Center for Nanomachines, California NanoSystems Institute, University of California, Los Angeles (UCLA), Los Angeles, CA, USA
- Department of Microbiology, Immunology, and Molecular Genetics, UCLA, Los Angeles, CA, USA
| | - Rui Zhao
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver Anschutz Medical Campus, Aurora, CO, USA
| | - Rhiju Das
- Biophysics Program, Stanford University, Stanford, CA, USA.
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA, USA.
- Department of Physics, Stanford University, Stanford, CA, USA.
| |
Collapse
|
42
|
Mori T, Kulik M, Miyashita O, Jung J, Tama F, Sugita Y. Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning. Structure 2018; 27:161-174.e3. [PMID: 30344106 DOI: 10.1016/j.str.2018.09.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2018] [Revised: 06/22/2018] [Accepted: 09/18/2018] [Indexed: 01/21/2023]
Abstract
Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca2+-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.
Collapse
Affiliation(s)
- Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Marta Kulik
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Osamu Miyashita
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Jaewoon Jung
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan; RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Florence Tama
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Department of Physics, Graduate School of Science, and Institute of Transformative Bio-Molecules, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan; RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; RIKEN Center for Biosystems Dynamics Research, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| |
Collapse
|
43
|
Gupta M, Sharma R, Kumar A. Docking techniques in pharmacology: How much promising? Comput Biol Chem 2018; 76:210-217. [PMID: 30067954 DOI: 10.1016/j.compbiolchem.2018.06.005] [Citation(s) in RCA: 86] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2017] [Revised: 02/21/2018] [Accepted: 06/30/2018] [Indexed: 01/01/2023]
|
44
|
|
45
|
Xu M, Chai X, Muthakana H, Liang X, Yang G, Zeev-Ben-Mordehai T, Xing EP. Deep learning-based subdivision approach for large scale macromolecules structure recovery from electron cryo tomograms. Bioinformatics 2018; 33:i13-i22. [PMID: 28881965 PMCID: PMC5946875 DOI: 10.1093/bioinformatics/btx230] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Motivation Cellular Electron CryoTomography (CECT) enables 3D visualization of cellular organization at near-native state and in sub-molecular resolution, making it a powerful tool for analyzing structures of macromolecular complexes and their spatial organizations inside single cells. However, high degree of structural complexity together with practical imaging limitations makes the systematic de novo discovery of structures within cells challenging. It would likely require averaging and classifying millions of subtomograms potentially containing hundreds of highly heterogeneous structural classes. Although it is no longer difficult to acquire CECT data containing such amount of subtomograms due to advances in data acquisition automation, existing computational approaches have very limited scalability or discrimination ability, making them incapable of processing such amount of data. Results To complement existing approaches, in this article we propose a new approach for subdividing subtomograms into smaller but relatively homogeneous subsets. The structures in these subsets can then be separately recovered using existing computation intensive methods. Our approach is based on supervised structural feature extraction using deep learning, in combination with unsupervised clustering and reference-free classification. Our experiments show that, compared with existing unsupervised rotation invariant feature and pose-normalization based approaches, our new approach achieves significant improvements in both discrimination ability and scalability. More importantly, our new approach is able to discover new structural classes and recover structures that do not exist in training data. Availability and Implementation Source code freely available at http://www.cs.cmu.edu/∼mxu1/software. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Min Xu
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xiaoqi Chai
- Biomedical Engineering Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Hariank Muthakana
- Computer Science Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xiaodan Liang
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Ge Yang
- Biomedical Engineering Department, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Tzviya Zeev-Ben-Mordehai
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK
| | - Eric P Xing
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, PA, USA
| |
Collapse
|
46
|
Xia K. Multiscale virtual particle based elastic network model (MVP-ENM) for normal mode analysis of large-sized biomolecules. Phys Chem Chem Phys 2018; 20:658-669. [PMID: 29227479 DOI: 10.1039/c7cp07177a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
In this paper, a multiscale virtual particle based elastic network model (MVP-ENM) is proposed for the normal mode analysis of large-sized biomolecules. The multiscale virtual particle (MVP) model is proposed for the discretization of biomolecular density data. With this model, large-sized biomolecular structures can be coarse-grained into virtual particles such that a balance between model accuracy and computational cost can be achieved. An elastic network is constructed by assuming "connections" between virtual particles. The connection is described by a special harmonic potential function, which considers the influence from both the mass distributions and distance relations of the virtual particles. Two independent models, i.e., the multiscale virtual particle based Gaussian network model (MVP-GNM) and the multiscale virtual particle based anisotropic network model (MVP-ANM), are proposed. It has been found that in the Debye-Waller factor (B-factor) prediction, the results from our MVP-GNM with a high resolution are as good as the ones from GNM. Even with low resolutions, our MVP-GNM can still capture the global behavior of the B-factor very well with mismatches predominantly from the regions with large B-factor values. Further, it has been demonstrated that the low-frequency eigenmodes from our MVP-ANM are highly consistent with the ones from ANM even with very low resolutions and a coarse grid. Finally, the great advantage of MVP-ANM model for large-sized biomolecules has been demonstrated by using two poliovirus virus structures. The paper ends with a conclusion.
Collapse
Affiliation(s)
- Kelin Xia
- Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371.
| |
Collapse
|
47
|
Harrer N, Schindler CEM, Bruetzel LK, Forné I, Ludwigsen J, Imhof A, Zacharias M, Lipfert J, Mueller-Planitz F. Structural Architecture of the Nucleosome Remodeler ISWI Determined from Cross-Linking, Mass Spectrometry, SAXS, and Modeling. Structure 2018; 26:282-294.e6. [PMID: 29395785 DOI: 10.1016/j.str.2017.12.015] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2017] [Revised: 10/25/2017] [Accepted: 12/27/2017] [Indexed: 11/17/2022]
Abstract
Chromatin remodeling factors assume critical roles by regulating access to nucleosomal DNA. To determine the architecture of the Drosophila ISWI remodeling enzyme, we developed an integrative structural approach that combines protein cross-linking, mass spectrometry, small-angle X-ray scattering, and computational modeling. The resulting structural model shows the ATPase module in a resting state with both ATPase lobes twisted against each other, providing support for a conformation that was recently trapped by crystallography. The autoinhibiting NegC region does not protrude from the ATPase module as suggested previously. The regulatory NTR domain is located near both ATPase lobes. The full-length enzyme is flexible and can adopt a compact structure in solution with the C-terminal HSS domain packing against the ATPase module. Our data imply a series of conformational changes upon activation of the enzyme and illustrate how the NTR, NegC, and HSS domains contribute to regulation of the ATPase module.
Collapse
Affiliation(s)
- Nadine Harrer
- Molecular Biology, Biomedical Center, Faculty of Medicine, LMU Munich, 82152 Martinsried, Germany
| | - Christina E M Schindler
- Physics Department (T38), Technical University of Munich, 85748 Garching, Germany; Center for Integrated Protein Science Munich, 81377 Munich, Germany
| | - Linda K Bruetzel
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, 80799 Munich, Germany
| | - Ignasi Forné
- Molecular Biology, Biomedical Center, Faculty of Medicine, LMU Munich, 82152 Martinsried, Germany
| | - Johanna Ludwigsen
- Molecular Biology, Biomedical Center, Faculty of Medicine, LMU Munich, 82152 Martinsried, Germany
| | - Axel Imhof
- Molecular Biology, Biomedical Center, Faculty of Medicine, LMU Munich, 82152 Martinsried, Germany
| | - Martin Zacharias
- Physics Department (T38), Technical University of Munich, 85748 Garching, Germany; Center for Integrated Protein Science Munich, 81377 Munich, Germany
| | - Jan Lipfert
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, LMU Munich, 80799 Munich, Germany.
| | - Felix Mueller-Planitz
- Molecular Biology, Biomedical Center, Faculty of Medicine, LMU Munich, 82152 Martinsried, Germany.
| |
Collapse
|
48
|
Al Nasr K, Yousef F, Jebril R, Jones C. Analytical Approaches to Improve Accuracy in Solving the Protein Topology Problem. Molecules 2018; 23:E28. [PMID: 29360779 PMCID: PMC6017786 DOI: 10.3390/molecules23020028] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Revised: 01/19/2018] [Accepted: 01/19/2018] [Indexed: 11/17/2022] Open
Abstract
To take advantage of recent advances in genomics and proteomics it is critical that the three-dimensional physical structure of biological macromolecules be determined. Cryo-Electron Microscopy (cryo-EM) is a promising and improving method for obtaining this data, however resolution is often not sufficient to directly determine the atomic scale structure. Despite this, information for secondary structure locations is detectable. De novo modeling is a computational approach to modeling these macromolecular structures based on cryo-EM derived data. During de novo modeling a mapping between detected secondary structures and the underlying amino acid sequence must be identified. DP-TOSS (Dynamic Programming for determining the Topology Of Secondary Structures) is one tool that attempts to automate the creation of this mapping. By treating the correspondence between the detected structures and the structures predicted from sequence data as a constraint graph problem DP-TOSS achieved good accuracy in its original iteration. In this paper, we propose modifications to the scoring methodology of DP-TOSS to improve its accuracy. Three scoring schemes were applied to DP-TOSS and tested: (i) a skeleton-based scoring function; (ii) a geometry-based analytical function; and (iii) a multi-well potential energy-based function. A test of 25 proteins shows that a combination of these schemes can improve the performance of DP-TOSS to solve the topology determination problem for macromolecule proteins.
Collapse
Affiliation(s)
- Kamal Al Nasr
- Department of Computer Science, Tennessee State University, Nashville, TN 37209, USA.
| | - Feras Yousef
- Department of Mathematics, The University of Jordan, Amman 11942, Jordan.
| | - Ruba Jebril
- Department of Computer Science, Tennessee State University, Nashville, TN 37209, USA.
| | - Christopher Jones
- Department of Computer Science, Tennessee State University, Nashville, TN 37209, USA.
| |
Collapse
|
49
|
Mei K, Li Y, Wang S, Shao G, Wang J, Ding Y, Luo G, Yue P, Liu JJ, Wang X, Dong MQ, Wang HW, Guo W. Cryo-EM structure of the exocyst complex. Nat Struct Mol Biol 2018; 25:139-146. [PMID: 29335562 PMCID: PMC5971111 DOI: 10.1038/s41594-017-0016-2] [Citation(s) in RCA: 93] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2017] [Accepted: 12/07/2017] [Indexed: 12/22/2022]
Abstract
The exocyst is an evolutionarily conserved octameric protein complex that mediates the tethering of post-Golgi secretory vesicles to the plasma membrane during exocytosis and is implicated in many cellular processes such as cell polarization, cytokinesis, ciliogenesis and tumor invasion. Using cryo-EM and chemical cross-linking MS (CXMS), we solved the structure of the Saccharomyces cerevisiae exocyst complex at an average resolution of 4.4 Å. Our model revealed the architecture of the exocyst and led to the identification of the helical bundles that mediate the assembly of the complex at its core. Sequence analysis suggests that these regions are evolutionarily conserved across eukaryotic systems. Additional cell biological data suggest a mechanism for exocyst assembly that leads to vesicle tethering at the plasma membrane.
Collapse
Affiliation(s)
- Kunrong Mei
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | - Yan Li
- Ministry of Education Key Laboratory of Protein Sciences, Tsinghua University, Beijing, China.,Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing, China.,School of Life Sciences, Tsinghua University, Beijing, China.,Tsinghua-Peking Joint Center for Life Sciences, Beijing, China
| | - Shaoxiao Wang
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | - Guangcan Shao
- National Institute of Biological Sciences, Beijing, China
| | - Jia Wang
- Ministry of Education Key Laboratory of Protein Sciences, Tsinghua University, Beijing, China.,Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing, China.,School of Life Sciences, Tsinghua University, Beijing, China.,Tsinghua-Peking Joint Center for Life Sciences, Beijing, China
| | - Yuehe Ding
- National Institute of Biological Sciences, Beijing, China
| | - Guangzuo Luo
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | - Peng Yue
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | - Jun-Jie Liu
- Ministry of Education Key Laboratory of Protein Sciences, Tsinghua University, Beijing, China.,Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing, China.,School of Life Sciences, Tsinghua University, Beijing, China.,Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Xinquan Wang
- Ministry of Education Key Laboratory of Protein Sciences, Tsinghua University, Beijing, China.,Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing, China.,School of Life Sciences, Tsinghua University, Beijing, China
| | - Meng-Qiu Dong
- National Institute of Biological Sciences, Beijing, China
| | - Hong-Wei Wang
- Ministry of Education Key Laboratory of Protein Sciences, Tsinghua University, Beijing, China. .,Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing, China. .,School of Life Sciences, Tsinghua University, Beijing, China. .,Tsinghua-Peking Joint Center for Life Sciences, Beijing, China.
| | - Wei Guo
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
50
|
Assessment of data-assisted prediction by inclusion of crosslinking/mass-spectrometry and small angle X-ray scattering data in the 12thCritical Assessment of protein Structure Prediction experiment. Proteins 2017; 86 Suppl 1:215-227. [DOI: 10.1002/prot.25442] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2017] [Revised: 11/16/2017] [Accepted: 12/10/2017] [Indexed: 12/26/2022]
|