Atomic Modeling of an Immature Retroviral Lattice Using Molecular Dynamics and Mutagenesis

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

All-Atom MD Simulations of the HBV Capsid Complexed with AT130 Reveal Secondary and Tertiary Structural Changes and Mechanisms of Allostery

The hepatitis B virus (HBV) capsid is an attractive drug target, relevant to combating viral hepatitis as a major public health concern. Among small molecules known to interfere with capsid assembly, the phenylpropenamides, including AT130, represent an important antiviral paradigm based on disrupting the timing of genome packaging. Here, all-atom molecular dynamics simulations of an intact AT130-bound HBV capsid reveal that the compound increases spike flexibility and improves recovery of helical secondary structure in the spike tips. Regions of the capsid-incorporated dimer that undergo correlated motion correspond to established sub-domains that pivot around the central chassis. AT130 alters patterns of correlated motion and other essential dynamics. A new conformational state of the dimer is identified, which can lead to dramatic opening of the intradimer interface and disruption of communication within the spike tip. A novel salt bridge is also discovered, which can mediate contact between the spike tip and fulcrum even in closed conformations, revealing a mechanism of direct communication across these sub-domains. Altogether, results describe a dynamical connection between the intra- and interdimer interfaces and enable mapping of allostery traversing the entire core protein dimer.

Design and exploration of C-3 benzoic acid bioisosteres and alkyl replacements in the context of GSK3532795 (BMS-955176) that exhibit broad spectrum HIV-1 maturation inhibition

GSK3532795 (formerly BMS-955176) is a second-generation HIV-1 maturation inhibitor that has shown broad spectrum antiviral activity and preclinical PK predictive of once-daily dosing in humans. Although efficacy was confirmed in clinical trials, the observation of gastrointestinal intolerability and the emergence of drug resistant virus in a Phase 2b clinical study led to the discontinuation of GSK3532795. As part of the effort to further map the maturation inhibitor pharmacophore and provide additional structural options, the evaluation of alternates to the C-3 phenyl substituent in this chemotype was pursued. A cyclohexene carboxylic acid provided exceptional inhibition of wild-type, V370A and ΔV370 mutant viruses in addition to a suitable PK profile following oral dosing to rats. In addition, a novel spiro[3.3]hept-5-ene was designed to extend the carboxylic acid further from the triterpenoid core while reducing side chain flexibility compared to the other alkyl substituents. This modification was shown to closely emulate the C-3 benzoic acid moiety of GSK3532795 from both a potency and PK perspective, providing a non-traditional, sp³-rich bioisostere of benzene. Herein, we detail additional modifications to the C-3 position of the triterpenoid core that offer effective replacements for the benzoic acid of GSK3532795 and capture the interplay between these new C-3 elements and C-17 modifications that contribute to enhanced polymorph coverage.

The impact of Gag non-cleavage site mutations on HIV-1 viral fitness from integrative modelling and simulations

The high mutation rate in retroviruses is one of the leading causes of drug resistance. In human immunodeficiency virus type-1 (HIV-1), synergistic mutations in its protease and the protease substrate - the Group-specific antigen (Gag) polyprotein - work together to confer drug resistance against protease inhibitors and compensate the mutations affecting viral fitness. Some Gag mutations can restore Gag-protease binding, yet most Gag-protease correlated mutations occur outside of the Gag cleavage site. To investigate the molecular basis for this, we now report multiscale modelling approaches to investigate various sequentially cleaved Gag products in the context of clinically relevant mutations that occur outside of the cleavage sites, including simulations of the largest Gag proteolytic product in its viral membrane-bound state. We found that some mutations, such as G123E and H219Q, involve direct interaction with cleavage site residues to influence their local environment, while certain mutations in the matrix domain lead to the enrichment of lipids important for Gag targeting and assembly. Collectively, our results reveal why non-cleavage site mutations have far-reaching implications outside of Gag proteolysis, with important consequences for drugging Gag maturation intermediates and tackling protease inhibitor resistance.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Structures of immature EIAV Gag lattices reveal a conserved role for IP6 in lentivirus assembly

