Probing electrostatic interactions and ligand binding in aspartyl-tRNA synthetase through site-directed mutagenesis and computer simulations

Uncovering substrate specificity determinants of class IIb aminoacyl-tRNA synthetases with machine learning

Specific amino acid (AA) binding by aminoacyl-tRNA synthetases (aaRSs) is necessary for correct translation of the genetic code. Sequence and structure analyses have revealed the main specificity determinants and allowed a partitioning of aaRSs into two classes and several subclasses. However, the information contributed by each determinant has not been precisely quantified, and other, minor determinants may still be unidentified. Growth of genomic data and development of machine learning classification methods allow us to revisit these questions. This work considered the subclass IIb, formed by the three enzymes aspartyl-, asparaginyl-, and lysyl-tRNA synthetase (LysRS). Over 35,000 sequences from the Pfam database were considered, and used to train a machine-learning model based on ensembles of decision trees. The model was trained to reproduce the existing classification of each sequence as AspRS, AsnRS, or LysRS, and to identify which sequence positions were most important for the classification. A few positions (5-8 depending on the AA substrate) sufficed for accurate classification. Most but not all of them were well-known specificity determinants. The machine learning models thus identified sets of mutations that distinguish the three subclass members, which might be targeted in engineering efforts to alter or swap the AA specificities for biotechnology applications.

Site-Specifically Studying Lysine Acetylation of Aminoacyl-tRNA Synthetases

Aminoacyl-tRNA synthetases (AARSs) charge their cognate tRNAs with corresponding amino acids, playing key roles in ribosomal protein synthesis. A series of proteomic studies have demonstrated that AARSs have levels of lysine acetylation much higher than those of other proteins in Escherichia coli. To study AARS acetylation, 25 site-specifically acetylated variants of four AARSs were generated by the genetic code expansion strategy. Kinetic analyses were performed to biochemically characterize the impact of site-specific acetylation on AARS functions, including amino acid activation, tRNA aminoacylation, and editing activities. The results showed that impacts of acetylation were different between class I and class II AARSs and also varied among the same class of AARSs. The results also showed that acetylation of threonyl-tRNA synthetase (ThrRS) could affect its editing function. Both in vivo and in vitro studies were further performed to explore the acetylation and deacetylation processes of ThrRS. Although nonenzymatic acetylation and CobB-dependent deacetylation were concluded, the results also indicated the existence of additional modifying enzymes or mechanisms for ThrRS acetylation and deacetylation.

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.

Redesigning the stereospecificity of tyrosyl-tRNA synthetase

D-Amino acids are largely excluded from protein synthesis, yet they are of great interest in biotechnology. Unnatural amino acids have been introduced into proteins using engineered aminoacyl-tRNA synthetases (aaRSs), and this strategy might be applicable to D-amino acids. Several aaRSs can aminoacylate their tRNA with a D-amino acid; of these, tyrosyl-tRNA synthetase (TyrRS) has the weakest stereospecificity. We use computational protein design to suggest active site mutations in Escherichia coli TyrRS that could increase its D-Tyr binding further, relative to L-Tyr. The mutations selected all modify one or more sidechain charges in the Tyr binding pocket. We test their effect by probing the aminoacyl-adenylation reaction through pyrophosphate exchange experiments. We also perform extensive alchemical free energy simulations to obtain L-Tyr/D-Tyr binding free energy differences. Agreement with experiment is good, validating the structural models and detailed thermodynamic predictions the simulations provide. The TyrRS stereospecificity proves hard to engineer through charge-altering mutations in the first and second coordination shells of the Tyr ammonium group. Of six mutants tested, two are active towards D-Tyr; one of these has an inverted stereospecificity, with a large preference for D-Tyr. However, its activity is low. Evidently, the TyrRS stereospecificity is robust towards charge rearrangements near the ligand. Future design may have to consider more distant and/or electrically neutral target mutations, and possibly design for binding of the transition state, whose structure however can only be modeled.

