Impact of bound ssRNA length on allostery in the Dengue Virus NS3 helicase.

The presence of ATP is known to stimulate helicase activity of the Dengue Virus Non-structural protein 3 helicase (NS3h), and the presence of RNA stimulates NS3h ATPase activity, however this coupling is still mechanistically unclear. Here we use atomistic models and molecular dynamics simulations to evaluate the single-stranded RNA (ssRNA)-length dependence of the NS3h-ssRNA binding affinity and its modulation by bound ATP. Considering complexes with 7, 11, 16 and 26 nucleotides (nts), we observe that both the binding affinity and its modulation by bound ATP are augmented with increased ssRNA lengths. In models with at least 11 nts bound, the binding of ATP results in a shift from a tightly bound to a weakly bound state. We find that the weakly bound state persists during both the ADP-Pi- and ADP-bound stages of the catalytic cycle. We obtain the equilibrium association constants for NS3h binding to an ssRNA 10-mer in vitro, both in the absence and presence of ADP, which further support the alternation between tightly and weakly bound states during the catalytic cycle. The length of bound ssRNA is critical for understanding the NS3h-RNA interaction as well as how it is modulated during the catalytic cycle.

Thermodynamic and mechanistic analysis of the functional properties of dengue virus NS3 helicase.

The Dengue Virus (DENV) non-structural protein 3 (NS3) is a multi-functional protein critical in the viral life cycle. The DENV NS3 is comprised of a serine protease domain and a helicase domain. The helicase domain itself acts as a molecular motor, either translocating in a unidirectional manner along single-stranded RNA or unwinding double-stranded RNA, processes fueled by the hydrolysis of nucleoside triphosphates. In this brief review, we summarize our contributions and ongoing efforts to uncover the thermodynamic and mechanistic functional properties of the DENV NS3 as an NTPase and helicase.

Insights into the product release mechanism of dengue virus NS3 helicase.

The non-structural protein 3 helicase (NS3h) is a multifunctional protein that is critical in RNA replication and other stages in the flavivirus life cycle. NS3h uses energy from ATP hydrolysis to translocate along single stranded nucleic acid and to unwind double stranded RNA. Here we present a detailed mechanistic analysis of the product release stage in the catalytic cycle of the dengue virus (DENV) NS3h. This study is based on a combined experimental and computational approach of product-inhibition studies and free energy calculations. Our results support a model in which the catalytic cycle of ATP hydrolysis proceeds through an ordered sequential mechanism that includes a ternary complex intermediate (NS3h-Pi-ADP), which evolves releasing the first product, phosphate (Pi), and subsequently ADP. Our results indicate that in the product release stage of the DENV NS3h a novel open-loop conformation plays an important role that may be conserved in NS3 proteins of other flaviviruses as well.

Atomistic insight into the essential binding event of ACE2-derived peptides to the SARS-CoV-2 spike protein.

The pathogenic agent of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters into human cells through the interaction between the receptor binding domain (RBD) of its spike glycoprotein and the angiotensin-converting enzyme 2 (ACE2) receptor. Efforts have been made towards finding antivirals that block this interaction, therefore preventing infection. Here, we determined the binding affinity of ACE2-derived peptides to the RBD of SARS-CoV-2 experimentally and performed MD simulations in order to understand key characteristics of their interaction. One of the peptides, p6, binds to the RBD of SARS-CoV-2 with nM affinity. Although the ACE2-derived peptides retain conformational flexibility when bound to SARS-CoV-2 RBD, we identified residues T27 and K353 as critical anchors mediating the interaction. New ACE2-derived peptides were developed based on the p6-RBD interface analysis and expecting the native conformation of the ACE2 to be maintained. Furthermore, we found a correlation between the helicity in trifluoroethanol and the binding affinity to RBD of the new peptides. Under the hypothesis that the conservation of peptide secondary structure is decisive to the binding affinity, we developed a cyclized version of p6 which had more helicity than p6 and approximately half of its KD value.

Nucleotide-dependent dynamics of the Dengue NS3 helicase.

