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

Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation

Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π^∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.

A computational study of the interface interaction between SARS-CoV-2 RBD and ACE2 from human, cat, dog, and ferret

The total impact of the worldwide COVID‐19 pandemic is still emerging, changing all relationships as a result, including those with pet animals. In the infection process, the use of angiotensin‐converting enzyme 2 (ACE2) as a cellular receptor to the spike protein of the new coronavirus is a fundamental step. In this sense, understanding which residue plays what role in the interaction between SARS‐CoV‐2 spike glycoprotein and ACE2 from cats, dogs, and ferrets is an important guide for helping to choose which animal model can be used to study the pathology of COVID‐19, and if there are differences between these interactions and those occurring in the human system. To help answer these questions, we performed classical molecular dynamics simulations to evaluate, from an atomistic point of view, the interactions in these systems. Our results show that there are significant differences in the interacting residues between the systems from different animal species, and the role of ACE2 key residues are different in each system, and can assist in the search for different inhibitors for each animal.

The Role of Capsid in the Early Steps of HIV-1 Infection: New Insights into the Core of the Matter

In recent years, major advances in research and experimental approaches have significantly increased our knowledge on the role of the HIV-1 capsid in the virus life cycle, from reverse transcription to integration and gene expression. This makes the capsid protein a good pharmacological target to inhibit HIV-1 replication. This review covers our current understanding of the role of the viral capsid in the HIV-1 life cycle and its interaction with different host factors that enable reverse transcription, trafficking towards the nucleus, nuclear import and integration into host chromosomes. It also describes different promising small molecules, some of them in clinical trials, as potential targets for HIV-1 therapy.

Algorithm to minimize MPI communications in the parallelized fast multipole method combined with molecular dynamics calculations

In the era of exascale supercomputers, large-scale, and long-time molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)-communication in the MPI-parallelized fast multipole method (FMM) combined with MD calculations under three-dimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimum-transferred data (MTD) method, should enable large-scale and long-time MD calculations to be calculated efficiently, under the condition of massive MPI-parallelization on exascale supercomputers.

NO● Represses the Oxygenation of Arachidonoyl PE by 15LOX/PEBP1: Mechanism and Role in Ferroptosis

We recently discovered an anti-ferroptotic mechanism inherent to M1 macrophages whereby high levels of NO^● suppressed ferroptosis via inhibition of hydroperoxy-eicosatetraenoyl-phosphatidylethanolamine (HpETE-PE) production by 15-lipoxygenase (15LOX) complexed with PE-binding protein 1 (PEBP1). However, the mechanism of NO^● interference with 15LOX/PEBP1 activity remained unclear. Here, we use a biochemical model of recombinant 15LOX-2 complexed with PEBP1, LC-MS redox lipidomics, and structure-based modeling and simulations to uncover the mechanism through which NO^● suppresses ETE-PE oxidation. Our study reveals that O₂ and NO^● use the same entry pores and channels connecting to 15LOX-2 catalytic site, resulting in a competition for the catalytic site. We identified residues that direct O₂ and NO^● to the catalytic site, as well as those stabilizing the esterified ETE-PE phospholipid tail. The functional significance of these residues is supported by in silico saturation mutagenesis. We detected nitrosylated PE species in a biochemical system consisting of 15LOX-2/PEBP1 and NO^● donor and in RAW264.7 M2 macrophages treated with ferroptosis-inducer RSL3 in the presence of NO^●, in further support of the ability of NO^● to diffuse to, and react at, the 15LOX-2 catalytic site. The results provide first insights into the molecular mechanism of repression of the ferroptotic Hp-ETE-PE production by NO^●.

Molecular determinants of Ebola nucleocapsid stability from molecular dynamics simulations

Ebola virus (EBOV) is a human pathogen with the ability to cause hemorrhagic fever and bleeding diathesis in hosts. The life cycle of EBOV depends on its nucleocapsid. The Ebola nucleocapsid consists of a helical assembly of nucleoproteins (NPs) encapsidating single-stranded viral RNA (ssRNA). Knowledge of the molecular determinants of Ebola nucleocapsid stability is essential for the development of therapeutics against EBOV. However, large degrees of freedom associated with the Ebola nucleocapsid helical assembly pose a computational challenge, thereby limiting the previous simulation studies to the level of monomers. In the present work, we have performed all atom molecular dynamics (MD) simulations of the helical assembly of EBOV nucleoproteins in the absence and presence of ssRNA. We found that ssRNA is essential for maintaining structural integrity of the nucleocapsid. Other molecular determinants observed to stabilize the nucleocapsid include NP-RNA and NP-NP interactions and ion distributions. Additionally, the structural and dynamical behavior of the nucleocapsid monomer depends on its position in the helical assembly. NP monomers present on the longitudinal edges of the helical tube are more exposed, flexible, and have weaker NP-NP interactions than those residing in the center. This work provides key structural features stabilizing the nucleocapsid that may serve as therapeutic targets.