Evaluation of residue-residue contact prediction in CASP10

We present the results of the assessment of the intramolecular residue-residue contact predictions from 26 prediction groups participating in the 10th round of the CASP experiment. The most recently developed direct coupling analysis methods did not take part in the experiment likely because they require a very deep sequence alignment not available for any of the 114 CASP10 targets. The performance of contact prediction methods was evaluated with the measures used in previous CASPs (i.e., prediction accuracy and the difference between the distribution of the predicted contacts and that of all pairs of residues in the target protein), as well as new measures, such as the Matthews correlation coefficient, the area under the precision-recall curve and the ranks of the first correctly and incorrectly predicted contact. We also evaluated the ability to detect interdomain contacts and tested whether the difficulty of predicting contacts depends upon the protein length and the depth of the family sequence alignment. The analyses were carried out on the target domains for which structural homologs did not exist or were difficult to identify. The evaluation was performed for all types of contacts (short, medium, and long-range), with emphasis placed on long-range contacts, i.e. those involving residues separated by at least 24 residues along the sequence. The assessment suggests that the best CASP10 contact prediction methods perform at approximately the same level, and comparably to those participating in CASP9.

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.

Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.

Species specificity of the complement inhibitor compstatin investigated by all-atom molecular dynamics simulations

The development of compounds to regulate the activation of the complement system in non-primate species is of profound interest because it can provide models for human diseases. The peptide compstatin inhibits protein C3 in primate mammals and is a potential therapeutic agent against unregulated activation of complement in humans but is inactive against nonprimate species. Here, we elucidate this species specificity of compstatin by molecular dynamics simulations of complexes between the most potent natural compstatin analog and human or rat C3. The results are compared against an experimental conformation of the human complex, determined recently by X-ray diffraction at 2.4-A resolution. The human complex simulations provide information on the relative contributions to stability of specific C3 and compstatin residues. In the rat simulations, the protein undergoes reproducible conformational changes, which eliminate or weaken specific interactions and reduce the complex stability. The simulation insights can be used to design improved compstatin-based inhibitors for human C3 and active inhibitors against lower mammals.

Computational design of protein-ligand binding: modifying the specificity of asparaginyl-tRNA synthetase

A method for computational design of protein-ligand interactions is implemented and tested on the asparaginyl- and aspartyl-tRNA synthetase enzymes (AsnRS, AspRS). The substrate specificity of these enzymes is crucial for the accurate translation of the genetic code. The method relies on a molecular mechanics energy function and a simple, continuum electrostatic, implicit solvent model. As test calculations, we first compute AspRS-substrate binding free energy changes due to nine point mutations, for which experimental data are available; we also perform large-scale redesign of the entire active site of each enzyme (40 amino acids) and compare to experimental sequences. We then apply the method to engineer an increased binding of aspartyl-adenylate (AspAMP) into AsnRS. Mutants are obtained using several directed evolution protocols, where four or five amino acid positions in the active site are randomized. Promising mutants are subjected to molecular dynamics simulations; Poisson-Boltzmann calculations provide an estimate of the corresponding, AspAMP, binding free energy changes, relative to the native AsnRS. Several of the mutants are predicted to have an inverted binding specificity, preferring to bind AspAMP rather than the natural substrate, AsnAMP. The computed binding affinities are significantly weaker than the native, AsnRS:AsnAMP affinity, and in most cases, the active site structure is significantly changed, compared to the native complex. This almost certainly precludes catalytic activity. One of the designed sequences has a higher affinity and more native-like structure and may represent a valid candidate for Asp activity.

Alchemical free energy simulations for biological complexes: powerful but temperamental

