Mature HIV-1 capsid structure by cryo-electron microscopy and all-atom molecular dynamics

GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models

GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.

The nuclear localization signal of CPSF6 governs post-nuclear import steps of HIV-1 infection

The early stages of HIV-1 infection include the trafficking of the viral core into the nucleus of infected cells. However, much remains to be understood about how HIV-1 accomplishes nuclear import and the consequences of the import pathways utilized on nuclear events. The host factor cleavage and polyadenylation specificity factor 6 (CPSF6) assists HIV-1 nuclear localization and post-entry integration targeting. Here, we used a CPSF6 truncation mutant lacking a functional nuclear localization signal (NLS), CPSF6-358, and appended heterologous NLSs to rescue nuclear localization. We show that some, but not all, NLSs drive CPSF6-358 into the nucleus. Interestingly, we found that some nuclear localized CPSF6-NLS chimeras supported inefficient HIV-1 infection. We found that HIV-1 still enters the nucleus in these cell lines but fails to traffic to speckle-associated domains (SPADs). Additionally, we show that HIV-1 fails to efficiently integrate in these cell lines. Collectively, our results demonstrate that the NLS of CPSF6 facilitates steps of HIV-1 infection subsequent to nuclear import and additionally identify the ability of canonical NLS sequences to influence cargo localization in the nucleus following nuclear import.

Author Summary

During HIV-1 infection, the viral capsid, which encloses the viral genome and accessory proteins required for reverse transcription (RT) and integration, traffics towards the nucleus and enters through the nuclear pore complex (NPC). Following entry into the nucleus, RT is completed and viral capsid disassembles releasing the preintegration complex (PIC) to integrate with the host chromosome. In this study, we investigated the early HIV-1 host factor CPSF6, and specifically focused on the C-terminal short amino acid nuclear localization signal (NLS) in CPSF6, in mediating viral nuclear entry and subsequent gene expression. Altering the NLS in CPSF6 with NLS from other proteins, significantly impacted HIV-1's ability to infect those cells. We further showed this defect in infection occurred at the level of viral integration. This study highlights the importance of the NLS in CPSF6 in dictating the NPC it associates with and its effect on HIV-1 infection. Moreover, our study emphasizes the function of NLS in targeting host cargos to different nuclear entry pathways.

Comparative analysis of retroviral Gag-host cell interactions: focus on the nuclear interactome

Retroviruses exploit host proteins to assemble and release virions from infected cells. Previously, most studies focused on interacting partners of retroviral Gag proteins that localize to the cytoplasm or plasma membrane. Given that several full-length Gag proteins have been found in the nucleus, identifying the Gag-nuclear interactome has high potential for novel findings involving previously unknown host processes. Here we systematically compared nuclear factors identified in published HIV-1 proteomic studies and performed our own mass spectrometry analysis using affinity-tagged HIV-1 and RSV Gag proteins mixed with nuclear extracts. We identified 57 nuclear proteins in common between HIV-1 and RSV Gag, and a set of nuclear proteins present in our analysis and ≥ 1 of the published HIV-1 datasets. Many proteins were associated with nuclear processes which could have functional consequences for viral replication, including transcription initiation/elongation/termination, RNA processing, splicing, and chromatin remodeling. Examples include facilitating chromatin remodeling to expose the integrated provirus, promoting expression of viral genes, repressing the transcription of antagonistic cellular genes, preventing splicing of viral RNA, altering splicing of cellular RNAs, or influencing viral or host RNA folding or RNA nuclear export. Many proteins in our pulldowns common to RSV and HIV-1 Gag are critical for transcription, including PolR2B, the second largest subunit of RNA polymerase II (RNAPII), and LEO1, a PAF1C complex member that regulates transcriptional elongation, supporting the possibility that Gag influences the host transcription profile to aid the virus. Through the interaction of RSV and HIV-1 Gag with splicing-related proteins CBLL1, HNRNPH3, TRA2B, PTBP1 and U2AF1, we speculate that Gag could enhance unspliced viral RNA production for translation and packaging. To validate one putative hit, we demonstrated an interaction of RSV Gag with Mediator complex member Med26, required for RNA polymerase II-mediated transcription. Although 57 host proteins interacted with both Gag proteins, unique host proteins belonging to each interactome dataset were identified. These results provide a strong premise for future functional studies to investigate roles for these nuclear host factors that may have shared functions in the biology of both retroviruses, as well as functions specific to RSV and HIV-1, given their distinctive hosts and molecular pathology.

Cyclophilin A Facilitates HIV-1 DNA Integration

Cyclophilin A (CypA) promotes HIV-1 infection by facilitating reverse transcription, nuclear entry and by countering the antiviral activity of TRIM5α. These multifunctional roles of CypA are driven by its binding to the viral capsid. Interestingly, recent studies suggest that the HIV-1 capsid lattice enters the nucleus of an infected cell and uncoats just before integration. Therefore, we tested whether CypA-capsid interaction regulates post-nuclear entry steps of infection, particularly integration. First, we challenged CypA-expressing (CypA +/+ ) and CypA-depleted (CypA -/- ) cells with HIV-1 particles and quantified the resulting levels of provirus. Surprisingly, CypA-depletion significantly reduced integration, an effect that was independent of CypA's effect on reverse transcription, nuclear entry, and the presence or absence of TRIM5α. Additionally, cyclosporin A, an inhibitor that disrupts CypA-capsid binding, inhibited HIV-1 integration in CypA +/+ cells but not in CypA -/- cells. Accordingly, HIV-1 capsid mutants (G89V and P90A) deficient in CypA binding were also blocked at integration in CypA +/+ cells but not in CypA -/- cells. Then, to understand the mechanism, we assessed the integration activity of HIV-1 preintegration complexes (PICs) extracted from infected cells. The PICs from CypA -/- cells had lower activity in vitro compared to those from CypA +/+ cells. PICs from cells depleted for CypA and TRIM5α also had lower activity, suggesting that CypA's effect on PIC activity is independent of TRIM5α. Finally, addition of CypA protein significantly stimulated the integration activity of PICs extracted from both CypA +/+ and CypA -/- cells. Collectively, these results suggest that CypA promotes HIV-1 integration, a previously unknown role of this host factor.

