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

Conformational dynamics of adenylate kinase in crystals

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

A pathway for assembling [4Fe-4S]2+ clusters in mitochondrial iron-sulfur protein biogenesis

During its late steps, the mitochondrial iron-sulfur cluster (ISC) assembly machinery leads to the formation of [4Fe-4S] clusters. In vivo studies revealed that several proteins are implicated in the biosynthesis and trafficking of [4Fe-4S] clusters in mitochondria. However, they do not provide a clear picture into how these proteins cooperate. Here, we showed that three late-acting components of the mitochondrial ISC assembly machinery (GLRX5, BOLA3, and NFU1) are part of a ISC assembly pathway leading to the synthesis of a [4Fe-4S]²⁺ cluster on NFU1. We showed that the [2Fe-2S]²⁺ GLRX5-BOLA3 complex transfers its cluster to monomeric apo NFU1 to form, in the presence of a reductant, a [4Fe-4S]²⁺ cluster bound to dimeric NFU1. The cluster formation on NFU1 does not occur with [2Fe-2S]²⁺ GLRX5, and thus, the [4Fe-4S] cluster assembly pathway is activated only in the presence of BOLA3. These results define NFU1 as an 'assembler' of [4Fe-4S] clusters, that is, a protein able of converting two [2Fe-2S]²⁺ clusters into a [4Fe-4S]²⁺ cluster. Finally, we found that the [4Fe-4S]²⁺ cluster bound to NFU1 has a coordination site which is easily accessible to sulfur-containing ligands, as is typically observed in metallochaperones. This finding supports a role for NFU1 in promoting rapid and controlled cluster-exchange reaction.

Enzymatic characteristics of two adenylate kinases, AdkA and AdkB, from Myxococcus xanthus

Adenylate kinase (Adk) plays a critical role in energy metabolism and adaptation of bacteria to environmental stresses. We have previously shown that Myxococcus xanthus expresses polyphosphate kinase 1 (Ppk1) that also has Adk activity in the absence of polyphosphates. In this study, we investigated the Adk activity of the other two M. xanthus enzymes, AdkA and AdkB. The activity of AdkA was increased by dithiothreitol (DTT), which also enhanced enzyme stability. Site-directed mutagenesis of three cysteine residues (C130, C150, and C153) present in the LID domain of AdkA revealed that the Adk activity and stability of C150S and C153S mutants were not affected by DTT addition, suggesting formation of a disulfide bond between C150 and C153 in AdkA. The Km of AdkA for AMP was 8 and 17 times lower than that for ADP and ATP, respectively. AdkB is a polyphosphate kinase 2 (Ppk2) homolog lacking the Ppk2 middle region and, consequently, Ppk activity. According to our analysis, AdkB also had Adk activity and its affinity for substrates was higher than that of AdkA. Thus, M. xanthus expresses three enzymes, AdkA, AdkB, and Ppk1, with Adk activity, which may function to support energy metabolism of the bacteria in different environmental conditions.

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of [Formula: see text]N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of [Formula: see text]N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2-4[Formula: see text] due to relatively low "compressibility" of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of [Formula: see text]/[Formula: see text] is sampled in a way similar to NUS of [Formula: see text]/[Formula: see text] domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional [Formula: see text]N HSQC spectra.

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

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

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Characterizing Active Site Conformational Heterogeneity along the Trajectory of an Enzymatic Phosphoryl Transfer Reaction

States along the phosphoryl transfer reaction catalyzed by the nucleoside monophosphate kinase UmpK were captured and changes in the conformational heterogeneity of conserved active site arginine side‐chains were quantified by NMR spin‐relaxation methods. In addition to apo and ligand‐bound UmpK, a transition state analog (TSA) complex was utilized to evaluate the extent to which active site conformational entropy contributes to the transition state free energy. The catalytically essential arginine side‐chain guanidino groups were found to be remarkably rigid in the TSA complex, indicating that the enzyme has evolved to restrict the conformational freedom along its reaction path over the energy landscape, which in turn allows the phosphoryl transfer to occur selectively by avoiding side reactions.

