Mechanism of polymorphism and curvature of HIV capsid assemblies probed by 3D simulations with a novel coarse grain model

IP6, PF74 affect HIV-1 capsid stability through modulation of hexamer-hexamer tilt angle preference

The HIV-1 capsid is an irregularly shaped protein complex containing the viral genome and several proteins needed for integration into the host cell genome. Small molecules, such as the drug-like compound PF-3450074 (PF74) and the anionic sugar inositolhexakisphosphate (IP6), are known to impact capsid stability, although the mechanisms through which they do so remain unknown. In this study, we employed atomistic molecular dynamics simulations to study the impact of molecules bound to hexamers at the central pore (IP6) and the FG-binding site (PF74) on the interface between capsid oligomers. We found that the IP6 cofactor stabilizes a pair of neighboring hexamers in their flattest configurations, whereas PF74 introduces a strong preference for intermediate tilt angles. These results suggest that the tilt angle between neighboring hexamers is a primary mechanism for the modulation of capsid stability. In addition, hexamer-pentamer interfaces were highly stable, suggesting that pentamers are likely not the locus of disassembly.

Effect of Pregenomic RNA on the Mechanical Stability of HBV Capsid by Coarse-Grained Molecular Simulations

Hepatitis B virus (HBV) is a double-stranded DNA virus, but its life cycle involves an intermediate stage, during which pregenomic RNA (pgRNA) is encapsulated in the capsid and then reverse-transcribed into the minus DNA strand. These immature HBV virions are the key target for antiviral drug discovery. In this study, we investigate the flexibility and mechanical stability of the HBV capsid containing pgRNA by employing residue-resolved coarse-grained molecular dynamics simulations. The results showed that the presence of pgRNA tends to decrease the overall flexibility of the capsid. In addition, the symmetrically arranged subunits of the capsid show asymmetry in the dominant modes of the conformational fluctuations with or without the presence of pgRNA. Furthermore, the simulations revealed that the presence of pgRNA enhances the overall mechanical stability of the virion particle. Electrostatic interactions between the disordered CTD of capsid and pgRNA were found to play a crucial role in modulating viral mechanical stability. Decreasing the electrostatic interactions by CTD phosphorylation or high salt concentration significantly reduces the mechanical stability of the HBV capsid containing pgRNA. Finally, the 2-fold symmetric sites have been proposed to be the most vulnerable to rupture during the initial stages of capsid disassembly. These findings could enhance our understanding of the physical basis of viral invasion and provide valuable insights into the development of antiviral drugs.

ASAP: An automatic sequential assignment program for congested multidimensional solid state NMR spectra

Accurate signal assignments can be challenging for congested solid-state NMR (ssNMR) spectra. We describe an automatic sequential assignment program (ASAP) to partially overcome this challenge. ASAP takes three input files: the residue type assignments (RTAs) determined from the better-resolved NCACX spectrum, the full peak list of the NCOCX spectrum, and the protein sequence. It integrates our auto-residue type assignment strategy (ARTIST) with the Monte Carlo simulated annealing (MCSA) algorithm to overcome the hurdle for accurate signal assignments caused by incomplete side-chain resonances and spectral congestion. Combined, ASAP demonstrates robust performance and accelerates signal assignments of large proteins (>200 residues) that lack crystalline order.

IP6 and PF74 affect HIV-1 Capsid Stability through Modulation of Hexamer-Hexamer Tilt Angle Preference

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.

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

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.

Molecular Dynamics Simulations in Protein-Protein Docking

Concerted interactions between all the cell components form the basis of biological processes. Protein-protein interactions (PPIs) constitute a tremendous part of this interaction network. Deeper insight into PPIs can help us better understand numerous diseases and lead to the development of new diagnostic and therapeutic strategies. PPI interfaces, until recently, were considered undruggable. However, it is now believed that the interfaces contain "hot spots," which could be targeted by small molecules. Such a strategy would require high-quality structural data of PPIs, which are difficult to obtain experimentally. Therefore, in silico modeling can complement or be an alternative to in vitro approaches. There are several computational methods for analyzing the structural data of the binding partners and modeling of the protein-protein dimer/oligomer structure. The major problem with in silico structure prediction of protein assemblies is obtaining sufficient sampling of protein dynamics. One of the methods that can take protein flexibility and the effects of the environment into account is Molecular Dynamics (MD). While sampling of the whole protein-protein association process with plain MD would be computationally expensive, there are several strategies to harness the method to PPI studies while maintaining reasonable use of resources. This chapter reviews known applications of MD in the PPI investigation workflows.