Importance

HIV-1 capsid interaction with host cellular factors is essential for establishing a productive infection. However, the molecular details of such virus-host interactions are not fully understood. Cyclophilin A (CypA) is the first host protein identified to specifically bind to the HIV-1 capsid. Now it is established that CypA promotes reverse transcription and nuclear entry steps of HIV-1 infection. In this report, we show that CypA promotes HIV-1 integration by binding to the viral capsid. Specifically, our results demonstrate that CypA promotes HIV-1 integration by stimulating the activity of the viral preintegration complex and identifies a novel role of CypA during HIV-1 infection. This new knowledge is important because recent reports suggest that an operationally intact HIV-1 capsid enters the nucleus of an infected cell.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Three prime repair exonuclease 1 preferentially degrades the integration-incompetent HIV-1 DNA through favorable kinetics, thermodynamic, structural, and conformational properties

HIV-1 integration into the human genome is dependent on 3'-processing of the viral DNA. Recently, we reported that the cellular Three Prime Repair Exonuclease 1 (TREX1) enhances HIV-1 integration by degrading the unprocessed viral DNA, while the integration-competent 3'-processed DNA remained resistant. Here, we describe the mechanism by which the 3'-processed HIV-1 DNA resists TREX1-mediated degradation. Our kinetic studies revealed that the rate of cleavage (kcat) of the 3'-processed DNA was significantly lower (approximately 2-2.5-fold) than the unprocessed HIV-1 DNA by TREX1. The kcat values of human TREX1 for the processed U5 and U3 DNA substrates were 3.8 s-1 and 4.5 s-1, respectively. In contrast, the unprocessed U5 and U3 substrates were cleaved at 10.2 s-1 and 9.8 s-1, respectively. The efficiency of degradation (kcat/Km) of the 3'-processed DNA (U5-70.2 and U3-28.05 pM-1s-1) was also significantly lower than the unprocessed DNA (U5-103.1 and U3-65.3 pM-1s-1). Furthermore, the binding affinity (Kd) of TREX1 was markedly lower (∼2-fold) for the 3'-processed DNA than the unprocessed DNA. Molecular docking and dynamics studies revealed distinct conformational binding modes of TREX1 with the 3'-processed and unprocessed HIV-1 DNA. Particularly, the unprocessed DNA was favorably positioned in the active site with polar interactions with the catalytic residues of TREX1. Additionally, a stable complex was formed between TREX1 and the unprocessed DNA compared the 3'-processed DNA. These results pinpoint the mechanism by which TREX1 preferentially degrades the integration-incompetent HIV-1 DNA and reveal the unique structural and conformational properties of the integration-competent 3'-processed HIV-1 DNA.

WEDAP: A Python Package for Streamlined Plotting of Molecular Simulation Data

Given the growing interest in path sampling methods for extending the timescales of molecular dynamics (MD) simulations, there has been great interest in software tools that streamline the generation of plots for monitoring the progress of large-scale simulations. Here, we present the WEDAP Python package for simplifying the analysis of data generated from either conventional MD simulations or the weighted ensemble (WE) path sampling method, as implemented in the widely used WESTPA software package. WEDAP facilitates (i) the parsing of WE simulation data stored in highly compressed, hierarchical HDF5 files, and (ii) incorporates trajectory weights from WE simulations into all generated plots. Our Python package consists of multiple user-friendly interfaces: a command-line interface, a graphical user interface, and a Python application programming interface. We demonstrate the plotting features of WEDAP through a series of examples using data from WE and conventional MD simulations that focus on the HIV-1 capsid protein C-terminal domain dimer as a showcase system. The source code for WEDAP is freely available on GitHub at https://github.com/chonglab-pitt/wedap .

TomoNet: A streamlined cryogenic electron tomography software pipeline with automatic particle picking on flexible lattices

Cryogenic electron tomography (cryoET) is capable of determining in situ biological structures of molecular complexes at near-atomic resolution by averaging half a million subtomograms. While abundant complexes/particles are often clustered in arrays, precisely locating and seamlessly averaging such particles across many tomograms present major challenges. Here, we developed TomoNet, a software package with a modern graphical user interface to carry out the entire pipeline of cryoET and subtomogram averaging to achieve high resolution. TomoNet features built-in automatic particle picking and three-dimensional (3D) classification functions and integrates commonly used packages to streamline high-resolution subtomogram averaging for structures in 1D, 2D, or 3D arrays. Automatic particle picking is accomplished in two complementary ways: one based on template matching and the other using deep learning. TomoNet's hierarchical file organization and visual display facilitate efficient data management as required for large cryoET datasets. Applications of TomoNet to three types of datasets demonstrate its capability of efficient and accurate particle picking on flexible and imperfect lattices to obtain high-resolution 3D biological structures: virus-like particles, bacterial surface layers within cellular lamellae, and membranes decorated with nuclear egress protein complexes. These results demonstrate TomoNet's potential for broad applications to various cryoET projects targeting high-resolution in situ structures.

HIV-1 capsid stability and reverse transcription are finely balanced to minimize sensing of reverse transcription products via the cGAS-STING pathway

A critical determinant for early post-entry events, the HIV-1 capsid (CA) protein forms the conical core when it rearranges around the dimeric RNA genome and associated viral proteins. Although mutations in CA have been reported to alter innate immune sensing of HIV-1, a direct link between core stability and sensing of HIV-1 nucleic acids has not been established. Herein, we assessed how manipulating the stability of the CA lattice through chemical and genetic approaches affects innate immune recognition of HIV-1. We found that destabilization of the CA lattice resulted in potent sensing of reverse transcription products when destabilization per se does not completely block reverse transcription. Surprisingly, due to the combined effects of enhanced reverse transcription and defects in nuclear entry, two separate CA mutants that form hyperstable cores induced innate immune sensing more potently than destabilizing CA mutations. At low concentrations that allowed the accumulation of reverse transcription products, CA-targeting compounds GS-CA1 and lenacapavir measurably impacted CA lattice stability in cells and modestly enhanced innate immune sensing of HIV. Interestingly, innate immune activation observed with viruses containing unstable cores was abolished by low doses of lenacapavir. Innate immune activation observed with both hyperstable and unstable CA mutants was dependent on the cGAS-STING DNA-sensing pathway and reverse transcription. Overall, our findings demonstrate that CA lattice stability and reverse transcription are finely balanced to support reverse transcription and minimize cGAS-STING-mediated sensing of the resulting viral DNA.

