RNA/Peptide Binding Driven by Electrostatics-Insight from Bidirectional Pulling Simulations

RNA/protein interactions play crucial roles in controlling gene expression. They are becoming important targets for pharmaceutical applications. Due to RNA flexibility and to the strength of electrostatic interactions, standard docking methods are insufficient. We here present a computational method which allows studying the binding of RNA molecules and charged peptides with atomistic, explicit-solvent molecular dynamics. In our method, a suitable estimate of the electrostatic interaction is used as an order parameter (collective variable) which is then accelerated using bidirectional pulling simulations. Since the electrostatic interaction is only used to enhance the sampling, the approximations used to compute it do not affect the final accuracy. The method is employed to characterize the binding of TAR RNA from HIV-1 and a small cyclic peptide. Our simulation protocol allows blindly predicting the binding pocket and pose as well as the binding affinity. The method is general and could be applied to study other electrostatics-driven binding events.

GFT projection NMR for efficient (1)H/ (13)C sugar spin system identification in nucleic acids

A newly implemented G-matrix Fourier transform (GFT) (4,3)D HC(C)CH experiment is presented in conjunction with (4,3)D HCCH to efficiently identify (1)H/(13)C sugar spin systems in (13)C labeled nucleic acids. This experiment enables rapid collection of highly resolved relay 4D HC(C)CH spectral information, that is, shift correlations of (13)C-(1)H groups separated by two carbon bonds. For RNA, (4,3)D HC(C)CH takes advantage of the comparably favorable 1'- and 3'-CH signal dispersion for complete spin system identification including 5'-CH. The (4,3)D HC(C)CH/HCCH based strategy is exemplified for the 30-nucleotide 3'-untranslated region of the pre-mRNA of human U1A protein.

Preparation of RNA samples with narrow line widths for solid state NMR investigations

Solid state NMR can provide detailed structural and dynamic information on biological systems that cannot be studied under solution conditions, and can investigate motions which occur with rates that cannot be fully studied by solution NMR. This approach has successfully been used to study proteins, but the application of multidimensional solid state NMR to RNA has been limited because reported line widths have been too broad to execute most multidimensional experiments successfully. A reliable method to generate spectra with narrow line widths is necessary to apply the full range of solid state NMR spectroscopic approaches to RNA. Using the HIV-1 transactivation response (TAR) RNA as a model, we present an approach based on precipitation with polyethylene glycol that improves the line width of (13)C signals in TAR from >6 ppm to about 1 ppm, making solid state 2D NMR studies of selectively enriched RNAs feasible at ambient temperature.

Aminoglycoside antibiotics bind to the influenza A virus RNA promoter

Aminoglycosides bind to the influenza A virus promoter (vRNA) at submicromolar concentration. The complex structure between the vRNA and neomycin illustrates that binding of neomycin causes a conformational change which would affect further transcription processes. Thus, aminoglycosides represent lead compounds for the discovery of antiviral therapeutics against influenza A virus.

Independent alignment of RNA for dynamic studies using residual dipolar couplings

Molecular motion and dynamics play an essential role in the biological function of many RNAs. An important source of information on biomolecular motion can be found in residual dipolar couplings which contain dynamics information over the entire ms-ps timescale. However, these methods are not fully applicable to RNA because nucleic acid molecules tend to align in a highly collinear manner in different alignment media. As a consequence, information on dynamics that can be obtained with this method is limited. In order to overcome this limitation, we have generated a chimeric RNA containing both the wild type TAR RNA, the target of our investigation of dynamics, as well as the binding site for U1A protein. When U1A protein was bound to the portion of the chimeric RNA containing its binding site, we obtained independent alignment of TAR by exploiting the physical chemical characteristics of this protein. This technique can allow the extraction of new information on RNA dynamics, which is particularly important for time scales not covered by relaxation methods where important RNA motions occur.

Counterion Redistribution upon Binding of a Tat-Protein Mimic to HIV-1 TAR RNA

Binding of proteins and small molecules to RNA involves many electrostatic interactions, which may alter the distribution of ions around the RNA molecule. Here, we use molecular dynamics simulations to investigate how binding of a cyclic peptide mimic of the HIV-1 Tat protein affects the ionic distribution around the HIV-1 TAR RNA element. The calculations reproduce the structural properties observed in NMR studies of TAR and its complex. They also provide insight into the rearrangement of counterions during the molecular recognition events leading to the formation of the protein/RNA complex.

The H/ACA RNP assembly factor SHQ1 functions as an RNA mimic

