Residue-level resolution of alphavirus envelope protein interactions in pH-dependent fusion

Resolving coupled pH titrations using alchemical free energy calculations

In a protein, nearby titratable sites can be coupled: the (de)protonation of one may affect the other. The degree of this interaction depends on several factors and can influence the measured p K a . Here, we derive a formalism based on double free energy differences ( Δ Δ G ) for quantifying the individual site p K a values of coupled residues. As Δ Δ G values can be obtained by means of alchemical free energy calculations, the presented approach allows for a convenient estimation of coupled residue p K a s in practice. We demonstrate that our approach and a previously proposed microscopic p K a formalism, can be combined with alchemical free energy calculations to resolve pH-dependent protein p K a values. Toy models and both, regular and constant-pH molecular dynamics simulations, alongside experimental data, are used to validate this approach. Our results highlight the insights gleaned when coupling and microstate probabilities are analyzed and suggest extensions to more complex enzymatic contexts. Furthermore, we find that naïvely computed p K a values that ignore coupling, can be significantly improved when coupling is accounted for, in some cases reducing the error by half. In short, alchemical free energy methods can resolve the p K a values of both uncoupled and coupled residues.

In silico identification and in vitro antiviral validation of potential inhibitors against Chikungunya virus

The Chikungunya virus (CHIKV) has become endemic in the Africa, Asia and Indian subcontinent, with its continuous re-emergence causing a significant public health crisis. The unavailability of specific antivirals and vaccines against the virus has highlighted an urgent need for novel therapeutics. In the present study, we have identified small molecule inhibitors targeting the envelope proteins of the CHIKV to interfere with the fusion process, eventually inhibiting the cell entry of the virus particles. We employed high throughput computational screening of large datasets against two different binding sites in the E1-E2 dimer to identify potential candidate inhibitors. Among them, four high affinity inhibitors were selected to confirm their anti-CHIKV activity in the in vitro assay. Quercetin derivatives, Taxifolin and Rutin, binds to the E1-E2 dimer at different sites and display inhibition of CHIKV infection with EC₅₀ values 3.6 μM and 87.67 μM, respectively. Another potential inhibitor with ID ChemDiv 8015-3006 binds at both the target sites and shows anti-CHIKV activity at EC₅₀ = 41 μM. The results show dose-dependent inhibitory effects of Taxifolin, Rutin and ChemDiv 8015-3006 against the CHIKV with minimal cytotoxicity. In addition, molecular dynamics studies revealed the structural stability of these inhibitors at their respective binding sites in the E1-E2 protein. In conclusion, our study reports Taxifolin, Rutin and ChemDiv 8015-3006 as potential inhibitors of the CHIKV entry. Also, this study suggests a few potential candidate inhibitors which could serve as a template to design envelope protein specific CHIKV entry inhibitors.

Cooperative Chikungunya Virus Membrane Fusion and Its Substoichiometric Inhibition by CHK-152 Antibody

Chikungunya virus (CHIKV) presents a major burden on healthcare systems worldwide, but specific treatment remains unavailable. Attachment and fusion of CHIKV to the host cell membrane is mediated by the E1/E2 protein spikes. We used an in vitro single-particle fusion assay to study the effect of the potent, neutralizing antibody CHK-152 on CHIKV binding and fusion. We find that CHK-152 shields the virions, inhibiting interaction with the target membrane and inhibiting fusion. The analysis of the ratio of bound antibodies to epitopes implied that CHIKV fusion is a highly cooperative process. Further, dissociation of the antibody at lower pH results in a finely balanced kinetic competition between inhibition and fusion, suggesting a window of opportunity for the spike proteins to act and mediate fusion, even in the presence of the antibody.

pH Effects and Cooperativity among Key Titratable Residues for Escherichia coli Glycinamide Ribonucleotide Transformylase

