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.     

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.     

Communication between the nucleotide site and the main molecular hinge of 3-phosphoglycerate kinase

3-Phosphoglycerate kinase is a hinge-bending enzyme with substrate-assisted domain closure. However, the closure mechanism has not been described in terms of structural details. Here we present experimental evidence of the participation of individual substrate binding side chains in the operation of the main hinge which is distant from the substrate binding sites. The combined mutational, kinetic, and structural (DSC and SAXS) data for human 3-phosphoglycerate kinase have shown that catalytic residue R38, which also binds the substrate 3-phosphoglycerate, is essential in inducing domain closure. Similarly, residues K219, N336, and E343 which interact with the nucleotide substrates are involved in the process of domain closure. The other catalytic residue, K215, covers a large distance during catalysis but has no direct role in domain closure. The transmission path of the nucleotide effect toward the main hinge of PGK is described for the first time at the level of interactions existing in the tertiary structure.     

Biochemical and X-ray crystallographic studies on shikimate kinase: the important structural role of the P-loop lysine

Shikimate kinase, despite low sequence identity, has been shown to be structurally a member of the nucleoside monophosphate (NMP) kinase family, which includes adenylate kinase. In this paper we have explored the roles of residues in the P-loop of shikimate kinase, which forms the binding site for nucleotides and is one of the most conserved structural features in proteins. In common with many members of the P-loop family, shikimate kinase contains a cysteine residue 2 amino acids upstream of the essential lysine residue; the side chains of these residues are shown to form an ion pair. The C13S mutant of shikimate kinase was found to be enzymatically active, whereas the K15M mutant was inactive. However, the latter mutant had both increased thermostability and affinity for ATP when compared to the wild-type enzyme. The structure of the K15M mutant protein has been determined at 1.8 A, and shows that the organization of the P-loop and flanking regions is heavily disturbed. This indicates that, besides its role in catalysis, the P-loop lysine also has an important structural role. The structure of the K15M mutant also reveals that the formation of an additional arginine/aspartate ion pair is the most likely reason for its increased thermostability. From studies of ligand binding it appears that, like adenylate kinase, shikimate kinase binds substrates randomly and in a synergistic fashion, indicating that the two enzymes have similar catalytic mechanisms.     

The structure of the periplasmic ligand-binding domain of the sensor kinase CitA reveals the first extracellular PAS domain

The integral membrane sensor kinase CitA of Klebsiella pneumoniae is part of a two-component signal transduction system that regulates the transport and metabolism of citrate in response to its environmental concentration. Two-component systems are widely used by bacteria for such adaptive processes, but the stereochemistry of periplasmic ligand binding and the mechanism of signal transduction across the membrane remain poorly understood. The crystal structure of the CitAP periplasmic sensor domain in complex with citrate reveals a PAS fold, a versatile ligand-binding structural motif that has not previously been observed outside the cytoplasm or implicated in the transduction of conformational signals across the membrane. Citrate is bound in a pocket that is shared among many PAS domains but that shows structural variation according to the nature of the bound ligand. In CitAP, some of the citrate contact residues are located in the final strand of the central beta-sheet, which is connected to the C-terminal transmembrane helix. These secondary structure elements thus provide a potential conformational link between the periplasmic ligand binding site and the cytoplasmic signaling domains of the receptor.     

Synthesis and evaluation of heteroaryl substituted diazaspirocycles as scaffolds to probe the ATP-binding site of protein kinases

With the success of protein kinase inhibitors as drugs to target cancer, there is a continued need for new kinase inhibitor scaffolds. We have investigated the synthesis and kinase inhibition of new heteroaryl-substituted diazaspirocyclic compounds that mimic ATP. Versatile syntheses of substituted diazaspirocycles through ring-closing metathesis were demonstrated. Diazaspirocycles directly linked to heteroaromatic hinge binder groups provided ligand efficient inhibitors of multiple kinases, suitable as starting points for further optimization. The binding modes of representative diazaspirocyclic motifs were confirmed by protein crystallography. Selectivity profiles were influenced by the hinge binder group and the interactions of basic nitrogen atoms in the scaffold with acidic side-chains of residues in the ATP pocket. The introduction of more complex substitution to the diazaspirocycles increased potency and varied the selectivity profiles of these initial hits through engagement of the P-loop and changes to the spirocycle conformation, demonstrating the potential of these core scaffolds for future application to kinase inhibitor discovery.     

