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

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.

Crystal structure of adenylate kinase from an extremophilic archaeon Aeropyrum pernix with ATP and AMP

The crystal structure of an adenylate kinase from an extremophilic archaeon Aeropyrum pernix was determined in complex with full ligands, ATP-Mg2+ and AMP, at a resolution of 2.0 Å. The protein forms a trimer as found for other adenylate kinases from archaea. Interestingly, the reacting three atoms, two phosphorus and one oxygen atoms, were located almost in line, supporting the SN2 nucleophilic substitution reaction mechanism. Based on the crystal structure obtained, the reaction coordinate was estimated by the quantum mechanics calculations combined with molecular dynamics. It was found that the reaction undergoes two energy barriers; the steps for breaking the bond between the oxygen and γ-phosphorus atoms of ATP to produce a phosphoryl fragment and creating the bond between the phosphoryl fragment and the oxygen atom of the β-phosphate group of ADP. The reaction coordinate analysis also suggested the role of amino-acid residues for the catalysis of adenylate kinase.

Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations

Proteins exist as dynamic ensembles of molecules, implying that protein amino acid sequences evolved to code for both the ground state structure as well as the entire energy landscape of excited states. Accumulating theoretical and experimental evidence suggests that enzymes use such conformational fluctuations to facilitate allosteric processes important for substrate binding and possibly catalysis. This phenomenon can be clearly demonstrated in Escherichia coli adenylate kinase, where experimentally observed local unfolding of the LID subdomain, as opposed to a more commonly postulated rigid-body opening motion, is related to substrate binding. Because "entropy promoting" glycine mutations designed to increase specifically the local unfolding of the LID domain also affect substrate binding, changes in the excited energy landscape effectively tune the function of this enzyme without changing the ground state structure or the catalytic site. Thus, additional thermodynamic information, above and beyond the single folded structure of an enzyme-substrate complex, is likely required for a full and quantitative understanding of how enzymes work.

Rational modulation of conformational fluctuations in adenylate kinase reveals a local unfolding mechanism for allostery and functional adaptation in proteins

Elucidating the complex interplay between protein structure and dynamics is a prerequisite to an understanding of both function and adaptation in proteins. Unfortunately, it has been difficult to experimentally decouple these effects because it is challenging to rationally design mutations that will either affect the structure but not the dynamics, or that will affect the dynamics but not the structure. Here we adopt a mutation approach that is based on a thermal adaptation strategy observed in nature, and we use it to study the binding interaction of Escherichia coli adenylate kinase (AK). We rationally design several single-site, surface-exposed glycine mutations to selectively perturb the excited state conformational repertoire, leaving the ground-state X-ray crystallographic structure unaffected. The results not only demonstrate that the conformational ensemble of AK is significantly populated by a locally unfolded state that is depopulated upon binding, but also that the excited-state conformational ensemble can be manipulated through mutation, independent of perturbations of the ground-state structures. The implications of these results are twofold. First, they indicate that it is possible to rationally design dynamic allosteric mutations, which do not propagate through a pathway of structural distortions connecting the mutated and the functional sites. Secondly and equally as important, the results reveal a general strategy for thermal adaptation that allows enzymes to modulate binding affinity by controlling the amount of local unfolding in the native-state ensemble. These findings open new avenues for rational protein design and fundamentally illuminate the role of local unfolding in function and adaptation.

Substrate induced structural and dynamics changes in human phosphomevalonate kinase and implications for mechanism