Human glycinamide ribonucleotide transformylase (GAR Tfase) is a regulatory enzyme in the de novo purine biosynthesis pathway that has been extensively studied as an anticancer target. To some extent, inhibition of GAR Tfase selectively targets cancer cells over normal cells and inhibits purine formation and DNA replication. In this study, we investigated E. coli GAR Tfase, which shares high sequence similarity with the human GAR Tfase, and most functional residues are conserved. Herein, we aim to predict the pH-activity curve through a computational approach. We carried out pH-replica exchange molecular dynamics (pH-REMD) simulations to investigate pH-dependent functions such as structural changes, ligand binding, and catalytic activity. To compute the pH-activity curve, we identified the catalytic residues in specific protonation states, referred to as the catalytic competent protonation states (CCPS), which maintain the structure, keep ligands bound, and facilitate catalysis. Our computed population of CCPS with respect to pH matches well with the experimental pH-activity curve. To compute the microscopic pK_a values in the catalytically active conformation, we devised a thermodynamic model that considers the coupling between protonation states of CCPS residues and conformational states. These results allow us to correctly identify the general acid and base catalysts and interpret the pH-activity curve at an atomistic level.

The Phosphatidylserine Receptor TIM-1 Enhances Authentic Chikungunya Virus Cell Entry

Chikungunya virus (CHIKV) is a re-emerging, mosquito-transmitted, enveloped positive stranded RNA virus. Chikungunya fever is characterized by acute and chronic debilitating arthritis. Although multiple host factors have been shown to enhance CHIKV infection, the molecular mechanisms of cell entry and entry factors remain poorly understood. The phosphatidylserine-dependent receptors, T-cell immunoglobulin and mucin domain 1 (TIM-1) and Axl receptor tyrosine kinase (Axl), are transmembrane proteins that can serve as entry factors for enveloped viruses. Previous studies used pseudoviruses to delineate the role of TIM-1 and Axl in CHIKV entry. Conversely, here, we use the authentic CHIKV and cells ectopically expressing TIM-1 or Axl and demonstrate a role for TIM-1 in CHIKV infection. To further characterize TIM-1-dependent CHIKV infection, we generated cells expressing domain mutants of TIM-1. We show that point mutations in the phosphatidylserine binding site of TIM-1 lead to reduced cell binding, entry, and infection of CHIKV. Ectopic expression of TIM-1 renders immortalized keratinocytes permissive to CHIKV, whereas silencing of endogenously expressed TIM-1 in human hepatoma cells reduces CHIKV infection. Altogether, our findings indicate that, unlike Axl, TIM-1 readily promotes the productive entry of authentic CHIKV into target cells.

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.

Acidic pH-Induced Conformational Changes in Chikungunya Virus Fusion Protein E1: a Spring-Twisted Region in the Domain I-III Linker Acts as a Hinge Point for Swiveling Motion of Domains

Chikungunya virus (CHIKV), a mosquito-transmitted alphavirus, enters a cell through endocytosis, followed by viral and cell membrane fusion. The fusion protein, E1, undergoes an acid pH-induced pre- to postfusion conformation change during membrane fusion. As part of the conformation change, E1 dissociates from the receptor-binding protein, E2, and swivels its domains I and II over domain III to form an extended intermediate and then eventually to form a postfusion hairpin homotrimer. In this study, we tested if the domain I-III linker acts as a "hinge" for the swiveling motion of E1 domains. We found a conserved spring-twisted structure in the linker, stabilized by a salt bridge between a conserved arginine-aspartic acid pair, as a "hinge point" for domain swiveling. Molecular dynamics (MD) simulation of the CHIKV E1 or E2-E1 structure predicted that the spring-twisted region untwists at pH 5.5. Corroborating the prediction, introduction of a "cystine staple" at the hinge point, replacing the conserved arginine-aspartic acid pair with cysteine residues, resulted in loss of fusion activity of E1. MD simulation also predicted domain I-III swiveling at acidic pH. We tested if breaking the His 331-Lys 16 H bond between domains I and III, seen only in the prefusion conformation, is important for domain swiveling. When domains I and III are "stapled" by introducing a disulfide bond in between, E1 showed loss of fusion activity, implying that domain I and III dissociation is a critical acid pH-induced step in membrane fusion. However, replacement of His 331 with an acidic residue did not affect the pH threshold for fusion, suggesting His 331 is not an acid-sensing residue.IMPORTANCE Aedes mosquito-transmitted viruses such as the Zika, dengue, and chikungunya viruses have spread globally. CHIKV, similar to many other enveloped viruses, enters cells in sequential steps: step 1 involves receptor binding followed by endocytosis, and step 2 involves viral-cell membrane fusion in the endocytic vesicle. The viral envelope surface protein, E1, performs membrane fusion. E1 is triggered to undergo conformational changes by acidic pH of the maturing endosome. Different domains of E1 rearrange during the pre- to postfusion conformation change. Using in silico analysis of the E1 structure and different biochemical experiments, we explained a structural mechanism of key conformational changes in E1 triggered by acidic pH. We noted two important structural changes in E1 at acidic pH. In the first, a spring-twisted region in a loop connecting two domains (I and III) untwists, bringing a swiveling motion of domains on each other. In the second, breaking of interactions between domains I and III and domain separation are required for membrane fusion. This knowledge will help devise new therapeutic strategies to block conformation changes in E1 and thus viral entry.

