Mg2+-free Bacillus stearothermophilus tryptophanyl-tRNA synthetase retains a major fraction of the overall rate enhancement for tryptophan activation

Domain acquisition by class I aminoacyl-tRNA synthetase urzymes coordinated the catalytic functions of HVGH and KMSKS motifs

Leucyl-tRNA synthetase (LeuRS) is a Class I aminoacyl-tRNA synthetase (aaRS) that synthesizes leucyl-tRNAleu for codon-directed protein synthesis. Two signature sequences, HxGH and KMSKS help stabilize transition-states for amino acid activation and tRNA aminoacylation by all Class I aaRS. Separate alanine mutants of each signature, together with the double mutant, behave in opposite ways in Pyrococcus horikoshii LeuRS and the 129-residue urzyme ancestral model generated from it (LeuAC). Free energy coupling terms, Δ(ΔG‡), for both reactions are large and favourable for LeuRS, but unfavourable for LeuAC. Single turnover assays with 32Pα-ATP show correspondingly different internal products. These results implicate domain motion in catalysis by full-length LeuRS. The distributed thermodynamic cycle of mutational changes authenticates LeuAC urzyme catalysis far more convincingly than do single point mutations. Most importantly, the evolutionary gain of function induced by acquiring the anticodon-binding (ABD) and multiple insertion modules in the catalytic domain appears to be to coordinate the catalytic function of the HxGH and KMSKS signature sequences. The implication that backbone elements of secondary structures achieve a major portion of the overall transition-state stabilization by LeuAC is also consistent with coevolution of the genetic code and metabolic pathways necessary to produce histidine and lysine sidechains.

Microcalorimetry reveals multi-state thermal denaturation of G. stearothermophilus tryptophanyl-tRNA synthetase

Mechanistic studies of Geobacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) afford an unusually detailed description-the escapement mechanism-for the distinct steps coupling catalysis to domain motion, efficiently converting the free energy of ATP hydrolysis into biologically useful alternative forms of information and work. Further elucidation of the escapement mechanism requires understanding thermodynamic linkages between domain configuration and conformational stability. To that end, we compare experimental thermal melting of fully liganded and apo TrpRS with a computational simulation of the melting of its fully liganded form. The simulation also provides important structural cameos at successively higher temperatures, enabling more confident interpretation. Experimental and simulated melting both proceed through a succession of three transitions at successively higher temperature. The low-temperature transition occurs at approximately the growth temperature of the organism and so may be functionally relevant but remains too subtle to characterize structurally. Structural metrics from the simulation imply that the two higher-temperature transitions entail forming a molten globular state followed by unfolding of secondary structures. Ligands that stabilize the enzyme in a pre-transition (PreTS) state compress the temperature range over which these transitions occur and sharpen the transitions to the molten globule and fully denatured states, while broadening the low-temperature transition. The experimental enthalpy changes provide a key parameter necessary to convert changes in melting temperature of combinatorial mutants into mutationally induced conformational free energy changes. The TrpRS urzyme, an excerpted model representing an early ancestral form, containing virtually the entire catalytic apparatus, remains largely intact at the highest simulated temperatures.

High-throughput thermal denaturation of tryptophanyl-tRNA synthetase combinatorial mutants reveals high-order energetic coupling determinants of conformational stability

Landscape descriptions provide a framework for identifying functionally significant dynamic linkages in proteins but cannot supply details. Rate measurements of combinatorial mutations can implicate dynamic linkages in catalysis. A major difficulty is filtering dynamic linkages from the vastly more numerous static interactions that stabilize domain folding. The Geobacillus stearothermophilus (TrpRS) D1 switch is such a dynamic packing motif; it links domain movement to catalysis and specificity. We describe Thermofluor and far UV circular dichroism melting curves for all 16 D1 switch variants to determine their higher-order impact on unliganded TrpRS stability. A prominent transition at intermediate temperatures in TrpRS thermal denaturation is molten globule formation. Combinatorial analysis of thermal melting transcends the protein landscape in four significant respects: (i) bioinformatic methods identify dynamic linkages from coordinates of multiple conformational states. (ii) Relative mutant melting temperatures, δTM, are proportional to free energy changes. (iii) Structural analysis of thermal melting implicates unexpected coupling between the D1 switch packing and regions of high local frustration. Those segments develop molten globular characteristics at the point of greatest complementarity to the chemical transition state and are the first TrpRS structures to melt. (iv) Residue F37 stabilizes both native and molten globular states; its higher-order interactions modify the relative intrinsic impacts of mutations to other D1 switch residues from those estimated for single point mutants. The D1 switch is a central component of an escapement mechanism essential to free energy transduction. These conclusions begin to relate the escapement mechanism to differential TrpRS conformational stabilities.

