Letter to the Editor: 1H, 13C and 15N backbone resonance assignment of Escherichia coli adenylate kinase, a 23.6 kDa protein

1H, 13C, 15N backbone resonance assignment of Escherichia coli adenylate kinase

Adenylate kinase reversibly catalyzes the conversion of ATP plus AMP to two ADPs. This essential catalyst is present in every cell, and the Escherichia coli protein is often employed as a model enzyme. Our aim is to use the E. coli enzyme to understand dry protein structure and protection. Here, we report the expression, purification, steady-state assay, NMR conditions and 1H, 13C, 15N backbone resonance NMR assignments of its C77S variant. These data will also help others utilize this prototypical enzyme.

Tracking the ATP-binding response in adenylate kinase in real time

The biological function of proteins is critically dependent on dynamics inherent to the native structure. Such structural dynamics obey a predefined order and temporal timing to execute the specific reaction. Determination of the cooperativity of key structural rearrangements requires monitoring protein reactions in real time. In this work, we used time-resolved x-ray solution scattering (TR-XSS) to visualize structural changes in the Escherichia coli adenylate kinase (AdK) enzyme upon laser-induced activation of a protected ATP substrate. A 4.3-ms transient intermediate showed partial closing of both the ATP- and AMP-binding domains, which indicates a cooperative closing mechanism. The ATP-binding domain also showed local unfolding and breaking of an Arg¹³¹-Asp¹⁴⁶ salt bridge. Nuclear magnetic resonance spectroscopy data identified similar unfolding in an Arg¹³¹Ala AdK mutant, which refolded in a closed, substrate-binding conformation. The observed structural dynamics agree with a “cracking mechanism” proposed to underlie global structural transformation, such as allostery, in proteins.

Structural basis for catalytically restrictive dynamics of a high-energy enzyme state

An emerging paradigm in enzymology is that transient high-energy structural states play crucial roles in enzymatic reaction cycles. Generally, these high-energy or ‘invisible' states cannot be studied directly at atomic resolution using existing structural and spectroscopic techniques owing to their low populations or short residence times. Here we report the direct NMR-based detection of the molecular topology and conformational dynamics of a catalytically indispensable high-energy state of an adenylate kinase variant. On the basis of matching energy barriers for conformational dynamics and catalytic turnover, it was found that the enzyme's catalytic activity is governed by its dynamic interconversion between the high-energy state and a ground state structure that was determined by X-ray crystallography. Our results show that it is possible to rationally tune enzymes' conformational dynamics and hence their catalytic power—a key aspect in rational design of enzymes catalysing novel reactions.

Adenylate kinase (AdK) plays a key role in cellular energy homeostasis by catalysing the reversible magnesium-dependent formation of ADP from AMP and ATP. Here the authors present a detailed analysis of adenylate kinase's conformational dynamics and characterize a high-energy state of AdK indispensable for catalysis.

Evidence for domain motion in proteins affecting global diffusion properties: a nuclear magnetic resonance study

The rotational diffusion of proteins is an important hydrodynamic property. Compact protein structures were found previously to exhibit hydration layer viscosity, etaloc, higher than the viscosity of bulk water, eta. This implies an apparent activation energy for rotational diffusion higher than the activation energy of water viscosity, Eeta=15.4+/-0.3 kJ/mol. In this study we examine etaloc of internally mobile proteins using 15N spin relaxation methods. We also examine the activation enthalpy, DeltaH#, and activation entropy, DeltaS#, for rotational diffusion. Of particular relevance are internally mobile ligand-free forms and compact ligand-bound forms of multidomain proteins. Adenylate kinase (AKeco) and Ca2+-calmodulin (Ca2+-CaM) are typical examples. For AKeco (Ca2+-CaM) we find that DeltaH# is 14.5+/-0.5 (15.7+/-0.4) kJ/mol. For the complex of AKeco with the inhibitor AP5A (the complex of Ca2+-CaM with the peptide smMLCKp), we find that DeltaH# is 18.1+/-0.7 (18.2+/-0.5) kJ/mol. The internally mobile outer surface protein A has DeltaH#=12.6+/-0.8 kJ/mol, and the compact protein Staphylococcal nuclease has DeltaH#=18.8+/-0.6 kJ/mol. For the internally mobile and compact proteins studied, <|DeltaS(|> equals 62+/-7 J/(mol K) and 44+/-5 J/(mol K), respectively. The fact is that etaloc>eta (DeltaH#>Eeta) for compact proteins was ascribed previously to electrostatic interactions between surface sites and water rigidifying the hydration layer. We find herein that obliteration of these interactions by domain motion leads to etaloc approximately eta, DeltaH# approximately Eeta, and large activation entropy for internally mobile protein structures.

