Investigating the binding specificity of U1A-RNA by computational mutagenesis.

The mammalian spliceosomal protein U1A binds a hairpin RNA with picomolar affinity. To examine the origin of this binding specificity, we carried out computational mutagenesis on protein and RNA residues in the U1A-RNA binding interface. Our computational mutagenesis methods calculate the relative binding affinity between mutant and wild-type as the sum of molecular mechanical energies and solvation free energies estimated with a continuum solvent model. We obtained good agreement with experimental studies and we verified mutations that abolish and improve binding. Therefore, we offer these methods as computationally inexpensive tools for investigating and predicting the effects of site-specific mutagenesis.

[1]  Conrad C. Huang,et al.  The MIDAS display system , 1988 .

[2]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[3]  P. Evans,et al.  Crystal structure of the RNA-binding domain of the U1 small nuclear ribonucleoprotein A , 1990, Nature.

[4]  J. Keene,et al.  RNA recognition: towards identifying determinants of specificity. , 1991, Trends in biochemical sciences.

[5]  C. Oubridge,et al.  Identification of molecular contacts between the U1 A small nuclear ribonucleoprotein and U1 RNA. , 1991, The EMBO journal.

[6]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[7]  Nobutoshi Ito,et al.  Crystal structure at 1.92 Å resolution of the RNA-binding domain of the U1A spliceosomal protein complexed with an RNA hairpin , 1994, Nature.

[8]  K. Hall,et al.  Interaction of RNA hairpins with the human U1A N-terminal RNA binding domain. , 1994, Biochemistry.

[9]  B. Honig,et al.  Evaluation of the conformational free energies of loops in proteins , 1994, Proteins.

[10]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[11]  P. Evans,et al.  The RNP domain: a sequence-specific RNA-binding domain involved in processing and transport of RNA. , 1995, Trends in biochemical sciences.

[12]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules J. Am. Chem. Soc. 1995, 117, 5179−5197 , 1996 .

[13]  Jian Shen,et al.  Calculation of binding energy differences for receptor–ligand systems using the Poisson‐Boltzmann method , 1995, J. Comput. Chem..

[14]  J. Belasco,et al.  Analysis of RNA-binding proteins by in vitro genetic selection: identification of an amino acid residue important for locking U1A onto its RNA target. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[15]  K. Hall,et al.  [12] Thermodynamics and mutations in RNA—Protein interactions , 1995 .

[16]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[17]  K. Hall,et al.  Contribution of the tyrosines to the structure and function of the human U1A N‐terminal RNA binding domain , 1996, Protein science : a publication of the Protein Society.

[18]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[19]  G. Varani,et al.  Solution structure of the N-terminal RNP domain of U1A protein: the role of C-terminal residues in structure stability and RNA binding. , 1996, Journal of molecular biology.

[20]  John J. Wendoloski,et al.  Electrostatic binding energy calculation using the finite difference solution to the linearized Poisson‐Boltzmann equation: Assessment of its accuracy , 1996 .

[21]  G. Varani,et al.  Specificity of ribonucleoprotein interaction determined by RNA folding during complex formation , 1996, Nature.

[22]  D E Koshland,et al.  Computational method for relative binding energies of enzyme‐substrate complexes , 1996, Protein science : a publication of the Protein Society.

[23]  K. Hall,et al.  RNA hairpins with non-nucleotide spacers bind efficiently to the human U1A protein. , 1996, Journal of molecular biology.

[24]  B Honig,et al.  On the calculation of binding free energies using continuum methods: Application to MHC class I protein‐peptide interactions , 1997, Protein science : a publication of the Protein Society.

[25]  G. Varani,et al.  Structural basis of the RNA‐binding specificity of human U1A protein , 1997, The EMBO journal.

[26]  Contribution of the C-terminal tail of U1A RBD1 to RNA recognition and protein stability. , 1997, RNA.

[27]  G. Varani,et al.  RNA recognition by RNP proteins during RNA processing. , 1998, Annual review of biophysics and biomolecular structure.

[28]  P A Kollman,et al.  Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution. , 1998, Journal of biomolecular structure & dynamics.

[29]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[30]  K. Hall,et al.  RNA binding mediates the local cooperativity between the beta-sheet and the C-terminal tail of the human U1A RBD1 protein. , 1998, Journal of molecular biology.

[31]  P A Kollman,et al.  Continuum solvent studies of the stability of RNA hairpin loops and helices. , 1998, Journal of biomolecular structure & dynamics.

[32]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[33]  K. Hall,et al.  RNA recognition by the human U1A protein is mediated by a network of local cooperative interactions that create the optimal binding surface. , 1999, Journal of molecular biology.

[34]  P. Kollman,et al.  Molecular dynamics studies of U1A-RNA complexes. , 1999, RNA.