Mechanical fatigue testing in silico: Dynamic evolution of material properties of nanoscale biological particles

Biological particles have evolved to possess mechanical characteristics necessary to carry out their functions. We developed a computational approach to "fatigue testing in silico", in which constant-amplitude cyclic loading is applied to a particle to explore its mechanobiology. We used this approach to describe dynamic evolution of nanomaterial properties and low-cycle fatigue in the thin spherical encapsulin shell, thick spherical Cowpea Chlorotic Mottle Virus (CCMV) capsid, and thick cylindrical microtubule (MT) fragment over 20 cycles of deformation. Changing structures and force-deformation curves enabled us to describe their damage-dependent biomechanics (strength, deformability, stiffness), thermodynamics (released and dissipated energies, enthalpy, and entropy) and material properties (toughness). Thick CCMV and MT particles experience material fatigue due to slow recovery and damage accumulation over 3-5 loading cycles; thin encapsulin shells show little fatigue due to rapid remodeling and limited damage. The results obtained challenge the existing paradigm: damage in biological particles is partially reversible owing to particle's partial recovery; fatigue crack may or may not grow with each loading cycle and may heal; and particles adapt to deformation amplitude and frequency to minimize the energy dissipated. Using crack size to quantitate damage is problematic as several cracks might form simultaneously in a particle. Dynamic evolution of strength, deformability, and stiffness, can be predicted by analyzing the cycle number (N) dependent damage, [Formula: see text] , where α is a power law and Nf is fatigue life. Fatigue testing in silico can now be used to explore damage-induced changes in the material properties of other biological particles. STATEMENT OF SIGNIFICANCE: Biological particles possess mechanical characteristics necessary to perform their functions. We developed "fatigue testing in silico" approach, which employes Langevin Dynamics simulations of constant-amplitude cyclic loading of nanoscale biological particles, to explore dynamic evolution of the mechanical, energetic, and material properties of the thin and thick spherical particles of encapsulin and Cowpea Chlorotic Mottle Virus, and the microtubule filament fragment. Our study of damage growth and fatigue development challenge the existing paradigm. Damage in biological particles is partially reversible as fatigue crack might heal with each loading cycle. Particles adapt to deformation amplitude and frequency to minimize energy dissipation. The evolution of strength, deformability, and stiffness, can be accurately predicted by analyzing the damage growth in particle structure.

Understanding Virus Structure and Dynamics through Molecular Simulations

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.

Curvature of the Retroviral Capsid Assembly Is Modulated by a Molecular Switch

During the maturation step, the retroviral capsid proteins (CAs) assemble into polymorphic capsids. Their acute curvature is largely determined by 12 pentamers inserted into the hexameric lattice. However, how the CA switches its conformation to control assembly curvature remains unclear. We report the high-resolution structural model of the Rous sarcoma virus (RSV) CA T = 1 capsid, established by molecular dynamics simulations combining solid-state NMR and prior cryoelectron tomography restraints. Comparing this with our previous model of the RSV CA tubular assembly, we identify the key residues for dictating the incorporation of acute curvatures. These residues undergo large torsion angle changes, resulting in a 34° rotation of the C-terminal domain relative to its N-terminal domain around the flexible interdomain linker, without substantial changes of either the conformation of individual domains or the assembly contact interfaces. This knowledge provides new insights to help decipher the mechanism of the retroviral capsid assembly.

Visualizing HIV-1 Capsid and Its Interactions with Antivirals and Host Factors

