Simulation study of the ability of a computationally‐designed peptide to recognize target tRNALys3 and other decoy tRNAs

In this paper, we investigate the ability of our computationally‐designed peptide, Pept10 (PNWNGNRWLNNCLRG), to recognize the anticodon stem and loop (ASL) domain of the hypermodified tRNALys3 (mcm5s2U34,ms2t6A37), a reverse transcription primer of HIV replication. Five other ASLs, the singly modified ASLLys3(ms2t6A37), ASLLys3(s2U34), ASLLys3(Ψ39), ASLLys1,2(t6A37), and ASLGlu(s2U34), were used as decoys. Explicit‐solvent atomistic molecular dynamics simulations were performed to examine the process of binding of Pept10 with the target ASLLys3(mcm5s2U34,ms2t6A37) and the decoy ASLs. Simulation results demonstrated that Pept10 is capable of recognizing the target ASLLys3(mcm5s2U34,ms2t6A37) as well as one of the decoys, ASLLys3(Ψ39), but screens out the other four decoy ASLs. The interchain van der Waals (VDW) and charge–charge (ELE + EGB) energies for the two best complexes were evaluated to shed light on the molecular recognition mechanism between Pept10 and ASLs. The results indicated that Pept10 recognizes and binds to the target ASLLys3 (mcm5s2U34,ms2t6A37) through residues W3 and R7 which interact with the nucleotides mcm5s2U34, U35, and ms2t6A37 via the interchain VDW energy. Pept10 also recognizes the decoy ASLLys3(Ψ39) through residue R14 which contacts the nucleotide U36 via the interchain VDW energy. Regardless of the type of ASL, the positively charged arginines on Pept10 are attracted to the negatively charged phosphate linkages on the ASL via the interchain ELE + EGB energy, thereby enhancing the binding affinity.

[1]  Joshua N. Leonard,et al.  Adding energy minimization strategy to peptide‐design algorithm enables better search for RNA‐binding peptides: Redesigned λ N peptide binds boxB RNA , 2016, J. Comput. Chem..

[2]  Xingqing Xiao,et al.  Introducing folding stability into the score function for computational design of RNA‐binding peptides boosts the probability of success , 2016, Proteins.

[3]  C. Hall,et al.  Designing peptide sequences in flexible chain conformations to bind RNA: a search algorithm combining Monte Carlo, self-consistent mean field and concerted rotation techniques. , 2015, Journal of chemical theory and computation.

[4]  C. Hall,et al.  Molecular recognition mechanism of peptide chain bound to the tRNALys3 anticodon loop in silico , 2015, Journal of biomolecular structure & dynamics.

[5]  C. Hall,et al.  The design of a peptide sequence to inhibit HIV replication: a search algorithm combining Monte Carlo and self-consistent mean field techniques , 2014, Journal of biomolecular structure & dynamics.

[6]  C. Hall,et al.  Amino Acid Signature Enables Proteins to Recognize Modified tRNA , 2014, Biochemistry.

[7]  J. Su,et al.  The interactions and recognition of cyclic peptide mimetics of Tat with HIV-1 TAR RNA: a molecular dynamics simulation study , 2013, Journal of biomolecular structure & dynamics.

[8]  Christian Kramer,et al.  MM/GBSA Binding Energy Prediction on the PDBbind Data Set: Successes, Failures, and Directions for Further Improvement , 2013, J. Chem. Inf. Model..

[9]  Franck A. P. Vendeix,et al.  Human tRNA(Lys3)(UUU) is pre-structured by natural modifications for cognate and wobble codon binding through keto-enol tautomerism. , 2012, Journal of molecular biology.

[10]  Julian Tirado-Rives,et al.  Improving MM-GB/SA Scoring through the Application of the Variable Dielectric Model. , 2011, Journal of chemical theory and computation.

[11]  William D. Graham,et al.  Functional recognition of the modified human tRNALys3(UUU) anticodon domain by HIV's nucleocapsid protein and a peptide mimic. , 2011, Journal of molecular biology.

[12]  Yunjie Zhao,et al.  Improvements of the Hierarchical Approach for Predicting RNA Tertiary Structure , 2011, Journal of biomolecular structure & dynamics.

[13]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[14]  Celia A Schiffer,et al.  Decomposing the energetic impact of drug resistant mutations in HIV-1 protease on binding DRV. , 2010, Journal of chemical theory and computation.