Fighting viruses with computers, right now

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

Integrative structural biology of HIV-1 capsid protein assemblies: combining experiment and computation

HIV-1 is the causative agent of acquired immunodeficiency syndrome (AIDS), a global pandemic that has claimed 32.7 million lives since 1981. Despite decades of research, there is no cure for the disease, with 38 million people currently infected with HIV. Attractive therapeutic targets for drug development are mature HIV-1 capsids, immature Gag polyprotein assemblies, and Gag maturation intermediates, although their complex architectures, pleomorphism, and dynamics render these assemblies challenging for structural biology. The recent development of integrative approaches, combining experimental and computational methods has enabled atomic-level characterization of structures and dynamics of capsid and Gag assemblies, and revealed their interactions with small-molecule inhibitors and host factors. These structures provide important insights that will guide the development of capsid and maturation inhibitors.

Atomistic dynamics of a viral infection process: Release of membrane lytic peptides from a non-enveloped virus

Molecular simulations have played an instrumental role in uncovering the structural dynamics and physical properties of virus capsids. In this work, we move beyond equilibrium physicochemical characterization of a virus system to study a stage of the infection process that is required for viral proliferation. Despite many biochemical and functional studies, the molecular mechanism of host cell entry by non-enveloped viruses remains largely unresolved. Flock House virus (FHV) is a model system for non-enveloped viruses and is the subject of the current study. FHV infects through the acid-dependent endocytic pathway, where low pH triggers externalization of membrane-disrupting (γ) peptides from the capsid interior. Using all-atom equilibrium and enhanced sampling simulations, the mechanism and energetics of γ peptide liberation and the effect of pH on this process are investigated. Our computations agree with experimental findings and reveal nanoscopic details regarding the pH control mechanism, which are not readily accessible in experiments.

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.

Vpx enhances innate immune responses independently of SAMHD1 during HIV-1 infection

Background

The genomes of HIV-2 and some SIV strains contain the accessory gene vpx, which carries out several functions during infection, including the downregulation of SAMHD1. Vpx is also commonly used in experiments to increase HIV-1 infection efficiency in myeloid cells, particularly in studies that investigate the activation of antiviral pathways. However, the potential effects of Vpx on cellular innate immune signaling is not completely understood. We investigated whether and how Vpx affects ISG responses in monocytic cell lines and MDMs during HIV-1 infection.

Results

HIV-1 infection at excessively high virus doses can induce ISG activation, although at the expense of high levels of cell death. At equal infection levels, the ISG response is potentiated by the presence of Vpx and requires the initiation of reverse transcription. The interaction of Vpx with the DCAF1 adaptor protein is important for the enhanced response, implicating Vpx-mediated degradation of a host factor. Cells lacking SAMHD1 show similarly augmented responses, suggesting an effect that is independent of SAMHD1 degradation. Overcoming SAMHD1 restriction in MDMs to reach equal infection levels with viruses containing and lacking Vpx reveals a novel function of Vpx in elevating innate immune responses.

Conclusions

Vpx likely has as yet undefined roles in infected cells. Our results demonstrate that Vpx enhances ISG responses in myeloid cell lines and primary cells independently of its ability to degrade SAMHD1. These findings have implications for innate immunity studies in myeloid cells that use Vpx delivery with HIV-1 infection.

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

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

Wrapping Up Viruses at Multiscale Resolution: Optimizing PACKMOL and SIRAH Execution for Simulating the Zika Virus