SHQ1 is an essential assembly factor for H/ACA ribonucleoproteins (RNPs) required for ribosome biogenesis, pre-mRNA splicing, and telomere maintenance. SHQ1 binds dyskerin/NAP57, the catalytic subunit of human H/ACA RNPs, and this interaction is modulated by mutations causing X-linked dyskeratosis congenita. We report the crystal structure of the C-terminal domain of yeast SHQ1, Shq1p, and its complex with yeast dyskerin/NAP57, Cbf5p, lacking its catalytic domain. The C-terminal domain of Shq1p interacts with the RNA-binding domain of Cbf5p and, through structural mimicry, uses the RNA-protein-binding sites to achieve a specific protein-protein interface. We propose that Shq1p operates as a Cbf5p chaperone during RNP assembly by acting as an RNA placeholder, thereby preventing Cbf5p from nonspecific RNA binding before association with an H/ACA RNA and the other core RNP proteins.

Theory of nonrigid rotational motion applied to NMR relaxation in RNA

Solution NMR spectroscopy can elucidate many features of the structure and dynamics of macromolecules, yet relaxation measurements, the most common source of experimental information on dynamics, can sample only certain ranges of dynamic rates. A complete characterization of motion of a macromolecule thus requires the introduction of complementary experimental approaches. Solid-state NMR spectroscopy successfully probes the time scale of nanoseconds to microseconds, a dynamic window where solution NMR results have been deficient, and probes conditions where the averaging effects of rotational diffusion of the molecule are absent. Combining the results of the two distinct techniques within a single framework provides greater insight into dynamics, but this task requires the common interpretation of results recorded under very different experimental conditions. Herein, we provide a unified description of dynamics that is robust to the presence of large-scale conformational exchange, where the diffusion tensor of the molecule varies on a time scale comparable to rotational diffusion in solution. We apply this methodology to the HIV-1 TAR RNA molecule, where conformational rearrangements are both substantial and functionally important. The formalism described herein is of greater generality than earlier combined solid-state/solution NMR interpretations, if detailed molecular structures are available, and can offer a more complete description of RNA dynamics than either solution or solid-state NMR spectroscopy alone.

Interactions of protein side chains with RNA defined with REDOR solid state NMR

Formation of the complex between human immunodeficiency virus type-1 Tat protein and the transactivation response region (TAR) RNA is vital for transcriptional elongation, yet the structure of the Tat-TAR complex remains to be established. The NMR structures of free TAR, and TAR bound to Tat-derived peptides have been obtained by solution NMR, but only a small number of intermolecular NOEs could be identified unambiguously, preventing the determination of a complete structure. Here we show that a combination of multiple solid state NMR REDOR experiments can be used to obtain multiple distance constraints from (15)N to (13)C spins within the backbone and side chain guanidinium groups of arginine in a Tat-derived peptide, using (19)F spins incorporated into the base of U23 in TAR and (31)P spins in the P22 and P23 phosphate groups. Distances between the side chain of Arg52 and the base and phosphodiester backbone near U23 measured by REDOR NMR are comparable to distances observed in solution NMR-derived structural models, indicating that interactions of TAR RNA with key amino acid side chains in Tat are the same in the amorphous solid state as in solution. This method is generally applicable to other protein-RNA complexes where crystallization or solution NMR has failed to provide high resolution structural information.

Comparative analysis of glutaredoxin domains from bacterial opportunistic pathogens

Glutaredoxin proteins (GLXRs) are essential components of the glutathione system that reductively detoxify substances such as arsenic and peroxides and are important in the synthesis of DNA via ribonucleotide reductases. NMR solution structures of glutaredoxin domains from two Gram-negative opportunistic pathogens, Brucella melitensis and Bartonella henselae, are presented. These domains lack the N-terminal helix that is frequently present in eukaryotic GLXRs. The conserved active-site cysteines adopt canonical proline/tyrosine-stabilized geometries. A difference in the angle of α-helix 2 relative to the β-sheet surface and the presence of an extended loop in the human sequence suggests potential regulatory regions and/or protein-protein interaction motifs. This observation is consistent with mutations in this region that suppress defects in GLXR-ribonucleotide reductase interactions. These differences between the human and bacterial forms are adjacent to the dithiol active site and may permit species-selective drug design.

NMR structure of an acyl-carrier protein from Borrelia burgdorferi