Activation Energy of Catalysis-Related Domain Motion in E. coli Adenylate Kinase

Adenylate kinase from E. coli (AKeco), folded into domains CORE, AMPbd, and LID, catalyzes the reaction AMP + ATP <--> 2ADP. Previous X-ray crystallography and optical solution methods showed that the domains AMPbd and LID, and the conserved P-loop, execute large-amplitude catalysis-related motions. We used (15)N NMR spin relaxation methods to find that the simplified model-free (MF) analysis does not, whereas our general Slowly Relaxing Local Structure analysis does, detect catalytic domain motion. SRLS set for the first time the correlation time for domain motion at tau(L)perpendicular = 8.2 ns, to be compared with tau(m) = 15.1 ns for global tumbling. These results were obtained at 303 K. Herein we conduct a temperature-dependent investigation of tau(L)perpendicular and tau(m) in the range of 288-310 K. We found that the activation energy for global tumbling is Ea = 16.9 +/- 0.5 kJ/mol, the hydrodynamic volume of hydrated AKeco is 65.6 +/- 2.1 nm3, its radius is 2.50 +/- 0.03 nm, and the number of hydration layers is 1.77. The average tau(L)perpendicular value decreases from 11 ns at 288 K to 4 ns at 310 K, with activation energies of 29.7 +/- 3.3, 32.1 +/- 4.3, and 30.4 +/- 4.3 kJ/mol for the domains AMPbd and LID, and the catalytic P-loop, respectively. These values are two-to-three times smaller than typical activation energies of enzymatic reactions. Hence kinase catalysis appears not to be controlled by domain motion in the ligand-free enzyme. However, the latter process clearly facilitates important mechanical aspects such as steric recognition and capturing of the AMP and ATP substrates, their proper positioning for phosphorylation, and the release of the ADP product.

Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair

A fundamental question is how enzymes can accelerate chemical reactions. Catalysis is not only defined by actual chemical steps, but also by enzyme structure and dynamics. To investigate the role of protein dynamics in enzymatic turnover, we measured residue-specific protein dynamics in hyperthermophilic and mesophilic homologs of adenylate kinase during catalysis. A dynamic process, the opening of the nucleotide-binding lids, was found to be rate-limiting for both enzymes as measured by NMR relaxation. Moreover, we found that the reduced catalytic activity of the hyperthermophilic enzyme at ambient temperatures is caused solely by a slower lid-opening rate. This comparative and quantitative study of activity, structure and dynamics revealed a close link between protein dynamics and catalytic turnover.

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.

Domain flexibility in ligand-free and inhibitor-bound Escherichia coli adenylate kinase based on a mode-coupling analysis of 15N spin relaxation

Adenylate kinase from Escherichia coli (AKeco), consisting of a 23.6-kDa polypeptide chain folded into domains CORE, AMPbd, and LID catalyzes the reaction AMP + ATP <--> 2ADP. The domains AMPbd and LID execute large-amplitude movements during catalysis. Backbone dynamics of ligand-free and AP(5)A-inhibitor-bound AKeco is studied with slowly relaxing local structure (SRLS) (15)N relaxation, an approach particularly suited when the global (tau(m)) and the local (tau) motions are likely to be coupled. For AKeco tau(m) = 15.1 ns, whereas for AKeco*AP(5)A tau(m) = 11.6 ns. The CORE domain of AKeco features an average squared order parameter, <S(2)>, of 0.84 and correlation times tau(f) = 5-130 ps. Most of the AKeco*AP(5)A backbone features <S(2)> = 0.90 and tau(f) = 33-193 ps. These data are indicative of relative rigidity. Domains AMPbd and LID of AKeco, and loops beta(1)/alpha(1), alpha(2)/alpha(3), alpha(4)/beta(3), alpha(5)/beta(4), and beta(8)/alpha(7) of AKeco*AP(5)A, feature a novel type of protein flexibility consisting of nanosecond peptide plane reorientation about the C(i-1)(alpha)-C(i)(alpha) axis, with correlation time tau(perpendicular) = 5.6-11.3 ns. The other microdynamic parameters underlying this dynamic model include S(2) = 0.13-0.5, tau(parallel) on the ps time scale, and a diffusion tilt beta(MD) ranging from 12 to 21 degrees. For the ligand-free enzyme the tau(perpendicular) mode was shown to represent segmental domain motion, accompanied by conformational exchange contributions R(ex) < or = 4.4 s(-1). Loop alpha(4)/beta(3) and alpha(5)/beta(4) dynamics in AKeco*AP(5)A is related to the "energetic counter-balancing of substrate binding" effect apparently driving kinase catalysis. The other flexible AKeco*AP(5)A loops may relate to domain motion toward product release.

