Evolvability of the mode of peptide binding by an RNA

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

A Review: The Antiviral Activity of Cyclic Peptides

In the design and development of therapeutic agents, macromolecules with restricted structures have stronger competitive edges than linear biological entities since cyclization can overcome the limitations of linear structures. The common issues of linear peptides include susceptibility to degradation of the peptidase enzyme, off-target effects, and necessity of routine dosing, leading to instability and ineffectiveness. The unique conformational constraint of cyclic peptides provides a larger surface area to interact with the target at the same time, improving the membrane permeability and in vivo stability compared to their linear counterparts. Currently, cyclic peptides have been reported to possess various activities, such as antifungal, antiviral and antimicrobial activities. To date, there is emerging interest in cyclic peptide therapeutics, and increasing numbers of clinically approved cyclic peptide drugs are available on the market. In this review, the medical significance of cyclic peptides in the defence against viral infections will be highlighted. Except for chikungunya virus, which lacks specific antiviral treatment, all the viral diseases targeted in this review are those with effective treatments yet with certain limitations to date. Thus, strategies and approaches to optimise the antiviral effect of cyclic peptides will be discussed along with their respective outcomes. Apart from isolated naturally occurring cyclic peptides, chemically synthesized or modified cyclic peptides with antiviral activities targeting coronavirus, herpes simplex viruses, human immunodeficiency virus, Ebola virus, influenza virus, dengue virus, five main hepatitis viruses, termed as type A, B, C, D and E and chikungunya virus will be reviewed herein.

Graphical Abstract

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Sequence and Functional Variation in the HIV-1 Rev Regulatory Axis

BACKGROUND

To complete its replication cycle, HIV-1 requires the nucleocytoplasmic export of intron-containing viral mRNAs. This process is ordinarily restricted by the cell, but HIV overcomes the block by means of a viral protein, Rev, and an RNA secondary structure found in all unspliced and incompletely spliced viral mRNAs called the Rev Response Element (RRE). In vivo activity of the Rev-RRE axis requires Rev binding to the RRE, oligomerization of Rev to form a competent ribonucleoprotein complex, and recruitment of cellular factors including Crm1 and RanGTP in order to export the targeted transcript. Sequence variability is observed among primary isolates in both Rev and the RRE, and the activity of both can be modulated through relatively small sequence changes. Primary isolates show differences in Rev-RRE activity and a few studies have found a correlation between lower Rev-RRE activity and slower progression of clinical disease. Lower Rev-RRE activity has also been associated with the evasion of cytotoxic T lymphocyte mediated killing.

CONCLUSION

The HIV-1 Rev-RRE regulatory axis is an understudied mechanism by which viral adaptation to diverse immune milieus may take place. There is evidence that this adaptation plays a role in HIV pathogenesis, particularly in immune evasion and latency, but further studies with larger sample sizes are warranted.

Native mass spectrometry reveals the initial binding events of HIV-1 rev to RRE stem II RNA

Nuclear export complexes composed of rev response element (RRE) ribonucleic acid (RNA) and multiple molecules of rev protein are promising targets for the development of therapeutic strategies against human immunodeficiency virus type 1 (HIV-1), but their assembly remains poorly understood. Using native mass spectrometry, we show here that rev initially binds to the upper stem of RRE IIB, from where it is relayed to binding sites that allow for rev dimerization. The newly discovered binding region implies initial rev recognition by nucleotides that are not part of the internal loop of RRE stem IIB RNA, which was previously identified as the preferred binding region. Our study highlights the unique capability of native mass spectrometry to separately study the binding interfaces of RNA/protein complexes of different stoichiometry, and provides a detailed understanding of the mechanism of RRE/rev association with implications for the rational design of potential drugs against HIV-1 infection.

Cyclic Peptide Inhibitors of the Tsg101 UEV Protein Interactions Refined through Global Docking and Gaussian Accelerated Molecular Dynamics Simulations

