The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase. A highly stable salt bridge D118-K136 at the LID-CORE interface, contributing substantially to the total nonbonded LID-CORE interactions, was identified as a major factor that stabilizes the open conformation. Alternative transition pathways, such as AMPbd opening following LID opening, seem unlikely, as full transition events were not observed along this pathway. The simulation data indicate a high enthalpic penalty, possibly obstructing transitions along this route.     

Conformational transitions in adenylate kinase. Allosteric communication reduces misligation

Large conformational changes in the LID and NMP domains of adenylate kinase (AKE) are known to be key to ligand binding and catalysis, yet the order of binding events and domain motion is not well understood. Combining the multiple available structures for AKE with the energy landscape theory for protein folding, a theoretical model was developed for allostery, order of binding events, and efficient catalysis. Coarse-grained models and nonlinear normal mode analysis were used to infer that intrinsic structural fluctuations dominate LID motion, whereas ligand-protein interactions and cracking (local unfolding) are more important during NMP motion. In addition, LID-NMP domain interactions are indispensable for efficient catalysis. LID domain motion precedes NMP domain motion, during both opening and closing. These findings provide a mechanistic explanation for the observed 1:1:1 correspondence between LID domain closure, NMP domain closure, and substrate turnover. This catalytic cycle has likely evolved to reduce misligation, and thus inhibition, of AKE. The separation of allosteric motion into intrinsic structural fluctuations and ligand-induced contributions can be generalized to further our understanding of allosteric transitions in other proteins.     

Many Local Motions Cooperate to Produce the Adenylate Kinase Conformational Transition

Conformational transitions are functionally important in many proteins. In the enzyme adenylate kinase (AK), two small domains (LID and NMP) close over the larger CORE domain; the reverse (opening) motion limits the rate of catalytic turnover. Here, using double-well Gō simulations of Escherichia coli AK, we elaborate on previous investigations of the AK transition mechanism by characterizing the contributions of rigid-body (Cartesian), backbone dihedral, and contact motions to transition-state (TS) properties. In addition, we compare an apo simulation to a pseudo-ligand-bound simulation to reveal insights into allostery. In Cartesian space, LID closure precedes NMP closure in the bound simulation, consistent with prior coarse-grained models of the AK transition. However, NMP-first closure is preferred in the apo simulation. In backbone dihedral space, we find that, as expected, backbone fluctuations are reduced in the O/C transition in parts of all three domains. Among these “quenching” residues, most in the CORE, especially residues 11–13, are rigidified in the TS of the bound simulation, while residues 42–44 in the NMP are flexible in the TS. In contact space, in both apo and bound simulations, one nucleus of closed-state contacts includes parts of the NMP and CORE; CORE–LID contacts are absent in the TS of the apo simulation but formed in the TS of the bound simulation. From these results, we predict mutations that will perturb the opening and/or closing transition rates by changing the entropy of dihedrals and/or the enthalpy of contacts. Furthermore, regarding allostery, the fully closed structure is populated in the apo simulation, but our contact results imply that ligand binding shifts the preferred O/C transition pathway, thus precluding a simple conformational selection mechanism. Finally, the analytical approach and the insights derived from this work may inform the rational design of flexibility and allostery in proteins.     

Conformational Dynamics of a Ligand-Free Adenylate Kinase

Adenylate kinase (AdK) is a phosphoryl-transfer enzyme with important physiological functions. Based on a ligand-free open structure and a ligand-bound closed structure solved by crystallography, here we use molecular dynamics simulations to examine the stability and dynamics of AdK conformations in the absence of ligands. We first perform multiple simulations starting from the open or the closed structure, and observe their free evolutions during a simulation time of 100 or 200 nanoseconds. In all seven simulations starting from the open structure, AdK remained stable near the initial conformation. The eight simulations initiated from the closed structure, in contrast, exhibited large variation in the subsequent evolutions, with most (seven) undergoing large-scale spontaneous conformational changes and approaching or reaching the open state. To characterize the thermodynamics of the transition, we propose and apply a new sampling method that employs a series of restrained simulations to calculate a one-dimensional free energy along a curved pathway in the high-dimensional conformational space. Our calculated free energy profile features a single minimum at the open conformation, and indicates that the closed state, with a high (∼13 kcal/mol) free energy, is not metastable, consistent with the observed behaviors of the unrestrained simulations. Collectively, our simulations suggest that it is energetically unfavorable for the ligand-free AdK to access the closed conformation, and imply that ligand binding may precede the closure of the enzyme.     

Crystal structure of human adenylate kinase 4 (L171P) suggests the role of hinge region in protein domain motion