Solution structure and dynamics of an open beta-sheet, glycolytic enzyme, monomeric 23.7 kDa phosphoglycerate mutase from Schizosaccharomyces pombe

The structure and backbone dynamics of a double labelled (15N,13C) monomeric, 23.7 kD phosphoglycerate mutase (PGAM) from Schizosaccharomyces pombe have been investigated in solution using NMR spectroscopy. A set of 3125 NOE-derived distance restraints, 148 restraints representing inferred hydrogen bonds and 149 values of (3)J(HNHalpha) were used in the structure calculation. The mean rmsd from the average structure for all backbone atoms from residues 6-205 in the best 21 calculated structures was 0.59 A. The core of the enzyme includes an open, twisted, six-stranded beta-sheet flanked by four alpha-helices and a short 3(10)-helix. An additional smaller domain contains two short antiparallel beta-strands and a further pair of alpha-helices. The C(alpha) atoms of the S. pombe PGAM may be superimposed on their equivalents in one of the four identical subunits of Saccharomyces cerevisiae PGAM with an rmsd of 1.34 A (0.92 A if only the beta-sheet is considered). Small differences between the two structures are attributable partly to the deletion in the S. pombe sequence of a 25 residue loop involved in stabilising the S. cerevisiae tetramer. Analysis of 15N relaxation parameters indicates that PGAM tumbles isotropically with a rotational correlation time of 8.7 ns and displays a range of dynamic features. Of 178 residues analysed, only 77 could be fitted without invoking terms for fast internal motion or chemical exchange, and out of the remainder, 77 required a chemical exchange term. Significantly, 46 of the slowly exchanging (milli- to microsecond) residues lie in helices, and these account for two-thirds of all analysed helix residues. On the contrary, only one beta-sheet residue required an exchange term. In contrast to other analyses of backbone dynamics reported previously, residues in slow exchange appeared to correlate with architectural features of the enzyme rather than congregating close to ligand binding sites.     

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

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

Phosphorylation- and Nucleotide-Binding-Induced Changes to the Stability and Hydrogen Exchange Patterns of JNK1β1 Provide Insight into Its Mechanisms of Activation

Many studies have characterized how changes to the stability and internal motions of a protein during activation can contribute to their catalytic function, even when structural changes cannot be observed. Here, unfolding studies and hydrogen–deuterium exchange (HX) mass spectrometry were used to investigate the changes to the stability and conformation/conformational dynamics of JNK1β1 induced by phosphorylative activation. Equivalent studies were also employed to determine the effects of nucleotide binding on both inactive and active JNK1β1 using the ATP analogue, 5?-adenylyl-imidodiphosphate (AMP-PNP). JNK1β1 phosphorylation alters HX in regions involved in catalysis and substrate binding, changes that can be ascribed to functional modifications in either structure and/or backbone flexibility. Increased HX in the hinge between the N- and C-terminal domains implied that it acquires enhanced flexibility upon phosphorylation that may be a prerequisite for interdomain closure. In combination with the finding that nucleotide binding destabilizes the kinase, the patterns of solvent protection by AMP-PNP were consistent with a novel mode of nucleotide binding to the C-terminal domain of a destabilized and open domain conformation of inactive JNK1β1. Solvent protection by AMP-PNP of both N- and C-terminal domains in active JNK1β1 revealed that the domains close around nucleotide upon phosphorylation, concomitantly stabilizing the kinase. This suggests that phosphorylation activates JNK1β1 in part by increasing hinge flexibility to facilitate interdomain closure and the creation of a functional active site. By uncovering the complex interplay that occurs between nucleotide binding and phosphorylation, we present new insight into the unique mechanisms by which JNK1β1 is regulated.     

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.     

Solution NMR study of the monomeric form of p13suc1 protein sheds light on the hinge region determining the affinity for a phosphorylated substrate.