Retrovirus assembly is driven by the multidomain structural protein Gag. Interactions between the capsid domains (CA) of Gag result in Gag multimerization, leading to an immature virus particle that is formed by a protein lattice based on dimeric, trimeric, and hexameric protein contacts. Among retroviruses the inter- and intra-hexamer contacts differ, especially in the N-terminal sub-domain of CA (CANTD). For HIV-1 the cellular molecule inositol hexakisphosphate (IP6) interacts with and stabilizes the immature hexamer, and is required for production of infectious virus particles. We have used in vitro assembly, cryo-electron tomography and subtomogram averaging, atomistic molecular dynamics simulations and mutational analyses to study the HIV-related lentivirus equine infectious anemia virus (EIAV). In particular, we sought to understand the structural conservation of the immature lentivirus lattice and the role of IP6 in EIAV assembly. Similar to HIV-1, IP6 strongly promoted in vitro assembly of EIAV Gag proteins into virus-like particles (VLPs), which took three morphologically highly distinct forms: narrow tubes, wide tubes, and spheres. Structural characterization of these VLPs to sub-4Å resolution unexpectedly showed that all three morphologies are based on an immature lattice with preserved key structural components, highlighting the structural versatility of CA to form immature assemblies. A direct comparison between EIAV and HIV revealed that both lentiviruses maintain similar immature interfaces, which are established by both conserved and non-conserved residues. In both EIAV and HIV-1, IP6 regulates immature assembly via conserved lysine residues within the CACTD and SP. Lastly, we demonstrate that IP6 stimulates in vitro assembly of immature particles of several other retroviruses in the lentivirus genus, suggesting a conserved role for IP6 in lentiviral assembly.

Multiscale modelling and simulation of viruses

In recent years, advances in structural biology, integrative modelling, and simulation approaches have allowed us to gain unprecedented insights into viral structure and dynamics. In this article we survey recent studies utilizing this wealth of structural information to build computational models of partial or complete viruses and to elucidate mechanisms of viral function. Additionally, the close interplay of viral pathogens with host factors - such as cellular and intracellular membranes, receptors, antibodies, and other host proteins - makes accurate models of viral interactions and dynamics essential. As viruses continue to pose severe challenges in prevention and treatment, enhancing our mechanistic understanding of viral infection is vital to enable the development of novel therapeutic strategies.

Structural analysis of pleomorphic and asymmetric viruses using cryo-electron tomography and subtomogram averaging

Describing the protein interactions that form pleomorphic and asymmetric viruses represents a considerable challenge to most structural biology techniques, including X-ray crystallography and single particle cryo-electron microscopy. Obtaining a detailed understanding of these interactions is nevertheless important, considering the number of relevant human pathogens that do not follow strict icosahedral or helical symmetry. Cryo-electron tomography and subtomogram averaging methods provide structural insights into complex biological environments and are well suited to go beyond structures of perfectly symmetric viruses. This chapter discusses recent developments showing that cryo-ET and subtomogram averaging can provide high-resolution insights into hitherto unknown structural features of pleomorphic and asymmetric virus particles. It also describes how these methods have significantly added to our understanding of retrovirus capsid assemblies in immature and mature viruses. Additional examples of irregular viruses and their associated proteins, whose structures have been studied via cryo-ET and subtomogram averaging, further support the versatility of these methods.

Single-molecule studies of flavivirus envelope dynamics: Experiment and computation

