The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations

pH induced structural and conformational changes in nucleocapsid protein leads to intermediate like conformation: a biophysical and computational approach

Nucleocapsid protein (N) of SARS-CoV-2 is a multivalent protein, which is responsible for viral replication, assembly, packaging and modulates host immune response. In this study, we report conformational measurements of N protein at different pH by observing transition in secondary and tertiary structural contents by biophysical and computational approaches. Spectroscopic measurements revealed that N protein loses its secondary and tertiary structure at extreme acidic pH while maintaining its native conformation at mild acidic and alkaline pH. Molecular dynamics simulation studies validated spectroscopic findings. Secondary structure estimation confirmed circular dichroism (CD) findings that participation of total number of average residues in formation of native structure is higher at physiological pH, and coil percentage is higher at acidic pH. In molten globule (MG) state, secondary structure is conserved but here, CD data reveal more random structure at low pH. In pre-MG, ANS (8-anilino-1-napthalene sulfonate) binds weakly to protein as compared to MG but here, ANS binds strongly to protein. All the above-mentioned findings suggested formation of intermediary like state at low pH, which can be attributed to an off-pathway species. Unravelling structural characteristics of N protein will help understand phase-separation, protein-protein interaction and host-immune response modulation behaviour, which will eventually help in designing novel therapeutic target against COVID-19.

Titr-DMD-A Rapid, Coarse-Grained Quasi-All-Atom Constant pH Molecular Dynamics Framework

The pH-dependence of enzyme fold stability and catalytic activity is a fundamentally dynamic, structural property which is difficult to study. The challenges and expense of investigating dynamic, atomic scale behavior experimentally means that computational methods, particularly constant pH molecular dynamics (CpHMD), are well situated tools for this. However, these methods often struggle with affordable sampling of sufficiently long time scales while also obtaining accurate pKa prediction and verifying the structures they generate. We introduce Titr-DMD, an affordable CpHMD method that combines the quasi-all-atom coarse-grained discrete molecular dynamics (DMD) method for conformational sampling with Propka for pKa prediction, to circumvent these issues. The combination enables rapid sampling on limited computational resources, while simulations are still performed on the atomic scale. We benchmark the method on a set of proteins with experimentally attested pKa and on the pH triggered conformational change in a staphylococcal nuclease mutant, a rare experimental study of such behavior. Our results show Titr-DMD to be an effective and inexpensive method to study pH-coupled protein dynamics.

Computational modelling of flavivirus dynamics: The ins and outs

Enveloped viruses such as the flaviviruses represent a significant burden to human health around the world, with hundreds of millions of people each year affected by dengue alone. In an effort to improve our understanding of the molecular basis for the infective mechanisms of these viruses, extensive computational modelling approaches have been applied to elucidate their conformational dynamics. Multiscale protocols have been developed to simulate flavivirus envelopes in close accordance with biophysical data, in particular derived from cryo-electron microscopy, enabling high-resolution refinement of their structures and elucidation of the conformational changes associated with adaptation both to host environments and to immunological factors such as antibodies. Likewise, integrative modelling efforts combining data from biophysical experiments and from genome sequencing with chemical modification are providing unparalleled insights into the architecture of the previously unresolved nucleocapsid complex. Collectively, this work provides the basis for the future rational design of new antiviral therapeutics and vaccine development strategies targeting enveloped viruses.

Molecular Changes in Dengue Envelope Protein Domain III upon Interaction with Glycosaminoglycans

Dengue fever is a rapidly emerging vector-borne viral disease with a growing global burden of approximately 390 million new infections per annum. The Dengue virus (DENV) is a flavivirus spread by female mosquitos of the aedes genus, but the mechanism of viral endocytosis is poorly understood at a molecular level, preventing the development of effective transmission blocking vaccines (TBVs). Recently, glycosaminoglycans (GAGs) have been identified as playing a role during initial viral attachment through interaction with the third domain of the viral envelope protein (EDIII). Here, we report a systematic study investigating the effect of a range of biologically relevant GAGs on the structure and oligomeric state of recombinantly generated EDIII. We provide novel in situ biophysical evidence that heparin and chondroitin sulphate C induce conformational changes in EDIII at the secondary structure level. Furthermore, we report the ability of chondroitin sulphate C to bind EDIII and induce higher-order dynamic molecular changes at the tertiary and quaternary structure levels which are dependent on pH, GAG species, and the GAG sulphation state. Lastly, we conducted ab initio modelling of Small Angle Neutron Scattering (SANS) data to visualise the induced oligomeric state of EDIII caused by interaction with chondroitin sulphate C, which may aid in TBV development.