Simulating huge biomolecular complexes of million atoms at relevant biological time scales is becoming accessible to the broad scientific community. That proves to be crucial for urgent responses against emergent diseases in real time. Yet, there are still issues to sort regarding the system setup so that molecular dynamics (MD) simulations can be run in a simple and standard way. Here, we introduce an optimized pipeline for building and simulating enveloped virus-like particles (VLP). First, the membrane packing problem is tackled with new features and optimized options in PACKMOL. This allows preparing accurate membrane models of thousands of lipids in the context of a VLP within a few hours using a single CPU. Then, the assembly of the VLP system is done within the multiscale framework of the coarse-grained SIRAH force field. Finally, the equilibration protocol provides a system ready for production MD simulations within a few days on broadly accessible GPU resources. The pipeline is applied to study the Zika virus as a test case for large biomolecular systems. The VLP stabilizes at approximately 0.5 μs of MD simulation, reproducing correlations greater than 0.90 against experimental density maps from cryo-electron microscopy. Detailed structural analysis of the protein envelope also shows very good agreement in root-mean-square deviations and B-factors with the experimental data. The level of details attained shows for the first time a possible role for anionic phospholipids in stabilizing the envelope. Combining an efficient and reliable setup procedure with an accurate coarse-grained force field provides a valuable pipeline for simulating arbitrary viral systems or subcellular compartments, paving the way toward whole-cell simulations.

Physics of viral dynamics

Viral capsids are often regarded as inert structural units, but in actuality they display fascinating dynamics during different stages of their life cycle. With the advent of single-particle approaches and high-resolution techniques, it is now possible to scrutinize viral dynamics during and after their assembly and during the subsequent development pathway into infectious viruses. In this Review, the focus is on the dynamical properties of viruses, the different physical virology techniques that are being used to study them, and the physical concepts that have been developed to describe viral dynamics.

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.

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP₆) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP₆ and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

This study shows that the HIV-1 capsid protein, in addition to its structural role, regulates reverse transcription, an essential metabolic process of the virus, by mediating the import of nucleotides. In addition, host cell metabolites such as inositol phosphates are recruited by the capsid to regulate viral DNA synthesis.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS

The introduction of accelerator devices such as graphics processing units (GPUs) has had profound impact on molecular dynamics simulations and has enabled order-of-magnitude performance advances using commodity hardware. To fully reap these benefits, it has been necessary to reformulate some of the most fundamental algorithms, including the Verlet list, pair searching, and cutoffs. Here, we present the heterogeneous parallelization and acceleration design of molecular dynamics implemented in the GROMACS codebase over the last decade. The setup involves a general cluster-based approach to pair lists and non-bonded pair interactions that utilizes both GPU and central processing unit (CPU) single instruction, multiple data acceleration efficiently, including the ability to load-balance tasks between CPUs and GPUs. The algorithm work efficiency is tuned for each type of hardware, and to use accelerators more efficiently, we introduce dual pair lists with rolling pruning updates. Combined with new direct GPU-GPU communication and GPU integration, this enables excellent performance from single GPU simulations through strong scaling across multiple GPUs and efficient multi-node parallelization.

Atomic-scale characterization of mature HIV-1 capsid stabilization by inositol hexakisphosphate (IP6)

Inositol hexakisphosphates (IP6) are cellular cofactors that promote the assembly of mature capsids of HIV. These negatively charged molecules coordinate an electropositive ring of arginines at the center of pores distributed throughout the capsid surface. Kinetic studies indicate that the binding of IP6 increases the stable lifetimes of the capsid by several orders of magnitude from minutes to hours. Using all-atom molecular dynamics simulations, we uncover the mechanisms that underlie the unusually high stability of mature capsids in complex with IP6 We find that capsid hexamers and pentamers have differential binding modes for IP6 Ligand density calculations show three sites of interaction with IP6 including at a known capsid inhibitor binding pocket. Free energy calculations demonstrate that IP6 preferentially stabilizes pentamers over hexamers to enhance fullerene modes of assembly. These results elucidate the molecular role of IP6 in stabilizing and assembling the retroviral capsid.

Modeling and scale-bridging using machine learning: nanoconfinement effects in porous media

