Intrinsic domain and loop dynamics commensurate with catalytic turnover in an induced-fit enzyme

Structure of McsB, a protein kinase for regulated arginine phosphorylation

Protein phosphorylation regulates key processes in all organisms. In Gram-positive bacteria, protein arginine phosphorylation plays a central role in protein quality control by regulating transcription factors and marking aberrant proteins for degradation. Here, we report structural, biochemical, and in vivo data of the responsible kinase, McsB, the founding member of an arginine-specific class of protein kinases. McsB differs in structure and mechanism from protein kinases that act on serine, threonine, and tyrosine residues and instead has a catalytic domain related to that of phosphagen kinases (PhKs), metabolic enzymes that phosphorylate small guanidino compounds. In McsB, the PhK-like phosphotransferase domain is structurally adapted to target protein substrates and is accompanied by a novel phosphoarginine (pArg)-binding domain that allosterically controls protein kinase activity. The identification of distinct pArg reader domains in this study points to a remarkably complex signaling system, thus challenging simplistic views of bacterial protein phosphorylation.

Elevated μs-ms timescale backbone dynamics in the transition state analog form of arginine kinase

Arginine kinase catalyzes reversible phosphoryl transfer between arginine and ATP. Crystal structures of arginine kinase in an open, substrate-free form and closed, transition state analog (TSA) complex indicate that the enzyme undergoes substantial domain and loop rearrangements required for substrate binding, catalysis, and product release. Nuclear magnetic resonance (NMR) has shown that substrate-free arginine kinase is rigid on the ps-ns timescale (average S²=0.84±0.08) yet quite dynamic on the µs-ms timescale (35 residues with R_ex, 12%), and that movements of the N-terminal domain and the loop comprising residues I182-G209 are rate-limiting on catalysis. Here, NMR of the TSA-bound enzyme shows similar rigidity on the ps-ns timescale (average S²=0.91±0.05) and substantially increased μs-ms timescale dynamics (77 residues; 22%). Many of the residues displaying μs-ms dynamics in NMR Carr-Purcell-Meiboom-Gill (CPMG) ¹⁵N backbone relaxation dispersion experiments of the TSA complex are also dynamic in substrate-free enzyme. However, the presence of additional dynamic residues in the TSA-bound form suggests that dynamics extend through much of the C-terminal domain, which indicates that in the closed form, a larger fraction of the protein takes part in conformational transitions to the excited state(s). Conformational exchange rate constants (k_ex) of the TSA complex are all approximately 2500s^-1, higher than any observed in the substrate-free enzyme (800-1900s^-1). Elevated μs-ms timescale protein dynamics in the TSA-bound enzyme is more consistent with recently postulated catalytic networks involving multiple interconnected states at each step of the reaction, rather than a classical single stabilized transition state.

The Michaelis Complex of Arginine Kinase Samples the Transition State at a Frequency That Matches the Catalytic Rate

Arginine kinase (AK), which is a member of the phosphagen kinase family, serves as a model system for studying the structural and dynamic determinants of biomolecular enzyme catalysis of all major states involved of the enzymatic cycle. These states are the apo state (substrate free), the Michaelis complex analogue AK:Arg:Mg·AMPPNP (MCA), a product complex analogue AK:pAIE:Mg·ADP (PCA), and the transition state analogue AK:Arg:Mg·ADP:NO₃^- (TSA). The conformational dynamics of these states have been studied by NMR relaxation dispersion measurements of the methyl groups of the Ile, Leu, and Val residues at two static magnetic fields. Although all states undergo significant amounts of μs-ms time scale dynamics, only the MCA samples a dominant excited state that resembles the TSA, as evidenced by the strong correlation between the relaxation dispersion derived chemical shift differences Δω and the equilibrium chemical shift differences Δδ of these states. The average lifetime of the MCA is 36 ms and the free energy difference to the TSA-like form is 8.5 kJ/mol. It is shown that the conformational energy landscape of the Michaelis complex analogue is shaped in a way that at room temperature it channels passage to the transition state, thereby determining the rate-limiting step of the phosphorylation reaction of arginine. Conversely, relaxation dispersion experiments of the TSA reveal that it samples the structures of the Michaelis complex analogue or the apo state as its dominant excited state. This reciprocal behavior shows that the free energy of the TSA, with all ligands bound, is lower by only about 8.9 kJ/mol than that of the Michaelis or apo complex conformations with the TSA ligands present.