Nearly complete resonance assignment and the high-resolution NMR structure of the acyl-carrier protein from Borrelia burgdorferi, a target of the Seattle Structural Genomics Center for Infectious Disease (SSGCID) structure-determination pipeline, are reported. This protein was chosen as a potential target for drug-discovery efforts because of its involvement in fatty-acid biosynthesis, an essential metabolic pathway, in bacteria. It was possible to assign >98% of backbone resonances and >92% of side-chain resonances using multidimensional NMR spectroscopy. The NMR structure was determined to a backbone r.m.s.d. of 0.4 Å and contained four α-helices and two 3(10)-helices. A structure-homology search revealed that this protein is highly similar to the acyl-carrier protein from Aquifex aeolicus.

Structural genomics of infectious disease drug targets: the SSGCID

The Seattle Structural Genomics Center for Infectious Disease (SSGCID) is a consortium of researchers at Seattle BioMed, Emerald BioStructures, the University of Washington and Pacific Northwest National Laboratory that was established to apply structural genomics approaches to drug targets from infectious disease organisms. The SSGCID is currently funded over a five-year period by the National Institute of Allergy and Infectious Diseases (NIAID) to determine the three-dimensional structures of 400 proteins from a variety of Category A, B and C pathogens. Target selection engages the infectious disease research and drug-therapy communities to identify drug targets, essential enzymes, virulence factors and vaccine candidates of biomedical relevance to combat infectious diseases. The protein-expression systems, purified proteins, ligand screens and three-dimensional structures produced by SSGCID constitute a valuable resource for drug-discovery research, all of which is made freely available to the greater scientific community. This issue of Acta Crystallographica Section F, entirely devoted to the work of the SSGCID, covers the details of the high-throughput pipeline and presents a series of structures from a broad array of pathogenic organisms. Here, a background is provided on the structural genomics of infectious disease, the essential components of the SSGCID pipeline are discussed and a survey of progress to date is presented.

Examining the relationship between RNA function and motion using nuclear magnetic resonance

The biological function of proteins and nucleic acids relies on their complex structures, yet dynamics provides an additional layer of functional adaptability. Numerous studies have demonstrated that RNA is only able to perform the multitude of functions for which it is responsible by readily changing its conformation in response to binding of proteins or small molecules. Examination of RNA dynamics is therefore essential to understanding its biological function. Nuclear magnetic resonance (NMR) has emerged as a leading technique for the examination of RNA motion and conformational transitions. It can examine domain motions as well as motion with atomic level resolution over a wide range of time scales. This review examines how NMR spectroscopy can be applied to examine the relationship between function and dynamics in RNA.

γ-AApeptides bind to RNA by mimicking RNA-binding proteins

The interactions between proteins and RNAs are of vital importance for many cellular processes, including transcription and processing of RNA, translation, and viral infections. Here we report an γ-AApeptide that can mimic HIV-1 Tat protein and bind to TAR RNAs of HIV and BIV with nanomolar affinity, comparable to that of the RNA-binding fragment of Tat (amino acids 49-58). The interaction is resistant to the presence of a large excess of tRNA. With resistance to proteolytic hydrolysis and limitless potential for diversification, γ-AApeptides may emerge as a new class of peptidomimetics to modulate RNA-protein interactions.

Finding the missing code of RNA recognition by PUF proteins

Pumilio and FBF homology (PUF) proteins represent highly promising candidates for engineering sequence-specific RNA recognition, but were only known to recognize G, A, and U, significantly limiting applications. Two groups (Filipovska et al., 2011; Dong et al., 2011) have now reported the discovery of the cytosine-recognition code for PUF proteins.

Inhibition of both HIV-1 reverse transcription and gene expression by a cyclic peptide that binds the Tat-transactivating response element (TAR) RNA

The RNA response element TAR plays a critical role in HIV replication by providing a binding site for the recruitment of the viral transactivator protein Tat. Using a structure-guided approach, we have developed a series of conformationally-constrained cyclic peptides that act as structural mimics of the Tat RNA binding region and block Tat-TAR interactions at nanomolar concentrations in vitro. Here we show that these compounds block Tat-dependent transcription in cell-free systems and in cell-based reporter assays. The compounds are also cell permeable, have low toxicity, and inhibit replication of diverse HIV-1 strains, including both CXCR4-tropic and CCR5-tropic primary HIV-1 isolates of the divergent subtypes A, B, C, D and CRF01_AE. In human peripheral blood mononuclear cells, the cyclic peptidomimetic L50 exhibited an IC₅₀ ∼250 nM. Surprisingly, inhibition of LTR-driven HIV-1 transcription could not account for the full antiviral activity. Timed drug-addition experiments revealed that L-50 has a bi-phasic inhibition curve with the first phase occurring after HIV-1 entry into the host cell and during the initiation of HIV-1 reverse transcription. The second phase coincides with inhibition of HIV-1 transcription. Reconstituted reverse transcription assays confirm that HIV-1 (−) strand strong stop DNA synthesis is blocked by L50-TAR RNA interactions in-vitro. These findings are consistent with genetic evidence that TAR plays critical roles both during reverse transcription and during HIV gene expression. Our results suggest that antiviral drugs targeting TAR RNA might be highly effective due to a dual inhibitory mechanism.

