1
|
Davis CM, Gruebele M. Cellular Sticking Can Strongly Reduce Complex Binding by Speeding Dissociation. J Phys Chem B 2021; 125:3815-3823. [PMID: 33826329 DOI: 10.1021/acs.jpcb.1c00950] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
While extensive studies have been carried out to determine protein-RNA binding affinities, mechanisms, and dynamics in vitro, such studies do not take into consideration the effect of the many weak nonspecific interactions in a cell filled with potential binding partners. Here we experimentally tested the role of the cellular environment on affinity and binding dynamics between a protein and RNA in living U-2 OS cells. Our model system is the spliceosomal protein U1A and its binding partner SL2 of the U1 snRNA. The binding equilibrium was perturbed by a laser-induced temperature jump and monitored by Förster resonance energy transfer. The apparent binding affinity in live cells was reduced by up to 2 orders of magnitude compared to in vitro. The measured in-cell dissociation rate coefficients were up to 2 orders of magnitude larger, whereas no change in the measured association rate coefficient was observed. The latter is not what would be anticipated due to macromolecular crowding or nonspecific sticking of the uncomplexed U1A and SL2 in the cell. A quantitative model fits our experimental results, with the major cellular effect being that U1A and SL2 sticking to cellular components are capable of binding, just not as strongly as the free complex. This observation suggests that high binding affinities measured or designed in vitro are necessary for proper binding in vivo, where competition with many nonspecific interactions exists, especially for strongly interacting species with high charge or large hydrophobic surface areas.
Collapse
|
2
|
Shao Q, Gong W, Li C. A study on allosteric communication in U1A-snRNA binding interactions: network analysis combined with molecular dynamics data. Biophys Chem 2020; 264:106393. [PMID: 32653695 DOI: 10.1016/j.bpc.2020.106393] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2019] [Revised: 05/02/2020] [Accepted: 05/03/2020] [Indexed: 01/09/2023]
Abstract
The allosteric regulation during the binding interactions between small nuclear RNAs (snRNAs) and the associated protein factors is critical to the function of spliceosomes in alternative RNA splicing. Although network models combined with molecular dynamics simulations have shown to be powerful tools for the analysis of protein allostery, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems. In this work, we use a residual network model combined with a coarse-grained Gaussian network model (GNM) to investigate the binding interactions between the snRNA and the human U1A protein which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle, and to identify the residues that play an important role in the allosteric communication in U1A during this process. We also utilize the Girvan-Newman method to detect the structural organization in U1A-snRNA recognition and interactions. Our results reveal that: (Ι) not only the residues at the binding sites that are traditionally considered to play a major role in U1A-snRNA association, but those residues that are far away from the RNA binding interface participate in the U1A's allosteric signal transmission induced by the RNA binding; (Π) the structure of U1A protein is well organized with different communities acting different roles for its RNA binding and allosteric regulation. The study demonstrates that the combination of the residual network and elastic network models is an effective and efficient method which can be readily extended to the investigation of the allosteric communication for other macromolecular interaction systems.
Collapse
Affiliation(s)
- Qi Shao
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Weikang Gong
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China.
| |
Collapse
|
3
|
Sharma M, Sharma S, Alawada A. Understanding the binding specificities of mRNA targets by the mammalian Quaking protein. Nucleic Acids Res 2020; 47:10564-10579. [PMID: 31602485 PMCID: PMC6847458 DOI: 10.1093/nar/gkz877] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Revised: 09/18/2019] [Accepted: 10/08/2019] [Indexed: 01/01/2023] Open
Abstract
Mammalian Quaking (QKI) protein, a member of STAR family of proteins is a mRNA-binding protein, which post-transcriptionally modulates the target RNA. QKI protein possesses a maxi-KH domain composed of single heterogeneous nuclear ribonucleoprotein K homology (KH) domain and C-terminal QUA2 domain, that binds a sequence-specific QKI RNA recognition element (QRE), CUAAC. To understand the binding specificities for different mRNA sequences of the KH-QUA2 domain of QKI protein, we introduced point mutations at different positions in the QRE resulting in twelve different mRNA sequences with single nucleotide change. We carried out long unbiased molecular dynamics simulations using two different sets of recently updated forcefield parameters: AMBERff14SB+RNAχOL3 and CHARMM36 (with CMAP correction). We analyzed the changes in intermolecular dynamics as a result of mutation. Our results show that AMBER forcefields performed better to model the interactions between mRNA and protein. We also calculated the binding affinities of different mRNA sequences and found that the relative order correlates to the reported experimental studies. Our study shows that the favorable binding with the formation of stable complex will occur when there is an increase of the total intermolecular contacts between mRNA and protein, but without the loss of native contacts within the KH-QUA domain.
Collapse
Affiliation(s)
- Monika Sharma
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Sector 81, Knowledge City, SAS Nagar, Punjab, India
| | - Shakshi Sharma
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Sector 81, Knowledge City, SAS Nagar, Punjab, India
| | - Apoorv Alawada
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Sector 81, Knowledge City, SAS Nagar, Punjab, India
| |
Collapse
|
4
|
Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model. Biophys J 2019; 116:1625-1636. [PMID: 30975455 DOI: 10.1016/j.bpj.2019.03.008] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Revised: 03/04/2019] [Accepted: 03/12/2019] [Indexed: 12/14/2022] Open
Abstract
The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.
Collapse
|
5
|
Sharma M, Anirudh CR. Mechanism of mRNA-STAR domain interaction: Molecular dynamics simulations of Mammalian Quaking STAR protein. Sci Rep 2017; 7:12567. [PMID: 28974714 PMCID: PMC5626755 DOI: 10.1038/s41598-017-12930-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Accepted: 09/20/2017] [Indexed: 01/08/2023] Open
Abstract
STAR proteins are evolutionary conserved mRNA-binding proteins that post-transcriptionally regulate gene expression at all stages of RNA metabolism. These proteins possess conserved STAR domain that recognizes identical RNA regulatory elements as YUAAY. Recently reported crystal structures show that STAR domain is composed of N-terminal QUA1, K-homology domain (KH) and C-terminal QUA2, and mRNA binding is mediated by KH-QUA2 domain. Here, we present simulation studies done to investigate binding of mRNA to STAR protein, mammalian Quaking protein (QKI). We carried out conventional MD simulations of STAR domain in presence and absence of mRNA, and studied the impact of mRNA on the stability, dynamics and underlying allosteric mechanism of STAR domain. Our unbiased simulations results show that presence of mRNA stabilizes the overall STAR domain by reducing the structural deviations, correlating the ‘within-domain’ motions, and maintaining the native contacts information. Absence of mRNA not only influenced the essential modes of motion of STAR domain, but also affected the connectivity of networks within STAR domain. We further explored the dissociation of mRNA from STAR domain using umbrella sampling simulations, and the results suggest that mRNA binding to STAR domain occurs in multi-step: first conformational selection of mRNA backbone conformations, followed by induced fit mechanism as nucleobases interact with STAR domain.
Collapse
Affiliation(s)
- Monika Sharma
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Sector 81, Knowledge City, SAS Nagar, Punjab, India.
| | - C R Anirudh
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER), Sector 81, Knowledge City, SAS Nagar, Punjab, India
| |
Collapse
|
6
|
A compare-and-contrast NMR dynamics study of two related RRMs: U1A and SNF. Biophys J 2015; 107:208-19. [PMID: 24988355 DOI: 10.1016/j.bpj.2014.05.026] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2013] [Revised: 05/14/2014] [Accepted: 05/19/2014] [Indexed: 11/21/2022] Open
Abstract
The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/β sandwich topology, and these RRMs use their β-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the β-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 β-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.
Collapse
|
7
|
Williams SG, Harms MJ, Hall KB. Resurrection of an Urbilaterian U1A/U2B″/SNF protein. J Mol Biol 2013; 425:3846-62. [PMID: 23796518 DOI: 10.1016/j.jmb.2013.05.031] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2013] [Revised: 05/06/2013] [Accepted: 05/08/2013] [Indexed: 11/20/2022]
Abstract
The U1A/U2B″/SNF family of proteins found in the U1 and U2 spliceosomal small nuclear ribonucleoproteins is highly conserved. In spite of the high degree of sequence and structural conservation, modern members of this protein family have unique RNA binding properties. These differences have necessarily resulted from evolutionary processes, and therefore, we reconstructed the protein phylogeny in order to understand how and when divergence occurred and how protein function has been modulated. Contrary to the conventional understanding of an ancient human U1A/U2B″ gene duplication, we show that the last common ancestor of bilaterians contained a single ancestral protein (URB). The gene for URB was synthesized, the protein was overexpressed and purified, and we assessed RNA binding to modern snRNA sequences. We find that URB binds human and Drosophila U1 snRNA SLII and U2 snRNA SLIV with higher affinity than do modern homologs, suggesting that both Drosophila SNF and human U1A/U2B″ have evolved into weaker binders of one RNA or both RNAs.
Collapse
MESH Headings
- Amino Acid Sequence
- Animals
- Base Sequence
- Drosophila
- Drosophila Proteins/chemistry
- Drosophila Proteins/genetics
- Drosophila Proteins/metabolism
- Evolution, Molecular
- Gene Duplication
- Humans
- Inverted Repeat Sequences
- Models, Molecular
- Molecular Sequence Data
- Nucleic Acid Conformation
- Phylogeny
- Protein Binding
- Protein Conformation
- RNA, Small Nuclear/chemistry
- RNA, Small Nuclear/genetics
- RNA, Small Nuclear/metabolism
- Ribonucleoprotein, U1 Small Nuclear/chemistry
- Ribonucleoprotein, U1 Small Nuclear/genetics
- Ribonucleoprotein, U1 Small Nuclear/metabolism
- Ribonucleoproteins, Small Nuclear/chemistry
- Ribonucleoproteins, Small Nuclear/genetics
- Ribonucleoproteins, Small Nuclear/metabolism
- Sequence Alignment
- Spliceosomes/metabolism
- snRNP Core Proteins/chemistry
- snRNP Core Proteins/genetics
- snRNP Core Proteins/metabolism
Collapse
Affiliation(s)
- Sandra G Williams
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63108, USA
| | | | | |
Collapse
|
8
|
Rau M, Stump WT, Hall KB. Intrinsic flexibility of snRNA hairpin loops facilitates protein binding. RNA (NEW YORK, N.Y.) 2012; 18:1984-1995. [PMID: 23012481 PMCID: PMC3479389 DOI: 10.1261/rna.035006.112] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2012] [Accepted: 08/02/2012] [Indexed: 06/01/2023]
Abstract
Stem-loop II of U1 snRNA and Stem-loop IV of U2 snRNA typically have 10 or 11 nucleotides in their loops. The fluorescent nucleobase 2-aminopurine was used as a substitute for the adenines in each loop to probe the local and global structures and dynamics of these unusually long loops. Using steady-state and time-resolved fluorescence, we find that, while the bases in the loops are stacked, they are able to undergo significant local motion on the picosecond/nanosecond timescale. In addition, the loops have a global conformational change at low temperatures that occurs on the microsecond timescale, as determined using laser T-jump experiments. Nucleobase and loop motions are present at temperatures far below the melting temperature of the hairpin stem, which may facilitate the conformational change required for specific protein binding to these RNA loops.
Collapse
Affiliation(s)
- Michael Rau
- Department of Biochemistry and Molecular Biophysics, Washington University Medical School, St. Louis, Missouri 63110, USA
| | - W. Tom Stump
- Department of Biochemistry and Molecular Biophysics, Washington University Medical School, St. Louis, Missouri 63110, USA
| | - Kathleen B. Hall
- Department of Biochemistry and Molecular Biophysics, Washington University Medical School, St. Louis, Missouri 63110, USA
| |
Collapse
|
9
|
Bothe JR, Nikolova EN, Eichhorn CD, Chugh J, Hansen AL, Al-Hashimi HM. Characterizing RNA dynamics at atomic resolution using solution-state NMR spectroscopy. Nat Methods 2011; 8:919-31. [PMID: 22036746 PMCID: PMC3320163 DOI: 10.1038/nmeth.1735] [Citation(s) in RCA: 112] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Many recently discovered noncoding RNAs do not fold into a single native conformation but sample many different conformations along their free-energy landscape to carry out their biological function. Here we review solution-state NMR techniques that measure the structural, kinetic and thermodynamic characteristics of RNA motions spanning picosecond to second timescales at atomic resolution, allowing unprecedented insights into the RNA dynamic structure landscape. From these studies a basic description of the RNA dynamic structure landscape is emerging, bringing new insights into how RNA structures change to carry out their function as well as applications in RNA-targeted drug discovery and RNA bioengineering.
Collapse
Affiliation(s)
- Jameson R. Bothe
- Department of Chemistry, The University of Michigan, Ann Arbor, Michigan, USA
| | - Evgenia N. Nikolova
- Chemical Biology Doctoral Program, The University of Michigan, Ann Arbor, Michigan, USA
| | - Catherine D. Eichhorn
- Chemical Biology Doctoral Program, The University of Michigan, Ann Arbor, Michigan, USA
| | - Jeetender Chugh
- Department of Biophysics, The University of Michigan, Ann Arbor, Michigan, USA
| | - Alexandar L. Hansen
- Department of Chemistry, The University of Toronto, Toronto, Ontario, Canada
- Department of Biochemistry, The University of Toronto, Toronto, Ontario, Canada
- Department of Molecular Genetics, The University of Toronto, Toronto, Ontario, Canada
| | - Hashim M. Al-Hashimi
- Department of Chemistry, The University of Michigan, Ann Arbor, Michigan, USA
- Department of Biophysics, The University of Michigan, Ann Arbor, Michigan, USA
| |
Collapse
|
10
|
Kormos BL, Pieniazek SN, Beveridge DL, Baranger AM. U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations. Biopolymers 2011; 95:591-606. [PMID: 21384338 DOI: 10.1002/bip.21616] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2010] [Revised: 01/01/2011] [Accepted: 02/07/2011] [Indexed: 01/02/2023]
Abstract
Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.
Collapse
Affiliation(s)
- Bethany L Kormos
- Chemistry Department and Molecular Biophysics Program, Wesleyan University, Middletown, CT 06459, USA
| | | | | | | |
Collapse
|
11
|
Wostenberg C, Noid WG, Showalter SA. MD simulations of the dsRBP DGCR8 reveal correlated motions that may aid pri-miRNA binding. Biophys J 2010; 99:248-56. [PMID: 20655853 DOI: 10.1016/j.bpj.2010.04.010] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2009] [Revised: 03/12/2010] [Accepted: 04/01/2010] [Indexed: 12/13/2022] Open
Abstract
Over the past decade, microRNAs (miRNAs) have been shown to affect gene regulation by basepairing with messenger RNA, and their misregulation has been directly linked with cancer. DGCR8, a protein that contains two dsRNA-binding domains (dsRBDs) in tandem, is vital for nuclear maturation of primary miRNAs (pri-miRNAs) in connection with the RNase III enzyme Drosha. The crystal structure of the DGCR8 Core (493-720) shows a unique, well-ordered structure of the linker region between the two dsRBDs that differs from the flexible linker connecting the two dsRBDs in the antiviral response protein, PKR. To better understand the interfacial interactions between the two dsRBDs, we ran extensive MD simulations of isolated dsRBDs (505-583 and 614-691) and the Core. The simulations reveal correlated reorientations of the two domains relative to one another, with the well-ordered linker and C-terminus serving as a pivot. The results demonstrate that motions at the domain interface dynamically impact the conformation of the RNA-binding surface and may provide an adaptive separation distance that is necessary to allow interactions with a variety of different pri-miRNAs with heterogeneous structures. These results thus provide an entry point for further in vitro studies of the potentially unique RNA-binding mode of DGCR8.
Collapse
Affiliation(s)
- Christopher Wostenberg
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania, USA
| | | | | |
Collapse
|
12
|
Qin F, Chen Y, Wu M, Li Y, Zhang J, Chen HF. Induced fit or conformational selection for RNA/U1A folding. RNA (NEW YORK, N.Y.) 2010; 16:1053-1061. [PMID: 20354153 PMCID: PMC2856877 DOI: 10.1261/rna.2008110] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2009] [Accepted: 02/07/2010] [Indexed: 05/29/2023]
Abstract
The hairpin II of U1 snRNA can bind U1A protein with high affinity and specificity. NMR spectra suggest that the loop region of apo-RNA is largely unstructured and undergoes a transition from unstructured to well-folded upon U1Abinding. However, the mechanism that RNA folding coupled protein binding is poorly understood. To get an insight into the mechanism, we have performed explicit-solvent molecular dynamics (MD) to study the folding kinetics of bound RNA and apo-RNA. Room-temperature MD simulations suggest that the conformation of bound RNA has significant adjustment and becomes more stable upon U1A binding. Kinetic analysis of high-temperature MD simulations shows that bound RNA and apo-RNA unfold via a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound RNA folds in the order of RNA contracting, U1A binding, and tertiary folding. The predicted Phi-values suggest that A8, C10, A11, and G16 are key bases for bound RNA folding. Mutant Arg52Gln analysis shows that electrostatic interaction and hydrogen bonds between RNA and U1A (Arg52Gln) decrease. These results are in qualitative agreement with experiments. Furthermore, this method could be used in other studies about biomolecule folding upon receptor binding.
Collapse
Affiliation(s)
- Fang Qin
- College of Life Sciences and Biotechnology, Shanghai Jiaotong University, Shanghai, 200240, China
| | | | | | | | | | | |
Collapse
|
13
|
Musselman C, Zhang Q, Al-Hashimi H, Andricioaei I. Referencing strategy for the direct comparison of nuclear magnetic resonance and molecular dynamics motional parameters in RNA. J Phys Chem B 2010; 114:929-39. [PMID: 20039757 PMCID: PMC4287414 DOI: 10.1021/jp905286h] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations are both techniques that can be used to characterize the structural dynamics of biomolecules and their underlying time scales. Comparison of relaxation parameters obtained through each methodology allows for cross validation of techniques and for complementarity in the analysis of dynamics. Here we present a combined NMR/MD study of the dynamics of HIV-1 transactivation response (TAR) RNA. We compute relaxation constants (R(1), R(2), and NOE) and model-free parameters (S(2) and tau) from a 65 ns molecular dynamics (MD) trajectory and compare them with the respective parameters measured in a domain-elongation NMR experiment. Using the elongated domain as the frame of reference for all computed parameters allows for a direct comparison between experiment and simulation. We see good agreement for many parameters and gain further insight into the nature of the local and global dynamics of TAR, which are found to be quite complex, spanning multiple time scales. For the few cases where agreement is poor, comparison of the dynamical parameters provides insight into the limits of each technique. We suggest a frequency-matching procedure that yields an upper bound for the time scale of dynamics to which the NMR relaxation experiment is sensitive.
Collapse
Affiliation(s)
- Catherine Musselman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | | | | |
Collapse
|
14
|
Interactions between PTB RRMs induce slow motions and increase RNA binding affinity. J Mol Biol 2010; 397:260-77. [PMID: 20080103 DOI: 10.1016/j.jmb.2009.12.051] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2009] [Revised: 11/19/2009] [Accepted: 12/24/2009] [Indexed: 11/23/2022]
Abstract
Polypyrimidine tract binding protein (PTB) participates in a variety of functions in eukaryotic cells, including alternative splicing, mRNA stabilization, and internal ribosomal entry site-mediated translation initiation. Its mechanism of RNA recognition is determined in part by the novel geometry of its two C-terminal RNA recognition motifs (RRM3 and RRM4), which interact with each other to form a stable complex (PTB1:34). This complex itself is unusual among RRMs, suggesting that it performs a specific function for the protein. In order to understand the advantage it provides to PTB, the fundamental properties of PTB1:34 are examined here as a comparative study of the complex and its two constituent RRMs. Both RRM3 and RRM4 adopt folded structures that NMR data show to be similar to their structure in PRB1:34. The RNA binding properties of the domains differ dramatically. The affinity of each separate RRM for polypyrimidine tracts is far weaker than that of PTB1:34, and simply mixing the two RRMs does not create an equivalent binding platform. (15)N NMR relaxation experiments show that PTB1:34 has slow, microsecond motions throughout both RRMs including the interdomain linker. This is in contrast to the individual domains, RRM3 and RRM4, where only a few backbone amides are flexible on this time scale. The slow backbone dynamics of PTB1:34, induced by packing of RRM3 and RRM4, could be essential for high-affinity binding to a flexible polypyrimidine tract RNA and also provide entropic compensation for its own formation.
Collapse
|
15
|
Julien KR, Sumita M, Chen PH, Laird-Offringa IA, Hoogstraten CG. Conformationally restricted nucleotides as a probe of structure-function relationships in RNA. RNA (NEW YORK, N.Y.) 2008; 14:1632-1643. [PMID: 18596252 PMCID: PMC2491483 DOI: 10.1261/rna.866408] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2008] [Accepted: 04/17/2008] [Indexed: 05/26/2023]
Abstract
We introduce the use of commercially available locked nucleic acids (LNAs) as a functional probe in RNA. LNA nucleotides contain a covalent linkage that restricts the pseudorotation phase of the ribose to C3'-endo (A-form). Introduction of an LNA at a single site thus allows the role of ribose structure and dynamics in RNA function to be assessed. We apply LNA probing at multiple sites to analyze self-cleavage in the lead-dependent ribozyme (leadzyme), thermodynamic stability in the UUCG tetraloop, and the kinetics of recognition of U1A protein by U1 snRNA hairpin II. In the leadzyme, locking a single guanosine residue into the C3'-endo pucker increases the catalytic rate by a factor of 20, despite the fact that X-ray crystallographic and NMR structures of the leadzyme ground state reported a C2'-endo conformation at this site. These results strongly suggest that a conformational change at this position is critical for catalytic function. Functional insights obtained in all three systems demonstrate the highly general applicability of LNA probing in analysis of the role of ribose orientation in RNA structure, dynamics, and function.
Collapse
Affiliation(s)
- Kristine R Julien
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
| | | | | | | | | |
Collapse
|
16
|
Pressure effects on the ensemble dynamics of ubiquitin inspected with molecular dynamics simulations and isotropic reorientational eigenmode dynamics. Biophys J 2008; 95:3943-55. [PMID: 18621814 DOI: 10.1529/biophysj.108.133702] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
According to NMR chemical shift data, the ensemble of ubiquitin is a mixture of "open" and "closed" conformations at rapid equilibrium. Pressure perturbations provide the means to study the transition between the two conformers by imposing an additional constraint on the system's partial molar volume. Here we use nanosecond-timescale molecular dynamics simulations to characterize the network of correlated motions accessible to the conformers at low- and high-pressure conditions. Using the isotropic reorientational eigenmode dynamics formalism to analyze our simulation trajectories, we reproduce NMR relaxation data without fitting any parameters of our model. Comparative analysis of our results suggests that the two conformations behave very differently. The dynamics of the "closed" conformation are almost unaffected by pressure and are dominated by large-amplitude correlated motions of residues 23-34 in the extended alpha-helix. The "open" conformation under conditions of normal pressure displays increased mobility, focused on the loop residues 17-20, 46-55, and 58-59 at the bottom of the core of the structure, as well as the C-terminal residues 69-76, that directly participate in key protein-protein interactions. For the same conformation, a pressure increase induces a loss of separability between molecular tumbling and internal dynamics, while motions between different backbone sites become uncorrelated.
Collapse
|
17
|
Anunciado D, Agumeh M, Kormos BL, Beveridge DL, Knee JL, Baranger AM. Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements. J Phys Chem B 2008; 112:6122-30. [PMID: 18293956 DOI: 10.1021/jp076896c] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.
Collapse
Affiliation(s)
- Divina Anunciado
- Department of Chemistry and Molecular Biophysics Program, Wesleyan University, Middletown, Connecticut 06459, USA
| | | | | | | | | | | |
Collapse
|
18
|
Mackerell AD, Nilsson L. Molecular dynamics simulations of nucleic acid-protein complexes. Curr Opin Struct Biol 2008; 18:194-9. [PMID: 18281210 DOI: 10.1016/j.sbi.2007.12.012] [Citation(s) in RCA: 130] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2007] [Revised: 12/17/2007] [Accepted: 12/21/2007] [Indexed: 10/22/2022]
Abstract
Molecular dynamics simulation studies of protein-nucleic acid complexes are more complicated than studies of either component alone-the force field has to be properly balanced, the systems tend to become very large, and a careful treatment of solvent and of electrostatic interactions is necessary. Recent investigations into several protein-DNA and protein-RNA systems have shown the feasibility of the simulation approach, yielding results of biological interest not readily accessible to experimental methods.
Collapse
Affiliation(s)
- Alexander D Mackerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA.
| | | |
Collapse
|
19
|
Recognition of essential purines by the U1A protein. BMC BIOCHEMISTRY 2007; 8:22. [PMID: 17980039 PMCID: PMC2203988 DOI: 10.1186/1471-2091-8-22] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/16/2007] [Accepted: 11/02/2007] [Indexed: 11/10/2022]
Abstract
Background The RNA recognition motif (RRM) is one of the largest families of RNA binding domains. The RRM is modulated so that individual proteins containing RRMs can specifically recognize RNA targets with diverse sequences and structures. Understanding the principles governing this specificity will be important for the rational modification and design of RRM-RNA complexes. Results In this paper we have investigated the origins of specificity of the N terminal RRM of the U1A protein for stem loop 2 (SL2) of U1 snRNA by substituting modified bases for essential purines in SL2 RNA. In one series of modified bases, hydrogen bond donors and acceptors were replaced by aliphatic groups to probe the importance of these functional groups to binding. In a second series of modified bases, hydrogen bond donors and acceptors were incorrectly placed on the purine bases to analyze the origins of discrimination between cognate and non-cognate RNA. The results of these experiments show that three different approaches are used by the U1A protein to gain specificity for purines. Specificity for the first base in the loop, A1, is based primarily on discrimination against RNA containing the incorrect base, specificity for the fourth base in the loop, G4, is based largely on recognition of the donors and acceptors of G4, while specificity for the sixth base in the loop, A6, results from a combination of direct recognition of the base and discrimination against incorrectly placed functional groups. Conclusion These investigations identify different roles that hydrogen bond donors and acceptors on bases in both cognate and non-cognate RNA play in the specific recognition of RNA by the U1A protein. Taken together with investigations of other RNA-RRM complexes, the results contribute to a general understanding of the origins of RNA-RRM specificity and highlight, in particular, the contribution of steric and electrostatic repulsion to binding specificity.
Collapse
|
20
|
Ghosh A, Vishveshwara S. A study of communication pathways in methionyl- tRNA synthetase by molecular dynamics simulations and structure network analysis. Proc Natl Acad Sci U S A 2007; 104:15711-6. [PMID: 17898174 PMCID: PMC2000407 DOI: 10.1073/pnas.0704459104] [Citation(s) in RCA: 186] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by approximately 70 A in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.
Collapse
Affiliation(s)
- Amit Ghosh
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
| | - Saraswathi Vishveshwara
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India
- *To whom correspondence may be addressed. E-mail:
| |
Collapse
|
21
|
Kormos BL, Baranger AM, Beveridge DL. Do collective atomic fluctuations account for cooperative effects? Molecular dynamics studies of the U1A-RNA complex. J Am Chem Soc 2007; 128:8992-3. [PMID: 16834346 PMCID: PMC2603296 DOI: 10.1021/ja0606071] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A complete understanding of gene expression relies on a comprehensive understanding of the protein-RNA recognition process. However, the study of protein-RNA recognition is complicated by many factors that contribute to both binding affinity and specificity, including structure, energetics, dynamical motions, and cooperative interactions. Several recent studies have suggested that energetic coupling between residues contributes to formation of the complex between the U1A protein and stem loop 2 of U1 snRNA as a consequence of a cooperative network of interactions. We have performed molecular dynamics simulations on the U1A-RNA complex, including explicit water and counterions, and analyzed the results based on the calculated positional cross-correlations of atomic fluctuations. The results indicate that cross-correlations calculated on a per residue basis agree well with the observed inter-residue cooperativity and predict that the networks identified to date may also be coupled into an extensive hyper-network that reflects the intrinsic rigidity of the RNA recognition motif. In addition, we report a comparison of the MD calculated correlations with the results of a positional covariance analysis based on the sequences of 330 RNA recognition motifs, including U1A. The calculated inter-residue cross-correlations agree very well with the results of the sites exhibiting positional covariance. Collectively, these results strongly support the hypothesis that collective fluctuations contribute to cooperativity and the corresponding observed thermodynamic coupling. Predictions of additional sites in U1A that may be involved in cooperative networks are advanced.
Collapse
|
22
|
Kormos BL, Benitex Y, Baranger AM, Beveridge DL. Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model. J Mol Biol 2007; 371:1405-19. [PMID: 17603075 PMCID: PMC2034351 DOI: 10.1016/j.jmb.2007.06.003] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2007] [Revised: 06/02/2007] [Accepted: 06/04/2007] [Indexed: 11/26/2022]
Abstract
An MM-GBSA computational protocol was used to investigate wild-type U1A-RNA and F56 U1A mutant experimental binding free energies. The trend in mutant binding free energies compared to wild-type is well-reproduced. Following application of a linear-response-like equation to scale the various energy components, the binding free energies agree quantitatively with observed experimental values. Conformational adaptation contributes to the binding free energy for both the protein and the RNA in these systems. Small differences in DeltaGs are the result of different and sometimes quite large relative contributions from various energetic components. Residual free energy decomposition indicates differences not only at the site of mutation, but throughout the entire protein. MM-GBSA and ab initio calculations performed on model systems suggest that stacking interactions may nearly, but not completely, account for observed differences in mutant binding affinities. This study indicates that there may be different underlying causes of ostensibly similar experimentally observed binding affinities of different mutants, and thus recommends caution in the interpretation of binding affinities and specificities purely by inspection.
Collapse
Affiliation(s)
- Bethany L Kormos
- Chemistry Department and Molecular Biophysics Program, Wesleyan University, Middletown, CT 06459, USA.
| | | | | | | |
Collapse
|
23
|
Musselman C, Al-Hashimi HM, Andricioaei I. iRED analysis of TAR RNA reveals motional coupling, long-range correlations, and a dynamical hinge. Biophys J 2007; 93:411-22. [PMID: 17449677 PMCID: PMC1896250 DOI: 10.1529/biophysj.107.104620] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The HIV-1 transactivation response RNA element (TAR), which is essential to the lifecycle of the virus, has been suggested, based on NMR and hydrodynamic measurements, to undergo substantial, collective, structural dynamics that are important for its function. To deal with the significant coupling between overall diffusional rotation and internal motion expected to exist in TAR, here we utilize an isotropic reorientational eigenmode dynamics analysis of simulated molecular trajectories to obtain a detailed description of TAR dynamics and an accurately quantified pattern of dynamical correlations. The analysis demonstrates the inseparability of internal and overall motional modes, confirms the existence and reveals the nature of collective domain dynamics, and additionally reveals that the hinge for these motions is centered on residues U23, C24, and C41. Results also indicate the existence of long-range communication between the loop and the core of the RNA, and between the loop and the bulge. Additionally, the isotropic reorientational eigenmode dynamics analysis explains, from a dynamical perspective, several existing biochemical mutational studies and suggests new mutations for future structural dynamics studies.
Collapse
Affiliation(s)
- Catherine Musselman
- Department of Chemistry and The Center for Computational Medicine and Biology, University of Michigan, Ann Arbor, Michigan, USA
| | | | | |
Collapse
|
24
|
McDowell SE, Špačková N, Šponer J, Walter NG. Molecular dynamics simulations of RNA: an in silico single molecule approach. Biopolymers 2007; 85:169-84. [PMID: 17080418 PMCID: PMC2018183 DOI: 10.1002/bip.20620] [Citation(s) in RCA: 123] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .
Collapse
Affiliation(s)
- S. Elizabeth McDowell
- Biophysics Research Division, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| | - Nad'a Špačková
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| |
Collapse
|
25
|
Eldho NV, Dayie KT. Internal Bulge and Tetraloop of the Catalytic Domain 5 of a Group II Intron Ribozyme Are Flexible: Implications for Catalysis. J Mol Biol 2007; 365:930-44. [PMID: 17098254 DOI: 10.1016/j.jmb.2006.10.037] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2006] [Revised: 09/28/2006] [Accepted: 10/09/2006] [Indexed: 01/27/2023]
Abstract
RNA molecules have an inherent flexibility that enables recognition of other interacting partners through potential disorder-order transitions, yet studies to quantify such motional dynamics remain few. With an increasing database of three-dimensional structures of biologically important RNA molecules, quantifying such motions becomes important to link structural deformations with function. One such system studied intensely is domain 5 (D5) from the self-splicing group II introns, which is at the heart of its catalytic machinery. We report the dynamics of a 36 nucleotide D5 from the Pylaiella littoralis group II intron in the presence and absence of magnesium ions, and at a range of temperatures (298K-318 K). Using high-resolution NMR experiments of heteronuclear nuclear Overhauser enhancement (NOE), spin-lattice (R(1)), and spin-spin (R(2)) (13)C relaxation rates, we determined the rotational diffusion tensor of D5 using the ROTDIF program modified for RNA dynamic analysis (ROTDIF_RNA). The D5 rotational diffusion tensor has an axial symmetric ratio (D(||)/D(perpendicular)) of 1.7+/-0.3, consistent with an estimated overall rotational correlation time of tau(m)=(2D(||)+4D(perpendicular))(-1) of 6.1(+/-0.3) ns at 298 K and 4.1(+/-0.2) ns at 318 K. The measured relaxation data were analyzed with the reduced spectral density mapping formalism using assumed values of the chemical shift anisotropy of the (13)C spins. Both the relaxation data and the values of the spectral density function reveal that the functional groups in D5 implicated in magnesium ion binding and catalysis (catalytic triad, internal bulge, and tetraloop regions) exhibit thermally induced motion on a wide variety of timescales. Because these motions parallel those observed in the intramolecular stem-loop of the U6 element within the spliceosome, we hypothesize that such extensive dynamic disorder likely facilitates D5 engaging both binding and catalytic regions of the ribozyme, and these may be a conserved feature of the catalytic machinery essential for catalysis.
Collapse
Affiliation(s)
- Nadukkudy V Eldho
- Department of Molecular Genetics and Center for Structural Biology, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio 44195, USA
| | | |
Collapse
|
26
|
Kormos BL, Baranger AM, Beveridge DL. A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations. J Struct Biol 2006; 157:500-13. [PMID: 17194603 PMCID: PMC1994251 DOI: 10.1016/j.jsb.2006.10.022] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2006] [Revised: 10/02/2006] [Accepted: 10/04/2006] [Indexed: 11/21/2022]
Abstract
Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.
Collapse
Affiliation(s)
- Bethany L Kormos
- Chemistry Department and Molecular Biophysics Program, Wesleyan University, 237 Church St., Middletown, CT 06459, USA.
| | | | | |
Collapse
|
27
|
Law MJ, Linde ME, Chambers EJ, Oubridge C, Katsamba PS, Nilsson L, Haworth IS, Laird-Offringa IA. The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA. Nucleic Acids Res 2006; 34:275-85. [PMID: 16407334 PMCID: PMC1326249 DOI: 10.1093/nar/gkj436] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.
Collapse
Affiliation(s)
- Michael J. Law
- Department of Biochemistry and Molecular Biology, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- Department of Surgery, Keck School of Medicine, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
| | - Michael E. Linde
- Department of Biochemistry and Molecular Biology, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- Department of Surgery, Keck School of Medicine, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
| | - Eric J. Chambers
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
| | - Chris Oubridge
- MRC Laboratory of Molecular BiologyHills Road, Cambridge CB2 2QH, UK
| | - Phinikoula S. Katsamba
- Department of Biochemistry and Molecular Biology, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- Department of Surgery, Keck School of Medicine, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences at NovumSE-141 57 Huddinge, Sweden
| | - Ian S. Haworth
- Department of Biochemistry and Molecular Biology, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
| | - Ite A. Laird-Offringa
- Department of Biochemistry and Molecular Biology, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- Department of Surgery, Keck School of Medicine, University of Southern CaliforniaLos Angeles, CA 90089-9176, USA
- To whom correspondence should be addressed. Tel: +1 323 865 0655; Fax: +1 323 865 0158;
| |
Collapse
|