[15]  Shigeyuki Yokoyama,et al.  Tertiary structure checkpoint at anticodon loop modification in tRNA functional maturation , 2009, Nature Structural &Molecular Biology.

[16]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[17]  Z. Luthey-Schulten,et al.  Dynamical networks in tRNA:protein complexes , 2009, Proceedings of the National Academy of Sciences.

[18]  Paul F Agris,et al.  Bringing order to translation: the contributions of transfer RNA anticodon‐domain modifications , 2008, EMBO reports.

[19]  Kai Zhu,et al.  Improved Methods for Side Chain and Loop Predictions via the Protein Local Optimization Program:  Variable Dielectric Model for Implicitly Improving the Treatment of Polarization Effects. , 2007, Journal of chemical theory and computation.

[20]  Jianyin Shao,et al.  Clustering Molecular Dynamics Trajectories: 1. Characterizing the Performance of Different Clustering Algorithms. , 2007, Journal of chemical theory and computation.

[21]  S. Vishveshwara,et al.  A study of communication pathways in methionyl- tRNA synthetase by molecular dynamics simulations and structure network analysis , 2007, Proceedings of the National Academy of Sciences.

[22]  Jacquelyn S Fetrow,et al.  Using molecular dynamics to map interaction networks in an aminoacyl‐tRNA synthetase , 2007, Proteins.

[23]  John SantaLucia,et al.  AMBER Force Field Parameters for the Naturally Occurring Modified Nucleosides in RNA. , 2007, Journal of chemical theory and computation.

[24]  Lennart Nilsson,et al.  Molecular Dynamics Study of Intrinsic Stability in Six RNA Terminal Loop Motifs , 2007, Journal of biomolecular structure & dynamics.

[25]  B. Sproat,et al.  Specificity of Phage Display Selected Peptides for Modified Anticodon Stem and Loop Domains of tRNA , 2007, The protein journal.

[26]  D. Davis,et al.  Structural effects of hypermodified nucleosides in the Escherichia coli and human tRNALys anticodon loop: the effect of nucleosides s2U, mcm5U, mcm5s2U, mnm5s2U, t6A, and ms2t6A. , 2005, Biochemistry.

[27]  C. Tisné.  Structural bases of the annealing of primer tRNA(3Lys) to the HIV-1 viral RNA. , 2005, Current HIV research.

[28]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[29]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[30]  B. Roques,et al.  Specific recognition of primer tRNA Lys 3 by HIV-1 nucleocapsid protein: involvement of the zinc fingers and the N-terminal basic extension. , 2003, Biochimie.

[31]  D. Myszka,et al.  An RNA complex of the HIV-1 A-loop and tRNA(Lys,3) is stabilized by nucleoside modifications. , 2002, Journal of the American Chemical Society.

[32]  C. Ehresmann,et al.  The crystal structure of HIV reverse-transcription primer tRNA(Lys,3) shows a canonical anticodon loop. , 2000, RNA.

[33]  D. Davis,et al.  Hypermodified nucleosides in the anticodon of tRNALys stabilize a canonical U-turn structure. , 2000, Biochemistry.

[34]  M. Karplus,et al.  Analysis of Calculated Normal Modes of a Set of Native and Partially Unfolded Proteins , 1999 .

[35]  C. Veltri,et al.  An RNA model system for investigation of pseudouridine stabilization of the codon-anticodon interaction in tRNALys, tRNAHis and tRNATyr. , 1998, Journal of biomolecular structure & dynamics.

[36]  P. Koehl,et al.  Application of a self-consistent mean field theory to predict protein side-chains conformation and estimate their conformational entropy. , 1994, Journal of molecular biology.

[37]  D. Theodorou,et al.  A concerted rotation algorithm for atomistic Monte Carlo simulation of polymer melts and glasses , 1993 .

[38]  Derek W. Robinson,et al.  Mean entropy of states in classical statistical mechanics , 1967 .

[39]  Kentaro Shimizu,et al.  Mechanism of the difference in the binding affinity of E. coli tRNAGln to glutaminyl-tRNA synthetase caused by noninterface nucleotides in variable loop. , 2007, Biophysical journal.

[40]  C. Ehresmann,et al.  tRNAs as primer of reverse transcriptases. , 1995, Biochimie.