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

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.

Distal Mutations in the β-Clamp of DNA Polymerase III* Disrupt DNA Orientation and Affect Exonuclease Activity

DNA polymerases are responsible for the replication and repair of DNA found in all DNA-based organisms. DNA Polymerase III is the main replicative polymerase of E. coli and is composed of over 10 proteins. A subset of these proteins (Pol III*) includes the polymerase (α), exonuclease (ϵ), clamp (β), and accessory protein (θ). Mutations of residues in, or around the active site of the catalytic subunits (α and ϵ), can have a significant impact on catalysis. However, the effects of distal mutations in noncatalytic subunits on the activity of catalytic subunits are less well-characterized. Here, we investigate the effects of two Pol III* variants, β-L82E/L82'E and β-L82D/L82'D, on the proofreading reaction catalyzed by ϵ. MD simulations reveal major changes in the dynamics of Pol III*, which extend throughout the complex. These changes are mostly induced by a shift in the position of the DNA substrate inside the β-clamp, although no major structural changes are observed in the protein complex. Quantum mechanics/molecular mechanics (QM/MM) calculations indicate that the β-L82D/L82'D variant has reduced catalytic proficiency due to highly endoergic reaction energies resulting from structural changes in the active site and differences in the electric field at the active site arising from the protein and substrate. Conversely, the β-L82E/L82'E variant is predicted to maintain proofreading activity, exhibiting a similar reaction barrier for nucleotide excision compared with the WT system. However, significant differences in the reaction mechanism are obtained due to the changes induced by the mutations on the β-clamp.

Reciprocally-Coupled Gating: Strange Loops in Bioenergetics, Genetics, and Catalysis

Bioenergetics, genetic coding, and catalysis are all difficult to imagine emerging without pre-existing historical context. That context is often posed as a "Chicken and Egg" problem; its resolution is concisely described by de Grasse Tyson: "The egg was laid by a bird that was not a chicken". The concision and generality of that answer furnish no details-only an appropriate framework from which to examine detailed paradigms that might illuminate paradoxes underlying these three life-defining biomolecular processes. We examine experimental aspects here of five examples that all conform to the same paradigm. In each example, a paradox is resolved by coupling "if, and only if" conditions for reciprocal transitions between levels, such that the consequent of the first test is the antecedent for the second. Each condition thus restricts fluxes through, or "gates" the other. Reciprocally-coupled gating, in which two gated processes constrain one another, is self-referential, hence maps onto the formal structure of "strange loops". That mapping uncovers two different kinds of forces that may help unite the axioms underlying three phenomena that distinguish biology from chemistry. As a physical analog for Gödel's logic, biomolecular strange-loops provide a natural metaphor around which to organize a large body of experimental data, linking biology to information, free energy, and the second law of thermodynamics.

Escapement mechanisms: Efficient free energy transduction by reciprocally-coupled gating

Conversion of the free energy of NTP hydrolysis efficiently into mechanical work and/or information by transducing enzymes sustains living systems far from equilibrium, and so has been of interest for many decades. Detailed molecular mechanisms, however, remain puzzling and incomplete. We previously reported that catalysis of tryptophan activation by tryptophanyl-tRNA synthetase, TrpRS, requires relative domain motion to re-position the catalytic Mg²⁺ ion, noting the analogy between that conditional hydrolysis of ATP and the escapement mechanism of a mechanical clock. The escapement allows the time-keeping mechanism to advance discretely, one gear at a time, if and only if the pendulum swings, thereby converting energy from the weight driving the pendulum into rotation of the hands. Coupling of catalysis to domain motion, however, mimics only half of the escapement mechanism, suggesting that domain motion may also be reciprocally coupled to catalysis, completing the escapement metaphor. Computational studies of the free energy surface restraining the domain motion later confirmed that reciprocal coupling: the catalytic domain motion is thermodynamically unfavorable unless the PPi product is released from the active site. These two conditional phenomena-demonstrated together only for the TrpRS mechanism-function as reciprocally-coupled gates. As we and others have noted, such an escapement mechanism is essential to the efficient transduction of NTP hydrolysis free energy into other useful forms of mechanical or chemical work and/or information. Some implementation of both gating mechanisms-catalysis by domain motion and domain motion by catalysis-will thus likely be found in many other systems.

Characterizing the relation of functional and Early Folding Residues in protein structures using the example of aminoacyl-tRNA synthetases

