Molecular simulations to delineate functional conformational transitions in the HCV polymerase.

Hepatitis C virus (HCV) is a global health concern for which there is no vaccine available. The HCV polymerase is responsible for the critical function of replicating the RNA genome of the virus. Transitions between at least two conformations (open and closed) are necessary to allow the enzyme to replicate RNA. In this study, molecular dynamic simulations were initiated from multiple crystal structures to understand the free energy landscape (FEL) explored by the enzyme as it interconverts between these conformations. Our studies reveal the location of distinct states within the FEL as well as the molecular interactions associated with these states. Specific hydrogen bonds appear to play a key role in modulating conformational transitions. This knowledge is essential to elucidate the role of these conformations in replication and may also be valuable in understanding the basis by which this enzyme is inhibited by small molecules. © 2016 Wiley Periodicals, Inc.

Allosteric inhibitors have distinct effects, but also common modes of action, in the HCV polymerase.

The RNA-dependent RNA polymerase from the Hepatitis C Virus (gene product NS5B) is a validated drug target because of its critical role in genome replication. There are at least four distinct allosteric sites on the polymerase to which several small molecule inhibitors bind. In addition, numerous crystal structures have been solved with different allosteric inhibitors bound to the polymerase. However, the molecular mechanisms by which these small molecules inhibit the enzyme have not been fully elucidated. There is evidence that allosteric inhibitors alter the intrinsic motions and distribution of conformations sampled by the enzyme. In this study we use molecular dynamics simulations to understand the structural and dynamic changes that result when inhibitors are bound at three different allosteric binding sites on the enzyme. We observe that ligand binding at each site alters the structure and dynamics of NS5B in a distinct manner. Nonetheless, our studies also highlight commonalities in the mechanisms of action of the different inhibitors. Each inhibitor alters the conformational states sampled by the enzyme, either by rigidifying the enzyme and preventing transitions between functional conformational states or by destabilizing the enzyme and preventing functionally relevant conformations from being adequately sampled. By illuminating the molecular mechanisms of allosteric inhibition, these studies delineate the intrinsic functional properties of the enzyme and pave the way for designing novel and more effective polymerase inhibitors. This information may also be important to understand how allosteric regulation occurs in related viral polymerases and other enzymes.

Dual Allosteric Inhibitors Jointly Modulate Protein Structure and Dynamics in the Hepatitis C Virus Polymerase.

The hepatitis C virus (HCV) infects close to 200 million people globally, resulting in a significant need for effective HCV therapies. The HCV polymerase (gene product NS5B) is a valuable target for therapeutics because of its role in replicating the viral genome. Various studies have identified inhibitors for this enzyme, including non-nucleoside inhibitors (NNIs) that bind distal to the enzyme active site. Recently, it has been shown that simultaneously challenging the enzyme with two NNIs results in enhanced inhibition relative to that observed after challenge with individual inhibitors, suggesting that employing multiple NNIs might be the basis of more effective therapeutics. Nevertheless, the molecular mechanisms responsible for this enhanced inhibition remain unclear. We employ molecular dynamics simulations to determine the origin of enhanced inhibition when two NNIs bind to NS5B. Our results suggest that nonoverlapping NNI sites are compatible with simultaneous binding of dual NNIs. We observe that both inhibitors act in concert to induce novel enzyme conformations and dynamics, allowing us to identify molecular mechanisms underlying enhanced inhibition of NS5B. This knowledge will be useful in optimizing combinations of NNIs to target NS5B, helping to prevent the acquisition of viral resistance that remains a significant barrier to the development of HCV therapeutics.

Inhibitors for the hepatitis C virus RNA polymerase explored by SAR with advanced machine learning methods.