IMPORTANCE

In HIV-1 particles, the dimeric RNA genome and associated viral proteins and enzymes are encased in a proteinaceous lattice composed of the viral capsid protein. Herein, we assessed how altering the stability of this capsid lattice through orthogonal genetic and chemical approaches impacts the induction of innate immune responses. Specifically, we found that decreasing capsid lattice stability results in more potent sensing of viral reverse transcription products, but not the genomic RNA, in a cGAS-STING-dependent manner. The recently developed capsid inhibitors lenacapavir and GS-CA1 enhanced the innate immune sensing of HIV-1. Unexpectedly, due to increased levels of reverse transcription and cytosolic accumulation of the resulting viral cDNA, capsid mutants with hyperstable cores also resulted in the potent induction of type I interferon-mediated innate immunity. Our findings suggest that HIV-1 capsid lattice stability and reverse transcription are finely balanced to minimize exposure of reverse transcription products in the cytosol of host cells.

Scalable computation of anisotropic vibrations for large macromolecular assemblies

The Normal Mode Analysis (NMA) is a standard approach to elucidate the anisotropic vibrations of macromolecules at their folded states, where low-frequency collective motions can reveal rearrangements of domains and changes in the exposed surface of macromolecules. Recent advances in structural biology have enabled the resolution of megascale macromolecules with millions of atoms. However, the calculation of their vibrational modes remains elusive due to the prohibitive cost associated with constructing and diagonalizing the underlying eigenproblem and the current approaches to NMA are not readily adaptable for efficient parallel computing on graphic processing unit (GPU). Here, we present eigenproblem construction and diagonalization approach that implements level-structure bandwidth-reducing algorithms to transform the sparse computation in NMA to a globally-sparse-yet-locally-dense computation, allowing batched tensor products to be most efficiently executed on GPU. We map, optimize, and compare several low-complexity Krylov-subspace eigensolvers, supplemented by techniques such as Chebyshev filtering, sum decomposition, external explicit deflation and shift-and-inverse, to allow fast GPU-resident calculations. The method allows accurate calculation of the first 1000 vibrational modes of some largest structures in PDB ( > 2.4 million atoms) at least 250 times faster than existing methods.

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

The capsid revolution

Lenacapavir, targeting the human immunodeficiency virus type-1 (HIV-1) capsid, is the first-in-class antiretroviral drug recently approved for clinical use. The development of Lenacapavir is attributed to the remarkable progress in our understanding of the capsid protein made during the last few years. Considered little more than a component of the virus shell to be shed early during infection, the capsid has been found to be a key player in the HIV-1 life cycle by interacting with multiple host factors, entering the nucleus, and directing integration. Here, we describe the key advances that led to this 'capsid revolution'.

Spatiotemporal binding of cyclophilin A and CPSF6 to capsid regulates HIV-1 nuclear entry and integration

Human immunodeficiency virus type 1 (HIV-1) capsid, which is the target of the antiviral lenacapavir, protects the viral genome and binds multiple host proteins to influence intracellular trafficking, nuclear import, and integration. Previously, we showed that capsid binding to cleavage and polyadenylation specificity factor 6 (CPSF6) in the cytoplasm is competitively inhibited by cyclophilin A (CypA) binding and regulates capsid trafficking, nuclear import, and infection. Here we determined that a capsid mutant with increased CypA binding affinity had significantly reduced nuclear entry and mislocalized integration. However, disruption of CypA binding to the mutant capsid restored nuclear entry, integration, and infection in a CPSF6-dependent manner. Furthermore, relocalization of CypA expression from the cell cytoplasm to the nucleus failed to restore mutant HIV-1 infection. Our results clarify that sequential binding of CypA and CPSF6 to HIV-1 capsid is required for optimal nuclear entry and integration targeting, informing antiretroviral therapies that contain lenacapavir.

Biomolecular dynamics with machine-learned quantum-mechanical force fields trained on diverse chemical fragments

The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.

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.

SPANA: Spatial decomposition analysis for cellular-scale molecular dynamics simulations

The rapid increase in computational power with the latest supercomputers has enabled atomistic molecular dynamics (MDs) simulations of biomolecules in biological membrane, cytoplasm, and other cellular environments. These environments often contain a million or more atoms to be simulated simultaneously. Therefore, their trajectory analyses involve heavy computations that can become a bottleneck in the computational studies. Spatial decomposition analysis (SPANA) is a set of analysis tools in the Generalized-Ensemble Simulation System (GENESIS) software package that can carry out MD trajectory analyses of large-scale biological simulations using multiple CPU cores in parallel. SPANA applies the spatial decomposition of a large biological system to distribute structural and dynamical analyses into individual CPU cores, which reduces the computational time and the memory size, significantly. SPANA opens new possibilities for detailed atomistic analyses of biomacromolecules as well as solvent water molecules, ions, and metabolites in MD simulation trajectories of very large biological systems containing more than millions of atoms in cellular environments.

Three Prime Repair Exonuclease 1 preferentially degrades the integration-incompetent HIV-1 DNA through favorable kinetics, thermodynamic, structural and conformational properties