Proteins are chains of amino acids which adopt a three-dimensional structure and are then able to catalyze chemical reactions or propagate signals in organisms. Without external influence, many proteins fold into their native structure, and a small number of Early Folding Residues (EFR) have previously been shown to initiate the formation of secondary structure elements and guide their respective assembly. Using the two diverse superfamilies of aminoacyl-tRNA synthetases (aaRS), it is shown that the position of EFR is preserved over the course of evolution even when the corresponding sequence conservation is small. Folding initiation sites are positioned in the center of secondary structure elements, independent of aaRS class. In class I, the predicted position of EFR resembles an ancient structural packing motif present in many seemingly unrelated proteins. Furthermore, it is shown that EFR and functionally relevant residues in aaRS are almost entirely disjoint sets of residues. The Start2Fold database is used to investigate whether this separation of EFR and functional residues can be observed for other proteins. EFR are found to constitute crucial connectors of protein regions which are distant at sequence level. Especially, these residues exhibit a high number of non-covalent residue-residue contacts such as hydrogen bonds and hydrophobic interactions. This tendency also manifests as energetically stable local regions, as substantiated by a knowledge-based potential. Despite profound differences regarding how EFR and functional residues are embedded in protein structures, a strict separation of structurally and functionally relevant residues cannot be observed for a more general collection of proteins.

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.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNA synthetases (aaRSs) are modular enzymes globally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation. Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g., in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show huge structural plasticity related to function and limited idiosyncrasies that are kingdom or even species specific (e.g., the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS). Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably between distant groups such as Gram-positive and Gram-negative Bacteria. The review focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation, and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulated in last two decades is reviewed, showing how the field moved from essentially reductionist biology towards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRS paralogs (e.g., during cell wall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointed throughout the review and distinctive characteristics of bacterium-like synthetases from organelles are outlined.

A modified PATH algorithm rapidly generates transition states comparable to those found by other well established algorithms

PATH rapidly computes a path and a transition state between crystal structures by minimizing the Onsager-Machlup action. It requires input parameters whose range of values can generate different transition-state structures that cannot be uniquely compared with those generated by other methods. We outline modifications to estimate these input parameters to circumvent these difficulties and validate the PATH transition states by showing consistency between transition-states derived by different algorithms for unrelated protein systems. Although functional protein conformational change trajectories are to a degree stochastic, they nonetheless pass through a well-defined transition state whose detailed structural properties can rapidly be identified using PATH.

Urzymology: experimental access to a key transition in the appearance of enzymes

Urzymes are catalysts derived from invariant cores of protein superfamilies. Urzymes from both aminoacyl-tRNA synthetase classes possess sophisticated catalytic mechanisms: pre-steady state bursts, significant transition-state stabilization of both amino acid activation, and tRNA acylation. However, they have insufficient specificity to ensure a fully developed genetic code, suggesting that they participated in synthesizing statistical proteins. They represent a robust experimental platform from which to articulate and test hypotheses both about their own ancestors and about how they, in turn, evolved into modern enzymes. They help reshape numerous paradigms from the RNA World hypothesis to protein structure databases and allostery.

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.

Full implementation of the genetic code by tryptophanyl-tRNA synthetase requires intermodular coupling

Tryptophanyl-tRNA Synthetase (TrpRS) Urzyme (fragments A and C), a 130-residue construct containing only secondary structures positioning the HIGH and KMSKS active site signatures and the specificity helix, accelerates tRNA(Trp) aminoacylation with ∼10-fold specificity toward tryptophan, relative to structurally related tyrosine. We proposed that including the 76-residue connecting peptide 1 insertion (Fragment B) might enhance tryptophan affinity and hence amino acid specificity, because that subdomain constrains the orientation of the specificity helix. We test that hypothesis by characterizing two new constructs: the catalytic domain (fragments A-C) and the Urzyme supplemented with the anticodon-binding domain (fragments A, C, and D). The three constructs, together with the full-length enzyme (fragments A-D), comprise a factorial experiment from which we deduce individual and combined contributions of the two modules to the steady-state kinetics parameters for tryptophan-dependent (32)PPi exchange, specificity for tryptophan versus tyrosine, and aminoacylation of tRNA(Trp). Factorial design directly measures the energetic coupling between the two more recent modules in the contemporary enzyme and demonstrates its functionality. Combining the TrpRS Urzyme individually in cis with each module affords an analysis of long term evolution of amino acid specificity and tRNA aminoacylation, both essential for expanding the genetic code. Either module significantly enhances tryptophan activation but unexpectedly eliminates amino acid specificity for tryptophan, relative to tyrosine, and significantly reduces tRNA aminoacylation. Exclusive dependence of both enhanced functionalities of full-length TrpRS on interdomain coupling energies between the two new modules argues that independent recruitment of connecting peptide 1 and the anticodon-binding domain during evolutionary development of Urzymes would have entailed significant losses of fitness.

Evolutionary and mechanistic insights into substrate and product accommodation of CTP:phosphocholine cytidylyltransferase from Plasmodium falciparum

The enzyme CTP:phosphocholine cytidylyltransferase (CCT) is essential in the lipid biosynthesis of Plasmodia (Haemosporida), presenting a promising antimalarial target. Here, we identified two independent gene duplication events of CCT within Apicomplexa and characterized a truncated construct of Plasmodium falciparum CCT that forms a dimer resembling the molecular architecture of CCT enzymes from other sources. Based on biophysical and enzyme kinetics methods, our data show that the CDP-choline product of the CCT enzymatic reaction binds to the enzyme considerably stronger than either substrate (CTP or choline phosphate). Interestingly, in the presence of Mg²⁺ , considered to be a cofactor of the enzyme, the binding of the CTP substrate is attenuated by a factor of 5. The weaker binding of CTP:Mg²⁺ , similarly to the related enzyme family of aminoacyl tRNA synthetases, suggests that, with lack of Mg²⁺ , positively charged side chain(s) of CCT may contribute to CTP accommodation. Thermodynamic investigations by isothermal titration calorimetry and fluorescent spectroscopy studies indicate that accommodation of the choline phosphate moiety in the CCT active site is different when it appears on its own as one of the substrates or when it is linked to the CDP-choline product. A tryptophan residue within the active site is identified as a useful internal fluorescence sensor of enzyme-ligand binding. Results indicate that the catalytic mechanism of Plasmodium falciparum CCT may involve conformational changes affecting the choline subsite of the enzyme.

Aminoacyl-tRNA Synthetases in the Bacterial World

Aminoacyl-tRNAsynthetases (aaRSs) are modular enzymesglobally conserved in the three kingdoms of life. All catalyze the same two-step reaction, i.e., the attachment of a proteinogenic amino acid on their cognate tRNAs, thereby mediating the correct expression of the genetic code. In addition, some aaRSs acquired other functions beyond this key role in translation.Genomics and X-ray crystallography have revealed great structural diversity in aaRSs (e.g.,in oligomery and modularity, in ranking into two distinct groups each subdivided in 3 subgroups, by additional domains appended on the catalytic modules). AaRSs show hugestructural plasticity related to function andlimited idiosyncrasies that are kingdom or even speciesspecific (e.g.,the presence in many Bacteria of non discriminating aaRSs compensating for the absence of one or two specific aaRSs, notably AsnRS and/or GlnRS).Diversity, as well, occurs in the mechanisms of aaRS gene regulation that are not conserved in evolution, notably betweendistant groups such as Gram-positive and Gram-negative Bacteria.Thereview focuses on bacterial aaRSs (and their paralogs) and covers their structure, function, regulation,and evolution. Structure/function relationships are emphasized, notably the enzymology of tRNA aminoacylation and the editing mechanisms for correction of activation and charging errors. The huge amount of genomic and structural data that accumulatedin last two decades is reviewed,showing how thefield moved from essentially reductionist biologytowards more global and integrated approaches. Likewise, the alternative functions of aaRSs and those of aaRSparalogs (e.g., during cellwall biogenesis and other metabolic processes in or outside protein synthesis) are reviewed. Since aaRS phylogenies present promiscuous bacterial, archaeal, and eukaryal features, similarities and differences in the properties of aaRSs from the three kingdoms of life are pinpointedthroughout the reviewand distinctive characteristics of bacterium-like synthetases from organelles are outlined.

Structural diversity and protein engineering of the aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases (aaRS) are the enzymes that ensure faithful transmission of genetic information in all living cells, and are central to the developing technologies for expanding the capacity of the translation apparatus to incorporate nonstandard amino acids into proteins in vivo. The 24 known aaRS families are divided into two classes that exhibit functional evolutionary convergence. Each class features an active site domain with a common fold that binds ATP, the amino acid, and the 3'-terminus of tRNA, embellished by idiosyncratic further domains that bind distal portions of the tRNA and enhance specificity. Fidelity in the expression of the genetic code requires that the aaRS be selective for both amino acids and tRNAs, a substantial challenge given the presence of structurally very similar noncognate substrates of both types. Here we comprehensively review central themes concerning the architectures of the protein structures and the remarkable dual-substrate selectivities, with a view toward discerning the most important issues that still substantially limit our capacity for rational protein engineering. A suggested general approach to rational design is presented, which should yield insight into the identities of the protein-RNA motifs at the heart of the genetic code, while also offering a basis for improving the catalytic properties of engineered tRNA synthetases emerging from genetic selections.

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.

Search for identical octapeptides in unrelated proteins: Structural plasticity revisited

Since proteins are dynamic in nature, they can alter their local structure in response to changes in their environment factors such as temperature, pH, phosphorylation, and binding of other small molecules. These conformational changes are extremely important for the correct folding and functioning of proteins. There are also a number of diseases associated with protein conformational change such as amyloid diseases. To stimulate research into the above factors which specify one conformation over another, different theoretical models have been proposed and tested against sequence similar distant structure protein fragments. In order to simplify the computational complexity of identifying conformational changes in proteins, various local sequence search algorithms were employed and the structural plasticity in unrelated proteins was examined by various research groups. In the present work, we revisit the mechanism of structural plasticity in unrelated proteins with increased number of structures in Protein Data Bank by comparing identical octapeptides in unrelated proteins with dictionary of protein secondary structure extracted from existing experimental data. Our goal is to bring out the influence of hydrophobic residues, hydrophilic residues, flanking residues, difference in secondary structural propensities of surrounding residues, difference in phi-psi angles and local and nonlocal interactions in identical octapeptides adopting different conformations. Also we have used surrounding hydrophobicity, environment dependent interaction energy, atomic mean force potential, structural unit contacts and difference profiles models to explore the factors which cause structural plasticity. The results discussed here may provide insights into protein folding, design and function.

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.

Six Rossmannoid folds, including the Class I aminoacyl-tRNA synthetases, share a partial core with the anti-codon-binding domain of a Class II aminoacyl-tRNA synthetase

MOTIVATION

Similarities in core residue packing provide evidence for divergence or convergence not reported using other methods.

RESULTS

We apply a new method for rapid structure comparison based on Simplicial Neighborhood Analysis of Protein Packing (SNAPP) to the diverse structural classification of proteins (SCOP) alpha/beta-class of protein folds. The procedure identifies inter-residue packing motifs shared by protein pairs from different folds. A threshold of 0.67 A RMSD for all atoms of corresponding residues ensures inclusion of only highly significant similarities comparable with those observed for identical catalytic residues in homologues. Many tertiary packing motifs are shared among the three classical Rossmannoid folds, as well as thousands of other motifs that occur in at least two distinct folds. Merging of neighboring packing motifs facilitated recognition of larger, recurrent substructures or cores. The anti-codon-binding domain of an archeal aminoacyl-tRNA synthetase (aaRS) was discovered to possess a packed core in which eight identical amino acid residues are within 0.55 A RMSD of the comparable structure in the FixJ receiver, a member of the Rossmannoid family that also includes the CheY signaling protein and flavodoxin-like proteins. Further investigation identified close variants of this core in five other Rossmannoid folds, including a functionally relevant core in Class Ia aminoacyl-tRNA synthetases. Although it is possible that the two essentially identical cores in the ProRS anti-codon-binding domain and the FixJ receiver converged to the same structure, the consensus core obtained from the structural and sequence alignments suggests that all the implicated protein folds descended from a simpler ancestral protein in which this core provided nucleotide binding and proto-allosteric functions.

AVAILABILITY

Programs are available at http://staff.vbi.vt.edu/cammer/snapp/download/

IMPLEMENTATION

Programs were written in Perl and c and run under Linux.

CONTACT

cammer@vbi.vt.edu.

Crystal structures of Saccharomyces cerevisiae tryptophanyl-tRNA synthetase: new insights into the mechanism of tryptophan activation and implications for anti-fungal drug design

Specific activation of amino acids by aminoacyl-tRNA synthetases is essential for maintaining translational fidelity. Here, we present crystal structures of Saccharomyces cerevisiae tryptophanyl-tRNA synthetase (sTrpRS) in apo form and in complexes with various ligands. In each complex, there is a sulfate ion bound at the active site which mimics the α- or β-phosphate group of ATP during tryptophan activation. In particular, in one monomer of the sTrpRS–TrpNH₂O complex, the sulfate ion appears to capture a snapshot of the α-phosphate of ATP during its movement towards tryptophan. Simulation study of a human TrpRS–Trp–ATP model shows that during the catalytic process the α-phosphate of ATP is driven to an intermediate position equivalent to that of the sulfate ion, then moves further and eventually fluctuates at around 2 Å from the nucleophile. A conserved Arg may interact with the oxygen in the scissile bond at the transition state, indicating its critical role in the nucleophilic substitution. Taken together, eukaryotic TrpRSs may adopt an associative mechanism for tryptophan activation in contrast to a dissociative mechanism proposed for bacterial TrpRSs. In addition, structural analysis of the apo sTrpRS reveals a unique feature of fungal TrpRSs, which could be exploited in rational antifungal drug design.

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.

Crystal structure of a mammalian CTP: phosphocholine cytidylyltransferase catalytic domain reveals novel active site residues within a highly conserved nucleotidyltransferase fold

CTP:phosphocholine cytidylyltransferase (CCT) is the key regulatory enzyme in the synthesis of phosphatidylcholine, the most abundant phospholipid in eukaryotic cell membranes. The CCT-catalyzed transfer of a cytidylyl group from CTP to phosphocholine to form CDP-choline is regulated by a membrane lipid-dependent mechanism imparted by its C-terminal membrane binding domain. We present the first analysis of a crystal structure of a eukaryotic CCT. A deletion construct of rat CCTalpha spanning residues 1-236 (CCT236) lacks the regulatory domain and as a result displays constitutive activity. The 2.2-A structure reveals a CCT236 homodimer in complex with the reaction product, CDP-choline. Each chain is composed of a complete catalytic domain with an intimately associated N-terminal extension, which together with the catalytic domain contributes to the dimer interface. Although the CCT236 structure reveals elements involved in binding cytidine that are conserved with other members of the cytidylyltransferase superfamily, it also features nonconserved active site residues, His-168 and Tyr-173, that make key interactions with the beta-phosphate of CDP-choline. Mutagenesis and kinetic analyses confirmed their role in phosphocholine binding and catalysis. These results demonstrate structural and mechanistic differences in a broadly conserved protein fold across the cytidylyltransferase family. Comparison of the CCT236 structure with those of other nucleotidyltransferases provides evidence for substrate-induced active site loop movements and a disorder-to-order transition of a loop element in the catalytic mechanism.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

Reaction mechanism of the epsilon subunit of E. coli DNA polymerase III: insights into active site metal coordination and catalytically significant residues

The 28 kDa epsilon subunit of Escherichia coli DNA polymerase III is the exonucleotidic proofreader responsible for editing polymerase insertion errors. Here, we study the mechanism by which epsilon carries out the exonuclease activity. We performed quantum mechanics/molecular mechanics calculations on the N-terminal domain containing the exonuclease activity. Both the free-epsilon and a complex epsilon bound to a theta homologue (HOT) were studied. For the epsilon-HOT complex Mg(2+) or Mn(2+) were investigated as the essential divalent metal cofactors, while only Mg(2+) was used for free-epsilon. In all calculations a water molecule bound to the catalytic metal acts as the nucleophile for hydrolysis of the phosphate bond. Initially, a direct proton transfer to H162 is observed. Subsequently, the nucleophilic attack takes place followed by a second proton transfer to E14. Our results show that the reaction catalyzed with Mn(2+) is faster than that with Mg(2+), in agreement with experiment. In addition, the epsilon-HOT complex shows a slightly lower energy barrier compared to free-epsilon. In all cases the catalytic metal is observed to be pentacoordinated. Charge and frontier orbital analyses suggest that charge transfer may stabilize the pentacoordination. Energy decomposition analysis to study the contribution of each residue to catalysis suggests that there are several important residues. Among these, H98, D103, D129, and D146 have been implicated in catalysis by mutagenesis studies. Some of these residues were found to be structurally conserved on human TREX1, the exonuclease domains from E. coli DNA-Pol I, and the DNA polymerase of bacteriophage RB69.

Independent saturation of three TrpRS subsites generates a partially assembled state similar to those observed in molecular simulations

Two new crystal structures of Bacillus stearothermophilus tryptophanyl-tRNA synthetase (TrpRS) afford evidence that a closed interdomain hinge angle requires a covalent bond between AMP and an occupant of either pyrophosphate or tryptophan subsite. They also are within experimental error of a cluster of structures observed in a nonequilibrium molecular dynamics simulation showing partial active-site assembly. Further, the highest energy structure in a minimum action pathway computed by using elastic network models for Open and Pretransition state (PreTS) conformations for the fully liganded TrpRS monomer is intermediate between that simulated structure and a partially disassembled structure from a nonequilibrium molecular dynamics trajectory for the unliganded PreTS. These mutual consistencies provide unexpected validation of inferences drawn from molecular simulations.

HIV-1 protease function and structure studies with the simplicial neighborhood analysis of protein packing method

The Simplicial Neighborhood Analysis of Protein Packing (SNAPP) method was used to predict the effect of mutagenesis on the enzymatic activity of the HIV-1 protease (HIVP). SNAPP relies on a four-body statistical scoring function derived from the analysis of spatially nearest neighbor residue compositional preferences in a diverse and representative subset of protein structures from the Protein Data Bank. The method was applied to the analysis of HIVP mutants with residue substitutions in the hydrophobic core as well as at the interface between the two protease monomers. Both wild-type and tethered structures were employed in the calculations. We obtained a strong correlation, with R(2) as high as 0.96, between DeltaSNAPP score (i.e., the difference in SNAPP scores between wild-type and mutant proteins) and the protease catalytic activity for tethered structures. However, a weaker but significant correlation was obtained for nontethered structures. Our analysis identified residues both in the hydrophobic core and at the dimeric interface that are very important for the protease function. This study demonstrates a potential utility of the SNAPP method for rational design of mutagenesis studies and protein engineering.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.

Catalytic mechanism of the tryptophan activation reaction revealed by crystal structures of human tryptophanyl-tRNA synthetase in different enzymatic states

Human tryptophanyl-tRNA synthetase (hTrpRS) differs from its bacterial counterpart at several key positions of the catalytic active site and has an extra N-terminal domain, implying possibly a different catalytic mechanism. We report here the crystal structures of hTrpRS in complexes with Trp, tryptophanamide and ATP and tryptophanyl-AMP, respectively, which represent three different enzymatic states of the Trp activation reaction. Analyses of these structures reveal the molecular basis of the mechanisms of the substrate recognition and the activation reaction. The dimeric hTrpRS is structurally and functionally asymmetric with half-of-the-sites reactivity. Recognition of Trp is by an induced-fit mechanism involving conformational change of the AIDQ motif that creates a perfect pocket for the binding and activation of Trp and causes coupled movements of the N-terminal and C-terminal domains. The KMSAS loop appears to have an inherent flexibility and the binding of ATP stabilizes it in a closed conformation that secures the position of ATP for catalysis. Our structural data indicate that the catalytic mechanism of the Trp activation reaction by hTrpRS involves more moderate conformational changes of the structural elements at the active site to recognize and bind the substrates, which is more complex and fine-tuned than that of bacterial TrpRS.

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

Few experimental data are available for rates of enzymatic phosphoryl-transfer reactions in the absence of the divalent metal ions associated with such reactions. Such data are of interest for amino acid activation by class Ic aminoacyl-tRNA synthetases, for which there is substantial evidence that binding energy of ATP may account for a major fraction of the overall rate enhancement, and it is crucial to know if these effects themselves depend on the divalent metal ion. We describe a nested, nonlinear model for the sum of metal-free and metal-catalyzed activities and its use in determining metal-free enzyme activity jointly with transition-state metal binding affinity, by fitting observed values obtained from Mg2+-depleted assays with increasing [EDTA] at known [Mg2+]total. Tryptophan activation by Bacillus stearothermophilus tryptophanyl-tRNA synthetase falls asymptotically to a plateau value 5 orders of magnitude below that observed for the Mg2+-supplemented enzyme at EDTA concentrations that reduce the free metal concentration to <1 pmolar. The fitted regression model parameters yield a relative rate acceleration of 9.3 x 10(4) attributable to the catalytic effect of Mg2+ and an enhanced (K(E)(double dagger) = 1.15 x 10(-7) M) transition-state binding of Mg2+. Factorial analysis indicates that 80% of the reduction in free energy of activation effected by TrpRS arises from protein-ligand interactions.