Mechanism of adenylate kinase. The "essential lysine" helps to orient the phosphates and the active site residues to proper conformations

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.

Mechanistic Basis for a Connection between the Catalytic Step and Slow Opening Dynamics of Adenylate Kinase

Escherichia coli adenylate kinase (AdK) is a small, monomeric enzyme that synchronizes the catalytic step with the enzyme's conformational dynamics to optimize a phosphoryl transfer reaction and the subsequent release of the product. Guided by experimental measurements of low catalytic activity in seven single-point mutation AdK variants (K13Q, R36A, R88A, R123A, R156K, R167A, and D158A), we utilized classical mechanical simulations to probe mutant dynamics linked to product release, and quantum mechanical and molecular mechanical calculations to compute a free energy barrier for the catalytic event. The goal was to establish a mechanistic connection between the two activities. Our calculations of the free energy barriers in AdK variants were in line with those from experiments, and conformational dynamics consistently demonstrated an enhanced tendency toward enzyme opening. This indicates that the catalytic residues in the wild-type AdK serve a dual role in this enzyme's function─one to lower the energy barrier for the phosphoryl transfer reaction and another to delay enzyme opening, maintaining it in a catalytically active, closed conformation for long enough to enable the subsequent chemical step. Our study also discovers that while each catalytic residue individually contributes to facilitating the catalysis, R36, R123, R156, R167, and D158 are organized in a tightly coordinated interaction network and collectively modulate AdK's conformational transitions. Unlike the existing notion of product release being rate-limiting, our results suggest a mechanistic interconnection between the chemical step and the enzyme's conformational dynamics acting as the bottleneck of the catalytic process. Our results also suggest that the enzyme's active site has evolved to optimize the chemical reaction step while slowing down the overall opening dynamics of the enzyme.

Proteasome-dependent protein quality control of the peroxisomal membrane protein Pxa1p

Peroxisomes are eukaryotic organelles that function in numerous metabolic pathways and defects in peroxisome function can cause serious developmental brain disorders such as adrenoleukodystrophy (ALD). Peroxisomal membrane proteins (PMPs) play a crucial role in regulating peroxisome function. Therefore, PMP homeostasis is vital for peroxisome function. Recently, we established that certain PMPs are degraded by the Ubiquitin Proteasome System yet little is known about how faulty/non-functional PMPs undergo quality control. Here we have investigated the degradation of Pxa1p, a fatty acid transporter in the yeast Saccharomyces cerevisiae. Pxa1p is a homologue of the human protein ALDP and mutations in ALDP result in the severe disorder ALD. By introducing two corresponding ALDP mutations into Pxa1p (Pxa1^MUT), fused to mGFP, we show that Pxa1^MUT-mGFP is rapidly degraded from peroxisomes in a proteasome-dependent manner, while wild type Pxa1-mGFP remains relatively stable. Furthermore, we identify a role for the ubiquitin ligase Ufd4p in Pxa1^MUT-mGFP degradation. Finally, we establish that inhibiting Pxa1^MUT-mGFP degradation results in a partial rescue of Pxa1p activity in cells. Together, our data demonstrate that faulty PMPs can undergo proteasome-dependent quality control. Furthermore, our observations may provide new insights into the role of ALDP degradation in ALD.

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.

Energetics and structural characterization of the large-scale functional motion of adenylate kinase

Adenylate Kinase (AK) is a signal transducing protein that regulates cellular energy homeostasis balancing between different conformations. An alteration of its activity can lead to severe pathologies such as heart failure, cancer and neurodegenerative diseases. A comprehensive elucidation of the large-scale conformational motions that rule the functional mechanism of this enzyme is of great value to guide rationally the development of new medications. Here using a metadynamics-based computational protocol we elucidate the thermodynamics and structural properties underlying the AK functional transitions. The free energy estimation of the conformational motions of the enzyme allows characterizing the sequence of events that regulate its action. We reveal the atomistic details of the most relevant enzyme states, identifying residues such as Arg119 and Lys13, which play a key role during the conformational transitions and represent druggable spots to design enzyme inhibitors. Our study offers tools that open new areas of investigation on large-scale motion in proteins.

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.

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.

Spectroscopy study on the noncovalent interactions in the binary and ternary systems of L-lysine, adenosine 5'-triphosphate and magnesium ions

Intermolecular interactions of adenosine 5'-triphosphate (ATP) with Lysine (Lys) and Mg(2+) were studied in aqueous solution by (1)H and (31)P NMR spectra. In the metal-free system, the N-1 atom of the purine ring of ATP and carboxyl group of Lys are the interaction sites at low pH conditions. With increasing pH, the interaction efficiency between the phosphate group of ATP and the protonated ammonium group of Lys increased significantly, while that with carboxyl group in Lys decreased. In the Mg(2+)-Lys-ATP system, multi-interactions, such as coordination, cations (Mg(2+), NH(3)(+))-π, hydrogen bonding, ion-pairing interactions and electrostatic interactions co-existed. In addition, the recognition of ATP by the amino acid cation (Lys) was significantly promoted by the addition of magnesium ion, which led to the coordination competition between Lys and ATP.

Elucidating the site of protein-ATP binding by top-down mass spectrometry

