Influence of RNA Binding on the Structure and Dynamics of the Lassa Virus Nucleoprotein

Lassa virus protects its viral genome through the formation of a ribonucleoprotein complex in which the nucleoprotein (NP) encapsidates the single-stranded RNA genome. Crystal structures provide evidence that a conformational change must occur to allow for RNA binding. In this study, the mechanism by which NP binds to RNA and how the conformational changes in NP are achieved was investigated with molecular-dynamics simulations. NP was structurally characterized in an open configuration when bound to RNA and in a closed form in the absence of RNA. Our results show that when NP is bound to RNA, the protein is highly dynamic and the system undergoes spontaneous deviations away from the open-state configuration. The equilibrium simulations are supported by free-energy calculations that quantify the influence of RNA on the free-energy surface, which governs NP dynamics. We predict that the globally stable states are qualitatively in agreement with the observed crystal structures, but that both open and closed conformations are thermally accessible in the presence of RNA. The free-energy calculations also provide a prediction of the location of the transition state for RNA binding and identify an intermediate metastable state that exhibits correlated motions that could promote RNA binding.     

Conformational Transition Pathway in the Inhibitor Binding Process of Human Monoacylglycerol Lipase

Human monoacylglycerol lipase (MGL) catalyzes the hydrolysis of 2-arachidonoylglycerol to arachidonic and glycerol, which plays a pivotal role in the normal biological processes of brain. Co-crystal structure of the MGL in complex with its inhibitor, compound 1, shows that the helix α4 undergoes large-scale conformational changes in response to the compound 1 binding compared to the apo MGL. However, the detailed conformational transition pathway of the helix α4 in the inhibitor binding process of MGL has remained unclear. Here, conventional molecular dynamics (MD) and nudged elastic band (NEB) simulations were performed to explore the conformational transition pathway of the helix α4. Conventional MD simulations unveiled that the compound 1 induced the closed conformation of the active site of MGL, reduced the conformational flexibility of the helix α4, and elicited the large-scale conformational rearrangement of the helix α4, leading to the complete folding of the helix α4. Moreover, NEB simulations revealed that the conformational transition pathway of helix α4 underwent an almost 180° counter-clockwise rotation of the helix α4. Our computational results advance the structural and mechanistic understanding of the inhibitory mechanism.     

RNA binding by a novel helical fold of b2 protein from wuhan nodavirus mediates the suppression of RNA interference and promotes b2 dimerization

Wuhan nodavirus (WhNV) is a newly identified member of the Nodaviridae family with a bipartite genome of positive-sense RNAs. A nonstructural protein encoded by subgenomic RNA3 of nodaviruses, B2, serves as a potent RNA silencing suppressor (RSS) by sequestering RNA duplexes. We have previously demonstrated that WhNV B2 blocks RNA silencing in cultured Drosophila cells. However, the molecular mechanism by which WhNV B2 functions remains unknown. Here, we successfully established an RNA silencing system in cells derived from Pieris rapae, a natural host of WhNV, by introducing into these cells double-stranded RNA (dsRNA)-expressing plasmids or chemically synthesized small interfering RNAs (siRNAs). Using this system, we revealed that the WhNV B2 protein inhibited Dicer-mediated dsRNA cleavage and the incorporation of siRNA into the RNA-induced silencing complex (RISC) by sequestering dsRNA and siRNA. Based on the modeled B2 3-dimensional structure, serial single alanine replacement mutations and N-terminal deletion analyses showed that the RNA-binding domain of B2 is formed by its helices α2 and α3, while helix α1 mediates B2 dimerization. Furthermore, positive feedback between RNA binding and B2 dimerization was uncovered by gel shift assay and far-Western blotting, revealing that B2 dimerization is required for its binding to RNA, whereas RNA binding to B2 in turn promotes its dimerization. All together, our findings uncovered a novel RNA-binding mode of WhNV B2 and provided evidence that the promotion effect of RNA binding on dimerization exists in a viral RSS protein.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Targeting human DNA polymerase α for the inhibition of keratinocyte proliferation. Part 1. Homology model,active site architecture and ligand binding