HIV-1 integration into the human genome is dependent on 3'-processing of the reverse transcribed viral DNA. Recently, we reported that the cellular Three Prime Repair Exonuclease 1 (TREX1) enhances HIV-1 integration by degrading the unprocessed viral DNA, while the integration-competent 3'-processed DNA remained resistant. Here, we describe the mechanism by which the 3'-processed HIV-1 DNA resists TREX1-mediated degradation. Our kinetic studies revealed that the rate of cleavage (kcat) of the 3'-processed DNA was significantly lower than the unprocessed HIV-1 DNA by TREX1. The efficiency of degradation (kcat/KM) of the 3'-processed DNA was also significantly lower than the unprocessed DNA. Furthermore, the binding affinity (Kd) of TREX1 was markedly lower to the 3'-processed DNA compared to the unprocessed DNA. Molecular docking and dynamics studies revealed distinct conformational binding modes of TREX1 with the 3'-processed and unprocessed HIV-1 DNA. Particularly, the unprocessed DNA was favorably positioned in the active site with polar interactions with the catalytic residues of TREX1. Additionally, a stable complex was formed between TREX1 and the unprocessed DNA compared the 3'-processed DNA. These results pinpoint the biochemical mechanism by which TREX1 preferentially degrades the integration-incompetent HIV-1 DNA and reveal the unique structural and conformational properties of the integration-competent 3'-processed HIV-1 DNA.

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.

All-Atom Biomolecular Simulation in the Exascale Era

Exascale supercomputers have opened the door to dynamic simulations, facilitated by AI/ML techniques, that model biomolecular motions over unprecedented length and time scales. This new capability holds the potential to revolutionize our understanding of fundamental biological processes. Here we report on some of the major advances that were discussed at a recent CECAM workshop in Pisa, Italy, on the topic with a primary focus on atomic-level simulations. First, we highlight examples of current large-scale biomolecular simulations and the future possibilities enabled by crossing the exascale threshold. Next, we discuss challenges to be overcome in optimizing the usage of these powerful resources. Finally, we close by listing several grand challenge problems that could be investigated with this new computer architecture.

Integrative approaches for characterizing protein dynamics: NMR, CryoEM, and computer simulations

Proteins are inherently dynamic and their internal motions are essential for biological function. Protein motions cover a broad range of timescales: 10-14-10 s, spanning from sub-picosecond vibrational motions of atoms via microsecond loop conformational rearrangements to millisecond large amplitude domain reorientations. Observing protein dynamics over all timescales and connecting motions and structure to biological mechanisms requires integration of multiple experimental and computational techniques. This review reports on state-of-the-art approaches for assessing dynamics in biological systems using recent examples of virus assemblies, enzymes, and molecular machines. By integrating NMR spectroscopy in solution and the solid state, cryo electron microscopy, and molecular dynamics simulations, atomistic pictures of protein motions are obtained, not accessible from any single method in isolation. This information provides fundamental insights into protein behavior that can guide the development of future therapeutics.

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.

HIV-1-induced translocation of CPSF6 to biomolecular condensates

Cleavage and polyadenylation specificity factor subunit 6 (CPSF6, also known as CFIm68) is a 68 kDa component of the mammalian cleavage factor I (CFIm) complex that modulates mRNA alternative polyadenylation (APA) and determines 3' untranslated region (UTR) length, an important gene expression control mechanism. CPSF6 directly interacts with the HIV-1 core during infection, suggesting involvement in HIV-1 replication. Here, we review the contributions of CPSF6 to every stage of the HIV-1 replication cycle. Recently, several groups described the ability of HIV-1 infection to induce CPSF6 translocation to nuclear speckles, which are biomolecular condensates. We discuss the implications for CPSF6 localization in condensates and the potential role of condensate-localized CPSF6 in the ability of HIV-1 to control the protein expression pattern of the cell.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

HIV-1 capsid and viral DNA integration

IMPORTANCE

HIV-1 capsid protein (CA)-independently or by recruiting host factors-mediates several key steps of virus replication in the cytoplasm and nucleus of the target cell. Research in the recent years have established that CA is multifunctional and genetically fragile of all the HIV-1 proteins. Accordingly, CA has emerged as a validated and high priority therapeutic target, and the first CA-targeting antiviral drug was recently approved for treating multi-drug resistant HIV-1 infection. However, development of next generation CA inhibitors depends on a better understanding of CA's known roles, as well as probing of CA's novel roles, in HIV-1 replication. In this timely review, we present an updated overview of the current state of our understanding of CA's multifunctional role in HIV-1 replication-with a special emphasis on CA's newfound post-nuclear roles, highlight the pressing knowledge gaps, and discuss directions for future research.

Capsid-host interactions for HIV-1 ingress

The HIV-1 capsid, composed of approximately 1,200 copies of the capsid protein, encases genomic RNA alongside viral nucleocapsid, reverse transcriptase, and integrase proteins. After cell entry, the capsid interacts with a myriad of host factors to traverse the cell cytoplasm, pass through the nuclear pore complex (NPC), and then traffic to chromosomal sites for viral DNA integration. Integration may very well require the dissolution of the capsid, but where and when this uncoating event occurs remains hotly debated. Based on size constraints, a long-prevailing view was that uncoating preceded nuclear transport, but recent research has indicated that the capsid may remain largely intact during nuclear import, with perhaps some structural remodeling required for NPC traversal. Completion of reverse transcription in the nucleus may further aid capsid uncoating. One canonical type of host factor, typified by CPSF6, leverages a Phe-Gly (FG) motif to bind capsid. Recent research has shown these peptides reside amid prion-like domains (PrLDs), which are stretches of protein sequence devoid of charged residues. Intermolecular PrLD interactions along the exterior of the capsid shell impart avid host factor binding for productive HIV-1 infection. Herein we overview capsid-host interactions implicated in HIV-1 ingress and discuss important research questions moving forward. Highlighting clinical relevance, the long-acting ultrapotent inhibitor lenacapavir, which engages the same capsid binding pocket as FG host factors, was recently approved to treat people living with HIV.

Emerging role of cyclophilin A in HIV-1 infection: from producer cell to the target cell nucleus

