Identification of specific interactions that drive ligand-induced closure in five enzymes with classic domain movements

In order to better understand ligand-induced closure in domain enzymes, open unliganded X-ray structures and closed liganded X-ray structures have been studied in five enzymes: adenylate kinase, aspartate aminotransferase, citrate synthase, liver alcohol dehydrogenase, and the catalytic subunit of cAMP-dependent protein kinase. A sequential model of ligand binding and domain closure was used to test the hypothesis that the ligand actively drives closure from an open conformation. The analysis supports the assumption that each enzyme has a dedicated binding domain to which the ligand binds first and a closing domain. In every case, a small number of residues are identified to interact with the ligand to initiate and drive domain closure. In all cases except adenylate kinase, the backbone of residues located in an interdomain-bending region (hinge site) is identified to interact with the ligand to aid in driving closure. In adenylate kinase, the side-chain of a residue located directly adjacent to a bending region drives closure. It is thought that by binding near a hinge site the ligand is able to get within interaction range of residues when the enzyme is in the open conformation. Interdomain bending regions not involved in inducing closure are involved in control, helping to determine the location of the hinge axis. Similarities have been discovered between aspartate aminotransferase and citrate synthase that only come to light in the context of their dynamical behaviour in response to binding their substrate. Similarity also exists between liver alcohol dehydrogenase and cAMP-dependent protein kinase whereby groups on NAD and ATP, respectively, mimic the backbone of a single amino acid residue in a process where a three residue segment located at the terminus of a beta-sheet, moves to form hydrogen bonds with the mimic that resemble those found in a parallel beta-sheet. This interaction helps to drive domain closure in a process that has analogy to protein folding.     

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.     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.     

Investigating the accessibility of the closed domain conformation of citrate synthase using essential dynamics sampling

A molecular dynamics study of pig heart citrate synthase is presented that aims to directly address the question of whether, for this enzyme, the ligand-induced closed domain conformation is accessible to the open unliganded enzyme. The approach utilises the technique of essential dynamics sampling, which is used in two modes. In exploring mode, the enzyme is encouraged to explore domain conformations it might not normally sample in free molecular dynamics simulation. In targeting mode, the enzyme is encouraged to adopt the domain conformation of a target structure. Using both modes extensively, it has been found that when the enzyme is prepared from a crystallographic open-domain structure and is in the unliganded state, it is unable to adopt the crystallographic closed-domain conformation of the liganded enzyme. Likewise, when the enzyme is prepared from the crystallographic closed liganded conformation with the ligands removed, it is unable to adopt the crystallographic open domain conformation. Structural investigations point to a common structural difference that is the source of this energy barrier; namely, the shift of alpha-helix 328-341 along its own axis relative to the large domain. Without this shift, the domains are unable to close or open fully. The charged substrate, oxaloacetate, binds near the base of this helix in the large domain and the interaction of Arg329 at the base of the helix with oxaloacetate is one that is consistent with the shift of this helix in going from the crystallographic open to closed structure. Therefore, the results suggest that without the substrate the enzyme remains in a partially open conformation ready to receive the substrate. In this way, the efficiency of the enzyme should be increased over one that is closed part of the time, with its binding site inaccessible to the substrate.     

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.     

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

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

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.     

Crystal structure of unligated guanylate kinase from yeast reveals GMP-induced conformational changes