Reticular dysgenesis-associated AK2 protects hematopoietic stem and progenitor cell development from oxidative stress

Rissone et al. demonstrate that adenylate kinase AK2, an enzyme mutated in reticular dysgenesis (RD) in humans, prevents oxidative stress during hematopoiesis. Using a zebrafish model, as well as induced pluripotent stem cells derived from an RD patient, they find that AK2 deficiency affects hematopoietic stem and progenitor development with increased oxidative stress. Antioxidant treatment rescues the hematopoietic phenotypes.

Adenylate kinases (AKs) are phosphotransferases that regulate the cellular adenine nucleotide composition and play a critical role in the energy homeostasis of all tissues. The AK2 isoenzyme is expressed in the mitochondrial intermembrane space and is mutated in reticular dysgenesis (RD), a rare form of severe combined immunodeficiency (SCID) in humans. RD is characterized by a maturation arrest in the myeloid and lymphoid lineages, leading to early onset, recurrent, and overwhelming infections. To gain insight into the pathophysiology of RD, we studied the effects of AK2 deficiency using the zebrafish model and induced pluripotent stem cells (iPSCs) derived from fibroblasts of an RD patient. In zebrafish, Ak2 deficiency affected hematopoietic stem and progenitor cell (HSPC) development with increased oxidative stress and apoptosis. AK2-deficient iPSCs recapitulated the characteristic myeloid maturation arrest at the promyelocyte stage and demonstrated an increased AMP/ADP ratio, indicative of an energy-depleted adenine nucleotide profile. Antioxidant treatment rescued the hematopoietic phenotypes in vivo in ak2 mutant zebrafish and restored differentiation of AK2-deficient iPSCs into mature granulocytes. Our results link hematopoietic cell fate in AK2 deficiency to cellular energy depletion and increased oxidative stress. This points to the potential use of antioxidants as a supportive therapeutic modality for patients with RD.

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

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

On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase

We characterized the conformational change of adenylate kinase (AK) between open and closed forms by conducting five all-atom molecular-dynamics simulations, each of 100 ns duration. Different initial structures and substrate binding configurations were used to probe the pathways of AK conformational change in explicit solvent, and no bias potential was applied. A complete closed-to-open and a partial open-to-closed transition were observed, demonstrating the direct impact of substrate-mediated interactions on shifting protein conformation. The sampled configurations suggest two possible pathways for connecting the open and closed structures of AK, affirming the prediction made based on available x-ray structures and earlier works of coarse-grained modeling. The trajectories of the all-atom molecular-dynamics simulations revealed the complexity of protein dynamics and the coupling between different domains during conformational change. Calculations of solvent density and density fluctuations surrounding AK did not show prominent variation during the transition between closed and open forms. Finally, we characterized the effects of local unfolding of an important hinge near Pro(177) on the closed-to-open transition of AK and identified a novel mechanism by which hinge unfolding modulates protein conformational change. The local unfolding of Pro(177) hinge induces alternative tertiary contacts that stabilize the closed structure and prevent the opening transition.

Strategies for the thermodynamic characterization of linked binding/local folding reactions within the native state application to the LID domain of adenylate kinase from Escherichia coli

Conformational fluctuations in proteins have emerged as an important aspect of biological function, having been linked to processes ranging from molecular recognition and catalysis to allostery and signal transduction. In spite of the realization of their importance, however, the connections between fluctuations and function have largely been empirical, even when they have been quantitative. Part of the problem in understanding the role of fluctuations in function is the fact that the mere existence of fluctuations complicates the interpretation of classic mutagenesis approaches. Namely, mutagenesis, which is typically targeted to an internal position (to elicit an effect), will change the fluctuations as well as the structure of the native state. Decoupling these effects is essential to an unambiguous understanding of the role of fluctuations in function. Here, we use a mutation strategy that targets surface-exposed sites in flexible parts of the molecule for mutation to glycine. Such mutations leave the ground-state structure unaffected. As a result, we can assess the nature of the fluctuations, develop a quantitative model relating fluctuations to function (in this case, molecular recognition), and unambiguously resolve the probabilities of the fluctuating states. We show that when this approach is applied to Escherichia coli adenylate kinase (AK), unique thermodynamic and structural insights are obtained, even when classic mutagenesis approaches targeted to the same region yield ambiguous results.