Phosphomevalonate kinase (PMK) catalyzes an essential step in the mevalonate pathway, which is the only pathway for synthesis of isoprenoids and steroids in humans. PMK catalyzes transfer of the gamma-phosphate of ATP to mevalonate 5-phosphate (M5P) to form mevalonate 5-diphosphate. Bringing these phosphate groups in proximity to react is especially challenging, given the high negative charge density on the four phosphate groups in the active site. As such, conformational and dynamics changes needed to form the Michaelis complex are of mechanistic interest. Herein, we report the characterization of substrate induced changes (Mg-ADP, M5P, and the ternary complex) in PMK using NMR-based dynamics and chemical shift perturbation measurements. Mg-ADP and M5P K(d)'s were 6-60 microM in all complexes, consistent with there being little binding synergy. Binding of M5P causes the PMK structure to compress (tau(c) = 13.5 nsec), whereas subsequent binding of Mg-ADP opens the structure up (tau(c) = 15.6 nsec). The overall complex seems to stay very rigid on the psec-nsec timescale with an average NMR order parameter of S(2) approximately 0.88. Data are consistent with addition of M5P causing movement around a hinge region to permit domain closure, which would bring the M5P domain close to ATP to permit catalysis. Dynamics data identify potential hinge residues as H55 and R93, based on their low order parameters and their location in extended regions that connect the M5P and ATP domains in the PMK homology model. Likewise, D163 may be a hinge residue for the lid region that is homologous to the adenylate kinase lid, covering the "Walker-A" catalytic loop. Binding of ATP or ADP appears to cause similar conformational changes; however, these observations do not indicate an obvious role for gamma-phosphate binding interactions. Indeed, the role of gamma-phosphate interactions may be more subtle than suggested by ATP/ADP comparisons, because the conservative O to NH substitution in the beta-gamma bridge of ATP causes a dramatic decrease in affinity and induces few chemical shift perturbations. In terms of positioning of catalytic residues, binding of M5P induces a rigidification of Gly21 (adjacent to the catalytically important Lys22), although exchange broadening in the ternary complex suggests some motion on a slower timescale does still occur. Finally, the first nine residues of the N-terminus are highly disordered, suggesting that they may be part of a cleavable signal or regulatory peptide sequence.

Apo adenylate kinase encodes its holo form: a principal component and varimax analysis

Adenylate kinase undergoes large-scale motions of its LID and AMP-binding (AMPbd) domains when its apo, open form closes over its substrates, AMP and Mg2+-ATP. It may be an example of an enzyme that provides an ensemble of conformations in its apo state from which its substrates can select and bind to produce catalytically competent conformations. In this work, the fluctuations of the enzyme apo Escherichia coli adenylate kinase (AKE) are obtained with molecular dynamics. The resulting trajectory is analyzed with principal component analysis (PCA) that decomposes the atom motions into orthogonal modes ordered by their decreasing contributions to the total protein fluctuation. In apo AKE, a small set of the PCA modes describes the bulk of the fluctuations. Identification of the atom motions that are important contributors to these modes is improved with the use of a varimax rotation method that rotates the PCA modes to a new mode set that concentrates the atom contributions to a smaller set of atoms in these new modes. In this way, the nature of the important motions of the LID and AMPbd domains are clarified. The dominant PCA modes are used to investigate if apo AKE can fluctuate to conformations that are holo-like, even though the apo trajectory is mainly confined to a region around the initial apo structure. This is accomplished by expressing the difference between the protein coordinates, obtained from the holo and apo crystal structures, using as a basis the PCA modes from the apo AKE trajectory. The coherent motion described by a small set of the apo PCA modes is shown to be able to produce protein conformations that are quite similar to the holo conformation of the protein. In this sense, apo AKE does encode in its fluctuations information about holo-like conformations.

Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism

Large-scale conformational changes in proteins are often associated with the binding of a substrate. Because conformational changes may be related to the function of an enzyme, understanding the kinetics and energetics of these motions is very important. We have delineated the atomically detailed conformational transition pathway of the phosphotransferase enzyme adenylate kinase (AdK) in the absence and presence of an inhibitor. The computed free energy profiles associated with conformational transitions offer detailed mechanistic insights into, as well as kinetic information on, the ligand binding mechanism. Specifically, potential of mean force calculations reveal that in the ligand-free state, there is no significant barrier separating the open and closed conformations of AdK. The enzyme samples near closed conformations, even in the absence of its substrate. The ligand binding event occurs late, toward the closed state, and transforms the free energy landscape. In the ligand-bound state, the closed conformation is energetically most favored with a large barrier to opening. These results emphasize the underlying dynamic nature of the enzyme and indicate that the conformational transitions in AdK are more intricate than a mere two-state jump between the crystal-bound and -unbound states. Based on the existence of the multiple conformations of the enzyme in the open and closed states, a different viewpoint of ligand binding is presented. Our estimated activation energy barrier for the conformational transition is also in reasonable accord with the experimental findings.

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.

Molecular Dynamics of Apo-Adenylate Kinase: A Principal Component Analysis

Adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes large-scale motions of its Lid and AMP-binding domains when its open form closes over its substrates, AMP and Mg2+-ATP. The third domain, the Core, is relatively stable during closing. The resulting trajectory is analyzed with a principal component analysis method that decomposes the atom motions into modes ordered by their decreasing contributions to the total protein fluctuation. Simulations at 303 K (normal T) and 500 K (high T) reveal that at both temperatures the first three modes account for 70% of the total fluctuation. The residues that contribute the most to these three modes are concentrated in the Lid and AMP-binding domains. Analysis of the normal T modes indicates that the Lid and AMP-binding domains sample a broad distribution of conformations indicating that AKE is designed to provide its substrates with a large set of conformations. The high T results show that the Lid initially closes toward the Core. Subsequently, the Lid rotates to a new stable conformation that is different from what is observed in the substrate-bound AKE. These results are discussed in the context of experimental data that indicate that adenylate kinases do sample more than one conformational state in solution and that each of these conformational states undergoes substantial fluctuations. A pair of residues is suggested for labeling that would be useful for monitoring distance fluctuations by energy transfer experiments.

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.

Structural characterization of protein kinase A as a function of nucleotide binding. Hydrogen-deuterium exchange studies using matrix-assisted laser desorption ionization-time of flight mass spectrometry detection

Transient state kinetic studies indicate that substrate phosphorylation in protein kinase A is partially rate-limited by conformational changes, some of which may be associated with nucleotide binding (Shaffer, J., and Adams, J. A. (1999) Biochemistry 38, 12072-12079). To assess whether specific structural changes are associated with the binding of nucleotides, hydrogen-deuterium exchange experiments were performed on the enzyme in the absence and presence of ADP. Four regions of the protein are protected from exchange in the presence of ADP. Two regions encompass the catalytic and glycine-rich loops and are integral parts of the active site. Conversely, protection of probes in the C terminus is consistent with nucleotide-induced domain closure. One protected probe encompasses a portion of helix C, a secondary structural element that does not make any direct contacts with the nucleotide but has been reported to undergo segmental motion upon the activation of some protein kinases. The combined data suggest that binding of the nucleotide has distal structural effects that may include stabilizing the closed state of the enzyme and altering the position of a critical helix outside the active site. The latter represents the first evidence that the nucleotide alone can induce changes in helix C in solution.

Interim report on genomics of Escherichia coli

We present a summary of recent progress in understanding Escherichia coli K-12 gene and protein functions. New information has come both from classical biological experimentation and from using the analytical tools of functional genomics. The content of the E. coli genome can clearly be seen to contain elements acquired by horizontal transfer. Nevertheless, there is probably a large, stable core of >3500 genes that are shared among all E. coli strains. The gene-enzyme relationship is examined, and, in many cases, it exhibits complexity beyond a simple one-to-one relationship. Also, the E. coli genome can now be seen to contain many multiple enzymes that carry out the same or closely similar reactions. Some are similar in sequence and may share common ancestry; some are not. We discuss the concept of a minimal genome as being variable among organisms and obligatorily linked to their life styles and defined environmental conditions. We also address classification of functions of gene products and avenues of insight into the history of protein evolution.