The HIV-1 transactivator protein (Tat), together with the elongation factor P-TEFb binds to an HIV-1 RNA secondary structure in the 5′-UTRs of nascent viral mRNAs (TAR) and promotes transcription elongation. This process has been an attractive target for drug development but previous inhibitors that bind either Tat or TAR have been plagued by poor inhibition of virus replication, limited cell penetration, and off-target effects. In this article, we describe a series of rationally designed cyclic peptides that block Tat-TAR interactions. L50, the most potent of these compounds, inhibits a wide range of HIV-1 strains from around the world. Remarkably, L50 inhibits two distinct steps in the HIV-1 lifecycle. As expected, L50 inhibits Tat-dependent HIV-1 transcription, but the majority of its anti-HIV activity is due to a block in reverse transcription, i.e. synthesis of the proviral DNA from the RNA genome. L50 inhibition of reverse transcription reveals an important role for TAR RNA during reverse transcription as well as providing one of first examples of a drug with a dual mechanism of action.

13C/15N-19F intermolecular REDOR NMR study of the interaction of TAR RNA with Tat peptides

The complex of the HIV TAR RNA with the viral regulatory protein Tat is of considerable interest, but the plasticity of this interaction has made it impossible so far to establish the structure of that complex. In order to explore a new approach to obtain structural information on protein-RNA complexes, we performed (13)C/(15)N-(19)F REDOR NMR experiments in the solid state on TAR bound to a peptide comprising the RNA-binding section of Tat. A critical arginine in the peptide was uniformly (13)C and (15)N labeled, and 5-fluorouridine was incorporated at the U23 position of TAR. REDOR irradiation resulted in dephasing of the (13)C and (15)N resonances, indicating the proximity of the U23(5F)-C and U23(5F)-N spin pairs. Best fits to the REDOR data show the U23(5F)-C distances and the U23(5F)-N distances are in good agreement with the distances obtained from solution NMR structures of partial complexes of Tat with TAR. These results demonstrate that it is possible to study protein-RNA complexes using solid-state REDOR NMR measurements, adding to a growing list of solid state techniques for studying protein-nucleic acid complexes.

Slow exchange model of nonrigid rotational motion in RNA for combined solid-state and solution NMR studies

Functional RNA molecules are conformationally dynamic and sample a multitude of dynamic modes over a wide range of frequencies. Thus, a comprehensive description of RNA dynamics requires the inclusion of a broad range of motions across multiple dynamic rates which must be derived from multiple spectroscopies. Here we describe a slow conformational exchange theoretical approach to combining the description of local motions in RNA that occur in the nanosecond to microsecond window and are detected by solid-state NMR with nonrigid rotational motion of the HIV-1 transactivation response element (TAR) RNA in solution as observed by solution NMR. This theoretical model unifies the experimental results generated by solution and solid-state NMR and provides a comprehensive view of the dynamics of HIV-1 TAR RNA, a well-known paradigm of an RNA where function requires extensive conformational rearrangements. This methodology provides a quantitative atomic level view of the amplitudes and rates of the local and collective displacements of the TAR RNA molecule and provides directly motional parameters for the conformational capture hypothesis of this classical RNA-ligand interaction.

Essential structural requirements for specific recognition of HIV TAR RNA by peptide mimetics of Tat protein

The pharmacological disruption of the interaction between the HIV Tat protein and its cognate transactivation response RNA (TAR) would generate novel anti-viral drugs with a low susceptibility to drug resistance, but efforts to discover ligands with sufficient potency to warrant pharmaceutical development have been unsuccessful. We have previously described a family of structurally constrained β-hairpin peptides that potently inhibits viral growth in HIV-infected cells. The nuclear magnetic resonance (NMR) structure of an inhibitory complex revealed that the peptide makes intimate contacts with the 3-nt bulge and the upper helix of the RNA hairpin, but that a single residue contacts the apical loop where recruitment of the essential cellular co-factor cyclin T₁ occurs. Attempting to extend the peptide to form more interactions with the RNA loop, we examined a library of longer peptides and achieved >6-fold improvement in affinity. The structure of TAR bound to one of the extended peptides reveals that the peptide slides down the major groove of the RNA, relative to our design, in order to maintain critical interactions with TAR. These conserved contacts involve three amino acid side chains and identify critical interaction points required for potent and specific binding to TAR RNA. They constitute a template of essential interactions required for inhibition of this RNA.

