1
|
Niblo JK, Swartley JR, Zhang Z, DuBay KH. 2D capsid formation within an oscillatory energy landscape: orderly self-assembly depends on the interplay between a dynamic potential and intrinsic relaxation times. SOFT MATTER 2024; 20:6702-6713. [PMID: 39046256 DOI: 10.1039/d4sm00455h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
Multiple dissipative self-assembly protocols designed to create novel structures or to reduce kinetic traps have recently emerged. Specifically, temporal oscillations of particle interactions have been shown effective at both aims, but investigations thus far have focused on systems of simple colloids or their binary mixtures. In this work, we expand our understanding of the effect of temporally oscillating interactions to a two-dimensional coarse-grained viral capsid-like model that undergoes a self-limited assembly. This model includes multiple intrinsic relaxation times due to the internal structure of the capsid subunits and, under certain interaction regimes, proceeds via a two-step nucleation mechanism. We find that oscillations much faster than the local intrinsic relaxation times can be described via a time averaged inter-particle potential across a wide range of interaction strengths, while oscillations much slower than these relaxation times result in structures that adapt to the attraction strength of the current half-cycle. Interestingly, oscillation periods similar to these relaxation times shift the interaction window over which orderly assembly occurs by enabling error correction during the half-cycles with weaker attractions. Our results provide fundamental insights to non-equilibrium self-assembly on temporally variant energy landscapes.
Collapse
Affiliation(s)
- Jessica K Niblo
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| | - Jacob R Swartley
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| | - Zhongmin Zhang
- Department of Chemistry, University of North Carolina at Chapel Hill, Campus Box 3290, Chapel Hill, NC 27599-3290, USA
| | - Kateri H DuBay
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| |
Collapse
|
2
|
Garza CM, Holcomb M, Santos-Martins D, Torbett BE, Forli S. IP6 and PF74 affect HIV-1 Capsid Stability through Modulation of Hexamer-Hexamer Tilt Angle Preference. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.11.584513. [PMID: 38559213 PMCID: PMC10979974 DOI: 10.1101/2024.03.11.584513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
The HIV-1 capsid is an irregularly shaped complex of about 1200 protein chains containing the viral genome and several viral proteins. Together, these components are the key to unlocking passage into the nucleus, allowing for permanent integration of the viral genome into the host cell genome. Recent interest into the role of the capsid in viral replication has been driven by the approval of the first-in-class drug lenacapavir, which marks the first drug approved to target a non-enzymatic HIV-1 viral protein. In addition to lenacapavir, other small molecules such as the drug-like compound PF74, and the anionic sugar inositolhexakisphosphate (IP6), are known to impact capsid stability, and although this is widely accepted as a therapeutic effect, the mechanisms through which they do so remain unknown. In this study, we employed a systematic atomistic simulation approach to study the impact of molecules bound to hexamers at the central pore (IP6) and the FG-binding site (PF74) on capsid oligomer dynamics, compared to apo hexamers and pentamers. We found that neither small molecule had a sizeable impact on the free energy of binding of the interface between neighboring hexamers but that both had impacts on the free energy profiles of performing angular deformations to the pair of oligomers akin to the variations in curvature along the irregular surface of the capsid. The IP6 cofactor, on one hand, stabilizes a pair of neighboring hexamers in their flattest configurations, whereas without IP6, the hexamers prefer a high tilt angle between them. On the other hand, having PF74 bound introduces a strong preference for intermediate tilt angles. These results suggest that structural instability is a natural feature of the HIV-1 capsid which is modulated by molecules bound in either the central pore or the FG-binding site. Such modulators, despite sharing many of the same effects on non-bonded interactions at the various protein-protein interfaces, have decidedly different effects on the flexibility of the complex. This study provides a detailed model of the HIV-1 capsid and its interactions with small molecules, informing structure-based drug design, as well as experimental design and interpretation.
Collapse
|
3
|
Grossman D, Katzav E, Sharon E. Packing of stiff rods on ellipsoids: Geometry. Phys Rev E 2021; 103:013001. [PMID: 33601586 DOI: 10.1103/physreve.103.013001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 12/17/2020] [Indexed: 06/12/2023]
Abstract
We suggest a geometrical mechanism for the ordering of slender filaments inside nonisotropic containers, using cortical microtubules in plant cells and the packing of viral genetic material inside capsids as concrete examples. We show analytically how the shape of the cell affects the ordering of phantom elastic rods that are not self-avoiding (i.e., self-crossing is allowed). We find that for oblate cells, the preferred orientation is along the equator, while for prolate spheroids with an aspect ratio close to 1, the orientation is along the principal (long axis). Surprisingly, at a high enough aspect ratio, a configurational phase transition occurs and the rods no longer point along the principal axis, but at an angle to it, due to high curvature at the poles. We discuss some of the possible effects of self-avoidance using energy considerations. These results are relevant to other packing problems as well, such as the spooling of filament in the industry or spider silk inside water droplets.
Collapse
Affiliation(s)
- Doron Grossman
- Racah Institute of Physics, Hebrew University, Jerusalem 9190401, Israel
| | - Eytan Katzav
- Racah Institute of Physics, Hebrew University, Jerusalem 9190401, Israel
| | - Eran Sharon
- Racah Institute of Physics, Hebrew University, Jerusalem 9190401, Israel
| |
Collapse
|
4
|
Jefferys EE, Sansom MSP. Computational Virology: Molecular Simulations of Virus Dynamics and Interactions. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2019; 1215:201-233. [DOI: 10.1007/978-3-030-14741-9_10] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
5
|
Abstract
An implicit aim in cellular infection biology is to understand the mechanisms how viruses, microbes, eukaryotic parasites, and fungi usurp the functions of host cells and cause disease. Mechanistic insight is a deep understanding of the biophysical and biochemical processes that give rise to an observable phenomenon. It is typically subject to falsification, that is, it is accessible to experimentation and empirical data acquisition. This is different from logic and mathematics, which are not empirical, but built on systems of inherently consistent axioms. Here, we argue that modeling and computer simulation, combined with mechanistic insights, yields unprecedented deep understanding of phenomena in biology and especially in virus infections by providing a way of showing sufficiency of a hypothetical mechanism. This ideally complements the necessity statements accessible to empirical falsification by additional positive evidence. We discuss how computational implementations of mathematical models can assist and enhance the quantitative measurements of infection dynamics of enveloped and non-enveloped viruses and thereby help generating causal insights into virus infection biology.
Collapse
|
6
|
Okumura H, Higashi M, Yoshida Y, Sato H, Akiyama R. Theoretical approaches for dynamical ordering of biomolecular systems. Biochim Biophys Acta Gen Subj 2017; 1862:212-228. [PMID: 28988931 DOI: 10.1016/j.bbagen.2017.10.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 09/30/2017] [Accepted: 10/04/2017] [Indexed: 01/21/2023]
Abstract
BACKGROUND Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies. SCOPE OF REVIEW We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified. MAJOR CONCLUSIONS The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function. GENERAL SIGNIFICANCE These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.
Collapse
Affiliation(s)
- Hisashi Okumura
- Research Center for Computational Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan; Department of Structural Molecular Science, The Graduate University for Advanced Studies, Okazaki, Aichi 444-8585, Japan.
| | - Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
| | - Yuichiro Yoshida
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan; Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| | - Ryo Akiyama
- Department of Chemistry, Kyushu University, Fukuoka 819-0395, Japan
| |
Collapse
|
7
|
Kobayashi C, Jung J, Matsunaga Y, Mori T, Ando T, Tamura K, Kamiya M, Sugita Y. GENESIS 1.1: A hybrid-parallel molecular dynamics simulator with enhanced sampling algorithms on multiple computational platforms. J Comput Chem 2017; 38:2193-2206. [PMID: 28718930 DOI: 10.1002/jcc.24874] [Citation(s) in RCA: 111] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Revised: 06/08/2017] [Accepted: 06/09/2017] [Indexed: 01/09/2023]
Abstract
GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Yasuhiro Matsunaga
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,JST PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Tadashi Ando
- Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center Computational Biology Research Core, 1-6-5 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Department of Applied Electronics, Faculty of Industrial Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan.,Water Frontier Science and Technology Research Center, Research Institute for Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan.,Research Division of Multiscale Interfacial Thermofluid Dynamics, Research Institute for Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan
| | - Koichi Tamura
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Motoshi Kamiya
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan.,Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center Computational Biology Research Core, 1-6-5 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| |
Collapse
|
8
|
Sha H, Zhu F. Parameter Optimization for Interaction between C-Terminal Domains of HIV-1 Capsid Protein. J Chem Inf Model 2017; 57:1134-1141. [PMID: 28426204 DOI: 10.1021/acs.jcim.7b00011] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
HIV-1 capsid proteins (CAs) assemble into a capsid that encloses the viral RNA. The binding between a pair of C-terminal domains (CTDs) constitutes a major interface in both the CA dimers and the large CA assemblies. Here, we attempt to use a general residue-level coarse-grained model to describe the interaction between two isolated CTDs in Monte Carlo simulations. With the standard parameters that depend only on the residue types, the model predicts a much weaker binding in comparison to the experiments. Detailed analysis reveals that some Lennard-Jones parameters are not compatible with the experimental CTD dimer structure, thus resulting in an unfavorable interaction energy. To improve the model for the CTD binding, we introduce ad hoc modifications to a small number of Lennard-Jones parameters for some specific pairs of residues at the binding interface. Through a series of extensive Monte Carlo simulations, we identify the optimal parameters for the CTD-CTD interactions. With the refined model parameters, both the binding affinity (with a dissociation constant of 13 ± 2 μM) and the binding mode are in good agreement with the experimental data. This study demonstrates that the general interaction model based on the Lennard-Jones potential, with some modest adjustment of the parameters for key residues, could correctly reproduce the reversible protein binding, thus potentially applicable for simulating the thermodynamics of the CA assemblies.
Collapse
Affiliation(s)
- Hao Sha
- Department of Physics, Indiana University - Purdue University Indianapolis , 402 N. Blackford, LD 154, Indianapolis, Indiana 46202, United States
| | - Fangqiang Zhu
- Department of Physics, Indiana University - Purdue University Indianapolis , 402 N. Blackford, LD 154, Indianapolis, Indiana 46202, United States
| |
Collapse
|
9
|
Sadiq SK. Reaction-diffusion basis of retroviral infectivity. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2016; 374:rsta.2016.0148. [PMID: 27698042 PMCID: PMC5052732 DOI: 10.1098/rsta.2016.0148] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 08/22/2016] [Indexed: 05/27/2023]
Abstract
Retrovirus particle (virion) infectivity requires diffusion and clustering of multiple transmembrane envelope proteins (Env3) on the virion exterior, yet is triggered by protease-dependent degradation of a partially occluding, membrane-bound Gag polyprotein lattice on the virion interior. The physical mechanism underlying such coupling is unclear and only indirectly accessible via experiment. Modelling stands to provide insight but the required spatio-temporal range far exceeds current accessibility by all-atom or even coarse-grained molecular dynamics simulations. Nor do such approaches account for chemical reactions, while conversely, reaction kinetics approaches handle neither diffusion nor clustering. Here, a recently developed multiscale approach is considered that applies an ultra-coarse-graining scheme to treat entire proteins at near-single particle resolution, but which also couples chemical reactions with diffusion and interactions. A model is developed of Env3 molecules embedded in a truncated Gag lattice composed of membrane-bound matrix proteins linked to capsid subunits, with freely diffusing protease molecules. Simulations suggest that in the presence of Gag but in the absence of lateral lattice-forming interactions, Env3 diffuses comparably to Gag-absent Env3 Initial immobility of Env3 is conferred through lateral caging by matrix trimers vertically coupled to the underlying hexameric capsid layer. Gag cleavage by protease vertically decouples the matrix and capsid layers, induces both matrix and Env3 diffusion, and permits Env3 clustering. Spreading across the entire membrane surface reduces crowding, in turn, enhancing the effect and promoting infectivity.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- S Kashif Sadiq
- Infection Biology Unit, Universitat Pompeu Fabra, Barcelona Biomedical Research Park (PRBB), C/Doctor Aiguader 88, 08003 Barcelona, Spain Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| |
Collapse
|
10
|
Smith GR, Xie L, Schwartz R. Modeling Effects of RNA on Capsid Assembly Pathways via Coarse-Grained Stochastic Simulation. PLoS One 2016; 11:e0156547. [PMID: 27244559 PMCID: PMC4887116 DOI: 10.1371/journal.pone.0156547] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2015] [Accepted: 05/16/2016] [Indexed: 12/02/2022] Open
Abstract
The environment of a living cell is vastly different from that of an in vitro reaction system, an issue that presents great challenges to the use of in vitro models, or computer simulations based on them, for understanding biochemistry in vivo. Virus capsids make an excellent model system for such questions because they typically have few distinct components, making them amenable to in vitro and modeling studies, yet their assembly can involve complex networks of possible reactions that cannot be resolved in detail by any current experimental technology. We previously fit kinetic simulation parameters to bulk in vitro assembly data to yield a close match between simulated and real data, and then used the simulations to study features of assembly that cannot be monitored experimentally. The present work seeks to project how assembly in these simulations fit to in vitro data would be altered by computationally adding features of the cellular environment to the system, specifically the presence of nucleic acid about which many capsids assemble. The major challenge of such work is computational: simulating fine-scale assembly pathways on the scale and in the parameter domains of real viruses is far too computationally costly to allow for explicit models of nucleic acid interaction. We bypass that limitation by applying analytical models of nucleic acid effects to adjust kinetic rate parameters learned from in vitro data to see how these adjustments, singly or in combination, might affect fine-scale assembly progress. The resulting simulations exhibit surprising behavioral complexity, with distinct effects often acting synergistically to drive efficient assembly and alter pathways relative to the in vitro model. The work demonstrates how computer simulations can help us understand how assembly might differ between the in vitro and in vivo environments and what features of the cellular environment account for these differences.
Collapse
Affiliation(s)
- Gregory R. Smith
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America
| | - Lu Xie
- Joint Carnegie Mellon/University of Pittsburgh Ph.D. Program in Computational Biology, Pittsburgh, Pennsylvania, United States of America
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America
| | - Russell Schwartz
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America
- * E-mail:
| |
Collapse
|
11
|
Abstract
The HIV genome materials are encaged by a proteinaceous shell called the capsid, constructed from ∼1000-1500 copies of the capsid proteins. Because its stability and integrity are critical to the normal life cycle and infectivity of the virus, the HIV capsid is a promising antiviral drug target. In this paper, we review the studies shaping our understanding of the structure and dynamics of the capsid proteins and various forms of their assemblies, as well as the assembly mechanism.
Collapse
Affiliation(s)
- Bo Chen
- Department of Physics, University of Central Florida , Orlando, Florida 32816, United States
| |
Collapse
|
12
|
Qiao X, Jeon J, Weber J, Zhu F, Chen B. Mechanism of polymorphism and curvature of HIV capsid assemblies probed by 3D simulations with a novel coarse grain model. Biochim Biophys Acta Gen Subj 2015; 1850:2353-67. [PMID: 26318016 DOI: 10.1016/j.bbagen.2015.08.017] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2015] [Revised: 08/17/2015] [Accepted: 08/19/2015] [Indexed: 12/11/2022]
Abstract
BACKGROUND During the maturation process, HIV capsid proteins self-assemble into polymorphic capsids. The strong polymorphism precludes high resolution structural characterization under in vivo conditions. In spite of the determination of structural models for various in vitro assemblies of HIV capsid proteins, the assembly mechanism is still not well-understood. METHODS We report 3D simulations of HIV capsid proteins by a novel coarse grain model that captures the backbone of the rigid segments in the protein accurately. The effects of protein dynamics on assembly are emulated by a static ensemble of subunits in conformations derived from molecular dynamics simulation. RESULTS We show that HIV capsid proteins robustly assemble into hexameric lattices in a range of conditions where trimers of dimeric subunits are the dominant oligomeric intermediates. Variations of hexameric lattice curvatures are observed in simulations with subunits of variable inter-domain orientations mimicking the conformation distribution in solution. Simulations with subunits based on pentameric structural models lead to assembly of sharp curved structures resembling the tips of authentic HIV capsids, along a distinct pathway populated by tetramers and pentamers with the characteristic quasi-equivalency of viral capsids. CONCLUSIONS Our results suggest that the polymorphism assembly is triggered by the inter-domain dynamics of HIV capsid proteins in solution. The assembly of highly curved structures arises from proteins in conformation with a highly specific inter-domain orientation. SIGNIFICANCE Our work proposes a mechanism of HIV capsid assembly based on available structural data, which can be readily verified. Our model can be applied to other large biomolecular assemblies.
Collapse
Affiliation(s)
- Xin Qiao
- Department of Physics, University of Central Florida, 4000 Central Florida Blvd, Orlando, FL 32816, USA
| | - Jaekyun Jeon
- Department of Physics, University of Central Florida, 4000 Central Florida Blvd, Orlando, FL 32816, USA
| | - Jeff Weber
- Department of Physics, University of Central Florida, 4000 Central Florida Blvd, Orlando, FL 32816, USA
| | - Fangqiang Zhu
- Department of Physics, Indiana University - Purdue University Indianapolis, IN, USA
| | - Bo Chen
- Department of Physics, University of Central Florida, 4000 Central Florida Blvd, Orlando, FL 32816, USA.
| |
Collapse
|
13
|
Affiliation(s)
- Fangqiang Zhu
- Department
of Physics, Indiana University - Purdue University, Indianapolis, Indiana 46202, United States
| | - Bo Chen
- Department
of Physics, University of Central Florida, Orlando, Florida 32816, United States
| |
Collapse
|
14
|
Zhang Z. Systematic methods for defining coarse-grained maps in large biomolecules. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2015; 827:33-48. [PMID: 25387958 DOI: 10.1007/978-94-017-9245-5_4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/21/2023]
Abstract
Large biomolecules are involved in many important biological processes. It would be difficult to use large-scale atomistic molecular dynamics (MD) simulations to study the functional motions of these systems because of the computational expense. Therefore various coarse-grained (CG) approaches have attracted rapidly growing interest, which enable simulations of large biomolecules over longer effective timescales than all-atom MD simulations. The first issue in CG modeling is to construct CG maps from atomic structures. In this chapter, we review the recent development of a novel and systematic method for constructing CG representations of arbitrarily complex biomolecules, in order to preserve large-scale and functionally relevant essential dynamics (ED) at the CG level. In this ED-CG scheme, the essential dynamics can be characterized by principal component analysis (PCA) on a structural ensemble, or elastic network model (ENM) of a single atomic structure. Validation and applications of the method cover various biological systems, such as multi-domain proteins, protein complexes, and even biomolecular machines. The results demonstrate that the ED-CG method may serve as a very useful tool for identifying functional dynamics of large biomolecules at the CG level.
Collapse
Affiliation(s)
- Zhiyong Zhang
- Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China,
| |
Collapse
|
15
|
Zhang L, Lua LHL, Middelberg APJ, Sun Y, Connors NK. Biomolecular engineering of virus-like particles aided by computational chemistry methods. Chem Soc Rev 2015; 44:8608-18. [DOI: 10.1039/c5cs00526d] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.
Collapse
Affiliation(s)
- Lin Zhang
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education
- School of Chemical Engineering and Technology
- Tianjin University
- Tianjin 300072, People's Republic of China
| | - Linda H. L. Lua
- Protein Expression Facility
- The University of Queensland
- Brisbane, Australia
| | - Anton P. J. Middelberg
- Australian Institute for Bioengineering and Nanotechnology
- The University of Queensland
- Brisbane, Australia
| | - Yan Sun
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education
- School of Chemical Engineering and Technology
- Tianjin University
- Tianjin 300072, People's Republic of China
| | - Natalie K. Connors
- Australian Institute for Bioengineering and Nanotechnology
- The University of Queensland
- Brisbane, Australia
| |
Collapse
|
16
|
Smith GR, Xie L, Lee B, Schwartz R. Applying molecular crowding models to simulations of virus capsid assembly in vitro. Biophys J 2014; 106:310-20. [PMID: 24411263 DOI: 10.1016/j.bpj.2013.11.022] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Revised: 11/03/2013] [Accepted: 11/11/2013] [Indexed: 11/29/2022] Open
Abstract
Virus capsid assembly has been widely studied as a biophysical system, both for its biological and medical significance and as an important model for complex self-assembly processes. No current technology can monitor assembly in detail and what information we have on assembly kinetics comes exclusively from in vitro studies. There are many differences between the intracellular environment and that of an in vitro assembly assay, however, that might be expected to alter assembly pathways. Here, we explore one specific feature characteristic of the intracellular environment and known to have large effects on macromolecular assembly processes: molecular crowding. We combine prior particle simulation methods for estimating crowding effects with coarse-grained stochastic models of capsid assembly, using the crowding models to adjust kinetics of capsid simulations to examine possible effects of crowding on assembly pathways. Simulations suggest a striking difference depending on whether or not a system uses nucleation-limited assembly, with crowding tending to promote off-pathway growth in a nonnucleation-limited model but often enhancing assembly efficiency at high crowding levels even while impeding it at lower crowding levels in a nucleation-limited model. These models may help us understand how complicated assembly systems may have evolved to function with high efficiency and fidelity in the densely crowded environment of the cell.
Collapse
Affiliation(s)
- Gregory R Smith
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Lu Xie
- Joint Carnegie Mellon/University of Pittsburgh Ph.D. Program in Computational Biology, Carnegie Mellon University, Pittsburgh, Pennsylvania; Lane Center for Computational Biology, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Byoungkoo Lee
- Department of Mathematics and Statistics, Georgia State University, Atlanta, Georgia
| | - Russell Schwartz
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania; Lane Center for Computational Biology, Carnegie Mellon University, Pittsburgh, Pennsylvania.
| |
Collapse
|
17
|
Kar P, Feig M. Recent advances in transferable coarse-grained modeling of proteins. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 96:143-80. [PMID: 25443957 PMCID: PMC5366245 DOI: 10.1016/bs.apcsb.2014.06.005] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Computer simulations are indispensable tools for studying the structure and dynamics of biological macromolecules. Biochemical processes occur on different scales of length and time. Atomistic simulations cannot cover the relevant spatiotemporal scales at which the cellular processes occur. To address this challenge, coarse-grained (CG) modeling of the biological systems is employed. Over the last few years, many CG models for proteins continue to be developed. However, many of them are not transferable with respect to different systems and different environments. In this review, we discuss those CG protein models that are transferable and that retain chemical specificity. We restrict ourselves to CG models of soluble proteins only. We also briefly review recent progress made in the multiscale hybrid all-atom/CG simulations of proteins.
Collapse
Affiliation(s)
- Parimal Kar
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA; Department of Chemistry, Michigan State University, East Lansing, Michigan, USA.
| |
Collapse
|
18
|
Abstract
By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
19
|
Takahashi K, Oda T, Naruse K. Coarse-grained molecular dynamics simulations of biomolecules. AIMS BIOPHYSICS 2014. [DOI: 10.3934/biophy.2014.1.1] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
20
|
Abstract
All matter has to obey the general laws of physics and living matter is not an exception. Viruses have not only learnt how to cope with them, but have managed to use them for their own survival. In this chapter we will review some of the exciting physics behind viruses and discuss simple physical models that can shed some light on different aspects of the viral life cycle and viral properties. In particular, we will focus on how the structure and shape of the capsid, its assembly and stability, and the entry and exit of viral particles and their genomes can be understood using fundamental physics theories.
Collapse
Affiliation(s)
- Antoni Luque
- Department of Fundamental Physics, Universitat de Barcelona, c/Martí i Franquès 1, 08028, Barcelona, Spain
| | | |
Collapse
|
21
|
Grime JMA, Voth GA. Early stages of the HIV-1 capsid protein lattice formation. Biophys J 2012; 103:1774-83. [PMID: 23083721 DOI: 10.1016/j.bpj.2012.09.007] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2012] [Revised: 08/20/2012] [Accepted: 09/04/2012] [Indexed: 11/26/2022] Open
Abstract
The early stages in the formation of the HIV-1 capsid (CA) protein lattice are investigated. The underlying coarse-grained (CG) model is parameterized directly from experimental data and examined under various native contact interaction strengths, CA dimer interfacial configurations, and local surface curvatures. The mechanism of early contiguous mature-style CA p6 lattice formation is explored, and a trimer-of-dimers structure is found to be crucial for CA lattice production. Quasi-equivalent generation of both the pentamer and hexamer components of the HIV-1 viral CA is also demonstrated, and the formation of pentamers is shown to be highly sensitive to local curvature, supporting the view that such inclusions in high-curvature regions allow closure of the viral CA surface. The complicated behavior of CA lattice self-assembly is shown to be reducible to a relatively simple function of the trimer-of-dimers behavior.
Collapse
Affiliation(s)
- John M A Grime
- Department of Chemistry, University of Chicago, Chicago, Illinois, USA
| | | |
Collapse
|
22
|
Ode H, Nakashima M, Kitamura S, Sugiura W, Sato H. Molecular dynamics simulation in virus research. Front Microbiol 2012; 3:258. [PMID: 22833741 PMCID: PMC3400276 DOI: 10.3389/fmicb.2012.00258] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2012] [Accepted: 07/02/2012] [Indexed: 01/24/2023] Open
Abstract
Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus–host interactions and evolution, and present future perspectives on this technique.
Collapse
Affiliation(s)
- Hirotaka Ode
- Clinical Research Center, National Hospital Organization Nagoya Medical Center Nagoya, Aichi, Japan
| | | | | | | | | |
Collapse
|
23
|
Chaudhri A, Zarraga IE, Kamerzell TJ, Brandt JP, Patapoff TW, Shire SJ, Voth GA. Coarse-grained modeling of the self-association of therapeutic monoclonal antibodies. J Phys Chem B 2012; 116:8045-57. [PMID: 22694284 DOI: 10.1021/jp301140u] [Citation(s) in RCA: 84] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Coarse-grained computational models of two therapeutic monoclonal antibodies are constructed to understand the effect of domain-level charge-charge electrostatics on the self-association phenomena at high protein concentrations. The coarse-grained representations of the individual antibodies are constructed using an elastic network normal-mode analysis. Two different models are constructed for each antibody for a compact Y-shaped and an extended Y-shaped configuration. The resulting simulations of these coarse-grained antibodies that interact through screened electrostatics are done at six different concentrations. It is observed that a particular monoclonal antibody (hereafter referred to as MAb1) forms three-dimensional heterogeneous structures with dense regions or clusters compared to a different monoclonal antibody (hereafter referred to as MAb2) that forms more homogeneous structures (no clusters). These structures, together with the potential mean force (PMF) and radial distribution functions (RDF) between pairs of coarse-grained regions on the MAbs, are qualitatively consistent with the experimental observation that MAb1 has a significantly higher viscosity compared to MAb2, especially at concentrations >50 mg/mL, even though the only difference between the MAbs lies with a few amino acids at the antigen-binding loops (CDRs). It is also observed that the structures in MAb1 are formed due to stronger Fab-Fab interactions in corroboration with experimental observations. Evidence is also shown that Fab-Fc interactions can be equally important in addition to Fab-Fab interactions. The coarse-grained representations are effective in picking up differences based on local charge distributions of domains and make predictions on the self-association characteristics of these protein solutions. This is the first computational study of its kind to show that there are differences in structures formed by two different monoclonal antibodies at high concentrations.
Collapse
Affiliation(s)
- Anuj Chaudhri
- Department of Chemistry, James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| | | | | | | | | | | | | |
Collapse
|
24
|
Tsiang M, Niedziela-Majka A, Hung M, Jin D, Hu E, Yant S, Samuel D, Liu X, Sakowicz R. A trimer of dimers is the basic building block for human immunodeficiency virus-1 capsid assembly. Biochemistry 2012; 51:4416-28. [PMID: 22564075 DOI: 10.1021/bi300052h] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Human immunodeficiency virus-1 (HIV-1) capsid protein (CA) has become a target of antiviral drug design in recent years. The recognition that binding of small molecules to the CA protein can result in the perturbation of capsid assembly or disassembly has led to mathematical modeling of the process. Although a number of capsid assembly models have been developed using biophysical parameters of the CA protein obtained experimentally, there is currently no model of CA polymerization that can be practically used to analyze in vitro CA polymerization data to facilitate drug discovery. Herein, we describe an equilibrium model of CA polymerization for the kinetic analysis of in vitro assembly of CA into polymer tubes. This new mathematical model has been used to assess whether a triangular trimer of dimers rather than a hexagonal hexamer can be the basic capsomere building block of CA polymer. The model allowed us to quantify for the first time the affinity for each of the four crucial interfaces involved in the polymerization process and indicated that the trimerization of CA dimers is a relatively slow step in CA polymerization in vitro. For wild-type CA, these four interfaces include the interface between two monomers of a CA dimer (K(D) = 6.6 μM), the interface between any two dimers within a CA trimer of dimers (K(D) = 32 nM), and two types of interfaces between neighboring trimers of dimers, either within the same ring around the perimeter of the polymer tube (K(D) = 438 nM) or from two adjacent rings (K(D) = 147 nM). A comparative analysis of the interface dissociation constants between wild-type and two mutant CA proteins, cross-linked hexamer (A14C/E45C/W184A/M185A) and A14C/E45C, yielded results that are consistent with the trimer of dimers with a triangular geometry being the capsomere building block involved in CA polymer growth. This work provides additional insights into the mechanism of HIV-1 CA assembly and may prove useful in elucidating how small molecule CA binding agents may disturb this essential step in the HIV-1 life cycle.
Collapse
Affiliation(s)
- Manuel Tsiang
- Gilead Sciences, Foster City, California 94404, United States.
| | | | | | | | | | | | | | | | | |
Collapse
|
25
|
Larini L, Shea JE. Coarse-Grained Modeling of Simple Molecules at Different Resolutions in the Absence of Good Sampling. J Phys Chem B 2012; 116:8337-49. [DOI: 10.1021/jp2097263] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Luca Larini
- Department of Chemistry and Biochemistry and of Physics, University of California, Santa Barbara, Santa Barbara,
California
93106, United States
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry and of Physics, University of California, Santa Barbara, Santa Barbara,
California
93106, United States
| |
Collapse
|
26
|
Pereverzev YV, Prezhdo E, Sokurenko EV. The two-pathway model of the biological catch-bond as a limit of the allosteric model. Biophys J 2012; 101:2026-36. [PMID: 22004757 DOI: 10.1016/j.bpj.2011.09.005] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2011] [Accepted: 08/11/2011] [Indexed: 01/25/2023] Open
Abstract
Catch-binding is a counterintuitive phenomenon in which the lifetime of a receptor/ligand bond increases when a force is applied to break the bond. Several mechanisms have been proposed to rationalize catch-binding. In the two-pathway model, the force drives the system away from its native dissociation pathway into an alternative pathway involving a higher energy barrier. Here, we analyze an allosteric model suggesting that a force applied to the complex alters the distribution of receptor conformations, and as a result, induces changes in the ligand-binding site. The model assumes explicitly that the allosteric transitions govern the properties of the ligand site. We demonstrate that the dynamics of the ligand is described by two relaxation times, one of which arises from the allosteric site. Therefore, we argue that one can characterize the allosteric transitions by studying the receptor/ligand binding. We show that the allosteric description reduces to the two-pathway model in the limit when the allosteric transitions are faster than the bond dissociation. The formal results are illustrated with two systems, P-selectin/PSGL-1 and FimH/mannose, subjected to both constant and time-dependent forces. The report advances our understanding of catch-binding by combining alternative physical models into a unified description and makes the problem more tractable for the bond mechanics community.
Collapse
Affiliation(s)
- Yuriy V Pereverzev
- Department of Chemistry, University of Rochester, Rochester, New York, USA
| | | | | |
Collapse
|
27
|
Bezkorovaynaya O, Lukyanov A, Kremer K, Peter C. Multiscale simulation of small peptides: Consistent conformational sampling in atomistic and coarse-grained models. J Comput Chem 2012; 33:937-49. [DOI: 10.1002/jcc.22915] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2011] [Accepted: 12/02/2011] [Indexed: 11/07/2022]
|
28
|
Saunders MG, Voth GA. Coarse-graining of multiprotein assemblies. Curr Opin Struct Biol 2012; 22:144-50. [PMID: 22277168 DOI: 10.1016/j.sbi.2012.01.003] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2011] [Revised: 01/02/2012] [Accepted: 01/04/2012] [Indexed: 11/24/2022]
Abstract
Multiscale models are important tools to elucidate how small changes in local subunit conformations may propagate to affect the properties of macromolecular complexes. We review recent advances in coarse-graining methods for poly-protein assemblies, systems that are composed of many copies of relatively few components, with a particular focus on viral capsids and cytoskeletal filaments. These methods are grouped into two broad categories-mapping methods, which use information from one scale of representation to parameterize a lower resolution model, and bridging methods, which repeatedly connect different scales during simulation-and we provide examples of both classes at different levels of complexity. Collectively, these models illustrate the numerous approaches to information transfer between scales and demonstrate that the complexity required of the model depends in general on the nature of the information sought.
Collapse
Affiliation(s)
- Marissa G Saunders
- Department of Chemistry, Institute for Biophysical Dynamics, James Franck Institute, and Computation Institute, University of Chicago, Chicago, IL 60637, United States
| | | |
Collapse
|
29
|
Chen B, Tycko R. Simulated self-assembly of the HIV-1 capsid: protein shape and native contacts are sufficient for two-dimensional lattice formation. Biophys J 2011; 100:3035-44. [PMID: 21689538 DOI: 10.1016/j.bpj.2011.05.025] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2010] [Revised: 05/09/2011] [Accepted: 05/12/2011] [Indexed: 11/16/2022] Open
Abstract
We report Monte Carlo simulations of the initial stages of self-assembly of the HIV-1 capsid protein (CA), using a coarse-grained representation that mimics the CA backbone structure and intermolecular contacts observed experimentally. A simple representation of N-terminal domain/N-terminal domain and N-terminal domain/C-terminal domain interactions, coupled with the correct protein shape, is sufficient to drive formation of an ordered lattice with the correct hexagonal symmetry in two dimensions. We derive an approximate concentration/temperature phase diagram for lattice formation, and we investigate the pathway by which the lattice develops from initially separated CA dimers. Within this model, lattice formation occurs in two stages: 1), condensation of CA dimers into disordered clusters; and 2), nucleation of the lattice by the appearance of one hexamer unit within a cluster. Trimers of CA dimers are important early intermediates, and pentamers are metastable within clusters. Introduction of a preformed hexamer at the beginning of a Monte Carlo run does not directly seed lattice formation, but does facilitate the formation of large clusters. We discuss possible connections between these simulations and experimental observations concerning CA assembly within HIV-1 and in vitro.
Collapse
Affiliation(s)
- Bo Chen
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, USA
| | | |
Collapse
|
30
|
Krishna V, Larini L. A generalized mean field theory of coarse-graining. J Chem Phys 2011; 135:124103. [DOI: 10.1063/1.3638044] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
|
31
|
Thorpe IF, Goldenberg DP, Voth GA. Exploration of Transferability in Multiscale Coarse-Grained Peptide Models. J Phys Chem B 2011; 115:11911-26. [DOI: 10.1021/jp204455g] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
| | | | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, University of Chicago, Chicago, Illinois, 60637, United States
| |
Collapse
|
32
|
López CS, Eccles JD, Still A, Sloan RE, Barklis RL, Tsagli SM, Barklis E. Determinants of the HIV-1 core assembly pathway. Virology 2011; 417:137-46. [PMID: 21676426 DOI: 10.1016/j.virol.2011.05.011] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2011] [Revised: 05/20/2011] [Accepted: 05/21/2011] [Indexed: 11/28/2022]
Abstract
Based on structural information, we have analyzed the mechanism of mature HIV-1 core assembly and the contributions of structural elements to the assembly process. Through the use of several in vitro assembly assay systems, we have examined details of how capsid (CA) protein helix 1, ß-hairpin and cyclophilin loop elements impact assembly-dependent protein interactions, and we present evidence for a contribution of CA helix 6 to the mature assembly-competent conformation of CA. Additional experiments with mixtures of proteins in assembly reactions provide novel analyses of the mature core assembly mechanism. Our results support a model in which initial assembly products serve as scaffolds for further assembly by converting incoming subunits to assembly proficient conformations, while mutant subunits increase the probability of assembly termination events.
Collapse
Affiliation(s)
- Claudia S López
- Department of Molecular Microbiology and Immunology and Vollum Institute, School of Medicine, Oregon Health & Science University, Portland, OR 97239, USA.
| | | | | | | | | | | | | |
Collapse
|
33
|
Ayton GS, Voth GA. Multiscale computer simulation of the immature HIV-1 virion. Biophys J 2011; 99:2757-65. [PMID: 21044572 DOI: 10.1016/j.bpj.2010.08.018] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2010] [Revised: 07/08/2010] [Accepted: 08/10/2010] [Indexed: 02/06/2023] Open
Abstract
Multiscale computer simulations, employing a combination of experimental data and coarse-graining methods, are used to explore the structure of the immature HIV-1 virion. A coarse-grained (CG) representation is developed for the virion membrane shell and Gag polypeptides using molecular level information. Building on the results from electron cryotomography experiments, the simulations under certain conditions reveal the existence of an incomplete p6 hexameric lattice formed from hexameric bundles of the Gag CA domains. In particular, the formation and stability of the immature Gag lattice at the CG level requires enhanced interfacial interactions of the CA protein C-terminal domains (CTDs). An exact mapping of the CG representation back to the molecular level then allows for detailed atomistic molecular dynamics studies to confirm the existence of these enhanced CA(CTD) interactions and to probe their possible origin. The multiscale simulations further provide insight into potential CA(CTD) mutations that may disrupt or modify the Gag immature lattice assembly process in the immature HIV-1 virion.
Collapse
Affiliation(s)
- Gary S Ayton
- Center for Biophysical Modeling and Simulation and Department of Biochemistry, University of Utah, Salt Lake City, Utah, USA
| | | |
Collapse
|
34
|
Larini L, Lu L, Voth GA. The multiscale coarse-graining method. VI. Implementation of three-body coarse-grained potentials. J Chem Phys 2010; 132:164107. [PMID: 20441258 DOI: 10.1063/1.3394863] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Many methodologies have been proposed to build reliable and computationally fast coarse-grained potentials. Typically, these force fields rely on the assumption that the relevant properties of the system under examination can be reproduced using a pairwise decomposition of the effective coarse-grained forces. In this work it is shown that an extension of the multiscale coarse-graining technique can be employed to parameterize a certain class of two-body and three-body force fields from atomistic configurations. The use of explicit three-body potentials greatly improves the results over the more commonly used two-body approximation. The method proposed here is applied to develop accurate one-site coarse-grained water models.
Collapse
Affiliation(s)
- Luca Larini
- Department of Chemistry and Center for Biophysical Modeling and Simulation, University of Utah, 315 S 1400 E, Salt Lake City, Utah 84112, USA
| | | | | |
Collapse
|