A molecular understanding of alphavirus entry

Alphaviruses cause severe human illnesses including persistent arthritis and fatal encephalitis. As alphavirus entry into target cells is the first step in infection, intensive research efforts have focused on elucidating aspects of this pathway, including attachment, internalization, and fusion. Herein, we review recent developments in the molecular understanding of alphavirus entry both in vitro and in vivo and how these advances might enable the design of therapeutics targeting this critical step in the alphavirus life cycle.

Novel Molecular Signatures of Chikungunya Virus in Puerto Rico

OBJECTIVE

The chikungunya virus (CHIKV) is an arthropod-borne Alphavirus transmitted to humans, primarily via Aedes mosquitoes. In Puerto Rico, the first locally transmitted infections were reported in May 2014. Although the virus strain in Puerto Rico is related to the Asian/American lineage, many autochthonous cases have emerged recently in the Caribbean region (including Puerto Rico), raising the question of how CHIKV will evolve and adapt in PR. Taking the role of the envelope glycoprotein (E1) in viral evolution and transmission as a given, we analyzed the genetic diversity of the Puerto Rican (PR) E1 gene sequences and the phylogenetic relationships between those sequences and sequences from other parts of the world.

MATERIALS AND METHODS

To analyze the overall genetic variation, 772 nucleotide sequences of the E1 gene were obtained from the Virus Pathogen Resource (ViPR). A maximum-likelihood analysis was performed to determine the phylogenetic relationships between the PR sequences and sequences from 48 countries around the world.

RESULTS

The analysis of the E1 gene identified variations at 4 nucleotide positions, which included synonymous and nonsynonymous mutations. In addition, 2 nonsynonymous amino acid changes, T207M and S120L, were unique to the PR CHIKV sequences, and T155I was found to be shared by the PR (n = 3) and Colombia (n = 1) strains.

CONCLUSION

Our analysis of the E1 gene revealed new molecular signatures in PR CHIKV sequences, 1 of which was also found in Colombia. While studies have shown possible relationships between T98A and A226V with viral adaptation and spread, no other PR sequence contained these vector-adaptive mutations. Thus, constant monitoring of the virus remains an essential factor in the establishment of control strategies to track viral spread.

Proton-Coupled Conformational Allostery Modulates the Inhibitor Selectivity for β-Secretase

Many important pharmaceutical targets, such as aspartyl proteases and kinases, exhibit pH-dependent dynamics, functions and inhibition. Accurate prediction of their binding free energies is challenging because current computational techniques neglect the effects of pH. Here we combine free energy perturbation calculations with continuous constant pH molecular dynamics to explore the selectivity of a small-molecule inhibitor for β-secretase (BACE1), an important drug target for Alzheimer's disease. The calculations predicted identical affinity for BACE1 and the closely related cathepsin D at high pH; however, at pH 4.6 the inhibitor is selective for BACE1 by 1.3 kcal/mol, in excellent agreement with experiment. Surprisingly, the pH-dependent selectivity can be attributed to the protonation of His45, which allosterically modulates a loop-inhibitor interaction. Allosteric regulation induced by proton binding is likely common in biology; considering such allosteric sites could lead to exciting new opportunities in drug design.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Insights into Adenovirus Uncoating from Interactions with Integrins and Mediators of Host Immunity

Human adenoviruses are large (150 MDa) nonenveloped double-stranded DNA (dsDNA) viruses that cause acute respiratory, gastrointestinal and ocular infections. Despite these disease associations, adenovirus has aided basic and clinical research efforts through studies of its association with cells and as a target of host antiviral responses. This review highlights the knowledge of adenovirus disassembly and nuclear transport gleaned from structural, biophysical and functional analyses of adenovirus interactions with soluble and membrane-associated host molecules.

Efficient implementation of constant pH molecular dynamics on modern graphics processors

