Using ¹⁵N-ammonium to characterise and map potassium binding sites in proteins by NMR spectroscopy

Detecting Bound Ions in Ion Channels by Solid-State NMR Experiments on 15N-Labelled Ammonium Ions

Solid-state NMR allows for the study of membrane proteins under physiological conditions. Here we describe a method for detection of bound ions in the selectivity filter of ion channels using solid-state NMR. This method employs standard 1H-detected solid-state NMR setup and experiment types, which is enabled by using 15N-labelled ammonium ions to mimic potassium ions.

Dissecting the Thermodynamics of ATP Binding to GroEL One Nucleotide at a Time

Variable-temperature electrospray ionization (vT-ESI) native mass spectrometry (nMS) is used to determine the thermodynamics for stepwise binding of up to 14 ATP molecules to the 801 kDa GroEL tetradecamer chaperonin complex. Detailed analysis reveals strong enthalpy-entropy compensation (EEC) for the ATP binding events leading to formation of GroEL-ATP₇ and GroEL-ATP₁₄ complexes. The observed variations in EEC and stepwise free energy changes of specific ATP binding are consistent with the well-established nested cooperativity model describing GroEL-ATP interactions, viz., intraring positive cooperativity and inter-ring negative cooperativity (Dyachenko A.; Proc. Natl. Acad. Sci. U.S.A.2013, 110, 7235-7239). Entropy-driven ATP binding is to be expected for ligand-induced conformational changes of the GroEL tetradecamer, though the magnitude of the entropy change suggests that reorganization of GroEL-hydrating water molecules and/or expulsion of water from the GroEL cavity may also play key roles. The capability for determining complete thermodynamic signatures (ΔG, ΔH, and -TΔS) for individual ligand binding reactions for the large, nearly megadalton GroEL complex expands our fundamental view of chaperonin functional chemistry. Moreover, this work and related studies of protein-ligand interactions illustrate important new capabilities of vT-ESI-nMS for thermodynamic studies of protein interactions with ligands and other molecules such as proteins and drugs.

Tuning the Sensitivity of Genetically Encoded Fluorescent Potassium Indicators through Structure-Guided and Genome Mining Strategies

Genetically encoded potassium indicators lack optimal binding affinity for monitoring intracellular dynamics in mammalian cells. Through structure-guided design and genome mining of potassium binding proteins, we developed green fluorescent potassium indicators with a broad range of binding affinities. KRaION1 (K+ ratiometric indicator for optical imaging based on mNeonGreen 1), based on the insertion of a potassium binding protein, Kbp, from E. coli (Ec-Kbp) into the fluorescent protein mNeonGreen, exhibits an isotonically measured Kd of 69 ± 10 mM (mean ± standard deviation used throughout). We identified Ec-Kbp's binding site using NMR spectroscopy to detect protein-thallium scalar couplings and refined the structure of Ec-Kbp in its potassium-bound state. Guided by this structure, we modified KRaION1, yielding KRaION1/D9N and KRaION2, which exhibit isotonically measured Kd's of 138 ± 21 and 96 ± 9 mM. We identified four Ec-Kbp homologues as potassium binding proteins, which yielded indicators with isotonically measured binding affinities in the 39-112 mM range. KRaIONs functioned in HeLa cells, but the Kd values differed from the isotonically measured case. We found that, by tuning the experimental conditions, Kd values could be obtained that were consistent in vitro and in vivo. We thus recommend characterizing potassium indicator Kd in the physiological context of interest before application.

Direct Detection of Bound Ammonium Ions in the Selectivity Filter of Ion Channels by Solid-State NMR

The flow of ions across cell membranes facilitated by ion channels is an important function for all living cells. Despite the huge amount of structural data provided by crystallography, elucidating the exact interactions between the selectivity filter atoms and bound ions is challenging. Here, we detect bound ¹⁵N-labeled ammonium ions as a mimic for potassium ions in ion channels using solid-state NMR under near-native conditions. The non-selective ion channel NaK showed two ammonium peaks corresponding to its two ion binding sites, while its potassium-selective mutant NaK2K that has a signature potassium-selective selectivity filter with four ion binding sites gave rise to four ammonium peaks. Ions bound in specific ion binding sites were identified based on magnetization transfer between the ions and carbon atoms in the selectivity filters. Magnetization transfer between bound ions and water molecules revealed that only one out of four ions in the selectivity filter of NaK2K is in close contact with water, which is in agreement with the direct knock-on ion conduction mechanism where ions are conducted through the channel by means of direct interactions without water molecules in between. Interestingly, the potassium-selective ion channels investigated here (NaK2K and, additionally, KcsA-Kv1.3) showed remarkably different chemical shifts for their bound ions, despite having identical amino acid sequences and crystal structures of their selectivity filters. Molecular dynamics simulations show similar ion binding and conduction behavior between ammonium and potassium ions and identify the origin of the differences between the investigated potassium channels.