A Fourier-transform ion cyclotron resonance (FT-ICR) top-down mass spectrometry strategy for determining the adenosine triphosphate (ATP)-binding site on chicken adenylate kinase is described. Noncovalent protein-ligand complexes are readily detected by electrospray ionization mass spectrometry (ESI-MS), but the ability to detect protein-ligand complexes depends on their stability in the gas phase. Previously, we showed that collisionally activated dissociation (CAD) of protein-nucleotide triphosphate complexes yield products from the dissociation of a covalent phosphate bond of the nucleotide with subsequent release of the nucleotide monophosphate (Yin, S. et al., J. Am. Soc. Mass Spectrom. 2008, 19, 1199-1208). The intrinsic stability of electrostatic interactions in the gas phase allows the diphosphate group to remain noncovalently bound to the protein. This feature is exploited to yield positional information on the site of ATP-binding on adenylate kinase. CAD and electron capture dissociation (ECD) of the adenylate kinase-ATP complex generate product ions bearing mono- and diphosphate groups from regions previously suggested as the ATP-binding pocket by NMR and crystallographic techniques. Top-down MS may be a viable tool to determine the ATP-binding sites on protein kinases and identify previously unknown protein kinases in a functional proteomics study.

Productive versus unproductive nucleotide binding in yeast guanylate kinase mutants: comparison of R41M with K14M by proton two dimensional transferred NOESY

The R41M and K14M mutant enzymes of yeast guanylate kinase (GKy) were studied to investigate the effects of these site-directed mutations on bound-substrate conformations. Published X-ray crystal structures of yeast guanylate kinase indicate that K14 is part of the "P" loop involved in ATP and ADP binding, while R41 is suggested as a hydrogen bonding partner for the phosphoryl moiety of GMP. Both of these residues might be involved in transition state stabilization. Adenosine conformations of ATP and ADP and guanosine conformations of GMP bound to R41M and K14M mutant yeast guanylate kinase in the complexes GKy.MgATP, GKy.MgADP, and GKy.MgADP.[u-(13)C]GMP were determined by two-dimensional transferred nuclear Overhauser effect (TRNOESY) measurements combined with molecular dynamics simulations, and these conformations were compared with previously published conformations for the wild type. In the fully constrained, two substrate complexes, GKy.MgADP.[u-(13)C]GMP, the guanyl glycosidic torsion angle, chi, is 51 +/- 5 degrees for R41M and 47 +/- 5 degrees for K14M. Both are similar to the published 50 +/- 5 degrees published for wild type. For R41M with adenyl nucleotides, the glycosidic torsion angle, chi, was 55 +/- 5 degrees with MgATP, and 47 +/- 5 degrees with MgADP, which compares well to the 54 +/- 5 degrees published for wild type. However, for K14M with adenyl nucleotides, the glycosidic torsion angle was 30 +/- 5 degrees with MgATP and 28 +/- 5 degrees with MgADP. The results indicate that bound adenyl-nucleotides have significantly different conformations in the wild-type and K14M mutant enzymes, suggesting that K14 plays an important role in orienting the triphosphate of MgATP for catalysis.

Role of CFTR's intrinsic adenylate kinase activity in gating of the Cl(-) channel

The cystic fibrosis transmembrane conductance regulator (CFTR) is a Cl(-)channel in the ATP-binding cassette (ABC) transporter protein family. CFTR features the modular design characteristic of ABC transporters, which includes two membrane-spanning domains forming the channel pore, and two ABC nucleotide-binding domains that interact with ATP and contain the enzymatic activity coupled to normal gating. Like other ABC transporters CFTR is an ATPase (ATP + H(2)O --> ADP + Pi). Recent work has shown that CFTR also possesses intrinsic adenylate kinase activity (ATP + AMP left arrow over right arrow ADP + ADP). This finding raises important questions: How does AMP influence CFTR gating? Why does ADP inhibit CFTR current? Which enzymatic activity gates CFTR in vivo? Are there implications for other ABC transporters? This minireview attempts to shed light on these questions by summarizing recent advances in our understanding of the role of the CFTR adenylate kinase activity for channel gating.

Erythrocyte adenylate kinase deficiency: characterization of recombinant mutant forms and relationship with nonspherocytic hemolytic anemia

OBJECTIVE

Red cell adenylate kinase (AK) deficiency is a rare hereditary erythroenzymopathy associated with moderate to severe nonspherocytic hemolytic anemia and, in some cases, with mental retardation and psychomotor impairment. To date, diagnosis of AK deficiency depends upon demonstration of low enzyme activity in red blood cells and detection of mutations in AK1 gene. To investigate the molecular bases of the AK deficiency, we characterized five variants of AK1 isoenzyme-bearing mutations (118G>A, 190G>A, 382C>T, 418-420del, and 491A>G) found in AK-deficient patients with chronic hemolytic anemia.

MATERIALS AND METHODS

The complete AK1 cDNA was obtained by standard procedures and using as template the reticulocyte RNA. The cDNA was cloned in a plasmid vector and the enzyme was expressed in Escherichia coli BL21(DE3)pLysS, and purified by standard protocols to homogeneity. DNA mutants bearing point mutations were obtained from the cloned wild-type cDNA using standard methods of site-directed mutagenesis, whereas the DNA mutant with deletion of codon 140 was obtained by a two-step method.

RESULTS

Four mutant enzymes (Gly40Arg, Gly64Arg, Arg128Trp, Asp140del) were severely affected in activity, displaying a catalytic efficiency of four orders of magnitude lower than the wild-type; one (Tyr164Cys) was grossly perturbed in protein stability.

CONCLUSIONS

The altered properties displayed by the mutant enzymes support the cause-effect relationship between AK1 mutations and hemolytic anemia.