Flaviviruses are simple enveloped viruses exhibiting complex structural and functional heterogeneities. Decades of research have provided crucial basic insights, antiviral medication and moderately successful gene therapy trials. The most infectious particle is, however, not always the most abundant one in a population, questioning the utility of classic ensemble-averaging virology approaches. Indeed, viral replication is often not particularly efficient, prone to errors or containing parallel routes. Here, we review different single-molecule sensitive fluorescence methods that are employed to investigate flaviviruses. In particular, we review how (i) time-resolved Förster resonance energy transfer (trFRET) was applied to probe dengue envelope conformations; (ii) FRET-fluorescence correlation spectroscopy to investigate dengue envelope intrinsic dynamics and (iii) single particle tracking to follow the path of dengue viruses in cells. We also discuss how such methods may be supported by molecular dynamics (MD) simulations over a range of spatio-temporal scales, to provide complementary data on the structure and dynamics of flaviviral systems. We describe recent improvements in multiscale MD approaches that allowed the simulation of dengue particle envelopes in near-atomic resolution. We hope this review is an incentive for setting up and applying similar single-molecule studies and combine them with MD simulations to investigate structural dynamics of entire flavivirus particles over the nanosecond-to-millisecond time-scale and follow viruses during infection in cells over milliseconds to minutes.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

All-atom molecular dynamics of the HBV capsid reveals insights into biological function and cryo-EM resolution limits

The hepatitis B virus capsid represents a promising therapeutic target. Experiments suggest the capsid must be flexible to function; however, capsid structure and dynamics have not been thoroughly characterized in the absence of icosahedral symmetry constraints. Here, all-atom molecular dynamics simulations are leveraged to investigate the capsid without symmetry bias, enabling study of capsid flexibility and its implications for biological function and cryo-EM resolution limits. Simulation results confirm flexibility and reveal a propensity for asymmetric distortion. The capsid’s influence on ionic species suggests a mechanism for modulating the display of cellular signals and implicates the capsid’s triangular pores as the location of signal exposure. A theoretical image reconstruction performed using simulated conformations indicates how capsid flexibility may limit the resolution of cryo-EM. Overall, the present work provides functional insight beyond what is accessible to experimental methods and raises important considerations regarding asymmetry in structural studies of icosahedral virus capsids.

Modeling the dynamics and kinetics of HIV-1 Gag during viral assembly

We report a computational model for the assembly of HIV-1 Gag into immature viral particles at the plasma membrane. To reproduce experimental structural and kinetic properties of assembly, a process occurring on the order of minutes, a coarse-grained representation consisting of a single particle per Gag molecule is developed. The model uses information relating the functional interfaces implicated in Gag assembly, results from cryo electron-tomography, and biophysical measurements from fluorescence microscopy, such as the dynamics of Gag assembly at single virions. These experimental constraints eliminated many classes of potential interactions, and narrowed the model to a single interaction scheme with two non-equivalent interfaces acting to form Gags into a hexamer, and a third interface acting to link hexamers together. This model was able to form into a hexameric structure with correct lattice spacing and reproduced biologically relevant growth rates. We explored the effect of genomic RNA seeding punctum growth, finding that RNA may be a factor in locally concentrating Gags to initiate assembly. The simulation results infer that completion of assembly cannot be governed simply by Gag binding kinetics. However the addition of membrane curvature suggests that budding of the virion from the plasma membrane could factor into slowing incorporation of Gag at an assembly site resulting in virions of the same size and number of Gag molecules independent of Gag concentration or the time taken to complete assembly. To corroborate the results of our simulation model, we developed an analytic model for Gag assembly finding good agreement with the simulation results.

Quenching protein dynamics interferes with HIV capsid maturation

Maturation of HIV-1 particles encompasses a complex morphological transformation of Gag via an orchestrated series of proteolytic cleavage events. A longstanding question concerns the structure of the C-terminal region of CA and the peptide SP1 (CA–SP1), which represents an intermediate during maturation of the HIV-1 virus. By integrating NMR, cryo-EM, and molecular dynamics simulations, we show that in CA–SP1 tubes assembled in vitro, which represent the features of an intermediate assembly state during maturation, the SP1 peptide exists in a dynamic helix–coil equilibrium, and that the addition of the maturation inhibitors Bevirimat and DFH-055 causes stabilization of a helical form of SP1. Moreover, the maturation-arresting SP1 mutation T8I also induces helical structure in SP1 and further global dynamical and conformational changes in CA. Overall, our results show that dynamics of CA and SP1 are critical for orderly HIV-1 maturation and that small molecules can inhibit maturation by perturbing molecular motions.