Understanding of the construction and function of the HIV capsid has advanced considerably in the last decade. This is due in large part to the development of more sophisticated structural techniques, particularly cryo-electron microscopy (cryoEM) and cryo-electron tomography (cryoET). The capsid is known to be a pleomorphic fullerene cone comprised of capsid protein monomers arranged into 200-250 hexamers and 12 pentamers. The latter of these induce high curvature necessary to close the cone at both ends. CryoEM/cryoET, NMR, and X-ray crystallography have collectively described these interactions to atomic or near-atomic resolutions. Further, these techniques have helped to clarify the role the HIV capsid plays in several parts of the viral life cycle, from reverse transcription to nuclear entry and integration into the host chromosome. This includes visualizing the capsid bound to host factors. Multiple proteins have been shown to interact with the capsid. Cyclophilin A, nucleoporins, and CPSF6 promote viral infectivity, while MxB and Trim5α diminish the viral infectivity. Finally, structural insights into the intra- and intermolecular interactions that govern capsid function have enabled development of small molecules, peptides, and truncated proteins to disrupt or stabilize the capsid to inhibit HIV replication. The most promising of these, GS6207, is now in clinical trial.

Switch from Polymorphic to Homogenous Self-Assembly of Virus-Like Particles of Simian Virus 40 through Double-Cysteine Substitution

Self-assembled virus-like particles (VLPs) hold great potential as natural nanomaterials for applications in many fields. For such purposes, monodisperse size distribution is a desirable property. However, the VLPs of simian virus 40 (SV40), a representative VLP platform, are characterized by polymorphism. In an attempt to eliminate the polymorphism, 15 mutants of the VLP subunit (VP1) are constructed through the substitution of double cysteines at the VP1 pentamer interfaces, generating a group of VLPs with altered size distributions. One of the mutants, SS2 (L102C/P300C), specifically forms homogenous T = 1-like tiny VLPs of 24 ± 3 nm in diameter. Moreover, the stability of the SS2 VLPs is markedly enhanced compared with that of wild-type VLPs. The homogeneous self-assembly and stability enhancement of SS2 VLPs can be attributed to the new disulfide bonds contributed by Cys102 and Cys300, which are identified by mass spectrometry and explored by molecular dynamics simulations. Endocytosis inhibition assays indicate that SS2 VLPs, like the polymorphic wild-type VLPs, preserve the multipathway feature of cellular uptake. SS2 VLPs may serve as an evolved version of SV40 VLPs in future studies and applications. The findings of this work would be useful for the design and fabrication of VLP-based materials and devices.

A Multiscale Model for the Self-Assembly of Coat Proteins in Bacteriophage MS2

The self-assembly of viral capsids is an essential step to the formation of infectious viruses. Elucidating the kinetic mechanisms of how a capsid or virus-like particle assembles could advance our knowledge about the viral lifecycle, as well as the general principles in self-assembly of biomaterials. However, current understanding of capsid assembly remains incomplete for many viruses due to the fact that the transient intermediates along the assembling pathways are experimentally difficult to be detected. In this paper, we constructed a new multiscale computational framework to simulate the self-assembly of virus-like particles. We applied our method to the coat proteins of bacteriophage MS2 as a specific model system. This virus-like particle of bacteriophage MS2 has a unique feature that its 90 sequence-identical dimers can be classified into two structurally various groups: one is the symmetric CC dimer, and the other is the asymmetric AB dimer. The homotypic interactions between AB dimers result in a 5-fold symmetric contact, while the heterotypic interactions between AB and CC dimers result in 6-fold symmetric contact. We found that the assembly can be described as a physical process of phase transition that is regulated by various factors such as concentration and specific stoichiometry between AB and CC dimers. Our simulations also demonstrate that heterotypic and homotypic interfaces play distinctive roles in modulating the assembling kinetics. The interaction between AB and CC dimers is much more dynamic than that between two AB dimers. We therefore suggest that the alternate growth of viral capsid through the heterotypic dimer interactions dominates the assembling pathways. This is, to the best of our knowledge, the first multiscale model to simulate the assembling process of coat proteins in bacteriophage MS2. The generality of this approach opens the door to its further applications in assembly of other viral capsids, virus-like particles, and novel drug delivery systems.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Off-Pathway Assembly: A Broad-Spectrum Mechanism of Action for Drugs That Undermine Controlled HIV-1 Viral Capsid Formation

