# 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

Violetta Weinreb, Li Li, Srinivas Niranj Chandrasekaran, Patrice Koehl, Marc Delarue, Charles Carter

► To cite this version:

Violetta Weinreb, Li Li, Srinivas Niranj Chandrasekaran, Patrice Koehl, Marc Delarue, et al. 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. Journal of Biological Chemistry, 2014, 289 (7), pp.4367-4376. 10.1074/jbc.M113.538660. pasteur-04079742

HAL Id: pasteur-04079742

https://pasteur.hal.science/pasteur-04079742

Submitted on 24 Apr 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# 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*§

Received for publication, November 27, 2013, and in revised form, December 27, 2013 Published, JBC Papers in Press, January 6, 2014, DOI 10.1074/jbc.M113.538660

Violetta Weinreb†, Li Li‡, Srinivas Niranj Chandrasekaran‡, Patrice Koehl§, Marc Delarue¶, and Charles W. Carter, Jr.†‡

From the ‡Department of Biochemistry and Biophysics, CB 7260, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599-7260, the §Department of Computer Science, University of California, Davis, California 95616, and the ¶Dynamique Structurale des Macromolecules, Unité de Mixte de Recherche 3528 du Centre National de la Recherche Scientifique, Pasteur Institute, 75015 Paris, France

**Background:** Amino acid selection by tryptophanyl-tRNA synthetase (TrpRS) requires intermodular coupling.

**Results:** Dynamic repacking of four side chains increases amino acid specificity 500-fold in contemporary TrpRS by reducing pocket size near the transition state.

**Conclusion:** An ancient tertiary packing motif not only activates the catalytic Mg<fnref n="2" />⁺ ion during catalysis, but also determines cognate amino acid specificity.

