RNA−Protein Intermolecular Recognition

mRNA interactions with disordered regions control protein activity

The cytoplasm is compartmentalized into different translation environments. mRNAs use their 3'UTRs to localize to distinct cytoplasmic compartments, including TIS granules (TGs). Many transcription factors, including MYC, are translated in TGs. It was shown that translation of proteins in TGs enables the formation of protein complexes that cannot be established when these proteins are translated in the cytosol, but the mechanism is poorly understood. Here we show that MYC protein complexes that involve binding to the intrinsically disordered region (IDR) of MYC are only formed when MYC is translated in TGs. TG-dependent protein complexes require TG-enriched mRNAs for assembly. These mRNAs bind to a new and widespread RNA-binding domain in neutral or negatively charged IDRs in several transcription factors, including MYC. RNA-IDR interaction changes the conformational ensemble of the IDR, enabling the formation of MYC protein complexes that act in the nucleus and control functions that cannot be accomplished by cytosolically-translated MYC. We propose that certain mRNAs have IDR chaperone activity as they control IDR conformations. In addition to post-translational modifications, we found a novel mode of protein activity regulation. Since RNA-IDR interactions are prevalent, we suggest that mRNA-dependent control of protein functional states is widespread.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

An efficient tRNA cleaver without additional co-reactants at physiological condition

A square-planar Cu(II) complex, [Cu(Me-Im)(gly-gly)]∙H₂O 1 (Me-Im = 1-methyl-imidazole, gly-gly = glycylglycinato), has been prepared and structurally characterized by single crystal X-ray. The complex 1 was tested for its ability to the transfer RNA by UV-vis spectroscopy, cyclic voltammetry (CV), capillary electrophoresis (CE). Comparative spectroscopic analysis shows a maximum fluorescence-quenching ratio of 0.41 of 1 upon binding to RNA, which gives a binding constant (K_b) of 1.24 × 10⁵ M^-1. Cyclic voltammograms of complex 1 attached on the mercaptoethanol (-SH) linked Au electrodes in phosphate buffer solution give a well-defined and quasi-reversible redox couple, indicate complex 1 can efficiently degrade the high-order structure of RNA in physiological conditions (pH 7.0 phosphate buffer solution at 37 °C) without additional co-reactants, yielding a digestion coefficient more than 90% within 113 h. This study targeting the genetic biomacromolecule degradation based on the strong binding of chemical nucleases paves an important way to the novel materials in the decontamination of microorganisms (e.g., bacteria and viruses) at mild condition.

Crystal structure of the ribonuclease-P-protein subunit from Staphylococcus aureus

Staphylococcus aureus ribonuclease-P-protein subunit (RnpA) is a promising antimicrobial target that is a key protein component for two essential cellular processes, RNA degradation and transfer-RNA (tRNA) maturation. The first crystal structure of RnpA from the pathogenic bacterial species, S. aureus, is reported at 2.0 Å resolution. The structure presented maintains key similarities with previously reported RnpA structures from bacteria and archaea, including the highly conserved RNR-box region and aromatic residues in the precursor-tRNA 5'-leader-binding domain. This structure will be instrumental in the pursuit of structure-based designed inhibitors targeting RnpA-mediated RNA processing as a novel therapeutic approach for treating S. aureus infections.

Determining Functional Aptamer-Protein Interaction by Biolayer Interferometry

Short single-stranded nucleic acids called aptamers are widely being explored as recognition molecules of high affinity and specificity for binding a wide range of target molecules, particularly protein targets. In biolayer interferometry (BLI), a simple Dip-and-Read approach in which the aptamer-coated biosensors are dipped into microplate wells is used to study the interactions between an aptamer and its target protein. Here we describe the protocol for the analysis of the interaction between a well-characterized anti-thrombin RNA aptamer with thrombin (Basic Protocol). We also report on the protocol for the affinity screening of a panel of anti-thrombin RNA aptamers with a single phosphorodithioate (PS2) modification, whereby the position of the modification along the RNA backbone is varied systematically (Support Protocol). The PS2 modification uses two sulfur atoms to replace two non-bridging oxygen atoms at an internucleotide phosphodiester backbone linkage. The PS2-modified RNAs are nuclease resistant and several in vitro and in vivo assays have demonstrated their biological activity. For example, combining the PS2 with the 2'-OMe modification affords increased loading of modified small interfering RNA (siRNA) duplexes into the RNA-induced silencing complex (RISC) as well as enhanced gene-silencing antitumor activity. © 2016 by John Wiley & Sons, Inc.