Experimental evolution of adenylate kinase reveals contrasting strategies toward protein thermostability

Success in evolution depends critically upon the ability of organisms to adapt, a property that is also true for the proteins that contribute to the fitness of an organism. Successful protein evolution is enhanced by mutational pathways that generate a wide range of physicochemical mechanisms to adaptation. In an earlier study, we used a weak-link method to favor changes to an essential but maladapted protein, adenylate kinase (AK), within a microbial population. Six AK mutants (a single mutant followed by five double mutants) had success within the population, revealing a diverse range of adaptive strategies that included changes in nonpolar packing, protein folding dynamics, and formation of new hydrogen bonds and electrostatic networks. The first mutation, AK(BSUB) Q199R, was essential in defining the structural context that facilitated subsequent mutations as revealed by a considerable mutational epistasis and, in one case, a very strong dependence upon the order of mutations. Namely, whereas the single mutation AK(BSUB) G213E decreases protein stability by >25 degrees C, the same mutation in the background of AK(BSUB) Q199R increases stability by 3.4 degrees C, demonstrating that the order of mutations can play a critical role in favoring particular molecular pathways to adaptation. In turn, protein folding kinetics shows that four of the five AK(BSUB) double mutants utilize a strategy in which an increase in the folding rate accompanied by a decrease in the unfolding rate results in additional stability. However, one mutant exhibited a dramatic increase in the folding relative to a modest increase in the unfolding rate, suggesting a different adaptive strategy for thermostability. In all cases, an increase in the folding rates for the double mutants appears to be the preferred mechanism in conferring additional stability and may be an important aspect of protein evolution. The range of overlapping as well as contrasting strategies for success illustrates both the power and subtlety of adaptation at even the smallest unit of change, a single amino acid.

Heterogeneous path ensembles for conformational transitions in semi-atomistic models of adenylate kinase

We performed "weighted ensemble" path-sampling simulations of adenylate kinase, using several semi-atomistic protein models. The models have an all-atom backbone with various levels of residue interactions. The primary result is that full statistically rigorous path sampling required only a few weeks of single-processor computing time with these models, indicating the addition of further chemical detail should be readily feasible. Our semi-atomistic path ensembles are consistent with previous biophysical findings: the presence of two distinct pathways, identification of intermediates, and symmetry of forward and reverse pathways.

Structural characterization of Burkholderia pseudomallei adenylate kinase (Adk): profound asymmetry in the crystal structure of the 'open' state

In all organisms adenylate kinases (Adks) play a vital role in cellular energy metabolism and nucleic acid synthesis. Due to differences in catalytic properties between the Adks found in prokaryotes and in the cytoplasm of eukaryotes, there is interest in targeting this enzyme for new drug therapies against infectious bacterial agents. Here we report the 2.1A resolution crystal structure for the 220-residue Adk from Burkholderia pseudomallei (BpAdk), the etiological agent responsible for the infectious disease melioidosis. The general structure of apo BpAdk is similar to other Adk structures, composed of a CORE subdomain with peripheral ATP-binding (ATP(bd)) and LID subdomains. The two molecules in the asymmetric unit have significantly different conformations, with a backbone RMSD of 1.46 A. These two BpAdk conformations may represent 'open' Adk sub-states along the preferential pathway to the 'closed' substrate-bound state.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

The effect of Asp54 phosphorylation on the energetics and dynamics in the response regulator protein Spo0F studied by molecular dynamics