**Significance:** Allosteric enforcement of specificity increases robustness to mutation.

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<fnref n="2" />⁺ ion. We report intrinsic and higher order interaction effects on the specificity ratio, $(k_{cat}/K_m)_{Trp}/(k_{cat}/K_m)_{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<fnref n="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.

A major unplumbed wonder of protein structure is the amazing extent to which different parts of a protein or enzyme

respond appropriately to what happens elsewhere. The nearly simultaneous behavior of multiple parts is called cooperativity, and is a quintessential enabling characteristic that makes proteins such flexible and powerful catalysts. The drive to understand what “cooperativity” actually means, starting with the work of Perutz (1, 2) on the allosteric behavior of hemoglobin, has driven the greatest advances in enzymology and protein science.

Allostery, and linkage are closely related, but distinct concepts. Allostery refers to the intramolecular communication responsible for cooperativity; linkage refers to the detailed energetic coupling underlying such mechanisms. Non-additivity of perturbations to linked modules or residues, A and B, measures the energetic coupling of the interaction, A × B.

*Bacillus stearothermophilus* tryptophanyl-tRNA synthetase (TrpRS)<fnref n="2" /> has recently opened important new windows on cooperative behavior. We began these studies hoping to identify specific mechanisms responsible for the very large non-additivity (~-6.0 kcal/mol) with which Mg<fnref n="2" />⁺ accelerates TrpRS-catalyzed tryptophan activation. The fact that the metal accelerates tryptophan activation 10⁵-fold in the presence of TrpRS, but less than 10-fold in water (3), implies that metal-protein coupling is responsible for virtually the entire catalytic effect of Mg<fnref n="2" />⁺. Yet, the catalytic metal ion interacts only with ATP and water molecules, and is not bound by any protein residues.

Combinatorial mutagenesis of active-site lysine residues physically coupled to the metal in crystal structures showed, by elimination, that this effect must arise from outside the active

<figure><img src="image_2.png" /></figure>

<footnote>* This work was supported, in whole or in part, by National Institutes of Health Grant 40906 from the NIGMS.</footnote>

<footnote>§ This article contains supplemental Movie S1.</footnote>

<footnote n="1">¹ To whom correspondence should be addressed. Tel.: 919-966-3263; Fax: 919-966-2852; E-mail: carter@med.unc.edu.</footnote>

<footnote n="2">² The abbreviations used are: TrpRS, tryptophanyl-tRNA synthetase; GlnRS, glutaminyl-tRNA synthetase; PheRS, phenylalanyl-tRNA synthetase; ABD, anticodon-binding domain; CP1, connecting peptide 1; PreTS, pre-transition state; aaRS, aminoacyl-tRNA synthetase.</footnote>
Allosteric Enhancement of Amino Acid Specificity in TrpRS

site (4). Instead, we identified a long range allosteric influence
(5) from a widely distributed and highly conserved core-pack-
ing motif (6) remote from the active site that accounted entirely
for the catalytic assist by Mg²⁺. This motif is the most extensive
of four discrete locations where Delaunay tessellation patterns
change during TrpRS catalysis, suggesting the name "D1
Switch" (7). Repacking within this motif mediates shear devel-
oped by relative movement of the essentially rigid Rossmann-
fold and anticodon-binding domains (ABD) during catalysis.

We recently showed that all three separate TrpRS functions:
catalysis of amino acid activation, specific amino acid recogni-
tion, and tRNA acylation, depend in quite similar ways on
dynamic interactions between the active site, the ABD, and an
insertion into the Rossmann-fold called Connecting Peptide 1
(CP1). We restored the CP1 and ABD modules individually to
the TrpRS Urzyme, a functional representation of ancestral
Class I aaRS, showing that the two restored domains are linked
by a quantitatively similar coupling energy (about −6 kcal/mol)
in both selection of tryptophan *versus* tyrosine and tRNA ami-
noacylation (8).

Previously (9) we had attempted to change the specificity
ratio to favor activation of tyrosine using mutagenesis of resi-
dues whose direct interaction with substrate tryptophan sug-
gested that they were responsible for discriminating against
tyrosine. That effort failed: mutations weakened catalytic pro-
ficiency, without favoring tyrosine activation.

Here we show instead that combinatorial mutagenesis of D1
Switch residues affects the specificity ratio, (kcat/Km)Trp/(kcat/
Km)Tyr, sufficiently to account for the increased specificity of
full-length TrpRS. Moreover, we illustrate the sequential struc-
tural changes involved in long range coupling by constructing a
minimum action path (10) to connect the three crystal struc-
tures that represent stable states along the path. Repacking of
the four mutated D1 Switch residues occurs twice, coincident
with the transition states for induced-fit and catalysis. Each
change is followed by reductions in the volume of the trypto-
phan binding pocket, consistent with the experimental speci-
ficity data.

The quantitative agreement of the simulated most probable
path with the modular (8) and point mutant coupling energies
extends our analysis of cooperativity, allostery, and linkage in
TrpRS enzymology to unprecedented levels of detail. As long
range interaction between side chains is responsible for amino
acid specificity, no obvious pattern of point mutations can actu-
ally change specificity unless they are also compatible with the
necessary domain movement, making specific amino acid
selection between competing canonical amino acids more
robust to mutation.

EXPERIMENTAL PROCEDURES

*Mutagenesis and Protein Purification*—Mutations were con-
structed using GeneTailor (Invitrogen). The expression plas-
mid pET11a containing WT TrpRS was methylated by DNA
methylase and used as template for mutagenic PCR. Primers
were designed for both directions, with a 15–20-nucleotide
overlap for efficient circularization. The PCR product was used
to transform DH5α *Escherichia coli* cells, and plated on LB

plates with ampicillin. Plasmids were sequenced to confirm the
mutations.

Mutant proteins were expressed in *E. coli* BL21(DE3)pLysS, with both ampicillin and chloramphenicol. Simultaneous 3-lister cultures of two mutants were resuspended in 20 mM HEPES, 0.3 M NaCl, 10 mM 2-mercaptoethanol, 30 mM imidazole (pH 7.6) and rapidly frozen. The cell paste was sonicated five times for 20 s with a Fisher Scientific Sonic Dismembrator and cleared by centrifugation at 27,800 × g for 30 min at 4 °C. Nickel-nitrilotriacetic acid beads (HisPur Ni-NTA Resin, Thermo Scientific; 1 ml) were added to the supernatant, followed by 1 h of shaking at 4 °C. After centrifugation the beads were washed with 20 ml of lysis buffer and shaken 10 min at 4 °C. Bound enzyme was eluted three times with lysis buffer plus 0.3 M imidazole. Purified mutant proteins were cleaved with a 1:10 ratio of tobacco etch virus protease overnight at 4 °C, diluted to reduce the imidazole concentration below 30 mM, and passed over the nickel resin again to bind uncleaved protein and tobacco etch virus. The eluate was concentrated using an Amicon PM10 Ultra membrane and stored at −20 °C in 50% glycerol. All proteins were >95% homogeneous as judged by gel electrophoresis.

*Active-site Titration*—The reaction mixture contained 50 mM HEPES (pH 7.5), 10 mM MgCl₂, 100 mM KCl, 1 mM DTT, 10 µM [γ-³²P]ATP, 2.5 mM Trp, 10 units/ml of inorganic pyrophosphatase. We added 3 µM enzyme and incubated the mixture at 37 °C. At specific time points the reaction was quenched by adding 3 µl of the reaction mixture to 6 µl of 400 mM sodium acetate (pH 5.0), 0.1% SDS. 3 µl of this mixture was spotted on a thin layer chromatography plate. To separate [³²P]ATP from [³²P]PPᵢ we pre-soaked and eluted the TLC plates with 750 mM KH₂PO₄ (pH 3.5), 4 M urea buffer (11). The TLC plate was dried, exposed 15–30 min, depending on the variant, and scanned on a Typhoon scanner and quantified with ImageJ software. All assays included data points well into the steady state (linear with time).

PreTS crystal structures for ~8 mutant proteins including the
quadruple mutant reveal very little structural variance from the
native structure, Protein Data Bank code 1MAU (12). Active titers
for the various mutant enzymes range from 0.3 to 0.95 with a mean
value of 0.51 ± 0.17. The structural and enzymatic evidence that
the mutant enzymes are properly folded is thus very strong.

*Michaelis-Menten Kinetics*—PPᵢ-exchange assays were done at 37 °C and initiated with 10 µl of enzyme to 190 µl of assay mixture: 0.1 M Tris-Cl, 0.01 M potassium fluoride, 5 mM MgCl₂, 2 mM ATP, 10 mM 2-mercaptoethanol (pH 8.0) plus 2 mM ³²P-PPᵢ at a specific radioactivity between 1 × 10⁵ and 2 × 10⁵ cpm/mol. Varying enzyme concentrations (100–400 µM) and incubation times (15 min for Trp, 20 min for Tyr) were used, depending on the activity level, which was determined by range finding for each variant. Michaelis-Menten kinetics were examined by varying the amino acid concentrations (0.01, 0.05, 0.5, 1.0, 2.0, 5.0, 10.0, 30, 60, 90, and 150 µM tryptophan and 0.005, 0.01, 0.05, 0.1, 0.15, 0.25, 0.5, 0.75, 1.0, 1.5, and 2.5 mM tyrosine). Assays were performed on 96-well plates with a Beckman-Coulter Biomek 3000 automatic liquid-handling work station, using Whatman filter plates and a Promega Vacman vacuum manifold for filtrations. [³²P]ATP was eluted from charcoal with pyridine as described (19). For metal substitution assays the mixtures were made without Mg²⁺ and treated with Chelex

<figure><img src="image_3.png" /></figure>
100 for 30 min at 4 °C to remove trace metals and supplemented with 0.9 mM MnCl₂. All experiments were repeated two or more times with four replicates, to improve the accuracy of regression coefficients.

*Multidimensional Thermodynamic Cycles*—The formal similarity between thermodynamic cycles (13) and factorial experimental design (14) can help clarify how this approach works. Briefly, as in a full factorial experiment, a multidimensional thermodynamic cycle measures the impact of each perturbation from a set in the context of all other perturbations in that set. This means that at least as many protein variants are examined as there are vertices in the corresponding thermodynamic cycle. In general, many or all of the variants are examined multiple times to estimate the experimental variance. For the fifth-dimensional thermodynamic cycle in Ref. 5, each of four sites was mutated to a single other amino acid, and the resulting 15 variants were assayed multiple times with both Mg²⁺ and Mn²⁺. We used the same mutants here.

Key to the relationship between these experimental results and the coupling energies is the *experimental design matrix*. Each experiment corresponds to a row in this matrix, giving the state of each perturbation and its higher-order interactions. For a double mutant cycle for sites A and B, the row corresponding to variant A has the form [0, 1, 0] because variant A is a mutant at site A (labeled 0), wild type at site B (labeled 1), and has 0 also in the third column, which is the product of columns A and B and corresponds to the A × B interaction.

Without experimental data, contributions of wild-type residues and their interaction energies to system function are unknown. Multiple column vectors correspond to experimental values for dependent variables. Examples are $\Delta G_{k_{\text{cat}}}$ and $\Delta G_{K_M}$ from Michaelis-Menten kinetics. Free energies are more useful than rate and equilibrium constants, because they are additive, and hence can be used to construct linear models from the state columns of the design matrix. These additional column vectors are pre-pended to the design matrix.

The final element of the analysis is a column vector of coefficients, which describe how much of a given predictor is necessary for optimal agreement between the observed data and corresponding results calculated from the linear model. Pre-multiplying the coefficient vector by the design matrix gives a set of simultaneous linear equations specifying the calculated results. If the number of experiments is greater than the number of unknowns in these simultaneous equations, then the solution and error estimates are found by least squares.

The algebraic formulation of the previous paragraph should clarify that the highest-order mutant is but one data point of the many that need to be summed together to estimate the coefficients, emphasizing the crucial difference between the experimental result for any given variant, and the coefficient accorded to the corresponding intrinsic or interaction effect by the linear equations. Just as the calculated data results all must be summed for each coefficient, coefficients must be estimated from all experimental points. In particular, it is overwhelmingly tempting to infer the coupling energy for the highest order interaction from the observed behavior of the highest order mutant (i.e. the quadruple mutant in a four-way experiment). The two are, in fact, almost totally unrelated.

*Minimum Action Path Simulation (10)*—A most probable trajectory connecting ground state crystal structures of *B. stearothermophilus* tryptophanyl-tRNA synthetase during induced-fit and subsequent catalysis was simulated from elastic network models (10) using a version of the program PATH (developed by P. Koehl) that finds the path with minimum action that connects three states. Spring constants for elastic network models of all states were derived using the MinAction-Path server.<fnref n="3" /> Total time intervals allocated to observe transitions, and energy differences between successive states were determined from the convergence surface.<fnref n="4" />

The PATH program carries out a rapid search for the transition state conformation between crystallographic states by assuring structural, velocity, and energy continuity to minimize the total Onsager-Machlup action of the path. Minimum action therefore assures the most probable path.

The sequence of structures was computed from the normal modes of the initial and final states of the two separate transitions, which are each given 100 steps, for a total of 200 steps. Thus the midpoint determines the transition from the induced-fit domain movement to that which takes place during the catalytic step. The two domain motions can be identified by the direction in which the ABD moves: during the induced-fit transition, it moves counterclockwise in the plane of the animation, during catalysis it rotates outward toward the viewer.

Colors used in supplemental trajectory Movie S1 are as follows: the Urzyme is slate blue, CP1 is forest green, the ABD is wheat. The β-strand forming the bottom of the tryptophan-binding pocket is yellow and the specificity helix forming the top is orange, as in Fig. 2A. Side chains forming the tryptophan binding pocket are hot pink; D1 Switch residue Ile-4 is yellow, and residues Phe-26, Tyr-33, and Phe-37 are deep blue.

## RESULTS

*Modular TrpRS Construction and the Conformational Cycle*—Crystallographic structures suggest that catalysis of amino acid activation entails three distinct TrpRS structural states: open, pre-transition state (PreTS), and products (Prod). Transition between the first two structures corresponds to induced-fit, between the second and third to catalysis, and the restorative step to product release stimulated by acyl-transfer to cognate tRNA<sup>Trp</sup>. Superposition of the three representative crystal structures (Protein Data Bank codes 1MB2, 1MAU, and 1I6L) reveals that the conformational cycle (7, 12, 16, 17) is well approximated by rigid-body movements of the ABD relative to the active site (Fig. 1).

Long range catalytic activation of the Mg²⁺ ion by the D1 Switch (5) suggested that it mediates coupling between catalytic and ABD domains. The TrpRS Urzyme (18, 19) enabled us to examine this intrinsic modularity directly (8). The Urzyme contains a functional active site without either CP1 or the ABD. It therefore likely represents an ancient Class I aminoacyl-tRNA synthetase (aaRS). The ABD and CP1 are presumably more modern modules (Fig. 2) that afford a mechanism to couple

<figure><img src="image_4.png" /></figure>

<footnote n="3">³ M. Delarue, MinActionPath, unpublished data.</footnote>

<footnote n="4">⁴ S. N. Chandrasekaran, P. Koehl, M. Delarue, S. Doniach, and C. W. Carter, Jr., manuscript in preparation.</footnote>
Allosteric Enhancement of Amino Acid Specificity in TrpRS

domain movement (Fig. 1) to structural changes within the active site (Fig. 2A). Energetic coupling (about −6 kcal/mol) between these two modules enhances all three TrpRS functions, Mg²⁺-dependent rate acceleration, specificity, and tRNA<sup>Trp</sup> aminoacylation to approximately the same degree (8).

The tryptophan binding site lies between two disjointed active site segments of the Urzyme; the N-terminal crossover and the specificity helix are separated by the CP1 module. Moreover, the CP1 module itself forms an exoskeleton that wraps around these two fragments, as suggested in Fig. 2A. Notably, the D1 Switch itself is a protoallosteric site contained entirely within the Urzyme (Fig. 2B). In this work, we extend the combinatorial mutagenesis of the D1 Switch and multidimensional thermodynamic cycles to the analysis of specificity during amino acid activation.

*Multidimensional Thermodynamic Analysis of TrpRS Specificity—*
Michaelis-Menten parameters for all experiments are summarized in Tables 1 and 2, and illustrated graphically in Fig. 3. The three plots show the correlation between observed free energies associated with the three parameters, *K*<sub>m</sub>, *k*<sub>cat</sub>, and *k*<sub>cat</sub>/K<sub>m</sub>, and those calculated using appropriate linear combinations of mutant sites, which are the independent variables responsible for systematic variation. *F*-ratios for the three models are highly significant: Δ*GK*<sub>M</sub>, *F* = 6.4 (*p* < 0.0001); Δ*Gk*<sub>cat</sub>, *F* = 188 (*p* < 0.0001); Δ*Gk*<sub>cat</sub>/K<sub>m</sub>, *F* = 33 (*p* < 0.0001). Distributions of Δ*GK*<sub>M</sub> for the two amino acids (Fig. 3A) are coextensive. However, plots for Δ*Gk*<sub>cat</sub> and Δ*Gk*<sub>cat</sub>/K<sub>m</sub> both clearly separate measurements for the two amino acids. They illustrate that specificity is driven largely by differences in *k*<sub>cat</sub> for the two amino acids (Fig. 3, *B* and *C*). Tryptophan, on average, is associated with a *k*<sub>cat</sub> value 220 times that observed with tyrosine (Δ*G*<sub>W</sub> = −3.2 kcal/mol).

The specificity ratios (Fig. 4) are computed from two differ-
ent experimental determinations, amplifying their errors. For
this reason, they are subject to higher experimental uncertain-
ties and require extensive replication to achieve suitable estimates
of variance. In view of this inherent experimental noise, we present
in Tables 1 and 2 statistical summaries of the Michaelis-Menten
parameters from which the specificity ratios in Table 3 were com-
puted. The important question can be stated succinctly: have we
replicated our measurements sufficiently to estimate confidently
the most important coupling energies?

<figure><img src="image_5.png" /></figure>

**FIGURE 1.** Relative domain motion during the TrpRS catalytic cycle. The Rossmann-fold and ABD move as rigid bodies relative to one another (12, 14, 17). Hinge bending and torsional motion of the anticodon-binding domain shown here reveal that the ABD assumes three distinct orientations relative to the Rossmann-fold domains that have been superimposed. These orientations are denoted LF (ligand-free, *wheat*), PreTS (pre-transition state; *sky blue*), and Products (*magenta*). Catalytic residues move with the ABD (11, 17).

<figure><img src="image_6.png" /></figure>

**FIGURE 2.** **Structural basis of allosteric behavior during amino acid activation by TrpRS.** *A*, modularity of the TrpRS monomer. The Urzyme contains three modules: the N-(yellow) and C (blue)-terminal crossover connections and the specificity helix (orange). It interacts snugly with both the CP1 (magenta) and anticodon-binding (teal) modules. Notably, the CP1 module circumscribes the Urzyme, constraining the relative movement of the specificity helix, relative to the N-terminal crossover connection of the Rossmann-fold, which determines the volume of the amino acid binding pocket. *B*, the main energetic communication within the Urzyme involves the D1 Switch. Four residues (Ile-4, Phe-26, Tyr-33, and Phe-37) were examined previously and are the subject of this work. The intact monomer is represented by the same colors as in *A*, except that the Urzyme is entirely *blue* and has a *transparent surface*, through which substrates tryptophan (W) and ATP (T), and D1 Switch residues are visible. The D1 Switch centroid is ~18 Å from the ATP α-phosphate and ~14 Å from the tryptophan α-carbon.
<table><caption>**TABLE 1**
Statistical properties of Michaelis-Menten parameters for TrpRS D1 Switch variants</caption><thead><tr><th></th><th>K<sub>m</sub></th><th>k<sub>cat</sub></th><th>k<sub>cat</sub>/K<sub>m</sub></th><th>ΔGk<sub>cat</sub>/K<sub>m</sub></th></tr><tr><th></th><th>M</th><th>S</th><th></th><th></th></tr></thead><tbody><tr><td><strong>Trp, Mg<sup>2+</sup></strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean</td><td>9.43E-06</td><td>5.45</td><td>1.09E+06</td><td>-8.23</td></tr><tr><td>S.D.</td><td>8.45E-06</td><td>3.46</td><td>1.00E+06</td><td>-8.18</td></tr><tr><td>Error<sub>Rel</sub></td><td>0.90</td><td>0.64</td><td>0.92</td><td>0.05</td></tr><tr><td>Min</td><td>4.24E-07</td><td>0.86</td><td>1.60E+05</td><td>-7.09</td></tr><tr><td>Max</td><td>3.65E-05</td><td>18.49</td><td>3.40E+06</td><td>-8.90</td></tr><tr><td>Ratio/difference</td><td>86.08</td><td>21.51</td><td>21</td><td>-1.81</td></tr><tr><td><strong>Trp, Mn<sup>2+</sup></strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean</td><td>5.27E-05</td><td>4.88</td><td>1.64E+05</td><td>-7.11</td></tr><tr><td>S.D.</td><td>4.01E-05</td><td>2.08</td><td>1.55E+05</td><td>-7.08</td></tr><tr><td>Error<sub>Rel</sub></td><td>0.76</td><td>0.43</td><td>0.95</td><td>0.03</td></tr><tr><td>Min</td><td>3.50E-06</td><td>1.03</td><td>2.09E+04</td><td>-5.89</td></tr><tr><td>Max</td><td>1.63E-04</td><td>8.16</td><td>7.41E+05</td><td>-8.00</td></tr><tr><td>Ratio/difference</td><td>46.57</td><td>7.90</td><td>35</td><td>-2.11</td></tr><tr><td><strong>Tyr, Mg<sup>2+</sup></strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean</td><td>4.35E-05</td><td>0.02</td><td>3.88E+03</td><td>-4.89</td></tr><tr><td>S.D.</td><td>7.21E-05</td><td>0.020</td><td>8.53E+03</td><td>-5.36</td></tr><tr><td>Error<sub>Rel</sub></td><td>1.66</td><td>0.94</td><td>2.20</td><td>-0.47</td></tr><tr><td>Min</td><td>1.30E-06</td><td>0.004</td><td>1.44E+02</td><td>-2.94</td></tr><tr><td>Max</td><td>3.57E-04</td><td>0.08</td><td>4.76E+04</td><td>-6.38</td></tr><tr><td>Ratio/difference</td><td>274.85</td><td>22.57</td><td>331</td><td>-3.43</td></tr><tr><td><strong>Tyr, Mn<sup>2+</sup></strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean</td><td>6.05E-05</td><td>0.042</td><td>3.06E+03</td><td>-4.75</td></tr><tr><td>S.D.</td><td>8.30E-05</td><td>0.040</td><td>5.20E+03</td><td>-5.07</td></tr><tr><td>Error<sub>Rel</sub></td><td>1.37</td><td>0.95</td><td>1.70</td><td>-0.31</td></tr><tr><td>Min</td><td>1.50E-06</td><td>0.009</td><td>1.01E+01</td><td>-1.37</td></tr><tr><td>Max</td><td>3.00E-04</td><td>0.238</td><td>2.62E+04</td><td>-6.02</td></tr><tr><td>Ratio/difference</td><td>200.00</td><td>26.72</td><td>2595</td><td>-4.65</td></tr></tbody></table>

<table><caption>**TABLE 2**
Statistical properties of specificity constants $k_{cat}/K_m(W)/k_{cat}/K_m(Y)$</caption><thead><tr><th></th><th>Specificity ratio</th><th>S.D.</th><th>Error<sub>rel</sub></th><th>ΔG<sub>spec</sub></th><th>ΔG<sub>error</sub></th></tr></thead><tbody><tr><td>Mean</td><td>926</td><td>830</td><td>0.8</td><td>3.5</td><td>-0.1</td></tr><tr><td>S.D.</td><td>1064</td><td>974</td><td>0.4</td><td>0.9</td><td>-0.6</td></tr><tr><td>Min</td><td>11</td><td>8</td><td>0.2</td><td>1.7</td><td>-1.0</td></tr><tr><td>Max</td><td>4192</td><td>3255</td><td>1.7</td><td>5.1</td><td>0.3</td></tr><tr><td>Ratio</td><td>388</td><td>430</td><td>9.6</td><td>3.0</td><td>1.3</td></tr></tbody></table>

The ultimate answer to this question lies in the statistics for the regression models. However, a simpler analysis of signal to noise from Tables 1 and 2 affords useful confirmation. Multiple determinations of the free energies for the apparent second-order rate constants in Table 1 are, on average, ~30 times the maximum relative error in multiple determinations of the same quantities. Similarly, the ratio of average specificity constant from Table 2 to the maximum relative error is 257; converted to free energies, this gives 8.6 for the range of free energies divided by the largest relative error for multiple estimates. These values imply that the range of specificities observed for the mutant proteins with Mg²⁺ and Mn²⁺ is substantially larger than the relative errors of the estimates, and hence that a significant signal can be extracted from the data, despite the increased relative error.

That conclusion is verified by regression modeling to estimate the strongest predictors of specificity from among the main and higher-order interaction effects (Table 4). Fifteen of 26 predictors have significant Student's *t* tests with *p* values <0.05. It should be noted that the *t*-statistic is the ratio of the coefficient to its standard error, so that *p* values are related to the impact of the experimental noise. In this way, the regression model apportions the experimental noise appropriately, thus highlighting the most important coefficients. Further evidence for the significance of these conclusions arises from the consistency of the implications summarized in the two following paragraphs with previously published results.

Two aspects of these coefficients stand out. First, by far the most important effect is the highest-order interaction between

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**FIGURE 3.** Linear models for Gibbs-free energies associated with Michaelis-Menten parameters summarized in Tables 1 and 2. Colors separate determinations for tryptophan (red) and tyrosine (blue). A, $ΔG_{K_m}$ distributions for the two amino acids overlap extensively. B, $ΔG_{k_{cat}}$ distributions for the two amino acids are entirely disjoint, because of the overwhelming effect of the amino acid. C, $ΔG_{k_{cat}/K_m}$ the distributions for $ΔG_{k_{cat}}$ are each broadened by the effect of $ΔG_{K_m}$ but remain dominated by the effect of the amino acid. Specificity is driven largely by $k_{cat}$.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**FIGURE 4.** Regression model for $ΔG_{spec}$, the free energy for selecting tryptophan versus tyrosine. Red open circles are for variants with the I4V mutation, blue closed circles are for variants with wild-type Ile-4.</figcaption></figure>

<figure><img src="image_11.png" /></figure>
# Allosteric Enhancement of Amino Acid Specificity in TrpRS

<table><caption>**TABLE 3**

Specificity ratios of TrpRS variants with Mg²⁺ and Mn²⁺</caption><thead><tr><th>Variant</th><th>k<sub>cat</sub>/K<sub>m</sub>(W)/k<sub>cat</sub>/K<sub>m</sub>(Y); Mg<sup>2+</sup></th><th>Variant</th><th>k<sub>cat</sub>/K<sub>m</sub>(W)/k<sub>cat</sub>/K<sub>m</sub>(Y); Mn<sup>2+</sup></th></tr></thead><tbody><tr><td>WT TrpRS</td><td>2192</td><td>WT TrpRS</td><td>1908</td></tr><tr><td>WT TrpRS</td><td>2289</td><td>WT TrpRS</td><td>587</td></tr><tr><td>WT TrpRS</td><td>4686</td><td>WT TrpRS</td><td>1996</td></tr><tr><td>WT TrpRS</td><td>4893</td><td>WT TrpRS</td><td>614</td></tr><tr><td>I4V</td><td>210</td><td>I4V</td><td>48</td></tr><tr><td>I4V</td><td>106</td><td>I4V</td><td>27</td></tr><tr><td>I4V</td><td>403</td><td>I4V</td><td>60</td></tr><tr><td>I4V</td><td>203</td><td>I4V</td><td>33</td></tr><tr><td>F26L</td><td>120</td><td>F26L</td><td>51</td></tr><tr><td>F26L</td><td>737</td><td>F26L</td><td>38</td></tr><tr><td>F26L</td><td>337</td><td>F26L</td><td>58</td></tr><tr><td>F26L</td><td>2065</td><td>F26L</td><td>43</td></tr><tr><td>Y33F</td><td>115</td><td>Y33F</td><td>62</td></tr><tr><td>Y33F</td><td>632</td><td>Y33F</td><td>50</td></tr><tr><td>Y33F</td><td>607</td><td>Y33F</td><td>278</td></tr><tr><td>Y33F</td><td>3346</td><td>Y33F</td><td>224</td></tr><tr><td>F37I</td><td>1629</td><td>Y33F</td><td>131</td></tr><tr><td>F37I</td><td>859</td><td>F37I</td><td>106</td></tr><tr><td>F37I</td><td>605</td><td>F37I</td><td>211</td></tr><tr><td>F37I</td><td>2831</td><td>F37I</td><td>188</td></tr><tr><td>F37I</td><td>1492</td><td>F37I</td><td>296</td></tr><tr><td>F37I</td><td>1051</td><td>F37I</td><td>263</td></tr><tr><td>I4V-F26L</td><td>4790</td><td>F37I</td><td>543</td></tr><tr><td>I4V-F26L</td><td>3714</td><td>F37I</td><td>483</td></tr><tr><td>I4V-F26L</td><td>772</td><td>I4V-F26L</td><td>260</td></tr><tr><td>I4V-F26L</td><td>8198</td><td>I4V-F26L</td><td>195</td></tr><tr><td>I4V-F26L</td><td>6356</td><td>I4V-F26L</td><td>288</td></tr><tr><td>I4V-F26L</td><td>1322</td><td>I4V-F26L</td><td>216</td></tr><tr><td>I4V-Y33F</td><td>2225</td><td>I4V-Y33F</td><td>72</td></tr><tr><td>I4V-Y33F</td><td>14</td><td>I4V-Y33F</td><td>48</td></tr><tr><td>I4V-Y33F</td><td>126</td><td>I4V-Y33F</td><td>54</td></tr><tr><td>I4V-Y33F</td><td>6916</td><td>I4V-Y33F</td><td>37</td></tr><tr><td>I4V-Y33F</td><td>43</td><td>I4V-Y33F</td><td>43</td></tr><tr><td>I4V-Y33F</td><td>393</td><td>I4V-Y33F</td><td>29</td></tr><tr><td>I4V-F37I</td><td>2025</td><td>I4V-F37I</td><td>481</td></tr><tr><td>I4V-F37I</td><td>389</td><td>I4V-F37I</td><td>63</td></tr><tr><td>I4V-F37I</td><td>6198</td><td>I4V-F37I</td><td>426</td></tr><tr><td>I4V-F37I</td><td>1192</td><td>I4V-F37I</td><td>56</td></tr><tr><td>F26L-Y33F</td><td>722</td><td>F26L-Y33F</td><td>145</td></tr><tr><td>F26L-Y33F</td><td>116</td><td>F26L-Y33F</td><td>268</td></tr><tr><td>F26L-Y33F</td><td>2714</td><td>F26L-Y33F</td><td>328</td></tr><tr><td>F26L-Y33F</td><td>128</td><td>F26L-Y33F</td><td>608</td></tr><tr><td>F26L-F37I</td><td>1290</td><td>F26L-Y33F</td><td>393</td></tr><tr><td>F26L-F37I</td><td>113</td><td>F26L-Y33F</td><td>728</td></tr><tr><td>F26L-F37I</td><td>2444</td><td>F26L-F37I</td><td>909</td></tr><tr><td>F26L-F37I</td><td>215</td><td>F26L-F37I</td><td>264</td></tr><tr><td>F26L-F37I</td><td></td><td>F26L-F37I</td><td>935</td></tr><tr><td>F26L-F37I</td><td></td><td>F26L-F37I</td><td>271</td></tr><tr><td>Y33F-F37I</td><td>1338</td><td>Y33F-F37I</td><td>2067</td></tr><tr><td>Y33F-F37I</td><td>111</td><td>Y33F-F37I</td><td>104</td></tr><tr><td>Y33F-F37I</td><td>757</td><td>Y33F-F37I</td><td>5487</td></tr><tr><td>Y33F-F37I</td><td>2647</td><td>Y33F-F37I</td><td>275</td></tr><tr><td>Y33F-F37I</td><td>220</td><td>Y33F-F37I</td><td>7883</td></tr><tr><td>Y33F-F37I</td><td>1499</td><td>Y33F-F37I</td><td>395</td></tr><tr><td>I4V-F26L-Y33F</td><td>52</td><td>I4V-F26L-Y33F</td><td>31</td></tr><tr><td>I4V-F26L-Y33F</td><td>61</td><td>I4V-F26L-Y33F</td><td>5</td></tr><tr><td>I4V-F26L-Y33F</td><td>144</td><td>I4V-F26L-Y33F</td><td>125</td></tr><tr><td>I4V-F26L-Y33F</td><td>170</td><td>I4V-F26L-Y33F</td><td>21</td></tr><tr><td>I4V-F26L-F37I</td><td>81</td><td>I4V-F26L-F37I</td><td>259</td></tr><tr><td>I4V-F26L-F37I</td><td>256</td><td>I4V-F26L-F37I</td><td>95</td></tr><tr><td>I4V-F26L-F37I</td><td>309</td><td>I4V-F26L-F37I</td><td>326</td></tr><tr><td>I4V-F26L-F37I</td><td>972</td><td>I4V-F26L-F37I</td><td>119</td></tr><tr><td>I4V-Y33F-F37I</td><td>136</td><td>I4V-Y33F-F37I</td><td>4</td></tr><tr><td>I4V-Y33F-F37I</td><td>450</td><td>I4V-Y33F-F37I</td><td>17</td></tr><tr><td>I4V-Y33F-F37I</td><td>1582</td><td>I4V-Y33F-F37I</td><td>4</td></tr><tr><td>I4V-Y33F-F37I</td><td>5231</td><td>I4V-Y33F-F37I</td><td>17</td></tr><tr><td>F26L-Y33F-F37I</td><td>57</td><td>F26L-Y33F-F37I</td><td>431</td></tr><tr><td>F26L-Y33F-F37I</td><td>45</td><td>F26L-Y33F-F37I</td><td>1988</td></tr><tr><td>F26L-Y33F-F37I</td><td>647</td><td>F26L-Y33F-F37I</td><td>858</td></tr><tr><td>F26L-Y33F-F37I</td><td>512</td><td>F26L-Y33F-F37I</td><td>3954</td></tr><tr><td>I4V-Y33F-F26L-F37I</td><td>69</td><td>I4V-F26L-Y33F-F37I</td><td>26</td></tr><tr><td>I4V-Y33F-F26L-F37I</td><td>58</td><td>I4V-F26L-Y33F-F37I</td><td>312</td></tr><tr><td>I4V-Y33F-F26L-F37I</td><td>447</td><td>I4V-F26L-Y33F-F37I</td><td>16</td></tr><tr><td>I4V-Y33F-F26L-F37I</td><td>375</td><td>I4V-F26L-Y33F-F37I</td><td>186</td></tr></tbody></table>

<figure><img src="image_12.png" /></figure>
TABLE 4

**Regression model for contributions to the free energy of the specificity ratio**

$R^2$ for this model is 0.65. The F-ratio is 8.6 with a *p* value of <0.0001. Coefficients are given in kcal/mol. The model is based on the 146 observations in Table 3, which afford 114 degrees of freedom for pure error. Most of the unexplained variance therefore arises because of the large experimental noise, as discussed under "Results."

<table><thead><tr><th>Predictor<sup>a</sup></th><th>Estimate</th><th>S.E.</th><th>t<br>Ratio</th><th>Prob>|t|</th></tr></thead><tbody><tr><td>Mg<sup>2+</sup></td><td>-0.73</td><td>0.11</td><td>-6.72</td><td>&lt;.0001</td></tr><tr><td>Ile-4</td><td>-0.71</td><td>0.11</td><td>-6.53</td><td>&lt;.0001</td></tr><tr><td>(Ile-4) × (Phe-26) × (Tyr-33) × (Phe-37)</td><td>-4.44</td><td>0.87</td><td>-5.10</td><td>&lt;.0001</td></tr><tr><td>Tyr-33</td><td>-0.44</td><td>0.11</td><td>-4.01</td><td>0.0001</td></tr><tr><td>(Mg<sup>2+</sup>) × (Ile-4)</td><td>0.86</td><td>0.22</td><td>3.97</td><td>0.0001</td></tr><tr><td>(Ile-4) × (Phe-26) × (Tyr-33)</td><td>-1.41</td><td>0.44</td><td>-3.23</td><td>0.0016</td></tr><tr><td>(Mg<sup>2+</sup>) × (Phe-26) × (Phe-37)</td><td>1.25</td><td>0.44</td><td>2.87</td><td>0.0049</td></tr><tr><td>(Mg<sup>2+</sup>) × (Phe-26)</td><td>-0.62</td><td>0.22</td><td>-2.87</td><td>0.0049</td></tr><tr><td>(Ile-4) × (Tyr-33)</td><td>0.61</td><td>0.22</td><td>2.80</td><td>0.0059</td></tr><tr><td>(Mg<sup>2+</sup>) × (Ile-4) × (Phe-37)</td><td>-1.13</td><td>0.43</td><td>-2.59</td><td>0.0108</td></tr><tr><td>(Mg<sup>2+</sup>) × (Ile-4) × (Tyr-33)</td><td>-1.08</td><td>0.43</td><td>-2.49</td><td>0.014</td></tr><tr><td>(Ile-4) × (Phe-26)</td><td>-0.52</td><td>0.22</td><td>-2.40</td><td>0.0179</td></tr><tr><td>(Mg<sup>2+</sup>) × (Phe-26) × (Tyr-33)</td><td>0.95</td><td>0.44</td><td>2.18</td><td>0.0316</td></tr><tr><td>(Mg<sup>2+</sup>) × (Ile-4) × (Tyr-33) × (Phe-37)</td><td>1.89</td><td>0.87</td><td>2.18</td><td>0.0316</td></tr><tr><td>(Ile-4) × (Phe-26) × (Phe-37)</td><td>-0.91</td><td>0.44</td><td>-2.09</td><td>0.0388</td></tr><tr><td>(Mg<sup>2+</sup>) × (Phe-37)</td><td>-0.40</td><td>0.22</td><td>-1.83</td><td>0.0704</td></tr><tr><td>Phe-37</td><td>0.18</td><td>0.11</td><td>1.67</td><td>0.0982</td></tr><tr><td>(Mg<sup>2+</sup>) × (Tyr-33)</td><td>-0.31</td><td>0.22</td><td>-1.44</td><td>0.1531</td></tr><tr><td>Phe-26</td><td>-0.15</td><td>0.11</td><td>-1.40</td><td>0.1631</td></tr><tr><td>(Ile-4) × (Tyr-33) × (Phe-37)</td><td>-0.49</td><td>0.44</td><td>-1.11</td><td>0.2675</td></tr><tr><td>(Tyr-33) × (Phe-37)</td><td>-0.21</td><td>0.22</td><td>-0.95</td><td>0.3462</td></tr><tr><td>(Phe-26) × (Tyr-33)</td><td>-0.11</td><td>0.22</td><td>-0.52</td><td>0.6046</td></tr><tr><td>(Mg<sup>2+</sup>) × (Tyr-33) × (Phe-37)</td><td>-0.17</td><td>0.43</td><td>-0.39</td><td>0.6958</td></tr><tr><td>(Phe-26) × (Tyr-33) × (Phe-37)</td><td>0.15</td><td>0.44</td><td>0.35</td><td>0.7288</td></tr><tr><td>(Phe-26) × (Phe-37)</td><td>0.07</td><td>0.22</td><td>0.33</td><td>0.7445</td></tr><tr><td>(Ile-4) × Phe-37</td><td>0.06</td><td>0.22</td><td>0.30</td><td>0.7675</td></tr></tbody></table>

the four mutated side chains (Fig. 5A). The magnitude of this interaction, -4.4 kcal/mol, is remarkably similar to the five-way coupling of about -6 kcal/mol associated with the catalytic activation of Mg²⁺, and to the interdomain coupling energies (~6 kcal/mol) estimated from the modular combinatorial analysis of contributions from CP1 and the ABD on amino acid activation, specificity, and tRNA<sup>Trp</sup> aminoacylation. Indeed, we propose that all effects represent the integrated, allosteric impact arising from domain motion of the ABD relative to the active site during catalysis of amino acid activation and acyl transfer.

Second, six of eight variants with wild-type Ile-4 contribute to higher specificity for tryptophan (*lightly shaded bars* in Fig. 5). Unlike the balanced contribution of the four D1 Switch residues to catalytic assist by Mg²⁺, the higher order interactions responsible for the allosteric enhancement of specificity suggest disproportionate participation by the wild-type Ile-4 residue. We explore implications of this observation under “Discussion.”

**The Role of the Metal**—Previously, we noted (4) that replacing Mg²⁺ with Mn²⁺ reduced the specificity for tryptophan, increasing the relative activation of tyrosine by ~3.5-fold. This value is almost precisely mirrored by the main effect of Mg²⁺ ion: -0.73 kcal/mol. Interactions involving Mg²⁺ actually lead to a net reduction of specificity: +2.9 kcal/mol. It appears, then, that unfavorable interactions between D1 residues and the metal modulate the overall net effect of D1 residues on the specificity ratio, about -8.3 kcal/mol, which implies a much higher tyrosine rejection than is observed.

DISCUSSION

Multidimensional thermodynamic cycles (13) are the “gold standard” for quantitative estimates of linkage-defining cou-

pling energies in allosteric and cooperativity. However, only rarely has the approach been extended to cycles of dimension >2, either in protein science or enzymology. Our results (5) and those of Sadovsky and Yifrach (20) illustrate two examples of functionally significant, higher order coupling interactions between individual residues. Our new paradigm, designed combinatorial mutagenesis of key dynamic residues connecting rigid bodies identified from structural studies, affords a path to new insight in other systems.

**Consistency between Modular, Point Mutant Thermodynamic Cycles**—We showed previously that catalytic assist by Mg²⁺ during tryptophan activation was the exclusive result of long range energetic coupling between the metal and the four D1 Switch residues mutated in this study. Implicating specific amino acids contributing to enhanced specificity in contemporary TrpRS, relative to that of the Urzyme uniquely complements that published recently (8) about the effects of adding the CP1 insertion or the ABD individually to the Urzyme. The surprising conclusion was that all three of the essential TrpRS functions, tryptophan activation, selection of tryptophan versus tyrosine, and tRNA<sup>Trp</sup> aminoacylation, benefit from approximately the same coupling energy, -6 kcal/mol between the two modules. This work presents remarkably consistent evidence that the functional intermodular communication is mediated by the D1 Switch.

**A Possible Role for Ile-4**—Differences between the higher order interactions illustrated in Fig. 5 and those illustrated in Fig. 4A of Ref. 5 merit comment. The highest order interaction is dominant by a significant margin in both distributions of coupling energies. In the allosteric coupling to catalytic assist by Mg²⁺ ion, none of the lower order interactions consistently involve the obvious, unidirectional affect of the same residue. In the present distribution, however, the Ile-4 residue is a prominent feature of the seven next largest effects (Fig. 5A), suggesting that this residue might serve an identifiable role in coordinating the effects of the D1 Switch in the selection of the correct amino acid. Ile-4 lies between the D1 Switch and the N-terminal β-strand that forms a “bottom” platform of the tryptophan-binding pocket, for which the specificity helix, residues 125–136, serves as the top. Thus, it occupies a unique position, relative to the tryptophan binding pocket. Fig. 2A illustrates that the top and bottom can be mechanically “pinched” together by the CP1 domain. Enhanced specificity arises primarily from differences in $k_{cat}$ (Fig. 3, B and C), suggesting that this effect, and the ensuing decision of whether or not to activate the bound amino acid, occur close to the transition state.

**The Observed Effects Are Associated with Changes in the Volume of the Tryptophan Binding Pocket**—TrpRS catalysis proceeds by a sequence of structural states represented by three distinct crystal structures (Protein Data Bank codes 1MAW, 1MAU, and 1I6L (7, 17)). We used the recently described algorithm POVME (21) to assess the tryptophan binding pocket volumes (287 Å³, 102 Å³, and 92 Å³) of these three states. These substantial changes suggest a possible mechanism to account for the association of specific substrate selection with domain movement, as well as the unique role of Ile-4: control over the tryptophan binding pocket volume by domain movement

<figure><img src="image_13.png" /></figure>

<footnote><sup>a</sup> Multidimensional thermodynamic analysis of TrpRS specificity.</footnote>
Allosteric Enhancement of Amino Acid Specificity in TrpRS

<figure><img src="image_14.png" /><figcaption>FIGURE 5. **Intrinsic and interaction effects responsible for specific selection of tryptophan versus tyrosine in B. stearothermophilus TrpRS.** A, intrinsic and interaction effects not involving the metal. B, values for terms involving coupling to Mg²⁺. Light bars denote terms that include wild-type residue Ile-4.</figcaption></figure>

might localize a stringent test for interactions specific for tryptophan to a point close to the chemical transition state.

Chemical transition states are so transient that their only experimental traces appear in kinetic measurements, such as those illustrated in Figs. 3 and 5. Computational methods, however, can potentially access the associated structural data. To that end, we have developed methods to construct the most probable pathways with minimum action (10, 17).⁴ These structural sequences illustrate semi-quantitatively (e.g. Fig. 6 and supplemental Movie S1) the successive states along the path (22), together with their energies.

Amino acid binding pocket volumes calculated for the trajectory using POVME (21) suggest two states along the path where these volumes assume minimum values (Fig. 6A), representing opportunities to validate the correct choice. Each follows closely after the corresponding conformational transition state encountered during induced-fit and catalysis (7).

Correlated movements of CP1 and the ABD impose subtle changes in the positions of both the N-terminal β-strand and the specificity helix, as suggested in Fig. 2A. These movements compress the volume by ~12% during the catalytic step (Fig. 6). The four D1 Switch residues rearrange in the conformational transition states that precede compression of pocket volumes (Fig. 6A and supplemental Movie S1), and so could couple specificity to these movements.

**Allosteric Enhancement Is Widely Observed in aaRS**—Both kinetic measurements (23) and indirect experiments (24) led to the early suggestion that amino acid specificity by aaRS involved the dominant role of *k*<sub>cat</sub>. Steady-state Michaelis-Menten kinetics of GlnRS mutants showed that its specific aminoacylation of tRNA<sup>Gln</sup> was primarily reflected in lower mutant *k*<sub>cat</sub> values (25, 26). Pre-steady state observations were made more recently for selectivity of GlnRS for glutamine *versus* glutamic acid (27–29). These results, and our own earlier failure to

improve specificity of TrpRS for tyrosine (9) can now be interpreted more clearly in terms of allosteric enhancement by domain motion. Similarly, the amino acid binding pocket size was implicated in specific recognition of phenylalanine by PheRS (30). The generality of these observations implies that allosteric enhancement of amino acid specificity by mechanisms analogous to those described here are a general property of all or most aaRS.

**Selective Advantage of Using Allostery to Enforce Specificity**—
The evident difficulty of changing aminoacyl-tRNA synthetase specificities by multiple obvious point mutations within the active site (9, 27–29) contrasts with the ease of creating orthogonal aaRS/tRNA pairs (31, 32). The apparent conflict merits two comments. First, the aaRSs derived for use in synthetic biology generally require sophisticated selection techniques involving suppression of nonsense mutations in essential genes by multisite mutants prepared from a “naive” amino acid binding pocket created by first mutating all relevant side chains in the amino acid binding pocket to alanine (33). Such a procedure may bypass the barrier to simplified mutagenesis, by selecting combinations that accommodate non-canonical amino acids without interfering overly much with the induced-fit and catalytic conformational changes. Second, whereas the engineered aaRS clearly accomplish their intended function of acylating orthogonal suppressor tRNAs with non-canonical amino acids, little has been done to define their actual specificity constants.

We are left with an important question: why off-load specific recognition to an allosteric mechanism? Several authors have suggested possible answers to this question. Cooperative transitions in proteins appear to be robust to mutation (15, 34). Thus, as the aminoacyl-tRNA synthetases depend strongly on retaining their high substrate specificities, it makes sense that allosteric mechanisms would be selected over time by the removal of point mutants that altered such specificity, whereas mutations leading to allos-

<figure><img src="image_15.png" /></figure>
Allosteric Enhancement of Amino Acid Specificity in TrpRS

<figure><img src="image_16.png" /></figure>

FIGURE 6. Trajectory parameters for modeling TrpRS interactions with substrate amino acids. A, a 200-step Minimum Action trajectory was computed to pass through three states provided by crystal structures (Protein Data Bank codes 1MAW, 1MAU, and 1I6L) that represent Open, Pre-transition, and Products states that demarcate the induced-fit and catalytic transitions, as indicated. Energies of states computed along the path are indicated by gray dots; tryptophan binding pocket volumes are given by black dots connected by the black line. Structures representing the first state with the new D1 side chain packing are indicated by large, filled circles, and coincide with the conformational transition states identified by PATH. Pocket volumes increase during both conformational transition states, and assume minimum values just thereafter, as indicated by the bold arrows. B, superposition of the tryptophan binding pocket residues from the minimum volumes following the conformational transition states associated with induced-fit (transparent) and catalysis as solid. The pocket (smudge green) is 12% smaller after the second transition state. The first β-strand and specificity helix are colored yellow and orange as described in the legend to Fig. 2.

teric mechanisms accumulated and began to function. Indeed, combinatorial mutagenesis of the D1 Switch residues in the TrpRS Urzyme will be a sensible extension of this work.

Acknowledgment—We thank H. Fried for helpful feedback and critical comments on the manuscript.

REFERENCES

1. Perutz, M. F., Wilkinson, A. J., Paoli, M., and Dodson, G. G. (1998) The stereochemical mechanism of the cooperative effects in hemoglobin revisited. *Annu. Rev. Biophys. Biomol. Struct.* **27**, 1–34

2. Perutz, M. F. (1970) Stereochemistry of cooperative effects of hemoglobin. *Nature* **228**, 726–739

3. Weinreb, V., and Carter, C. W., Jr. (2008) Mg²⁺-free *B. stearothermophil*
Allosteric Enhancement of Amino Acid Specificity in TrpRS

*lus* tryptophanyl-tRNA synthetase activates tryptophan with a major frac-
tion of the overall rate enhancement. *J. Am. Chem. Soc.* **130**, 1488–1494

4. Weinreb, V., Li, L., Campbell, C. L., Kaguni, L. S., and Carter, C. W., Jr. (2009) Mg²⁺-assisted catalysis by *B. stearothermophilus* TrpRS is promoted by allosteric effects. *Structure* **17**, 952–964

5. Weinreb, V., Li, L., and Carter, C. W., Jr. (2012) A master switch couples Mg²⁺-assisted catalysis to domain motion in *B. stearothermophilus* tryptophanyl-tRNA synthetase. *Structure* **20**, 128–138

6. Cammer, S., and Carter, C. W., Jr. (2010) Six rossmannoid folds, including the class I aminoacyl-tRNA synthetases, share a partial core with the anticodon-binding domain of a class II aminoacyl-tRNA synthetase. *Bioinformatics* **26**, 709–714

7. Kapustina, M., Weinreb, V., Li, L., Kuhlman, B., and Carter, C. W., Jr. (2007) A conformational transition state accompanies tryptophan activation by *B. stearothermophilus* tryptophanyl-tRNA synthetase. *Structure* **15**, 1272–1284

8. Li, L., and Carter, C. W., Jr. (2013) Full implementation of the genetic code by tryptophanyl-tRNA synthetase requires cooperativity between recently acquired structural modules. *J. Biol. Chem.* **288**, 34736–34745

9. Praetorius-Ibba, M., Stange-Thomann, N., Kitabatake, M., Ali, K., Söll, I., Carter, C. W., Jr., Ibba, M., and Söll, D. (2000) Ancient adaptation of the active site of tryptophanyl-tRNA synthetase for tryptophan binding. *Biochemistry* **39**, 13136–13143

10. Franklin, J., Koehl, P., Doniach, S., and Delarue, M. (2007) MinActionPath. Maximum likelihood trajectory for large-scale structural transitions in a coarse-grained locally harmonic energy landscape. *Nucleic Acids Res.* **35**, W477-W482

11. Francklyn, C. S., First, E. A., Perona, J. J., and Hou, Y.-M. (2008) Methods for kinetic and thermodynamic analysis of aminoacyl-tRNA synthetases. *Methods* **44**, 100–118

12. Retailleau, P., Hu, M., Vachette, P., and Carter, C. W., Jr. (2003) Interconversion of ATP binding and conformational free energies by tryptophanyl-tRNA synthetase. Structures of ATP bound to open and closed, pre-transition conformations. *J. Mol. Biol.* **325**, 39–63

13. Horovitz, A., and Fersht, A. R. (1990) Strategy for analyzing the co-operativity of intramolecular interactions in peptides and proteins. *J. Mol. Biol.* **214**, 613–617

14. Carter, C. W., Jr. (1990) Efficient factorial designs and the analysis of macromolecular crystal growth conditions. *Methods* **1**, 12–24

15. Weinkam, P., Pons, J., and Sali, A. (2012) Structure-based model of allosery predicts coupling between distant sites. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 4875–4880

16. Deleted in proof

17. Laowanapiban, P., Kapustina, M., Vonrhein, C., Delarue, M., Koehl, P., and Carter, C. W., Jr. (2009) Independent saturation of three TrpRS subsites generates a partially assembled state similar to those observed in molecular simulations. *Proc Natl. Acad. Sci. U.S.A.* **106**, 1790–1795

18. Pham, Y., Kuhlman, B., Butterfoss, G. L., Hu, H., Weinreb, V., and Carter, C. W., Jr. (2010) Tryptophanyl-tRNA synthetase Urzyme. A model to recapitulate molecular evolution and investigate intramolecular complementation. *J. Biol. Chem.* **285**, 38590–38601

19. Pham, Y., Li, L., Kim, A., Erdogan, O., Weinreb, V., Butterfoss, G. L., Kuhlman, B., and Carter, C. W., Jr. (2007) A minimal TrpRS catalytic domain supports sense/antisense ancestry of class I and II aminoacyl-tRNA synthetases. *Mol. Cell* **25**, 851–862

20. Sadovsky, E., and Yifrach, O. (2007) Principles underlying energetic coupling along an allosteric communication trajectory of a voltage-activated K channel. *Proc. Natl. Acad. Sci. U.S.A.* **104**, 19813–19818

21. Durrant, J. D., de Oliveira, C. A., and McCammon, J. A. (2011) POVME. An algorithm for measuring binding-pocket volumes. *J. Mol. Graph. Model* **29**, 773–776

22. Prevost, M. S., Sauguet, L., Nury, H., Van Renterghem, C., Huon, C., Poitevin, F., Baaden, M., Delarue, M., and Corringer, P.-J. (2012) A locally closed conformation of a bacterial pentameric proton-gated ion channel. *Nat. Struct. Mol. Biol.* **19**, 642–649

23. Lam, S. S., and Schimmel, P. R. (1975) Equilibrium measurements of cognate and noncognate interactions between aminoacyl transfer RNA synthetases and transfer RNA. *Biochemistry* **14**, 2775–2780

24. Ebel, J. P., Giegé, R., Bonnet, J., Kern, D., Befort, N., Bollack, C., Fasiolo, F., Gangloff, J., and Dirheimer, G. (1973) Factors determining the specificity of the tRNA aminoacylation reaction non-absolute specificity of tRNA-aminoacyl-tRNA synthetase recognition and particular importance of the maximal velocity. *Biochimie* **55**, 547–557

25. Rogers, M. J., Adachi, T., Inokuchi, H., and Söll, D. (1992) Switching tRNA<sup>Gln</sup> identity from glutamine to tryptophan. *Proc. Natl. Acad. Sci. U.S.A.* **89**, 3463–3467

26. Rogers, M. J., Weygand-Durasevi, I., Schwob, E., Sherman, J. M., Rogers, K. C., Adachi, T., Inokuchi, H., and Söll, D. (1993) Selectivity and specificity in the recognition of tRNA by *E. coli* by glutaminyl-tRNA synthetase. *Biochimie* **75**, 1083–1090

27. Bullock, T. L., Uter, N., Nissan, T. A., and Perona, J. J. (2003) Amino acid discrimination by a class I aminoacyl-tRNA synthetase specified by negative determinants. *J. Mol. Biol.* **328**, 395–408

28. Bullock, T. L., Rodríguez-Hernández, A., Corigliano, E. M., and Perona, J. J. (2008) A rationally engineered misacylating aminoacyl-tRNA synthetase. *Proc. Natl. Acad. Sci. U.S.A.* **105**, 7428–7433

29. Corigliano, E. M., and Perona, J. J. (2009) Architectural underpinnings of the genetic code for glutamine. *Biochemistry* **28**, 676–687

30. Ibba, M., Kast, P., and Hennecke, H. (1994) Substrate specificity is determined by amino acid binding pocket size in *Escherichia coli* phenylalanyl-tRNA synthetase. *Biochemistry* **33**, 7107–7112

31. Chin, J. W., Cropp, T. A., Anderson, J. C., Mukherji, M., Zhang, Z., and Schultz, P. G. (2003) An expanded eukaryotic genetic code. *Science* **301**, 964–967

32. Liu, C. C., and Schultz, P. G. (2010) Adding new chemistries to the genetic code. *Annu. Rev. Biochem.* **79**, 413–444

33. Wang, L., Brock, A., Herberich, B., and Schultz, P. G. (2001) Expanding the genetic code of *Escherichia coli*. *Science* **292**, 498–500

34. del Sol, A., Fujihashi, H., Amoros, D., and Nussinov, R. (2006) Residues crucial for maintaining short paths in network communication mediate signaling in proteins. *Mol. Syst. Biol.* **2**, 2006.0019

<figure><img src="image_17.png" /></figure>