Evoking picomolar binding in RNA by a single phosphorodithioate linkage

RNA aptamers are synthetic oligonucleotide-based affinity molecules that utilize unique three-dimensional structures for their affinity and specificity to a target such as a protein. They hold the promise of numerous advantages over biologically produced antibodies; however, the binding affinity and specificity of RNA aptamers are often insufficient for successful implementation in diagnostic assays or as therapeutic agents. Strong binding affinity is important to improve the downstream applications. We report here the use of the phosphorodithioate (PS2) substitution on a single nucleotide of RNA aptamers to dramatically improve target binding affinity by ∼1000-fold (from nanomolar to picomolar). An X-ray co-crystal structure of the α-thrombin:PS2-aptamer complex reveals a localized induced-fit rearrangement of the PS2-containing nucleotide which leads to enhanced target interaction. High-level quantum mechanical calculations for model systems that mimic the PS2 moiety and phenylalanine demonstrate that an edge-on interaction between sulfur and the aromatic ring is quite favorable, and also confirm that the sulfur analogs are much more polarizable than the corresponding phosphates. This favorable interaction involving the sulfur atom is likely even more significant in the full aptamer-protein complexes than in the model systems.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Synthetic RNA recognition motifs that selectively recognize HIV-1 trans-activation response element hairpin RNA

A multitude of RNA hairpins are directly implicated in human disease. Many of these RNAs are potentially valuable targets for drug discovery and basic research. However, very little is known about the molecular requirements for achieving sequence-selective recognition of a particular RNA sequence and structure. Although a relatively modest number of synthetic small to medium-sized RNA-binding molecules have been reported, rapid identification of sequence-selective RNA-binding molecules remains a daunting challenge. RNA recognition motif (RRM) domains may represent unique privileged scaffolds for the generation of synthetic proteins that selectively recognize structured disease-relevant RNAs, including RNA hairpins. As a demonstration of this potential, we mutated putative RNA-binding regions within the U1A RRM and a variant thereof and screened these synthetic proteins for affinity to HIV-1 trans-activation response (TAR) element hairpin RNA. Some of these U1A-derived proteins bind TAR with single-digit micromolar dissociation constants, and they do so preferentially over the native protein's original target RNA (U1hpII) and a DNA TAR variant. Binding affinity is not appreciably diminished by addition of 10 molar equivalents of cellular tRNAs from Escherichia coli. Taken together, our findings represent the first synthetic RRMs that selectively bind a disease-relevant RNA hairpin and may represent a general approach for achieving sequence-selective recognition of RNA hairpins, which are the focus of therapeutic discovery and basic research.

Analysis of protein-RNA complexes involving a RNA recognition motif engineered to bind hairpins with seven- and eight-nucleotide loops

U1A binds U1hpII, a hairpin RNA with a 10-nucleotide loop. A U1A mutant (ΔK50ΔM51) binds U1hpII-derived hairpins with shorter loops, making it an interesting scaffold for engineering or evolving proteins that bind similarly sized disease-related hairpin RNAs. However, a more detailed understanding of complexes involving ΔK50ΔM51 is likely a prerequisite to generating such proteins. Toward this end, we measured mutational effects for complexes involving U1A ΔK50ΔM51 and U1hpII-derived hairpin RNAs with seven- or eight-nucleotide loops and identified contacts that are critical to the stabilization of these complexes. Our data provide valuable insight into sequence-selective recognition of seven- or eight-nucleotide loop hairpins by an engineered RNA binding protein.

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

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.

Extracellular functional noncoding nucleic acid bioaptamers and angiotropin RNP ribokines in vascularization and self-tolerance

