Modeling HIV-1 viral capsid nucleation by dynamical systems

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.

A New Model System for Exploring Assembly Mechanisms of the HIV-1 Immature Capsid In Vivo

The assembly of the HIV-1 immature capsid (HIC) is an essential step in the virus life cycle. In vivo, the HIC is composed of [Formula: see text] hexameric building blocks, and it takes 5-6 min to complete the assembly process. The involvement of numerous building blocks and the rapid timecourse makes it difficult to understand the HIC assembly process. In this work, we study HIC assembly in vivo by using differential equations. We first obtain a full model with 420 differential equations. Then, we reduce six addition reactions for separate building blocks to a single complex reaction. This strategy reduces the full model to 70 equations. Subsequently, the theoretical analysis of the reduced model shows that it might not be an effective way to decrease the HIC concentration at the equilibrium state by decreasing the microscopic on-rate constants. Based on experimental data, we estimate that the nucleating structure is much smaller than the HIC. We also estimate that the microscopic on-rate constant for nucleation reactions is far less than that for elongation reactions. The parametric collinearity investigation testifies the reliability of these two characteristics, which might explain why free building blocks do not readily polymerize into higher-order polymers until their concentration reaches a threshold value. These results can provide further insight into the assembly mechanisms of the HIC in vivo.

Evolution of Coagulation-Fragmentation Stochastic Processes Using Accurate Chemical Master Equation Approach

Coagulation and fragmentation (CF) is a fundamental process in which smaller particles attach to each other to form larger clusters while existing clusters break up into smaller particles . It is a ubiquitous process that plays important roles in many physical and biological phenomena. CF is typically a stochastic process that often occurs in confined spaces with a limited number of available particles . Here, we study the CF process formulated with the discrete Chemical Master Equation (dCME). Using the newly developed Accurate Chemical Master Equation (ACME) method, we examine the time-dependent behavior of the CF system. We investigate the effects of a number of important factors that influence the overall behavior of the system, including the dimensionality, the ratio of attachment to detachment rates among clusters, and the initial conditions. By comparing CF in one and three dimensions, we conclude that systems in three dimensions are more likely to form large clusters. We also demonstrate how the ratio of the attachment to detachment rates affects the dynamics and the steady-state of the system. Finally, we demonstrate the relationship between the formation of large clusters and the initial condition.

Mathematical modeling and quantitative analysis of HIV-1 Gag trafficking and polymerization

Gag, as the major structural protein of HIV-1, is necessary for the assembly of the HIV-1 sphere shell. An in-depth understanding of its trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. We developed a mathematical model to simulate two biophysical processes, specifically Gag monomer and dimer transport in the cytoplasm and the polymerization of monomers to form a hexamer underneath the plasma membrane. Using experimental data, an optimization approach was utilized to identify the model parameters, and the identifiability and sensitivity of these parameters were then analyzed. Using our model, we analyzed the weight of the pathways involved in the polymerization reactions and concluded that the predominant pathways for the formation of a hexamer might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. We then deduced that the dimer and trimer intermediates might be crucial in hexamer formation. We also explored four theoretical combined methods for Gag suppression, and hypothesized that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. This work serves as a guide for future theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization, and might help accelerate the efficiency of anti-AIDS drug design.

The human immunodeficiency virus (HIV-1) is a retrovirus that causes acquired immunodeficiency syndrome (AIDS), an infectious disease with high annual mortality. Gag protein is the major structural protein of HIV-1 and can self-assemble into the HIV-1 sphere shell. Therefore, an in-depth understanding of Gag protein trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. Through mathematical modeling, optimization and quantitative analysis, we hypothesized the budding and release time of virus-like particles and revealed that the dimer and trimer intermediates might be crucial in hexamer formation. We also concluded that the predominant pathways in hexamer formation might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. Our analysis also suggested that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. These results serve as a guide for further theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization and could aid anti-AIDS drug design.

In vitro assembly of the Rous Sarcoma Virus capsid protein into hexamer tubes at physiological temperature

During a proteolytically-driven maturation process, the orthoretroviral capsid protein (CA) assembles to form the convex shell that surrounds the viral genome. In some orthoretroviruses, including Rous Sarcoma Virus (RSV), CA carries a short and hydrophobic spacer peptide (SP) at its C-terminus early in the maturation process, which is progressively removed as maturation proceeds. In this work, we show that RSV CA assembles in vitro at near-physiological temperatures, forming hexamer tubes that effectively model the mature capsid surface. Tube assembly is strongly influenced by electrostatic effects, and is a nucleated process that remains thermodynamically favored at lower temperatures, but is effectively arrested by the large Gibbs energy barrier associated with nucleation. RSV CA tubes are multi-layered, being formed by nested and concentric tubes of capsid hexamers. However the spacer peptide acts as a layering determinant during tube assembly. If only a minor fraction of CA-SP is present, multi-layered tube formation is blocked, and single-layered tubes predominate. This likely prevents formation of biologically aberrant multi-layered capsids in the virion. The generation of single-layered hexamer tubes facilitated 3D helical image reconstruction from cryo-electron microscopy data, revealing the basic tube architecture.

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.