Protein dynamics and function from solution state NMR spectroscopy

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

The Sampling of Conformational Dynamics in Ambient-Temperature Crystal Structures of Arginine Kinase

Arginine kinase provides a model for functional dynamics, studied through crystallography, enzymology, and nuclear magnetic resonance. Structures are now solved, at ambient temperature, for the transition state analog (TSA) complex. Analysis of quasi-rigid sub-domain displacements show that differences between the two TSA structures average about 5% of changes between substrate-free and TSA forms, and they are nearly co-linear. Small backbone hinge rotations map to sites that also flex on substrate binding. Anisotropic atomic displacement parameters (ADPs) are refined using rigid-body TLS constraints. Consistency between crystal forms shows that they reflect intrinsic molecular properties more than crystal lattice effects. In many regions, the favored directions of thermal/static displacement are appreciably correlated with movements on substrate binding. Correlation between ADPs and larger substrate-associated movements implies that the latter approximately follow paths of low-energy intrinsic motions.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Rapid Determination of Fast Protein Dynamics from NMR Chemical Exchange Saturation Transfer Data

Functional motions of (15)N-labeled proteins can be monitored by solution NMR spin relaxation experiments over a broad range of timescales. These experiments however typically take of the order of several days to a week per protein. Recently, NMR chemical exchange saturation transfer (CEST) experiments have emerged to probe slow millisecond motions complementing R1ρ and CPMG-type experiments. CEST also simultaneously reports on site-specific R1 and R2 parameters. It is shown here how CEST-derived R1 and R2 relaxation parameters can be measured within a few hours at an accuracy comparable to traditional relaxation experiments. Using a "lean" version of the model-free approach S(2) order parameters can be determined that match those from the standard model-free approach applied to (15)N R1, R2 , and {(1)H}-(15)N NOE data. The new methodology, which is demonstrated for ubiquitin and arginine kinase (42 kDa), should serve as an effective screening tool of protein dynamics from picosecond-to-millisecond timescales.

HdeB chaperone activity is coupled to its intrinsic dynamic properties

Enteric bacteria encounter extreme acidity when passing through hosts' stomach. Since the bacterial periplasmic space quickly equilibrates with outer environment, an efficient acid resistance mechanism is essential in preventing irreversible protein denaturation/aggregation and maintaining bacteria viability. HdeB, along with its homolog HdeA, was identified as a periplasmic acid-resistant chaperone. Both proteins exist as homodimers and share similar monomeric structures under neutral pH, while showing different dimeric packing interfaces. Previous investigations show that HdeA functions through an acid-induced dimer-to-monomer transition and partial unfolding at low pH (pH 2-3), resulting in exposure of hydrophobic surfaces that bind substrate proteins. In contrast, HdeB appears to have a much higher optimal activation pH (pH 4-5), under which condition the protein maintains a well-folded dimer and the mechanism for its chaperone activity remains elusive. Herein, we present an NMR study of HdeB to investigate its dynamic properties. Our results reveal that HdeB undergoes significant micro- to milli-second timescale conformational exchanges at neutral to near-neutral pH, under the later condition it exhibits optimal activity. The current study indicates that HdeB activation is coupled to its intrinsic dynamics instead of structural changes, and therefore its functional mechanism is apparently different from HdeA.

Parsimony in Protein Conformational Change