Cyclin-dependent kinase subunit (CKS) proteins bind to cyclin-dependent kinases and target various proteins to phosphorylation and proteolysis during cell division. Crystal structures showed that CKS can exist both in a closed monomeric conformation when bound to the kinase and in an inactive C-terminal beta-strand-exchanged conformation. With the exception of the hinge loop, however, both crystal structures are identical, and no new protein interface is formed in the dimer. Protein engineering studies have pinpointed the crucial role of the proline 90 residue of the p13(suc1) CKS protein from Schizosaccharomyces pombe in the monomer-dimer equilibrium and have led to the concept of a loaded molecular spring of the beta-hinge motif. Mutation of this hinge proline into an alanine stabilizes the protein and prevents the occurrence of swapping. However, other mutations further away from the hinge as well as ligand binding can equally shift the equilibrium between monomer and dimer. To address the question of differential affinity through relief of the strain, here we compare the ligand binding of the monomeric form of wild-type S. pombe p13(suc1) and its hinge mutant P90A in solution by NMR spectroscopy. We indeed observed a 5-fold difference in affinity with the wild-type protein being the most strongly binding. Our structural study further indicates that both wild-type and the P90A mutant proteins adopt in solution the closed conformation but display different dynamic properties in the C-terminal beta-sheet involved in domain swapping and protein interactions.     

Crystal structure of the complex of plasminogen activator inhibitor 2 with a peptide mimicking the reactive center loop

The structure of the serpin, plasminogen activator inhibitor type-2 (PAI-2), in a complex with a peptide mimicking its reactive center loop (RCL) has been determined at 1.6-A resolution. The structure shows the relaxed state serpin structure with a prominent six-stranded beta-sheet. Clear electron density is seen for all residues in the peptide. The P1 residue of the peptide binds to a well defined pocket at the base of PAI-2 that may be important in determining the specificity of protease inhibition. The stressed-to-relaxed state (S --> R) transition in PAI-2 can be modeled as the relative motion between a quasirigid core domain and a smaller segment comprising helix hF and beta-strands s1A, s2A, and s3A. A comparison of the Ramachandran plots of the stressed and relaxed state PAI-2 structures reveals the location of several hinge regions connecting these two domains. The hinge regions cluster in three locations on the structure, ensuring a cooperative S --> R transition. We hypothesize that the hinge formed by the conserved Gly(206) on beta-strand s3A in the breach region of PAI-2 effects the S --> R transition by altering its backbone torsion angles. This torsional change is due to the binding of the P14 threonine of the RCL to the open breach region of PAI-2.     

Thermodynamic analysis of substrate induced domain closure of 3-phosphoglycerate kinase

The energetic changes accompanying domain closure of 3-phosphoglycerate kinase, a typical hinge-bending enzyme, were assessed. Calorimetric titrations of the enzyme with each substrate, both in the absence and presence of the other one, provide information not only about the energetics of substrate binding, but of the associated conformational changes, including domain closure. Our results suggest that conformational rearrangements in the hinge generated by binding of both substrates provide the main driving force for domain closure overcoming the slightly unfavourable contact interactions between the domains.     

Substrate binding modifies the hinge bending characteristics of human 3‐phosphoglycerate kinase: A molecular dynamics study

3‐Phosphogycerate kinase (PGK) is a two domain enzyme, with a binding site of the 1,3‐bisphosphoglycerate on the N‐domain and of the ADP on the C‐domain. To transfer a phosphate group the enzyme has to undergo a hinge bending motion from open to closed conformation to bring the substrates to close proximity. Molecular dynamics simulation was used to elucidate the effect of ligand binding onto the domain motions of this enzyme. The simulation results of the apo form indicate a hinge bending motion in the ns timescale while the time period of the hinge bending motion of the complex form is clearly over the 20 ns simulation time. The apo form exhibits several hinge points that contribute to the hinge bending motion while upon binding the ligands, the hinge bending becomes strictly restrained with one dominant hinge point in the vicinity of the substrates. At the same time, ligand binding results in an enhanced correlation of internal domain motions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

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.     

Mutation of a highly conserved aspartate residue in subdomain IX abolishes Fer protein-tyrosine kinase activity