Aromatic side-chain flips orchestrate the conformational sampling of functional loops in human histone deacetylase 8

Human histone deacetylase 8 (HDAC8) is a key hydrolase in gene regulation and an important drug-target. High-resolution structures of HDAC8 in complex with substrates or inhibitors are available, which have provided insights into the bound state of HDAC8 and its function. Here, using long all-atom unbiased molecular dynamics simulations and Markov state modelling, we show a strong correlation between the conformation of aromatic side chains near the active site and opening and closing of the surrounding functional loops of HDAC8. We also investigated two mutants known to allosterically downregulate the enzymatic activity of HDAC8. Based on experimental data, we hypothesise that I19S-HDAC8 is unable to release the product, whereas both product release and substrate binding are impaired in the S39E-HDAC8 mutant. The presented results deliver detailed insights into the functional dynamics of HDAC8 and provide a mechanism for the substantial downregulation caused by allosteric mutations, including a disease causing one.

Nitrogen Assimilation and Transport by Ex Planta Nitrogen-Fixing Bradyrhizobium diazoefficiens Bacteroids Is Modulated by Oxygen, Bacteroid Density and l-Malate

Symbiotic nitrogen fixation requires the transfer of fixed organic nitrogen compounds from the symbiotic bacteria to a host plant, yet the chemical nature of the compounds is in question. Bradyrhizobium diazoefficiens bacteroids were isolated anaerobically from soybean nodules and assayed at varying densities, varying partial pressures of oxygen, and varying levels of l-malate. Ammonium was released at low bacteroid densities and high partial pressures of oxygen, but was apparently taken up at high bacteroid densities and low partial pressures of oxygen in the presence of l-malate; these later conditions were optimal for amino acid excretion. The ratio of partial pressure of oxygen/bacteroid density of apparent ammonium uptake and of alanine excretion displayed an inverse relationship. Ammonium uptake, alanine and branch chain amino acid release were all dependent on the concentration of l-malate displaying similar K0.5 values of 0.5 mM demonstrating concerted regulation. The hyperbolic kinetics of ammonium uptake and amino acid excretion suggests transport via a membrane carrier and also suggested that transport was rate limiting. Glutamate uptake displayed exponential kinetics implying transport via a channel. The chemical nature of the compounds released were dependent upon bacteroid density, partial pressure of oxygen and concentration of l-malate demonstrating an integrated metabolism.

A distal regulatory region of a class I human histone deacetylase

Histone deacetylases (HDACs) are key enzymes in epigenetics and important drug targets in cancer biology. Whilst it has been established that HDACs regulate many cellular processes, far less is known about the regulation of these enzymes themselves. Here, we show that HDAC8 is allosterically regulated by shifts in populations between exchanging states. An inactive state is identified, which is stabilised by a range of mutations and resembles a sparsely-populated state in equilibrium with active HDAC8. Computational models show that the inactive and active states differ by small changes in a regulatory region that extends up to 28 Å from the active site. The regulatory allosteric region identified here in HDAC8 corresponds to regions in other class I HDACs known to bind regulators, thus suggesting a general mechanism. The presented results pave the way for the development of allosteric HDAC inhibitors and regulators to improve the therapy for several disease states.

Enzyme-Substrate-Cofactor Dynamical Networks Revealed by High-Resolution Field Cycling Relaxometry

The remarkable power and specificity of enzyme catalysis rely on the dynamic alignment of the enzyme, substrates, and cofactors, yet the role of dynamics has usually been approached from the perspective of the protein. We have been using an underappreciated NMR technique, subtesla high-resolution field cycling 31P NMR relaxometry, to investigate the dynamics of the enzyme-bound substrates and cofactor on guanosine-5'-monophosphate reductase (GMPR). GMPR forms two dead end, yet catalytically competent, complexes that mimic distinct steps in the catalytic cycle: E·IMP·NADP+ undergoes a partial hydride transfer reaction, while E·GMP·NADP+ undergoes a partial deamination reaction. A different cofactor conformation is required for each partial reaction. Here we report the effects of mutations designed to perturb cofactor conformation and ammonia binding with the goal of identifying the structural features that contribute to the distinct dynamic signatures of the hydride transfer and deamination complexes. These experiments suggest that Asp129 is a central cog in a dynamic network required for both hydride transfer and deamination. In contrast, Lys77 modulates the conformation and mobility of substrates and cofactors in a reaction-specific manner. Thr105 and Tyr318 are part of a deamination-specific dynamic network that includes the 2'-OH of GMP. These residues have comparatively little effect on the dynamic properties of the hydride transfer complex. These results further illustrate the potential of high-resolution field cycling NMR relaxometry for the investigation of ligand dynamics. In addition, exchange experiments indicate that NH3/NH4+ has a high affinity for the deamination complex but a low affinity for the hydride transfer complex, suggesting that the movement of ammonia may gate the cofactor conformational change. Collectively, these experiments reinforce the view that the enzyme, substrates, and cofactor are linked in intricate, reaction-specific, dynamic networks and demonstrate that distal portions of the substrates and cofactors are critical features in these networks.