The treatment of pH sensitive ionization states for titratable residues in proteins is often omitted from molecular dynamics (MD) simulations. While static charge models can answer many questions regarding protein conformational equilibrium and protein-ligand interactions, pH-sensitive phenomena such as acid-activated chaperones and amyloidogenic protein aggregation are inaccessible to such models. Constant pH molecular dynamics (CPHMD) coupled with the Generalized Born with a Simple sWitching function (GBSW) implicit solvent model provide an accurate framework for simulating pH sensitive processes in biological systems. Although this combination has demonstrated success in predicting pKa values of protein structures, and in exploring dynamics of ionizable side-chains, its speed has been an impediment to routine application. The recent availability of low-cost graphics processing unit (GPU) chipsets with thousands of processing cores, together with the implementation of the accurate GBSW implicit solvent model on those chipsets (Arthur and Brooks, J. Comput. Chem. 2016, 37, 927), provide an opportunity to improve the speed of CPHMD and ionization modeling greatly. Here, we present a first implementation of GPU-enabled CPHMD within the CHARMM-OpenMM simulation package interface. Depending on the system size and nonbonded force cutoff parameters, we find speed increases of between one and three orders of magnitude. Additionally, the algorithm scales better with system size than the CPU-based algorithm, thus allowing for larger systems to be modeled in a cost effective manner. We anticipate that the improved performance of this methodology will open the door for broad-spread application of CPHMD in its modeling pH-mediated biological processes. © 2016 Wiley Periodicals, Inc.

Computational virology: From the inside out

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

Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations

The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.

All-Atom Molecular Dynamics of Virus Capsids as Drug Targets

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

Mayaro virus: a neglected arbovirus of the Americas

Mayaro virus is a neglected tropical arbovirus that causes a mild, self-limited febrile syndrome, sometimes accompanied by a highly incapacitating arthralgia. First isolated in Trinidad and Tobago in 1954, it was reported in several countries within the tropical regions of South and Central America. Human infections are accidental spillover of the enzootic cycle. Little epidemiological data are available due to inadequate surveillance and the generic nature of clinical manifestations resulting in the misdiagnosis with other viral fevers. Despite its restricted distribution, Mayaro fever may become a public health issue due to their urbanization potential. Accurate epidemiological data are urgently needed to access the real distribution of this virus guiding public health policies better.

Early Events in Chikungunya Virus Infection-From Virus Cell Binding to Membrane Fusion

Chikungunya virus (CHIKV) is a rapidly emerging mosquito-borne alphavirus causing millions of infections in the tropical and subtropical regions of the world. CHIKV infection often leads to an acute self-limited febrile illness with debilitating myalgia and arthralgia. A potential long-term complication of CHIKV infection is severe joint pain, which can last for months to years. There are no vaccines or specific therapeutics available to prevent or treat infection. This review describes the critical steps in CHIKV cell entry. We summarize the latest studies on the virus-cell tropism, virus-receptor binding, internalization, membrane fusion and review the molecules and compounds that have been described to interfere with virus cell entry. The aim of the review is to give the reader a state-of-the-art overview on CHIKV cell entry and to provide an outlook on potential new avenues in CHIKV research.

pKa values in proteins determined by electrostatics applied to molecular dynamics trajectories

For a benchmark set of 194 measured pKa values in 13 proteins, electrostatic energy computations are performed in which pKa values are computed by solving the Poisson-Boltzmann equation. In contrast to the previous approach of Karlsberg(+) (KB(+)) that essentially used protein crystal structures with variations in their side chain conformations, the present approach (KB2(+)MD) uses protein conformations from four molecular dynamics (MD) simulations of 10 ns each. These MD simulations are performed with different specific but fixed protonation patterns, selected to sample the conformational space for the different protonation patterns faithfully. The root-mean-square deviation between computed and measured pKa values (pKa RMSD) is shown to be reduced from 1.17 pH units using KB(+) to 0.96 pH units using KB2(+)MD. The pKa RMSD can be further reduced to 0.79 pH units, if each conformation is energy-minimized with a dielectric constant of εmin = 4 prior to calculating the electrostatic energy. The electrostatic energy expressions upon which the computations are based have been reformulated such that they do not involve terms that mix protein and solvent environment contributions and no thermodynamic cycle is needed. As a consequence, conformations of the titratable residues can be treated independently in the protein and solvent environments. In addition, the energy terms used here avoid the so-called intrinsic pKa and can therefore be interpreted without reference to arbitrary protonation states and conformations.