The process of HIV particle maturation involves complex molecular transitions. Here the authors combine NMR spectroscopy, cryo-EM, and molecular dynamics simulations to provide insight into the conformational equilibria in CA-SP1 assemblies relevant to HIV-1 maturation intermediates formation.

Molecular Mechanism of Flavin Photoprotection by Archaeal Dodecin: Photoinduced Electron Transfer and Mg2+-Promoted Proton Transfer

Photoinduced biochemical reactions are ubiquitously governed by derivatives of flavin, which is a key player in a manifold of cellular redox reactions. The photoreactivity of flavins is also one of their greatest disadvantages as the molecules are sensitive to photodegradation. To prevent this unfavorable reaction, UV-light-exposed archaea bacteria, such as Halobacterium salinarum, manage the task of protecting flavin derivatives by dodecin, a protein which stores flavins and efficiently photoprotects them. In this study, we shed light on the photoprotection mechanism, i.e., the excited state quenching mechanism by dodecin using computational methodology. Molecular dynamics (MD) simulations unraveled the hydrogen bond network in the flavin binding pocket as a starting point for proton transfer upon preceding electron transfer. Using high-level ab initio quantum chemical methods, different proton transfer channels have been investigated and an energetically feasible Mg²⁺-promoted channel has been identified fully explaining previous experimental observations. This is the first extensive theoretical study of archaeal dodecin, furthering the understanding of its photocycle and manipulation.

A Large Size Chimeric Highly Immunogenic Peptide Presents Multistage Plasmodium Antigens as a Vaccine Candidate System against Malaria

Rational strategies for obtaining malaria vaccine candidates should include not only a proper selection of target antigens for antibody stimulation, but also a versatile molecular design based on ordering the right pieces from the complex pathogen molecular puzzle towards more active and functional immunogens. Classical Plasmodium falciparum antigens regarded as vaccine candidates have been selected as model targets in this study. Among all possibilities we have chosen epitopes of PfCSP, STARP; MSA1 and Pf155/RESA from pre- and erythrocyte stages respectively for designing a large 82-residue chimeric immunogen. A number of options aimed at diminishing steric hindrance for synthetic procedures were assessed based on standard Fmoc chemistry such as building block orthogonal ligation; pseudo-proline and microwave-assisted procedures, therefore the large-chimeric target was produced, characterized and immunologically tested. Antigenicity and functional in vivo efficacy tests of the large-chimera formulations administered alone or as antigen mixtures have proven the stimulation of high antibody titers, showing strong correlation with protection and parasite clearance of vaccinated BALB/c mice after being lethally challenged with both P. berghei-ANKA and P. yoelii 17XL malaria strains. Besides, 3D structure features shown by the large-chimera encouraged as to propose using these rational designed large synthetic molecules as reliable vaccine candidate-presenting systems.

All-Atom Molecular Dynamics Simulations of Entire Virus Capsid Reveal the Role of Ion Distribution in Capsid's Stability

Present experimental methods do not have sufficient resolution to investigate all processes in virus particles at atomistic details. We report the results of molecular dynamics simulations and analyze the connection between the number of ions inside an empty capsid of PCV2 virus and its stability. We compare the crystallographic structures of the capsids with unresolved N-termini and without them in realistic conditions (room temperature and aqueous solution) and show that the structure is preserved. We find that the chloride ions play a key role in the stability of the capsid. A low number of chloride ions results in loss of the native icosahedral symmetry, while an optimal number of chloride ions create a neutralizing layer next to the positively charged inner surface of the capsid. Understanding the dependence of the capsid stability on the distribution of the ions will help clarify the details of the viral life cycle that is ultimately connected to the role of packaged viral genome inside the capsid.

Oligomeric Conjugated Polyelectrolytes Display Site-Preferential Binding to an MS2 Viral Capsid