Endogenous extracellular and circulating functional small noncoding nucleic acids (ncNAs; <200 nucleotides) and complexes with proteins (ribonucleoproteins; RNPs) make up varying biolibraries of molecular imprints of cellular histories. They are nascently formed upon cellular activation by extrinsic (environmental) factors, including mitogens, cell-mediated immune memory reactions (Landsteiner-Chase-Lawrence transfer factors), and metabolic (hypoxia) and (physical) shear stress forces. Those factors are conventional models for epigenetic (non-Mendelian) vascular remodeling variations directed rather to proteinaceous gene expression and regulation than genomic DNA sequence changes. Structurally defined ncNAs are described as small hairpin nc-shRNA bioaptamers in interaction with proteins forming functional (Cu,Ca,Na,K)-metalloregulated complexes (CuRNP; angiotropins). As nonmitogenic angiomorphogen cytokines (ribokines), they may reprogram confluent quiescent (contact-inhibited) endothelial cell types to migratory, phagokinetically active phenotypes in the morphogenesis of tolerated neovascular patterns. Their functions in organized and mess-chaotic vascular patterns were investigated with regard to master gene, information, epigenetic, vascular, and tumor factors. Some ncNAs feature three-dimensional codes (3D episcripts) for distinct protein conformer phases. They are suggested as being specific recognition types, the estimated repertoires of which are superior in diversity and specificity to conventional immune (glyco-)proteins. For episcription of phenotype variations, they may address and integrate information flow on molecular shapes to protein-mediated nucleic acid processing and [post-]translational modification mechanisms in ncNA-, redox, and metalloregulated conformation phase pathway-locked loops (CPLL). Several vascular and cancer epigenetic regulator proteins are shown to be entangled by sharing helix-nucleating structural (proteomic) domains for interaction with functional nc-shRNA, termed K/RxxxH (K/R3H, -xK/RxxxHx(7-9)h/xx(7-9)h/xx(5-20)K/Rx-). This would suggest a tolerated mess-chaotic tumor vascularization as a bioaptamer disorder in ncNA-switched proteinaceous genetic and epigenetic processes.

Biophysical analysis of nucleic acids

This overview unit provides a thorough overview of biophysical methods used for structure analysis, including X-ray diffraction, nuclear magnetic resonance, optical spectroscopy, theoretical and computational methods, and single-molecule methods. Advantages and disadvantages of the methods are compared.

Conserved nucleotides in an RNA essential for hepatitis B virus replication show distinct mobility patterns

The number of regulatory RNAs with identified non-canonical structures is increasing, and structural transitions often play a role in their biological function. This stimulates interest in internal motions of RNA, which can underlie structural transitions. Heteronuclear NMR relaxation measurements, which are commonly used to study internal motion, only report on local motions of few sites within the molecule. Here we have studied a 27-nt segment of the human hepatitis B virus (HBV) pregenomic RNA, which is essential for viral replication. We combined heteronuclear relaxation with the new off-resonance ROESY technique, which reports on internal motions of H,H contacts. Using off-resonance ROESY, we could for the first time detect motion of through-space H,H contacts, such as in intra-residue base-ribose contacts or inter-nucleotide contacts, both essential for NMR structure determination. Motions in non-canonical structure elements were found primarily on the sub-nanosecond timescale. Different patterns of mobility were observed among several mobile nucleotides. The most mobile nucleotides are highly conserved among different HBV strains, suggesting that their mobility patterns may be necessary for the RNA’s biological function.

A New Tool in Peptide Engineering: A Photoswitchable Stilbene-type β-Hairpin Mimetic

Peptide secondary structure mimetics are important tools in medicinal chemistry, as they provide analogues of endogenous peptides with new physicochemical and pharmacological properties. The development, synthesis, photochemical investigation, and conformational analysis of a stilbene-type beta-hairpin mimetic capable of light-triggered conformational changes have been achieved. In addition to standard spectroscopic techniques (nuclear Overhauser effects, amide temperature coefficients, circular dichroism spectroscopy), the applicability of self-diffusion measurements (longitudinal eddy current delay pulsed-field gradient spin echo (LED-PGSE) NMR technique) in conformational studies of oligopeptides is demonstrated. The title compound shows photoisomerization of the stilbene chromophore, resulting in a change in solution conformation between an unfolded structure and a folded beta-hairpin.

Understanding the Template Preorganization Step of an Artificial Arginine Receptor§

A biomimetic complex which mimics the arginine-phosphonate diester interaction of the arginine fork is investigated with respect to structure and energetics of stable configurations. Within this work, we provide knowledge on local minima of the isolated system obtained from first-principles calculations. Non-negligible solvation effects are studied in a microsolvation approach. The interactions which govern the structural patterns of molecular recognition in this tweezer-guest complex can be significantly modulated by the action of hydrogen bond accepting and donating solvent molecules, such as dimethyl sulfoxide or water, which were present in experimental investigations on this system. Different tweezer-guest structures are evaluated with respect to their temperature-dependent thermodynamical properties as products of the first association reaction step of the bisphosphonate tweezer template and the guanidinium moiety.

Investigating local conformations of double-stranded DNA by low-energy circular dichroism of pyrrolo-cytosine