In order to understand the binding modes of human DNA polymerase α (pol α) inhibitors on a molecular level, a 3D homology model of the active site of the enzyme was proposed based on the application of molecular modelling methods and molecular dynamic simulations using available crystal coordinates of pol α relatives. Docking results for a series of known nucleotide analogue inhibitors were consistent with reported experimental binding data and offered the possibility to elucidate structure-activity relationships via investigations of active site-inhibitor interactions. Furthermore, the study could explain, at least partially, the inhibitory effect of aphidicolin on pol α. In molecular dynamics simulations, aphidicolin occupied the catalytic centre, but acted in a not truly competitive manner with respect to nucleotides. It destabilized the replicating “closed” form of the pol α and transferred the enzyme into the inactive “open” conformation. This result is consistent with recent experiments on the binding mode of aphidicolin.     

The influenza A virus NS1 protein interacts with the nucleoprotein of viral ribonucleoprotein complexes

The influenza A virus genome consists of eight RNA segments that associate with the viral polymerase proteins (PB1, PB2, and PA) and nucleoprotein (NP) to form ribonucleoprotein complexes (RNPs). The viral NS1 protein was previously shown to associate with these complexes, although it was not clear which RNP component mediated the interaction. Using individual TAP (tandem affinity purification)-tagged PB1, PB2, PA, and NP, we demonstrated that the NS1 protein interacts specifically with NP and not the polymerase subunits. The region of NS1 that binds NP was mapped to the RNA-binding domain.     

RNA binding by the novel helical domain of the influenza virus NS1 protein requires its dimer structure and a small number of specific basic amino acids

The RNA-binding/dimerization domain of the NS1 protein of influenza A virus (73 amino acids in length) exhibits a novel dimeric six-helical fold. It is not known how this domain binds to its specific RNA targets, one of which is double-stranded RNA. To elucidate the mode of RNA binding, we introduced single alanine replacements into the NS1 RNA-binding domain at specific positions in the three-dimensional structure. Our results indicate that the dimer structure is essential for RNA binding, because any alanine replacement that causes disruption of the dimer also leads to the loss of RNA-binding activity. Surprisingly, the arginine side chain at position 38, which is in the second helix of each monomer, is the only amino-acid side chain that is absolutely required only for RNA binding and not for dimerization, indicating that this side chain probably interacts directly with the RNA target. This interaction is primarily electrostatic, because replacement of this arginine with lysine had no effect on RNA binding. A second basic amino acid, the lysine at position 41, which is also in helix 2, makes a strong contribution to the affinity of binding. We conclude that helix 2 and helix 2', which are antiparallel and next to each other in the dimer conformation, constitute the interaction face between the NS1 RNA-binding domain and its RNA targets, and that the arginine side chain at position 38 and possibly the lysine side chain at position 41 in each of these antiparallel helices contact the phosphate backbone of the RNA target.     

Minimum sequence requirements for the binding of paromomycin to the rRNA decoding site A

We have recently introduced a computational methodology that combines molecular dynamics (MD) simulations, free-energy calculations, and in vitro binding assays to predict the minimum RNA structural requirements for selective, high-affinity RNA binding to small-molecule ligands. Here, we show that this methodology can be applied to the conformationally flexible aminoglycoside antibiotic paromomycin. A RNA consisting of an 11-mer:10-mer duplex that contains one 16S ribosome RNA decoding A-site bound to paromomycin was simulated for 4 ns. The methodology predicts that the 11-mer:10-mer duplex binds to paromomycin with high affinity, whereas smaller RNA duplexes lose complex stability and the ability to bind paromomycin. The predicted high-affinity binding to paromomycin of the 11-mer:10-mer duplex was confirmed experimentally (EC(50) = 0.28 microM), as well as the inability of smaller complexes to bind. Our simulations show good agreement with experiment for dynamic and structural properties of the isolated A-site, including hydrogen-bonding networks and RNA structural rearrangements upon ligand binding. The results suggest that MD simulations can supplement in vitro methods as a tool for predicting minimum RNA-binding motifs for both small, rigid ligands, and large, flexible ligands when structural information is available.     