Tsg101 UEV domain proteins are potential targets for virus infection therapy, especially for HIV and Ebola viruses. Peptides are key in curbing virus transmission, and cyclic peptides have a greater survival time than their linear peptides. To date, the accurate prediction of cyclic peptide-protein receptors binding conformations still is challenging because of high peptide flexibility. Here, a useful approach combined the global peptide docking, Gaussian accelerated molecular dynamics (GaMD), two-dimensional (2D) potential of mean force (PMF), normal molecular dynamics (cMD), and solvated interaction energy (SIE) techniques. Then we used this approach to investigate the binding conformations of UEV domain proteins with three cyclic peptides inhibitors. We reported the possible cyclic peptide-UEV domain protein binding conformations via 2D PMF free energy profiles and SIE free energy calculations. The residues Trp145, Tyr147, and Trp148 of the native cyclic peptide (CP1) indeed play essential roles in the cyclic peptides-UEV domain proteins interactions. Our findings might increase the accuracy of cyclic peptide-protein conformational prediction, which may facilitate cyclic peptide inhibitor design. Our approach is expected to further aid in addressing the challenges in cyclic peptide inhibitor design.

Computational modeling of cyclic peptide inhibitor-MDM2/MDMX binding through global docking and Gaussian accelerated molecular dynamics simulations

MDM2 and MDMX are potential targets for p53-dependent cancer therapy. Peptides are key in cellular immunology and oncology, and cyclic peptides generally have higher half-life than their linear counterparts. However, prediction of cyclic peptide-protein binding is challenging with normal molecular simulation approaches because of high peptide flexibility. Here, we used global peptide docking, normal molecular dynamics, Gaussian accelerated molecular dynamics (GaMD), two-dimensional (2D) potential of mean force (PMF) profiles, and solvated interaction energy (SIE) techniques to investigate the interactions of MDM2/MDMX with three N-to-C-terminal cyclic peptide-based inhibitors. We determined the possible cyclic peptide-MDM2/MDMX complex structures via 2D PMF profiles and SIE calculations. Our findings increase the accuracy of peptide-protein structural prediction, which may facilitate cyclic peptide drug design. Advancements in the computational methods and computing power may further aid in addressing the challenges in cyclic peptide drug design. Communicated by Ramaswamy H. Sarma.

The new (dis)order in RNA regulation

RNA-binding proteins play a key role in the regulation of all aspects of RNA metabolism, from the synthesis of RNA to its decay. Protein-RNA interactions have been thought to be mostly mediated by canonical RNA-binding domains that form stable secondary and tertiary structures. However, a number of pioneering studies over the past decades, together with recent proteome-wide data, have challenged this view, revealing surprising roles for intrinsically disordered protein regions in RNA binding. Here, we discuss how disordered protein regions can mediate protein-RNA interactions, conceptually grouping these regions into RS-rich, RG-rich, and other basic sequences, that can mediate both specific and non-specific interactions with RNA. Disordered regions can also influence RNA metabolism through protein aggregation and hydrogel formation. Importantly, protein-RNA interactions mediated by disordered regions can influence nearly all aspects of co- and post-transcriptional RNA processes and, consequently, their disruption can cause disease. Despite growing interest in disordered protein regions and their roles in RNA biology, their mechanisms of binding, regulation, and physiological consequences remain poorly understood. In the coming years, the study of these unorthodox interactions will yield important insights into RNA regulation in cellular homeostasis and disease.

Diverse mutants of HIV RRE IIB recognize wild-type Rev ARM or Rev ARM R35G-N40V

The binding of human immunodeficiency virus Rev protein via its arginine-rich motif (ARM) to an internal loop in the Rev-response element region IIB (RRE IIB) is necessary for viral replication. Many variant RNAs and ARMs that bind Rev and RRE IIB have been found. Despite the essential role of Rev asparagine 40 in recognition, the Rev ARM double-mutant R35G-N40V functions well in a Rev-RRE IIB reporter assay, indicating R35G-N40V uses a distinct recognition strategy. To examine how RRE IIB may evolve specificity to wild-type Rev ARM and R35G-N40V, 10 RRE IIB libraries, each completely randomized in overlapping regions, were screened with wild-type Rev ARM and R35G-N40V using a reporter system based on bacteriophage λ N antitermination. Consistent with previous studies, a core element of RRE IIB did not vary, and substitutions occurred at conserved residues only in the presence of other substitutions. Notably, the groove-widening, non-canonical base-pair G48:G71 was mutable to U48:G71 without strong loss of binding to wild-type Rev ARM, suggesting U48:G71 performs the same role by adopting the nearly isosteric, reverse wobble base pair. Originating from RRE IIB, as few as one or two substitutions are sufficient to confer specificity to wild-type Rev or Rev R35G-N40. The diversity of RRE IIB mutants that maintain binding to wild-type Rev ARM and R35G-N40V supports neutral theories of evolution and illustrates paths by which viral RNA-protein interactions can evolve new specificities. Rev-RRE offers an excellent model with which to study the fine structure of how specificity evolves.