Local base stacking and conformational mobility play a major role in the structure and function of nucleic acids. We have recently shown that the low-energy CD spectrum of 2-aminopurine (2-AP), i.e., the CD spectral region above 300 nm, can be used to monitor conformational changes in polynucleotides at or near mono- and dinucleotide 2-AP residues that replace adenine residues in DNA and RNA. Here, we extend this technique to pyrrolo-cytosine (PC), a fluorescent analogue of cytosine. The low-energy CD spectrum of a PC dinucleotide in dsDNA exhibits an exciton couplet with two bands of opposite sign centered at 350 nm. This signal is characteristic of base stacking between adjacent PC residues in a right helical conformation. Isolated PC nucleotide residues inserted into polynucleotide chains also display chirality that reflects the asymmetric environment of their sequence context. Thus, we show that the low-energy CD spectra of C(PC)A and A(PC)C sequences in dsDNA have opposite signs. It appears that the measurement of the low-energy CD spectra of PC residues will usefully complement 2-AP measurements by serving to characterize the local conformations and dynamics of nucleic acids near C residues and G-C base pairs.

The double-stranded-RNA-binding motif: interference and much more

RNA duplexes have been catapulted into the spotlight by the discovery of RNA interference and related phenomena. But double-stranded and highly structured RNAs have long been recognized as key players in cell processes ranging from RNA maturation and localization to the antiviral response in higher organisms. Penetrating insights into the metabolism and functions of such RNAs have come from the identification and study of proteins that contain the double-stranded-RNA-binding motif.

Molecular phylogenetics and functional evolution of major RNA recognition domains of recently cloned and characterized autoimmune RNA-binding particle

Heterogeneous nuclear ribonucleoproteins (hnRNPs) are spliceosomal macromolecular assemblages and thus actively participate in pre-mRNA metabolism. They are composed of evolutionarily conserved and tandemly repeated motifs, where both RNA-binding and protein-protein recognition occur to achieve cellular activities. By yet unknown mechanisms, these ribonucleoprotein (RNP) particles are targeted by autoantibodies and hence play significant role in a variety of human systemic autoimmune diseases. This feature makes them important prognostic markers in terms of molecular epidemiology and pathogenesis of autoimmunity. Since RNP domain is one of the most conserved and widespread scaffolds, evolutionary analyses of these RNA-binding domains can provide further clues on disease-specific epitope formation. The study presented herein represents a sequence comparison of RNA-recognition regions of recently cloned and characterized human hnRNP A3 with those of other relevant hnRNP A/B-type proteins. Their implications in human autoimmunity are particularly emphasized.

Cyclic aromatic oligoamides as highly selective receptors for the guanidinium ion

A class of six-residue, shape-persistent aromatic oligoamide macrocycles bind the guanidinium ion with very high selectivity.

Validation of automated docking programs for docking and database screening against RNA drug targets

The increasing awareness of the essential role of RNA in controlling viral replication and in bacterial protein synthesis emphasizes the potential of ribonucleoproteins as targets for developing new antibacterial and antiviral drugs. RNA forms well defined three-dimensional structures with clefts and binding pockets reminiscent of the active sites of proteins. Furthermore, it precedes proteins in the translation pathway; inhibiting the function of a single RNA molecule would result in inhibition of multiple proteins. Thus, small molecules that bind RNA specifically would combine the advantages of antisense and RNAi strategies with the much more favorable medicinal chemistry of small-molecule therapeutics. The discovery of small-molecule inhibitors of RNA with attractive pharmacological potential would be facilitated if we had available effective computational tools of structure-based drug design. Here, we systematically test automated docking tools developed for proteins using existing three-dimensional structures of RNA-small molecule complexes. The results show that the native structures can generally be reproduced to within 2.5 angstroms more than 50-60% of the time. For more than half of the test complexes, the native ligand ranked among the top 10% compounds in a database-scoring test. Through this work, we provide parameters for the validated application of automated docking tools to the discovery of new inhibitors of RNA function.

Extracellular and circulating redox- and metalloregulated eRNA and eRNP: copper ion-structured RNA cytokines (angiotropin ribokines) and bioaptamer targets imparting RNA chaperone and novel biofunctions to S100-EF-hand and disease-associated proteins