Fine-scale models that represent first-principles physics are challenging to represent at larger scales of interest in many application areas. In nanoporous media such as tight-shale formations, where the typical pore size is less than 50 nm, confinement effects play a significant role in how fluids behave. At these scales, fluids are under confinement, affecting key properties such as density, viscosity, adsorption, etc. Pore-scale Lattice Boltzmann Methods (LBM) can simulate flow in complex pore structures relevant to predicting hydrocarbon production, but must be corrected to account for confinement effects. Molecular dynamics (MD) can model confinement effects but is computationally expensive in comparison. The hurdle to bridging MD with LBM is the computational expense of MD simulations needed to perform this correction. Here, we build a Machine Learning (ML) surrogate model that captures adsorption effects across a wide range of parameter space and bridges the MD and LBM scales using a relatively small number of MD calculations. The model computes upscaled adsorption parameters across varying density, temperature, and pore width. The ML model is 7 orders of magnitude faster than brute force MD. This workflow is agnostic to the physical system and could be generalized to further scale-bridging applications.

Intrinsic curvature of the HIV-1 CA hexamer underlies capsid topology and interaction with cyclophilin A

The mature retrovirus capsid consists of a variably curved lattice of capsid protein (CA) hexamers and pentamers. High-resolution structures of the curved assembly, or in complex with host factors, have not been available. By devising cryo-EM methodologies for exceedingly flexible and pleomorphic assemblies, we have determined cryo-EM structures of apo-CA hexamers and in complex with cyclophilin A (CypA) at near-atomic resolutions. The CA hexamers are intrinsically curved, flexible and asymmetric, revealing the capsomere and not the previously touted dimer or trimer interfaces as the key contributor to capsid curvature. CypA recognizes specific geometries of the curved lattice, simultaneously interacting with three CA protomers from adjacent hexamers via two noncanonical interfaces, thus stabilizing the capsid. By determining multiple structures from various helical symmetries, we further revealed the essential plasticity of the CA molecule, which allows formation of continuously curved conical capsids and the mechanism of capsid pattern sensing by CypA.

Coronavirus through Delaware's Computational Microscope

The Perilla/Hadden-Perilla research team at the University of Delaware presents an overview of computational structural biology, their efforts to model the SARS-CoV-2 viral particle, and their perspective on how their work and training endeavors can contribute to public health.

Experimentally-driven protein structure modeling

Revolutions in natural and exact sciences started at the dawn of last century have led to the explosion of theoretical, experimental, and computational approaches to determine structures of molecules, complexes, as well as their rich conformational dynamics. Since different experimental methods produce information that is attributed to specific time and length scales, corresponding computational methods have to be tailored to these scales and experiments. These methods can be then combined and integrated in scales, hence producing a fuller picture of molecular structure and motion from the "puzzle pieces" offered by various experiments. Here, we describe a number of computational approaches to utilize experimental data to glance into structure of proteins and understand their dynamics. We will also discuss the limitations and the resolution of the constraints-based modeling approaches. SIGNIFICANCE: Experimentally-driven computational structure modeling and determination is a rapidly evolving alternative to traditional approaches for molecular structure determination. These new hybrid experimental-computational approaches are proving to be a powerful microscope to glance into the structural features of intrinsically or partially disordered proteins, dynamics of molecules and complexes. In this review, we describe various approaches in the field of experimentally-driven computational structure modeling.

Molecular Basis of Differential Stability and Temperature Sensitivity of ZIKA versus Dengue Virus Protein Shells

Rapid spread of ZIKA virus (ZIKV) and its association with severe birth defects have raised worldwide concern. Recent studies have shown that ZIKV retains its infectivity and remains structurally stable at temperatures up to 40 °C, unlike dengue and other flaviviruses. In spite of recent cryo-EM structures that showed similar architecture of ZIKA and dengue virus (DENV) E protein shells, little is known that makes ZIKV so temperature insensitive. Here, we attempt to unravel the molecular basis of greater thermal stability of ZIKV over DENV2 by executing atomistic molecular dynamics (MD) simulations on the viral E protein shells at 37 °C. Our results suggest that ZIKA E protein shell retains its structural integrity through stronger inter-raft communications facilitated by a series of electrostatic and H-bonding interactions among multiple inter-raft residues. In comparison, the DENV2 E protein shell surface was loosly packed that exhibited holes at all 3-fold vertices, in close agreement with another EM structure solved at 37 °C. The residue-level information obtained from our study could pave way for designing small molecule inhibitors and specific antibodies to inhibit ZIKV E protein assembly and membrane fusion.

A Review on Applications of Computational Methods in Drug Screening and Design