Dengue represents a substantial public health burden, particularly in low-resource countries. Non-structural protein 3 (NS3) is a multifunctional protein critical in the virus life cycle and has been identified as a promising anti-viral drug target. Despite recent crystallographic studies of the NS3 helicase domain, only subtle structural nucleotide-dependent differences have been identified, such that its coupled ATPase and helicase activities remain mechanistically unclear. Here we use molecular dynamics simulations to explore the nucleotide-dependent conformational landscape of the Dengue virus NS3 helicase and identify substantial changes in the protein flexibility during the ATP hydrolysis cycle. We relate these changes to the RNA-protein interactions and proposed translocation models for other monomeric helicases. Furthermore, we report a novel open-loop conformation with a likely escape route for Pi after hydrolysis, providing new insight into the conformational changes that underlie the ATPase activity of NS3.

Lessons learned about steered molecular dynamics simulations and free energy calculations.

The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.

On the accurate estimation of free energies using the jarzynski equality.

The Jarzynski equality is one of the most widely celebrated and scrutinized nonequilibrium work theorems, relating free energy to the external work performed in nonequilibrium transitions. In practice, the required ensemble average of the Boltzmann weights of infinite nonequilibrium transitions is estimated as a finite sample average, resulting in the so-called Jarzynski estimator, ΔF^J . Alternatively, the second-order approximation of the Jarzynski equality, though seldom invoked, is exact for Gaussian distributions and gives rise to the Fluctuation-Dissipation estimator ΔF^FD . Here we derive the parametric maximum-likelihood estimator (MLE) of the free energy ΔF^ML considering unidirectional work distributions belonging to Gaussian or Gamma families, and compare this estimator to ΔF^J . We further consider bidirectional work distributions belonging to the same families, and compare the corresponding bidirectional ΔF^ML∗ to the Bennett acceptance ratio ( ΔF^BAR ) estimator. We show that, for Gaussian unidirectional work distributions, ΔF^FD is in fact the parametric MLE of the free energy, and as such, the most efficient estimator for this statistical family. We observe that ΔF^ML and ΔF^ML∗ perform better than ΔF^J and ΔF^BAR , for unidirectional and bidirectional distributions, respectively. These results illustrate that the characterization of the underlying work distribution permits an optimal use of the Jarzynski equality. © 2018 Wiley Periodicals, Inc.

Multiscale approach to the activation and phosphotransfer mechanism of CpxA histidine kinase reveals a tight coupling between conformational and chemical steps.

Sensor histidine kinases (SHKs) are an integral component of the molecular machinery that permits bacteria to adapt to widely changing environmental conditions. CpxA, an extensively studied SHK, is a multidomain homodimeric protein with each subunit consisting of a periplasmic sensor domain, a transmembrane domain, a signal-transducing HAMP domain, a dimerization and histidine phospho-acceptor sub-domain (DHp) and a catalytic and ATP-binding subdomain (CA). The key activation event involves the rearrangement of the HAMP-DHp helical core and translation of the CA towards the acceptor histidine, which presumably results in an autokinase-competent complex. In the present work we integrate coarse-grained, all-atom, and hybrid QM-MM computer simulations to probe the large-scale conformational reorganization that takes place from the inactive to the autokinase-competent state (conformational step), and evaluate its relation to the autokinase reaction itself (chemical step). Our results highlight a tight coupling between conformational and chemical steps, underscoring the advantage of CA walking along the DHp core, to favor a reactive tautomeric state of the phospho-acceptor histidine. The results not only represent an example of multiscale modelling, but also show how protein dynamics can promote catalysis.

Distribution of FMR1 and FMR2 Repeats in Argentinean Patients with Primary Ovarian Insufficiency.