Opportunistic bacteria and viruses are a worldwide health threat prompting the need to develop new targeting modalities. A class of novel synthetic poly(phenylene ethynylene) (PPE)-based oligomeric conjugated polyelectrolytes (OPEs) have demonstrated potent wide-spectrum biocidal activity. A subset of cationic OPEs display high antiviral activity against the MS2 bacteriophage. The oligomers have been found to inactivate the bacteriophage and perturb the morphology of the MS2 viral capsid. However, details of the initial binding and interactions between the OPEs and the viruses are not well understood. In this study, we use a multiscale computational approach, including random sampling, molecular dynamics, and electronic structure calculations, to gain an understanding of the molecular-level interactions of a series of OPEs that vary in length, charge, and functional groups with the MS2 capsid. Our results show that OPEs strongly bind to the MS2 capsid protein assembly with binding energies of up to -30 kcal/mol. Free-energy analysis shows that the binding is dominated by strong van der Waals interactions between the hydrophobic OPE backbone and the capsid surface and strong electrostatic free energy contributions between the OPE charged moieties and charged residues on the capsid surface. This knowledge provides molecular-level insight into how to tailor the OPEs to optimize viral capsid disruption and increase OPE efficacy to target amphiphilic protein coats of icosahedral-based viruses.

A chimeric protein-based malaria vaccine candidate induces robust T cell responses against Plasmodium vivax MSP119

The most widespread Plasmodium species, Plasmodium vivax, poses a significant public health threat. An effective vaccine is needed to reduce global malaria burden. Of the erythrocytic stage vaccine candidates, the 19 kDa fragment of the P. vivax Merozoite Surface Protein 1 (PvMSP1₁₉) is one of the most promising. Our group has previously defined several promiscuous T helper epitopes within the PvMSP1 protein, with features that allow them to bind multiple MHC class II alleles. We describe here a P. vivax recombinant modular chimera based on MSP1 (PvRMC-MSP1) that includes defined T cell epitopes genetically fused to PvMSP1₁₉. This vaccine candidate preserved structural elements of the native PvMSP1₁₉ and elicited cytophilic antibody responses, and CD4+ and CD8+ T cells capable of recognizing PvMSP1₁₉. Although CD8+ T cells that recognize blood stage antigens have been reported to control blood infection, CD8+ T cell responses induced by P. falciparum or P. vivax vaccine candidates based on MSP1₁₉ have not been reported. To our knowledge, this is the first time a protein based subunit vaccine has been able to induce CD8+ T cell against PvMSP1₁₉. The PvRMC-MSP1 protein was also recognized by naturally acquired antibodies from individuals living in malaria endemic areas with an antibody profile associated with protection from infection. These features make PvRMC-MSP1 a promising vaccine candidate.

Functional Equivalence of Retroviral MA Domains in Facilitating Psi RNA Binding Specificity by Gag

Retroviruses specifically package full-length, dimeric genomic RNA (gRNA) even in the presence of a vast excess of cellular RNA. The “psi” (Ψ) element within the 5′-untranslated region (5′UTR) of gRNA is critical for packaging through interaction with the nucleocapsid (NC) domain of Gag. However, in vitro Gag binding affinity for Ψ versus non-Ψ RNAs is not significantly different. Previous salt-titration binding assays revealed that human immunodeficiency virus type 1 (HIV-1) Gag bound to Ψ RNA with high specificity and relatively few charge interactions, whereas binding to non-Ψ RNA was less specific and involved more electrostatic interactions. The NC domain was critical for specific Ψ binding, but surprisingly, a Gag mutant lacking the matrix (MA) domain was less effective at discriminating Ψ from non-Ψ RNA. We now find that Rous sarcoma virus (RSV) Gag also effectively discriminates RSV Ψ from non-Ψ RNA in a MA-dependent manner. Interestingly, Gag chimeras, wherein the HIV-1 and RSV MA domains were swapped, maintained high binding specificity to cognate Ψ RNAs. Using Ψ RNA mutant constructs, determinants responsible for promoting high Gag binding specificity were identified in both systems. Taken together, these studies reveal the functional equivalence of HIV-1 and RSV MA domains in facilitating Ψ RNA selectivity by Gag, as well as Ψ elements that promote this selectivity.