Randomized codon mutagenesis reveals that the HIV Rev arginine-rich motif is robust to substitutions and that double substitution of two critical residues alters specificity

The binding of the arginine-rich motif (ARM) of HIV Rev protein to its high-affinity site in stem IIB in the Rev response element (RRE) initiates assembly of a ribonucleoprotein complex that mediates the export of essential, incompletely spliced viral transcripts. Many biochemical, genetic, and structural studies of Rev-RRE IIB have been published, yet the roles of many peptide residues in Rev ARM are unconfirmed by mutagenesis. Rev aptamer I (RAI) is an optimized RRE IIB that binds Rev with higher affinity and for which mutational data are sparse. Randomized-codon libraries of Rev ARM were assayed for their ability to bind RRE IIB and RAI using a bacterial reporter system based on bacteriophage λ N-nut antitermination. Most Rev ARM residues tolerated substitutions without strong loss of binding to RRE IIB, and all except arginine 39 tolerated substitution without strong loss of binding to RAI. The pattern of critical Rev residues is not the same for RRE IIB and RAI, suggesting important differences between the interactions. The results support and aid the interpretation of existing structural models. Observed clinical variation is consistent with additional constraints on Rev mutation. By chance, we found double mutants of two highly critical residues, arginine 35 (to glycine) and asparagine 40 (to valine or lysine), that bind RRE IIB well, but not RAI. That an apparently distinct binding mode occurs with only two mutations highlights the ability of ARMs to evolve new recognition strategies and supports the application of neutral theories of evolution to protein-RNA recognition.

Amino acid requirement for the high affinity binding of a selected arginine-rich peptide with the HIV Rev-response element RNA

The arginine-rich motif is a class of short arginine-rich peptides that bind to specific RNA structures that has been found to be a versatile framework for the design and selection of RNA-binding peptides. We previously identified novel peptides that bind to the Rev-response element (RRE) RNA of the HIV from an arginine-rich polypeptide library (ARPL) consisting of a polyarginine (15 mer) randomized at the N-terminal 10 positions. The selected peptides bound more strongly to the RRE than the natural binding partner, Rev, and contained glutamine residues that were assumed to be important for recognition of the G-A base pair. In addition, the peptides were predicted to bind to the RRE in an alpha-helical conformation. In this study, in order to understand the mechanism of the interaction between the RRE and the putative alpha-helical glutamine-containing peptides, the amino acid requirements for high affinity binding were analyzed by a combinatorial approach using a bacterial system for detecting RNA-peptide interactions. A consensus peptide, the DLA peptide, was elucidated, which consists of a single glutamine residue within a polyarginine context with the glutamine residue flanked at specific positions by three nonarginine residues, two of which appear to be important for alpha-helix stabilization. In addition, the DLA peptide was found to bind extremely tightly to the RRE with an affinity 50-fold higher than that of the Rev peptide as determined by a gel shift assay. A working model for the interaction of the DLA peptide to the RRE is proposed, which should aid in the development of peptide-based drugs that inhibit HIV replication, as well as in our understanding of polypeptide-RNA interactions.

Role of non-coding RNAs in neurodegeneration and stress response in Drosophila