The premutation state of FMR1 (Fragile X Mental Retardation 1) has been associated with primary ovarian insufficiency (POI), and is the most common known genetic cause for 46,XX patients. Nevertheless, very few studies have analyzed its frequency in Latin American populations. Additionally, a relationship between alleles carrying a cryptic microdeletion in the 5'UTR of FMR2 and the onset of POI has only been studied in one population. Our aim was to analyze the incidence of FMR1 premutations and putative microdeletions in exon 1 of FMR2 in a cohort of Argentinean women with POI. We studied 133 patients and 84 controls. Fluorescent PCR was performed, and the FMR2 exon 1 was further sequenced in samples presenting less than 11 repeats. We found the frequency of FMR1 premutations to be 6.7% and 2.9% for familial and sporadic patients, respectively. Among controls, 1/84 women presented a premutation. In addition, although we did not find microdeletions in FMR2, we observed a change (T >C) adjacent to the repeats in two sisters with POI. Given the repetitive nature of the sequence involved, we could not ascertain whether this represents a single nucleotide polymorphism (SNP) or a deletion. Therefore, a relationship between FMR2 and POI could not be established for our population.

Structural Insights into the HWE Histidine Kinase Family: The Brucella Blue Light-Activated Histidine Kinase Domain.

In response to light, as part of a two-component system, the Brucella blue light-activated histidine kinase (LOV-HK) increases its autophosphorylation, modulating the virulence of this microorganism. The Brucella histidine kinase (HK) domain belongs to the HWE family, for which there is no structural information. The HWE family is exclusively present in proteobacteria and usually coupled to a wide diversity of light sensor domains. This work reports the crystal structure of the Brucella HK domain, which presents two different dimeric assemblies in the asymmetric unit: one similar to the already described canonical parallel homodimers (C) and the other, an antiparallel non-canonical (NC) dimer, each with distinct relative subdomain orientations and dimerization interfaces. Contrary to these crystallographic structures and unlike other HKs, in solution, the Brucella HK domain is monomeric and still active, showing an astonishing instability of the dimeric interface. Despite this instability, using cross-linking experiments, we show that the C dimer is the functionally relevant species. Mutational analysis demonstrates that the autophosphorylation activity occurs in cis. The different relative subdomain orientations observed for the NC and C states highlight the large conformational flexibility of the HK domain. Through the analysis of these alternative conformations by means of molecular dynamics simulations, we also propose a catalytic mechanism for Brucella LOV-HK.

Accelerated adaptive integration method.

Conformational changes that occur upon ligand binding may be too slow to observe on the time scales routinely accessible using molecular dynamics simulations. The adaptive integration method (AIM) leverages the notion that when a ligand is either fully coupled or decoupled, according to λ, barrier heights may change, making some conformational transitions more accessible at certain λ values. AIM adaptively changes the value of λ in a single simulation so that conformations sampled at one value of λ seed the conformational space sampled at another λ value. Adapting the value of λ throughout a simulation, however, does not resolve issues in sampling when barriers remain high regardless of the λ value. In this work, we introduce a new method, called Accelerated AIM (AcclAIM), in which the potential energy function is flattened at intermediate values of λ, promoting the exploration of conformational space as the ligand is decoupled from its receptor. We show, with both a simple model system (Bromocyclohexane) and the more complex biomolecule Thrombin, that AcclAIM is a promising approach to overcome high barriers in the calculation of free energies, without the need for any statistical reweighting or additional processors.

Inactivating mutation in histone deacetylase 3 stabilizes its active conformation.

Histone deacetylases (HDACs), together with histone acetyltransferases (HATs), regulate gene expression by modulating the acetylation level of chromatin. HDAC3 is implicated in many important cellular processes, particularly in cancer cell proliferation and metastasis, making inhibition of HDAC3 a promising epigenetic treatment for certain cancers. HDAC3 is activated upon complex formation with both inositol tetraphosphate (IP4) and the deacetylase-activating domain (DAD) of multi-protein nuclear receptor corepressor complexes. In previous studies, we have shown that binding of DAD and IP4 to HDAC3 significantly restricts its conformational space towards its stable ternary complex conformation, and suggest this to be the active conformation. Here, we report a single mutation of HDAC3 that is capable of mimicking the stabilizing effects of DAD and IP4, without the presence of either. This mutation, however, results in a total loss of deacetylase activity, prompting a closer evaluation of our understanding of the activation of HDAC3.

w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations.

Free energy governs the equilibrium extent of many biological processes. High barriers separating free energy minima often limit the sampling in molecular dynamics (MD) simulations, leading to inaccurate free energies. Here, we demonstrate enhanced sampling and improved free energy calculations, relative to conventional MD, using windowed accelerated MD within a Hamiltonian replica exchange framework (w-REXAMD). We show that for a case in which multiple conformations are separated by large free energy barriers, w-REXAMD is a useful enhanced sampling technique, without any necessary reweighting.