Putting amino acids onto tRNAs: The aminoacyl-tRNA synthetases as catalysts

In this chapter we consider the catalytic approaches used by aminoacyl-tRNA synthetase (AARS) enzymes to synthesize aminoacyl-tRNA from cognate amino acid and tRNA. This ligase reaction proceeds through an activated aminoacyl-adenylate (aa-AMP). Common themes among AARSs include use of induced fit to drive catalysis and transition state stabilization by class-conserved sequence and structure motifs. Active site metal ions contribute to the amino acid activation step, while amino acid transfer to tRNA is generally a substrate-assisted concerted mechanism. A distinction between classes is the rate-limiting step for aminoacylation. We present some examples for each aspect of aminoacylation catalysis, including the experimental approaches developed to address questions of AARS chemistry.

Simple yet functional phosphate-loop proteins

Abundant and essential motifs, such as phosphate-binding loops (P-loops), are presumed to be the seeds of modern enzymes. The Walker-A P-loop is absolutely essential in modern NTPase enzymes, in mediating binding, and transfer of the terminal phosphate groups of NTPs. However, NTPase function depends on many additional active-site residues placed throughout the protein's scaffold. Can motifs such as P-loops confer function in a simpler context? We applied a phylogenetic analysis that yielded a sequence logo of the putative ancestral Walker-A P-loop element: a β-strand connected to an α-helix via the P-loop. Computational design incorporated this element into de novo designed β-α repeat proteins with relatively few sequence modifications. We obtained soluble, stable proteins that unlike modern P-loop NTPases bound ATP in a magnesium-independent manner. Foremost, these simple P-loop proteins avidly bound polynucleotides, RNA, and single-strand DNA, and mutations in the P-loop's key residues abolished binding. Binding appears to be facilitated by the structural plasticity of these proteins, including quaternary structure polymorphism that promotes a combined action of multiple P-loops. Accordingly, oligomerization enabled a 55-aa protein carrying a single P-loop to confer avid polynucleotide binding. Overall, our results show that the P-loop Walker-A motif can be implemented in small and simple β-α repeat proteins, primarily as a polynucleotide binding motif.

Augmenting the anisotropic network model with torsional potentials improves PATH performance, enabling detailed comparison with experimental rate data

PATH algorithms for identifying conformational transition states provide computational parameters-time to the transition state, conformational free energy differences, and transition state activation energies-for comparison to experimental data and can be carried out sufficiently rapidly to use in the "high throughput" mode. These advantages are especially useful for interpreting results from combinatorial mutagenesis experiments. This report updates the previously published algorithm with enhancements that improve correlations between PATH convergence parameters derived from virtual variant structures generated by RosettaBackrub and previously published kinetic data for a complete, four-way combinatorial mutagenesis of a conformational switch in Tryptophanyl-tRNA synthetase.

Combining multi-mutant and modular thermodynamic cycles to measure energetic coupling networks in enzyme catalysis

We measured and cross-validated the energetics of networks in Bacillus stearothermophilus Tryptophanyl-tRNA synthetase (TrpRS) using both multi-mutant and modular thermodynamic cycles. Multi-dimensional combinatorial mutagenesis showed that four side chains from this "molecular switch" move coordinately with the active-site Mg2+ ion as the active site preorganizes to stabilize the transition state for amino acid activation. A modular thermodynamic cycle consisting of full-length TrpRS, its Urzyme, and the Urzyme plus each of the two domains deleted in the Urzyme gives similar energetics. These dynamic linkages, although unlikely to stabilize the transition-state directly, consign the active-site preorganization to domain motion, assuring coupled vectorial behavior.