The crystal structure of guanylate kinase (GK) from yeast (Saccharomyces cerevisiae) with a non-acetylated N terminus has been determined in its unligated form (apo-GK) as well as in complex with GMP (GK.GMP). The structure of apo-GK was solved with multiwavelength anomalous diffraction data and refined to an R-factor of 0.164 (R(free)=0.199) at 2.3 A resolution. The structure of GK.GMP was determined using the crystal structure of GK with an acetylated N terminus as the search model and refined to an R-factor of 0.156 (R(free)=0.245) at 1.9 A. GK belongs to the family of nucleoside monophosphate (NMP) kinases and catalyzes the reversible phosphoryl transfer from ATP to GMP. Like other NMP kinases, GK consists of three dynamic domains: the CORE, LID, and NMP-binding domains. Dramatic movements of the GMP-binding domain and smaller but significant movements of the LID domain have been revealed by comparing the structures of apo-GK and GK.GMP. apo-GK has a much more open conformation than the GK.GMP complex. Systematic analysis of the domain movements using the program DynDom shows that the large movements of the GMP-binding domain involve a rotation around an effective hinge axis approximately parallel with helix 3, which connects the GMP-binding and CORE domains. The C-terminal portion of helix 3, which connects to the CORE domain, has strikingly higher temperature factors in GK.GMP than in apo-GK, indicating that these residues become more mobile upon GMP binding. The results suggest that helix 3 plays an important role in domain movement. Unlike the GMP-binding domain, which moves toward the active center of the enzyme upon GMP binding, the LID domain moves away from the active center and makes the presumed ATP-binding site more open. Therefore, the LID domain movement may facilitate the binding of MgATP. The structure of the recombinant GK.GMP complex superimposes very well with that of the native GK.GMP complex, indicating that N-terminal acetylation does not have significant impact on the three-dimensional structure of GK.     

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.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

Mechanisms of domain closure in proteins

Certain enzymes respond to the binding of substrates and coenzymes by the closure of an active site that lies in a cleft between two domains. We have examined the mechanism of the domain closure in citrate synthase, for which atomic co-ordinates are available for "open" and "closed" forms. We show that the mechanism of domain closure involves small shifts and rotations of packed helices within the two domains and at their interface. Large motions of distant segments of the structure are the cumulative effect of the small relative shifts in intervening pairs of packed segments. These shifts are accommodated not by changes in packing but rather by small conformational changes in side-chains. We call this the helix interface shear mechanism of domain closure. The relative movements of packed helices follow the principles suggested by our recent study of insulin. This mechanism of domain closure is quite different from the hinge mechanisms that allow the rigid body movements of domains in immunoglobulins. The large interface between the domains of citrate synthase precludes a simple hinge mechanism for its conformational change. The helix interface shear mechanism of conformational change occurs in other enzymes that contain extensive domain-domain interfaces.     

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.     

Fine structure of conformational ensembles in adenylate kinase

Adenylate kinase (ADK) catalyzes the reversible Mg²⁺‐dependent phosphoryl transfer reaction Mg²⁺+2ADP ?Mg²⁺+ATP + AMP in essential cellular systems. This reaction is a major player in cellular energy homeostasis and the isoform network of ADK plays an important role in AMP metabolic signaling circuits. ADK has 3 domains, the LID, NMP, and CORE domains, that undergo large conformational rearrangements during ADK's catalytic cycle. In spite of extensive experimental and computational studies, details of the conformational pathway from open to closed forms remain uncertain. In this paper we explore this pathway using coarse‐grained molecular dynamics (MD) trajectories of ADK calculated by GROMACS using a SMOG model and classify the conformations within the resultant trajectories by K‐means clustering. ADK conformations segregate naturally into open; intermediate; and closed forms with long‐term residence in the intermediate state. Structural clustering divides the intermediate conformation into 3 sub‐states that are distinguished from one another on the basis of differences in both structure and dynamics. These distinctions are defined on the basis of a number of different metrics including radius of gyration, dihedral angle fluctuation, and fluctuations of interatomic pair distances. Furthermore, differences in the sub‐states appear to correspond to the distinct ways each sub‐state contributes to the molecular mechanism of catalysis: One sub‐state acts as a gate‐way to the open conformation; one sub‐state a gate‐way to the closed conformation. A third intermediate sub‐state appears to represent a metastable off‐pathway structure that is nevertheless frequently visited during the passage from open to closed state.     

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

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

Modulation of functionally significant conformational equilibria in adenylate kinase by high concentrations of trimethylamine oxide attributed to volume exclusion