Before the structure of cAMP-dependent protein kinase had been solved, sequence alignments had already suggested that several highly conserved peptide motifs described as kinase subdomains I through XI might play some functional role in catalysis. Crystal structures of several members of the protein kinase superfamily have suggested that the nearly invariant aspartate residue within subdomain IX contributes to the conformational stability of the catalytic loop by forming hydrogen bonds with backbone amides within subdomain VI. However, substitution of this aspartate with alanine or threonine in some protein kinases have indicated that these interactions are not essential for activity. In contrast, we show here that conversion of this aspartate to arginine abolished the catalytic activity of the Fer protein-tyrosine kinase when expressed either in mammalian cells or in bacteria. Structural modeling predicted that the catalytic loop of the FerD743R mutant was disrupted by van der Waal's repulsion between the side chains of the substituted arginine residue in subdomain IX and histidine-683 in subdomain VI. The FerD743R mutant model predicted a shift in the peptide backbone of the catalytic loop, and an outward rotation of histidine-683 and arginine-684 side chains. However, the position and orientation of the presumptive catalytic base, aspartate-685, was not substantially changed. The proposed model explains how substitutions of some, but not all residues could be tolerated at this nearly invariant aspartate in kinase subdomain IX.     

Structural characterization of a proline-driven conformational switch within the Itk SH2 domain

Interleukin-2 tyrosine kinase (Itk) is a T cell-specific kinase required for a proper immune response following T cell receptor engagement. In addition to the kinase domain, Itk is composed of several noncatalytic regulatory domains, including a Src homology 2 (SH2) domain that contains a conformationally heterogeneous Pro residue. Cis-trans isomerization of a single prolyl imide bond within the SH2 domain mediates conformer-specific ligand recognition that may have functional implications in T cell signaling. To better understand the mechanism by which a proline switch regulates ligand binding, we have used NMR spectroscopy to determine two structures of Itk SH2 corresponding to the cis and trans imide bond-containing conformers. The structures indicate that the heterogeneous Pro residue acts as a hinge that modulates ligand recognition by controlling the relative orientation of protein-binding surfaces.     

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.     

Crystal structure of diaminopelargonic acid synthase: evolutionary relationships between pyridoxal-5'-phosphate-dependent enzymes.

The three-dimensional structure of diaminopelargonic acid synthase, a vitamin B6-dependent enzyme in the pathway of the biosynthesis of biotin, has been determined to 1.8 A resolution by X-ray crystallography. The structure was solved by multi-wavelength anomalous diffraction techniques using a crystal derivatized with mercury ions. The protein model has been refined to a crystallographic R -value of 17.5% (R -free 22.6%). Each enzyme subunit consists of two domains, a large domain (residues 50-329) containing a seven-stranded predominantly parallel beta-sheet, surrounded by alpha-helices, and a small domain comprising residues 1-49 and 330-429. Two subunits, related by a non-crystallographic dyad in the crystals, form the homodimeric molecule, which contains two equal active sites. Pyridoxal-5'-phosphate is bound in a cleft formed by both domains of one subunit and the large domain of the second subunit. The cofactor is anchored to the enzyme by a covalent linkage to the side-chain of the invariant residue Lys274. The phosphate group interacts with main-chain nitrogen atoms and the side-chain of Ser113, located at the N terminus of an alpha-helix. The pyridine nitrogen forms a hydrogen bond to the side-chain of the invariant residue Asp245. Electron density corresponding to a metal ion, most likely Na(+), was found in a tight turn at the surface of the enzyme. Structure analysis reveals that diaminopelargonic acid synthase belongs to the family of vitamin B6-dependent aminotransferases with the same fold as originally observed in aspartate aminotransferase. A multiple structure alignment of enzymes in this family indicated that they form at least six different subclasses. Striking differences in the fold of the N-terminal part of the polypeptide chain are one of the hallmarks of these subclasses. Diaminopelargonic acid synthase is a member of the aminotransferase subclass III. From the structure of the non-productive complex of the holoenzyme with the substrate 7-keto-8-aminopelargonic acid the location of the active site and residues involved in substrate binding have been identified.