Hepatitis C virus (HCV) is a global health challenge, affecting approximately 200 million people worldwide. In this study we developed SAR models with advanced machine learning classifiers Random Forest and k Nearest Neighbor Simulated Annealing for 679 small molecules with measured inhibition activity for NS5B genotype 1b. The activity was expressed as a binary value (active/inactive), where actives were considered molecules with IC50 ≤0.95 µM. We applied our SAR models to various drug-like databases and identified novel chemical scaffolds for NS5B inhibitors. Subsequent in vitro antiviral assays suggested a new activity for an existing prodrug, Candesartan cilexetil, which is currently used to treat hypertension and heart failure but has not been previously tested for anti-HCV activity. We also identified NS5B inhibitors with two novel non-nucleoside chemical motifs.

Intramolecular hydrogen bonding in ortho-substituted arylamide oligomers: a computational and experimental study of ortho-fluoro- and ortho-chloro-N-methylbenzamides.

As a part of our systematic study of foldamer structural elements, we analyze and quantify the conformational behavior of two model compounds based on a frequently used class of aromatic oligoamide building blocks. Combining computational and NMR approaches, we investigate ortho-fluoro- and ortho-chloro-N-methylbenzamide. Our results indicate that the -F substituent in an ortho position can be used to fine-tune the rigidity of the oligomer backbone. It provides a measurably attenuated but still considerably strong hydrogen bond (H-bond) to the peptide group proton when compared to the -OCH3 substituent in the same position. On the other hand, the ortho-Cl substituent does not impose significant restrictions on the flexibility of the backbone. Its effect on the final shape of an oligomer is likely governed by its size rather than by noncovalent intramolecular interactions. Furthermore, the effect of solvent on the conformational preferences of these building blocks has been quantified. The number of intramolecularly H-bonded conformations decreases significantly when going from nonprotic to protic environments. This study will facilitate rational design of novel arylamide foldamers.

Computational predictions suggest that structural similarity in viral polymerases may lead to comparable allosteric binding sites.

The identification of ligand-binding sites is often the first step in drug targeting and design. To date there are numerous computational tools available to predict ligand binding sites. These tools can guide or mitigate the need for experimental methods to identify binding sites, which often require significant resources and time. Here, we evaluate four ligand-binding site predictor (LBSP) tools for their ability to predict allosteric sites within the Hepatitis C Virus (HCV) polymerase. Our results show that the LISE LBSP is able to identify all three target allosteric sites within the HCV polymerase as well as a known allosteric site in the Coxsackievirus polymerase. LISE was then employed to identify novel binding sites within the polymerases of the Dengue, West Nile, and Foot-and-mouth Disease viruses. Our results suggest that all three viral polymerases have putative sites that share structural or chemical similarities with allosteric pockets of the HCV polymerase. Thus, these binding locations may represent an evolutionarily conserved structural feature of several viral polymerases that could be exploited for the development of small molecule therapeutics.

Amino, Ammonio and Aminioethenes: A Theoretical Study of their Structure and Energetics.

We have performed high level ab initio quantum mechanical calculations for aminoethene and the three isomeric 1,1- (Z)- or (E)-1,2-diaminoethenes as well as their singly and doubly charged cations derived by loss of electrons and/or upon protonation. Gas phase molecular structures were computed at the MP2/6-311+G(3df,2p) level. Standard molar enthalpies of formation in the gas phase, at T = 298.15 K were estimated using the G3 composite method and atomization, isodesmic and homodesmotic reactions. Other energetic parameters were also calculated at the G3 level: proton affinities, basicities and adiabatic ionization enthalpies. Theoretical and experimental data are compared. The reported experimental data refer only to aminoethene wherein the standard molar enthalpy of formation has a considerable uncertainty, although the molecular structure is well-established. There are no such data, neither structural nor thermochemical, for any of the three isomeric diaminoethenes. Isoelectronic comparisons are made. For example, the diprotonated diaminoethenes are isoelectronic to isobutene and (Z)- and (E)-butene, while the doubly ionized diaminoethenes are likewise related to trimethylenemethane and 1,3-butadiene.