Protein conformational change is analyzed by finding the minimalist backbone torsion angle rotations that superpose crystal structures within experimental error. Of several approaches for enforcing parsimony during flexible least-squares superposition, an ℓ(1)-norm restraint provided greatest consistency with independent indications of flexibility from nuclear magnetic resonance relaxation dispersion and chemical shift perturbation in arginine kinase and four previously studied systems. Crystallographic cross-validation shows that the dihedral parameterization describes conformational change more accurately than rigid-group approaches. The rotations that superpose the principal elements of structure constitute a small fraction of the raw (φ, ψ) differences that also reflect local conformation and experimental error. Substantial long-range displacements can be mediated by modest dihedral rotations, accommodated even within α helices and β sheets without disruption of hydrogen bonding at the hinges. Consistency between ligand-associated and intrinsic motions (in the unliganded state) implies that induced changes tend to follow low-barrier paths between conformational sub-states that are in intrinsic dynamic equilibrium.

Identifying sequential substrate binding at the single-molecule level by enzyme mechanical stabilization

Enzyme-substrate binding is a dynamic process intimately coupled to protein structural changes, which in turn changes the unfolding energy landscape. By the use of single-molecule force spectroscopy (SMFS), we characterize the open-to-closed conformational transition experienced by the hyperthermophilic adenine diphosphate (ADP)-dependent glucokinase from Thermococcus litoralis triggered by the sequential binding of substrates. In the absence of substrates, the mechanical unfolding of TlGK shows an intermediate 1, which is stabilized in the presence of Mg·ADP(-), the first substrate to bind to the enzyme. However, in the presence of this substrate, an additional unfolding event is observed, intermediate 1*. Finally, in the presence of both substrates, the unfolding force of intermediates 1 and 1* increases as a consequence of the domain closure. These results show that SMFS can be used as a powerful experimental tool to investigate binding mechanisms of different enzymes with more than one ligand, expanding the repertoire of protocols traditionally used in enzymology.

The substrate-free and -bound crystal structures of the duplicated taurocyamine kinase from the human parasite Schistosoma mansoni

The taurocyamine kinase from the blood fluke Schistosoma mansoni (SmTK) belongs to the phosphagen kinase (PK) family and catalyzes the reversible Mg(2+)-dependent transfer of a phosphoryl group between ATP and taurocyamine. SmTK is derived from gene duplication, as are all known trematode TKs. Our crystallographic study of SmTK reveals the first atomic structure of both a TK and a PK with a bilobal structure. The two unliganded lobes present a canonical open conformation and interact via their respective C- and N-terminal domains at a helix-mediated interface. This spatial arrangement differs from that observed in true dimeric PKs, in which both N-terminal domains make contact. Our structures of SmTK complexed with taurocyamine or l-arginine compounds explain the mechanism by which an arginine residue of the phosphagen specificity loop is crucial for substrate specificity. An SmTK crystal was soaked with the dead end transition state analog (TSA) components taurocyamine-NO3 (2-)-MgADP. One SmTK monomer was observed with two bound TSAs and an asymmetric conformation, with the first lobe semiclosed and the second closed. However, isothermal titration calorimetry and enzyme kinetics experiments showed that the two lobes function independently. A small angle x-ray scattering model of SmTK-TSA in solution with two closed active sites was generated.

Enzyme dynamics from NMR spectroscopy

Biological activities of enzymes, including regulation or coordination of mechanistic stages preceding or following the chemical step, may depend upon kinetic or equilibrium changes in protein conformations. Exchange of more open or flexible conformational states with more closed or constrained states can influence inhibition, allosteric regulation, substrate recognition, formation of the Michaelis complex, side reactions, and product release. NMR spectroscopy has long been applied to the study of conformational dynamic processes in enzymes because these phenomena can be characterized over multiple time scales with atomic site resolution. Laboratory-frame spin-relaxation measurements, sensitive to reorientational motions on picosecond–nanosecond time scales, and rotating-frame relaxation-dispersion measurements, sensitive to chemical exchange processes on microsecond–millisecond time scales, provide information on both conformational distributions and kinetics. This Account reviews NMR spin relaxation studies of the enzymes ribonuclease HI from mesophilic (Escherichia coli) and thermophilic (Thermus thermophilus) bacteria, E. coli AlkB, and Saccharomyces cerevisiae triosephosphate isomerase to illustrate the contributions of conformational flexibility and dynamics to diverse steps in enzyme mechanism.