High-Dimensional Mutant and Modular Thermodynamic Cycles, Molecular Switching, and Free Energy Transduction

Understanding how distinct parts of proteins produce coordinated behavior has driven and continues to drive advances in protein science and enzymology. However, despite consensus about the conceptual basis for allostery, the idiosyncratic nature of allosteric mechanisms resists general approaches. Computational methods can identify conformational transition states from structural changes, revealing common switching mechanisms that impose multistate behavior. Thermodynamic cycles use factorial perturbations to measure coupling energies between side chains in molecular switches that mediate shear during domain motion. Such cycles have now been complemented by modular cycles that measure energetic coupling between separable domains. For one model system, energetic coupling between domains has been shown to be quantitatively equivalent to that between dynamic side chains. Linkages between domain motion, switching residues, and catalysis make nucleoside triphosphate hydrolysis conditional on domain movement, confirming an essential yet neglected aspect of free energy transduction and suggesting the potential generality of these studies.

Coding of Class I and II Aminoacyl-tRNA Synthetases

The aminoacyl-tRNA synthetases and their cognate transfer RNAs translate the universal genetic code. The twenty canonical amino acids are sufficiently diverse to create a selective advantage for dividing amino acid activation between two distinct, apparently unrelated superfamilies of synthetases, Class I amino acids being generally larger and less polar, Class II amino acids smaller and more polar. Biochemical, bioinformatic, and protein engineering experiments support the hypothesis that the two Classes descended from opposite strands of the same ancestral gene. Parallel experimental deconstructions of Class I and II synthetases reveal parallel losses in catalytic proficiency at two novel modular levels-protozymes and Urzymes-associated with the evolution of catalytic activity. Bi-directional coding supports an important unification of the proteome; affords a genetic relatedness metric-middle base-pairing frequencies in sense/antisense alignments-that probes more deeply into the evolutionary history of translation than do single multiple sequence alignments; and has facilitated the analysis of hitherto unknown coding relationships in tRNA sequences. Reconstruction of native synthetases by modular thermodynamic cycles facilitated by domain engineering emphasizes the subtlety associated with achieving high specificity, shedding new light on allosteric relationships in contemporary synthetases. Synthetase Urzyme structural biology suggests that they are catalytically-active molten globules, broadening the potential manifold of polypeptide catalysts accessible to primitive genetic coding and motivating revisions of the origins of catalysis. Finally, bi-directional genetic coding of some of the oldest genes in the proteome places major limitations on the likelihood that any RNA World preceded the origins of coded proteins.

An Ancestral Tryptophanyl-tRNA Synthetase Precursor Achieves High Catalytic Rate Enhancement without Ordered Ground-State Tertiary Structures

Urzymes-short, active core modules derived from enzyme superfamilies-prepared from the two aminoacyl-tRNA synthetase (aaRS) classes contain only the modules shared by all related family members. They have been described as models for ancestral forms. Understanding them currently depends on inferences drawn from the crystal structures of the full-length enzymes. As aaRS Urzymes lack much of the mass of modern aaRS's, retaining only a small portion of the hydrophobic cores of the full-length enzymes, it is desirable to characterize their structures. We report preliminary characterization of (15)N tryptophanyl-tRNA synthetase Urzyme by heteronuclear single quantum coherence (HSQC) NMR spectroscopy supplemented by circular dichroism, thermal melting, and induced fluorescence of bound dye. The limited dispersion of (1)H chemical shifts (0.5 ppm) is inconsistent with a narrow ensemble of well-packed structures in either free or substrate-bound forms, although the number of resonances from the bound state increases, indicating a modest, ligand-dependent gain in structure. Circular dichroism spectroscopy shows the presence of helices and evidence of cold denaturation, and all ligation states induce Sypro Orange fluorescence at ambient temperatures. Although the term "molten globule" is difficult to define precisely, these characteristics are consistent with most such definitions. Active-site titration shows that a majority of molecules retain ∼60% of the transition state stabilization free energy observed in modern synthetases. In contrast to the conventional view that enzymes require stable tertiary structures, we conclude that a highly flexible ground-state ensemble can nevertheless bind tightly to the transition state for amino acid activation.

