Abstract:The pyrrolysyl-tRNA synthetase (PylRS)/tRNACUA pair is one of the most widely used systems for the incorporation of noncanonical amino acids (ncAAs) into proteins at specific positions. Although directed evolution of PylRS have enabled over 300 ncAAs to be incorporated into proteins, most of the ncAA-containing proteins are expressed in a limited yield due to low activities of PylRS variants. Here, we applied machine learning (ML) to engineer the tRNA-binding domain of PylRS with a fast Fourier transform-partial least square regression (FFT-PLSR) model and three zero-shot prediction ML models. FFT-PLSR was first applied to explore a sequence space composed of pairwise combinations of 12 single mutations, and the best variant, Com1-IFRS, showed an 11-fold increase in activity compared to IFRS, a PylRS variant. The deep learning models ESM-1v, Mutcompute, and ProRefiner were then used to identify new mutation sites impacting the activity of Com1-IFRS. FFT-PLSR was used again to identify a variant, Com2-IFRS, from a sequence space containing 11520 mutations, which showed a 30-fold increase in activity. Com2-IFRS also enhanced enzyme activity against 12 other ncAAs by up to 3944.8-fold. Transplantation of the evolved mutations into 7 other PylRS-derived synthetases improved yields of proteins containing six types of ncAAs, including derivatives of Phe, Tyr, Trp, Cys, His and Lys, by up to 1149.7-fold. Molecular dynamics simulations revealed that mutations reshaped the hydrogen bond network between tRNA and protein, which increased tRNA binding affinity, shortened the reaction distance between tRNA and ncAA, and even enhanced the dynamics correlation network. This paper offers new PylRS variants that increase the utility of the orthogonal translation system and provide a machine learning framework for identifying optimized multiple-point combinatorial mutations in a vast sequence space.

tRNA shape is an identity element for an archaeal pyrrolysyl-tRNA synthetase from the human gut

Ancestral archaea expanded the genetic code with pyrrolysine

tRNAPyl: Structure, function, and applications

Pyrrolo-C As a Molecular Probe for Monitoring Conformations of the Trna 3' End.

Coexistence of bacterial leucyl-tRNA synthetases with archaeal tRNA binding domains that distinguish tRNA(Leu) in the archaeal mode.

The tRNA discriminator base defines the mutual orthogonality of two distinct pyrrolysyl-tRNA synthetase/tRNAPyl pairs in the same organism

Chimeric design of pyrrolysyl-tRNA synthetase/tRNA pairs and canonical synthetase/tRNA pairs for genetic code expansion

Pyrrolysyl-tRNA synthetase: an ordinary enzyme but an outstanding genetic code expansion tool.

The Pyrrolysyl-tRNA Synthetase Activity Can Be Improved by a P188 Mutation That Stabilizes the Full-Length Enzyme.

Aminoacylation of tRNA 2'- or 3'-hydroxyl by phosphoseryl- and pyrrolysyl-tRNA synthetases

Methanomethylophilus alvus Mx1201 Provides Basis for Mutual Orthogonal Pyrrolysyl tRNA/Aminoacyl-tRNA Synthetase Pairs in Mammalian Cells

Directed Evolution of Candidatus Methanomethylophilus alvus Pyrrolysyl-tRNA Synthetase for the Genetic Incorporation of Two Different Noncanonical Amino Acids in One Protein

Structural basis for aminoacylation of cellular modified tRNA by human lysyl-tRNA synthetase

Structural basis for shape-selective recognition and aminoacylation of a D-armless human mitochondrial tRNA

X-shaped structure of bacterial heterotetrameric tRNA synthetase suggests cryptic prokaryote functions and a rationale for synthetase classifications

Distinct Ways of G:U Recognition by Conserved Trna Binding Motifs.

A Designed, Highly Efficient Pyrrolysyl-tRNA Synthetase Mutant Binds O-Chlorophenylalanine Using Two Halogen Bonds

ThrRS-Mediated Translation Regulation of the RNA Polymerase III Subunit RPC10 Occurs through an Element with Similarity to Cognate tRNA ASL and Affects tRNA Levels

Machine Learning-Guided Evolution of Pyrrolysyl-Trna Synthetase for Improved Incorporation Efficiency of Diverse Noncanonical Amino Acids

The Yin and Yang of Trna: Proper Binding of Acceptor End Determines the Catalytic Balance of Editing and Aminoacylation

Amino acid activation analysis of primitive aminoacyl-tRNA synthetases encoded by both strands of a single gene using the malachite green assay