Spin relaxation measurements and molecular dynamics (MD) simulations of the bacterial ribonuclease H enzymes show that the handle region, one of three loop regions that interact with substrates, interconverts between two conformations. Comparison of these conformations with the structure of the complex between Homo sapiens ribonuclease H and a DNA:RNA substrate suggests that the more closed state is inhibitory to binding. The large population of the closed conformation in T. thermophilus ribonuclease H contributes to the increased Michaelis constant compared with the E. coli enzyme. NMR spin relaxation and fluorescence spectroscopy have characterized a conformational transition in AlkB between an open state, in which the side chains of methionine residues in the active site are disordered, and a closed state, in which these residues are ordered. The open state is highly populated in the AlkB/Zn(II) complex, and the closed state is highly populated in the AlkB/Zn(II)/2OG/substrate complex, in which 2OG is the 2-oxoglutarate cosubstrate and the substrate is a methylated DNA oligonucleotide. The equilibrium is shifted to approximately equal populations of the two conformations in the AlkB/Zn(II)/2OG complex. The conformational shift induced by 2OG ensures that 2OG binds to AlkB/Zn(II) prior to the substrate. In addition, the opening rate of the closed conformation limits premature release of substrate, preventing generation of toxic side products by reaction with water. Closure of active site loop 6 in triosephosphate isomerase is critical for forming the Michaelis complex, but reopening of the loop after the reaction is (partially) rate limiting. NMR spin relaxation and MD simulations of triosephosphate isomerase in complex with glycerol 3-phosphate demonstrate that closure of loop 6 is a highly correlated rigid-body motion. The MD simulations also indicate that motions of Gly173 in the most flexible region of loop 6 contribute to opening of the active site loop for product release. Considered together, these three enzyme systems illustrate the power of NMR spin relaxation investigations in providing global insights into the role of conformational dynamic processes in the mechanisms of enzymes from initial activation to final product release.

Backbone resonance assignments of the 42 kDa enzyme arginine kinase in the transition state analogue form

Nearly complete backbone resonance assignments for the 357 residue, 42 kDa enzyme arginine kinase in a transition state analogue (TSA) complex are presented. The TSA is a quaternary complex of arginine kinase, MgADP, arginine, and nitrate. About 93% (320 of 344) of the non-proline backbone amides were assigned using an enzyme enriched with (2)H, (13)C, and (15)N in combination with three enzyme samples prepared with a single (15)N-labeled amino acid (K, L, and R). The amide assignments will provide the foundation for investigating the dynamics of arginine kinase when in a TSA complex.

Crystal structures of arginine kinase in complex with ADP, nitrate, and various phosphagen analogs

Arginine kinase catalyzes the reversible transfer of a phosphoryl group between ATP and l-arginine and is a monomeric homolog of the human enzyme creatine kinase. Arginine and creatine kinases belongs to the phosphagen kinase family of enzymes, which consists of eight known members, each of which is specific for its own phosphagen. Here, the source of phosphagen specificity in arginine kinase is investigated through the use of phosphagen analogs. Crystal structures have been determined for Limulus polyphemus arginine kinase with one of four arginine analogs bound in a transition state analog complex: l-ornithine, l-citrulline, imino-l-ornithine, and d-arginine. In all complexes, the enzyme achieves a closed conformation very similar to that of the cognate transition state analog complex, but differences are observed in the configurations of bound ligands. Arginine kinase exhibits no detectable activity towards ornithine, citrulline, or imino-l-ornithine, and only trace activity towards d-arginine. The crystal structures presented here demonstrate that phosphagen specificity is derived neither from a lock-and-key mechanism nor a modulation of induced-fit conformational changes, but potentially from subtle distortions in bound substrate configurations.