It is well known that motion of LID and NMP-binding (NMP_bind) domains in adenylate kinase (AK) is important in ligand binding and catalysis. However, the nature of such domain motions is poorly characterized. One of the critical hinge regions is hinge IV, which connects the CORE and LID domains. In addition, the hinge IV contains a strictly conserved residue, L171, in the AK family. To investigate the role of hinge IV, crystal structure of human adenylate kinase 4 (AK4) L171P mutant was determined. This mutation dramatically changes the orientation of the LID domain, which could be described as a novel twisted-and-closed conformation in contrast to the open and closed conformations in other AKs. This mutant provides a new example of domain motions in AK family.     

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7–15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4–3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its “parents” up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5–2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.     

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, , 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.     

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.     

Atomistic View of the Conformational Activation of Src Kinase Using the String Method with Swarms-of-Trajectories

The inactive-to-active conformational transition of the catalytic domain of human c-Src tyrosine kinase is characterized using the string method with swarms-of-trajectories with all-atom explicit solvent molecular dynamics simulations. The activation process occurs in two main steps in which the activation loop (A-loop) opens first, followed by the rotation of the αC helix. The computed potential of mean force energy along the activation pathway displays a local minimum, which allows the identification of an intermediate state. These results show that the string method with swarms-of-trajectories is an effective technique to characterize complex and slow conformational transitions in large biomolecular systems.     

New insights into FAK signaling and localization based on detection of a FAT domain folding intermediate

Mounting evidence suggests that the focal adhesion targeting (FAT) domain, an antiparallel four-helix bundle, exists in alternative conformations that may modulate phosphorylation, ligand binding, and the subcellular localization of focal adhesion kinase (FAK). In order to characterize the conformational dynamics of the FAT domain, we have developed a novel method for reconstructing the folding pathway of the FAT domain by using discrete molecular dynamics (DMD) simulations, with free energy constraints derived from NMR hydrogen exchange data. The DMD simulations detect a folding intermediate, in which a cooperative unfolding event causes helix 1 to lose helical character while separating from the helix bundle. The conformational dynamic features of helix 1 in the intermediate state of the FAT domain are likely to facilitate Y926 phosphorylation, yet interfere with paxillin binding. The presence of this intermediate state in vivo may promote FAK signaling via the ERK/MAPK pathway and by release of FAK from focal adhesions.     

How maltose influences structural changes to bind to maltose‐binding protein: Results from umbrella sampling simulation

A well‐studied periplasmic‐binding protein involved in the abstraction of maltose is maltose‐binding protein (MBP), which undergoes a ligand‐induced conformational transition from an open (ligand‐free) to a closed (ligand‐bound) state. Umbrella sampling simulations have been us to estimate the free energy of binding of maltose to MBP and to trace the potential of mean force of the unbinding event using the center‐of‐mass distance between the protein and ligand as the reaction coordinate. The free energy thus obtained compares nicely with the experimentally measured value justifying our theoretical basis. Measurement of the domain angle (N‐terminal‐domain – hinge – C‐terminal‐domain) along the unbinding pathway established the existence of three different states. Starting from a closed state, the protein shifts to an open conformation during the initial unbinding event of the ligand then resides in a semi‐open conformation and later resides predominantly in an open‐state. These transitions along the ligand unbinding pathway have been captured in greater depth using principal component analysis. It is proposed that in mixed‐model, both conformational selection and an induced‐fit mechanism combine to the ligand recognition process in MBP. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

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.     

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

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

Stability domains, substrate-induced conformational changes, and hinge-bending motions in a psychrophilic phosphoglycerate kinase. A microcalorimetric study

The cold-active phosphoglycerate kinase from the Antarctic bacterium Pseudomonas sp. TACII18 exhibits two distinct stability domains in the free, open conformation. It is shown that these stability domains do not match the structural N- and C-domains as the heat-stable domain corresponds to about 80 residues of the C-domain, including the nucleotide binding site, whereas the remaining of the protein contributes to the main heat-labile domain. This was demonstrated by spectroscopic and microcalorimetric analyses of the native enzyme, of its mutants, and of the isolated recombinant structural domains. It is proposed that the heat-stable domain provides a compact structure improving the binding affinity of the nucleotide, therefore increasing the catalytic efficiency at low temperatures. Upon substrate binding, the enzyme adopts a uniformly more stable closed conformation. Substrate-induced stability changes suggest that the free energy of ligand binding is converted into an increased conformational stability used to drive the hinge-bending motions and domain closure.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

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.     

Conformational dynamics and thermodynamics of protein-ligand binding studied by NMR relaxation