Systematic comparison of catalytic mechanisms of hydrolysis and transfer reactions classified in the EzCatDB database

Catalytic mechanisms of 270 enzymes from 131 superfamilies, mainly hydrolases and transferases, were analyzed based on their enzyme structures. A method of systematic comparison and classification of the catalytic reactions was developed. Hydrolysis and transfer reactions closely resemble one another, displaying common mechanisms, single displacement, and double displacement. These displacement mechanisms might be further subclassified according to the type of catalytic factors and nucleophilic substitution involved. Several types of catalytic factors exist: nucleophile, acid, base, stabilizer, modulator, cofactors. Nucleophilic substitution might be categorized as S(N)1/S(N)2 (or dissociative/associative) reactions. The classification indicates that some mechanisms favor particular types of catalytic factors. In hydrolyses of amide bonds and phosphoric ester bonds, mechanisms with single displacement tend to use inorganic cofactors such as zinc and magnesium ions as important catalysts, whereas those with double displacement frequently do not use such cofactors. In contrast, hydrolyses of O-glycoside bond rarely use such cofactors, with one exception. The trypsin-like hydrolytic reaction, which is catalyzed by the classic catalytic triad comprising serine/histidine/aspartate, can be considered as a "super-reaction" because it is observed in at least three nonhomologous enzymes, whereas most reactions are singlets without any nonhomologous enzymes. By dividing complex reactions into several reactions, correlations between active site structures and catalytic functions can be suggested. This classification method is applicable to other reactions such as elimination and isomerization. Furthermore, it will facilitate annotation of enzyme functions from 3D patterns of enzyme active sites. The classification is available at http://mbs.cbrc.jp/EzCatDB/RLCP/index.html.

Essential dynamics sampling study of adenylate kinase: Comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.

Molecular dynamics of apo-adenylate kinase: a distance replica exchange method for the free energy of conformational fluctuations

A large domain motion in adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes a large-scale rearrangement of its lid and AMP-binding domains when the open form closes over its substrates, AMP, and Mg2+-ATP, whereby the AMP-binding and lid domains come closer to the core. The third domain, the core, is relatively stable during this motion. A reaction coordinate that monitors the distance between the AMP-binding and core domains is selected to be able to compare with the results of energy transfer experiments. Sampling along this reaction coordinate is carried out by using a distance replica exchange method (DREM), where systems that differ by a restraint potential enforcing different reaction coordinate values are independently simulated with periodic attempts at exchange of these systems. Several methods are used to study the efficiency and convergence properties of the DREM simulation and compared with an analogous non-DREM simulation. The DREM greatly accelerates the rate and extent of configurational sampling and leads to equilibrium sampling as measured by monitoring collective modes obtained from a principal coordinate analysis. The potential of mean force along the reaction coordinate reveals a rather flat region for distances from the open to a relatively closed AKE conformation. The potential of mean force for smaller distances has a distinct minimum that is quite close to that found in the closed form X-ray structure. In concert with a decrease in the reaction coordinate distance (AMP-binding-to-core distance) the lid-to-core distance of AKE also decreases. Therefore, apo AKE can fluctuate from its open form to conformations that are quite similar to its closed form X-ray structure, even in the absence of its substrates.

Analyses of co-operative transitions in Plasmodium falciparum beta-ketoacyl acyl carrier protein reductase upon co-factor and acyl carrier protein binding

The type II fatty acid synthase pathway of Plasmodium falciparum is a validated unique target for developing novel antimalarials because of its intrinsic differences from the type I pathway operating in humans. beta-Ketoacyl-acyl carrier protein reductase is the only enzyme of this pathway that has no isoforms and thus selective inhibitors can be developed for this player of the pathway. We report here intensive studies on the direct interactions of Plasmodiumbeta-ketoacyl-acyl carrier protein reductase with its cofactor, NADPH, acyl carrier protein, acetoacetyl-coenzyme A and other ligands in solution, by monitoring the intrinsic fluorescence (lambdamax 334 nM) of the protein as a result of its lone tryptophan, as well as the fluorescence of NADPH (lambdamax 450 nM) upon binding to the enzyme. Binding of the reduced cofactor makes the enzyme catalytically efficient, as it increases the binding affinity of the substrate, acetoacetyl-coenzyme A, by 16-fold. The binding affinity of acyl carrier protein to the enzyme also increases by approximately threefold upon NADPH binding. Plasmodiumbeta-ketoacyl-acyl carrier protein reductase exhibits negative, homotropic co-operative binding for NADPH, which is enhanced in the presence of acyl carrier protein. Acyl carrier protein increases the accessibility of NADPH to beta-ketoacyl-acyl carrier protein reductase, as evident from the increase in the accessibility of the tryptophan of beta-ketoacyl-acyl carrier protein reductase to acrylamide, from 81 to 98%. In the presence of NADP+, the reaction proceeds in the reverse direction (Ka=23.17 microM-1). These findings provide impetus for exploring the influence of ligands on the structure-activity relationship of Plasmodiumbeta-ketoacyl-acyl carrier protein reductase.

Structures of S. aureus thymidylate kinase reveal an atypical active site configuration and an intermediate conformational state upon substrate binding