Dynamical properties of the loop 320s of substrate-free and substrate-bound muscle creatine kinase by NMR: evidence for independent subunits

Muscle creatine kinase (MCK; EC2.7.3.2) is a 86 kDa homodimer that belongs to the family of guanidino kinases. MCK has been intensively studied for several decades, but it is still not known why it is a dimer because this quaternary structure does not translate into obvious structural or functional advantages over the homologous monomeric arginine kinase. In particular, it remains to be demonstrated whether MCK subunits are independent. Here, we describe NMR chemical-shift perturbation and relaxation experiments designed to study the active site 320s flexible loop of this enzyme. The analysis was performed with the enzyme in its ligand-free and MgADP-complexed forms, as well as with the transition-state analogue abortive complex (MCK-Mg-ADP-creatine-nitrate ion). Our data indicate that each subunit can bind substrates independently.

Understanding the molecular mechanism of enzyme dynamics of ribonuclease A through protonation/deprotonation of HIS48

Molecular dynamics simulation is carried out to investigate the enzyme dynamics of RNase A with the HIS48 in three different states (HIP48 (protonated), HID48 (deprotonated), and H48A mutant). Insights derived from the current theoretical study, combined with the available experimental observations, enabled us to provide a microscopic picture for the efficient enzyme dynamics. Specifically, in the "closed" state or HIP48, the N-terminal hinge loop is intact and the enzyme remains in a relatively stable conformation which is preferred for catalytic reaction. Deprotonation of HIS48 induces the denaturing of this hinge-loop into a 3(10)-helix, causing it to break the original interaction network around the loop-1 and drive the partial unfolding of the N-terminal. The enhanced dynamic motion of the N-terminal helix facilitates the release of the catalytic product (the rate limiting step) and speeds up the overall catalytic process. The current study established that HIS49 acts as a modulator for the transformation of conformational states through the perturbing of hydrogen bond networks across loop-1, the N-terminal helix, and other residues nearby. Our study suggests that HIS48 may also serve to transport loop-1's kinetic energy to the reaction center.

Rate-limiting domain and loop motions in arginine kinase

Arginine kinase catalyzes the reversible transfer of a phosphoryl group between ATP and arginine. It is the arthropod homologue of creatine kinase, buffering cellular ATP levels. Crystal structures of arginine kinase, in substrate-free and substrate-bound forms, have revealed large conformational changes associated with the catalytic cycle. Recent nuclear magnetic resonance identified movements of the N-terminal domain and a loop comprising residues I182--G209 with conformational exchange rates in the substrate-free enzyme similar to the turnover rate. Here, to understand whether these motions might be rate-limiting, we determined activation barriers for both the intrinsic dynamics and enzyme turnover using measurements over a temperature range of 15-30 °C. (15)N transverse relaxation dispersion yields activation barriers of 46 ± 8 and 34 ± 12 kJ/mol for the N-terminal domain and I182--G209 loop, respectively. An activation barrier of 34 ± 13 kJ/mol was obtained for enzyme turnover from steady-state kinetics. The similarity between the activation barriers is indeed consistent with turnover being limited by backbone conformational dynamics and pinpoints the locations of potentially rate-limiting motions.

The structure of lombricine kinase: implications for phosphagen kinase conformational changes