ATPase activity of GroEL is dependent on GroES and it is response for environmental stress in Riemerella anatipestifer

Riemerella anatipestifer (Ra) is a serious gram-negative pathogen of birds and can cause considerable economic losses. The survival mechanisms of R. anatipestifer in the host and environment remain largely unknown. Previous results have demonstrated that GroEL is a molecular chaperone and an important component of the response to various stresses in most bacteria. This study focused on whether GroEL is implicated in this process in R. anatipestifer. The 1629 bp groEL is highly conserved among other gram-negative bacteria (levels of sequence similarity > 60%). A structural analysis and ATPase activity assay revealed that RaGroEL had weak ATPase activity and that the enzyme activity was temperature and ion dependent. GroES partially enhanced the GroEL ATPase activity in the same temperature range. In addition, we studied the mRNA expression of groEL under abiotic stresses caused by heat shock, pH, salt and hydrogen peroxide. These stresses increased the transcription of groEL to varying degrees. In R. anatipestifer, the ATPase activity of GroEL is dependent on GroES and temperature. The expression of groEL was strongly induced by heat, pH, hydrogen peroxide and salt stress. This study is the first to show that GroEL in R. anatipestifer might play a major role in response to environmental stress.

Computer-Aided Drug Design in Epigenetics

Epigenetic dysfunction has been widely implicated in several diseases especially cancers thus highlights the therapeutic potential for chemical interventions in this field. With rapid development of computational methodologies and high-performance computational resources, computer-aided drug design has emerged as a promising strategy to speed up epigenetic drug discovery. Herein, we make a brief overview of major computational methods reported in the literature including druggability prediction, virtual screening, homology modeling, scaffold hopping, pharmacophore modeling, molecular dynamics simulations, quantum chemistry calculation, and 3D quantitative structure activity relationship that have been successfully applied in the design and discovery of epi-drugs and epi-probes. Finally, we discuss about major limitations of current virtual drug design strategies in epigenetics drug discovery and future directions in this field.

Insights into structural features of HDAC1 and its selectivity inhibition elucidated by Molecular dynamic simulation and Molecular Docking

Histone deacetylases (HDACs) are a family of proteins whose main function is the removal of acetyl groups from lysine residues located on histone and non-histone substrates, which regulates gene transcription and other activities in cells. HDAC1 dysfunction has been implicated in cancer development and progression; thus, its inhibition has emerged as a new therapeutic strategy. Two additional metal binding sites (Site 1 and Site 2) in HDACs have been described that are primarily occupied by potassium ions, suggesting a possible structural role that affects HDAC activity. In this work, we explored the structural role of potassium ions in Site 1 and Site 2 and how they affect the interactions of compounds with high affinities for HDAC1 (AC1OCG0B, Chlamydocin, Dacinostat and Quisinostat) and SAHA (a pan-inhibitor) using molecular docking and molecular dynamics (MD) simulations in concert with a Molecular-Mechanics-Generalized-Born-Surface-Area (MMGBSA) approach. Four models were generated: one with a potassium ion (K⁺) in both sites (HDAC1_k), a second with K⁺ only at site 1 (HDAC1_ks1), a third with K⁺ only at site 2 (HDAC1_ks2) and a fourth with no K⁺ (HDAC1_wk). We found that the presence or absence of K⁺ not only impacted the structural flexibility of HDAC1, but also its molecular recognition, consistent with experimental findings. These results could therefore be useful for further structure-based drug design studies addressing new HDAC1 inhibitors.

Measurement of 15N longitudinal relaxation rates in 15NH4+ spin systems to characterise rotational correlation times and chemical exchange

Many chemical and biological processes rely on the movement of monovalent cations and an understanding of such processes can therefore only be achieved by characterising the dynamics of the involved ions. It has recently been shown that 15N-ammonium can be used as a proxy for potassium to probe potassium binding in bio-molecules such as DNA quadruplexes and enzymes. Moreover, equations have been derived to describe the time-evolution of 15N-based spin density operator elements of 15NH4+ spin systems. Herein NMR pulse sequences are derived to select specific spin density matrix elements of the 15NH4+ spin system and to measure their longitudinal relaxation in order to characterise the rotational correlation time of the 15NH4+ ion as well as report on chemical exchange events of the 15NH4+ ion. Applications to 15NH4+ in acidic aqueous solutions are used to cross-validate the developed pulse sequence while measurements of spin-relaxation rates of 15NH4+ bound to a 41kDa domain of the bacterial Hsp70 homologue DnaK are presented to show the general applicability of the derived pulse sequence. The rotational correlation time obtained for 15N-ammonium bound to DnaK is similar to the correlation time that describes the rotation about the threefold axis of a methyl group. The methodology presented here provides, together with the previous theoretical framework, an important step towards characterising the motional properties of cations in macromolecular systems.