The response regulator protein Spo0F acts as an intermediate phospho-messenger in the signal transduction pathway that controls initiation of the differentiation process of sporculation in the bacterium Bacillus subtilis. The regulatory domain of Spo0F contains a triad of three conserved aspartate residues, whereof one aspartate (Asp54) is phosphorylated. Using molecular dynamics simulations, we have studied the changes in flexibility induced by phosphorylation and estimated the free energy cost of introducing a phosphate group at this position using alchemical free energy calculations. The deduced conformational flexibility compares well with experimental NMR results. We find that the apo-conformation of the protein explores a rough energy landscape resulting in a broad population of conformational substates. Phosphorylation of Spo0F reduces protein flexibility, and in particular, the so-called anchor and recognition regions exhibit lower mobility relative to the apo-conformation. Phosphorylation of Asp54 (P-Asp54), in which the apo-structure coordinates to the magnesium ion, results in extension of the sidechain, and depending on which carboxylate oxygen is phosphorylated, distinct interactions between P-Asp54 and magnesium ion as well as residues in its proximity are established. However, phosphorylation does not affect the coordination number of the magnesium ion yielding, within the statistical uncertainties, the same free energy change.

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.

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.

Direct Mg(2+) binding activates adenylate kinase from Escherichia coli

We report evidence that adenylate kinase (AK) from Escherichia coli can be activated by the direct binding of a magnesium ion to the enzyme, in addition to ATP-complexed Mg(2+). By systematically varying the concentrations of AMP, ATP, and magnesium in kinetic experiments, we found that the apparent substrate inhibition of AK, formerly attributed to AMP, was suppressed at low magnesium concentrations and enhanced at high magnesium concentrations. This previously unreported magnesium dependence can be accounted for by a modified random bi-bi model in which Mg(2+) can bind to AK directly prior to AMP binding. A new kinetic model is proposed to replace the conventional random bi-bi mechanism with substrate inhibition and is able to describe the kinetic data over a physiologically relevant range of magnesium concentrations. According to this model, the magnesium-activated AK exhibits a 23- +/- 3-fold increase in its forward reaction rate compared with the unactivated form. The findings imply that Mg(2+) could be an important affecter in the energy signaling network in cells.

Small- and large-scale conformational changes of adenylate kinase: a molecular dynamics study of the subdomain motion and mechanics

Adenylate kinase, an enzyme that catalyzes the phosphoryl transfer between ATP and AMP, can interconvert between the open and catalytically potent (closed) forms even without binding ligands. Several aspects of the enzyme elasticity and internal dynamics are analyzed here by atomistic molecular dynamics simulations covering a total time span of 100 ns. This duration is sufficiently long to reveal a partial conversion of the enzyme that proceeds through jumps between structurally different substates. The intra- and intersubstates contributions to the enzyme's structural fluctuations are analyzed and compared both in magnitude and directionality. It is found that, despite the structural heterogeneity of the visited conformers, the generalized directions accounting for conformational fluctuations within and across the substates are mutually consistent and can be described by a limited set of collective modes. The functional-oriented nature of the consensus modes is suggested by their good overlap with the deformation vector bridging the open and closed crystal structures. The consistency of adenylate kinase's internal dynamics over timescales wide enough to capture intra- and intersubstates fluctuations adds elements in favor of the recent proposal that the free (apo) enzyme possesses an innate ability to sustain the open/close conformational changes.

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.

Structure of the yeast SR protein Npl3 and Interaction with mRNA 3'-end processing signals

Yeast Npl3 is homologous to SR proteins in higher eukaryotes, a family of RNA-binding proteins that have multiple essential roles in RNA metabolism. This protein competes with 3'-end processing factors for binding to the nascent RNA, protecting the transcript from premature termination and coordinating transcription termination and the packaging of the fully processed transcript for export. The NMR structure of its RNA-binding domain shows two unusually compact RNA recognition motifs (RRMs), and identifies the RNA recognition surface in Npl3. Biochemical and NMR studies identify a class of G+U-rich RNA sequences with high specificity for this protein. The protein binds to RNA and forms a single globular structure, but the two RRMs of Npl3 are not equivalent, with the second domain forming much stronger interactions with G+U-rich RNA sequences that occur independently of the interaction of the first RRM. The specific binding to G+U-rich RNAs observed for the two RRMs of Npl3 is masked in the full-length protein by a much stronger but non-sequence-specific RNA-binding activity residing outside of its RRMs. The preference of Npl3 for G+U-rich sequences supports the model for its function in regulating recognition of 3'-end processing sites through competition with the Rna15 (yeast analog of human CstF-64 protein) subunit of the processing complex.