Enhanced amino acid selection in fully evolved tryptophanyl-tRNA synthetase, relative to its urzyme, requires domain motion sensed by the D1 switch, a remote dynamic packing motif

We previously showed (Li, L., and Carter, C. W., Jr. (2013) J. Biol. Chem. 288, 34736-34745) that increased specificity for tryptophan versus tyrosine by contemporary Bacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) over that of TrpRS Urzyme results entirely from coupling between the anticodon-binding domain and an insertion into the Rossmann-fold known as Connecting Peptide 1. We show that this effect is closely related to a long range catalytic effect, in which side chain repacking in a region called the D1 Switch, accounts fully for the entire catalytic contribution of the catalytic Mg(2+) ion. We report intrinsic and higher order interaction effects on the specificity ratio, (kcat/Km)Trp/(kcat/Km)Tyr, of 15 combinatorial mutants from a previous study (Weinreb, V., Li, L., and Carter, C. W., Jr. (2012) Structure 20, 128-138) of the catalytic role of the D1 Switch. Unexpectedly, the same four-way interaction both activates catalytic assist by Mg(2+) ion and contributes -4.4 kcal/mol to the free energy of the specificity ratio. A minimum action path computed for the induced-fit and catalytic conformation changes shows that repacking of the four residues precedes a decrease in the volume of the tryptophan-binding pocket. We suggest that previous efforts to alter amino acid specificities of TrpRS and glutaminyl-tRNA synthetase (GlnRS) by mutagenesis without extensive, modular substitution failed because mutations were incompatible with interdomain motions required for catalysis.

Synthetic and editing mechanisms of aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases (aaRS) ensure the faithful transmission of genetic information in all living cells. The 24 known aaRS families are divided into 2 structurally distinct classes (class I and class II), each featuring a catalytic domain with a common fold that binds ATP, amino acid, and the 3'-terminus of tRNA. In a common two-step reaction, each aaRS first uses the energy stored in ATP to synthesize an activated aminoacyl adenylate intermediate. In the second step, either the 2'- or 3'-hydroxyl oxygen atom of the 3'-A76 tRNA nucleotide functions as a nucleophile in synthesis of aminoacyl-tRNA. Ten of the 24 aaRS families are unable to distinguish cognate from noncognate amino acids in the synthetic reactions alone. These enzymes possess additional editing activities for hydrolysis of misactivated amino acids and misacylated tRNAs, with clearance of the latter species accomplished in spatially separate post-transfer editing domains. A distinct class of trans-acting proteins that are homologous to class II editing domains also perform hydrolytic editing of some misacylated tRNAs. Here we review essential themes in catalysis with a view toward integrating the kinetic, stereochemical, and structural mechanisms of the enzymes. Although the aaRS have now been the subject of investigation for many decades, it will be seen that a significant number of questions regarding fundamental catalytic functioning still remain unresolved.

A master switch couples Mg²⁺-assisted catalysis to domain motion in B. stearothermophilus tryptophanyl-tRNA Synthetase

We demonstrate how tryptophanyl-tRNA synthetase uses conformation-dependent Mg(2+) activation to couple catalysis of tryptophan activation to specific, functional domain movements. Rate acceleration by Mg(2+) requires ∼-6.0 kcal/mol in protein⋅Mg(2+) interaction energy, none of which arises from the active site. A highly cooperative interaction between Mg(2+) and four residues from a remote, conserved motif that mediates the shear of domain movement (1) destabilizes the pretransition state conformation, thereby (2) inducing the Mg(2+) to stabilize the transition state for k(cat) by ∼-5.0 kcal/mol. Cooperative, long-range conformational effects on the metal therefore convert an inactive Mg(2+) coordination into one that can stabilize the transition state if, and only if, domain motion occurs. Transient, conformation-dependent Mg(2+) activation, analogous to the escapement in mechanical clocks, explains vectorial coupling.