The HIV-1 genome encodes a small number of proteins with structural, enzymatic, regulatory, and accessory functions. These viral proteins interact with a number of host factors to promote the early and late stages of HIV-1 infection. During the early stages of infection, interactions between the viral proteins and host factors enable HIV-1 to enter the target cell, traverse the cytosol, dock at the nuclear pore, gain access to the nucleus, and integrate into the host genome. Similarly, the viral proteins recruit another set of host factors during the late stages of infection to orchestrate HIV-1 transcription, translation, assembly, and release of progeny virions. Among the host factors implicated in HIV-1 infection, Cyclophilin A (CypA) was identified as the first host factor to be packaged within HIV-1 particles. It is now well established that CypA promotes HIV-1 infection by directly binding to the viral capsid. Mechanistic models to pinpoint CypA's role have spanned from an effect in the producer cell to the early steps of infection in the target cell. In this review, we will describe our understanding of the role(s) of CypA in HIV-1 infection, highlight the current knowledge gaps, and discuss the potential role of this host factor in the post-nuclear entry steps of HIV-1 infection.

Building the next generation of virtual cells to understand cellular biology

Cell science has made significant progress by focusing on understanding individual cellular processes through reductionist approaches. However, the sheer volume of knowledge collected presents challenges in integrating this information across different scales of space and time to comprehend cellular behaviors, as well as making the data and methods more accessible for the community to tackle complex biological questions. This perspective proposes the creation of next-generation virtual cells, which are dynamic 3D models that integrate information from diverse sources, including simulations, biophysical models, image-based models, and evidence-based knowledge graphs. These virtual cells would provide statistically accurate and holistic views of real cells, bridging the gap between theoretical concepts and experimental data, and facilitating productive new collaborations among researchers across related fields.

Multidisciplinary studies with mutated HIV-1 capsid proteins reveal structural mechanisms of lattice stabilization

HIV-1 capsid (CA) stability is important for viral replication. E45A and P38A mutations enhance and reduce core stability, thus impairing infectivity. Second-site mutations R132T and T216I rescue infectivity. Capsid lattice stability was studied by solving seven crystal structures (in native background), including P38A, P38A/T216I, E45A, E45A/R132T CA, using molecular dynamics simulations of lattices, cryo-electron microscopy of assemblies, time-resolved imaging of uncoating, biophysical and biochemical characterization of assembly and stability. We report pronounced and subtle, short- and long-range rearrangements: (1) A38 destabilized hexamers by loosening interactions between flanking CA protomers in P38A but not P38A/T216I structures. (2) Two E45A structures showed unexpected stabilizing CANTD-CANTD inter-hexamer interactions, variable R18-ring pore sizes, and flipped N-terminal β-hairpin. (3) Altered conformations of E45Aa α9-helices compared to WT, E45A/R132T, WTPF74, WTNup153, and WTCPSF6 decreased PF74, CPSF6, and Nup153 binding, and was reversed in E45A/R132T. (4) An environmentally sensitive electrostatic repulsion between E45 and D51 affected lattice stability, flexibility, ion and water permeabilities, electrostatics, and recognition of host factors.

Lipid nanodiscs as a template for high-resolution cryo-EM structures of peripheral membrane proteins

Peripheral membrane proteins are ubiquitous throughout cell biology and are required for a variety of cellular processes such as signal transduction, membrane trafficking, and autophagy. Transient binding to the membrane has a profound impact on protein function, serving to induce conformational changes and alter biochemical and biophysical parameters by increasing the local concentration of factors and restricting diffusion to two dimensions. Despite the centrality of the membrane in serving as a template for cell biology, there are few reported high-resolution structures of peripheral membrane proteins bound to the membrane. We analyzed the utility of lipid nanodiscs to serve as a template for cryo-EM analysis of peripheral membrane proteins. We tested a variety of nanodiscs and we report a 3.3 Å structure of the AP2 clathrin adaptor complex bound to a 17-nm nanodisc, with sufficient resolution to visualize a bound lipid head group. Our data demonstrate that lipid nanodiscs are amenable to high-resolution structure determination of peripheral membrane proteins and provide a framework for extending this analysis to other systems.

Cryo-electron tomography to study viral infection

Developments in cryo-electron microscopy (cryo-EM) have been interwoven with the study of viruses ever since its first applications to biological systems. Following the success of single particle cryo-EM in the last decade, cryo-electron tomography (cryo-ET) is now rapidly maturing as a technology and catalysing great advancement in structural virology as its application broadens. In this review, we provide an overview of the use of cryo-ET to study viral infection biology, discussing the key workflows and strategies used in the field. We highlight the vast body of studies performed on purified viruses and virus-like particles (VLPs), as well as discussing how cryo-ET can characterise host-virus interactions and membrane fusion events. We further discuss the importance of in situ cellular imaging in revealing previously unattainable details of infection and highlight the need for validation of high-resolution findings from purified ex situ systems. We give perspectives for future developments to achieve the full potential of cryo-ET to characterise the molecular processes of viral infection.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Sub-stoichiometric Modulation of Viral Targets-Potent Antiviral Agents That Exploit Target Vulnerability

The modulation of oligomeric viral targets at sub-stoichiometric ratios of drug to target has been advocated for its efficacy and potency, but there are only a limited number of documented examples. In this Viewpoint, we summarize the invention of the HIV-1 maturation inhibitor fipravirimat and discuss the emerging details around the mode of action of this class of drug that reflects inhibition of a protein composed of 1,300-1,600 monomers that interact in a cooperative fashion. Similarly, the HCV NS5A inhibitor daclatasvir has been shown to act in a highly sub-stoichiometric fashion, inhibiting viral replication at concentrations that are ∼23,500 lower than that of the protein target.

HIV-2 Immature Particle Morphology Provides Insights into Gag Lattice Stability and Virus Maturation