A novel view of domain flexibility in E. coli adenylate kinase based on structural mode-coupling (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. In the ligand-free enzyme the domains AMPbd and LID execute large-amplitude movements controlling substrate binding and product release during catalysis. Domain flexibility is investigated herein with the slowly relaxing local structure (SRLS) model for (15)N relaxation. SRLS accounts rigorously for coupling between the global and local N-H motions through a local ordering potential exerted by the protein structure at the N-H bond. The latter reorients with respect to its protein surroundings, which reorient on the slower time scale associated with the global protein tumbling. AKeco diffuses globally with correlation time tau(m)=15.1 ns, while locally two different dynamic cases prevail. The domain CORE features ordering about the equilibrium N-H bond orientation with order parameters, S(2), of 0.8-0.9 and local motional correlation times, tau, mainly between 5-130 ps. This represents a conventional rigid protein structure with rapid small-amplitude N-H fluctuations. The domains AMPbd and LID feature small parallel (Z(M)) ordering of S(2)=0.2-0.5 which can be reinterpreted as high perpendicular (Y(M)) ordering. M denotes the local ordering/local diffusion frame. Local motion about Z(M) is given by tau( parallel) approximately 5 ps and local motion of the effective Z(M) axis about Y(M) by tau( perpendicular)=6-11 ns. Z(M) is tilted at approximately 20 degrees from the N-H bond. The orientation of the Y(M) axis may be considered parallel to the C(alpha)(i-1)-C(alpha)(i) axis. The tau( perpendicular) mode reflects collective nanosecond peptide-plane motions, interpretable as domain motion. A powerful new model of protein flexibility/domain motion has been established. Conformational exchange (R(ex)) processes accompany the tau( perpendicular) mode. The SRLS analysis is compared with the conventional model-free analysis.

Multinuclear magnetic resonance studies of Escherichia coli adenylate kinase in free and bound forms. Resonance assignment, secondary structure and ligand binding

The crystal structure of Escherichia coli adenylate kinase (AKe) revealed three main components: a CORE domain, composed of a five-stranded parallel beta-sheet surrounded by alpha-helices, and two peripheral domains involved in covering the ATP in the active site (LID) and binding of the AMP (NMPbind). We initiated a long-term NMR study aiming to characterize the solution structure, binding mechanism and internal dynamics of the various domains. Using single (15N) and double-labeled (13C and 15N) samples and double- and triple-resonance NMR experiments we assigned 97% of the 1H, 13C and 15N backbone resonances, and proton and 13Cbeta resonances for more than 40% of the side chains in the free protein. Analysis of a 15N-labeled enzyme in complex with the bi-substrate analogue [P1,P5-bis(5'-adenosine)-pentaphosphate] (Ap5A) resulted in the assignment of 90% of the backbone 1H and 15N resonances and 42% of the side chain resonances. Based on short-range NOEs and 1H and 13C secondary chemical shifts, we identified the elements of secondary structure and the topology of the beta-strands in the unliganded form. The alpha-helices and the beta-strands of the parallel beta-sheet in solution have the same limits (+/- 1 residue) as those observed in the crystal. The first helix (alpha1) appears to have a frayed N-terminal side. Significant differences relative to the crystal were noticed in the LID domain, which in solution exhibits four antiparallel beta-strands. The secondary structure of the nucleoside-bound form, as deduced from intramolecular NOEs and the 1Halpha chemical shifts, is similar to that of the free enzyme. The largest chemical shift differences allowed us to map the regions of protein-ligand contacts. 1H/2H exchange experiments performed on free and Ap5A-bound enzymes showed a general decrease of the structural flexibility in the complex which is accompanied by a local increased flexibility on the N-side of the parallel beta-sheet.