Tryptophanyl-tRNA synthetase Urzyme: a model to recapitulate molecular evolution and investigate intramolecular complementation

We substantiate our preliminary description of the class I tryptophanyl-tRNA synthetase minimal catalytic domain with details of its construction, structure, and steady-state kinetic parameters. Generating that active fragment involved deleting 65% of the contemporary enzyme, including the anticodon-binding domain and connecting peptide 1, CP1, a 74-residue internal segment from within the Rossmann fold. We used protein design (Rosetta), rather than phylogenetic sequence alignments, to identify mutations to compensate for the severe loss of modularity, thus restoring stability, as evidenced by renaturation described previously and by 70-ns molecular dynamics simulations. Sufficient solubility to enable biochemical studies was achieved by expressing the redesigned Urzyme as a maltose-binding protein fusion. Michaelis-Menten kinetic parameters from amino acid activation assays showed that, compared with the native full-length enzyme, TrpRS Urzyme binds ATP with similar affinity. This suggests that neither of the two deleted structural modules has a strong influence on ground-state ATP binding. However, tryptophan has 10(3) lower affinity, and the Urzyme has comparably reduced specificity relative to the related amino acid, tyrosine. Molecular dynamics simulations revealed how CP1 may contribute significantly to cognate amino acid specificity. As class Ia editing domains are nested within the CP1, this finding suggests that this module enhanced amino acid specificity continuously, throughout their evolution. We call this type of reconstructed protein catalyst an Urzyme (Ur prefix indicates original, primitive, or earliest). It establishes a model for recapitulating very early steps in molecular evolution in which fitness may have been enhanced by accumulating entire modules, rather than by discrete amino acid sequence changes.

Mg2+-assisted catalysis by B. stearothermophilus TrpRS is promoted by allosteric effects

Mn(2+)-assisted catalysis by B. stearothermophilus TrpRS parallels that in polymerases and reduces specificity in amino acid activation. As predicted by nonequilibrium molecular dynamics simulations, multivariant thermodynamic cycles with [ATP]-dependent Michaelis-Menten kinetics and Mn(2+) for Mg(2+) substitution demonstrate energetic coupling of ATP affinities to the metal; to lysines K111 and K192, which interact via the PPi leaving group; and to K195, which couples differently to the metal via the alpha-phosphate. However, net coupling to the metal opposes catalysis in both ground (K(M)) and transition (k(cat)) states. The 10(5)-fold rate acceleration by Mg(2+)-protein interactions therefore requires additional favorable protein-metal couplings. Examples include longer-range, i.e., allosteric, interactions previously illustrated by the remote F37I mutation, which both reduces k(cat) and enhances catalytic assist by Mn(2+), relative to that by Mg(2+). These data support a model linking metal-assisted phosphoryl transfer catalysis to domain movement, and hence to free-energy transduction in a broad range of enzymes.

A conformational transition state accompanies tryptophan activation by B. stearothermophilus tryptophanyl-tRNA synthetase

B. stearothermophilus tryptophanyl-tRNA synthetase catalysis proceeds via high-energy protein conformations. Unliganded MD trajectories of the pretransition-state complex with Mg(2+)ATP and the (post) transition-state analog complex with adenosine tetraphosphate relax rapidly in opposite directions, the former regressing, the latter progressing along the structural reaction coordinate. The two crystal structures (rmsd 0.7 A) therefore lie on opposite sides of a conformational free-energy maximum as the chemical transition state forms. SNAPP analysis illustrates the complexity of the associated long-range conformational coupling. Switching interactions in four nonpolar core regions are locally isoenergetic throughout the transition. Different configurations, however, propagate their effects to unfavorable, longer-range interactions at the molecular surface. Designed mutation shows that switching interactions enhance the rate, perhaps by destabilizing the ground state immediately before the transition state and limiting nonproductive diffusion before and after the chemical transition state, thereby reducing the activation entropy. This paradigm may apply broadly to energy-transducing enzymes.