Retrovirus immature particle morphology consists of a membrane enclosed, pleomorphic, spherical and incomplete lattice of Gag hexamers. Previously, we demonstrated that human immunodeficiency virus type 2 (HIV-2) immature particles possess a distinct and extensive Gag lattice morphology. To better understand the nature of the continuously curved hexagonal Gag lattice, we have used the single particle cryo-electron microscopy method to determine the HIV-2 Gag lattice structure for immature virions. The reconstruction map at 5.5 Å resolution revealed a stable, wineglass-shaped Gag hexamer structure with structural features consistent with other lentiviral immature Gag lattice structures. Cryo-electron tomography provided evidence for nearly complete ordered Gag lattice structures in HIV-2 immature particles. We also solved a 1.98 Å resolution crystal structure of the carboxyl-terminal domain (CTD) of the HIV-2 capsid (CA) protein that identified a structured helix 12 supported via an interaction of helix 10 in the absence of the SP1 region of Gag. Residues at the helix 10-12 interface proved critical in maintaining HIV-2 particle release and infectivity. Taken together, our findings provide the first 3D organization of HIV-2 immature Gag lattice and important insights into both HIV Gag lattice stabilization and virus maturation.

How hydrophobicity shapes the architecture of protein assemblies

The interactions that give rise to protein self-assembly are basically electrical and hydrophobic in origin. The electrical interactions are approached in this study as the interaction between electrostatic dipoles originated by the asymmetric distribution of their charged amino acids. However, hydrophobicity is not easily derivable from basic physicochemical principles. Its treatment is carried out here considering a hydrophobic force field originated by "hydrophobic charges". These charges are indices obtained experimentally from the free energies of transferring amino acids from polar to hydrophobic media. Hydrophobic dipole moments are used here in a manner analogous to electric dipole moments, and an empirical expression of interaction energy between hydrophobic dipoles is derived. This methodology is used with two examples of self-assembly systems of different complexity. It was found that the hydrophobic dipole moments of proteins tend to interact in such a way that they align parallel to each other in a completely analogous way to how phospholipids are oriented in biological membranes to form the well-known double layer. In this biological membrane model (BM model), proteins tend to interact in a similar way, although in this case this alignment is modulated by the tendency of the corresponding electrostatic dipoles to counter-align. Helical conformation of influenza virus PDBid: 6Z5L. Two monomers are shown in cyan and green. The corresponding dipole moment vectors are shown in red (electric dipoles) and blue (hydrophobic dipoles). From the inset figure, it can be seen that the growth of the helix is due to electrical attraction of the monomers, overcoming a hydrophobic repulsion (see text).

Unconventional stabilization of the human T-cell leukemia virus type 1 immature Gag lattice

Human T-cell leukemia virus type 1 (HTLV-1) has an atypical immature particle morphology compared to other retroviruses. This indicates that these particles are formed in a way that is unique. Here we report the results of cryo-electron tomography (cryo-ET) studies of HTLV-1 virus-like particles (VLPs) assembled in vitro, as well as derived from cells. This work shows that HTLV-1 employs an unconventional mechanism of Gag-Gag interactions to form the immature viral lattice. Analysis of high-resolution structural information from immature CA tubular arrays reveals that the primary stabilizing component in HTLV-1 is CA-NTD. Mutagenesis and biophysical analysis support this observation. This distinguishes HTLV-1 from other retroviruses, in which the stabilization is provided primarily by the CA-CTD. These results are the first to provide structural details of the quaternary arrangement of Gag for an immature deltaretrovirus, and this helps explain why HTLV-1 particles are morphologically distinct.

Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories

The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.

HIV-1 capsid stability enables inositol phosphate-independent infection of target cells and promotes integration into genes

The mature HIV-1 capsid is stabilized by host and viral determinants. The capsid protein CA binds to the cellular metabolites inositol hexakisphosphate (IP6) and its precursor inositol (1, 3, 4, 5, 6) pentakisphosphate (IP5) to stabilize the mature capsid. In target cells, capsid destabilization by the antiviral compounds lenacapavir and PF74 reveals a HIV-1 infectivity defect due to IP5/IP6 (IP5/6) depletion. To test whether intrinsic HIV-1 capsid stability and/or host factor binding determines HIV-1 insensitivity to IP5/6 depletion, a panel of CA mutants was assayed for infection of IP5/6-depleted T cells and wildtype cells. Four CA mutants with unstable capsids exhibited dependence on host IP5/6 for infection and reverse transcription (RTN). Adaptation of one such mutant, Q219A, by spread in culture resulted in Vpu truncation and a capsid three-fold interface mutation, T200I. T200I increased intrinsic capsid stability as determined by in vitro uncoating of purified cores and partially reversed the IP5/6-dependence in target cells for each of the four CA mutants. T200I further rescued the changes to lenacapavir sensitivity associated with the parental mutation. The premature dissolution of the capsid caused by the IP5/6-dependent mutations imparted a unique defect in integration targeting that was rescued by T200I. Collectively, these results demonstrate that T200I restored other capsid functions after RTN for the panel of mutants. Thus, the hyperstable T200I mutation stabilized the instability defects imparted by the parental IP5/6-dependent CA mutation. The contribution of Vpu truncation to mutant adaptation was linked to BST-2 antagonization, suggesting that cell-to-cell transfer promoted replication of the mutants. We conclude that interactions at the three-fold interface are adaptable, key mediators of capsid stability in target cells and are able to antagonize even severe capsid instability to promote infection.

Heterologous Assembly of Pleomorphic Bacterial Microcompartment Shell Architectures Spanning the Nano- to Microscale

Many bacteria use protein-based organelles known as bacterial microcompartments (BMCs) to organize and sequester sequential enzymatic reactions. Regardless of their specialized metabolic function, all BMCs are delimited by a shell made of multiple structurally redundant, yet functionally diverse, hexameric (BMC-H), pseudohexameric/trimeric (BMC-T), or pentameric (BMC-P) shell protein paralogs. When expressed without their native cargo, shell proteins have been shown to self-assemble into 2D sheets, open-ended nanotubes, and closed shells of ≈40 nm diameter that are being developed as scaffolds and nanocontainers for applications in biotechnology. Here, by leveraging a strategy for affinity-based purification, it is demonstrated that a wide range of empty synthetic shells, many differing in end-cap structures, can be derived from a glycyl radical enzyme-associated microcompartment. The range of pleomorphic shells observed, which span ≈2 orders of magnitude in size from ≈25 nm to ≈1.8 µm, reveal the remarkable plasticity of BMC-based biomaterials. In addition, new capped nanotube and nanocone morphologies are observed that are consistent with a multicomponent geometric model in which architectural principles are shared among asymmetric carbon, viral protein, and BMC-based structures.

