Coarse-grained protein-protein stiffnesses and dynamics from all-atom simulations

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.

Antiviral compounds modulate elasticity, strength and material fatigue of a virus capsid framework

This study investigates whether the biochemical and antiviral effects of organic compounds that bind different sites in the mature human immunodeficiency virus capsid may be related to the modulation of different mechanical properties of the protein lattice from which the capsid is built. Mechanical force was used as a probe to quantify, in atomic force microscopy experiments at physiological pH and ionic strength, ligand-mediated changes in capsid lattice elasticity, breathing, strength against local dislocation by mechanical stress, and resistance to material fatigue. The results indicate that the effects of the tested compounds on assembly or biochemical stability can be linked, from a physics-based perspective, to their interference with the mechanical behavior of the viral capsid framework. The antivirals CAP-1 and CAI-55 increased the intrinsic elasticity and breathing of the capsid protein lattice and may entropically decrease the probability of the capsid protein to assemble into a functionally competent conformation. Antiviral PF74 increased the resistance of the capsid protein lattice to disruption by mechanical stress and material fatigue and may enthalpically strengthen the basal capsid lattice against breakage and disintegration. This study provides proof of concept that the interrogation of the mechanical properties of the nanostructured protein material that makes a virus capsid may provide fundamental insights into the biophysical action of capsid-binding antiviral agents. The implications for drug design by specifically targeting the biomechanics of viruses are discussed.

A Genetically Engineered, Chain Mail-Like Nanostructured Protein Material with Increased Fatigue Resistance and Enhanced Self-Healing

Protein-based nanostructured materials are being developed for many biomedical and nanotechnological applications. Despite their many desirable features, protein materials are highly susceptible to disruption by mechanical stress and fatigue. This study is aimed to increase fatigue resistance and enhance self-healing of a natural protein-based supramolecular nanomaterial through permanent genetic modification. The authors envisage the conversion of a model nanosheet, formed by a regular array of noncovalently bound human immunodeficiency virus capsid protein molecules, into a supramolecular "chain mail." Rationally engineered mutations allow the formation of a regular network of disulfide bridges in the protein lattice. This network links each molecule in the lattice to each adjacent molecule through one covalent bond, analogous to the rivetting of interlinked iron rings in the chain mail of a medieval knight. The engineered protein nanosheet shows greatly increased thermostability and resistance to mechanical stress and fatigue in particular, as well as enhanced self-healing, without undesirable stiffening compared to the original material. The results provide proof of concept for a genetic design to permanently increase fatigue resistance and enhance self-healing of protein-based nanostructured materials. They also provide insights into the molecular basis for fatigue of protein materials.

Elastic continuum stiffness of contractile tail sheaths from molecular dynamics simulations

Contractile tails are key components of the biological nanomachinery involved in cell membrane puncturing, where they provide a means to deliver molecules and ions inside cells. Two intriguing examples of contractile tails are those from bacteriophage T4 and R2-pyocin. Although the two systems are different in terms of biological activity, they share a fascinatingly similar injection mechanism, during which the tail sheaths of both systems contract from a so-called extended state to around half of their length (the contracted state), accompanied by release of elastic energy originally stored in the sheath. Despite the great prevalence and biomedical importance of contractile delivery systems, many fundamental details of their injection machinery and dynamics are still unknown. In this work, we calculate the bending and torsional stiffness constants of a helical tail sheath strand of bacteriophage T4 and R2-pyocin, in both extended and contracted states, using molecular dynamics simulations of about one-sixth of the entire sheath. Differences in stiffness constants between the two systems are rationalized by comparing their all-atom monomer structures, changes in sheath architecture on contraction, and differences in interstrand interactions. The calculated coefficients indicate that the T4 strand is stiffer for both bending and torsion than the corresponding R2-pyocin strands in both extended and contracted conformations. The sheath strands also have greater stiffness in the contracted state for both systems. As the main application of this study, we describe how the stiffness constants can be incorporated in a model to simulate the dynamics of contractile nanoinjection machineries.

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.

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.

Quantification and modification of the equilibrium dynamics and mechanics of a viral capsid lattice self-assembled as a protein nanocoating

Self-assembling, protein-based bidimensional lattices are being developed as functionalizable, highly ordered biocoatings for multiple applications in nanotechnology and nanomedicine. Unfortunately, protein assemblies are soft materials that may be too sensitive to mechanical disruption, and their intrinsic conformational dynamism may also influence their applicability. Thus, it may be critically important to characterize, understand and manipulate the mechanical features and dynamic behavior of protein assemblies in order to improve their suitability as nanomaterials. In this study, the capsid protein of the human immunodeficiency virus was induced to self-assemble as a continuous, single layered, ordered nanocoating onto an inorganic substrate. Atomic force microscopy (AFM) was used to quantify the mechanical behavior and the equilibrium dynamics ("breathing") of this virus-based, self-assembled protein lattice in close to physiological conditions. The results uniquely provided: (i) evidence that AFM can be used to directly visualize in real time and quantify slow breathing motions leading to dynamic disorder in protein nanocoatings and viral capsid lattices; (ii) characterization of the dynamics and mechanics of a viral capsid lattice and protein-based nanocoating, including flexibility, mechanical strength and remarkable self-repair capacity after mechanical damage; (iii) proof of principle that chemical additives can modify the dynamics and mechanics of a viral capsid lattice or protein-based nanocoating, and improve their applied potential by increasing their mechanical strength and elasticity. We discuss the implications for the development of mechanically resistant and compliant biocoatings precisely organized at the nanoscale, and of novel antiviral agents acting on fundamental physical properties of viruses.

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

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.

Modeling Viral Capsid Assembly

I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Structural systems biology and multiscale signaling models

We review current advances in experimental as well as computational modeling and simulation approaches to structural systems biology, whose overall aim is to build quantitative models of signaling networks while retaining the crucial elements of molecular specificity. We briefly discuss the current and emerging experimental and computational methods, particularly focusing on hybrid and multiscale methods, and highlight several applications in cell signaling with quantitative and predictive capabilities. The scope of such models range from delineating protein-protein interactions to describing clinical implications.

Coarse-graining of multiprotein assemblies

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.