Bioassays for cellular differentiation and tissue morphogenesis were used to design methods for isolation of bioactive redox- and metalloregulated nucleic acids and copper ion complexes with proteins from extracellular, circulating, wound, and supernatant fluids of cultured cells. In extracellular biospheres, diversities of nucleic acids were found to be secreted by cells upon activation. They may reflect nucleic acid biolibraries with molecular imprints of cellular history. After removal of protein components, eRNA prototypes exuded by activated cells were sequenced. They are small, endogenous, highly modified and edited, redox- and metalloregulated 5'-end phosphorylated extracellular eRNA (approximately 2-200 bases) with cellular, enzymic, and bioaptamer functions. Fenton-type OH* radical redox reactions may form modified nucleotides in RNA as wobbles eRNA per se, or as copper ion-complex with protein (e.g., S100A12-EF-hand protein, angiotropin-related protein, calgranulin-C, hippocampal neurite differentiation factor) are shown to be bioactive in vivo and in vitro as cytokines (ribokines) and as nonmitogenic angiomorphogens for endothelial cell differentiation in the formation of organoid supracellular capillary structures. As bioaptamers, copper ion-structured eRNA imparts novel biofunctions to proteins that they do not have on their own. The origin of extracellular RNA and intermediate precursors (up to 500 bases) was traced to intracellular parent nucleic acids. Intermediate precursors with and without partial homology were found. This suggests that bioaptamers are not directly retranslatable gene products. Metalloregulated eRNA bioaptamer function was investigated by domains (e.g. 5'...CUG...3' hairpin loop) for folding, bioactivity, and binding of protein with copper, calcium, and alkali metal ion affinity. Vice versa, metalloregulated nucleic acid-binding domains (K3H, R3H) in proteins were identified. Interaction of protein and eRNA docking potentials were visualized by 3D-rapid prototyping of accurate molecular image models based on crystallographic or NMR data. For S100A12-homologous proteins, receptor- and metalloregulated RNA chaperone-shaped protein assemblies were investigated. They suggest insight into signaling cascades as to how eRNA transmits its cytokine (ribokine) bioinformation from the extracellular RNA biosphere into cells. Proteomics of the extracellular RNA biosphere demonstrate the presence of nucleic acid-binding domain homologies in defense-, aging-, and disease-associated neuronal and other proteins as targets for RNA orphans. By structural relationships found to transmissible processes, proteinaceous transfer ("infectivity") and feedback of bioinformation beyond the central dogma of molecular biology are considered in terms of metalloregulated RNA bioaptamer function, nucleic acid-binding domains, and protein conformation.

Kinetic characterization of TAR RNA-Tat peptide and neomycin interactions by acoustic wave biosensor

The kinetics of binding of short Tat peptides and an aminoglycoside molecule to the human immunodeficiency virus-type 1(HIV-1) TAR RNA and to a bulge mutant analogue (MTAR) is studied in a biosensor format by monitoring the time course of the response in a series resonance frequency, using an acoustic wave biosensor. Association and dissociation rate constants are evaluated by fitting the experimental data to a simple 1:1 (Langmuir) model. Kinetic rate and equilibrium dissociation constants show that MTAR-peptide complexes are subject to a higher dissociation rate and are less stable compared to the corresponding TAR-peptide complexes. In addition, longer peptides display enhanced discrimination ability than a shorter peptide according to the equilibrium dissociation constants evaluated using this technique. K(D) values for TAR-Tat vs. MTAR-Tat complexes are 2.6 vs. 3.8 microM for Tat-12, 0.87 vs. 4.3 microM for Tat-18 and 0.93 vs. 1.6 microM for Tat-20. The equilibrium dissociation constant for TAR-neomycin complex is 12.4 microM and it is comparable to the values obtained from non-biosensor type assays. These findings are in parallel with those cited in the literature and the results from this study underline the potential of the acoustic wave sensor for detailed biophysical analysis of nucleic acid-ligand binding.

Computational analysis of hydrogen bonds in protein-RNA complexes for interaction patterns

Structural analysis of protein-RNA complexes is labor-intensive, yet provides insight into the interaction patterns between a protein and RNA. As the number of protein-RNA complex structures reported has increased substantially in the last few years, a systematic method is required for automatically identifying interaction patterns. This paper presents a computational analysis of the hydrogen bonds in the most representative set of protein-RNA complexes. The analysis revealed several interesting interaction patterns. (1) While residues in the beta-sheets favored unpaired nucleotides, residues in the helices showed no preference and residues in turns favored paired nucleotides. (2) The backbone hydrogen bonds were more dominant than the base hydrogen bonds in the paired nucleotides, but the reverse was observed in the unpaired nucleotides. (3) The protein-RNA complexes contained more paired nucleotides than unpaired nucleotides, but the unpaired nucleotides were observed more frequently interacting with the proteins. And (4) Arg-U, Thr-A, Lys-A, and Asn-U were the most frequently observed pairs. The interaction patterns discovered from the analysis will provide us with useful information in predicting the structure of the RNA binding protein and the structure of the protein binding RNA.

