The closed conformation of a highly flexible protein: the structure of E. coli adenylate kinase with bound AMP and AMPPNP

Nucleotide binding to the ATP-cone in anaerobic ribonucleotide reductases allosterically regulates activity by modulating substrate binding

A small, nucleotide-binding domain, the ATP-cone, is found at the N-terminus of most ribonucleotide reductase (RNR) catalytic subunits. By binding adenosine triphosphate (ATP) or deoxyadenosine triphosphate (dATP) it regulates the enzyme activity of all classes of RNR. Functional and structural work on aerobic RNRs has revealed a plethora of ways in which dATP inhibits activity by inducing oligomerisation and preventing a productive radical transfer from one subunit to the active site in the other. Anaerobic RNRs, on the other hand, store a stable glycyl radical next to the active site and the basis for their dATP-dependent inhibition is completely unknown. We present biochemical, biophysical, and structural information on the effects of ATP and dATP binding to the anaerobic RNR from Prevotella copri. The enzyme exists in a dimer-tetramer equilibrium biased towards dimers when two ATP molecules are bound to the ATP-cone and tetramers when two dATP molecules are bound. In the presence of ATP, P. copri NrdD is active and has a fully ordered glycyl radical domain (GRD) in one monomer of the dimer. Binding of dATP to the ATP-cone results in loss of activity and increased dynamics of the GRD, such that it cannot be detected in the cryo-EM structures. The glycyl radical is formed even in the dATP-bound form, but the substrate does not bind. The structures implicate a complex network of interactions in activity regulation that involve the GRD more than 30 Å away from the dATP molecules, the allosteric substrate specificity site and a conserved but previously unseen flap over the active site. Taken together, the results suggest that dATP inhibition in anaerobic RNRs acts by increasing the flexibility of the flap and GRD, thereby preventing both substrate binding and radical mobilisation.

AK7-deficiency reversal inhibits ccRCC progression and boosts anti-PD1 immunotherapy sensitivity

Clear cell renal cell carcinoma (ccRCC) is a common kidney cancer with subtle early symptoms, high recurrence rates, and low sensitivity to traditional treatments like radiotherapy and chemotherapy. Identifying novel therapeutic targets is critical. The expression level of adenylate kinases 7 (AK7) in ccRCC was examined by the TCGAportal and UALCAN databases. The effect of AK7 on proliferation, invasion and migration of ccRCC cell lines was evaluated by cell assay. The correlation between AK7 expression and prognosis, as well as its direct relationship with immunotherapy efficacy, was analyzed using CANCERTOOL and Kaplan-Meier plotter data. Moreover, the TISIDB database was used to study the relationship between AK7 and immune markers. The effect of overexpressed AK7 combined with PD1 monoclonal antibody on ccRCC was evaluated in animal experiments. The results showed that low level of AK7 expression was observed in ccRCC tissues. The expression of AK7 can regulate the proliferation, invasion, and migration of human ccRCC cell lines. The level of AK7 expression was associated with OS of ccRCC patients. This was potentially due to the negative connection between AK7 expression and CD8+ T cell depletion, indicating that immunotherapy might be less effective in individuals with low AK7 expression. Conversely, augmenting AK7 demonstrated an enhanced effectiveness of anti-PD1 therapy. The findings of our research strongly indicated that AK7 could serve as both a prognostic indicator and therapeutic target for patients with ccRCC. Moreover, the overexpression of AK7 combined with anti-PD1 held promising potential as a therapeutic approach for treating ccRCC.

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Four-dimensional imaging for cryo-electron microscopy experiments using molecular simulations and manifold learning

Elucidating protein conformational changes is essential because conformational changes are closely related to the functions of proteins. Cryo-electron microscopy (cryo-EM) experiment can be used to reconstruct protein conformational changes via a method that involves using the experimental data (two-dimensional protein images). In this study, a reconstruction method, referred to as the "four-dimensional imaging," was proposed. In our four-dimensional imaging technique, the protein conformational change was obtained using the two-dimensional protein images (the three-dimensional electron density maps used in previously proposed techniques were not used). The protein conformation for each two-dimensional protein image was obtained using our original protocol with molecular dynamics simulations. Using a manifold-learning technique and two-dimensional protein images, the protein conformations were arranged according to the conformational change of the protein. By arranging the protein conformations according to the arrangement of the protein images, four-dimensional imaging is constructed. A simulation for a cryo-EM experiment demonstrated the validity of our four-dimensional imaging technique.

