Recognition mechanism of siRNA by viral p19 suppressor of RNA silencing: a molecular dynamics study

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.

Structural insights into interactions between viral suppressor of RNA silencing protein p19 mutants and small RNAs

Viral suppressors of RNA silencing (VSRSs) are a diverse group of viral proteins that have evolved to disrupt eukaryotic RNA silencing pathways, thereby contributing to viral pathogenicity. The p19 protein is a VSRS that selectively binds to short interfering RNAs (siRNAs) over microRNAs (miRNAs). Mutational analysis has identified single amino acid substitutions that reverse this selectivity through new high-affinity interactions with human miR-122. Herein, we report crystal structures of complexed p19-T111S (2.6 Å), p19-T111H (2.3 Å) and wild-type p19 protein (2.2 Å) from the Carnation Italian ringspot virus with small interfering RNA (siRNA) ligands. Structural comparisons reveal that these mutations do not lead to major changes in p19 architecture, but instead promote subtle rearrangement of residues and solvent molecules along the p19 midline. These observations suggest p19 uses many small interactions to distinguish siRNAs from miRNAs and perturbing these interactions can create p19 variants with novel RNA-recognition properties. DATABASE: Model data are deposited in the PDB database under the accession numbers 6BJG, 6BJH and 6BJV.

Exploration of HIV-1 fusion peptide–antibody VRC34.01 binding reveals fundamental neutralization sites

VRC34.01 antibody binding to a vulnerable site of HIV envelope glycoprotein (Env), the gp41 fusion peptide, renders robust HIV neutralization, but several critical mutations decrease binding affinity and result in unbinding.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Interacting protein partners of Arabidopsis RNA-binding protein AtRBP45b

