1
|
Gien H, Morse M, McCauley MJ, Rouzina I, Gorelick RJ, Williams MC. Cationic Residues of the HIV-1 Nucleocapsid Protein Enable DNA Condensation to Maintain Viral Core Particle Stability during Reverse Transcription. Viruses 2024; 16:872. [PMID: 38932164 PMCID: PMC11209390 DOI: 10.3390/v16060872] [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: 03/29/2024] [Revised: 04/30/2024] [Accepted: 05/23/2024] [Indexed: 06/28/2024] Open
Abstract
The HIV-1 nucleocapsid protein (NC) is a multifunctional viral protein necessary for HIV-1 replication. Recent studies have demonstrated that reverse transcription (RT) completes in the intact viral capsid, and the timing of RT and uncoating are correlated. How the small viral core stably contains the ~10 kbp double stranded (ds) DNA product of RT, and the role of NC in this process, are not well understood. We showed previously that NC binds and saturates dsDNA in a non-specific electrostatic binding mode that triggers uniform DNA self-attraction, condensing dsDNA into a tight globule against extending forces up to 10 pN. In this study, we use optical tweezers and atomic force microscopy to characterize the role of NC's basic residues in dsDNA condensation. Basic residue mutations of NC lead to defective interaction with the dsDNA substrate, with the constant force plateau condensation observed with wild-type (WT) NC missing or diminished. These results suggest that NC's high positive charge is essential to its dsDNA condensing activity, and electrostatic interactions involving NC's basic residues are responsible in large part for the conformation, size, and stability of the dsDNA-protein complex inside the viral core. We observe DNA re-solubilization and charge reversal in the presence of excess NC, consistent with the electrostatic nature of NC-induced DNA condensation. Previous studies of HIV-1 replication in the presence of the same cationic residue mutations in NC showed significant defects in both single- and multiple-round viral infectivity. Although NC participates in many stages of viral replication, our results are consistent with the hypothesis that cationic residue mutations inhibit genomic DNA condensation, resulting in increased premature capsid uncoating and contributing to viral replication defects.
Collapse
Affiliation(s)
- Helena Gien
- Department of Physics, Northeastern University, Boston, MA 02115, USA; (H.G.); (M.M.); (M.J.M.)
| | - Michael Morse
- Department of Physics, Northeastern University, Boston, MA 02115, USA; (H.G.); (M.M.); (M.J.M.)
| | - Micah J. McCauley
- Department of Physics, Northeastern University, Boston, MA 02115, USA; (H.G.); (M.M.); (M.J.M.)
| | - Ioulia Rouzina
- Department of Chemistry and Biochemistry, Center for Retroviral Research and Center for RNA Biology, The Ohio State University, Columbus, OH 43210, USA;
| | - Robert J. Gorelick
- AIDS and Cancer Virus Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA;
| | - Mark C. Williams
- Department of Physics, Northeastern University, Boston, MA 02115, USA; (H.G.); (M.M.); (M.J.M.)
| |
Collapse
|
2
|
Chen H, Yan G, Wen MH, Brooks KN, Zhang Y, Huang PS, Chen TY. Advancements and Practical Considerations for Biophysical Research: Navigating the Challenges and Future of Super-resolution Microscopy. CHEMICAL & BIOMEDICAL IMAGING 2024; 2:331-344. [PMID: 38817319 PMCID: PMC11134610 DOI: 10.1021/cbmi.4c00019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/10/2024] [Revised: 04/06/2024] [Accepted: 04/10/2024] [Indexed: 06/01/2024]
Abstract
The introduction of super-resolution microscopy (SRM) has significantly advanced our understanding of cellular and molecular dynamics, offering a detailed view previously beyond our reach. Implementing SRM in biophysical research, however, presents numerous challenges. This review addresses the crucial aspects of utilizing SRM effectively, from selecting appropriate fluorophores and preparing samples to analyzing complex data sets. We explore recent technological advancements and methodological improvements that enhance the capabilities of SRM. Emphasizing the integration of SRM with other analytical methods, we aim to overcome inherent limitations and expand the scope of biological insights achievable. By providing a comprehensive guide for choosing the most suitable SRM methods based on specific research objectives, we aim to empower researchers to explore complex biological processes with enhanced precision and clarity, thereby advancing the frontiers of biophysical research.
Collapse
Affiliation(s)
- Huanhuan Chen
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Guangjie Yan
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Meng-Hsuan Wen
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Kameron N. Brooks
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Yuteng Zhang
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Pei-San Huang
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| | - Tai-Yen Chen
- Department of Chemistry, University of Houston, Houston, Texas 77204, United States
| |
Collapse
|
3
|
Taylor IA, Fassati A. The capsid revolution. J Mol Cell Biol 2024; 15:mjad076. [PMID: 38037430 PMCID: PMC11193064 DOI: 10.1093/jmcb/mjad076] [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: 07/25/2023] [Revised: 11/23/2023] [Accepted: 11/29/2023] [Indexed: 12/02/2023] Open
Abstract
Lenacapavir, targeting the human immunodeficiency virus type-1 (HIV-1) capsid, is the first-in-class antiretroviral drug recently approved for clinical use. The development of Lenacapavir is attributed to the remarkable progress in our understanding of the capsid protein made during the last few years. Considered little more than a component of the virus shell to be shed early during infection, the capsid has been found to be a key player in the HIV-1 life cycle by interacting with multiple host factors, entering the nucleus, and directing integration. Here, we describe the key advances that led to this 'capsid revolution'.
Collapse
Affiliation(s)
- Ian A Taylor
- Macromolecular Structure Laboratory, The Francis Crick Institute, London NW1 1AT, UK
| | - Ariberto Fassati
- Division of Infection and Immunity, University College London, London WC1E 6JF, UK
- Institute of Immunity and Transplantation, University College London, London NW3 2PP, UK
| |
Collapse
|
4
|
Coshic K, Maffeo C, Winogradoff D, Aksimentiev A. The structure and physical properties of a packaged bacteriophage particle. Nature 2024; 627:905-914. [PMID: 38448589 PMCID: PMC11196859 DOI: 10.1038/s41586-024-07150-4] [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: 09/01/2023] [Accepted: 02/01/2024] [Indexed: 03/08/2024]
Abstract
A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.
Collapse
Affiliation(s)
- Kush Coshic
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Christopher Maffeo
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - David Winogradoff
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Aleksei Aksimentiev
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
5
|
Shrivastav G, Borkotoky S, Dey D, Singh B, Malhotra N, Azad K, Jayaram B, Agarwal M, Banerjee M. Structure and energetics guide dynamic behaviour in a T = 3 icosahedral virus capsid. Biophys Chem 2024; 305:107152. [PMID: 38113782 DOI: 10.1016/j.bpc.2023.107152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 12/03/2023] [Accepted: 12/04/2023] [Indexed: 12/21/2023]
Abstract
Although virus capsids appear as rigid, symmetric particles in experimentally determined structures; biochemical studies suggest a significant degree of structural flexibility in the particles. We carried out all-atom simulations on the icosahedral capsid of an insect virus, Flock House Virus, which show intriguing differences in the degree of flexibility of quasi-equivalent capsid subunits consistent with previously described biological behaviour. The flexibility of all the β and γ subunits of the protein and RNA fragments is analysed and compared. Both γA subunit and RNA fragment exhibit higher flexibility than the γB and γC subunits. The capsid shell is permeable to the bidirectional movement of water molecules, and the movement is heavily influenced by the geometry of the capsid shell along specific symmetry axes. In comparison to the symmetry axes along I5 and I3, the I2 axis exhibits a slightly higher water content. This enriched water environment along I2 could play a pivotal role in facilitating the structural transitions necessary for RNA release, shedding some light on the intricate and dynamic processes underlying the viral life cycle. Our study suggests that the physical characterization of whole virus capsids is the key to identifying biologically relevant transition states in the virus life cycle and understanding the basis of virus infectivity.
Collapse
Affiliation(s)
- Gourav Shrivastav
- Department of Chemical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.
| | - Subhomoi Borkotoky
- Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Debajit Dey
- Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Bhumika Singh
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Nidhi Malhotra
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Kimi Azad
- Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - B Jayaram
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Manish Agarwal
- Computer Services Centre, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.
| | - Manidipa Banerjee
- Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.
| |
Collapse
|
6
|
Sha H, Zhu F. Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model. J Phys Chem B 2024; 128:960-972. [PMID: 38251836 DOI: 10.1021/acs.jpcb.3c06881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2024]
Abstract
HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.
Collapse
Affiliation(s)
- Hao Sha
- Department of Physics, Indiana University─Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Fangqiang Zhu
- Department of Physics, Indiana University─Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
7
|
Antolínez S, Jones PE, Phillips JC, Hadden-Perilla JA. AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD. J Chem Inf Model 2024; 64:543-554. [PMID: 38176097 PMCID: PMC10806814 DOI: 10.1021/acs.jcim.3c01648] [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: 10/10/2023] [Revised: 12/18/2023] [Accepted: 12/18/2023] [Indexed: 01/06/2024]
Abstract
All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.
Collapse
Affiliation(s)
- Santiago Antolínez
- Department
of Chemistry and Biochemistry, University
of Delaware, Newark, Delaware 19716, United States
| | - Peter Eugene Jones
- Department
of Chemistry and Biochemistry, University
of Delaware, Newark, Delaware 19716, United States
| | - James C. Phillips
- National
Center for Supercomputing Applications, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Jodi A. Hadden-Perilla
- Department
of Chemistry and Biochemistry, University
of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
8
|
Brown C, Agarwal A, Luque A. pyCapsid: identifying dominant dynamics and quasi-rigid mechanical units in protein shells. Bioinformatics 2024; 40:btad761. [PMID: 38113434 PMCID: PMC10786678 DOI: 10.1093/bioinformatics/btad761] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 11/01/2023] [Accepted: 12/15/2023] [Indexed: 12/21/2023] Open
Abstract
SUMMARY pyCapsid is a Python package developed to facilitate the characterization of the dynamics and quasi-rigid mechanical units of protein shells and other protein complexes. The package was developed in response to the rapid increase of high-resolution structures, particularly capsids of viruses, requiring multiscale biophysical analyses. Given a protein shell, pyCapsid generates the collective vibrations of its amino-acid residues, identifies quasi-rigid mechanical regions associated with the disassembly of the structure, and maps the results back to the input proteins for interpretation. pyCapsid summarizes the main results in a report that includes publication-quality figures. AVAILABILITY AND IMPLEMENTATION pyCapsid's source code is available under MIT License on GitHub. It is compatible with Python 3.8-3.10 and has been deployed in two leading Python package-management systems, PIP and Conda. Installation instructions and tutorials are available in the online documentation and in the pyCapsid's YouTube playlist. In addition, a cloud-based implementation of pyCapsid is available as a Google Colab notebook. pyCapsid Colab does not require installation and generates the same report and outputs as the installable version. Users can post issues regarding pyCapsid in the repository's issues section.
Collapse
Affiliation(s)
- Colin Brown
- Viral Information Institute, San Diego State University, San Diego, CA 92116, United States
- Department of Physics, San Diego State University, San Diego, CA 92116, United States
| | - Anuradha Agarwal
- Viral Information Institute, San Diego State University, San Diego, CA 92116, United States
- Computational Science Research Center, San Diego State University, San Diego, CA 92116, United States
| | - Antoni Luque
- Viral Information Institute, San Diego State University, San Diego, CA 92116, United States
- Computational Science Research Center, San Diego State University, San Diego, CA 92116, United States
- Department of Mathematics and Statistics, San Diego State University, San Diego, CA 92116, United States
- Department of Biology, University of Miami, Coral Gables, FL 33146, United States
| |
Collapse
|
9
|
Bogetti A, Leung JMG, Chong LT. LPATH: A Semiautomated Python Tool for Clustering Molecular Pathways. J Chem Inf Model 2023; 63:7610-7616. [PMID: 38048485 PMCID: PMC10751797 DOI: 10.1021/acs.jcim.3c01318] [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: 08/18/2023] [Revised: 10/14/2023] [Accepted: 11/09/2023] [Indexed: 12/06/2023]
Abstract
The pathways by which a molecular process transitions to a target state are highly sought-after as direct views of a transition mechanism. While great strides have been made in the physics-based simulation of such pathways, the analysis of these pathways can be a major challenge due to their diversity and variable lengths. Here, we present the LPATH Python tool, which implements a semiautomated method for linguistics-assisted clustering of pathways into distinct classes (or routes). This method involves three steps: 1) discretizing the configurational space into key states, 2) extracting a text-string sequence of key visited states for each pathway, and 3) pairwise matching of pathways based on a text-string similarity score. To circumvent the prohibitive memory requirements of the first step, we have implemented a general two-stage method for clustering conformational states that exploits machine learning. LPATH is primarily designed for use with the WESTPA software for weighted ensemble simulations; however, the tool can also be applied to conventional simulations. As demonstrated for the C7eq to C7ax conformational transition of the alanine dipeptide, LPATH provides physically reasonable classes of pathways and corresponding probabilities.
Collapse
Affiliation(s)
- Anthony
T. Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Jeremy M. G. Leung
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Lillian T. Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
10
|
Ahmed M, Maldonado AM, Durrant JD. From Byte to Bench to Bedside: Molecular Dynamics Simulations and Drug Discovery. ARXIV 2023:arXiv:2311.16946v1. [PMID: 38076508 PMCID: PMC10705576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Molecular dynamics (MD) simulations and computer-aided drug design (CADD) have advanced substantially over the past two decades, thanks to continuous computer hardware and software improvements. Given these advancements, MD simulations are poised to become even more powerful tools for investigating the dynamic interactions between potential small-molecule drugs and their target proteins, with significant implications for pharmacological research.
Collapse
Affiliation(s)
- Mayar Ahmed
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Alex M. Maldonado
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Jacob D. Durrant
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| |
Collapse
|
11
|
Bogetti AT, Leung JMG, Chong LT. LPATH: A semi-automated Python tool for clustering molecular pathways. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.17.553774. [PMID: 37645995 PMCID: PMC10462149 DOI: 10.1101/2023.08.17.553774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
The pathways by which a molecular process transitions to a target state are highly sought-after as direct views of a transition mechanism. While great strides have been made in the physics-based simulation of such pathways, the analysis of these pathways can be a major challenge due to their diversity and variable lengths. Here we present the LPATH Python tool, which implements a semi-automated method for linguistics-assisted clustering of pathways into distinct classes (or routes). This method involves three steps: 1) discretizing the configurational space into key states, 2) extracting a text-string sequence of key visited states for each pathway, and 3) pairwise matching of pathways based on a text-string similarity score. To circumvent the prohibitive memory requirements of the first step, we have implemented a general two-stage method for clustering conformational states that exploits machine learning. LPATH is primarily designed for use with the WESTPA software for weighted ensemble simulations; however, the tool can also be applied to conventional simulations. As demonstrated for the C7eq to C7ax conformational transition of alanine dipeptide, LPATH provides physically reasonable classes of pathways and corresponding probabilities.
Collapse
Affiliation(s)
- Anthony T. Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - Jeremy M. G. Leung
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| | - Lillian T. Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260
| |
Collapse
|
12
|
Callil-Soares PH, Biasi LCK, Pessoa Filho PDA. Effect of preprocessing and simulation parameters on the performance of molecular docking studies. J Mol Model 2023; 29:251. [PMID: 37452150 DOI: 10.1007/s00894-023-05637-x] [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: 02/13/2023] [Accepted: 06/26/2023] [Indexed: 07/18/2023]
Abstract
CONTEXT Molecular docking is an important and rapid tool that provides a comprehensive view of different molecular mechanisms. It is often used to verify the binding interactions of many pairs of molecules and is much faster than more rigorous approaches. However, its application requires carefully preprocessing each molecule and selecting a series of simulation parameters, which is not always done correctly. We show how preprocessing and simulation parameters can positively or negatively impact molecular docking performance. For example, the inclusion of hydrogen atoms leads to better redocking scores, but molecular dynamics simulations must be performed under certain constraints; otherwise, it may worsen performance rather than improve it. This study clarifies the importance and influence of these different parameters in the simulation results. METHODS We analyzed the influence of different parameters on the predictive ability of molecular docking techniques using two software packages: AutoDock Vina and AutoDock-GPU. Thus, 90 receptor-ligand complexes were redocked, evaluating the root mean square deviation (RMSD) between the original position of the ligand (receptor-ligand complex obtained experimentally) and that obtained by the software for every analysis. We investigated the influence of hydrogen atoms (on the receptor and on the receptor-ligand complex), partial charges (QEq, QTPIE, EEM, EEM2015ha, MMFF94, Gasteiger-Marsili, and no charge), search boxes (size and exhaustiveness), ligand characteristics (size and number of torsions), and the use of molecular dynamics (of the receptor or the receptor-ligand complex) before docking analyses.
Collapse
Affiliation(s)
- Pedro Henrique Callil-Soares
- Chemical Engineering Department, Polytechnic School of the University of São Paulo, Av. Lineu Prestes, 580, São Paulo, 05508-000, Brazil
| | - Lilian Caroline Kramer Biasi
- Chemical Engineering Department, Polytechnic School of the University of São Paulo, Av. Lineu Prestes, 580, São Paulo, 05508-000, Brazil.
| | - Pedro de Alcântara Pessoa Filho
- Chemical Engineering Department, Polytechnic School of the University of São Paulo, Av. Lineu Prestes, 580, São Paulo, 05508-000, Brazil
| |
Collapse
|
13
|
Jana AK, Sharawy M, May ER. Non-equilibrium virus particle dynamics: Microsecond MD simulations of the complete Flock House virus capsid under different conditions. J Struct Biol 2023; 215:107964. [PMID: 37105277 PMCID: PMC10205670 DOI: 10.1016/j.jsb.2023.107964] [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/16/2023] [Revised: 03/22/2023] [Accepted: 04/14/2023] [Indexed: 04/29/2023]
Abstract
Flock House virus (FHV) is an animal virus and considered a model system for non-enveloped viruses. It has a small, icosahedral capsid (T=3) and a bipartite positive-sense RNA genome. We present an extensive study of the FHV capsid dynamics from all-atom molecular dynamics simulations of the complete capsid. The simulations explore different biologically relevant conditions (neutral/low pH, with/without RNA in the capsid) using the CHARMM force field. The results show that low pH destabilizes the capsid, causing radial expansion, and RNA stabilizes the capsid. The finding of low pH destabilization is biologically relevant because the capsid is exposed to low pH in the endosome, where conformational changes occur leading to genome release. We also observe structural changes at the fivefold and twofold symmetry axes that likely relate to the externalization of membrane active γ peptides through the fivefold vertex and extrusion of RNA at the twofold axis. Simulations using the Amber force field at neutral pH are also performed and display similar characteristics to the CHARMM simulations.
Collapse
Affiliation(s)
- Asis K Jana
- DepartmentofMolecularandCellBiology, UniversityofConnecticut, Storrs, CT06269-3125, USA; Department of Microbiology and Biotechnology, Sister Nivedita University, New Town, West Bengal 700156, India
| | - Mahmoud Sharawy
- DepartmentofMolecularandCellBiology, UniversityofConnecticut, Storrs, CT06269-3125, USA
| | - Eric R May
- DepartmentofMolecularandCellBiology, UniversityofConnecticut, Storrs, CT06269-3125, USA.
| |
Collapse
|
14
|
Lynch DL, Pavlova A, Fan Z, Gumbart JC. Understanding Virus Structure and Dynamics through Molecular Simulations. J Chem Theory Comput 2023. [PMID: 37192279 DOI: 10.1021/acs.jctc.3c00116] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.
Collapse
Affiliation(s)
- Diane L Lynch
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Anna Pavlova
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Zixing Fan
- Interdisciplinary Bioengineering Graduate Program, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| |
Collapse
|
15
|
Highland CM, Tan A, Ricaña CL, Briggs JAG, Dick RA. Structural insights into HIV-1 polyanion-dependent capsid lattice formation revealed by single particle cryo-EM. Proc Natl Acad Sci U S A 2023; 120:e2220545120. [PMID: 37094124 PMCID: PMC10160977 DOI: 10.1073/pnas.2220545120] [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: 12/02/2022] [Accepted: 03/12/2023] [Indexed: 04/26/2023] Open
Abstract
The HIV-1 capsid houses the viral genome and interacts extensively with host cell proteins throughout the viral life cycle. It is composed of capsid protein (CA), which assembles into a conical fullerene lattice composed of roughly 200 CA hexamers and 12 CA pentamers. Previous structural analyses of individual CA hexamers and pentamers have provided valuable insight into capsid structure and function, but detailed structural information about these assemblies in the broader context of the capsid lattice is lacking. In this study, we combined cryoelectron tomography and single particle analysis (SPA) cryoelectron microscopy to determine structures of continuous regions of the capsid lattice containing both hexamers and pentamers. We also developed a method of liposome scaffold-based in vitro lattice assembly ("lattice templating") that enabled us to directly study the lattice under a wider range of conditions than has previously been possible. Using this approach, we identified a critical role for inositol hexakisphosphate in pentamer formation and determined the structure of the CA lattice bound to the capsid-targeting antiretroviral drug GS-6207 (lenacapavir). Our work reveals key structural details of the mature HIV-1 CA lattice and establishes the combination of lattice templating and SPA as a robust strategy for studying retroviral capsid structure and capsid interactions with host proteins and antiviral compounds.
Collapse
Affiliation(s)
- Carolyn M. Highland
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY14853
- Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY14853
| | - Aaron Tan
- Structural Studies Division, Medical Research Council Laboratory of Molecular Biology, CambridgeCB2 0QH, UK
| | - Clifton L. Ricaña
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY14853
| | - John A. G. Briggs
- Structural Studies Division, Medical Research Council Laboratory of Molecular Biology, CambridgeCB2 0QH, UK
- Department of Cell and Virus Structure, Max Planck Institute of Biochemistry, Munich82512, Germany
| | - Robert A. Dick
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY14853
| |
Collapse
|
16
|
Bryer AJ, Rey JS, Perilla JR. Performance efficient macromolecular mechanics via sub-nanometer shape based coarse graining. Nat Commun 2023; 14:2014. [PMID: 37037809 PMCID: PMC10086035 DOI: 10.1038/s41467-023-37801-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 03/30/2023] [Indexed: 04/12/2023] Open
Abstract
Dimensionality reduction via coarse grain modeling is a valuable tool in biomolecular research. For large assemblies, ultra coarse models are often knowledge-based, relying on a priori information to parameterize models thus hindering general predictive capability. Here, we present substantial advances to the shape based coarse graining (SBCG) method, which we refer to as SBCG2. SBCG2 utilizes a revitalized formulation of the topology representing network which makes high-granularity modeling possible, preserving atomistic details that maintain assembly characteristics. Further, we present a method of granularity selection based on charge density Fourier Shell Correlation and have additionally developed a refinement method to optimize, adjust and validate high-granularity models. We demonstrate our approach with the conical HIV-1 capsid and heteromultimeric cofilin-2 bound actin filaments. Our approach is available in the Visual Molecular Dynamics (VMD) software suite, and employs a CHARMM-compatible Hamiltonian that enables high-performance simulation in the GPU-resident NAMD3 molecular dynamics engine.
Collapse
Affiliation(s)
- Alexander J Bryer
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA
| | - Juan S Rey
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA.
| |
Collapse
|
17
|
Shen Q, Kumari S, Xu C, Jang S, Shi J, Burdick RC, Levintov L, Xiong Q, Wu C, Devarkar SC, Tian T, Tripler TN, Hu Y, Yuan S, Temple J, Feng Q, Lusk CP, Aiken C, Engelman AN, Perilla JR, Pathak VK, Lin C, Xiong Y. The capsid lattice engages a bipartite NUP153 motif to mediate nuclear entry of HIV-1 cores. Proc Natl Acad Sci U S A 2023; 120:e2202815120. [PMID: 36943880 PMCID: PMC10068764 DOI: 10.1073/pnas.2202815120] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 01/30/2023] [Indexed: 03/23/2023] Open
Abstract
Increasing evidence has suggested that the HIV-1 capsid enters the nucleus in a largely assembled, intact form. However, not much is known about how the cone-shaped capsid interacts with the nucleoporins (NUPs) in the nuclear pore for crossing the nuclear pore complex. Here, we elucidate how NUP153 binds HIV-1 capsid by engaging the assembled capsid protein (CA) lattice. A bipartite motif containing both canonical and noncanonical interaction modules was identified at the C-terminal tail region of NUP153. The canonical cargo-targeting phenylalanine-glycine (FG) motif engaged the CA hexamer. By contrast, a previously unidentified triple-arginine (RRR) motif in NUP153 targeted HIV-1 capsid at the CA tri-hexamer interface in the capsid. HIV-1 infection studies indicated that both FG- and RRR-motifs were important for the nuclear import of HIV-1 cores. Moreover, the presence of NUP153 stabilized tubular CA assemblies in vitro. Our results provide molecular-level mechanistic evidence that NUP153 contributes to the entry of the intact capsid into the nucleus.
Collapse
Affiliation(s)
- Qi Shen
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
- Department of Cell Biology, Yale School of Medicine, New Haven, CT06520
- Nanobiology Institute, Yale University, West Haven, CT06516
| | - Sushila Kumari
- HIV Dynamics and Replication Program, Center for Cancer Research, National Cancer Institute at Frederick, Frederick, MD21702
| | - Chaoyi Xu
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE19716
| | - Sooin Jang
- Department of Cancer Immunology and Virology, Dana-Farber Cancer Institute, Boston, MA02215
- Department of Medicine, Harvard Medical School, Boston, MA02115
| | - Jiong Shi
- Department of Pathology, Microbiology and Immunology, Vanderbilt University Medical Center, Nashville, TN37232
| | - Ryan C. Burdick
- HIV Dynamics and Replication Program, Center for Cancer Research, National Cancer Institute at Frederick, Frederick, MD21702
| | - Lev Levintov
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE19716
| | - Qiancheng Xiong
- Department of Cell Biology, Yale School of Medicine, New Haven, CT06520
- Nanobiology Institute, Yale University, West Haven, CT06516
| | - Chunxiang Wu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Swapnil C. Devarkar
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Taoran Tian
- Nanobiology Institute, Yale University, West Haven, CT06516
| | - Therese N. Tripler
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Yingxia Hu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Shuai Yuan
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Joshua Temple
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| | - Qingzhou Feng
- Department of Cell Biology, Yale School of Medicine, New Haven, CT06520
- Nanobiology Institute, Yale University, West Haven, CT06516
| | - C. Patrick Lusk
- Department of Cell Biology, Yale School of Medicine, New Haven, CT06520
| | - Christopher Aiken
- Department of Pathology, Microbiology and Immunology, Vanderbilt University Medical Center, Nashville, TN37232
| | - Alan N. Engelman
- Department of Cancer Immunology and Virology, Dana-Farber Cancer Institute, Boston, MA02215
- Department of Medicine, Harvard Medical School, Boston, MA02115
| | - Juan R. Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE19716
| | - Vinay K. Pathak
- HIV Dynamics and Replication Program, Center for Cancer Research, National Cancer Institute at Frederick, Frederick, MD21702
| | - Chenxiang Lin
- Department of Cell Biology, Yale School of Medicine, New Haven, CT06520
- Nanobiology Institute, Yale University, West Haven, CT06516
- Department of Biomedical Engineering, Yale University, New Haven, CT06511
| | - Yong Xiong
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT06511
| |
Collapse
|
18
|
Pezeshkian W, Grünewald F, Narykov O, Lu S, Arkhipova V, Solodovnikov A, Wassenaar TA, Marrink SJ, Korkin D. Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling. Structure 2023; 31:492-503.e7. [PMID: 36870335 DOI: 10.1016/j.str.2023.02.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 11/15/2022] [Accepted: 02/07/2023] [Indexed: 03/06/2023]
Abstract
Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.
Collapse
Affiliation(s)
- Weria Pezeshkian
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
| | - Fabian Grünewald
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands
| | - Oleksandr Narykov
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | - Senbao Lu
- Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA
| | | | | | - Tsjerk A Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands; Institute for Life Science and Technology, Hanze University of Applied Sciences, 9747AS Groningen, the Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747AG Groningen, the Netherlands.
| | - Dmitry Korkin
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, USA; Bioinformatics and Computational Biology Program, Worcester Polytechnic Institute, Worcester, MA 01609, USA.
| |
Collapse
|
19
|
Gupta M, Pak AJ, Voth GA. Critical mechanistic features of HIV-1 viral capsid assembly. SCIENCE ADVANCES 2023; 9:eadd7434. [PMID: 36608139 PMCID: PMC9821859 DOI: 10.1126/sciadv.add7434 10.1126/sciadv.add7434] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 12/08/2022] [Indexed: 11/04/2023]
Abstract
The maturation of HIV-1 capsid protein (CA) into a cone-shaped lattice capsid is critical for viral infectivity. CA can self-assemble into a range of capsid morphologies made of ~175 to 250 hexamers and 12 pentamers. The cellular polyanion inositol hexakisphosphate (IP6) has recently been demonstrated to facilitate conical capsid formation by coordinating a ring of arginine residues within the central cavity of capsid hexamers and pentamers. However, the kinetic interplay of events during IP6 and CA coassembly is unclear. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanism of capsid formation, including the role played by IP6. We show that IP6, in small quantities at first, promotes curvature generation by trapping pentameric defects in the growing lattice and shifts assembly behavior toward kinetically favored outcomes. Our analysis also suggests that IP6 can stabilize metastable capsid intermediates and can induce structural pleomorphism in mature capsids.
Collapse
Affiliation(s)
- Manish Gupta
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| | | | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| |
Collapse
|
20
|
Gupta M, Pak AJ, Voth GA. Critical mechanistic features of HIV-1 viral capsid assembly. SCIENCE ADVANCES 2023; 9:eadd7434. [PMID: 36608139 PMCID: PMC9821859 DOI: 10.1126/sciadv.add7434+10.1126/sciadv.add7434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 12/08/2022] [Indexed: 01/21/2024]
Abstract
The maturation of HIV-1 capsid protein (CA) into a cone-shaped lattice capsid is critical for viral infectivity. CA can self-assemble into a range of capsid morphologies made of ~175 to 250 hexamers and 12 pentamers. The cellular polyanion inositol hexakisphosphate (IP6) has recently been demonstrated to facilitate conical capsid formation by coordinating a ring of arginine residues within the central cavity of capsid hexamers and pentamers. However, the kinetic interplay of events during IP6 and CA coassembly is unclear. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanism of capsid formation, including the role played by IP6. We show that IP6, in small quantities at first, promotes curvature generation by trapping pentameric defects in the growing lattice and shifts assembly behavior toward kinetically favored outcomes. Our analysis also suggests that IP6 can stabilize metastable capsid intermediates and can induce structural pleomorphism in mature capsids.
Collapse
Affiliation(s)
- Manish Gupta
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| | | | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| |
Collapse
|
21
|
Gupta M, Pak AJ, Voth GA. Critical mechanistic features of HIV-1 viral capsid assembly. SCIENCE ADVANCES 2023; 9:eadd7434. [PMID: 36608139 PMCID: PMC9821859 DOI: 10.1126/sciadv.add7434] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 12/08/2022] [Indexed: 05/29/2023]
Abstract
The maturation of HIV-1 capsid protein (CA) into a cone-shaped lattice capsid is critical for viral infectivity. CA can self-assemble into a range of capsid morphologies made of ~175 to 250 hexamers and 12 pentamers. The cellular polyanion inositol hexakisphosphate (IP6) has recently been demonstrated to facilitate conical capsid formation by coordinating a ring of arginine residues within the central cavity of capsid hexamers and pentamers. However, the kinetic interplay of events during IP6 and CA coassembly is unclear. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanism of capsid formation, including the role played by IP6. We show that IP6, in small quantities at first, promotes curvature generation by trapping pentameric defects in the growing lattice and shifts assembly behavior toward kinetically favored outcomes. Our analysis also suggests that IP6 can stabilize metastable capsid intermediates and can induce structural pleomorphism in mature capsids.
Collapse
|
22
|
Casalino L, Seitz C, Lederhofer J, Tsybovsky Y, Wilson IA, Kanekiyo M, Amaro RE. Breathing and Tilting: Mesoscale Simulations Illuminate Influenza Glycoprotein Vulnerabilities. ACS CENTRAL SCIENCE 2022; 8:1646-1663. [PMID: 36589893 PMCID: PMC9801513 DOI: 10.1021/acscentsci.2c00981] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/19/2022] [Indexed: 05/28/2023]
Abstract
Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from convalescent human donor, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.
Collapse
Affiliation(s)
- Lorenzo Casalino
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California92093, United States
| | - Christian Seitz
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California92093, United States
| | - Julia Lederhofer
- Vaccine
Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, Maryland20892, United States
| | - Yaroslav Tsybovsky
- Electron
Microscopy Laboratory, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research
Sponsored by the National Cancer Institute, Frederick, Maryland21702, United States
| | - Ian A. Wilson
- Department
of Integrative Structural and Computational Biology and the Skaggs
Institute for Chemical Biology, The Scripps
Research Institute, La Jolla, California92037, United States
| | - Masaru Kanekiyo
- Vaccine
Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, Maryland20892, United States
| | - Rommie E. Amaro
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
23
|
Pathak AK, Bandyopadhyay T. Heat-induced transitions of an empty minute virus of mice capsid in explicit water: all-atom MD simulation. J Biomol Struct Dyn 2022; 40:11900-11913. [PMID: 34459706 DOI: 10.1080/07391102.2021.1969283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The capsid-like structure of the virus-based protein nanoparticles (NPs) can serve as bionanomaterials, with applications in biomedicines and nanotechnology. Release of packaged material from these nanocontainers is associated with subtle conformational changes of the NP structure, which in vitro, is readily accomplished by heating. Characterizing the structural changes as a function of temperature may provide fresh insights into nanomaterial/antiviral strategies. Here, we have calculated heat induced changes in the properties of an empty minute virus of mice particle using large-scale ≈ 3.0 × 106 all-atom molecular dynamics simulations. We focus on two heat induced structural changes of the NP, namely, dynamical transition (DT) and breathing transition (BT), both characterized by sudden and sharp change of measured parameters at temperatures, TDT and TBT, respectively. While DT is assessed by mean-square fluctuation of hydrogen atoms of the NP, BT is monitored through internal volume and permeation rate of water molecules through the NP. Both the transitions, resulting primarily from collective atomistic motion, are found to occur at temperatures widely separated from one another (TBT>TDT). The breathing motions, responsible for the translocation events of the packaged materials through the NP to kick off, are further probed by computing atomic resolution stresses from NVE simulations. Distribution of equilibrium atomistic stresses on the NP reveals a largely asymmetric nature and suggests structural breathing may actually represent large dynamic changes in the hotspot regions, far from the NP pores, which is in remarkable resemblance with recently conducted hydrogen-deuterium exchange coupled to mass spectrometry experiment. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Arup Kumar Pathak
- Theoretical Chemistry Section, Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India.,Homi Bhabha National Institute, Mumbai, India
| | - Tusar Bandyopadhyay
- Theoretical Chemistry Section, Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India.,Homi Bhabha National Institute, Mumbai, India
| |
Collapse
|
24
|
Kumar P, Bhardwaj T, Kumar A, Garg N, Giri R. One microsecond MD simulations of the SARS-CoV-2 main protease and hydroxychloroquine complex reveal the intricate nature of binding. J Biomol Struct Dyn 2022; 40:10763-10770. [PMID: 34320905 DOI: 10.1080/07391102.2021.1948447] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Currently, several vaccines and antivirals across the globe are in clinical trials. Hydroxychloroquine (HCQ) was reported to inhibit the SARS-CoV-2 virus in antiviral assays. Here, it raises the curiosity about the molecular target of HCQ inside the cell. It may inhibit some of the viral targets, or some other complex mechanisms must be at disposal towards action mechanisms. In some of the viruses, proteases are experimentally reported to be a potential target of HCQ. However, no in-depth investigations are available in the literature yet. Henceforth, we have carried out extensive, one-microsecond long molecular dynamics simulations of the bound complex of hydroxychloroquine with main protease (Mpro) of SARS-CoV-2. Our analysis found that HCQ binds within the catalytic pocket of Mpro and remains stable upto one-third of simulation time but further causes increased fluctuations in simulation parameters. In the end, the HCQ does not possess any pre-formed hydrogen bond, other non-covalent interactions with Mpro, ultimately showing the unsteadiness in binding at catalytic binding pocket and may suggest that HCQ may not inhibit the Mpro. In the future, this study would require experimental validation on enzyme assays against Mpro, and that may be the final say. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Prateek Kumar
- School of Basic Sciences, VPO Kamand, Indian Institute of Technology, Mandi, Mandi, India
| | - Taniya Bhardwaj
- School of Basic Sciences, VPO Kamand, Indian Institute of Technology, Mandi, Mandi, India
| | - Ankur Kumar
- School of Basic Sciences, VPO Kamand, Indian Institute of Technology, Mandi, Mandi, India
| | - Neha Garg
- Faculty of Ayurveda, Department of Medicinal Chemistry, Institute of Medical Sciences, Banaras Hindu University, Varanasi, India
| | - Rajanish Giri
- School of Basic Sciences, VPO Kamand, Indian Institute of Technology, Mandi, Mandi, India
| |
Collapse
|
25
|
Apostolidou D, Zhang P, Yang W, Marszalek PE. Mechanical Unfolding and Refolding of NanoLuc via Single-Molecule Force Spectroscopy and Computer Simulations. Biomacromolecules 2022; 23:5164-5178. [PMID: 36350253 DOI: 10.1021/acs.biomac.2c00997] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A highly bioluminescent protein, NanoLuc (Nluc), has seen numerous applications in biological assays since its creation. We recently engineered a NanoLuc polyprotein that showed high bioluminescence but displayed a strong misfolding propensity after mechanical unfolding. Here, we present our single-molecule force spectroscopy (SMFS) studies by atomic force microscopy (AFM) and steered molecular dynamics (SMD) simulations on two new hybrid protein constructs comprised of Nluc and I91 titin domains, I91-I91-Nluc-I91-I91-I91-I91 (I912-Nluc-I914) and I91-Nluc-I91-Nluc-I91-Nluc-I91, to characterize the unfolding behavior of Nluc in detail and to further investigate its misfolding properties that we observed earlier for the I912-Nluc3-I912 construct. Our SMFS results confirm that Nluc's unfolding proceeds similarly in all constructs; however, Nluc's refolding differs in these constructs, and its misfolding is minimized when Nluc is monomeric or separated by I91 domains. Our simulations on monomeric Nluc, Nluc dyads, and Nluc triads pinpointed the origin of its mechanical stability and captured interesting unfolding intermediates, which we also observed experimentally.
Collapse
Affiliation(s)
- Dimitra Apostolidou
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina27708, United States
| | - Pan Zhang
- Department of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina27708, United States
| | - Piotr E Marszalek
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina27708, United States
| |
Collapse
|
26
|
Kumar G, Duggisetty SC, Srivastava A. A Review of Mechanics-Based Mesoscopic Membrane Remodeling Methods: Capturing Both the Physics and the Chemical Diversity. J Membr Biol 2022; 255:757-777. [PMID: 36197492 DOI: 10.1007/s00232-022-00268-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 08/29/2022] [Indexed: 12/24/2022]
Abstract
Specialized classes of proteins, working together in a tightly orchestrated manner, induce and maintain highly curved cellular and organelles membrane morphology. Due to the various experimental constraints, including the resolution limits of imaging techniques, it is non-trivial to accurately elucidate interactions among the various components involved in membrane deformation. The spatial and temporal scales of the systems also make it formidable to investigate them using simulations with molecular details. Interestingly, mechanics-based mesoscopic models have been used with great success in recapitulating the membrane deformations observed in experiments. In this review, we collate together and discuss the various mechanics-based mesoscopic models for protein-mediated membrane deformation studies. In particular, we provide an elaborate description of a mesoscopic model where the membrane is modeled as a triangulated sheet and proteins are represented as either nematics or filaments. This representation allows us to explore the various aspects of protein-protein and protein-membrane interactions as well as examine the underlying mechanistic pathways for emergent behavior such as curvature-mediated protein localization and membrane deformation. We also put forward current efforts in the field towards back-mapping these mesoscopic models to finer-grained particle-based models-a framework that could be used to explore how molecular interactions propagate to physical scales and vice-versa. We end the review with an integrative-modeling-based road map where experimental imaging micrograph and biochemical data are combined with mesoscopic and molecular simulations methods in a theoretically consistent manner to faithfully recapitulate the multiple length and time scales in the membrane remodeling processes.
Collapse
Affiliation(s)
- Gaurav Kumar
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, C. V. Raman Road, Bangalore, Karnataka, 560012, India
| | - Satya Chaithanya Duggisetty
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, C. V. Raman Road, Bangalore, Karnataka, 560012, India
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, C. V. Raman Road, Bangalore, Karnataka, 560012, India.
| |
Collapse
|
27
|
Hendam A, Al-Sadek AF, Hefny HA. Molecular Dynamic Simulation of Neurexin1α Mutations Associated with Mental Disorder. J Mol Neurosci 2022; 72:2252-2272. [PMID: 36197641 PMCID: PMC9532826 DOI: 10.1007/s12031-022-02072-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 09/22/2022] [Indexed: 12/14/2022]
Abstract
Neurexin1 gene is essential for formulating synaptic cell adhesion to establish synapses. In a previous work, 38 SNPs in Neurexin1 recoded in mental disorder patients have been collected. Five computational prediction tools have been used to predict the effect of SNPs on protein function and stability. Only four SNPs in Neurexin1α have deleterious prediction results from at least four tools. The current work aims to use molecular dynamic simulation (MD) to study the effects of the four mutations on Neurexin1α both on the whole protein as well as identifying affected domains by mutations. A protein model that consists of five domains out of six domains in the real protein was used; missing residues were added, and model was tested for quality. The MD experiment has last for 1.5 μs where four parameters have been used for studying the whole protein in addition to three more parameters for the domain analysis. The whole protein study has shown that two mutations E427I for Autism and R525C for non-syndromic intellectual disability (NSID) have distinctive behavior across the four used parameters. Domain study has confirmed the previous results where the five domains of R525C have acted differently from wild type (WT), while E427I has acted differently for four domains from wild type. The other two mutations D104H and G379E have three domains that only acted differently from wild type. The fourth domain of all mutations has an obvious distinctive behavior from wild type. Further study of E427I and R525C mutations can lead to better understanding of autism and NSID.
Collapse
Affiliation(s)
| | | | - Hesham Ahmed Hefny
- Faculty of Graduate Studies for Statistical Research, Cairo University, Giza, Egypt
| |
Collapse
|
28
|
Kulke M, Vermaas JV. Reversible Unwrapping Algorithm for Constant-Pressure Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:6161-6171. [PMID: 36129782 DOI: 10.1021/acs.jctc.2c00327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Molecular simulation technologies have afforded researchers a unique look into the nanoscale interactions driving physical processes. However, a limitation for molecular dynamics (MD) simulations is that they must be performed on finite-sized systems in order to map onto computational resources. To minimize artifacts arising from finite-sized simulation systems, it is common practice for MD simulations to be performed with periodic boundary conditions (PBCs). However, in order to calculate specific physical properties, such as mean square displacements to calculate diffusion coefficients, continuous particle trajectories where the atomic movements are continuous and do not jump between cell faces are required. In these cases, modifying atomic coordinates through unwrapping schemes is an essential post-processing tool to remove these jumps. Here, two established trajectory unwrapping schemes are applied to 1 μs wrapped trajectories for a small water box and lysozyme in water. The existing schemes can result in spurious diffusion coefficients, long bonds within unwrapped molecules, and inconsistent atomic coordinates when coordinates are rewrapped after unwrapping. We determine that prior unwrapping schemes do not account for changing periodic box dimensions and introduce an additional correction term to the existing displacement unwrapping scheme to correct for these artifacts. We also demonstrate that the resulting algorithm is a hybrid between the existing heuristic and displacement unwrapping schemes. After treatment using this new unwrapping scheme, molecular geometries are correct even after long simulations. In anticipation for longer MD trajectories, we develop implementations for this new scheme in multiple PBC handling tools.
Collapse
Affiliation(s)
- Martin Kulke
- MSU-DOE Plant Research Laboratory and Department of Biochemistry and Molecular Biology, Michigan State University, 612 Wilson Rd, East Lansing, Michigan 48824, United States
| | - Josh V Vermaas
- MSU-DOE Plant Research Laboratory and Department of Biochemistry and Molecular Biology, Michigan State University, 612 Wilson Rd, East Lansing, Michigan 48824, United States
| |
Collapse
|
29
|
Wang Q, Kuci D, Bhattacharya S, Hadden‐Perilla JA, Gupta R. Dynamic regulation of Zn(II) sequestration by calgranulin C. Protein Sci 2022; 31:e4403. [PMID: 36367084 PMCID: PMC9650546 DOI: 10.1002/pro.4403] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 07/09/2022] [Accepted: 07/17/2022] [Indexed: 12/13/2022]
Abstract
Calgranulin C performs antimicrobial activity in the human immune response by sequestering Zn(II). This biological function is afforded with the aid of two structurally distinct Ca(II)-binding EF hand motifs, wherein one of which bears an unusual amino acid sequence. Here, we utilize solution state NMR relaxation measurements to investigate the mechanism of Ca(II)-modulated enhancement of Zn(II) sequestration by calgranulin C. Using C13 /N15 CPMG dispersion experiments we have measured pH-dependent major and minor state populations exchanging on micro-to-millisecond timescale. This conformational exchange takes place exclusively in the Ca(II)-bound state and can be mapped to residues located in the EF-I loop and the linker between the tandem EF hands. Molecular dynamics (MD) simulations spanning nano-to-microsecond timescale offer insights into the role of pH-dependent electrostatic interactions in EF-hand dynamics. Our results suggest a pH-regulated dynamic equilibrium of conformations that explore a range of "closed" and partially "open" sidechain configurations within the Zn(II) binding site. We propose a novel mechanism by which Ca(II) binding to a non-canonical EF loop regulates its flexibility and tunes the antimicrobial activity of calgranulin C.
Collapse
Affiliation(s)
- Qian Wang
- Department of Chemistry, College of Staten IslandCity University of New YorkNew YorkNew YorkUSA
| | - Deniz Kuci
- Department of Chemistry, College of Staten IslandCity University of New YorkNew YorkNew YorkUSA
| | | | | | - Rupal Gupta
- Department of Chemistry, College of Staten IslandCity University of New YorkNew YorkNew YorkUSA
- Ph.D. Programs in Biochemistry and ChemistryThe Graduate Center of the City University of New YorkNew YorkNew YorkUSA
| |
Collapse
|
30
|
Integrative modeling of the cell. Acta Biochim Biophys Sin (Shanghai) 2022; 54:1213-1221. [PMID: 36017893 PMCID: PMC9909318 DOI: 10.3724/abbs.2022115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
A whole-cell model represents certain aspects of the cell structure and/or function. Due to the high complexity of the cell, an integrative modeling approach is often taken to utilize all available information including experimental data, prior knowledge and prior models. In this review, we summarize an emerging workflow of whole-cell modeling into five steps: (i) gather information; (ii) represent the modeled system into modules; (iii) translate input information into scoring function; (iv) sample the whole-cell model; (v) validate and interpret the model. In particular, we propose the integrative modeling of the cell by combining available (whole-cell) models to maximize the accuracy, precision, and completeness. In addition, we list quantitative predictions of various aspects of cell biology from existing whole-cell models. Moreover, we discuss the remaining challenges and future directions, and highlight the opportunity to establish an integrative spatiotemporal multi-scale whole-cell model based on a community approach.
Collapse
|
31
|
Hubman A, Urbic T. Structure and thermodynamics of a 2D Lennard-Jones hexagonal fluid. MOLECULAR SIMULATION 2022; 48:1435-1444. [PMID: 37727614 PMCID: PMC10508885 DOI: 10.1080/08927022.2022.2096219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Accepted: 06/20/2022] [Indexed: 10/17/2022]
Abstract
The thermodynamic and structural properties of the 2D hexagonal soft-sites fluid are examined by integral equation theory benchmarked against extensive Monte Carlo simulations. Hexamers are built of six equal Lennard-Jones segments. Site-site integral equation theory is used to compute site-site correlation functions, excess internal energies and isotherms over a wide range of conditions and compared with results obtained from Monte Carlo simulations. Various approaches for computing the pressure are discussed as well. Satisfactory qualitative agreement between theory and simulations is found with details depending on the applied closure relation.
Collapse
Affiliation(s)
- Anže Hubman
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
| |
Collapse
|
32
|
Tasneem N, Szyszka TN, Jenner EN, Lau YH. How Pore Architecture Regulates the Function of Nanoscale Protein Compartments. ACS NANO 2022; 16:8540-8556. [PMID: 35583458 DOI: 10.1021/acsnano.2c02178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Self-assembling proteins can form porous compartments that adopt well-defined architectures at the nanoscale. In nature, protein compartments act as semipermeable barriers to enable spatial separation and organization of complex biochemical processes. The compartment pores play a key role in their overall function by selectively controlling the influx and efflux of important biomolecular species. By engineering the pores, the functionality of compartments can be tuned to facilitate non-native applications, such as artificial nanoreactors for catalysis. In this review, we analyze how protein structure determines the porosity and impacts the function of both native and engineered compartments, highlighting the wealth of structural data recently obtained by cryo-EM and X-ray crystallography. Through this analysis, we offer perspectives on how current structural insights can inform future studies into the design of artificial protein compartments as nanoreactors with tunable porosity and function.
Collapse
Affiliation(s)
- Nuren Tasneem
- School of Chemistry, The University of Sydney, Eastern Avenue, Camperdown, New South Wales 2006, Australia
| | - Taylor N Szyszka
- School of Chemistry, The University of Sydney, Eastern Avenue, Camperdown, New South Wales 2006, Australia
- University of Sydney Nano Institute, Camperdown, New South Wales 2006, Australia
| | - Eric N Jenner
- School of Chemistry, The University of Sydney, Eastern Avenue, Camperdown, New South Wales 2006, Australia
| | - Yu Heng Lau
- School of Chemistry, The University of Sydney, Eastern Avenue, Camperdown, New South Wales 2006, Australia
- University of Sydney Nano Institute, Camperdown, New South Wales 2006, Australia
| |
Collapse
|
33
|
Nicastro G, Lucci M, Oregioni A, Kelly G, Frenkiel TA, Taylor IA. CP-MAS and solution NMR studies of allosteric communication in CA-assemblies of HIV-1. J Mol Biol 2022; 434:167691. [PMID: 35738429 DOI: 10.1016/j.jmb.2022.167691] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 06/14/2022] [Accepted: 06/15/2022] [Indexed: 10/18/2022]
Abstract
Solution and solid-state NMR spectroscopy are highly complementary techniques for studying structure and dynamics in very high molecular weight systems. Here we have analysed the dynamics of HIV-1 capsid (CA) assemblies in presence of the cofactors IP6 and ATPγS and the host-factor CPSF6 using a combination of solution state and cross polarisation magic angle spinning (CP-MAS) solid-state NMR. In particular, dynamical effects on ns to µs and µs to ms timescales are observed revealing diverse motions in assembled CA. Using CP-MAS NMR, we exploited the sensitivity of the amide/Cα-Cβ backbone chemical shifts in DARR and NCA spectra to observe the plasticity of the HIV-1 CA tubular assemblies and also map the binding of cofactors and the dynamics of cofactor-CA complexes. In solution, we measured how the addition of host- and co-factors to CA -hexamers perturbed the chemical shifts and relaxation properties of CA-Ile and -Met methyl groups using transverse-relaxation-optimized NMR spectroscopy to exploit the sensitivity of methyl groups as probes in high-molecular weight proteins. These data show how dynamics of the CA protein assembly over a range of spatial and temporal scales play a critical role in CA function. Moreover, we show that binding of IP6, ATPγS and CPSF6 results in local chemical shift as well as dynamic changes for a significant, contiguous portion of CA, highlighting how allosteric pathways communicate ligand interactions between adjacent CA protomers.
Collapse
Affiliation(s)
- Giuseppe Nicastro
- Macromolecular Structure Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK
| | - Massimo Lucci
- CIRMMP, University of Florence, Via L. Sacconi, 6 50019 Sesto Fiorentino (FI), Italy
| | - Alain Oregioni
- The Medical Research Council Biomedical NMR Centre, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK
| | - Geoff Kelly
- The Medical Research Council Biomedical NMR Centre, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK
| | - Tom A Frenkiel
- The Medical Research Council Biomedical NMR Centre, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK
| | - Ian A Taylor
- Macromolecular Structure Laboratory, The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK.
| |
Collapse
|
34
|
Bayarri G, Andrio P, Hospital A, Orozco M, Gelpí JL. BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations. Nucleic Acids Res 2022; 50:W99-W107. [PMID: 35639735 PMCID: PMC9252775 DOI: 10.1093/nar/gkac380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 04/19/2022] [Accepted: 05/02/2022] [Indexed: 11/15/2022] Open
Abstract
We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users’ own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users’ projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.
Collapse
Affiliation(s)
- Genís Bayarri
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Pau Andrio
- Barcelona Supercomputing Center (BSC), Jordi Girona 29, 08034, Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona). The Barcelona Institute of Science and Technology. Baldiri Reixac 10-12, 08028 Barcelona, Spain.,Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| | - Josep Lluís Gelpí
- Barcelona Supercomputing Center (BSC), Jordi Girona 29, 08034, Barcelona, Spain.,Department of Biochemistry and Molecular Biology, University of Barcelona, 08028 Barcelona, Spain
| |
Collapse
|
35
|
Yu A, Lee EMY, Briggs JAG, Ganser-Pornillos BK, Pornillos O, Voth GA. Strain and rupture of HIV-1 capsids during uncoating. Proc Natl Acad Sci U S A 2022; 119:e2117781119. [PMID: 35238630 PMCID: PMC8915963 DOI: 10.1073/pnas.2117781119] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 01/25/2022] [Indexed: 12/18/2022] Open
Abstract
SignificanceThe mature capsids of HIV-1 are transiently stable complexes that self-assemble around the viral genome during maturation, and uncoat to release preintegration complexes that archive a double-stranded DNA copy of the virus in the host cell genome. However, a detailed view of how HIV cores rupture remains lacking. Here, we elucidate the physical properties involved in capsid rupture using a combination of large-scale all-atom molecular dynamics simulations and cryo-electron tomography. We find that intrinsic strain on the capsid forms highly correlated patterns along the capsid surface, along which cracks propagate. Capsid rigidity also increases with high strain. Our findings provide fundamental insight into viral capsid uncoating.
Collapse
Affiliation(s)
- Alvin Yu
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Elizabeth M. Y. Lee
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, IL 60637
| | - John A. G. Briggs
- Department of Cell and Virus Structure, Max Planck Institute of Biochemistry, 82152 Martinsried, Germany
| | - Barbie K. Ganser-Pornillos
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA 22903
| | - Owen Pornillos
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA 22903
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
36
|
Gupta C, Sarkar D, Tieleman DP, Singharoy A. The ugly, bad, and good stories of large-scale biomolecular simulations. Curr Opin Struct Biol 2022; 73:102338. [PMID: 35245737 DOI: 10.1016/j.sbi.2022.102338] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 11/29/2021] [Accepted: 01/24/2022] [Indexed: 12/20/2022]
Abstract
Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.
Collapse
Affiliation(s)
- Chitrak Gupta
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA. https://twitter.com/ChitrakGupta2
| | - Daipayan Sarkar
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; MSU-DOE Plant Research Laboratory, Michigan State University, East Lansing, MI, 48824-1319, USA. https://twitter.com/17Dsarkar
| | - D Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada.
| | - Abhishek Singharoy
- School of Molecular Sciences, Center for Applied Structural Discovery, Arizona State University at Tempe, Tempe, AZ, 85282, USA; Biodesign Institute, Tempe, AZ, 85281, USA.
| |
Collapse
|
37
|
Exploring cryo-electron microscopy with molecular dynamics. Biochem Soc Trans 2022; 50:569-581. [PMID: 35212361 DOI: 10.1042/bst20210485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 01/31/2022] [Accepted: 02/02/2022] [Indexed: 11/17/2022]
Abstract
Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.
Collapse
|
38
|
Wong H, Crowet JM, Dauchez M, Ricard-Blum S, Baud S, Belloy N. Multiscale modelling of the extracellular matrix. Matrix Biol Plus 2022; 13:100096. [PMID: 35072037 PMCID: PMC8763633 DOI: 10.1016/j.mbplus.2021.100096] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2021] [Revised: 12/06/2021] [Accepted: 12/07/2021] [Indexed: 12/29/2022] Open
Abstract
The extracellular matrix is a complex three-dimensional network of molecules that provides cells with a complex microenvironment. The major constituents of the extracellular matrix such as collagen, elastin and associated proteins form supramolecular assemblies contributing to its physicochemical properties and organization. The structure of proteins and their supramolecular assemblies such as fibrils have been studied at the atomic level (e.g., by X-ray crystallography, Nuclear Magnetic Resonance and cryo-Electron Microscopy) or at the microscopic scale. However, many protein complexes are too large to be studied at the atomic level and too small to be studied by microscopy. Most extracellular matrix components fall into this intermediate scale, so-called the mesoscopic scale, preventing their detailed characterization. Simulation and modelling are some of the few powerful and promising approaches that can deepen our understanding of mesoscale systems. We have developed a set of modelling tools to study the self-organization of the extracellular matrix and large motion of macromolecules at the mesoscale level by taking advantage of the dynamics of articulated rigid bodies as a mean to study a larger range of motions at the cost of atomic resolution.
Collapse
Key Words
- Basement membrane
- CG, coarse-grained
- Cryo-EM, cryogenic electron microscopy
- DOF, degrees of freedom
- ECM, extracellular matrix
- EGF, epidermal growth factor
- Extracellular matrix
- FEM, finite element method
- MD, molecular dynamics
- Mesoscopic scale
- Modelling
- NC, non-collagenous
- NMR, nuclear magnetic resonance
- Rigid bodies
- SAXS, small-angle X-ray scattering
- Simulation
Collapse
Affiliation(s)
- Hua Wong
- Université de Reims Champagne Ardenne, CNRS, MEDyC UMR 7369, 51097 Reims, France
| | - Jean-Marc Crowet
- Université de Reims Champagne Ardenne, CNRS, MEDyC UMR 7369, 51097 Reims, France
| | - Manuel Dauchez
- Université de Reims Champagne Ardenne, CNRS, MEDyC UMR 7369, 51097 Reims, France
| | - Sylvie Ricard-Blum
- Univ. Lyon, University Claude Bernard Lyon 1, ICBMS, UMR 5246 CNRS, 69622 Villeurbanne Cedex, France
| | - Stéphanie Baud
- Université de Reims Champagne Ardenne, CNRS, MEDyC UMR 7369, 51097 Reims, France
- Université de Reims Champagne Ardenne, Plateau de Modélisation Moléculaire Multi-Echelle (P3M), Maison de la Simulation de Champagne-Ardenne (MaSCA), 51097 Reims, France
| | - Nicolas Belloy
- Université de Reims Champagne Ardenne, CNRS, MEDyC UMR 7369, 51097 Reims, France
- Université de Reims Champagne Ardenne, Plateau de Modélisation Moléculaire Multi-Echelle (P3M), Maison de la Simulation de Champagne-Ardenne (MaSCA), 51097 Reims, France
| |
Collapse
|
39
|
Nie T, Meng F, Lu F, Sun J, Bie X, Lu Z, Lu Y. Molecular dynamics insight of novel Enzybiotic Salmcide-p1 lysis peptidoglycan to inhibit Salmonella Typhimurium. Food Control 2022. [DOI: 10.1016/j.foodcont.2021.108564] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
40
|
Russo JD, Zhang S, Leung JMG, Bogetti AT, Thompson JP, DeGrave AJ, Torrillo PA, Pratt AJ, Wong KF, Xia J, Copperman J, Adelman JL, Zwier MC, LeBard DN, Zuckerman DM, Chong LT. WESTPA 2.0: High-Performance Upgrades for Weighted Ensemble Simulations and Analysis of Longer-Timescale Applications. J Chem Theory Comput 2022; 18:638-649. [PMID: 35043623 DOI: 10.1021/acs.jctc.1c01154] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The weighted ensemble (WE) family of methods is one of several statistical mechanics-based path sampling strategies that can provide estimates of key observables (rate constants and pathways) using a fraction of the time required by direct simulation methods such as molecular dynamics or discrete-state stochastic algorithms. WE methods oversee numerous parallel trajectories using intermittent overhead operations at fixed time intervals, enabling facile interoperability with any dynamics engine. Here, we report on the major upgrades to the WESTPA software package, an open-source, high-performance framework that implements both basic and recently developed WE methods. These upgrades offer substantial improvements over traditional WE methods. The key features of the new WESTPA 2.0 software enhance the efficiency and ease of use: an adaptive binning scheme for more efficient surmounting of large free energy barriers, streamlined handling of large simulation data sets, exponentially improved analysis of kinetics, and developer-friendly tools for creating new WE methods, including a Python API and resampler module for implementing both binned and "binless" WE strategies.
Collapse
Affiliation(s)
- John D Russo
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - She Zhang
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Jeremy M G Leung
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Anthony T Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Jeff P Thompson
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Alex J DeGrave
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Paul A Torrillo
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - A J Pratt
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kim F Wong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Junchao Xia
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Jeremy Copperman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - Joshua L Adelman
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Matthew C Zwier
- Department of Chemistry, Drake University, Des Moines, Iowa 50311-4505, United States
| | - David N LeBard
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - Lillian T Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
41
|
de Oliveira AB, Contessoto VG, Hassan A, Byju S, Wang A, Wang Y, Dodero‐Rojas E, Mohanty U, Noel JK, Onuchic JN, Whitford PC. SMOG 2 and OpenSMOG: Extending the limits of structure-based models. Protein Sci 2022; 31:158-172. [PMID: 34655449 PMCID: PMC8740843 DOI: 10.1002/pro.4209] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 10/12/2021] [Accepted: 10/12/2021] [Indexed: 01/03/2023]
Abstract
Applying simulations with structure-based G o ¯ - like models has proven to be an effective strategy for investigating the factors that control biomolecular dynamics. The common element of these models is that some (or all) of the intra/inter-molecular interactions are explicitly defined to stabilize an experimentally determined structure. To facilitate the development and application of this broad class of models, we previously released the SMOG 2 software package. This suite allows one to easily customize and distribute structure-based (i.e., SMOG) models for any type of polymer-ligand system. The force fields generated by SMOG 2 may then be used to perform simulations in highly optimized MD packages, such as Gromacs, NAMD, LAMMPS, and OpenMM. Here, we describe extensions to the software and demonstrate the capabilities of the most recent version (SMOG v2.4.2). Changes include new tools that aid user-defined customization of force fields, as well as an interface with the OpenMM simulation libraries (OpenSMOG v1.1.0). The OpenSMOG module allows for arbitrary user-defined contact potentials and non-bonded potentials to be employed in SMOG models, without source-code modifications. To illustrate the utility of these advances, we present applications to systems with millions of atoms, long polymers and explicit ions, as well as models that include non-structure-based (e.g., AMBER-based) energetic terms. Examples include large-scale rearrangements of the SARS-CoV-2 Spike protein, the HIV-1 capsid with explicit ions, and crystallographic lattices of ribosomes and proteins. In summary, SMOG 2 and OpenSMOG provide robust support for researchers who seek to develop and apply structure-based models to large and/or intricate biomolecular systems.
Collapse
Affiliation(s)
| | | | - Asem Hassan
- Department of PhysicsNortheastern University, Dana Research CenterBostonMassachusettsUSA
- Center for Theoretical Biological PhysicsNortheastern UniversityBostonMassachusettsUSA
| | - Sandra Byju
- Department of PhysicsNortheastern University, Dana Research CenterBostonMassachusettsUSA
- Center for Theoretical Biological PhysicsNortheastern UniversityBostonMassachusettsUSA
| | - Ailun Wang
- Center for Theoretical Biological PhysicsNortheastern UniversityBostonMassachusettsUSA
| | - Yang Wang
- Department of ChemistryBoston CollegeChestnut HillMassachusettsUSA
| | | | - Udayan Mohanty
- Department of ChemistryBoston CollegeChestnut HillMassachusettsUSA
| | - Jeffrey K. Noel
- CrystallographyMax Delbrück Center for Molecular MedicineBerlinGermany
- Present address:
Electric Ant Lab, Science Park 106AmsterdamThe Netherlands
| | - Jose N. Onuchic
- Center for Theoretical Biological PhysicsRice UniversityHoustonTexasUSA
- Department of Physics & AstronomyRice UniversityHoustonTexasUSA
- Department of ChemistryRice UniversityHoustonTexasUSA
- Department of BiosciencesRice UniversityHoustonTexasUSA
| | - Paul C. Whitford
- Department of PhysicsNortheastern University, Dana Research CenterBostonMassachusettsUSA
- Center for Theoretical Biological PhysicsNortheastern UniversityBostonMassachusettsUSA
| |
Collapse
|
42
|
Ni T, Zhu Y, Yang Z, Xu C, Chaban Y, Nesterova T, Ning J, Böcking T, Parker MW, Monnie C, Ahn J, Perilla JR, Zhang P. Structure of native HIV-1 cores and their interactions with IP6 and CypA. SCIENCE ADVANCES 2021; 7:eabj5715. [PMID: 34797722 PMCID: PMC8604400 DOI: 10.1126/sciadv.abj5715] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 10/01/2021] [Indexed: 05/24/2023]
Abstract
The viral capsid plays essential roles in HIV replication and is a major platform engaging host factors. To overcome challenges in study native capsid structure, we used the perfringolysin O to perforate the membrane of HIV-1 particles, thus allowing host proteins and small molecules to access the native capsid while improving cryo–electron microscopy image quality. Using cryo–electron tomography and subtomogram averaging, we determined the structures of native capsomers in the presence and absence of inositol hexakisphosphate (IP6) and cyclophilin A and constructed an all-atom model of a complete HIV-1 capsid. Our structures reveal two IP6 binding sites and modes of cyclophilin A interactions. Free energy calculations substantiate the two binding sites at R18 and K25 and further show a prohibitive energy barrier for IP6 to pass through the pentamer. Our results demonstrate that perfringolysin O perforation is a valuable tool for structural analyses of enveloped virus capsids and interactions with host cell factors.
Collapse
Affiliation(s)
- Tao Ni
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Yanan Zhu
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Zhengyi Yang
- Electron Bio-Imaging Centre, Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, UK
| | - Chaoyi Xu
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA
| | - Yuriy Chaban
- Electron Bio-Imaging Centre, Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, UK
| | - Tanya Nesterova
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA
| | - Jiying Ning
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Till Böcking
- EMBL Australia Node in Single Molecule Science and ARC Centre of Excellence in Advanced Molecular Imaging, School of Medical Sciences, UNSW, Sydney, Australia
| | - Michael W. Parker
- Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, Victoria 3010, Australia
- St. Vincent’s Institute of Medical Research, Fitzroy, Victoria 3065, Australia
| | - Christina Monnie
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Jinwoo Ahn
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Juan R. Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, USA
| | - Peijun Zhang
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Electron Bio-Imaging Centre, Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, UK
- Chinese Academy of Medical Sciences Oxford Institute, University of Oxford, Oxford OX3 7BN, UK
| |
Collapse
|
43
|
Ingram Z, Fischer DK, Ambrose Z. Disassembling the Nature of Capsid: Biochemical, Genetic, and Imaging Approaches to Assess HIV-1 Capsid Functions. Viruses 2021; 13:v13112237. [PMID: 34835043 PMCID: PMC8618418 DOI: 10.3390/v13112237] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 10/31/2021] [Accepted: 11/01/2021] [Indexed: 12/20/2022] Open
Abstract
The human immunodeficiency virus type 1 (HIV-1) capsid and its disassembly, or capsid uncoating, has remained an active area of study over the past several decades. Our understanding of the HIV-1 capsid as solely a protective shell has since shifted with discoveries linking it to other complex replication events. The interplay of the HIV-1 capsid with reverse transcription, nuclear import, and integration has led to an expansion of knowledge of capsid functionality. Coincident with advances in microscopy, cell, and biochemistry assays, several models of capsid disassembly have been proposed, in which it occurs in either the cytoplasmic, nuclear envelope, or nuclear regions of the cell. Here, we discuss how the understanding of the HIV-1 capsid has evolved and the key methods that made these discoveries possible.
Collapse
Affiliation(s)
- Zachary Ingram
- Department of Microbiology and Molecular Genetics, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15219, USA; (Z.I.); (D.K.F.)
- Pittsburgh Center for HIV Protein Interactions, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | - Douglas K. Fischer
- Department of Microbiology and Molecular Genetics, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15219, USA; (Z.I.); (D.K.F.)
- Pittsburgh Center for HIV Protein Interactions, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | - Zandrea Ambrose
- Department of Microbiology and Molecular Genetics, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15219, USA; (Z.I.); (D.K.F.)
- Pittsburgh Center for HIV Protein Interactions, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA
- Correspondence:
| |
Collapse
|
44
|
Perišić O, Wriggers W. Mechanism for the Unfolding of the TOP7 Protein in Steered Molecular Dynamics Simulations as Revealed by Mutual Information Analysis. Front Mol Biosci 2021; 8:696609. [PMID: 34660691 PMCID: PMC8516001 DOI: 10.3389/fmolb.2021.696609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2021] [Accepted: 08/30/2021] [Indexed: 12/03/2022] Open
Abstract
We employed mutual information (MI) analysis to detect motions affecting the mechanical resistance of the human-engineered protein Top7. The results are based on the MI analysis of pair contact correlations measured in steered molecular dynamics (SMD) trajectories and their statistical dependence on global unfolding. This study is the first application of the MI analysis to SMD forced unfolding, and we furnish specific SMD recommendations for the utility of parameters and options in the TimeScapes package. The MI analysis provided a global overview of the effect of perturbation on the stability of the protein. We also employed a more conventional trajectory analysis for a detailed description of the mechanical resistance of Top7. Specifically, we investigated 1) the hydropathy of the interactions of structural segments, 2) the H2O concentration near residues relevant for unfolding, and 3) the changing hydrogen bonding patterns and main chain dihedral angles. The results show that the application of MI in the study of protein mechanical resistance can be useful for the engineering of more resistant mutants when combined with conventional analysis. We propose a novel mutation design based on the hydropathy of residues that would stabilize the unfolding region by mimicking its more stable symmetry mate. The proposed design process does not involve the introduction of covalent crosslinks, so it has the potential to preserve the conformational space and unfolding pathway of the protein.
Collapse
Affiliation(s)
| | - Willy Wriggers
- Department of Mechanical and Aerospace Engineering, Old Dominion University, Norfolk, VA, United States
| |
Collapse
|
45
|
Fujimoto K, Yamaguchi Y, Urano R, Shinoda W, Ishikawa T, Omagari K, Tanaka Y, Nakagawa A, Okazaki S. All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution. J Chem Phys 2021; 155:145101. [PMID: 34654297 DOI: 10.1063/5.0065765] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Immature hepatitis B virus (HBV) captures nucleotides in its capsid for reverse transcription. The nucleotides and nucleotide analog drugs, which are triphosphorylated and negatively charged in the cell, approach the capsid via diffusion and are absorbed into it. In this study, we performed a long-time molecular dynamics calculation of the entire HBV capsid containing pregenome RNA to investigate the interactions between the capsid and negatively charged substances. Electric field analysis demonstrated that negatively charged substances can approach the HBV capsid by thermal motion, avoiding spikes. The substances then migrate all over the floor of the HBV capsid. Finally, they find pores through which they can pass through the HBV capsid shell. Free energy profiles were calculated along these pores for small ions to understand their permeability through the pores. Anions (Cl-) showed higher free energy barriers than cations (Na+ and K+) through all pores, and the permeation rate of Cl- was eight times slower than that of K+ or Na+. Furthermore, the ions were more stable in the capsid than in the bulk water. Thus, the HBV capsid exerts ion selectivity for uptake and provides an environment for ions, such as nucleotides and nucleotide analog drugs, to be stabilized within the capsid.
Collapse
Affiliation(s)
- Kazushi Fujimoto
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Youhei Yamaguchi
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Ryo Urano
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Tetsuya Ishikawa
- Department of Integrated Health Sciences, Nagoya University, Nagoya, Japan
| | | | - Yasuhito Tanaka
- Department of Virology and Liver Unit, Nagoya City University, Nagoya, Japan
| | | | - Susumu Okazaki
- Department of Advanced Materials Science, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
46
|
Jones PE, Pérez-Segura C, Bryer AJ, Perilla JR, Hadden-Perilla JA. Molecular dynamics of the viral life cycle: progress and prospects. Curr Opin Virol 2021; 50:128-138. [PMID: 34464843 PMCID: PMC8651149 DOI: 10.1016/j.coviro.2021.08.003] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Revised: 08/09/2021] [Accepted: 08/09/2021] [Indexed: 01/29/2023]
Abstract
Molecular dynamics (MD) simulations across spatiotemporal resolutions are widely applied to study viruses and represent the central technique uniting the field of computational virology. We discuss the progress of MD in elucidating the dynamics of the viral life cycle, including the status of modeling intact extracellular virions and leveraging advanced simulations to mimic active life cycle processes. We further remark on the prospects of MD for continued contributions to the basic science characterization of viruses, especially given the increasing availability of high-quality experimental data and supercomputing power. Overall, integrative computational methods that are closely guided by experiments are unmatched in the level of detail they provide, enabling-now and in the future-new discoveries relevant to thwarting viral infection.
Collapse
Affiliation(s)
- Peter Eugene Jones
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Carolina Pérez-Segura
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Alexander J Bryer
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Jodi A Hadden-Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States.
| |
Collapse
|
47
|
Schultze S, Grubmüller H. Time-Lagged Independent Component Analysis of Random Walks and Protein Dynamics. J Chem Theory Comput 2021; 17:5766-5776. [PMID: 34449229 PMCID: PMC8444338 DOI: 10.1021/acs.jctc.1c00273] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
![]()
Time-lagged independent
component analysis (tICA) is a widely used
dimension reduction method for the analysis of molecular dynamics
(MD) trajectories and has proven particularly useful for the construction
of protein dynamics Markov models. It identifies those “slow”
collective degrees of freedom onto which the projections of a given
trajectory show maximal autocorrelation for a given lag time. Here
we ask how much information on the actual protein dynamics and, in
particular, the free energy landscape that governs these dynamics
the tICA-projections of MD-trajectories contain, as opposed to noise
due to the inherently stochastic nature of each trajectory. To answer
this question, we have analyzed the tICA-projections of high dimensional
random walks using a combination of analytical and numerical methods.
We find that the projections resemble cosine functions and strongly
depend on the lag time, exhibiting strikingly complex behavior. In
particular, and contrary to previous studies of principal component
projections, the projections change noncontinuously with increasing
lag time. The tICA-projections of selected 1 μs protein trajectories
and those of random walks are strikingly similar, particularly for
larger proteins, suggesting that these trajectories contain only little
information on the energy landscape that governs the actual protein
dynamics. Further the tICA-projections of random walks show clusters
very similar to those observed for the protein trajectories, suggesting
that clusters in the tICA-projections of protein trajectories do not
necessarily reflect local minima in the free energy landscape. We
also conclude that, in addition to the previous finding that certain
ensemble properties of nonconverged protein trajectories resemble
those of random walks; this is also true for their time correlations.
Collapse
Affiliation(s)
- Steffen Schultze
- Max Planck Institute for Biophysical Chemistry, Göttingen, 37077, Germany
| | - Helmut Grubmüller
- Max Planck Institute for Biophysical Chemistry, Göttingen, 37077, Germany
| |
Collapse
|
48
|
Abstract
The disassembly of a viral capsid leading to the release of its genetic material into the host cell is a fundamental step in viral infection. In hepatitis B virus (HBV), the capsid consists of identical protein monomers that dimerize and then arrange themselves into pentamers or hexamers on the capsid surface. By applying atomistic molecular dynamics simulation to an entire solvated HBV capsid subjected to a uniform mechanical stress protocol, we monitor the capsid-disassembly process and analyze the process down to the level of individual amino acids in 20 independent simulation replicas. The strain of an isotropic external force, combined with structural fluctuations, causes structurally heterogeneous cracks to appear in the HBV capsid. Analysis of the monomer-monomer interfaces reveals that, in contrast to the expectation from purely mechanical considerations, the cracks mainly occur within hexameric sites, whereas pentameric sites remain largely intact. Only a small subset of the capsid protein monomers, different in each simulation, are engaged in each instance of disassembly. We identify specific residues whose interactions are most readily lost during disassembly; R127, I139, Y132, N136, A137, and V149 are among the hot spots at the interfaces between dimers that lie within hexamers, leading to disassembly. The majority of these hot-spot residues are conserved by evolution, hinting to their importance for disassembly by avoiding overstabilization of capsids.
Collapse
|
49
|
Wilson E, Vant J, Layton J, Boyd R, Lee H, Turilli M, Hernández B, Wilkinson S, Jha S, Gupta C, Sarkar D, Singharoy A. Large-Scale Molecular Dynamics Simulations of Cellular Compartments. Methods Mol Biol 2021; 2302:335-356. [PMID: 33877636 DOI: 10.1007/978-1-0716-1394-8_18] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.
Collapse
Affiliation(s)
- Eric Wilson
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA
| | - John Vant
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA
| | - Jacob Layton
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA
| | - Ryan Boyd
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA
| | - Hyungro Lee
- RADICAL, ECE, Rutgers University, Piscataway, NJ, USA
| | | | | | | | - Shantenu Jha
- RADICAL, ECE, Rutgers University, Piscataway, NJ, USA.,Brookhaven National Laboratory, Upton, NY, USA
| | - Chitrak Gupta
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA.
| | - Daipayan Sarkar
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA. .,Department of Biological Sciences, Purdue University, West Lafayette, IN, USA.
| | - Abhishek Singharoy
- The School of Molecular Sciences, Arizona State University, Tempe, AZ, USA.
| |
Collapse
|
50
|
Dongre AV, Das S, Bellur A, Kumar S, Chandrashekarmath A, Karmakar T, Balaram P, Balasubramanian S, Balaram H. Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation. Biophys J 2021; 120:3732-3746. [PMID: 34302792 DOI: 10.1016/j.bpj.2021.07.014] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Revised: 06/18/2021] [Accepted: 07/19/2021] [Indexed: 10/20/2022] Open
Abstract
Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.
Collapse
Affiliation(s)
- Aparna Vilas Dongre
- Molecular Biology and Genetics Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India
| | - Sudip Das
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India
| | - Asutosh Bellur
- Molecular Biology and Genetics Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India
| | - Sanjeev Kumar
- Molecular Biology and Genetics Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India; National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India
| | - Anusha Chandrashekarmath
- Molecular Biology and Genetics Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India
| | - Tarak Karmakar
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India; Department of Chemistry and Applied Biosciences, ETH Zurich, Lugano, Ticino, Switzerland; Facoltà di Informatica, Istituto di Scienze Computationali, Università della Svizzera Italiana, Lugano, Ticino, Switzerland
| | - Padmanabhan Balaram
- National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India; Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India
| | - Sundaram Balasubramanian
- Chemistry and Physics of Materials Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India.
| | - Hemalatha Balaram
- Molecular Biology and Genetics Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bangalore, India.
| |
Collapse
|