Drug development is one of the most significant processes in the pharmaceutical industry. Various computational methods have dramatically reduced the time and cost of drug discovery. In this review, we firstly discussed roles of multiscale biomolecular simulations in identifying drug binding sites on the target macromolecule and elucidating drug action mechanisms. Then, virtual screening methods (e.g., molecular docking, pharmacophore modeling, and QSAR) as well as structure- and ligand-based classical/de novo drug design were introduced and discussed. Last, we explored the development of machine learning methods and their applications in aforementioned computational methods to speed up the drug discovery process. Also, several application examples of combining various methods was discussed. A combination of different methods to jointly solve the tough problem at different scales and dimensions will be an inevitable trend in drug screening and design.

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.

A new analogue of islet neogenesis associated protein with higher structural and plasma stability

Islet Neogenesis Associated Protein pentadecapeptide (INGAP-PP) increases β-cell mass and function in experimental animals. A short clinical trial also yielded promising results. However, HTD4010, a new peptide derived from INGAP-PP, was developed in order to optimize its specific effects by minimizing its side effects. To study and compare the tertiary structure, stability dynamics, and plasma stability of HTD4010, an INGAP-PP analogue. Both peptides were pre-incubated in human, rat and mouse plasma at 37 °C, and their presence was identified and quantified by high performance liquid chromatography at different time-points. GROMACS 2019 package and the Gromos 54A7 force field were used to evaluate overall correlated motion of the peptide molecule during molecular dynamics simulation by essential dynamics. HTD4010 exhibited significantly larger plasma stability than INGAP-PP, and its structural stability was almost 3.36-fold higher than INGAP-PP. These results suggest that HTD4010 may facilitate longer tissue interaction, thereby developing higher potential biological effects. If so, HTD4010 may become a promising therapeutic agent to treat people with diabetes. Communicated by Ramaswamy H. Sarma.

Computational Nanoscopy of Tight Junctions at the Blood-Brain Barrier Interface

The selectivity of the blood-brain barrier (BBB) is primarily maintained by tight junctions (TJs), which act as gatekeepers of the paracellular space by blocking blood-borne toxins, drugs, and pathogens from entering the brain. The BBB presents a significant challenge in designing neurotherapeutics, so a comprehensive understanding of the TJ architecture can aid in the design of novel therapeutics. Unraveling the intricacies of TJs with conventional experimental techniques alone is challenging, but recently developed computational tools can provide a valuable molecular-level understanding of TJ architecture. We employed the computational methods toolkit to investigate claudin-5, a highly expressed TJ protein at the BBB interface. Our approach started with the prediction of claudin-5 structure, evaluation of stable dimer conformations and nanoscale assemblies, followed by the impact of lipid environments, and posttranslational modifications on these claudin-5 assemblies. These led to the study of TJ pores and barriers and finally understanding of ion and small molecule transport through the TJs. Some of these in silico, molecular-level findings, will need to be corroborated by future experiments. The resulting understanding can be advantageous towards the eventual goal of drug delivery across the BBB. This review provides key insights gleaned from a series of state-of-the-art nanoscale simulations (or computational nanoscopy studies) performed on the TJ architecture.

High-Performance Analysis of Biomolecular Containers to Measure Small-Molecule Transport, Transbilayer Lipid Diffusion, and Protein Cavities

Compartmentalization is a central theme in biology. Cells are composed of numerous membrane-enclosed structures, evolved to facilitate specific biochemical processes; viruses act as containers of genetic material, optimized to drive infection. Molecular dynamics simulations provide a mechanism to study biomolecular containers and the influence they exert on their environments; however, trajectory analysis software generally lacks knowledge of container interior versus exterior. Further, many relevant container analyses involve large-scale particle tracking endeavors, which may become computationally prohibitive with increasing system size. Here, a novel method based on 3-D ray casting is presented, which rapidly classifies the space surrounding biomolecular containers of arbitrary shape, enabling fast determination of the identities and counts of particles (e.g., solvent molecules) found inside and outside. The method is broadly applicable to the study of containers and enables high-performance characterization of properties such as solvent density, small-molecule transport, transbilayer lipid diffusion, and topology of protein cavities. The method is implemented in VMD, a widely used simulation analysis tool that supports personal computers, clouds, and parallel supercomputers, including ORNL's Summit and Titan and NCSA's Blue Waters, where the method can be employed to efficiently analyze trajectories encompassing millions of particles. The ability to rapidly characterize the spatial relationships of particles relative to a biomolecular container over many trajectory frames, irrespective of large particle counts, enables analysis of containers on a scale that was previously unfeasible, at a level of accuracy that was previously unattainable.