The early and late stages of human immunodeficiency virus (HIV) replication are orchestrated by the capsid (CA) protein, which self-assembles into a conical protein shell during viral maturation. Small molecule drugs known as capsid inhibitors (CIs) impede the highly regulated activity of CA. Intriguingly, a few CIs, such as PF-3450074 (PF74) and GS-CA1, exhibit effects at multiple stages of the viral lifecycle at effective concentrations in the pM to nM regimes, while the majority of CIs target a single stage of the viral lifecycle and are effective at nM to μM concentrations. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanisms that enable CIs to have such curious broad-spectrum activity. Our quantitatively analyzed findings show that CIs can have a profound impact on the hierarchical self-assembly of CA by perturbing populations of small CA oligomers. The self-assembly process is accelerated by the emergence of alternative assembly pathways that favor the rapid incorporation of CA pentamers, and leads to increased structural pleomorphism in mature capsids. Two relevant phenotypes are observed: (1) eccentric capsid formation that may fail to encase the viral genome and (2) rapid disassembly of the capsid, which express at late and early stages of infection, respectively. Finally, our study emphasizes the importance of adopting a dynamical perspective on inhibitory mechanisms and provides a basis for the design of future therapeutics that are effective at low stoichiometric ratios of drug to protein.

Novel Intersubunit Interaction Critical for HIV-1 Core Assembly Defines a Potentially Targetable Inhibitor Binding Pocket

Precise assembly and disassembly of the HIV-1 capsid core are key to the success of viral replication. The forces that govern capsid core formation and dissociation involve intricate interactions between pentamers and hexamers formed by HIV-1 CA. We identified one particular interaction between E28 of one CA and K30′ of the adjacent CA that appears more frequently in pentamers than in hexamers and that is important for capsid assembly. Targeting the corresponding site could lead to the development of antivirals which disrupt this interaction and affect capsid assembly.

HIV-1 capsid protein (CA) plays critical roles in both early and late stages of the viral replication cycle. Mutagenesis and structural experiments have revealed that capsid core stability significantly affects uncoating and initiation of reverse transcription in host cells. This has led to efforts in developing antivirals targeting CA and its assembly, although none of the currently identified compounds are used in the clinic for treatment of HIV infection. A specific interaction that is primarily present in pentameric interfaces in the HIV-1 capsid core was identified and is reported to be important for CA assembly. This is shown by multidisciplinary characterization of CA site-directed mutants using biochemical analysis of virus-like particle formation, transmission electron microscopy of in vitro assembly, crystallographic studies, and molecular dynamic simulations. The data are consistent with a model where a hydrogen bond between CA residues E28 and K30′ from neighboring N-terminal domains (CA_NTDs) is important for CA pentamer interactions during core assembly. This pentamer-preferred interaction forms part of an N-terminal domain interface (NDI) pocket that is amenable to antiviral targeting.

Advances in coarse-grained modeling of macromolecular complexes

Recent progress in coarse-grained (CG) molecular modeling and simulation has facilitated an influx of computational studies on biological macromolecules and their complexes. Given the large separation of length-scales and time-scales that dictate macromolecular biophysics, CG modeling and simulation are well-suited to bridge the microscopic and mesoscopic or macroscopic details observed from all-atom molecular simulations and experiments, respectively. In this review, we first summarize recent innovations in the development of CG models, which broadly include structure-based, knowledge-based, and dynamics-based approaches. We then discuss recent applications of different classes of CG models to explore various macromolecular complexes. Finally, we conclude with an outlook for the future in this ever-growing field of biomolecular modeling.

Modeling the Self-Assembly of Protein Complexes through a Rigid-Body Rotational Reaction-Diffusion Algorithm