Probing of arginine residues in peptides and proteins using selective tagging and electrospray ionization mass spectrometry

A general labelling method is presented which allows the determination of the number of guanidine groups (related to arginine and homoarginine in peptides and proteins) by means of mass spectrometry. It implies a guanidine-selective derivatization step with 2,3-butanedione and an arylboronic acid under aqueous, alkaline conditions (pH 8-10). The reaction mixture is then directly analysed by electrospray ionization mass spectrometry without further sample pretreatment. Other amino acids are not affected by this reaction although it is demonstrated that lysine side-chains may be unambiguously identified when they are converted to homoarginine prior to derivatization. Guanidine functionalities, as e.g. in the amino acid arginine, are easily identified by the characteristic mass shift between underivatized and derivatized analyte. The tagging procedure is straightforward and selective for guanidine groups. The influence of several experimental parameters, especially the pH of the solution and the choice of reagents, is examined and the method is applied to various arginine-containing peptides and to lysozyme as a representative protein. Possible applications of this technique and its limitations are discussed.

The solution structure of YbcJ from Escherichia coli reveals a recently discovered alphaL motif involved in RNA binding

The structure of the recombinant Escherichia coli protein YbcJ, a representative of a conserved family of bacterial proteins (COG2501), was determined by nuclear magnetic resonance. The fold of YbcJ identified it as a member of the larger family of S4-like RNA binding domains. These domains bind to structured RNA, such as that found in tRNA, rRNA, and a pseudoknot of mRNA. The structure of YbcJ revealed a highly conserved patch of basic residues, comprising amino acids K26, K38, R55, K56, and K59, which likely participate in RNA binding.

Molecular dynamics simulation of the RNA complex of a double-stranded RNA-binding domain reveals dynamic features of the intermolecular interface and its hydration

The interaction between double-stranded RNA (dsRNA) and the third double-stranded domain (dsRBD) from Drosophila Staufen protein represents a paradigm to understand how the dsRBD protein family, one of the most common RNA-binding protein units, binds dsRNA. The nuclear magnetic resonance (NMR) structure of this complex and the x-ray structure of another family member revealed the stereochemical basis for recognition, but also raised new questions. Although the crystallographic studies revealed a highly ordered interface containing numerous water-mediated contacts, NMR suggested extensive residual motion at the interface. To address how interfacial motion contributes to molecular recognition in the dsRBD-dsRNA system, we conducted a 2-ns molecular dynamics simulation of the complex derived from Staufen protein and of the separate protein and RNA components. The results support the observation that a high degree of conformational flexibility is retained upon complex formation and that this involves interfacial residues that are critical for dsRBD-dsRNA binding. The structural origin of this residual flexibility is revealed by the analysis of the trajectory of motion. Individual basic side chains switch continuously from one RNA polar group to another with a residence time seldom exceeding 100 ps, while retaining favorable interaction with RNA throughout much of the simulation. Short-lived water molecules mediate some of these interactions for a large fraction of the trajectory studied here. This result indicates that water molecules are not statically associated with the interface, but continuously exchange with the bulk solvent on a 1-10-ps time scale. This work provides new insight into dsRBD-dsRNA recognition and builds upon a growing body of evidence, suggesting that short-lived dynamic interactions play important roles in protein-nucleic acid interactions.

Role of pyridine hydrogen-bonding sites in recognition of basic amino acid side chains

A series of three, new artificial receptors for guanidinium and ammonium guests has been synthesized. All three receptors have highly preorganized clefts bearing two carboxylate groups. They differ in the number of nitrogen atoms contained in their clefts, as follows: four N atoms in receptor 3, three N atoms in 4, and two nitrogens in 5. Crystallographic studies have produced the solid-state structures of the following guanidinium complexes of each receptor: 3.2CH(3)CH(2)NHC(NH(2))(2)(+), 4.2CH(3)NHC(NH(2))(2)(+), and 5.2C(NH(2))(3)(+). The conformations of the receptor molecules in all three complexes are very similar. N-Alkylguanidinium guests are bound in the clefts of 3 and 4 in similar manners, despite the loss of one hydrogen-bond acceptor nitrogen in 4 and the possible hindrance of the cavity by a CH group. In the guanidinium complex of 5, neither guest enters the cavity containing two CH groups. Complexation studies were conducted in methanol by (1)H NMR titration for several guanidinium and ammonium guests, including derivatives of the amino acids arginine and lysine. Receptor 5 binds all such guests weakly (K(s) < 4000), while 3 binds most guests very strongly (K(s) > 100 000). Receptor 3 is selective for arginine versus lysine, while 4 binds lysine better than does 3. The results generally underscore the importance of receptor preorganization and hydrogen-bonding complementarity in the design of receptors that can serve as probes for biomolecules.