Perspectives towards antiviral drug discovery against Ebola virus

Ebola virus disease (EVD), caused by Ebola viruses, resulted in more than 11 500 deaths according to a recent 2018 WHO report. With mortality rates up to 90%, it is nowadays one of the most deadly infectious diseases. However, no Food and Drug Administration‐approved Ebola drugs or vaccines are available yet with the mainstay of therapy being supportive care. The high fatality rate and absence of effective treatment or vaccination make Ebola virus a category‐A biothreat pathogen. Fortunately, a series of investigational countermeasures have been developed to control and prevent this global threat. This review summarizes the recent therapeutic advances and ongoing research progress from research and development to clinical trials in the development of small‐molecule antiviral drugs, small‐interference RNA molecules, phosphorodiamidate morpholino oligomers, full‐length monoclonal antibodies, and vaccines. Moreover, difficulties are highlighted in the search for effective countermeasures against EVD with additional focus on the interplay between available in silico prediction methods and their evidenced potential in antiviral drug discovery.

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

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

Mesoscale model of the assembly and cross-linking of HPV virus-like particles

We present a novel kinetic Monte Carlo model to simulate the real process time-scale of the assembly of Human Papillomavirus (HPV) virus-like particles (VLPs) incorporating the formation of intercapsomeric disulfide bonds. The objective was to develop insights into the underlying mechanisms of HPV VLP assembly and cross-linking during in vitro production of the HPV vaccine. The model integrates actual experimental data and detailed information of VLP geometrical structure in microscopic mechanistic steps. The principal novelty of this model is in the concurrent simulation of VLP assembly and cross-linking including a variable for spatial angular arrangement of capsomeres during their assembly that affects the overall rates of VLP assembly and cross-linking. The cross-linking modeled by using the mechanistic probability rules between involved cysteine residues. The model was utilized to better understand the actual process data and check on the hypothesis related to factors affecting the rates of HPV growth and maturation.

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

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

Maturation of retroviruses

During retrovirus maturation, cleavage of the precursor structural Gag polyprotein by the viral protease induces architectural rearrangement of the virus particle from an immature into a mature, infectious form. The structural rearrangement encapsidates the viral RNA genome in a fullerene capsid, producing a diffusible viral core that can initiate infection upon entry into the cytoplasm of a host cell. Maturation is an important therapeutic window against HIV-1. In this review, we highlight recent breakthroughs in understanding of the structures of retroviral immature and mature capsid lattices that define the boundary conditions of maturation and provide novel insights on capsid transformation. We also discuss emerging insights on encapsidation of the viral genome in the mature capsid, as well as remaining questions for further study.

Structure-based insights into the mechanism of nucleotide import by HIV-1 capsid

Nucleotide import is critical to the reverse transcription and thus a key step in the life cycle of HIV-1 virus. In this work, boosted by the recently determined atomic structure of the whole HIV-1 capsid and our understanding of its dynamics, we carry out a computational structure-based study of the mechanism of nucleotide import by HIV-1 capsid. Our analysis reveals a number of novel insights into the mechanism. Specifically, our structure analysis indicates that arginine 18 side chains use a bind-translocate-release mechanism to bind and guide the long nucleotide to pass through the arginine ring before releasing it into the capsid. We discover that LYS 25 plays a critical role in defining the directionality of the transport, in facilitating the disassociation of nucleotides from the arginine ring, and in relaying them further into the interior of the capsid. We predict that mutations of LYS 25 to other residues, assuming that they still result in viable viruses, should reverse the directionality of the transport and significantly lower the infectivity.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

MS2 bacteriophage capsid studied using all-atom molecular dynamics

The all-atom model of an MS2 bacteriophage particle without its genome (the capsid) was built using high-resolution cryo-electron microscopy (EM) measurements for initial conformation. The structural characteristics of the capsid and the dynamics of the surrounding solution were examined using molecular dynamics simulation. The model demonstrates the overall preservation of the cryo-EM structure of the capsid at physiological conditions (room temperature and ions composition). The formation of a dense anion layer near the inner surface and a diffuse cation layer near the outer surface of the capsid was detected. The flow of water molecules and ions across the capsid through its pores were quantified, which was considerable for water and substantial for ions.