Methicillin-resistant Staphylococcus aureus (MRSA) poses a major threat to human health, particularly through hospital acquired infection. The spread of MRSA means that novel targets are required to develop potential inhibitors to combat infections caused by such drug-resistant bacteria. Thymidylate kinase (TMK) is attractive as an antibacterial target as it is essential for providing components for DNA synthesis. Here, we report crystal structures of unliganded and thymidylate-bound forms of S. aureus thymidylate kinase (SaTMK). His-tagged and untagged SaTMK crystallize with differing lattice packing and show variations in conformational states for unliganded and thymidylate (TMP) bound forms. In addition to open and closed forms of SaTMK, an intermediate conformation in TMP binding is observed, in which the site is partially closed. Analysis of these structures indicates a sequence of events upon TMP binding, with helix alpha3 shifting position initially, followed by movement of alpha2 to close the substrate site. In addition, we observe significant conformational differences in the TMP-binding site in SaTMK as compared to available TMK structures from other bacterial species, Escherichia coli and Mycobacterium tuberculosis as well as human TMK. In SaTMK, Arg 48 is situated at the base of the TMP-binding site, close to the thymine ring, whereas a cis-proline occupies the equivalent position in other TMKs. The observed TMK structural differences mean that design of compounds highly specific for the S. aureus enzyme looks possible; such inhibitors could minimize the transfer of drug resistance between different bacterial species.

Identification, mutagenesis, and transcriptional analysis of the methanesulfonate transport operon of Methylosulfonomonas methylovora