Inhibition of the U1A–RNA complex by an aminoacridine derivative

The RNA recognition motif (RRM) is one of the most common RNA binding domains. There have been few investigations of small molecule inhibitors of RRM-RNA complexes, although these inhibitors could be valuable tools for probing biological processes involving RRM-RNA complexes and would have the potential to be effective drugs. In this paper, the inhibition by small molecules of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA has been investigated. An aminoacridine derivative has been found to promote dissociation of the U1A-stem loop 2 RNA complex with an IC(50) value of 1 microM. Fluorescence experiments indicate that two aminoacridine ligands bind to each RNA target site. RNase A footprinting suggests that one binding site may be near the base pair that closes the loop and the other may be in a more flexible region of the loop. The addition of the aminoacridine derivative to stem loop 2 RNA increases the susceptibility of other portions of the loop to digestion by RNase A, which implies that binding of the ligand changes the conformation or dynamics of the stem loop target site. Either direct binding to the RNA or indirect alteration of the structure or dynamics of the loop would be likely to inhibit binding of the U1A protein to this RNA.

How to detect internal motion by homonuclear NMR

NOESY and ROESY cross-peak intensities depend on internuclear distances and internal motion. Internal motion is usually ignored, and NOESY cross-peak intensities are interpreted in terms of internuclear distances only. Off-resonance ROESY experiments measure a weighted average of NOE and ROE. The weight can be described and experimentally set by an angle theta;. For large enough molecules, NOE and ROE have opposite signs. Therefore, each cross-peak intensity becomes zero for an angle theta;(0). For any sample, the maximum angle theta;(0) is determined by the overall motion of the molecule. Smaller theta;(0) values reflect the angular component of internal motions. Because individual cross-peaks are analyzed, the method evaluates internal motions of individual H,H vectors. The reduction of theta;(0) is largest for internal motions on a time scale of 100-300 ps. The sensitivity of theta;(0) for internal motions decreases with increasing molecular weight. We estimate that detecting internal motions will be practicable for molecules up to about 15 kDa. We describe a protocol to measure theta;(0) from a series of off-resonance ROESY spectra. For such a series, we describe the choice of experimental parameters, a procedure to extract theta;(0) from the raw data, and the interpretation of theta;(0) in terms of internal motions. In the small protein BPTI, we analyzed 75 cross-peaks. The precision of theta;(0) was 0.25 degrees, as compared to typical reductions of theta;(0) of 3 degrees. We found a well-defined maximum theta;(0) for cross-peaks in rigid parts of the molecule, which reflects the overall motion of the molecule. For BPTI, also many structurally important long-range cross-peaks appear rigid. The lower theta;(0) values of long-range contacts involving methyl groups are consistent with methyl rotation on the 25-ps time scale. The lower theta;(0) values of the flexible C-terminus and of flexible side chains translate into upper limits for the angular order parameter of 0.4 and 0.5-0.8, respectively. Off-resonance ROESY can monitor internal motions of H,H contacts that are used in a structure calculation. Because no isotope labeling is needed, off-resonance ROESY can be used to detect internal motions in a wide range of natural products.

Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions

Computer modeling has been developed and widely applied in studying molecules of biological interest. The force field is the cornerstone of computer simulations, and many force fields have been developed and successfully applied in these simulations. Two interesting areas are (a) studying enzyme catalytic mechanisms using a combination of quantum mechanics and molecular mechanics, and (b) studying macromolecular dynamics and interactions using molecular dynamics (MD) and free energy (FE) calculation methods. Enzyme catalysis involves forming and breaking of covalent bonds and requires the use of quantum mechanics. Noncovalent interactions appear ubiquitously in biology, but here we confine ourselves to review only noncovalent interactions between protein and protein, protein and ligand, and protein and nucleic acids.

Recent advances in RNA-protein recognition

The past few years have witnessed remarkable progress in knowledge of the structure and function of RNA-binding proteins and their RNA complexes. X-ray crystallography and NMR spectroscopy have provided structures for all major classes of RNA-binding proteins, both alone and complexed with RNA. New computational and experimental tools have provided unprecedented insight into the molecular basis of RNA recognition.

Binding of a cyclic BIV beta-Tat peptide with its TAR RNA construct