The dynamical interplay between a megadalton peptide nanocage and solutes probed by microsecond atomistic MD; implications for design

Better understanding of the dynamics of protein-based supramolecular capsids can be applied to synthetic biology and biotechnology.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

Dynamic regulation of HIV-1 capsid interaction with the restriction factor TRIM5α identified by magic-angle spinning NMR and molecular dynamics simulations

The host factor protein TRIM5α plays an important role in restricting the host range of HIV-1, interfering with the integrity of the HIV-1 capsid. TRIM5 triggers an antiviral innate immune response by functioning as a capsid pattern recognition receptor, although the precise mechanism by which the restriction is imposed is not completely understood. Here we used an integrated magic-angle spinning nuclear magnetic resonance and molecular dynamics simulations approach to characterize, at atomic resolution, the dynamics of the capsid's hexameric and pentameric building blocks, and the interactions with TRIM5α in the assembled capsid. Our data indicate that assemblies in the presence of the pentameric subunits are more rigid on the microsecond to millisecond timescales than tubes containing only hexamers. This feature may be of key importance for controlling the capsid's morphology and stability. In addition, we found that TRIM5α binding to capsid induces global rigidification and perturbs key intermolecular interfaces essential for higher-order capsid assembly, with structural and dynamic changes occurring throughout the entire CA polypeptide chain in the assembly, rather than being limited to a specific protein-protein interface. Taken together, our results suggest that TRIM5α uses several mechanisms to destabilize the capsid lattice, ultimately inducing its disassembly. Our findings add to a growing body of work indicating that dynamic allostery plays a pivotal role in capsid assembly and HIV-1 infectivity.

All-Atom Molecular Dynamics Simulations of Whole Viruses

Classical molecular dynamics modeling of whole viruses or their capsids in explicit water is discussed, and known examples from the literature are analyzed. Only works on all-atom modeling in explicit water are included. Physical chemistry of the whole system is the focus, which includes the structure and dynamics of the biomolecules as well as water and ion behavior in and around the virus particle. It was demonstrated that in most investigations molecular phenomena that currently can not be studied experimentally are successfully reproduced and explained by the simulations. These include, for example, transport and distribution of ions inside viruses that ultimately connected to their stability, the hydrodynamic pressure in the capsid related to viruses' elastic properties, the role of metal ions in virus swelling, and others. Current and future tendencies in the development of all-atom virus simulations are outlined.

All-atom normal mode dynamics of HIV-1 capsid

Dynamics of biomolecular assemblies offer invaluable insights into their functional mechanisms. For extremely large biomolecular systems, such as HIV-1 capsid that has nearly 5 millions atoms, obtaining its normal mode dynamics using even coarse-grained models can be a challenging task. In this work, we have successfully carried out a normal mode analysis of an entire HIV-1 capsid in full all-atom details. This is made possible through our newly developed BOSE (Block of Selected Elasticity) model that is founded on the principle of resonance discovered in our recent work. The resonance principle makes it possible to most efficiently compute the vibrations of a whole capsid at any given frequency by projecting the motions of component capsomeres into a narrow subspace. We have conducted also assessments of the quality of the BOSE modes by comparing them with benchmark modes obtained directly from the original Hessian matrix. Our all-atom normal mode dynamics study of the HIV-1 capsid reveals the dynamic role of the pentamers in stabilizing the capsid structure and is in agreement with experimental findings that suggest capsid disassembly and uncoating start when the pentamers become destabilized. Our results on the dynamics of hexamer pores suggest that nucleotide transport should take place mostly at hexamers near pentamers, especially at the larger hemispherical end.

Supramolecular assemblies are large biomolecular complexes composed of hundreds or even thousands of protein chains. They function as molecular machines or as large containers that store or facilitate the chemical reactions of other molecules. Whatever they do, their functional mechanisms are tightly linked to their structures and intrinsic dynamics. Recently, due to breakthroughs in experimental techniques, many supramolecular assemblies have been determined, such as the capsid of human immunodeficiency virus (HIV) that is composed of nearly 5 millions of atoms. Computational studies of these systems are challenging due to their extremely large sizes. In this work, we have successfully carried out a dynamics study of an entire HIV capsid in full all-atom details. Our study reveals new insights into the dynamics of the N-terminal loops, the stabilizing role of the pentamers, and where the nucleotide transport may take place.

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.