Hydrophobic effect drives oxygen uptake in myoglobin via histidine E7.

Since the elucidation of the myoglobin (Mb) structure, a histidine residue on the E helix (His-E7) has been proposed to act as a gate with an open or closed conformation controlling access to the active site. Although it is believed that at low pH, the His-E7 gate is in its open conformation, the full relationship between the His-E7 protonation state, its conformation, and ligand migration in Mb is hotly debated. We used molecular dynamics simulations to first address the effect of His-E7 protonation on its conformation. We observed the expected shift from the closed to the open conformation upon protonation, but more importantly, noted a significant difference between the conformations of the two neutral histidine tautomers. We further computed free energy profiles for oxygen migration in each of the possible His-E7 states as well as in two instructive Mb mutants: Ala-E7 and Trp-E7. Our results show that even in the closed conformation, the His-E7 gate does not create a large barrier to oxygen migration and permits oxygen entry with only a small rotation of the imidazole side chain and movement of the E helix. We identify, instead, a hydrophobic site in the E7 channel that can accommodate an apolar diatomic ligand and enhances ligand uptake particularly in the open His-E7 conformation. This rate enhancement is diminished in the closed conformation. Taken together, our results provide a new conceptual framework for the histidine gate hypothesis.

Structural insight into the separate roles of inositol tetraphosphate and deacetylase-activating domain in activation of histone deacetylase 3.

Histone deacetylases (HDACs) repress transcription by deacetylating acetyllysines on specific histone tails. HDAC3 is implicated in neurodegenerative diseases, certain leukemias, and even in disrupting HIV-1 latency. A recent crystal structure of HDAC3 in complex with the deacetylase-activating domain (DAD) of its corepressor complex revealed an inositol tetraphosphate (IP4) molecule at the protein-protein interface. IP4 was shown to play an important, yet mechanistically ambiguous, role in the activity of HDAC3. The goal of this article is to explore the conformational ensemble of HDAC3 in its inactive apo state and in the presence of each or both of DAD and IP4. Using triplicate, 100 ns molecular dynamic simulations, we study the apo, ternary, and intermediate DAD-bound or IP4-bound HDAC3 states. We find that a population-shift effect is induced by the presence of each corepressor, and is most notable in the presence of both. Our results offer new insights into the change in dynamics necessary for the activation of HDAC3 and highlight the roles of IP4 and DAD in this process.

Substrate stereo-specificity in tryptophan dioxygenase and indoleamine 2,3-dioxygenase.

The first and rate-limiting step of the kynurenine pathway, in which tryptophan (Trp) is converted to N-formylkynurenine is catalyzed by two heme-containing proteins, Indoleamine 2,3-dioxygenase (IDO), and Tryptophan 2,3-dioxygenase (TDO). In mammals, TDO is found exclusively in liver tissue, IDO is found ubiquitously in all tissues. IDO has become increasingly popular in pharmaceutical research as it was found to be involved in many physiological situations, including immune escape of cancer. More importantly, small-molecule inhibitors of IDO are currently utilized in cancer therapy. One of the main concerns for the design of human IDO (hIDO) inhibitors is that they should be selective enough to avoid inhibition of TDO. In this work, we have used a combination of classical molecular dynamics (MD) and hybrid quantum-classical (QM/MM) methodologies to establish the structural basis that determine the differences in (a) the interactions of TDO and IDO with small ligands (CO/O(2)) and (b) the substrate stereo-specificity in hIDO and TDO. Our results indicate that the differences in small ligand bound structures of IDO and TDO arise from slight differences in the structure of the bound substrate complex. The results also show that substrate stereo-specificity of TDO is achieved by the perfect fit of L-Trp, but not D-Trp, which exhibits weaker interactions with the protein matrix. For hIDO, the presence of multiple stable binding conformations for L/D-Trp reveal the existence of a large and dynamic active site. Taken together, our data allow determination of key interactions useful for the future design of more potent hIDO-selective inhibitors.