Molecular dynamics-based refinement and validation for sub-5 Å cryo-electron microscopy maps

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

DOI:http://dx.doi.org/10.7554/eLife.16105.001

To understand the roles that proteins and other large molecules play inside cells, it is important to determine their structures. One of the techniques that researchers can use to do this is called cryo-electron microscopy (cryo-EM), which rapidly freezes molecules to fix them in position before imaging them in fine detail.

The cryo-EM images are like maps that show the approximate position of atoms. These images must then be processed in order to build a three-dimensional model of the protein that shows how its atoms are arranged relative to each other. One computational approach called Molecular Dynamics Flexible Fitting (MDFF) works by flexibly fitting possible atomic structures into cryo-EM maps. Although this approach works well with relatively undetailed (or ‘low resolution’) cryo-EM images, it struggles to handle the high-resolution cryo-EM maps now being generated.

Singharoy, Teo, McGreevy et al. have now developed two MDFF methods – called cascade MDFF and resolution exchange MDFF – that help to resolve atomic models of biological molecules from cryo-EM images. Each method can refine poorly guessed models into ones that are consistent with the high-resolution experimental images. The refinement is achieved by interpreting a range of images that starts with a ‘fuzzy’ image. The contrast of the image is then progressively improved until an image is produced that has a resolution that is good enough to almost distinguish individual atoms.

The method works because each cryo-EM image shows not just one, but a collection of atomic structures that the molecule can take on, with the fuzzier parts of the image representing the more flexible parts of the molecule. By taking into account this flexibility, the large-scale features of the protein structure can be determined first from the fuzzier images, and increasing the contrast of the images allows smaller-scale refinements to be made to the structure.

The MDFF tools have been designed to be easy to use and are available to researchers at low cost through cloud computing platforms. They can now be used to unravel the structure of many different proteins and protein complexes including those involved in photosynthesis, respiration and protein synthesis.

DOI:http://dx.doi.org/10.7554/eLife.16105.002

Computational virology: From the inside out

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Contributions of Charged Residues in Structurally Dynamic Capsid Surface Loops to Rous Sarcoma Virus Assembly

Extensive studies of orthoretroviral capsids have shown that many regions of the CA protein play unique roles at different points in the virus life cycle. The N-terminal domain (NTD) flexible-loop (FL) region is one such example: exposed on the outer capsid surface, it has been implicated in Gag-mediated particle assembly, capsid maturation, and early replication events. We have now defined the contributions of charged residues in the FL region of the Rous sarcoma virus (RSV) CA to particle assembly. Effects of mutations on assembly were assessed in vivo and in vitro and analyzed in light of new RSV Gag lattice models. Virus replication was strongly dependent on the preservation of charge at a few critical positions in Gag-Gag interfaces. In particular, a cluster of charges at the beginning of FL contributes to an extensive electrostatic network that is important for robust Gag assembly and subsequent capsid maturation. Second-site suppressor analysis suggests that one of these charged residues, D87, has distal influence on interhexamer interactions involving helix α7. Overall, the tolerance of FL to most mutations is consistent with current models of Gag lattice structures. However, the results support the interpretation that virus evolution has achieved a charge distribution across the capsid surface that (i) permits the packing of NTD domains in the outer layer of the Gag shell, (ii) directs the maturational rearrangements of the NTDs that yield a functional core structure, and (iii) supports capsid function during the early stages of virus infection.

IMPORTANCE

The production of infectious retrovirus particles is a complex process, a choreography of protein and nucleic acid that occurs in two distinct stages: formation and release from the cell of an immature particle followed by an extracellular maturation phase during which the virion proteins and nucleic acids undergo major rearrangements that activate the infectious potential of the virion. This study examines the contributions of charged amino acids on the surface of the Rous sarcoma virus capsid protein in the assembly of appropriately formed immature particles and the maturational transitions that create a functional virion. The results provide important biological evidence in support of recent structural models of the RSV immature virions and further suggest that immature particle assembly and virion maturation are controlled by an extensive network of electrostatic interactions and long-range communication across the capsid surface.