Associative mechanism for phosphoryl transfer: a molecular dynamics simulation of Escherichia coli adenylate kinase complexed with its substrates

The ternary complex of Escherichia coli adenylate kinase (ECAK) with its substrates adenosine monophosphate (AMP) and Mg-ATP, which catalyzes the reversible transfer of a phosphoryl group between adenosine triphosphate (ATP) and AMP, was studied using molecular dynamics. The starting structure for the simulation was assembled from the crystal structures of ECAK complexed with the bisubstrate analog diadenosine pentaphosphate (AP(5)A) and of Bacillus stearothermophilus adenylate kinase complexed with AP(5)A, Mg(2+), and 4 coordinated water molecules, and by deleting 1 phosphate group from AP(5)A. The interactions of ECAK residues with the various moieties of ATP and AMP were compared to those inferred from NMR, X-ray crystallography, site-directed mutagenesis, and enzyme kinetic studies. The simulation supports the hypothesis that hydrogen bonds between AMP's adenine and the protein are at the origin of the high nucleoside monophosphate (NMP) specificity of AK. The ATP adenine and ribose moieties are only loosely bound to the protein, while the ATP phosphates are strongly bound to surrounding residues. The coordination sphere of Mg(2+), consisting of 4 waters and oxygens of the ATP beta- and gamma-phosphates, stays approximately octahedral during the simulation. The important role of the conserved Lys13 in the P loop in stabilizing the active site by bridging the ATP and AMP phosphates is evident. The influence of Mg(2+), of its coordination waters, and of surrounding charged residues in maintaining the geometry and distances of the AMP alpha-phosphate and ATP beta- and gamma-phosphates is sufficient to support an associative reaction mechanism for phosphoryl transfer.

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.

Mapping of epitopes for autoantibodies to the type 1 diabetes autoantigen IA-2 by peptide phage display and molecular modeling: overlap of antibody and T cell determinants

IA-2 is a major target of autoimmunity in type 1 diabetes. IA-2 responsive T cells recognize determinants within regions represented by amino acids 787-817 and 841-869 of the molecule. Epitopes for IA-2 autoantibodies are largely conformational and not well defined. In this study, we used peptide phage display and homology modeling to characterize the epitope of a monoclonal IA-2 Ab (96/3) from a human type 1 diabetic patient. This Ab competes for IA-2 binding with Abs from the majority of patients with type 1 diabetes and therefore binds a region close to common autoantibody epitopes. Alignment of peptides obtained after screening phage-displayed peptide libraries with purified 96/3 identified a consensus binding sequence of Asn-x-Glu-x-x-(aromatic)-x-x-Gly. The predicted surface on a three-dimensional homology model of the tyrosine phosphatase domain of IA-2 was analyzed for clusters of Asn, Glu, and aromatic residues and amino acids contributing to the epitope investigated using site-directed mutagenesis. Mutation of each of amino acids Asn(858), Glu(836), and Trp(799) reduced 96/3 Ab binding by >45%. Mutations of these residues also inhibited binding of serum autoantibodies from IA-2 Ab-positive type 1 diabetic patients. This study identifies a region commonly recognized by autoantibodies in type 1 diabetes that overlaps with dominant T cell determinants.

Structures of thermophilic and mesophilic adenylate kinases from the genus Methanococcus