The reaction-diffusion equations provide a powerful framework for modeling nonequilibrium, cell-scale dynamics over the long time scales that are inaccessible by traditional molecular modeling approaches. Single-particle reaction-diffusion offers the highest resolution technique for tracking such dynamics, but it has not been applied to the study of protein self-assembly due to its treatment of reactive species as single-point particles. Here, we develop a relatively simple but accurate approach for building rigid structure and rotation into single-particle reaction-diffusion methods, providing a rate-based method for studying protein self-assembly. Our simplifying assumption is that reactive collisions can be evaluated purely on the basis of the separations between the sites, and not their orientations. The challenge of evaluating reaction probabilities can then be performed using well-known equations based on translational diffusion in both 3D and 2D, by employing an effective diffusion constant we derive here. We show how our approach reproduces both the kinetics of association, which is altered by rotational diffusion, and the equilibrium of reversible association, which is not. Importantly, the macroscopic kinetics of association can be predicted on the basis of the microscopic parameters of our structurally resolved model, allowing for critical comparisons with theory and other rate-based simulations. We demonstrate this method for efficient, rate-based simulations of self-assembly of clathrin trimers, highlighting how formation of regular lattices impacts the kinetics of association.

An isolate of Potato Virus X capsid protein from N. benthamiana: Insights from homology modeling and molecular dynamics simulation

Since Potato Virus X (PVX) is easily transmitted mechanically between their hosts, its control is difficult. We have previously reported new isolate of this virus (PVX-Iran, GenBank Accession number FJ461343). However, the molecular basis of resistance breaking activity and its relation to capsid protein structure are still not well-understood. SDS-PAGE, ELISA, Western blot and RT-PCR molecular examinations were performed on the inoculated plants Nicotiana benthamiana. The pathological symptoms were related to the PVX isolate. The capsid protein (CP) structure were modeled based on homology and subjected to three independent 80 ns molecular dynamics minimization (GROMACS, OPLS force field) in the SPC water box. The RMSD, RMSF, SASA, and electrostatic properties were retrieved from the trajectories. Flexibility and hydrophilic nature of the N-terminal residues (1-34) of solvated CP could be observed in conformational changes upon minimization. The obtained structure was then docked with NbPCIP1 using ClusPro 2.0. The strong binding affinity of these two proteins (≈-16.0 Kcal mol-1) represents the formation of inclusion body and hence appearance of the symptoms.

Physical properties of the HIV-1 capsid from all-atom molecular dynamics simulations

Human immunodeficiency virus type 1 (HIV-1) infection is highly dependent on its capsid. The capsid is a large container, made of ∼1,300 proteins with altogether 4 million atoms. Although the capsid proteins are all identical, they nevertheless arrange themselves into a largely asymmetric structure made of hexamers and pentamers. The large number of degrees of freedom and lack of symmetry pose a challenge to studying the chemical details of the HIV capsid. Simulations of over 64 million atoms for over 1 μs allow us to conduct a comprehensive study of the chemical–physical properties of an empty HIV-1 capsid, including its electrostatics, vibrational and acoustic properties, and the effects of solvent (ions and water) on the capsid. The simulations reveal critical details about the capsid with implications to biological function.

The large and complex HIV-1 capsid modulates several molecular events during HIV-1’s infective cycle. Here the authors use all-atom molecular dynamic simulations to probe the biophysical properties of the genome-free HIV-1 capsid.

Parameter Optimization for Interaction between C-Terminal Domains of HIV-1 Capsid Protein

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.

In vitro protease cleavage and computer simulations reveal the HIV-1 capsid maturation pathway

HIV-1 virions assemble as immature particles containing Gag polyproteins that are processed by the viral protease into individual components, resulting in the formation of mature infectious particles. There are two competing models for the process of forming the mature HIV-1 core: the disassembly and de novo reassembly model and the non-diffusional displacive model. To study the maturation pathway, we simulate HIV-1 maturation in vitro by digesting immature particles and assembled virus-like particles with recombinant HIV-1 protease and monitor the process with biochemical assays and cryoEM structural analysis in parallel. Processing of Gag in vitro is accurate and efficient and results in both soluble capsid protein and conical or tubular capsid assemblies, seemingly converted from immature Gag particles. Computer simulations further reveal probable assembly pathways of HIV-1 capsid formation. Combining the experimental data and computer simulations, our results suggest a sequential combination of both displacive and disassembly/reassembly processes for HIV-1 maturation.