Lombricine kinase is a member of the phosphagen kinase family and a homolog of creatine and arginine kinases, enzymes responsible for buffering cellular ATP levels. Structures of lombricine kinase from the marine worm Urechis caupo were determined by x-ray crystallography. One form was crystallized as a nucleotide complex, and the other was substrate-free. The two structures are similar to each other and more similar to the substrate-free forms of homologs than to the substrate-bound forms of the other phosphagen kinases. Active site specificity loop 309-317, which is disordered in substrate-free structures of homologs and is known from the NMR of arginine kinase to be inherently dynamic, is resolved in both lombricine kinase structures, providing an improved basis for understanding the loop dynamics. Phosphagen kinases undergo a segmented closing on substrate binding, but the lombricine kinase ADP complex is in the open form more typical of substrate-free homologs. Through a comparison with prior complexes of intermediate structure, a correlation was revealed between the overall enzyme conformation and the substrate interactions of His(178). Comparative modeling provides a rationale for the more relaxed specificity of these kinases, of which the natural substrates are among the largest of the phosphagen substrates.

Arginine kinase: joint crystallographic and NMR RDC analyses link substrate-associated motions to intrinsic flexibility

The phosphagen kinase family, including creatine and arginine kinases (AKs), catalyzes the reversible transfer of a "high-energy" phosphate between ATP and a phosphoguanidino substrate. They have become a model for the study of both substrate-induced conformational change and intrinsic protein dynamics. Prior crystallographic studies indicated large substrate-induced domain rotations, but differences among a recent set of AK structures were interpreted as a plastic deformation. Here, the structure of Limulus substrate-free AK is refined against high-resolution crystallographic data and compared quantitatively with NMR chemical shifts and residual dipolar couplings (RDCs). This demonstrates the feasibility of this type of RDC analysis of proteins that are large by NMR standards (42 kDa) and illuminates the solution structure, free from crystal-packing constraints. Detailed comparison of the 1.7 Å resolution substrate-free crystal structure against the 1.7 Å transition-state analog complex shows large substrate-induced domain motions that can be broken down into movements of smaller quasi-rigid bodies. The solution-state structure of substrate-free AK is most consistent with an equilibrium of substrate-free and substrate-bound structures, with the substrate-free form dominating, but with varying displacements of the quasi-rigid groups. Rigid-group rotations evident from the crystal structures are about axes previously associated with intrinsic millisecond dynamics using NMR relaxation dispersion. Thus, "substrate-induced" motions are along modes that are intrinsically flexible in the substrate-free enzyme and likely involve some degree of conformational selection.

Enzyme dynamics point to stepwise conformational selection in catalysis

Recent data increasingly reveal that conformational dynamics are indispensable to enzyme function throughout the catalytic cycle, in substrate recruiting, chemical transformation, and product release. Conformational transitions may involve conformational selection and induced fit, which can be viewed as a special case in the catalytic network. NMR, X-ray crystallography, single-molecule FRET, and simulations clearly demonstrate that the free enzyme dynamics already encompass all the conformations necessary for substrate binding, preorganization, transition-state stabilization, and product release. Conformational selection and substate population shift at each step of the catalytic turnover can accommodate enzyme specificity and efficiency. Within such a framework, entropy can have a larger role in conformational dynamics than in direct energy transfer in dynamically promoted catalysis.

Direct evidence of conformational heterogeneity in human pancreatic glucokinase from high-resolution nuclear magnetic resonance

High-resolution nuclear magnetic resonance is used to investigate the conformational dynamics of human glucokinase, a 52 kDa monomeric enzyme that displays kinetic cooperativity. (1)H-(15)N transverse relaxation optimized spectra of uniformly labeled glucokinase, recorded in the absence and presence of glucose, reveal significant cross-peak overlap and heterogeneous peak intensities that persist over a range of temperatures. (15)N-specific labeling of isoleucines and tryptophans, reporting on backbone and side chain dynamics, respectively, demonstrates that both unliganded and glucose-bound enzymes sample multiple conformations, although glucose stabilizes certain conformations. These results provide the first direct evidence of glucokinase conformational heterogeneity and hence shed light on the molecular basis of cooperativity.