Protein conformational dynamics can be critical for ligand binding in two ways that relate to kinetics and thermodynamics respectively. First, conformational transitions between different substates can control access to the binding site (kinetics). Secondly, differences between free and ligand-bound states in their conformational fluctuations contribute to the entropy of ligand binding (thermodynamics). In the present paper, I focus on the second topic, summarizing our recent results on the role of conformational entropy in ligand binding to Gal3C (the carbohydrate-recognition domain of galectin-3). NMR relaxation experiments provide a unique probe of conformational entropy by characterizing bond-vector fluctuations at atomic resolution. By monitoring differences between the free and ligand-bound states in their backbone and side chain order parameters, we have estimated the contributions from conformational entropy to the free energy of binding. Overall, the conformational entropy of Gal3C increases upon ligand binding, thereby contributing favourably to the binding affinity. Comparisons with the results from isothermal titration calorimetry indicate that the conformational entropy is comparable in magnitude to the enthalpy of binding. Furthermore, there are significant differences in the dynamic response to binding of different ligands, despite the fact that the protein structure is virtually identical in the different protein-ligand complexes. Thus both affinity and specificity of ligand binding to Gal3C appear to depend in part on subtle differences in the conformational fluctuations that reflect the complex interplay between structure, dynamics and ligand interactions.     

Interconversion of functional motions between mesophilic and thermophilic adenylate kinases

Dynamic properties are functionally important in many proteins, including the enzyme adenylate kinase (AK), for which the open/closed transition limits the rate of catalytic turnover. Here, we compare our previously published coarse-grained (double-well Gō) simulation of mesophilic AK from E. coli (AKmeso) to simulations of thermophilic AK from Aquifex aeolicus (AKthermo). In AKthermo, as with AKmeso, the LID domain prefers to close before the NMP domain in the presence of ligand, but LID rigid-body flexibility in the open (O) ensemble decreases significantly. Backbone foldedness in O and/or transition state (TS) ensembles increases significantly relative to AKmeso in some interdomain backbone hinges and within LID. In contact space, the TS of AKthermo has fewer contacts at the CORE-LID interface but a stronger contact network surrounding the CORE-NMP interface than the TS of AKmeso. A "heated" simulation of AKthermo at 375K slightly increases LID rigid-body flexibility in accordance with the "corresponding states" hypothesis. Furthermore, while computational mutation of 7 prolines in AKthermo to their AKmeso counterparts produces similar small perturbations, mutation of these sites, especially positions 8 and 155, to glycine is required to achieve LID rigid-body flexibility and hinge flexibilities comparable to AKmeso. Mutating the 7 sites to proline in AKmeso reduces some hinges' flexibilities, especially hinge 2, but does not reduce LID rigid-body flexibility, suggesting that these two types of motion are decoupled in AKmeso. In conclusion, our results suggest that hinge flexibility and global functional motions alike are correlated with but not exclusively determined by the hinge residues. This mutational framework can inform the rational design of functionally important flexibility and allostery in other proteins toward engineering novel biochemical pathways.     

Event detection and sub‐state discovery from biomolecular simulations using higher‐order statistics: Application to enzyme adenylate kinase

Biomolecular simulations at millisecond and longer time‐scales can provide vital insights into functional mechanisms. Because post‐simulation analyses of such large trajectory datasets can be a limiting factor in obtaining biological insights, there is an emerging need to identify key dynamical events and relating these events to the biological function online, that is, as simulations are progressing. Recently, we have introduced a novel computational technique, quasi‐anharmonic analysis (QAA) (Ramanathan et al., PLoS One 2011;6:e15827), for partitioning the conformational landscape into a hierarchy of functionally relevant sub‐states. The unique capabilities of QAA are enabled by exploiting anharmonicity in the form of fourth‐order statistics for characterizing atomic fluctuations. In this article, we extend QAA for analyzing long time‐scale simulations online. In particular, we present HOST4MD—a higher‐order statistical toolbox for molecular dynamics simulations, which (1) identifies key dynamical events as simulations are in progress, (2) explores potential sub‐states, and (3) identifies conformational transitions that enable the protein to access those sub‐states. We demonstrate HOST4MD on microsecond timescale simulations of the enzyme adenylate kinase in its apo state. HOST4MD identifies several conformational events in these simulations, revealing how the intrinsic coupling between the three subdomains (LID, CORE, and NMP) changes during the simulations. Further, it also identifies an inherent asymmetry in the opening/closing of the two binding sites. We anticipate that HOST4MD will provide a powerful and extensible framework for detecting biophysically relevant conformational coordinates from long time‐scale simulations. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

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, , of 0.84 and correlation times tau(f) = 5-130 ps. Most of the AKeco*AP(5)A backbone features = 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.