RNA binding proteins, important players in post-transcriptional gene regulation, usually exist in ribonuclear complexes. However, even in model systems like Arabidopsis characterisation of RBP associated proteins is limited. In this study, we investigated the interacting proteins of the Arabidopsis AtRBP45b, which is involved in stress signalling. In vivo localisation of AtRBP45b was conducted using 35S-GFP. FLAG-tagged AtRBP45b under control of the 35S promoter in the Atrbp45b-1 mutant background was used to pull down AtRBP45b interacting proteins. Yeast two-hybrid analysis, fluorescence energy resonance transfer assays were used to confirm the veracity of the AtRBP45b interacting proteins. In planta GFP-tagging indicated AtRBP45b is localised to the nucleus and the cytosol. AtRBP45b protein has a N-terminal proline-rich region and a C-terminal glutamine-rich domain that are usually involved in protein-protein interactions. Co-immunoprecipitation followed by mass spectrometry-based protein sequencing led to identification of 30 proteins that interacted with AtRBP45b. Using information from interactome databases (BIOGRID, INTACT and STRING), pull-down assays and localisation data, 12 putative interacting proteins were selected for yeast two-hybrid analysis. Cap-binding protein (CBP20, At5g44200) and polyA-binding protein (PAB8, At1g49760) were shown to interact with AtRBP45b. Based on its interacting partners we speculate that AtRBP45b may play an important role in RNA metabolism, especially in aspects related to mRNA stability and translation initiation during stress conditions in plants.

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13-26) with hydrogen bonds and π-π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16-20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM-PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (-43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.

Steered molecular dynamics identifies critical residues of the Nodamura virus B2 suppressor of RNAi

Nearly all RNA viruses produce double-stranded RNA (dsRNA) during their replication cycles--an important pathogen-associated molecular pattern recognized by the RNA interference (RNAi) pathway in invertebrates and plants. Nodamura virus (NoV) encodes a suppressor of RNA silencing termed B2, which binds to dsRNA and prevents the initiation of RNAi as well as the loading of silencing complexes. Using the published crystal structure of NoV-B2, we performed a series of molecular dynamics (MD) simulations to determine the relative electrostatic and van der Waals contributions of various residues in binding dsRNA, identifying four novel potential interactors: R56, E48, P68 and R69. Additionally, steered MD was used to simulate the binding affinity of NoV-B2 sequences bearing substitutions at positions F49, R56 or R59 to dsRNA, with F49S and R56L/R59L substitutions found to have a significant negative impact on the ability of NoV-B2 to bind dsRNA. NoV RNA1 variants were tested for self-directed replication in both vertebrate (RNAi⁻) and invertebrate (RNAi⁺) cultured cells. Consistent with a role in dsRNA binding, NoV replication in F49C and F49S variant constructs was affected negatively only in RNAi⁺ cells. Thus, we used a combination of MD simulations and experimental mutagenesis to further characterize residues important for NoV-dsRNA interactions.

pH-sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short-interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH-dependent residues that affect CIRV p19-siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue. At high pH, the deprotonation of Lys60, Lys67, Lys71, and Cys134 has the largest effect on the binding stability. Similarly, deprotonation of several acidic residues (Asp9, Glu12, Asp20, Glu35, and/or Glu41) at low pH results in a decrease in binding stability. At neutral pH, residues Glu17 and His132 provide a small increase in the binding stability and we find that the optimal pH range for siRNA binding is between 7.0 and 10.0. Overall, our findings further inform recent experiments and are in excellent agreement with data on the pH-dependent binding profile.

Large domain motions in Ago protein controlled by the guide DNA-strand seed region determine the Ago-DNA-mRNA complex recognition process

The recognition mechanism and cleavage activity of argonaute (Ago), miRNA, and mRNA complexes are the core processes to the small non-coding RNA world. The 5′ nucleation at the ‘seed’ region (position 2–8) of miRNA was believed to play a significant role in guiding the recognition of target mRNAs to the given miRNA family. In this paper, we have performed all-atom molecular dynamics simulations of the related and recently revealed Ago-DNA:mRNA ternary complexes to study the dynamics of the guide-target recognition and the effect of mutations by introducing “damaging” C·C mismatches at different positions in the seed region of the DNA-RNA duplex. Our simulations show that the A-form-like helix duplex gradually distorts as the number of seed mismatches increases and the complex can survive no more than two such mismatches. Severe distortions of the guide-target heteroduplex are observed in the ruinous 4-sites mismatch mutant, which give rise to a bending motion of the PAZ domain along the L1/L2 “hinge-like” connection segment, resulting in the opening of the nucleic-acid-binding channel. These long-range interactions between the seed region and PAZ domain, moderated by the L1/L2 segments, reveal the central role of the seed region in the guide-target strands recognition: it not only determines the guide-target heteroduplex’s nucleation and propagation, but also regulates the dynamic motions of Ago domains around the nucleic-acid-binding channel.

Studying the RNA silencing pathway with the p19 protein

The origins of the RNA silencing pathway are in defense against invading viruses and in response, viruses have evolved counter-measures to interfere with the host pathway. The p19 protein is expressed by tombusviruses as a suppressor of RNA silencing and functions to sequester small RNA duplexes, thereby preventing induction of the pathway. p19 exhibits size-specific and sequence-independent binding of its small RNA ligands, binding with high affinity to duplexes 20-22 nucleotides long. p19's binding specificity and its ability to sequester small RNAs has made it a unique protein-based tool for probing the molecular mechanisms of the highly complex RNA silencing pathway in a variety of systems. Furthermore, protein engineering of this 'molecular caliper' promises novel applications in biotechnology and medicine where small RNA molecules are of remarkable interest given their potent gene regulatory abilities.

Free-energy simulations reveal that both hydrophobic and polar interactions are important for influenza hemagglutinin antibody binding

Antibodies binding to conserved epitopes can provide a broad range of neutralization to existing influenza subtypes and may also prevent the propagation of potential pandemic viruses by fighting against emerging strands. Here we propose a computational framework to study structural binding patterns and detailed molecular mechanisms of viral surface glycoprotein hemagglutinin (HA) binding with a broad spectrum of neutralizing monoclonal antibody fragments (Fab). We used rigorous free-energy perturbation (FEP) methods to calculate the antigen-antibody binding affinities, with an aggregate underlying molecular-dynamics simulation time of several microseconds (∼2 μs) using all-atom, explicit-solvent models. We achieved a high accuracy in the validation of our FEP protocol against a series of known binding affinities for this complex system, with <0.5 kcal/mol errors on average. We then introduced what to our knowledge are novel mutations into the interfacial region to further study the binding mechanism. We found that the stacking interaction between Trp-21 in HA2 and Phe-55 in the CDR-H2 of Fab is crucial to the antibody-antigen association. A single mutation of either W21A or F55A can cause a binding affinity decrease of ΔΔG > 4.0 kcal/mol (equivalent to an ∼1000-fold increase in the dissociation constant K(d)). Moreover, for group 1 HA subtypes (which include both the H1N1 swine flu and the H5N1 bird flu), the relative binding affinities change only slightly (< ±1 kcal/mol) when nonpolar residues at the αA helix of HA mutate to conservative amino acids of similar size, which explains the broad neutralization capability of antibodies such as F10 and CR6261. Finally, we found that the hydrogen-bonding network between His-38 (in HA1) and Ser-30/Gln-64 (in Fab) is important for preserving the strong binding of Fab against group 1 HAs, whereas the lack of such hydrogen bonds with Asn-38 in most group 2 HAs may be responsible for the escape of antibody neutralization. These large-scale simulations may provide new insight into the antigen-antibody binding mechanism at the atomic level, which could be essential for designing more-effective vaccines for influenza.

Viral RNAi suppressor reversibly binds siRNA to outcompete Dicer and RISC via multiple turnover

RNA interference is a conserved gene regulatory mechanism employed by most eukaryotes as a key component of their innate immune response to viruses and retrotransposons. During viral infection, the RNase-III-type endonuclease Dicer cleaves viral double-stranded RNA into small interfering RNAs (siRNAs) 21-24 nucleotides in length and helps load them into the RNA-induced silencing complex (RISC) to guide the cleavage of complementary viral RNA. As a countermeasure, many viruses have evolved viral RNA silencing suppressors (RSS) that tightly, and presumably quantitatively, bind siRNAs to thwart RNA-interference-mediated degradation. Viral RSS proteins also act across kingdoms as potential immunosuppressors in gene therapeutic applications. Here we report fluorescence quenching and electrophoretic mobility shift assays that probe siRNA binding by the dimeric RSS p19 from Carnation Italian Ringspot Virus, as well as by human Dicer and RISC assembly complexes. We find that the siRNA:p19 interaction is readily reversible, characterized by rapid binding [(1.69 ± 0.07) × 10(8) M(-)(1) s(-1)] and marked dissociation (k(off)=0.062 ± 0.002 s(-1)). We also observe that p19 efficiently competes with recombinant Dicer and inhibits the formation of RISC-related assembly complexes found in human cell extract. Computational modeling based on these results provides evidence for the transient formation of a ternary complex between siRNA, human Dicer, and p19. An expanded model of RNA silencing indicates that multiple turnover by reversible binding of siRNAs potentiates the efficiency of the suppressor protein. Our predictive model is expected to be applicable to the dosing of p19 as a silencing suppressor in viral gene therapy.

Multiple functions of Rice dwarf phytoreovirus Pns10 in suppressing systemic RNA silencing

RNA silencing is a potent mechanism of antiviral defense response in plants and other organisms. For counterdefense, viruses have evolved a variety of suppressors of RNA silencing (VSRs) that can inhibit distinct steps of a silencing pathway. We previously identified Pns10 encoded by Rice dwarf phytoreovirus (RDV) as a VSR, the first of its kind from double-stranded RNA (dsRNA) viruses. In this study we investigated the mechanisms of Pns10 function in suppressing systemic RNA silencing in the widely used Nicotiana benthamiana model plant. We report that Pns10 suppresses local and systemic RNA silencing triggered by sense mRNA, enhances viral replication and/or viral RNA stability in inoculated leaves, accelerates the systemic spread of viral infection, and enables viral invasion of shoot apices. Mechanistically, Pns10 interferes with the perception of silencing signals in recipient tissues, binds double-stranded small interfering RNA (siRNAs) with two-nucleotide 3' overhangs, and causes the downregulated expression of RDR6. These results significantly deepen our mechanistic understanding of the VSR functions encoded by a dsRNA virus and contribute additional evidence that binding siRNAs and interfering with RDR6 expression are broad mechanisms of VSR functions encoded by diverse groups of viruses.