Adding energy minimization strategy to peptide‐design algorithm enables better search for RNA‐binding peptides: Redesigned λ N peptide binds boxB RNA

Our previously developed peptide‐design algorithm was improved by adding an energy minimization strategy which allows the amino acid sidechains to move in a broad configuration space during sequence evolution. In this work, the new algorithm was used to generate a library of 21‐mer peptides which could substitute for λ N peptide in binding to boxB RNA. Six potential peptides were obtained from the algorithm, all of which exhibited good binding capability with boxB RNA. Atomistic molecular dynamics simulations were then conducted to examine the ability of the λ N peptide and three best evolved peptides, viz. Pept01, Pept26, and Pept28, to bind to boxB RNA. Simulation results demonstrated that our evolved peptides are better at binding to boxB RNA than the λ N peptide. Sequence searches using the old (without energy minimization strategy) and new (with energy minimization strategy) algorithms confirm that the new algorithm is more effective at finding good RNA‐binding peptides than the old algorithm. © 2016 Wiley Periodicals, Inc.

[1]  J. Williamson,et al.  Analysis of bacteriophage N protein and peptide binding to boxB RNA using polyacrylamide gel coelectrophoresis (PACE). , 1997, RNA.

[2]  H. Heus,et al.  Structural features that give rise to the unusual stability of RNA hairpins containing GNRA loops. , 1991, Science.

[3]  D. Baker,et al.  Elicitation of structure-specific antibodies by epitope scaffolds , 2010, Proceedings of the National Academy of Sciences.

[4]  Bruce Tidor,et al.  Progress in computational protein design. , 2007, Current opinion in biotechnology.

[5]  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.

[6]  B. Stoddard,et al.  Computational Thermostabilization of an Enzyme , 2005, Science.

[7]  C. G. Broyden The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations , 1970 .

[8]  S. Phillips,et al.  The crystal structure of a high affinity RNA stem-loop complexed with the bacteriophage MS2 capsid: further challenges in the modeling of ligand-RNA interactions. , 2004, RNA.

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

[10]  Eric A. Althoff,et al.  De Novo Computational Design of Retro-Aldol Enzymes , 2008, Science.

[11]  D. Goldfarb A family of variable-metric methods derived by variational means , 1970 .

[12]  C. Barreau,et al.  Tethering of proteins to RNAs by bacteriophage proteins , 2008, Biology of the cell.

[13]  D. Shanno Conditioning of Quasi-Newton Methods for Function Minimization , 1970 .

[14]  B. Kuhlman,et al.  Computational design of a single amino acid sequence that can switch between two distinct protein folds. , 2006, Journal of the American Chemical Society.

[15]  J. Richardson,et al.  The penultimate rotamer library , 2000, Proteins.

[16]  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.

[17]  L. Liljas,et al.  Investigating the structural basis of purine specificity in the structures of MS2 coat protein RNA translational operator hairpins. , 2002, Nucleic acids research.

[18]  L. Stamatatos,et al.  Computational design of epitope-scaffolds allows induction of antibodies specific for a poorly immunogenic HIV vaccine epitope. , 2010, Structure.

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

[20]  Xavier Darzacq,et al.  Imaging gene expression in single living cells , 2004, Nature Reviews Molecular Cell Biology.

[21]  T. Xia,et al.  Conformational distributions at the N-peptide/boxB RNA interface studied using site-directed spin labeling. , 2010, RNA.

[22]  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.

[23]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[24]  J. Murray,et al.  The three-dimensional structures of two complexes between recombinant MS2 capsids and RNA operator fragments reveal sequence-specific protein-RNA interactions. , 1997, Journal of molecular biology.

[25]  Camille J. Delebecque,et al.  Designing and using RNA scaffolds to assemble proteins in vivo , 2012, Nature Protocols.

[26]  Christopher M. MacDermaid,et al.  Theoretical and computational protein design. , 2011, Annual review of physical chemistry.

[27]  A. Zewail,et al.  RNA-protein recognition: single-residue ultrafast dynamical control of structural specificity and function. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

[29]  Jinming Li,et al.  A novel delivery platform based on Bacteriophage MS2 virus-like particles , 2015, Virus Research.

[30]  D. Lazinski,et al.  Sequence-specific recognition of RNA hairpins by bacteriophage antiterminators requires a conserved arginine-rich motif , 1989, Cell.

[31]  J. Murray,et al.  Crystal structures of MS2 coat protein mutants in complex with wild-type RNA operator fragments. , 1998, Nucleic acids research.

[32]  R. Fletcher,et al.  A New Approach to Variable Metric Algorithms , 1970, Comput. J..

[33]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[34]  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.

[35]  Wei Wang,et al.  Computational de novo design and characterization of a four-helix bundle protein that selectively binds a nonbiological cofactor. , 2005, Journal of the American Chemical Society.

[36]  Jack Greenblatt,et al.  NMR Structure of the Bacteriophage λ N Peptide/boxB RNA Complex: Recognition of a GNRA Fold by an Arginine-Rich Motif , 1998, Cell.

[37]  M. Boehm,et al.  Antitermination in bacteriophage lambda. The structure of the N36 peptide-boxB RNA complex. , 2000, European journal of biochemistry.

[38]  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.

[39]  J. Murray,et al.  Crystallographic studies of RNA hairpins in complexes with recombinant MS2 capsids: implications for binding requirements. , 1999, RNA.

[40]  Mookyung Cheon,et al.  Extending the PRIME model for protein aggregation to all 20 amino acids , 2010, Proteins.

[41]  R. Friesner,et al.  The VSGB 2.0 model: A next generation energy model for high resolution protein structure modeling , 2011, Proteins.