Fragment-based discovery of novel thymidylate synthase leads by NMR screening and group epitope mapping

Solution-state nuclear magnetic resonance (NMR) is a versatile tool for the study of binding interactions between small molecules and macromolecular targets. We applied ligand-based NMR techniques to the study of human thymidylate synthase (hTS) using known nanomolar inhibitors and a library of small molecule fragments. Screening by NMR led to the rapid identification of ligand pairs that bind in proximal sites within the cofactor-binding pocket of hTS. Screening hits were used as search criteria within commercially available sources, and a subset of catalog analogs were tested for potency by in vitro assay and binding affinity by quantitative saturation transfer difference (STD)-NMR titration. Two compounds identified by this approach possess low micromolar affinity and potency, as well as excellent binding efficiency against hTS. Relative binding orientations for both leads were modeled using AutoDock, and the most likely bound conformations were validated using experimentally derived STD-NMR binding epitope data. These ligands represent novel starting points for fragment-based drug design of non-canonical TS inhibitors, and their binding epitopes highlight important and previously unexploited interactions with conserved residues in the cofactor-binding site.

Novel protein-protein contacts facilitate mRNA 3'-processing signal recognition by Rna15 and Hrp1

Precise 3'-end processing of mRNA is essential for correct gene expression, yet in yeast, 3'-processing signals consist of multiple ambiguous sequence elements. Two neighboring elements upstream of the cleavage site are particularly important for the accuracy (positioning element) and efficiency (efficiency element) of 3'-processing and are recognized by the RNA-binding proteins Rna15 and Hrp1, respectively. In vivo, these interactions are strengthened by the scaffolding protein Rna14 that stabilizes their association. The NMR structure of the 34 -kDa ternary complex of the RNA recognition motif (RRM) domains of Hrp1 and Rna15 bound to this pair of RNA elements was determined by residual dipolar coupling and paramagnetic relaxation experiments. It reveals how each of the proteins binds to RNA and introduces a novel class of protein-protein contact in regions of previously unknown function. These interdomain contacts had previously been overlooked in other multi-RRM structures, although a careful analysis suggests that they may be frequently present. Mutations in the regions of these contacts disrupt 3'-end processing, suggesting that they may structurally organize the ribonucleoprotein complexes responsible for RNA processing.

The Seattle Structural Genomics Center for Infectious Disease (SSGCID)

The NIAID-funded Seattle Structural Genomics Center for Infectious Disease (SSGCID) is a consortium established to apply structural genomics approaches to potential drug targets from NIAID priority organisms for biodefense and emerging and re-emerging diseases. The mission of the SSGCID is to determine approximately 400 protein structures over the next five years. In order to maximize biomedical impact, ligand-based drug-lead discovery campaigns will be pursued for a small number of high-impact targets. Here we review the center's target selection processes, which include pro-active engagement of the infectious disease research and drug therapy communities to identify drug targets, essential enzymes, virulence factors and vaccine candidates of biomedical relevance to combat infectious diseases. This is followed by a brief overview of the SSGCID structure determination pipeline and ligand screening methodology. Finally, specifics of our resources available to the scientific community are presented. Physical materials and data produced by SSGCID will be made available to the scientific community, with the aim that they will provide essential groundwork benefiting future research and drug discovery.

Intermediate rate atomic trajectories of RNA by solid-state NMR spectroscopy

Many RNAs undergo large conformational changes in response to the binding of proteins and small molecules. However, when RNA functional dynamics occur in the nanosecond-microsecond time scale, they become invisible to traditional solution NMR relaxation methods. Residual dipolar coupling methods have revealed the presence of extensive nanosecond-microsecond domain motions in HIV-1 TAR RNA, but this technique lacks information on the rates of motions. We have used solid-state deuterium NMR to quantitatively describe trajectories of key residues in TAR by exploiting the sensitivity of this technique to motions that occur in the nanosecond-microsecond regime. Deuterium line shape and relaxation data were used to model motions of residues within the TAR binding interface. The resulting motional models indicate two functionally essential bases within the single-stranded bulge sample both the free and Tat-bound conformations on the microsecond time scale in the complete absence of the protein. Thus, our results strongly support a conformational capture mechanism for recognition: the protein does not induce a new RNA structure, but instead captures an already-populated conformation.