The inherent limitations of genetic analysis in humans and other mammals as well as striking conservation of most genes controlling nervous system functioning in flies and mammals made Drosophila an attractive model to investigate various aspects of brain diseases. Since RNA research has made great progress in recent years here we present an overview of studies demonstrating the role of various non-coding RNAs in neurodegeneration and stress response in Drosophila as a model organism. We put special emphasis on the role of non-coding micro RNAs, hsr-omega transcripts, and artificial small highly structured RNAs as triggers of neuropathology including aggregates formation, cognitive abnormalities and other symptoms. Cellular stress is a conspicuous feature of many neurodegenerative diseases and the production of specialized proteins protects the nerve cells against aggregates formation. Therefore, herein we describe some data implicating various classes of non-coding RNAs in stress response in Drosophila. All these findings highlight Drosophila as an important model system to investigate various brain diseases potentially mediated by some non-coding RNAs including polyglutamine diseases, Alzheimer's disease, Huntigton's disease, and many others.

The RNA-binding domain of bacteriophage P22 N protein is highly mutable, and a single mutation relaxes specificity toward lambda

Antitermination in bacteriophage P22, a lambdoid phage, uses the arginine-rich domain of the N protein to recognize boxB RNAs in the nut site of two regulated transcripts. Using an antitermination reporter system, we screened libraries in which each nonconserved residue in the RNA-binding domain of P22 N was randomized. Mutants were assayed for the ability to complement N-deficient virus and for antitermination with P22 boxB(left) and boxB(right) reporters. Single amino acid substitutions complementing P22 N(-) virus were found at 12 of the 13 positions examined. We found evidence for defined structural roles for seven nonconserved residues, which was generally compatible with the nuclear magnetic resonance model. Interestingly, a histidine can be replaced by any other aromatic residue, although no planar partner is obvious. Few single substitutions showed bias between boxB(left) and boxB(right), suggesting that the two RNAs impose similar constraints on genetic drift. A separate library comprising only hybrids of the RNA-binding domains of P22, lambda, and phi21 N proteins produced mutants that displayed bias. P22 N(-) plaque size plotted against boxB(left) and boxB(right) reporter activities suggests that lytic viral fitness depends on balanced antitermination. A few N proteins were able to complement both lambda N- and P22 N-deficient viruses, but no proteins were found to complement both P22 N- and phi21 N-deficient viruses. A single tryptophan substitution allowed P22 N to complement both P22 and lambda N(-). The existence of relaxed-specificity mutants suggests that conformational plasticity provides evolutionary transitions between distinct modes of RNA-protein recognition.

Tailoring the peptide-binding specificity of an RNA by combinations of specificity-altering mutations

In this study, the ability to tailor the peptide-binding specificity of an RNA was investigated. First, variants of the Rev-response element (RRE) RNA with different specificities toward the natural binding partner, Rev, and two RRE-binding aptamers, the RSG-1.2 and the Kl peptides, were identified. Next, hybrid RRE mutants with combinations of two sets of specificity-altering substitutions were tested for peptide-binding specificity. It was shown that in most cases the results of the combination of individual mutations were of an additive nature, therefore providing a way to manipulate the peptide-binding specificity of an RNA in a predictable manner.

Non-coding RNA as a trigger of neuropathologic disorder phenotypes in transgenic Drosophila

At most, many protein-misfolding diseases develop as environmentally induced sporadic disorders. Recent studies indicate that the dynamic interplay between a wide repertoire of noncoding RNAs and the environment play an important role in brain development and pathogenesis of brain disorders. To elucidate this new issue, novel animal models which reproduce the most prominent disease manifestations are required. For this, transgenic Drosophila strains were constructed to express small highly structured, non-coding RNA under control of a heat shock promoter. Expression of the RNA induced formation of intracellular aggregates revealed by Thioflafin T in embryonic cell culture and Congo Red in the brain of transgenic flies. Also, this strongly perturbed the brain control of locomotion monitored by the parameters of sound production and memory retention of young 5-day-old males. This novel model demonstrates that expression of non-coding RNA alone is sufficient to trigger neuropathology.

Recognition of a flipped base in a hairpinloop DNA by a small peptide

