Mehrnoosh Arrar

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.

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/O2) 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. Proteins 2010;

Hydrophobic effect drives oxygen uptake in myoglobin via histidine E7.

Background: The distal histidine in myoglobin is thought to act as a gate regulating oxygen uptake. Results: Neither the open nor closed conformation hinders oxygen uptake; a hydrophobic site is more apparent in the open conformation. Conclusion: The driving force for ligand uptake is the hydrophobic effect. Significance: This amplifies our understanding of the mechanisms through which proteins regulate ligand uptake. 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.

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.

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.