All-Atom Molecular Dynamics of Virus Capsids as Drug Targets

Virus capsids are protein shells that package the viral genome. Although their morphology and biological functions can vary markedly, capsids often play critical roles in regulating viral infection pathways. A detailed knowledge of virus capsids, including their dynamic structure, interactions with cellular factors, and the specific roles that they play in the replication cycle, is imperative for the development of antiviral therapeutics. The following Perspective introduces an emerging area of computational biology that focuses on the dynamics of virus capsids and capsid-protein assemblies, with particular emphasis on the effects of small-molecule drug binding on capsid structure, stability, and allosteric pathways. When performed at chemical detail, molecular dynamics simulations can reveal subtle changes in virus capsids induced by drug molecules a fraction of their size. Here, the current challenges of performing all-atom capsid-drug simulations are discussed, along with an outlook on the applicability of virus capsid simulations to reveal novel drug targets.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Molecular Architecture of the Retroviral Capsid

Retroviral capsid cores are proteinaceous containers that self-assemble to encase the viral genome and a handful of proteins that promote infection. Their function is to protect and aid in the delivery of viral genes to the nucleus of the host, and, in many cases, infection pathways are influenced by capsid-cellular interactions. From a mathematical perspective, capsid cores are polyhedral cages and, as such, follow well-defined geometric rules. However, marked morphological differences in shapes exist, depending on virus type. Given the specific roles of capsid in the viral life cycle, the availability of detailed molecular structures, particularly at assembly interfaces, opens novel avenues for targeted drug development against these pathogens. Here, we summarize recent advances in the structure and understanding of retroviral capsid, with particular emphasis on assemblies and the capsid cores.

The Role of the Membrane in the Structure and Biophysical Robustness of the Dengue Virion Envelope

The dengue virion is surrounded by an envelope of membrane proteins surrounding a lipid bilayer. We have combined the cryoelectron microscopy structures of the membrane proteins (PDB: 3J27) with a lipid bilayer whose composition is based on lipidomics data for insect cell membranes, to obtain a near-atomic resolution computational model of the envelope of the dengue virion. A coarse-grained molecular dynamics simulation on the microsecond timescale enables analysis of key biophysical properties of the dengue outer envelope. Properties analyzed include area per lipid values (for a spherical virion with a mixed lipid composition), bilayer thickness, and lipid diffusion coefficients. Despite the absence of cholesterol from the lipid bilayer, the virion exhibits biophysical robustness (slow lipid diffusion alongside stable bilayer thickness, virion diameter, and shape) that matches the cholesterol-rich membrane of influenza A, with similarly anomalous diffusion of lipids. Biophysical robustness of the envelope may confer resilience to environmental perturbations.

•

The dengue virus envelope is a lipid bilayer plus an outer layer of membrane proteins

•

The structures of the proteins plus lipidomics data were used to model the envelope

•

Microsecond MD simulations revealed the dynamic behavior of the lipid bilayer

•

Protein interactions confer “raft-like” robustness on a cholesterol-free membrane

Membrane Binding of the Rous Sarcoma Virus Gag Protein Is Cooperative and Dependent on the Spacer Peptide Assembly Domain