Free energy simulations compare multiple ligand:receptor complexes by "alchemically" transforming one into another, yielding binding free energy differences. Since their introduction in the 1980s, many technical and theoretical obstacles were surmounted, and the method ("MDFE," since molecular dynamics are often used) has matured into a powerful tool. We describe its current status, its effectiveness, and the challenges it faces. MDFE has provided chemical accuracy for many systems but remains expensive, with significant human overhead costs. The bottlenecks have shifted, partly due to increased computer power. To study diverse sets of ligands, force field availability and accuracy can be a major difficulty. Another difficulty is the frequent need to consider multiple states, related to sidechain protonation or buried waters, for example. Sophisticated, automated methods to sample these states are maturing, such as constant pH simulations. Meanwhile, combinations of MDFE and simpler approaches, like continuum dielectric models, can be very effective. As illustrations, we show how, with careful force field parameterization, MDFE accurately predicts binding specificities between complex tetracycline ligands and their targets. We describe substrate binding to the aspartyl-tRNA synthetase enzyme, where many distinct electrostatic states play a role, and a histidine and a Mg(2+) ion act as coupled switches that help enforce a strict preference for the aspartate substrate, relative to several analogs. Overall, MDFE has achieved a predictive status, where novel ligands can be studied and molecular recognition elucidated in depth. It should play an increasing role in the analysis of complex cellular processes and biomolecular engineering.

Exit strategies for charged tRNA from GluRS

For several class I aminoacyl-tRNA synthetases (aaRSs), the rate-determining step in aminoacylation is the dissociation of charged tRNA from the enzyme. In this study, the following factors affecting the release of the charged tRNA from aaRSs are computationally explored: the protonation states of amino acids and substrates present in the active site, and the presence and the absence of AMP and elongation factor Tu. Through molecular modeling, internal pK(a) calculations, and molecular dynamics simulations, distinct, mechanistically relevant post-transfer states with charged tRNA bound to glutamyl-tRNA synthetase from Thermus thermophilus (Glu-tRNA(Glu)) are considered. The behavior of these nonequilibrium states is characterized as a function of time using dynamical network analysis, local energetics, and changes in free energies to estimate transitions that occur during the release of the tRNA. The hundreds of nanoseconds of simulation time reveal system characteristics that are consistent with recent experimental studies. Energetic and network results support the previously proposed mechanism in which the transfer of amino acid to tRNA is accompanied by the protonation of AMP to H-AMP. Subsequent migration of proton to water reduces the stability of the complex and loosens the interface both in the presence and in the absence of AMP. The subsequent undocking of AMP or tRNA then proceeds along thermodynamically competitive pathways. Release of the tRNA acceptor stem is further accelerated by the deprotonation of the alpha-ammonium group on the charging amino acid. The proposed general base is Glu41, a residue binding the alpha-ammonium group that is conserved in both structure and sequence across nearly all class I aaRSs. This universal handle is predicted through pK(a) calculations to be part of a proton relay system for destabilizing the bound charging amino acid following aminoacylation. Addition of elongation factor Tu to the aaRS.tRNA complex stimulates the dissociation of the tRNA core and the tRNA acceptor stem.

Computer simulations reveal a novel nucleotide-type binding orientation for ellipticine-based anticancer c-kit kinase inhibitors

Receptor tyrosine kinase (RTK) enzymes regulate cell signaling pathways and so are an important target for cancer chemotherapy. Current inhibitors of c-kit, a key RTK stem cell factor receptor, are inactive against the most common mutated variant Asp816Val, associated with highly malignant cancers. Recent combined experimental/simulation work has highlighted the utility of the ellipticine pharmacore in inhibiting mutant c-kit, and the present simulation study applies a combination of high-level simulation tools to probe further the binding of ellipticine-based derivatives to c-kit. We find a large preference for protonation of bound ellipticine, which stabilizes the negative protein residues that coordinated ADP.Mg (2+) in the native complex. The resulting ellipticine inhibitor binding mode resembles the native nucleotide complex and serves to explain some existing experimental data on binding specificities, indicating that functionalization at the C4/C5 sites of ellipticine derivatives may be important for the design of novel nucleotide analogues that inhibit mutant c-kit.