Two tiny hairpin DNAs, CORE (dAGGCTTCGGCCT) and AP2 (dAGGCTXCGGCCT; X: abasic nucleotide), fold into almost the same tetraloop hairpin structure with one exception, that is, the sixth thymine (T6) of CORE is exposed to the solvent water (Kawakami, J. et al., Chem. Lett. 2001, 258-259). In the present study, we selected small peptides that bind to CORE or AP2 from a combinatorial pentapeptide library with 2.5 x 10(6) variants. On the basis of the structural information, the selected peptide sequences should indicate the essential qualifications for recognition of the hairpin loop DNA with and without a flipped base. In the selected DNA binding peptides, aromatic amino acids such as histidine for CORE and glutamine/aspartic acid for AP2 were found to be abundant amino acids. This amino acid preference suggests that CORE-binding peptides use pi-pi stacking to recognize the target while hydrogen bonding is dominant for AP2-binding peptides. To investigate the binding properties of the selected peptide to the target, surface plasmon resonance was used. The binding constant of the interaction between CORE and a CORE-binding peptide (HWHHE) was about 1.1 x 10(6) M(-1) at 25 degrees C and the resulting binding free energy change at 25 degrees C (DeltaG degrees (25)) was -8.2 kcal mol(-1). The binding of the peptide to AP2 was also analyzed and the resulting binding constant and DeltaG degrees (25) were about 4.2 x 10(4) M(-1) and -6.3 kcal mol(-1), respectively. The difference in the binding free energy changes (DeltaDeltaG degrees (25)) of 1.9 kcal mol(-1) was comparable to the values reported in other systems and was considered a consequence of the loss of pi-pi stacking. Moreover, the stabilization effect by stacking affected the dissociation step as well as the association step. Our results suggest that the existence of an aromatic ring (T6 base) produces new dominant interactions between peptides and nucleic acids, although hydrogen bonding is the preferable mode of interaction in the absence of the flipping base. These findings regarding CORE and AP2 recognition are expected to give useful information in the design of novel artificial DNA binding peptides.

Bacteriophage P22 antitermination boxB sequence requirements are complex and overlap with those of lambda

Transcription antitermination in phages lambda and P22 uses N proteins that bind to similar boxB RNA hairpins in regulated transcripts. In contrast to the lambda N-boxB interaction, the P22 N-boxB interaction has not been extensively studied. A nuclear magnetic resonance structure of the P22 N peptide boxB(left) complex and limited mutagenesis have been reported but do not reveal a consensus sequence for boxB. We have used a plasmid-based antitermination system to screen boxBs with random loops and to test boxB mutants. We find that P22 N requires boxB to have a GNRA-like loop with no simple requirements on the remaining sequences in the loop or stem. U:A or A:U base pairs are strongly preferred adjacent to the loop and appear to modulate N binding in cooperation with the loop and distal stem. A few GNRA-like hexaloops have moderate activity. Some boxB mutants bind P22 and lambda N, indicating that the requirements imposed on boxB by P22 N overlap those imposed by lambda N. Point mutations can dramatically alter boxB specificity between P22 and lambda N. A boxB specific for P22 N can be mutated to lambda N specificity by a series of single mutations via a bifunctional intermediate, as predicted by neutral theories of evolution.

Specificity of Phage Display Selected Peptides for Modified Anticodon Stem and Loop Domains of tRNA

Protein recognition of RNA has been studied using Peptide Phage Display Libraries, but in the absence of RNA modifications. Peptides from two libraries, selected for binding the modified anticodon stem and loop (ASL) of human tRNA(LyS3) having 2-thiouridine (s(2)U34) and pseudouridine (psi39), bound the modified human ASL(Lys3)(s(2)U34;psi39) preferentially and had significant homology with RNA binding proteins. Selected peptides were narrowed to a manageable number using a less sensitive, but inexpensive assay before conducting intensive characterization. The affinity and specificity of the best binding peptide (with an N-terminal fluorescein) were characterized by fluorescence spectrophotometry. The peptide exhibited the highest binding affinity for ASL(LYS3)(s(2)U34; psi39), followed by the hypermodified ASL(Lys3) (mcm(5)s(2) U34; ms(2)t(6)A37) and the unmodified ASL(Lys3), but bound poorly to singly modified ASL(Lys3) constructs (psi39, ms(2)t(6)A37, s(2)34), ASL(Lys1,2) (t(6)A37) and Escherichia coli ASL(Glu) (s(2)U34). Thus, RNA modifications are potentially important recognition elements for proteins and can be targets for selective recognition by peptides.