Inhibition mechanism of SAHA in HDAC: a revisit

SAHA (vorinostat, Merck) is a famous clinical drug for zinc-containing histone deacetylase (HDAC) targets against cancer and several other human disorders, whose inhibition mechanism (namely the protonation mechanism) upon binding to HDAC has been debated for more than ten years. It is very challenging to verify experimentally and is still controversial theoretically. The popular "Class-dependent" (namely "Tyr-dependent") hypothesis is that the deprotonation of SAHA is mostly regulated by the conserved Tyr308 in class I HDAC while it is replaced by the His843 in class IIa HDAC. Herein, by elaborate QM(DFT)/MM MD simulations, we exclude the prevalent "Class-dependent" mechanism and advance a novel "Metal-dependent" mechanism, where the remote second metal site (K(+) in most HDAC and Ca(2+) in HDAC2) determines the protonation of SAHA. This proof-of-principle "Metal-dependent" mechanism opens up a new avenue to utilize the second metal site for isoform-selective inhibitor design.

You are lost without a map: Navigating the sea of protein structures

X-ray crystal structures propel biochemistry research like no other experimental method, since they answer many questions directly and inspire new hypotheses. Unfortunately, many users of crystallographic models mistake them for actual experimental data. Crystallographic models are interpretations, several steps removed from the experimental measurements, making it difficult for nonspecialists to assess the quality of the underlying data. Crystallographers mainly rely on "global" measures of data and model quality to build models. Robust validation procedures based on global measures now largely ensure that structures in the Protein Data Bank (PDB) are largely correct. However, global measures do not allow users of crystallographic models to judge the reliability of "local" features in a region of interest. Refinement of a model to fit into an electron density map requires interpretation of the data to produce a single "best" overall model. This process requires inclusion of most probable conformations in areas of poor density. Users who misunderstand this can be misled, especially in regions of the structure that are mobile, including active sites, surface residues, and especially ligands. This article aims to equip users of macromolecular models with tools to critically assess local model quality. Structure users should always check the agreement of the electron density map and the derived model in all areas of interest, even if the global statistics are good. We provide illustrated examples of interpreted electron density as a guide for those unaccustomed to viewing electron density.

Heteronuclear transverse and longitudinal relaxation in AX4 spin systems: application to (15)N relaxations in (15)NH4(+)

The equations that describe the time-evolution of transverse and longitudinal (15)N magnetisations in tetrahedral ammonium ions, (15)NH4(+), are derived from the Bloch-Wangsness-Redfield density operator relaxation theory. It is assumed that the relaxation of the spin-states is dominated by (1) the intra-molecular (15)N-(1)H and (1)H-(1)H dipole-dipole interactions and (2) interactions of the ammonium protons with remote spins, which also include the contribution to the relaxations that arise from the exchange of the ammonium protons with the bulk solvent. The dipole-dipole cross-correlated relaxation mechanisms between each of the (15)N-(1)H and (1)H-(1)H interactions are explicitly taken into account in the derivations. An application to (15)N-ammonium bound to a 41kDa domain of the protein DnaK is presented, where a comparison between experiments and simulations show that the ammonium ion rotates rapidly within its binding site with a local correlation time shorter than approximately 1ns. The theoretical framework provided here forms the basis for further investigations of dynamics of AX4 spin systems, with ammonium ions in solution and bound to proteins of particular interest.

Inhibition and mechanism of HDAC8 revisited

Histone deacetylases (HDACs) have found intense interest as drug targets for a variety of diseases, but there is disagreement about basic aspects of the inhibition and mechanism of HDACs. QM/MM calculations of HDAC8 including a large QM region provide a model that is consistent with the available crystal structures and structure-activity relationships of different HDAC inhibitors. The calculations support a spontaneous proton transfer from a hydroxamic acid to an active site histidine upon binding to the zinc. The role of the H142/D176 catalytic dyad as the general base of the reaction is elucidated. The reasons for the disagreements between previous proposals are discussed. The results provide detailed insights into the unique mechanism of HDACs, including the role of the two catalytic dyads and function of the potassium near the active site. They also have important implications for the design of novel inhibitors for a number of HDACs such as the class IIa HDACs.