The ability of RNA structures to adopt diverse yet complex tertiary structures has resulted in numerous fascinating RNA-protein recognition events. It was recently reported that a close relative of the HIV Rev peptide, namely a 17 residue Tat peptide from bovine immuno-deficiency virus (BIV), is able to bind to the 28 nucleotide BIV TAR RNA construct. Here we report that by simply converting the 17 residue beta-ribbon peptide structure to a 19 residue cyclopeptide, the binding affinity (Kd) of the resulting cyclopeptide to the TAR RNA target, observed by fluorescence binding study, was enhanced approximately 5-fold.

Recognition between flexible protein molecules: induced and assisted folding

This review focuses on a very important but little understood type of molecular recognition--the recognition between highly flexible molecular structures. The formation of a specific complex in this case is a dynamic process that can occur through sequential steps of mutual conformational adaptation. This allows modulation of specificity and affinity of interaction in extremely broad ranges. The interacting partners can interact together to form a complex with entirely new properties and produce conformational signal transduction at substantial distance. We show that this type of recognition is frequent in formation of different protein-protein and protein-nucleic acid complexes. It is also characteristic for self-assembly of protein molecules from their unfolded fragments as well as for interaction of molecular chaperones with their substrates and it can be the origin of 'protein misfolding' diseases. Thermodynamic and kinetic features of this type of dynamic recognition and the principles underlying their modeling and analysis are discussed.

RNA recognition by a Staufen double-stranded RNA-binding domain

The double-stranded RNA-binding domain (dsRBD) is a common RNA-binding motif found in many proteins involved in RNA maturation and localization. To determine how this domain recognizes RNA, we have studied the third dsRBD from Drosophila Staufen. The domain binds optimally to RNA stem-loops containing 12 uninterrupted base pairs, and we have identified the amino acids required for this interaction. By mutating these residues in a staufen transgene, we show that the RNA-binding activity of dsRBD3 is required in vivo for Staufen-dependent localization of bicoid and oskar mRNAs. Using high-resolution NMR, we have determined the structure of the complex between dsRBD3 and an RNA stem-loop. The dsRBD recognizes the shape of A-form dsRNA through interactions between conserved residues within loop 2 and the minor groove, and between loop 4 and the phosphodiester backbone across the adjacent major groove. In addition, helix alpha1 interacts with the single-stranded loop that caps the RNA helix. Interactions between helix alpha1 and single-stranded RNA may be important determinants of the specificity of dsRBD proteins.

Structure of Hsp15 reveals a novel RNA-binding motif

We have solved the crystal structure of the heat shock protein Hsp15, a newly isolated and very highly inducible heat shock protein that binds the ribosome. Comparison of its structure with those of two RNA-binding proteins, ribosomal protein S4 and threonyl-tRNA synthetase, reveals a novel RNA-binding motif. This newly recognized motif is remarkably common, present in at least eight different protein families that bind RNA. The motif's surface is populated by conserved, charged residues that define a likely RNA-binding site. An intriguing pattern emerges: stress proteins, ribosomal proteins and tRNA synthetases repeatedly share a conserved motif. This may imply a hitherto unrecognized functional similarity between these three protein classes.

The structure of ribonuclease P protein from Staphylococcus aureus reveals a unique binding site for single-stranded RNA

Ribonuclease P (RNaseP) catalyses the removal of the 5'-leader sequence from pre-tRNA to produce the mature 5' terminus. The prokaryotic RNaseP holoenzyme consists of a catalytic RNA component and a protein subunit (RNaseP protein), which plays an auxiliary but essential role in vivo by binding to the 5'-leader sequence and broadening the substrate specificity of the ribozyme. We determined the three-dimensional high-resolution structure of the RNaseP protein from Staphylococcus aureus (117 amino acid residues) by nuclear magnetic resonance (NMR) spectroscopy in solution. The protein has an alphabeta-fold, similar to the ribonucleoprotein domain. We used small nucleic acid molecules as a model for the 5'-leader sequence to probe the propensity for generic single-stranded RNA binding on the protein surface. The NMR results reveal a contiguous interaction site, which is identical with the previously identified leader sequence binding site in RNaseP holoenzyme. The conserved arginine-rich motif does not bind single-stranded RNA. It is likely that this peptide segment binds selectively to double-stranded sections of P RNA, which are conformationally more rigid. Given the essentiality of RNaseP for the viability of the organism, knowledge of the S. aureus protein structure and insight into its interaction with RNA will help us to develop RNaseP and RNaseP protein as targets for novel antibiotics against this pathogen.