The Effect of Inositol Hexakisphosphate on HIV-1 Particle Production and Infectivity can be Modulated by Mutations that Affect the Stability of the Immature Gag Lattice

The assembly of an HIV-1 particle begins with the construction of a spherical lattice composed of hexamer subunits of the Gag polyprotein. The cellular metabolite inositol hexakisphosphate (IP6) binds and stabilizes the immature Gag lattice via an interaction with the six-helix bundle (6HB), a crucial structural feature of Gag hexamers that modulates both virus assembly and infectivity. The 6HB must be stable enough to promote immature Gag lattice formation, but also flexible enough to be accessible to the viral protease, which cleaves the 6HB during particle maturation. 6HB cleavage liberates the capsid (CA) domain of Gag from the adjacent spacer peptide 1 (SP1) and IP6 from its binding site. This pool of IP6 molecules then promotes the assembly of CA into the mature conical capsid that is required for infection. Depletion of IP6 in virus-producer cells results in severe defects in assembly and infectivity of wild-type (WT) virions. Here we show that in an SP1 double mutant (M4L/T8I) with a hyperstable 6HB, IP6 can block virion infectivity by preventing CA-SP1 processing. Thus, depletion of IP6 in virus-producer cells markedly increases M4L/T8I CA-SP1 processing and infectivity. We also show that the introduction of the M4L/T8I mutations partially rescues the assembly and infectivity defects induced by IP6 depletion on WT virions, likely by increasing the affinity of the immature lattice for limiting IP6. These findings reinforce the importance of the 6HB in virus assembly, maturation, and infection and highlight the ability of IP6 to modulate 6HB stability.

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.

Structural insights into HIV-1 polyanion-dependent capsid lattice formation revealed by single particle cryo-EM

The HIV-1 capsid houses the viral genome and interacts extensively with host cell proteins throughout the viral life cycle. It is composed of capsid protein (CA), which assembles into a conical fullerene lattice composed of roughly 200 CA hexamers and 12 CA pentamers. Previous structural analyses of individual CA hexamers and pentamers have provided valuable insight into capsid structure and function, but detailed structural information about these assemblies in the broader context of the capsid lattice is lacking. In this study, we combined cryoelectron tomography and single particle analysis (SPA) cryoelectron microscopy to determine structures of continuous regions of the capsid lattice containing both hexamers and pentamers. We also developed a method of liposome scaffold-based in vitro lattice assembly ("lattice templating") that enabled us to directly study the lattice under a wider range of conditions than has previously been possible. Using this approach, we identified a critical role for inositol hexakisphosphate in pentamer formation and determined the structure of the CA lattice bound to the capsid-targeting antiretroviral drug GS-6207 (lenacapavir). Our work reveals key structural details of the mature HIV-1 CA lattice and establishes the combination of lattice templating and SPA as a robust strategy for studying retroviral capsid structure and capsid interactions with host proteins and antiviral compounds.

Electrostatic Screening, Acidic pH and Macromolecular Crowding Increase the Self-Assembly Efficiency of the Minute Virus of Mice Capsid In Vitro

The hollow protein capsids from a number of different viruses are being considered for multiple biomedical or nanotechnological applications. In order to improve the applied potential of a given viral capsid as a nanocarrier or nanocontainer, specific conditions must be found for achieving its faithful and efficient assembly in vitro. The small size, adequate physical properties and specialized biological functions of the capsids of parvoviruses such as the minute virus of mice (MVM) make them excellent choices as nanocarriers and nanocontainers. In this study we analyzed the effects of protein concentration, macromolecular crowding, temperature, pH, ionic strength, or a combination of some of those variables on the fidelity and efficiency of self-assembly of the MVM capsid in vitro. The results revealed that the in vitro reassembly of the MVM capsid is an efficient and faithful process. Under some conditions, up to ~40% of the starting virus capsids were reassembled in vitro as free, non aggregated, correctly assembled particles. These results open up the possibility of encapsidating different compounds in VP2-only capsids of MVM during its reassembly in vitro, and encourage the use of virus-like particles of MVM as nanocontainers.

HIV-1 Group M Capsid Amino Acid Variability: Implications for Sequence Quality Control of Genotypic Resistance Testing

BACKGROUND

With the approval of the HIV-1 capsid inhibitor, lenacapavir, capsid sequencing will be required for managing lenacapavir-experienced individuals with detectable viremia. Successful sequence interpretation will require examining new capsid sequences in the context of previously published sequence data.

METHODS

We analyzed published HIV-1 group M capsid sequences from 21,012 capsid-inhibitor naïve individuals to characterize amino acid variability at each position and influence of subtype and cytotoxic T lymphocyte (CTL) selection pressure. We determined the distributions of usual mutations, defined as amino acid differences from the group M consensus, with a prevalence ≥ 0.1%. Co-evolving mutations were identified using a phylogenetically-informed Bayesian graphical model method.

RESULTS

162 (70.1%) positions had no usual mutations (45.9%) or only conservative usual mutations with a positive BLOSUM62 score (24.2%). Variability correlated independently with subtype-specific amino acid occurrence (Spearman rho = 0.83; p < 1 × 10-9) and the number of times positions were reported to contain an HLA-associated polymorphism, an indicator of CTL pressure (rho = 0.43; p = 0.0002).

CONCLUSIONS

Knowing the distribution of usual capsid mutations is essential for sequence quality control. Comparing capsid sequences from lenacapavir-treated and lenacapavir-naïve individuals will enable the identification of additional mutations potentially associated with lenacapavir therapy.