Two competing models—disassembly/reassembly and displacive—have been proposed for how immature spherical HIV virions transform into mature particles with conical cores. Here the authors provide evidence that both disassembly/reassembly and displacive processes occur sequentially during the maturation process.

Reaction-diffusion basis of retroviral infectivity

Retrovirus particle (virion) infectivity requires diffusion and clustering of multiple transmembrane envelope proteins (Env₃) 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 Env₃ 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, Env₃ diffuses comparably to Gag-absent Env₃ Initial immobility of Env₃ 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 Env₃ diffusion, and permits Env₃ 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'.

Coarse-grained simulation reveals key features of HIV-1 capsid self-assembly

The maturation of HIV-1 viral particles is essential for viral infectivity. During maturation, many copies of the capsid protein (CA) self-assemble into a capsid shell to enclose the viral RNA. The mechanistic details of the initiation and early stages of capsid assembly remain to be delineated. We present coarse-grained simulations of capsid assembly under various conditions, considering not only capsid lattice self-assembly but also the potential disassembly of capsid upon delivery to the cytoplasm of a target cell. The effects of CA concentration, molecular crowding, and the conformational variability of CA are described, with results indicating that capsid nucleation and growth is a multi-stage process requiring well-defined metastable intermediates. Generation of the mature capsid lattice is sensitive to local conditions, with relatively subtle changes in CA concentration and molecular crowding influencing self-assembly and the ensemble of structural morphologies.

Significant morphological changes occur during the conversion of the immature HIV virion into a mature infectious form. Here the authors use coarse-grained molecular dynamics simulations to model HIV-1 capsid self-assembly and disassembly events that suggests several metastable capsid intermediates sensitive to local conditions.

HIV Capsid Assembly, Mechanism, and Structure

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.

Major Variations in HIV-1 Capsid Assembly Morphologies Involve Minor Variations in Molecular Structures of Structurally Ordered Protein Segments

We present the results of solid state nuclear magnetic resonance (NMR) experiments on HIV-1 capsid protein (CA) assemblies with three different morphologies, namely wild-type CA (WT-CA) tubes with 35-60 nm diameters, planar sheets formed by the Arg(18)-Leu mutant (R18L-CA), and R18L-CA spheres with 20-100 nm diameters. The experiments are intended to elucidate molecular structural variations that underlie these variations in CA assembly morphology. We find that multidimensional solid state NMR spectra of (15)N,(13)C-labeled CA assemblies are remarkably similar for the three morphologies, with only small differences in (15)N and (13)C chemical shifts, no significant differences in NMR line widths, and few differences in the number of detectable NMR cross-peaks. Thus, the pronounced differences in morphology do not involve major differences in the conformations and identities of structurally ordered protein segments. Instead, morphological variations are attributable to variations in conformational distributions within disordered segments, which do not contribute to the solid state NMR spectra. Variations in solid state NMR signals from certain amino acid side chains are also observed, suggesting differences in the intermolecular dimerization interface between curved and planar CA lattices, as well as possible differences in intramolecular helix-helix packing.

Construction of a novel coarse grain model for simulations of HIV capsid assembly to capture the backbone structure and inter-domain motions in solution

We show the construction of a novel coarse grain model for simulations of HIV capsid assembly based on four structural models of HIV capsid proteins: isolated hexamer 3H47.pdb, tubular assembly 3J34.pdb, isolated pentamer 3P05.pdb and C-terminus dimer 2KOD.pdb. The data demonstrates the derivation of inter-domain motions from all atom Molecular Dynamics simulations and comparison with the motions derived from the analysis of solution NMR results defined in 2M8L.pdb. Snapshots from a representative Monte Carlo simulation with 128 dimeric subunit proteins based on 3J34.pdb are shown in addition to the quantitative analysis of its assembly pathway. Movies of the assembly process are compiled with snapshots of representative simulations of four structural models. The methods and data in this article were utilized in Qiao et al. (in press) [1] to probe the mechanism of polymorphism and curvature control of HIV capsid assembly.