The effect of an inert small molecule osmolyte, trimethyl amine N-oxide (TMAO), upon the conformational equilibria of Escherichia coli adenylate kinase was studied using time-resolved FRET. The relative populations of open and closed clefts between the LID and the CORE domains were measured as functions of the concentrations of the substrate ATP over the concentration range 0–18 mM and TMAO over the concentration range 0–4 M. A model was constructed according to which the enzyme exists in equilibrium among four conformational states, corresponding to combinations of open and closed conformations of the LID-CORE and AMP-CORE clefts. ATP is assumed to bind only to those conformations with the closed LID-CORE cleft, and TMAO is assumed to be differentially excluded as a hard spherical particle from each of the four conformations in accordance with calculations based upon x-ray crystallographic structures. This model was found to describe quantitatively the dependence of the fraction of the closed LID-CORE cleft upon the concentrations of both ATP and TMAO over the entire range of concentrations with just five undetermined parameters.     

Substrate-induced double sided H-bond network as a means of domain closure in 3-phosphoglycerate kinase

Closure of the two domains of 3-phosphoglycerate kinase, upon substrate binding, is essential for the enzyme function. The available crystal structures cannot provide sufficient information about the mechanism of substrate assisted domain closure and about the requirement of only one or both substrates, since lattice forces may hinder the large scale domain movements. In this study the known X-ray data, obtained for the open and closed conformations, were probed by solution small-angle X-ray scattering experiments. The results prove that binding of both substrates is essential for domain closure. Molecular graphical analysis, indeed, reveals formation of a double-sided H-bond network, which affects substantially the shape of the main molecular hinge at beta-strand L, under the concerted action of both substrates.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulations of NAD+-induced domain closure in horse liver alcohol dehydrogenase

Horse liver alcohol dehydrogenase is a homodimer, the protomer having a coenzyme-binding domain and a catalytic domain. Using all available x-ray structures and 50 ns of molecular dynamics simulations, we investigated the mechanism of NAD+-induced domain closure. When the well-known loop at the domain interface was modeled to its conformation in the closed structure, the NAD+-induced domain closure from the open structure could be simulated with remarkable accuracy. Native interactions in the closed structure between Arg369, Arg47, His51, Ala317, Phe319, and NAD+ were seen to form at different stages during domain closure. Removal of the Arg369 side-chain charge resulted in the loss of the tendency to close, verifying that specific interactions do help drive the domains closed. Further simulations and a careful analysis of x-ray structures suggest that the loop prevents domain closure in the absence of NAD+, and a cooperative mechanism operates between the subunits for domain closure. This cooperative mechanism explains the role of the loop as a block to closure because in the absence of NAD+ it would prevent the occurrence of an unliganded closed subunit when the other subunit closes on NAD+. Simulations that started with one subunit open and one closed supported this.     

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.     

A spring-loaded release mechanism regulates domain movement and catalysis in phosphoglycerate kinase

Phosphoglycerate kinase (PGK) is the enzyme responsible for the first ATP-generating step of glycolysis and has been implicated extensively in oncogenesis and its development. Solution small angle x-ray scattering (SAXS) data, in combination with crystal structures of the enzyme in complex with substrate and product analogues, reveal a new conformation for the resting state of the enzyme and demonstrate the role of substrate binding in the preparation of the enzyme for domain closure. Comparison of the x-ray scattering curves of the enzyme in different states with crystal structures has allowed the complete reaction cycle to be resolved both structurally and temporally. The enzyme appears to spend most of its time in a fully open conformation with short periods of closure and catalysis, thereby allowing the rapid diffusion of substrates and products in and out of the binding sites. Analysis of the open apoenzyme structure, defined through deformable elastic network refinement against the SAXS data, suggests that interactions in a mostly buried hydrophobic region may favor the open conformation. This patch is exposed on domain closure, making the open conformation more thermodynamically stable. Ionic interactions act to maintain the closed conformation to allow catalysis. The short time PGK spends in the closed conformation and its strong tendency to rest in an open conformation imply a spring-loaded release mechanism to regulate domain movement, catalysis, and efficient product release.