Engineering excitonically coupled dimers in an artificial protein for light harvesting via computational modeling

In photosynthesis, pigment-protein complexes achieve outstanding photoinduced charge separation efficiencies through a set of strategies in which excited states delocalization over multiple pigments ("excitons") and charge-transfer states play key roles. These concepts, and their implementation in bioinspired artificial systems, are attracting increasing attention due to the vast potential that could be tapped by realizing efficient photochemical reactions. In particular, de novo designed proteins provide a diverse structural toolbox that can be used to manipulate the geometric and electronic properties of bound chromophore molecules. However, achieving excitonic and charge-transfer states requires closely spaced chromophores, a non-trivial aspect since a strong binding with the protein matrix needs to be maintained. Here, we show how a general-purpose artificial protein can be optimized via molecular dynamics simulations to improve its binding capacity of a chlorophyll derivative, achieving complexes in which chromophores form two closely spaced and strongly interacting dimers. Based on spectroscopy results and computational modeling, we demonstrate each dimer is excitonically coupled, and propose they display signatures of charge-transfer state mixing. This work could open new avenues for the rational design of chromophore-protein complexes with advanced functionalities.

Molecular size dependence on achievable resolution from XFEL single-particle 3D reconstruction

Single-particle analysis using x-ray free-electron lasers (XFELs) is a novel method for obtaining structural information of samples in a state close to nature. In particular, it is suitable for observing the inner structure of large biomolecules by taking advantage of the high transmittance of x-rays. However, systematic studies on the resolution achievable for large molecules are lacking. In this study, the molecular size dependence of the resolution of a three-dimensional (3D) structure resulting from XFEL single-particle reconstruction is evaluated using synthetic data. Evidently, 3D structures of larger molecules can be restored with higher detail (defined relative to the molecular sizes) than smaller ones; however, reconstruction with high absolute resolution (defined in nm^-1) is challenging. Our results provide useful information for the experimental design of 3D structure reconstruction using coherent x-ray diffraction patterns of single-particles.

Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike

Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.

Understanding Retroviral Life Cycle and its Genomic RNA Packaging

Members of the family Retroviridae are important animal and human pathogens. Being obligate parasites, their replication involves a series of steps during which the virus hijacks the cellular machinery. Additionally, many of the steps of retrovirus replication are unique among viruses, including reverse transcription, integration, and specific packaging of their genomic RNA (gRNA) as a dimer. Progress in retrovirology has helped identify several molecular mechanisms involved in each of these steps, but many are still unknown or remain controversial. This review summarizes our present understanding of the molecular mechanisms involved in various stages of retrovirus replication. Furthermore, it provides a comprehensive analysis of our current understanding of how different retroviruses package their gRNA into the assembling virions. RNA packaging in retroviruses holds a special interest because of the uniqueness of packaging a dimeric genome. Dimerization and packaging are highly regulated and interlinked events, critical for the virus to decide whether its unspliced RNA will be packaged as a "genome" or translated into proteins. Finally, some of the outstanding areas of exploration in the field of RNA packaging are highlighted, such as the role of epitranscriptomics, heterogeneity of transcript start sites, and the necessity of functional polyA sequences. An in-depth knowledge of mechanisms that interplay between viral and cellular factors during virus replication is critical in understanding not only the virus life cycle, but also its pathogenesis, and development of new antiretroviral compounds, vaccines, as well as retroviral-based vectors for human gene therapy.

Investigational drugs for HIV: trends, opportunities and key players

INTRODUCTION

Since the first antiretroviral drug was described, the field of HIV treatment and prevention has undergone two drug-based revolutions: the first one, enabled by the virtually concomitant discovery of non-nucleoside reverse transcriptase and protease inhibitors, was the inception of combined antiretroviral therapy. The second followed the creation of integrase strand-transfer inhibitors with improved safety, potency, and resistance profiles. Long-acting antiretroviral drugs, including broadly neutralizing antibodies, now offer the opportunity for a third transformational change in HIV management.

AREAS COVERED

Our review focused on HIV treatment and prevention with investigational drugs that offer the potential for infrequent dosing, including drugs not yet approved for clinical use. We also discussed approved drugs for which administration modalities or formulations are being optimized. We performed a literature search in published manuscripts, conference communications, and registered clinical trials.

EXPERT OPINION

While the field focuses on extending dosing intervals, we identify drug tissue penetration as an understudied opportunity to improve HIV care. We repeat that self-administration remains an essential milestone to reach the full potential of long-acting drugs. Treatments and prevention strategies based on broadly neutralizing antibodies require a deeper understanding of their antiretroviral properties.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Computing the Volume, Surface Area, Mean, and Gaussian Curvatures of Molecules and Their Derivatives

Geometry is crucial in our efforts to comprehend the structures and dynamics of biomolecules. For example, volume, surface area, and integrated mean and Gaussian curvature of the union of balls representing a molecule are used to quantify its interactions with the water surrounding it in the morphometric implicit solvent models. The Alpha Shape theory provides an accurate and reliable method for computing these geometric measures. In this paper, we derive homogeneous formulas for the expressions of these measures and their derivatives with respect to the atomic coordinates, and we provide algorithms that implement them into a new software package, AlphaMol. The only variables in these formulas are the interatomic distances, making them insensitive to translations and rotations. AlphaMol includes a sequential algorithm and a parallel algorithm. In the parallel version, we partition the atoms of the molecule of interest into 3D rectangular blocks, using a kd-tree algorithm. We then apply the sequential algorithm of AlphaMol to each block, augmented by a buffer zone to account for atoms whose ball representations may partially cover the block. The current parallel version of AlphaMol leads to a 20-fold speed-up compared to an independent serial implementation when using 32 processors. For instance, it takes 31 s to compute the geometric measures and derivatives of each atom in a viral capsid with more than 26 million atoms on 32 Intel processors running at 2.7 GHz. The presence of the buffer zones, however, leads to redundant computations, which ultimately limit the impact of using multiple processors. AlphaMol is available as an OpenSource software.