Multiple Virtual Screening Strategies for the Discovery of Novel Compounds Active Against Dengue Virus: A Hit Identification Study

Dengue infection is caused by a mosquito-borne virus, particularly in children, which may even cause death. No effective prevention or therapeutic agents to cure this disease are available up to now. The dengue viral envelope (E) protein was discovered to be a promising target for inhibition in several steps of viral infection. Structure-based virtual screening has become an important technique to identify first hits in a drug screening process, as it is possible to reduce the number of compounds to be assayed, allowing to save resources. In the present study, pharmacophore models were generated using the common hits approach (CHA), starting from trajectories obtained from molecular dynamics (MD) simulations of the E protein complexed with the active inhibitor, flavanone (FN5Y). Subsequently, compounds presented in various drug databases were screened using the LigandScout 4.2 program. The obtained hits were analyzed in more detail by molecular docking, followed by extensive MD simulations of the complexes. The highest-ranked compound from this procedure was then synthesized and tested on its inhibitory efficiency by experimental assays.

The structural basis of acid resistance in Mycobacterium tuberculosis: insights from multiple pH regime molecular dynamics simulations

The dormant Mycobacterium tuberculosis is evolved to develop the tolerance against the acidification of phagolysosome by the action of gamma interferon. The molecular mechanism responsible for the development of the resistance towards the acidic conditions in M. tuberculosis is not fully understood. Therefore, the current analysis was performed which studies the mechanism of acid tolerance by correlating the alteration in the protonation state of conserved residues in virulent proteins with changes in their folding states. The pH dependencies of proteins were studied using an efficient computational scheme which enables the understanding of their conformational behavior by molecular dynamics (MD) simulations. The adopted methodology involves cyclically updating of the ionization states of titrable residues in the studied proteins with conventional MD steps, which were applied to the newly generated ionization configuration. Significant pH-dependent protein structural stability parameters consistent with the changes of the protonation states of conserved residues were observed. Among the studied proteins, the peptidoglycan binding protein ompATB, carboxylesterase LipF and two-component systems' transcriptional regulator PhoP showed highest structural conservation in the observed acidic pH range throughout the course of MD simulations. The current study provides a better understanding of acid tolerance mechanisms present in M. tuberculosis and can facilitate the drug development strategies against the dormant protein targets.Communicated by Ramaswamy H. Sarma.

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

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

Study of the mechanism of protonated histidine-induced conformational changes in the Zika virus dimeric envelope protein using accelerated molecular dynamic simulations

The Zika virus has drawn worldwide attention because of the epidemic diseases it causes. It is a flavivirus that has an icosahedral protein shell constituted by an envelope glycoprotein (E-protein) and membrane protein (M-protein) in the mature virion. The multistep process of membrane fusion to infect the host cell is pH-induced. To understand the mechanism of the conformational changes in the (E-M)₂ protein homodimer embedded in the membrane, two 200-ns accelerated dynamic simulations were performed under different pH conditions. The low pH condition weakens the interactions and correlations in both E-protein monomers and in the E-M heterodimer. The highly conserved residues, His249, His288, His323 and His446, are protonated under low pH conditions and play key roles in driving the fusion process. The analysis and discussion in this study may provide some insight into the molecular mechanism of Zika virus infection.

Receptor-Guided De Novo Design of Dengue Envelope Protein Inhibitors

Inhibitor design associated with the dynamics of dengue envelope protein at pre-fusion stage is a prominent strategy to interfere fusion transition of dengue virus with the host cell membrane. Receptor-guided de novo inhibitors were designed based on the knowledge of co-crystallized detergent, β-octyl glucoside. Pharmacophore features distribution showed the preference of aromatic groups with H bonding features connected to aliphatic bulky group as the skeleton for inhibitor design. Molecular dynamic simulations revealed (2R)-2-[(6-amino-1-oxohexan-2-yl)amino]-4-[6-(4-phenylpiperidine-1-yl)-1,2-benzoxazol-3-yl]butanoate as the probable binder which developed extensive conservative interactions despite the local pocket residues movements especially from kl β-hairpin, the key structural unit for initiating conformational changes required for fusion transition. The electronic and hydrophobic potentials also indicated that butanoate molecule as the initial lead for envelope protein inhibitors.