Switching Go̅-Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions

Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and β2-adrenergic receptor (β2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Conformational dynamics of adenylate kinase in crystals

Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.

Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase

Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.

Accurate and Efficient Estimation of Lennard-Jones Interactions for Coarse-Grained Particles via a Potential Matching Method

The Lennard-Jones (LJ) energy functions are commonly used to describe the nonbonded interactions in bulk coarse-grained (CG) models, which contribute significantly to the stabilization of a local binding configuration or a self-assembly system. In many cases, systematic development of the LJ interaction parameters in a CG model requires a comprehensive sampling of the objective molecules at the all-atom (AA) level, which is therefore extremely time-consuming for large systems. Inspired by the concept of electrostatic potential (ESP), we define the LJ static potential (LJSP), by which the embedding potential energy surface can be constructed analytically. A semianalytic approach, namely, the LJSP matching method, is developed here to derive the CG parameters by minimizing the LJSP difference between the AA and the CG models, which provides a universal way to derive the CG LJ parameters from the AA models without doing presampling. The LJSP matching method is successful not only in deriving the LJ interaction energy landscape in the CG models for proteins, lipids, and DNA but also in reproducing the critical properties such as intermediate structures and enthalpy contributions as exemplified in simulating the self-assembly process of the dipalmitoylphosphatidylcholine (DPPC) lipids.

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Effects of antibiotics on enhanced biological phosphorus removal and its mechanisms

Many kinds of antibiotics are continuously discharged into wastewater and typically cause a great decrease in sewage treatment performance, whereas mechanisms of differences in the impacts of commonly used antibiotics on phosphate removal are still elusive. Thus, an enhanced biological phosphorus removal (EBPR) system, as an effective method of phosphate removal, was developed, and its performance in the treatment of artificial wastewater containing antibiotics at short- (8 h) and long-term (15 days) exposure was investigated. The results show that phosphorus removal was consistently inhibited by the addition of antibiotics with a significant difference (P < 0.05). To interpret the phenomena, mechanistic equations were developed, and the results indicate that for short-term tests, the difference was mainly caused by the suppression of polyhydroxyalkanoate (PHA) degradation and the activity of polyphosphate kinase (PPK), resulting in the different inhibition of the soluble orthophosphorus (SOP) uptake process. For long-term tests, the difference in SOP uptake was principally caused by the inhibition of PHA degradation and the activity of PPK, whereas the difference in SOP release resulted from the inhibition of activities of exopolyphosphatase (PPX) and adenylate kinase (ADK). Moreover, micro-mechanisms of such inhibition were identified from molecular docking and electrostatic potential.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Double-Well Ultra-Coarse-Grained Model to Describe Protein Conformational Transitions

The double-well model is usually used to describe the conformational transition between two states of a protein. Since conformational changes usually occur within a relatively large time scale, coarse-grained models are often used to accelerate the dynamic process due to their inexpensive computational cost. In this work, we develop a double-well ultra-coarse-grained (DW-UCG) model to describe the conformational transitions of the adenylate kinase, glutamine-binding protein, and lactoferrin. The coarse-grained simulation results show that the DW-UCG model of adenylate kinase captures the crucial intermediate states in the LID-closing and NMP-closing pathways, reflecting the key secondary structural changes in the conformational transition. A comparison of the different DW-UCG models of adenylate kinase indicates that an appropriate choice of bead resolution could generate the free energy landscape that is comparable to that from the residue-based model. The coarse-grained simulations for the glutamine-binding protein and lactoferrin also demonstrate that the DW-UCG model is valid in reproducing the correct two-state behavior for their functional study, which indicates the potential application of the DW-UCG model in investigating the mechanism of conformational changes of large proteins.

Structural Basis for GTP versus ATP Selectivity in the NMP Kinase AK3

ATP and GTP are exceptionally important molecules in biology with multiple, and often discrete, functions. Therefore, enzymes that bind to either of them must develop robust mechanisms to selectively utilize one or the other. Here, this specific problem is addressed by molecular studies of the human NMP kinase AK3, which uses GTP to phosphorylate AMP. AK3 plays an important role in the citric acid cycle, where it is responsible for GTP/GDP recycling. By combining a structural biology approach with functional experiments, we present a comprehensive structural and mechanistic understanding of the enzyme. We discovered that AK3 functions by recruitment of GTP to the active site, while ATP is rejected and nonproductively bound to the AMP binding site. Consequently, ATP acts as an inhibitor with respect to GTP and AMP. The overall features with specific recognition of the correct substrate and nonproductive binding by the incorrect substrate bear a strong similarity to previous findings for the ATP specific NMP kinase adenylate kinase. Taken together, we are now able to provide the fundamental principles for GTP and ATP selectivity in the large NMP kinase family. As a side-result originating from nonlinearity of chemical shifts in GTP and ATP titrations, we find that protein surfaces offer a general and weak binding affinity for both GTP and ATP. These nonspecific interactions likely act to lower the available intracellular GTP and ATP concentrations and may have driven evolution of the Michaelis constants of NMP kinases accordingly.

Molecular docking investigation of the amantadine binding to the enzymes upregulated or downregulated in Parkinson's disease

Parkinson's disease (PD) is a progressive neurodegenerative disease. Levodopa in combination with amantadine has a demonstrated efficacy in motility impairment. An extensive investigation of some enzymes described to be upregulated or downregulated in PD was made - adenylate kinase (AK), adenine phosphoribosyltransferase (APRT), ectonucleoside triphosphate diphosphohydrolase 1 (ENTPD1), nucleoside-diphosphate kinase 3 (NDK3), purine nucleoside phosphorylase 1 (PNP1), and ecto-5'-nucleotidase (NT5E). Also, creatine kinase (CK) was included in the study because it is one of the main enzymes involved in the regulation of the nucleotide ratio in energy metabolism. To date, there is no proven link between amantadine treatment of PD and these enzymes. Because there are many AKs isoforms modified in PD, the AK was the first investigated. The molecular docking experiments allow the analysis of the selective binding of amantadine - unionized (with -NH₂ group) and ionized form (with -NH₃ ⁺ group) - to the AKs' isoforms implicated in PD. Using available X-ray 3D structures of human AKs in closed-conformation, it was demonstrated that there are notable differences between the interactions of the two forms of amantadine for the zebrafish AK1 (5XZ2), human AK2 (2C9Y), human AK5 (2BWJ), and AK from B.stearothermophilus. The cytosolic human AK1 and human AK2 mostly interact with ionized amantadine by AMP binding residues. The human AK5 interaction with ionized amantadine does not involve the residues from the catalytic site. Among other enzymes tested in the present study, APRT revealed the best results in respect of binding amantadine ionized form. The results offer a new perspective for further investigation of the connections between amantadine treatment of PD and some enzymes involved in purine metabolism.

Optical Control of Cellular ATP Levels with a Photocaged Adenylate Kinase

We have developed a new tool for the optical control of cellular ATP concentrations with a photocaged adenylate kinase (Adk). The photocaged Adk is generated by substituting a catalytically essential lysine with a hydroxycoumarin-protected lysine through site-specific unnatural amino acid mutagenesis in both E. coli and mammalian cells. Caging of the critical lysine residue offers complete suppression of Adk's phosphotransferase activity and rapid restoration of its function both in vitro and in vivo upon optical stimulation. Light-activated Adk renders faster rescue of cell growth than chemically inducible expression of wild-type Adk in E. coli as well as rapid ATP depletion in mammalian cells. Thus, caging Adk provides a new tool for direct conditional perturbation of cellular ATP concentrations thereby enabling the investigation of ATP-coupled physiological events in temporally dynamic contexts.

Adenylate Kinase: A Ubiquitous Enzyme Correlated with Medical Conditions

Adenylate kinase is a small, usually monomeric, enzyme found in every living thing due to its crucial role in energetic metabolism. This paper outlines the most relevant data about adenylate kinases isoforms, and the connection between dysregulation or mutation of human adenylate kinase and medical conditions. The following datadases were consulted: National Centre for Biotechnology Information, Protein Data Bank, and Mouse Genomic Informatics. The SmartBLAST tool, EMBOSS Needle Program, and Clustal Omega Program were used to analyze the best protein match, and to perform pairwise sequence alignment and multiple sequence alignment. Human adenylate kinase genes are located on different chromosomes, six of them being on the chromosomes 1 and 9. The adenylate kinases' intracellular localization and organ distribution explain their dysregulation in many diseases. The cytosolic isoenzyme 1 and the mitochondrial isoenzyme 2 are the main adenylate kinases that are integrated in the vast network of inflammatory modulators. The cytosolic isoenzyme 5 is correlated with limbic encephalitis and Leu673Pro mutation of the isoenzyme 7 leads to primary male infertility due to impairment of the ciliary function. The impairment of the mitochondrial isoenzymes 2 and 4 is demonstrated in neuroblastoma or glioma. The adenylate kinases are disease modifier that can assess the risk of diseases where oxidative stress plays a crucial role in pathogenesis like metabolic syndrome or neurodegenerative diseases. Because adenylate kinases has ATP as substrate, they are integrated in the global network of energetic process of any organism therefore are valid target for new pharmaceutical compounds.

Determining Optimal Coarse-Grained Representation for Biomolecules Using Internal Cluster Validation Indexes

The development of ultracoarse-grained models for large biomolecules needs to derive the optimal number of coarse-grained (CG) sites to represent the targets. In this work, we propose to use the statistical internal cluster validation indexes to determine the optimal number of CG sites that are optimized based on the essential dynamics coarse-graining method. The calculated curves of Calinski-Harabasz and Silhouette Coefficient indexes exhibit the extrema corresponding to the similar CG numbers. The calculated ratios of the optimal CG numbers to the residue numbers of fine-grained models are in the range from 4 to 2. The comparison of the stability of index results indicates that Calinski-Harabasz index is the better choice to determine the optimal CG representation in coarse-graining. © 2019 Wiley Periodicals, Inc.

Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase

Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.

Quantifying the Intrinsic Conformation Energy Landscape Topography of Proteins with Large-Scale Open-Closed Transition

Large-scale conformational changes of proteins, including the open-closed transitions, are crucial for a variety of protein functions. These open-closed transitions are often associated with ligand binding. However, the understandings of the underlying mechanisms of the conformational changes within proteins during the open-closed transitions are still challenging at present. In this study, we quantified the intrinsic underlying conformational energy landscapes of five different proteins with large-scale open-closed transitions. This is realized by exploring the underlying density of states and the intrinsic conformational energy landscape topography measure Λ. Λ is a dimensionless ratio of conformational energy gap δE versus conformational energy roughness δE and configurational entropy S or size of the intrinsic conformational energy landscape. By quantifying the Λ of intrinsic open-closed conformational (Λ^oc) and intrinsic global folding (Λ^global) energy landscapes, we show that both intrinsic open-closed conformation energy and entropy landscapes are funneled toward the closed state. Furthermore, our results indicate the strong correlations between Λ and thermodynamics (conformational state transition temperature against trapping temperature) as well as between Λ and kinetics (open-closed kinetic time) of these proteins. This shows that the intrinsic conformational landscape topography determines both the conformational thermodynamic stability and kinetic speed of the conformational dynamics. Our investigations provide important insights for understanding the fundamental mechanisms of the protein conformational dynamics in a physical and global way.

Molecular Docking Evaluation of (E)-5-arylidene-2-thioxothiazolidin-4-one Derivatives as Selective Bacterial Adenylate Kinase Inhibitors

Multi-drug resistant microorganism infections with emerging problems that require not only a prevention strategy, but also the development of new inhibitory compounds. Six previously synthesized 5-arylidene-2-thioxothiazolidin-4-one derivatives 1a–f, were screened for inhibitory activity on adenylate kinases of different origins by molecular docking. The compounds 1c and 1d were the most efficient inhibitors of bacterial and some archean adenylate kinases. Hydrogen bond interactions were observed with the residues belonging to the ATP binding site. Moreover human adenylate kinases are poor targets, suggesting that this selectivity offers promising prospectives for refining the structure of our compounds.

Molecular mechanism of ATP versus GTP selectivity of adenylate kinase

Enzymatic substrate selectivity is critical for the precise control of metabolic pathways. In cases where chemically related substrates are present inside cells, robust mechanisms of substrate selectivity are required. Here, we report the mechanism utilized for catalytic ATP versus GTP selectivity during adenylate kinase (Adk) -mediated phosphorylation of AMP. Using NMR spectroscopy we found that while Adk adopts a catalytically competent and closed structural state in complex with ATP, the enzyme is arrested in a catalytically inhibited and open state in complex with GTP. X-ray crystallography experiments revealed that the interaction interfaces supporting ATP and GTP recognition, in part, are mediated by coinciding residues. The mechanism provides an atomic view on how the cellular GTP pool is protected from Adk turnover, which is important because GTP has many specialized cellular functions. In further support of this mechanism, a structure-function analysis enabled by synthesis of ATP analogs suggests that a hydrogen bond between the adenine moiety and the backbone of the enzyme is vital for ATP selectivity. The importance of the hydrogen bond for substrate selectivity is likely general given the conservation of its location and orientation across the family of eukaryotic protein kinases.

Effect of ligand binding on a protein with a complex folding landscape

Ligand binding to a protein can stabilize it significantly against unfolding. The variation of the folding free energy, ΔΔG⁰, due to ligand binding can be derived from a simple reaction scheme involving exclusive binding to the native state. One obtains the following expression: , where K_d is the ligand dissociation constant and L is its concentration, R is the universal gas constant and T is the temperature. This expression has been shown to correctly describe experimental results on multiple proteins. In the current work we studied the effect of ligand binding on the stability of the multi-domain protein adenylate kinase from E. coli (AKE). Unfolding experiments were conducted using single-molecule FRET spectroscopy, which allowed us to directly obtain the fraction of unfolded protein in a model-free way from FRET efficiency histograms. Surprisingly, it was found that the effect of two inhibitors (Ap₅A and AMPPNP) and a substrate (AMP) on the stability of AKE was much smaller than expected based on K_d values obtained independently using microscale thermophoresis. To shed light on this issue, we measured the K_d for Ap₅A over a range of chemical denaturant concentrations where the protein is still folded. It was found that K_d increases dramatically over this range, likely due to the population of folding intermediates, whose binding to the ligand is much weaker than that of the native state. We propose that binding to folding intermediates may dominate the effect of ligands on the stability of multi-domain proteins, and could therefore have a strong impact on protein homeostasis in vivo.

Effects of Catalytic Action and Ligand Binding on Conformational Ensembles of Adenylate Kinase

Crystal structures of adenylate kinase (AdK) from Escherichia coli capture two states: an "open" conformation (apo) obtained in the absence of ligands and a "closed" conformation in which ligands are bound. Other AdK crystal structures suggest intermediate conformations that may lie on the transition pathway between these two states. To characterize the transition from open to closed states in solution, X-ray solution scattering data were collected from AdK in the apo form and with progressively increasing concentrations of five different ligands. Scattering data from apo AdK are consistent with scattering predicted from the crystal structure of AdK in the open conformation. In contrast, data from AdK samples saturated with Ap5A do not agree with that calculated from AdK in the closed conformation. Using cluster analysis of available structures, we selected representative structures in five conformational states: open, partially open, intermediate, partially closed, and closed. We used these structures to estimate the relative abundances of these states for each experimental condition. X-ray solution scattering data obtained from AdK with AMP are dominated by scattering from AdK in the open conformation. For AdK in the presence of high concentrations of ATP and ADP, the conformational ensemble shifts to a mixture of partially open and closed states. Even when AdK is saturated with Ap5A, a significant proportion of AdK remains in a partially open conformation. These results are consistent with an induced-fit model in which the transition of AdK from an open state to a closed state is initiated by ATP binding.

Structural basis for ligand binding to an enzyme by a conformational selection pathway

Proteins can bind target molecules through either induced fit or conformational selection pathways. In the conformational selection model, a protein samples a scarcely populated high-energy state that resembles a target-bound conformation. In enzymatic catalysis, such high-energy states have been identified as crucial entities for activity and the dynamic interconversion between ground states and high-energy states can constitute the rate-limiting step for catalytic turnover. The transient nature of these states has precluded direct observation of their properties. Here, we present a molecular description of a high-energy enzyme state in a conformational selection pathway by an experimental strategy centered on NMR spectroscopy, protein engineering, and X-ray crystallography. Through the introduction of a disulfide bond, we succeeded in arresting the enzyme adenylate kinase in a closed high-energy conformation that is on-pathway for catalysis. A 1.9-Å X-ray structure of the arrested enzyme in complex with a transition state analog shows that catalytic sidechains are properly aligned for catalysis. We discovered that the structural sampling of the substrate free enzyme corresponds to the complete amplitude that is associated with formation of the closed and catalytically active state. In addition, we found that the trapped high-energy state displayed improved ligand binding affinity, compared with the wild-type enzyme, demonstrating that substrate binding to the high-energy state is not occluded by steric hindrance. Finally, we show that quenching of fast time scale motions observed upon ligand binding to adenylate kinase is dominated by enzyme-substrate interactions and not by intramolecular interactions resulting from the conformational change.

Two-bead polarizable water models combined with a two-bead multipole force field (TMFF) for coarse-grained simulation of proteins

(a) Four water molecules contained in the polarizable CG water models in (b) two-bead polarizable water model 1 (TPW1) and (c) two-bead polarizable water model 2 (TPW2).

Phosphoribosyl Diphosphate (PRPP): Biosynthesis, Enzymology, Utilization, and Metabolic Significance

Phosphoribosyl diphosphate (PRPP) is an important intermediate in cellular metabolism. PRPP is synthesized by PRPP synthase, as follows: ribose 5-phosphate + ATP → PRPP + AMP. PRPP is ubiquitously found in living organisms and is used in substitution reactions with the formation of glycosidic bonds. PRPP is utilized in the biosynthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, the cofactors NAD and tetrahydromethanopterin, arabinosyl monophosphodecaprenol, and certain aminoglycoside antibiotics. The participation of PRPP in each of these metabolic pathways is reviewed. Central to the metabolism of PRPP is PRPP synthase, which has been studied from all kingdoms of life by classical mechanistic procedures. The results of these analyses are unified with recent progress in molecular enzymology and the elucidation of the three-dimensional structures of PRPP synthases from eubacteria, archaea, and humans. The structures and mechanisms of catalysis of the five diphosphoryltransferases are compared, as are those of selected enzymes of diphosphoryl transfer, phosphoryl transfer, and nucleotidyl transfer reactions. PRPP is used as a substrate by a large number phosphoribosyltransferases. The protein structures and reaction mechanisms of these phosphoribosyltransferases vary and demonstrate the versatility of PRPP as an intermediate in cellular physiology. PRPP synthases appear to have originated from a phosphoribosyltransferase during evolution, as demonstrated by phylogenetic analysis. PRPP, furthermore, is an effector molecule of purine and pyrimidine nucleotide biosynthesis, either by binding to PurR or PyrR regulatory proteins or as an allosteric activator of carbamoylphosphate synthetase. Genetic analyses have disclosed a number of mutants altered in the PRPP synthase-specifying genes in humans as well as bacterial species.

Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case

Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.

Subnanometre enzyme mechanics probed by single-molecule force spectroscopy

Enzymes are molecular machines that bind substrates specifically, provide an adequate chemical environment for catalysis and exchange products rapidly, to ensure fast turnover rates. Direct information about the energetics that drive conformational changes is difficult to obtain. We used subnanometre single-molecule force spectroscopy to study the energetic drive of substrate-dependent lid closing in the enzyme adenylate kinase. Here we show that in the presence of the bisubstrate inhibitor diadenosine pentaphosphate (AP5A), closing and opening of both lids is cooperative and tightly coupled to inhibitor binding. Surprisingly, binding of the substrates ADP and ATP exhibits a much smaller energetic drive towards the fully closed state. Instead, we observe a new dominant energetic minimum with both lids half closed. Our results, combining experiment and molecular dynamics simulations, give detailed mechanical insights into how an enzyme can cope with the seemingly contradictory requirements of rapid substrate exchange and tight closing, to ensure efficient catalysis.

Adenylate kinase catalyses the interconversion of adenosine phosphates, and plays a crucial role in maintaining cellular energy homeostasis. Here, the authors use single molecule optical tweezers to understand how the enzyme's conformation dynamics modulates catalysis.

A new algorithm for construction of coarse-grained sites of large biomolecules

The development of coarse-grained (CG) models for large biomolecules remains a challenge in multiscale simulations, including a rigorous definition of CG representations for them. In this work, we proposed a new stepwise optimization imposed with the boundary-constraint (SOBC) algorithm to construct the CG sites of large biomolecules, based on the s cheme of essential dynamics CG. By means of SOBC, we can rigorously derive the CG representations of biomolecules with less computational cost. The SOBC is particularly efficient for the CG definition of large systems with thousands of residues. The resulted CG sites can be parameterized as a CG model using the normal mode analysis based fluctuation matching method. Through normal mode analysis, the obtained modes of CG model can accurately reflect the functionally related slow motions of biomolecules. The SOBC algorithm can be used for the construction of CG sites of large biomolecules such as F-actin and for the study of mechanical properties of biomaterials.

Nucleotides, Nucleosides, and Nucleobases

We review literature on the metabolism of ribo- and deoxyribonucleotides, nucleosides, and nucleobases in Escherichia coli and Salmonella,including biosynthesis, degradation, interconversion, and transport. Emphasis is placed on enzymology and regulation of the pathways, at both the level of gene expression and the control of enzyme activity. The paper begins with an overview of the reactions that form and break the N-glycosyl bond, which binds the nucleobase to the ribosyl moiety in nucleotides and nucleosides, and the enzymes involved in the interconversion of the different phosphorylated states of the nucleotides. Next, the de novo pathways for purine and pyrimidine nucleotide biosynthesis are discussed in detail.Finally, the conversion of nucleosides and nucleobases to nucleotides, i.e.,the salvage reactions, are described. The formation of deoxyribonucleotides is discussed, with emphasis on ribonucleotidereductase and pathways involved in fomation of dUMP. At the end, we discuss transport systems for nucleosides and nucleobases and also pathways for breakdown of the nucleobases.

Classification of projection images of proteins with structural polymorphism by manifold: a simulation study for x-ray free-electron laser diffraction imaging

Coherent x-ray diffraction imaging (CXDI) enables us to visualize noncrystalline sample particles with micrometer to submicrometer dimensions. Using x-ray free-electron laser (XFEL) sources, two-dimensional diffraction patterns are collected from fresh samples supplied to the irradiation area in the "diffraction-before-destruction" scheme. A recent significant increase in the intensity of the XFEL pulse is promising and will allow us to visualize the three-dimensional structures of proteins using XFEL-CXDI in the future. For the protocol proposed for molecular structure determination using future XFEL-CXDI [T. Oroguchi and M. Nakasako, Phys. Rev. E 87, 022712 (2013)10.1103/PhysRevE.87.022712], we require an algorithm that can classify the data in accordance with the structural polymorphism of proteins arising from their conformational dynamics. However, most of the algorithms proposed primarily require the numbers of conformational classes, and then the results are biased by the numbers. To improve this point, here we examine whether a method based on the manifold concept can classify simulated XFEL-CXDI data with respect to the structural polymorphism of a protein that predominantly adopts two states. After random sampling of the conformations of the two states and in-between states from the trajectories of molecular dynamics simulations, a diffraction pattern is calculated from each conformation. Classification was performed by using our custom-made program suite named enma, in which the diffusion map (DM) method developed based on the manifold concept was implemented. We successfully classify most of the projection electron density maps phase retrieved from diffraction patterns into each of the two states and in-between conformations without the knowledge of the number of conformational classes. We also examined the classification of the projection electron density maps of each of the three states with respect to the Euler angle. The present results suggest that the DM method is imperative for future applications of XFEL-CXDI experiments for proteins, and clarify issues to be taken care of in the future application.

Mutating the Conserved Q-loop Glutamine 1291 Selectively Disrupts Adenylate Kinase-dependent Channel Gating of the ATP-binding Cassette (ABC) Adenylate Kinase Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) and Reduces Channel Function in Primary Human Airway Epithelia