The principles underlying membrane binding and assembly of retroviral Gag proteins into a lattice are understood. However, little is known about how these processes are related. Using purified Rous sarcoma virus Gag and Gag truncations, we studied the interrelation of Gag-Gag interaction and Gag-membrane interaction. Both by liposome binding and by surface plasmon resonance on a supported bilayer, Gag bound to membranes much more tightly than did matrix (MA), the isolated membrane binding domain. In principle, this difference could be explained either by protein-protein interactions leading to cooperativity in membrane binding or by the simultaneous interaction of the N-terminal MA and the C-terminal nucleocapsid (NC) of Gag with the bilayer, since both are highly basic. However, we found that NC was not required for strong membrane binding. Instead, the spacer peptide assembly domain (SPA), a putative 24-residue helical sequence comprising the 12-residue SP segment of Gag and overlapping the capsid (CA) C terminus and the NC N terminus, was required. SPA is known to be critical for proper assembly of the immature Gag lattice. A single amino acid mutation in SPA that abrogates assembly in vitro dramatically reduced binding of Gag to liposomes. In vivo, plasma membrane localization was dependent on SPA. Disulfide cross-linking based on ectopic Cys residues showed that the contacts between Gag proteins on the membrane are similar to the known contacts in virus-like particles. Taken together, we interpret these results to mean that Gag membrane interaction is cooperative in that it depends on the ability of Gag to multimerize.

IMPORTANCE

The retroviral structural protein Gag has three major domains. The N-terminal MA domain interacts directly with the plasma membrane (PM) of cells. The central CA domain, together with immediately adjoining sequences, facilitates the assembly of thousands of Gag molecules into a lattice. The C-terminal NC domain interacts with the genome, resulting in packaging of viral RNA. For assembly in vitro with purified Gag, in the absence of membranes, binding of NC to nucleic acid somehow facilitates further Gag-Gag interactions that lead to formation of the Gag lattice. The contributions of MA-mediated membrane binding to virus particle assembly are not well understood. Here, we report that in the absence of nucleic acid, membranes provide a platform that facilitates Gag-Gag interactions. This study demonstrates that the binding of Gag, but not of MA, to membranes is cooperative and identifies SPA as a major factor that controls this cooperativity.

The Structure of Immature Virus-Like Rous Sarcoma Virus Gag Particles Reveals a Structural Role for the p10 Domain in Assembly

The polyprotein Gag is the primary structural component of retroviruses. Gag consists of independently folded domains connected by flexible linkers. Interactions between the conserved capsid (CA) domains of Gag mediate formation of hexameric protein lattices that drive assembly of immature virus particles. Proteolytic cleavage of Gag by the viral protease (PR) is required for maturation of retroviruses from an immature form into an infectious form. Within the assembled Gag lattices of HIV-1 and Mason-Pfizer monkey virus (M-PMV), the C-terminal domain of CA adopts similar quaternary arrangements, while the N-terminal domain of CA is packed in very different manners. Here, we have used cryo-electron tomography and subtomogram averaging to study in vitro-assembled, immature virus-like Rous sarcoma virus (RSV) Gag particles and have determined the structure of CA and the surrounding regions to a resolution of ∼8 Å. We found that the C-terminal domain of RSV CA is arranged similarly to HIV-1 and M-PMV, whereas the N-terminal domain of CA adopts a novel arrangement in which the upstream p10 domain folds back into the CA lattice. In this position the cleavage site between CA and p10 appears to be inaccessible to PR. Below CA, an extended density is consistent with the presence of a six-helix bundle formed by the spacer-peptide region. We have also assessed the affect of lattice assembly on proteolytic processing by exogenous PR. The cleavage between p10 and CA is indeed inhibited in the assembled lattice, a finding consistent with structural regulation of proteolytic maturation.

IMPORTANCE

Retroviruses first assemble into immature virus particles, requiring interactions between Gag proteins that form a protein layer under the viral membrane. Subsequently, Gag is cleaved by the viral protease enzyme into separate domains, leading to rearrangement of the virus into its infectious form. It is important to understand how Gag is arranged within immature retroviruses, in order to understand how virus assembly occurs, and how maturation takes place. We used the techniques cryo-electron tomography and subtomogram averaging to obtain a detailed structural picture of the CA domains in immature assembled Rous sarcoma virus Gag particles. We found that part of Gag next to CA, called p10, folds back and interacts with CA when Gag assembles. This arrangement is different from that seen in HIV-1 and Mason-Pfizer monkey virus, illustrating further structural diversity of retroviral structures. The structure provides new information on how the virus assembles and undergoes maturation.