Conformational and dynamics simulation study of antimicrobial peptide hedistin—heterogeneity of its helix–turn–helix motif

Hedistin is an antimicrobial peptide isolated from the coelomocytes of Nereis diversicolor, possessing activity against a large spectrum of bacteria including the methicillin resistant Staphylococcus aureus and Vibrio alginolyticus. The three-dimensional structure of hedistin in both aqueous solution and deuterated dodecylphosphocholine (DPC) micelles was examined using circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques. And, the early events of the antibacterial process of hedistin were simulated using palmitoyl-oleoyl-phophatidylcholine (POPC) lipid bilayers and molecular dynamics (MD) simulation methods. Hedistin lacks secondary structure in aqueous solution, however, in DPC micelles, it features with a heterogeneous helix–turn–helix moiety and exhibits obvious amphipathic nature. The turn region (residues Val9–Thr12) in the moiety is a four-residue hinge, lying in between the first N-terminal α-helix (residues Leu5–Lys8) and the second α-helix (residues Val13–Ala17) regions and causing an ～ 120° angle between the axes of the two helices. The segmental and nonlinear nature of hedistin structure is referred to as the heterogeneity of its helix–turn–helix motif which was found to be corresponding to a kind of discrete dynamics behavior, herein coined as its dynamical heterogeneity, at the early stage (0–50 ns) of the MD simulations. That is, the first helix segment, prior to (at 310 K) or following (at 363 K) the second helix, binds to the lipid head-group region and subsequently permeates into the hydrophobic lipid tail region, and the hinge is the last portion entering the lipid environment. This result implies that hedistin may adopt a “carpet” model action when disrupting bacterial membrane.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Identification of an RNA binding region within the N-terminal third of the influenza A virus nucleoprotein.

The influenza A virus nucleoprotein (NP) has been examined with regard to its RNA-binding characteristics. NP, purified from virions and devoid of RNA, bound synthetic RNAs in vitro and interacted with the ribonucleotide homopolymers poly(A), poly(G), poly(U), and poly(C) in a salt-dependent manner, showing higher binding affinity for polypyrimidine homopolymers. To map the NP regions involved in RNA binding, a series of deleted forms of the NP were prepared, and these truncated polypeptides were tested for their ability to bind poly(U) and poly(C) homopolymers linked to agarose beads. Proteins containing deletions at the N terminus of the NP molecule showed reduced RNA-binding activity, indicating that this part of the protein was required to bind RNA. To identify the NP region or regions which directly interact with RNA, proteins having the maltose-binding protein fused with various NP fragments were obtained and tested for binding to radioactively labeled RNAs in three different assays: (i) nitrocellulose filter binding assays, (ii) gel shift assays, and (iii) UV light-induced cross-linking experiments. A maltose-binding protein fusion containing the N-terminal 180 amino acids of NP behaved as an RNA-binding protein in the three assays, demonstrating that the N terminus of NP can directly interact with RNA. This NP region could be further subdivided into two smaller regions (amino acids 1 to 77 and 79 to 180) that also retained RNA-binding activity.     

Dominant-negative effects in prion diseases: insights from molecular dynamics simulations on mouse prion protein chimeras

Mutations in the prion protein (PrP) can cause spontaneous prion diseases in humans (Hu) and animals. In transgenic mice, mutations can determine the susceptibility to the infection of different prion strains. Some of these mutations also show a dominant-negative effect, thus halting the replication process by which wild type mouse (Mo) PrP is converted into Mo scrapie. Using all-atom molecular dynamics (MD) simulations, here we studied the structure of HuPrP, MoPrP, 10?Hu/MoPrP chimeras, and 1 Mo/sheepPrP chimera in explicit solvent. Overall, ～2?μs of MD were collected. Our findings suggest that the interactions between α1 helix and N-terminal of α3 helix are critical in prion propagation, whereas the β2–α2 loop conformation plays a role in the dominant-negative effect.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:4.