The ATP-binding cassette (ABC) transporter cystic fibrosis transmembrane conductance regulator (CFTR) and two other non-membrane-bound ABC proteins, Rad50 and a structural maintenance of chromosome (SMC) protein, exhibit adenylate kinase activity in the presence of physiologic concentrations of ATP and AMP or ADP (ATP + AMP ⇆ 2 ADP). The crystal structure of the nucleotide-binding domain of an SMC protein in complex with the adenylate kinase bisubstrate inhibitor P(1),P(5)-di(adenosine-5') pentaphosphate (Ap5A) suggests that AMP binds to the conserved Q-loop glutamine during the adenylate kinase reaction. Therefore, we hypothesized that mutating the corresponding residue in CFTR, Gln-1291, selectively disrupts adenylate kinase-dependent channel gating at physiologic nucleotide concentrations. We found that substituting Gln-1291 with bulky side-chain amino acids abolished the effects of Ap5A, AMP, and adenosine 5'-monophosphoramidate on CFTR channel function. 8-Azidoadenosine 5'-monophosphate photolabeling of the AMP-binding site and adenylate kinase activity were disrupted in Q1291F CFTR. The Gln-1291 mutations did not alter the potency of ATP at stimulating current or ATP-dependent gating when ATP was the only nucleotide present. However, when physiologic concentrations of ADP and AMP were added, adenylate kinase-deficient Q1291F channels opened significantly less than wild type. Consistent with this result, we found that Q1291F CFTR displayed significantly reduced Cl(-) channel function in well differentiated primary human airway epithelia. These results indicate that a highly conserved residue of an ABC transporter plays an important role in adenylate kinase-dependent CFTR gating. Furthermore, the results suggest that adenylate kinase activity is important for normal CFTR channel function in airway epithelia.

The energy landscape of adenylate kinase during catalysis

Kinases perform phosphoryl-transfer reactions in milliseconds; without enzymes, these reactions would take about 8,000 years under physiological conditions. Despite extensive studies, a comprehensive understanding of kinase energy landscapes, including both chemical and conformational steps, is lacking. Here we scrutinize the microscopic steps in the catalytic cycle of adenylate kinase, through a combination of NMR measurements during catalysis, pre-steady-state kinetics, molecular-dynamics simulations and crystallography of active complexes. We find that the Mg(2+) cofactor activates two distinct molecular events: phosphoryl transfer (>10(5)-fold) and lid opening (10(3)-fold). In contrast, mutation of an essential active site arginine decelerates phosphoryl transfer 10(3)-fold without substantially affecting lid opening. Our results highlight the importance of the entire energy landscape in catalysis and suggest that adenylate kinases have evolved to activate key processes simultaneously by precise placement of a single, charged and very abundant cofactor in a preorganized active site.

In the multi-domain protein adenylate kinase, domain insertion facilitates cooperative folding while accommodating function at domain interfaces

Having multiple domains in proteins can lead to partial folding and increased aggregation. Folding cooperativity, the all or nothing folding of a protein, can reduce this aggregation propensity. In agreement with bulk experiments, a coarse-grained structure-based model of the three-domain protein, E. coli Adenylate kinase (AKE), folds cooperatively. Domain interfaces have previously been implicated in the cooperative folding of multi-domain proteins. To understand their role in AKE folding, we computationally create mutants with deleted inter-domain interfaces and simulate their folding. We find that inter-domain interfaces play a minor role in the folding cooperativity of AKE. On further analysis, we find that unlike other multi-domain proteins whose folding has been studied, the domains of AKE are not singly-linked. Two of its domains have two linkers to the third one, i.e., they are inserted into the third one. We use circular permutation to modify AKE chain-connectivity and convert inserted-domains into singly-linked domains. We find that domain insertion in AKE achieves the following: (1) It facilitates folding cooperativity even when domains have different stabilities. Insertion constrains the N- and C-termini of inserted domains and stabilizes their folded states. Therefore, domains that perform conformational transitions can be smaller with fewer stabilizing interactions. (2) Inter-domain interactions are not needed to promote folding cooperativity and can be tuned for function. In AKE, these interactions help promote conformational dynamics limited catalysis. Finally, using structural bioinformatics, we suggest that domain insertion may also facilitate the cooperative folding of other multi-domain proteins.

Most individual protein domains fold in an all or nothing fashion. This cooperative folding is important because it reduces the existence of partially folded proteins which can stick to each other and create disease causing aggregates. However, numerous proteins have multiple domains, independent units of folding, stability and/or function. Several such proteins also fold cooperatively. It is thought that strong interactions between individual domains allow the folding to propagate from a nucleating domain to neighbouring ones and this enables cooperative folding in multi-domain proteins. Here, we computationally study the folding of the three-domain protein AKE and find instead that the topology of the protein, wherein the two less stable domains are inserted into the more stable one, promotes folding cooperativity. When the more stable domain is folded, the ends of the inserted domains are constrained and this allows them to fold easily. In such a protein topology, strong inter-domain interactions are not needed to promote folding cooperativity. Interface amino acids which would have been involved in ensuring that the domains fit together correctly can now be tuned for binding or catalysis or conformational transitions. Thus, inserted domains may be present in multi-domain proteins to promote both function and folding.

ATP and AMP mutually influence their interaction with the ATP-binding cassette (ABC) adenylate kinase cystic fibrosis transmembrane conductance regulator (CFTR) at separate binding sites

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel in the ATP-binding cassette (ABC) transporter protein family. In the presence of ATP and physiologically relevant concentrations of AMP, CFTR exhibits adenylate kinase activity (ATP + AMP ⇆ 2 ADP). Previous studies suggested that the interaction of nucleotide triphosphate with CFTR at ATP-binding site 2 is required for this activity. Two other ABC proteins, Rad50 and a structural maintenance of chromosome protein, also have adenylate kinase activity. All three ABC adenylate kinases bind and hydrolyze ATP in the absence of other nucleotides. However, little is known about how an ABC adenylate kinase interacts with ATP and AMP when both are present. Based on data from non-ABC adenylate kinases, we hypothesized that ATP and AMP mutually influence their interaction with CFTR at separate binding sites. We further hypothesized that only one of the two CFTR ATP-binding sites is involved in the adenylate kinase reaction. We found that 8-azidoadenosine 5′-triphosphate (8-N₃-ATP) and 8-azidoadenosine 5′-monophosphate (8-N₃-AMP) photolabeled separate sites in CFTR. Labeling of the AMP-binding site with 8-N₃-AMP required the presence of ATP. Conversely, AMP enhanced photolabeling with 8-N₃-ATP at ATP-binding site 2. The adenylate kinase active center probe P¹,P⁵-di(adenosine-5′) pentaphosphate interacted simultaneously with an AMP-binding site and ATP-binding site 2. These results show that ATP and AMP interact with separate binding sites but mutually influence their interaction with the ABC adenylate kinase CFTR. They further indicate that the active center of the adenylate kinase comprises ATP-binding site 2.

Molecular dynamics studies on the conformational transitions of adenylate kinase: a computational evidence for the conformational selection mechanism

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.

Structural topology and activation of an initial adenylate kinase-substrate complex

Enzymatic activity is ultimately defined by the structure, chemistry, and dynamics of the Michaelis complex. A large number of experimentally determined structures between enzymes and substrates, substrate analogues, or inhibitors exist. However, transient, short-lived encounter and equilibrium structures also play fundamental roles during enzymatic reaction cycles. Such structures are inherently difficult to study with conventional experimental techniques. The enzyme adenylate kinase undergoes major conformational rearrangements in response to binding of its substrates, ATP and AMP. ATP is sandwiched between two binding surfaces in the closed and active enzyme conformation. Thus, adenylate kinase harbors two spatially distant surfaces in the substrate free open conformation, of which one is responsible for the initial interaction with ATP. Here, we have performed primarily nuclear magnetic resonance experiments on Escherichia coli adenylate kinase (AK(eco)) variants that allowed identification of the site responsible for the initial ATP interaction. This allowed a characterization of the structural topology of an initial equilibrium complex between AK(eco) and ATP. On the basis of the results, we suggest that the ATP binding mechanism for AK(eco) is a mixture between "induced fit" and "conformational selection" models. It is shown that ATP is activated in the initial enzyme-bound complex because it displays an appreciable rate of nonproductive ATP hydrolysis. In summary, our results provide novel structural and functional insights into adenylate kinase catalysis.

Elucidation of Mg²⁺ binding activity of adenylate kinase from Mycobacterium tuberculosis H₃₇Rv using fluorescence studies

Adenylate kinase (AK) is a small ubiquitous enzyme that catalyzes the reversible transfer of the terminal phosphate group from adenine triphosphate (ATP): magnesium ion (Mg²⁺) to adenine monophosphate (AMP) to form two molecules of adenine diphosphate (ADP). AK thus maintains the homeostasis of adenine nucleotides in eukaryotes and prokaryotes. Because the [ATP]/[ADP] ratio is an important parameter in energy regulation in cells, Mg²⁺-activated AK has an important biological role, particularly in the case of bacteria, as imbalance in the ratio of [ATP]/[ADP] has been associated with alteration in its DNA supercoiling state. In the present study, magnesium-binding assays were carried out by systematically varying the concentrations of Mg²⁺, protein, AMP, ATP, and indicator in kinetic experiments. We report evidence that during magnesium-binding assay, the fluorescence level of the indicator "Mag-Indo-1" changes with protein concentration, suggesting that magnesium ions are binding to AK. The dual activity of AK both as nucleoside monophosphate and diphosphate kinases suggests that this enzyme may have a role in RNA and DNA biosynthesis in addition to its role in intracellular nucleotide metabolism. According to the proposed model, the magnesium-activated AK exhibits an increase in its forward reaction rate compared with the inactivated form. These findings imply that Mg²⁺ could be an important regulator in the energy signaling network in cells.

Exploring the Dynamic Functional Landscape of Adenylate Kinase Modulated by Substrates

Adenylate kinase (ADK) has been explored widely, through both experimental and theoretical studies. However, still less is known about how the functional dynamics of ADK is modulated explicitly by its natural substrates. Here, we report a quantitative study of the dynamic energy landscape for ADK responding to the substrate binding by integrating both experimental investigations and theoretical modeling. We make theoretical predictions which are in remarkable agreement with the single molecule experiments on the substrate-bound complex. With our combined models of ADK in its apo form, in the presence of AMP or ATP, and in complex with both substrates, we specifically address the following key questions: (1) Are there intermediate state(s) during their catalytic cycle and if so how many? (2) How many pathways are there along the open-to-closed transitions and what are their corresponding weights? (3) How do substrates influence the pathway weights and the stability of the intermediates? (4) Which lid's motion is rate-limiting along the turnover cycle, the NMP or the LID domain? Our models predict two major parallel stepwise pathways and two on-pathway intermediates which are denoted as IN (NMP domain open while LID domain closed) and IL (LID domain open and NMP domain closed), respectively. Further investigation of temperature effects suggests that the IN pathway is dominant at room temperature, but the IL pathway is dominant at the optimal temperature. This leads us to propose that the IL pathway is more dominant by entropy and IN pathway by enthalpy. Remarkably, our results show that even with maximum concentrations of natural substrates, ADK still fluctuates between multiple functional states, reflecting an intrinsic capability of large-scale conformational fluctuations which may be essential to its biological function. The results based on the dual-ligands model provide the theoretical validation of random bisubstrate biproducts (Bi-Bi) mechanism for the enzymatic reaction of ADK. Additionally, the pathway flux analysis strongly suggests that the motion of the NMP domain is the rate-determining step for the conformational cycle (opening and closing).

Demonstration of phosphoryl group transfer indicates that the ATP-binding cassette (ABC) transporter cystic fibrosis transmembrane conductance regulator (CFTR) exhibits adenylate kinase activity

Cystic fibrosis transmembrane conductance regulator (CFTR) is a membrane-spanning adenosine 5'-triphosphate (ATP)-binding cassette (ABC) transporter. ABC transporters and other nuclear and cytoplasmic ABC proteins have ATPase activity that is coupled to their biological function. Recent studies with CFTR and two nonmembrane-bound ABC proteins, the DNA repair enzyme Rad50 and a structural maintenance of chromosome (SMC) protein, challenge the model that the function of all ABC proteins depends solely on their associated ATPase activity. Patch clamp studies indicated that in the presence of physiologically relevant concentrations of adenosine 5'-monophosphate (AMP), CFTR Cl(-) channel function is coupled to adenylate kinase activity (ATP+AMP <==> 2 ADP). Work with Rad50 and SMC showed that these enzymes catalyze both ATPase and adenylate kinase reactions. However, despite the supportive electrophysiological results with CFTR, there are no biochemical data demonstrating intrinsic adenylate kinase activity of a membrane-bound ABC transporter. We developed a biochemical assay for adenylate kinase activity, in which the radioactive γ-phosphate of a nucleotide triphosphate could transfer to a photoactivatable AMP analog. UV irradiation could then trap the (32)P on the adenylate kinase. With this assay, we discovered phosphoryl group transfer that labeled CFTR, thereby demonstrating its adenylate kinase activity. Our results also suggested that the interaction of nucleotide triphosphate with CFTR at ATP-binding site 2 is required for adenylate kinase activity. These biochemical data complement earlier biophysical studies of CFTR and indicate that the ABC transporter CFTR can function as an adenylate kinase.

Deprotonated imidodiphosphate in AMPPNP-containing protein structures

Many different proteins utilize the chemical energy provided by the cofactor adenosine triphosphate (ATP) for their proper function. A number of structures in the Protein Data Bank (PDB) contain adenosine 5'-(β,γ-imido)triphosphate (AMPPNP), a nonhydrolysable analog of ATP in which the bridging O atom between the two terminal phosphate groups is substituted by the imido function. Under mild conditions imides do not have acidic properties and thus the imide nitrogen should be protonated. However, an analysis of protein structures containing AMPPNP reveals that the imide group is deprotonated in certain complexes if the negative charges of the phosphate moieties in AMPPNP are in part neutralized by coordinating divalent metals or a guanidinium group of an arginine.

hCINAP is an atypical mammalian nuclear adenylate kinase with an ATPase motif: structural and functional studies

Human coilin interacting nuclear ATPase protein (hCINAP) directly interacts with coilin, a marker protein of Cajal Bodies (CBs), nuclear organelles involved in the maturation of small nuclear ribonucleoproteins UsnRNPs and snoRNPs. hCINAP has previously been designated as an adenylate kinase (AK6), but is very atypical as it exhibits unusually broad substrate specificity, structural features characteristic of ATPase/GTPase proteins (Walker motifs A and B) and also intrinsic ATPase activity. Despite its intriguing structure, unique properties and cellular localization, the enzymatic mechanism and biological function of hCINAP have remained poorly characterized. Here, we offer the first high-resolution structure of hCINAP in complex with the substrate ADP (and dADP), the structure of hCINAP with a sulfate ion bound at the AMP binding site, and the structure of the ternary complex hCINAP-Mg(2+) ADP-Pi. Induced fit docking calculations are used to predict the structure of the hCINAP-Mg(2+) ATP-AMP ternary complex. Structural analysis suggested a functional role for His79 in the Walker B motif. Kinetic analysis of mutant hCINAP-H79G indicates that His79 affects both AK and ATPase catalytic efficiency and induces homodimer formation. Finally, we show that in vivo expression of hCINAP-H79G in human cells is toxic and drastically deregulates the number and appearance of CBs in the cell nucleus. Our findings suggest that hCINAP may not simply regulate nucleotide homeostasis, but may have broader functionality, including control of CB assembly and disassembly in the nucleus of human cells.

Modulation of functionally significant conformational equilibria in adenylate kinase by high concentrations of trimethylamine oxide attributed to volume exclusion

The effect of an inert small molecule osmolyte, trimethyl amine N-oxide (TMAO), upon the conformational equilibria of Escherichia coli adenylate kinase was studied using time-resolved FRET. The relative populations of open and closed clefts between the LID and the CORE domains were measured as functions of the concentrations of the substrate ATP over the concentration range 0-18 mM and TMAO over the concentration range 0-4 M. A model was constructed according to which the enzyme exists in equilibrium among four conformational states, corresponding to combinations of open and closed conformations of the LID-CORE and AMP-CORE clefts. ATP is assumed to bind only to those conformations with the closed LID-CORE cleft, and TMAO is assumed to be differentially excluded as a hard spherical particle from each of the four conformations in accordance with calculations based upon x-ray crystallographic structures. This model was found to describe quantitatively the dependence of the fraction of the closed LID-CORE cleft upon the concentrations of both ATP and TMAO over the entire range of concentrations with just five undetermined parameters.

Seeing the forest for the trees: fluorescence studies of single enzymes in the context of ensemble experiments

Enzymes are remarkable molecular machines that make many difficult biochemical reactions possible under mild biological conditions with incredible precision and efficiency. Our understanding of the working principles of enzymes, however, has not reached the level where one can readily deduce the mechanism and the catalytic rates from an enzyme's structure. Resolving the dynamics that relate the three-dimensional structure of an enzyme to its function has been identified as a key issue. While still challenging to implement, single-molecule techniques have emerged as one of the most useful methods for studying enzymes. We review enzymes studied using single-molecule fluorescent methods but placing them in the context of results from other complementary experimental work done on bulk samples. This review primarily covers three enzyme systems--flavoenzymes, dehydrofolate reductase, and adenylate kinase--with additional enzymes mentioned where appropriate. When the single-molecule experiments are discussed together with other methods aiming at the same scientific question, the weakness, strength, and unique contributions become clear.

Diversity of function-related conformational changes in proteins: coordinate uncertainty, fragment rigidity, and stability

It was found that the variety of function-related conformational changes ("movements") in proteins is beyond the earlier simple classifications. Here we offer biochemists a more comprehensive, transparent, and easy-to-use approach allowing a detailed and accurate interpretation of such conformational changes. It makes possible a more multifaceted characterization of protein flexibility via identification of rigidly and nonrigidly repositioned fragments, stable and nonstable fragments, and domain and nondomain repositioning. "Coordinate uncertainty thresholds" derived from computed differences between independently determined coordinates of the same molecules are used as the criteria for conformational identity. "Identical" rigid substructures are localized in the distance difference matrices (DDMs). A sequence of simple transformations determines whether a structural change occurs by rigid-body movements of fragments or largely through non-rigid-body deformations. We estimate the stability of protein fragments and compare stable and rigidly moving fragments. The motions computed with the coarse-grained elastic networks are also compared to those of their DDM analogues. We study and suggest a classification for 17 structural pairs, differing in their functional states. For five of the 17 proteins, conformational change cannot be accomplished by rigid-body transformations and requires significant non-rigid-body deformations. Stable fragments rarely coincide with rigidly moving fragments and often disagree with the CATH identifications of domains. Almost all monomeric apo chains, containing stable fragments and/or domains, indicate instability of the entire molecule, suggesting the importance of fragments and domains motions prior to stabilization by substrate binding or crystallization. Notably, kinases exhibit the greatest extent of nonrigidity among the proteins investigated.