Recently identified genes located downstream (3') of the msmEF (transport encoding) gene cluster, msmGH, and located 5' of the structural genes for methanesulfonate monooxygenase (MSAMO) are described from Methylosulfonomonas methylovora. Sequence analysis of the derived polypeptide sequences encoded by these genes revealed a high degree of identity to ABC-type transporters. MsmE showed similarity to a putative periplasmic substrate binding protein, MsmF resembled an integral membrane-associated protein, and MsmG was a putative ATP-binding enzyme. MsmH was thought to be the cognate permease component of the sulfonate transport system. The close association of these putative transport genes to the MSAMO structural genes msmABCD suggested a role for these genes in transport of methanesulfonic acid (MSA) into M. methylovora. msmEFGH and msmABCD constituted two operons for the coordinated expression of MSAMO and the MSA transporter systems. Reverse-transcription-PCR analysis of msmABCD and msmEFGH revealed differential expression of these genes during growth on MSA and methanol. The msmEFGH operon was constitutively expressed, whereas MSA induced expression of msmABCD. A mutant defective in msmE had considerably slower growth rates than the wild type, thus supporting the proposed role of MsmE in the transport of MSA into M. methylovora.

Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase

The plastic network model (PNM) is used to generate a conformational change pathway for Escherichia coli adenylate kinase based on two crystal structures, namely that of an open and a closed conformer. In this model, the energy basins corresponding to known conformers are connected at their lowest common energies. The results are used to evaluate and analyze the minimal energy pathways between these basins. The open to closed transition analysis provides an identification of hinges that is in agreement with the existing definitions based on the available X-ray structures. The elastic energy distribution and the C(alpha) pseudo-dihedral variation provide similar information on these hinges. The ensemble of the 45 published structures for this protein and closely related proteins is shown to always be within 3.0 A of the pathway, which corresponds to a conformational change between two end structures that differ by a C(alpha)-atom root-mean-squared deviation of 7.1A.

Structural and dynamic studies on ligand-free adenylate kinase from Mycobacterium tuberculosis revealed a closed conformation that can be related to the reduced catalytic activity

Tuberculosis is the leading cause of death worldwide from a single infectious disease. Search of new therapeutic tools requires the discovery and biochemical characterization of new potential targets among the bacterial proteins essential for the survival and virulence. Among them are the nucleoside monophosphate kinases, involved in the nucleotide biosynthesis. In this work, we determined the solution structure of adenylate kinase (AK) from Mycobacterium tuberculosis (AKmt), a protein of 181 residues that was found to be essential for bacterial survival. The structure was calculated by a simulated annealing protocol and energy minimization using experimental restraints, collected by nuclear magnetic resonance spectroscopy. The final, well-defined 20 NMR structures show an average root-mean-square deviation of 0.77 A for the backbone atoms in regular secondary structure segments. The protein has a central CORE domain, composed of a five-stranded parallel beta-sheet surrounded by seven alpha-helices, and two peripheral domains, AMPbd and LID. As compared to other crystallographic structures of free form AKs, AKmt is more compact, with the AMP(bd) domain closer to the CORE of the protein. Analysis of the (15)N relaxation data enabled us to obtain the global rotational correlation time (9.19 ns) and the generalized order parameters (S(2)) of amide vectors along the polypeptide sequence. The protein exhibits restricted movements on a picosecond to nanosecond time scale in the secondary structural regions with amplitudes characterized by an average S(2)() value of 0.87. The loops beta1/alpha1, beta2/alpha2, alpha2/alpha3, alpha3/alpha4, alpha4/beta3, beta3/alpha5, alpha6/alpha7 (LID), alpha7/alpha8, and beta5/alpha9 exhibit rapid fluctuations with enhanced amplitudes. These structural and dynamic features of AKmt may be related to its low catalytic activity that is 10-fold lower than in their eukaryote counterparts.

An intrinsic adenylate kinase activity regulates gating of the ABC transporter CFTR

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel in the ATP binding cassette (ABC) transporter family. Like other ABC transporters, it can hydrolyze ATP. Yet while ATP hydrolysis influences channel gating, it has long seemed puzzling that CFTR would require this reaction because anions flow passively through CFTR. Moreover, no other ion channel is known to require the large energy of ATP hydrolysis to gate. We found that CFTR also has adenylate kinase activity (ATP + AMP <=> ADP + ADP) that regulates gating. When functioning as an adenylate kinase, CFTR showed positive cooperativity for ATP suggesting its two nucleotide binding domains may dimerize. Thus, channel activity could be regulated by two different enzymatic reactions, ATPase and adenylate kinase, that share a common ATP binding site in the second nucleotide binding domain. At physiologic nucleotide concentrations, adenylate kinase activity, rather than ATPase activity may control gating, and therefore involve little energy consumption.

Mutational analysis defines the 5'-kinase and 3'-phosphatase active sites of T4 polynucleotide kinase

T4 polynucleotide kinase (Pnk) is a bifunctional 5'-kinase/3'-phosphatase that aids in the repair of broken termini in RNA by converting 3'-PO4/5'-OH ends into 3'-OH/5'-PO4 ends, which are then sealed by RNA ligase. Here we have employed site-directed mutagenesis (introducing 31 mutations at 16 positions) to locate candidate catalytic residues within the 301 amino acid Pnk polypeptide. We found that alanine substitutions for Arg38 and Arg126 inactivated the 5'-kinase, but spared the 3'-phosphatase activity. Conservative substitutions of lysine or glutamine for Arg38 and Arg126 did not restore 5'-kinase activity. These results, together with previous mutational studies, highlight a constellation of five amino acids (Lys15, Ser16, Asp35, Arg38 and Arg126) that likely comprise the 5'-kinase active site. Four of these residues are conserved at the active sites of adenylate kinases (Adk), suggesting that Pnk and Adk are structurally and mechanistically related. We found that alanine substitutions for Asp165, Asp167, Arg176, Arg213, Asp254 and Asp278 inactivated the 3'-phosphatase, but spared the 5'-kinase. Conservative substitutions of asparagine or glutamate for Asp165, Asp167 and Asp254 did not revive the 3'-phosphatase activity, nor did lysine substitutions for Arg176 and Arg213. Glutamate in lieu of Asp278 partially restored activity, whereas asparagine had no salutary effect. Alanine substitutions for Arg246 and Arg279 partially inactivated the 3'-phosphatase; the conservative R246K change restored activity, whereas R279K had no benefit. The essential phosphatase residues Asp165 and Asp167 are located within a 165DxDxT169 motif that defines a superfamily of phosphotransferases. Our data suggest that the 3'-phosphatase active site incorporates multiple additional functional groups.

Domain structure and mutational analysis of T4 polynucleotide kinase

T4 polynucleotide kinase (Pnk) is the founding member of a family of 5'-kinase/3'-phosphatase enzymes that heal broken termini in RNA or DNA by converting 3'-PO(4)/5'-OH ends into 3'-OH/5'-PO(4) ends, which are then suitable for sealing by RNA or DNA ligases. Here we employed site-directed mutagenesis and biochemical methods to dissect the domain structure of the homotetrameric T4 Pnk protein and to localize essential constituents of the apparently separate active sites for the 5'-kinase and 3'-phosphatase activities. We characterized deletion mutants Pnk(42-301) and Pnk(1-181), which correspond to domains defined by proteolysis with chymotrypsin. Pnk(1-181) is a monomer with no 3'-phosphatase and low residual 5'-kinase activity. Pnk(42-301) is a dimer with no 5'-kinase and low residual 3'-phosphatase activity. Four classes of missense mutational effects were observed. (i) Mutations K15A, S16A, and D35A inactivated the 5'-kinase but did not affect the 3'-phosphatase or the tetrameric quaternary structure of T4 Pnk. 5'-kinase activity was ablated by the conservative mutations K15R, K15Q, and D35N; however, kinase activity was restored by the S16T change. (ii) Mutation D167A inactivated the 3'-phosphatase without affecting the 5'-kinase or tetramerization. (iii) Mutation D85A caused a severe decrement in 5'-kinase activity and only a modest effect on the 3'-phosphatase; the nearby N87A mutation resulted in a significantly reduced 3'-phosphatase activity and slightly reduced 5'-kinase activity. D85A and N87A both affected the quaternary structure, resulting in a mixed population of tetramer and dimer species. (iv) Alanine mutations at 11 other conserved positions had no significant effect on either 5'-kinase or 3'-phosphatase activity.

Biochemical and X-ray crystallographic studies on shikimate kinase: the important structural role of the P-loop lysine

Shikimate kinase, despite low sequence identity, has been shown to be structurally a member of the nucleoside monophosphate (NMP) kinase family, which includes adenylate kinase. In this paper we have explored the roles of residues in the P-loop of shikimate kinase, which forms the binding site for nucleotides and is one of the most conserved structural features in proteins. In common with many members of the P-loop family, shikimate kinase contains a cysteine residue 2 amino acids upstream of the essential lysine residue; the side chains of these residues are shown to form an ion pair. The C13S mutant of shikimate kinase was found to be enzymatically active, whereas the K15M mutant was inactive. However, the latter mutant had both increased thermostability and affinity for ATP when compared to the wild-type enzyme. The structure of the K15M mutant protein has been determined at 1.8 A, and shows that the organization of the P-loop and flanking regions is heavily disturbed. This indicates that, besides its role in catalysis, the P-loop lysine also has an important structural role. The structure of the K15M mutant also reveals that the formation of an additional arginine/aspartate ion pair is the most likely reason for its increased thermostability. From studies of ligand binding it appears that, like adenylate kinase, shikimate kinase binds substrates randomly and in a synergistic fashion, indicating that the two enzymes have similar catalytic mechanisms.

Protein folding and function: the N-terminal fragment in adenylate kinase

Three-dimensional protein folds range from simple to highly complex architectures. In complex folds, some building block fragments are more important for correct protein folding than others. Such fragments are typically buried in the protein core and mediate interactions between other fragments. Here we present an automated, surface area-based algorithm that is able to indicate which, among all local elements of the structure, is critical for the formation of the native fold, and apply it to structurally well-characterized proteins. In particular, we focus on adenylate kinase. The fragment containing the phosphate binding, P-loop (the "giant anion hole") flanked by a beta-strand and an alpha-helix near the N-terminus, is identified as a critical building block. This building block shows a high degree of sequence and structural conservation in all adenylate kinases. The results of our molecular dynamics simulations are consistent with this identification. In its absence, the protein flips to a stable, non-native state. In this misfolded conformation, the other local elements of the structure are in their native-like conformations; however, their association is non-native. Furthermore, this element is critically important for the function of the enzyme, coupling folding, and function.

Phosphorylation of Wzc, a tyrosine autokinase, is essential for assembly of group 1 capsular polysaccharides in Escherichia coli

Wzc proteins are tyrosine autokinases. They are found in some important bacterial pathogens of humans and livestock as well as plant-associated bacteria, and are often encoded within gene clusters determining synthesis and assembly of capsular and extracellular polysaccharides. Autophosphorylation of Wzc(cps) is essential for assembly of the serotype K30 group 1 capsule in Escherichia coli O9a:K30, although a genetically unlinked Wzc(cps)-homologue (Etk) can also participate with low efficiency. While autophosphorylation of Wzc(cps) is required for assembly of high molecular weight K30 capsular polysaccharide, it is not essential for either the synthesis of the K30 repeat units or for activity of the K30 polymerase enzyme. Paradoxically, the cognate phosphotyrosine protein phosphatase for Wzc(cps), Wzb(cps), is also required for capsule expression. The tyrosine-rich domain at the C terminus of Wzc(cps) was identified as the site of phosphorylation and autophosphorylation of Wzc requires a functional Walker A motif. Intermolecular transphosphorylation of Wzc(cps) was detected in strains expressing a combination of mutant Wzc(cps) derivatives. The N- and C-terminal domains of Wzc(cps) were expressed independently to mimic the situation found naturally in Gram-positive bacteria. In this format, both domains were required for phosphorylation of the Wzc(cps) C terminus, and for capsule assembly. Regulation by a post-translational phosphorylation event represents a new dimension in the assembly of bacterial cell-surface polysaccharides.

Backbone dynamics of escherichia coli adenylate kinase at the extreme stages of the catalytic cycle studied by (15)N NMR relaxation

Adenylate kinase from Escherichia coli (AKeco), consisting of a single 23.6 kDa polypeptide chain folded into domains CORE, AMPbd, and LID, catalyzes the reaction AMP + ATP --> 2ADP. Domains LID and AMPbd execute large-scale movements during catalysis. Backbone dynamics of ligand-free and AP(5)A-inhibitor-bound AKeco were studied comparatively with (15)N NMR relaxation methods. Overall diffusion with correlation times of 15.05 (11.42) ns and anisotropy D(parallel)/D(perp) = 1.25 (1.10), and fast internal motions with correlation times up to 100 ps (50 ps), were determined for AKeco (AKecoAP(5)A). Fast internal motions affect 93% of the AKeco sites, with pronounced preference for domains AMPbd and LID, and 47% of the AKecoAP(5)A sites, with limited variability along the chain. The mean squared generalized order parameters, <S(2)>, of secondary structure elements and loops are affected by ligand binding differentially and in a domain-specific manner. Nanosecond motions predominate within AMPbd. Prominent exchange contributions, associated in particular with residue G10 of the nucleotide-binding P-loop motif, are interpreted to reflect hydrogen-bond dynamics at the inhibitor-binding site. The hypothesis of energetic counter balancing of substrate binding based on crystallographic data is strongly supported by the solution NMR results. Correlations between backbone dynamics and domain displacement are established.

wt p53 dependent expression of a membrane-associated isoform of adenylate kinase

Six novel p53-inducible transcripts were recently cloned from Val5, a murine cell line stably expressing a temperature-sensitive p53 allele. One of the isolated clones represented a novel isoform of cytosolic adenylate kinase (AK1), a highly conserved monomeric enzyme involved in cellular homeostasis of adenine nucleotides. The corresponding protein, which we named AK1beta, was specifically induced upon activation of wt p53 in Val5 cells. The AK1beta protein differs from cytoplasmic AK1 by having 18 extra amino acids at the N-terminus. The extra residues in AK1beta provide a consensus signal for N-terminal myristoylation; as expected, AK1beta was shown to localize to the plasma membrane. The human AK1 gene contains several consensus p53 binding sites and we report that p53-dependent induction of the alternative AK1beta transcript also occurs in human cells. By using antisense ablation experiments in Val5 fibroblasts we show that AK1beta plays a relevant role in the establishment of reversible cell-cycle arrest as induced by p53 in these cells. These findings suggest that within a p53-dependent genetic program, a specific isoform of adenylate kinase has a previously undescribed growth-regulatory function, which might not necessarily require its best characterized biochemical activity.

Role of lysine 240 in the mechanism of yeast pyruvate kinase catalysis

Site-directed mutagenesis was used to change Lys 240 of yeast pyruvate kinase (Lys 269 in muscle PK) to Met. K240M has an absolute requirement for FBP for catalysis. K240M is 100- and 1000-fold less active than wild-type YPK in the presence of Mn(2+) and Mg(2+), respectively. Steady-state fluorescence titration data suggest that the substrate PEP binds to K240M with the same affinity as it does to wild-type YPK. The rate of phosphoryl transfer in K240M has been decreased >1000-fold compared to wild-type YPK. The detritiation of 3-[(3)H]pyruvate catalyzed by YPK occurs at a rate significantly greater than the spontaneous rate. Detritiation of pyruvate by wild-type YPK occurs as a divalent metal- and FBP-dependent process requiring ATP. There is no detectable detritiation of pyruvate catalyzed by K240M. The solvent deuterium isotope effect on k(cat) is 2.7 +/- 0.2 and 1.6 +/- 0.1 for the wild type and for K240M YPK, respectively. This suggests that the isotope sensitive step in the PK reaction does not involve Lys 240 and that the enolpyruvate intermediate is still protonated by K240M. Isotope trapping was used to characterize enolpyruvate protonation by K240M. While there was enrichment of the methyl protons of pyruvate from labeled solvent formed by catalysis with muscle PK and wild-type YPK, only background levels of tritium were trapped with K240M. In K240M, the proton donor exchanges protons with the solvent at a higher rate relative to turnover than does the proton donor in wild-type YPK. The pH-rate profile of K240M exhibits the loss of a pK(a) value of 8. 8 observed with wild-type YPK. The above data and recent crystal structure data suggest that Lys 240 interacts with the phosphoryl group of phosphoenolpyruvate and helps to stabilize the pentavalent phosphate transition state during phosphoryl transfer. Phosphoryl transfer is highly coupled to proton transfer, or Lys 240 also affects enolate protonation.

Nucleoside monophosphate kinases: structure, mechanism, and substrate specificity

The catalytic mechanisms of adenylate kinase, guanylate kinase, uridylate kinase, and cytidylate kinase are reviewed in terms of kinetic and structural information that has been obtained in recent years. All four kinases share a highly related tertiary structure, characterized by a central five-stranded parallel beta-sheet with helices on both sides, as well as the three regions designated as the CORE, NMPbind, and LID domains. The catalytic mechanism continues to be refined to higher levels of resolution by iterative structure-function studies, and the strengths and limitations of site-directed mutagenesis are well illustrated in the case of adenylate kinase. The identity and roles of active site residues now appear to be resolved, and this review describes how specific site substitutions with unnatural amino acid side-chains have proven to be a major advance. Likewise, there is mounting evidence that phosphoryl transfer occurs by an associative transition state, based on (a) the stereochemical course of phosphoryl transfer, (b) geometric considerations, (c) examination of likely electronic distributions, (d) the orientation of the phosphoryl acceptor relative to the phosphoryl being transferred, (e) the most likely role of magnesium ion, (f) the lack of restricted access of solvent water, and (g) the results of oxygen-18 kinetic isotope. effect experiments.

3'-Phosphoadenosine 5'-phosphosulfate binding site of flavonol 3-sulfotransferase studied by affinity chromatography and 31P NMR

The function of Lys-59, Arg-141, and Arg-277 in PAPS binding and catalysis of the flavonol 3-sulfotransferase was investigated. Affinity chromatography of conservative mutants with PAPS analogues allowed us to determine that Lys-59 interacts with the 5' portion of the nucleotide, while Arg-141 interacts with the 3' portion, confirming assignments deduced from the crystal structure of mouse estrogen sulfotransferase [Kakuta, Y., Pedersen, L. G., Carter, C. W. , Negishi, M., and Pedersen, L. C. (1997) Nat. Struct. Biol. 4, 904-908]. The affinity chromatography method could be used to characterize site-directed mutants for other types of enzymes that bind nucleoside 3',5'- or 2',5'-diphosphates. 31P NMR spectra of enzyme-PAP complexes were recorded for the wild-type enzyme and K59R and K59A mutants. The results of these experiments suggest that Lys-59 is involved in the determination of the proper orientation of the phosphosulfate group for catalysis.

An ExeAB complex in the type II secretion pathway of Aeromonas hydrophila: effect of ATP-binding cassette mutations on complex formation and function

The energy-dependent secretion of aerolysin by Aeromonas hydrophila requires the ExeA and ExeB proteins. An 85 kDa complex containing the two proteins was identified in wild-type cells but not in cells producing either protein alone. Radiolabelling followed by cross-linking, immunoprecipitation and then reduction of the cross-links confirmed the presence of the two proteins in the same complex. The complex could also be extracted intact from cell membranes with non-ionic detergents. A G229D substitution in the kinase-3a motif of ExeA strongly reduced the level of aerolysin secretion, as did the replacement of the invariant Lys of the kinase-1a motif (K56) with Arg. The G229D mutant contained very little of the ExeA-ExeB complex, but overexpression of the mutant complex until wild-type levels were achieved allowed normal secretion. In contrast, the K56R mutation had no effect on complex formation, but normal secretion levels occurred only when there was a far greater amount of the complex present. These results are consistent with a model in which binding of ATP by ExeA is required for ExeA-ExeB complex formation, while hydrolysis is required for its function in secretion once established.

Direct demonstration of the catalytic role of binding interactions in an enzymatic reaction

It has been suggested that the fundamental feature that distinguishes enzymes from simple chemical catalysts is the ability of enzymes to use binding interactions for catalysis. Results with the Tetrahymena group I RNA enzyme described herein directly demonstrate the catalytic contributions of binding interactions. With wild-type ribozyme, specific functional groups at a distance from the site of chemical transformation facilitate substrate binding without accelerating reaction of bound substrate; with modified ribozymes, these functional groups provide the same overall energetic effect but instead accelerate reaction of bound substrate without increasing binding. These observations are quantitatively described by a structural framework that was established by previous results. The P1 duplex between the substrate and the ribozyme's recognition sequence exists in two states, the open complex, in which the substrate is localized to the ribozyme solely by base-pairing interactions, or the closed complex, in which the duplex is docked into tertiary interactions and positioned with respect to the catalytic groups in the active site. In the absence of sufficient binding energy to ensure stable docking in the ground state, added P1 functional groups accelerate reaction of the bound substrate by helping to overcome the energetic barrier for docking into the reactive, closed complex. When the functional groups present on the P1 duplex are sufficient to ensure stable docking in the closed complex, added functional groups give stronger binding without accelerating reaction of the bound substrate. This behavior is a manifestation of the inextricable link between binding interactions and catalysis. The conclusions also have implications for interpreting the effects of site-directed mutagenesis and for the evolution of active site interactions.

The adenylation domain of tyrocidine synthetase 1--structural and functional role of the interdomain linker region and the (S/T)GT(T/S)GXPKG core sequence

Sequence analysis of peptide synthetases revealed extensive structure similarity with firefly luciferase, whose crystal structure has recently become available, providing evidence for the localization of the active site at the interface between two subdomains separated by a distorted linker region [Conti, E., Franks, N. P. & Brick, P. (1996) Structure 4, 287-298]. The functional importance of two flexible loops, corresponding to the linker region of firefly luciferase and the highly conserved (S/T)GT(T/S)GXPKG core sequence, has been studied in view of the proposed conformational changes by the use of mutant analysis, limited proteolysis and chemical modification of tyrocidine synthetase 1. Substitution of the highly conserved Arg416, residing in the loop separating the subdomains of the adenylation domain, resulted in profound loss of activity. Limited proteolysis of the mutant suggested significant structural changes as manifested by lack of protection to degradation in the presence of substrates, revealing a probable disturbance of the induced-fit mechanism regulating the transformation from an open to a closed conformation. Mutants, obtained by replacement of the conserved Lys186 from the (S/T)GT(T/S)GXPKG core sequence, displayed only minor differences in substrate-binding affinity despite significant reduction of catalytic efficiency. Residue Lys186 appears to play an important role in either stabilization of the bound substrate through charge-charge-interactions, and/or fixing of the loop for maintainance of the active-site conformation.

Modelling the 2-kinase domain of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase on adenylate kinase

Simultaneous multiple alignment of available sequences of the bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase revealed several segments of conserved residues in the 2-kinase domain. The sequence of the kinase domain was also compared with proteins of known three-dimensional structure. No similarity was found between the kinase domain of 6-phosphofructo-2-kinase and 6-phosphofructo-1-kinase. This questions the modelling of the 2-kinase domain on bacterial 6-phosphofructo-1-kinase that has previously been proposed [Bazan, Fletterick and Pilkis (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 9642-9646]. However, sequence similarities were found between the 2-kinase domain and several nucleotide-binding proteins, the most similar being adenylate kinase. A structural model of the 2-kinase domain based on adenylate kinase is proposed. It accommodates all the results of site-directed mutagenesis studies carried out to date on residues in the 2-kinase domain. It also allows residues potentially involved in catalysis and/or substrate binding to be predicted.

Identification of amino acid residues critical for catalysis and cosubstrate binding in the flavonol 3-sulfotransferase

The comparison of the deduced amino acid sequences of plant and animal sulfotransferases (ST) has allowed the identification of four well conserved regions, and previous experimental evidence suggested that regions I and IV might be involved in the binding of the cosubstrate, 3'-phosphoadenosine 5'-phosphosulfate (PAPS). Moreover, region IV is homologous to the glycine-rich phosphate binding loop (P-loop) motif known to be involved in nucleotide phosphate binding in several protein families. In this study, the function of amino acid residues within these two regions was investigated by site-directed mutagenesis of the plant flavonol 3-ST. In region I, our results identify Lys59 as critical for catalysis, since replacement of this residue with alanine resulted in a 300-fold decrease in specific activity, while a 15-fold reduction was observed after the conservative replacement with arginine. Photoaffinity labeling of K59R and K59A with [35S]PAPS revealed that Lys59 is not required for cosubstrate binding. However, the K59A mutant had a reduced affinity for 3'-phosphoadenosine 5'-phosphate (PAP)-agarose, suggesting that Lys59 may participate in the stabilization of an intermediate during the reaction. In region IV, all substitutions of Arg276 resulted in a marked decrease in specific activity. Conservative and unconservative replacements of Arg276 resulted in weak photoaffinity labeling with [35S]PAPS and the R276A/T73A and R276E enzymes displayed reduced affinities for PAP-agarose, suggesting that the Arg276 side chain is required to bind the cosubstrate. The analysis of the kinetic constants of mutant enzymes at residues Lys277, Gly281, and Lys284 allowed to confirm that region IV is involved in cosubstrate binding.