The crystal structures of adenylate kinases from the thermophile Methanococcus thermolithotrophicus and the mesophile Methanococcus voltae have been solved to resolutions of 2.8A and 2.5A, respectively. The structures of the enzymes are similar to that of the adenylate kinase from archaeal Sulfolobus acidocaldarius in many respects such as the extended central beta-sheets, the short LID domain, and the trimeric state. The analysis of unligated and AMP-bound subunits of M.voltae suggests that movements of two mobile domains are not independent of each other. The methanococcal structures are examined with respect to their lack of the "invariant" Lys residue within the phosphate-binding loop, and two Arg residues in the LID domain are proposed as substituting residues based on their conservation among archaeal adenylate kinases and mobility within the structures. Since S.acidocaldarius adenylate kinase has the invariant Lys residue as well as the two Arg residues, its phosphate-binding loop is examined and compared with those of other adenylate kinases. On the basis of the comparison and other available biochemical data, the unusual conformation of the Lys residue in S.acidocaldarius adenylate kinase is explained. Despite possessing 78% sequence identity, the methanococcal enzymes exhibit significantly different thermal stabilities. To study the determinants of thermostability, several structural features including salt-links, hydrogen bonds, packing density, surface to volume ratio and buried surface area are compared between the enzymes. From their difference in apolar buried surface area, hydrophobic interaction is proposed to be a basis for the disparate thermostabilities, and the corresponding free energy difference is also estimated. Results of previous mutational studies are interpreted in terms of the crystal structures, and support the importance of hydrophobic interactions in thermostability.

Essential motions in a fungal lipase with bound substrate, covalently attached inhibitor and product

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics. This 'essential dynamics' analysis reveals different modes of response along the pathway suggesting that binding, catalysis and product clearance occur along different energy surfaces. Motions in the enzyme with a covalently attached ligand are more complex and occur along several eigenvectors. The magnitudes of the fluctuations in these individual subspaces are significantly smaller than those observed for the substrate and product molecules, indicating that the energy surface is shallow and that a relatively large number of conformational substates are accessible. On the other hand, substrate binding and product release occur at distinct modes of the protein flexibility suggesting that these processes occur along rough energy surfaces with only a few minima. Detailed energetic analyses along the trajectories indicated that in all cases binding is dominated by van der Waals interactions. The carboxylate form of the product is stabilized by a tight hydrogen bond network involving in particular Ser82, which may be a potential cause of product inhibition. Considerations such as these should aid the understanding of mechanisms of substrate, inhibitor or product recognition and could become of importance in the design of new substrates or inhibitors for enzymes.

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.

Conformational and functional significance of residue proline 17 in chicken muscle adenylate kinase

The effect of mutation proline 17 on the multiple conformations and catalytic function in chicken muscle adenylate kinase (AK) has been studied. The substitution of proline 17 with glycine or valine altered the distribution of multiple conformations. Compared with the wild-type enzyme, the P17G and P17V mutants contained decreased fraction of minor conformer from 18% to 9% and 11%, respectively. Due to the mutation, the enzyme showed lower secondary structural content, poorer affinity to substrates or substrate analogues, and reduced catalytic efficiency. The results revealed the significance of proline 17 in the conformation and function of AK.

Functional dynamics in the active site of the ribonuclease binase

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-micros time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a beta-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum k(cat) (1,500 s(-1)) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

A structural mode-coupling approach to 15N NMR relaxation in proteins

The two-body Slowly Relaxing Local Structure (SRLS) model was applied to (15)N NMR spin relaxation in proteins and compared with the commonly used original and extended model-free (MF) approaches. In MF, the dynamic modes are assumed to be decoupled, local ordering at the N-H sites is represented by generalized order parameters, and internal motions are described by effective correlation times. SRLS accounts for dynamical coupling between the global diffusion of the protein and the internal motion of the N-H bond vector. The local ordering associated with the coupling potential and the internal N-H diffusion are tensors with orientations that may be tilted relative to the global diffusion and magnetic frames. SRLS generates spectral density functions that differ from the MF formulas. The MF spectral densities can be regarded as limiting cases of the SRLS spectral density. SRLS-based model-fitting and model-selection schemes similar to the currently used MF-based ones were devised, and a correspondence between analogous SRLS and model-free parameters was established. It was found that experimental NMR data are sensitive to the presence of mixed modes. Our results showed that MF can significantly overestimate order parameters and underestimate local motion correlation times in proteins. The extent of these digressions in the derived microdynamic parameters is estimated in the various parameter ranges, and correlated with the time scale separation between local and global motions. The SRLS-based analysis was tested extensively on (15)N relaxation data from several isotropically tumbling proteins. The results of SRLS-based fitting are illustrated with RNase H from E. coli, a protein extensively studied previously with MF.