Molecular dynamics simulations of Ago silencing complexes reveal a large repertoire of admissible 'seed-less' targets

Probing the conformational dynamics of an Ago-RNA complex in water/methanol solution

Argonaute (Ago) proteins mediate target recognition guiding miRNA to bind complementary mRNA primarily in the seed region. However, additional pairing can occur beyond the seed, forming a supplementary duplex that can contribute to the guide-target affinity. In order to shed light on the connection, between protein-RNA interactions and miRNA-mRNA seed and supplementary duplex mobility, we carried out molecular dynamics simulations at the microsecond time-scale using a different approach compared to the ones normally used. Until now, theoretical investigations with classical MD on Ago-RNA complexes have been focused primarily on pure water solvent, which mimics the natural environment of biological molecules. Here, we explored the conformational space of a human Ago2 (hAgo2) bound to the seed + supplementary miRNA-mRNA duplex, using the solvent environment as a molecular probe. MD simulations have been performed in a mixture of water/MeOH at a molar ratio of 70 : 30 as well as in pure water for comparison. Our findings revealed that the mixed solvent promotes protein RNA association, principally enhancing salt-linkages between basic amino acid side-chains and acidic phosphates of the sugar-phosphate backbone. The primary effect registered was the restriction of supplementary duplex flexibility and the stabilization of the miRNA 3' terminus. Interestingly, we observed that the influence of the solvent appears to have almost no impact on the conformation of the seed duplex.

Novel Inhibitory Role of Fenofibric Acid by Targeting Cryptic Site on the RBD of SARS-CoV-2

The emergence of the recent pandemic causing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has created an alarming situation worldwide. It also prompted extensive research on drug repurposing to find a potential treatment for SARS-CoV-2 infection. An active metabolite of the hyperlipidemic drug fenofibrate (also called fenofibric acid or FA) was found to destabilize the receptor-binding domain (RBD) of the viral spike protein and therefore inhibit its binding to human angiotensin-converting enzyme 2 (hACE2) receptor. Despite being considered as a potential drug candidate for SARS-CoV-2, FA's inhibitory mechanism remains to be elucidated. We used molecular dynamics (MD) simulations to investigate the binding of FA to the RBD of the SARS-CoV-2 spike protein and revealed a potential cryptic FA binding site. Free energy calculations were performed for different FA-bound RBD complexes. The results suggest that the interaction of FA with the cryptic binding site of RBD alters the conformation of the binding loop of RBD and effectively reduces its binding affinity towards ACE2. Our study provides new insights for the design of SARS-CoV-2 inhibitors targeting cryptic sites on the RBD of SARS-CoV-2.

Critical role of backbone coordination in the mRNA recognition by RNA induced silencing complex

Despite its functional importance, the molecular mechanism underlying target mRNA recognition by Argonaute (Ago) remains largely elusive. Based on extensive all-atom molecular dynamics simulations, we constructed quasi-Markov State Model (qMSM) to reveal the dynamics during recognition at position 6-7 in the seed region of human Argonaute 2 (hAgo2). Interestingly, we found that the slowest mode of motion therein is not the gRNA-target base-pairing, but the coordination of the target phosphate groups with a set of positively charged residues of hAgo2. Moreover, the ability of Helix-7 to approach the PIWI and MID domains was found to reduce the effective volume accessible to the target mRNA and therefore facilitate both the backbone coordination and base-pair formation. Further mutant simulations revealed that alanine mutation of the D358 residue on Helix-7 enhanced a trap state to slow down the loading of target mRNA. Similar trap state was also observed when wobble pairs were introduced in g6 and g7, indicating the role of Helix-7 in suppressing non-canonical base-paring. Our study pointed to a general mechanism for mRNA recognition by eukaryotic Agos and demonstrated the promise of qMSM in investigating complex conformational changes of biomolecular systems.

Length-Dependent Structural Transformations of Huntingtin PolyQ Domain Upon Binding to 2D-Nanomaterials

There is a strong negative correlation between the polyglutamine (polyQ) domain length (Q-length) in the intrinsically disordered Huntingtin protein (Htt) exon-1 and the age of onset of Huntington's disease (HD). PolyQ of Q-length longer than 40 has the propensity of forming very compact aggregate structures, leading to HD at full penetrance. Recent advances in nanobiotechnology provided a new platform for the development of novel diagnosis and therapeutics. Here, we explore the possibility of utilizing 2D-nanomaterials to inhibit the formation of supercompact polyQ structures through the so-called “folding-upon-binding” where the protein structure is dependent on the binding substrate. Using molecular dynamics simulations, we characterize two polyQ peptides with Q-length of 22 (Q22, normal length) and 46 (Q46, typical length causing HD) binding to both graphene and molybdenum disulfide (MoS₂) nanosheets, which have been applied as antibacterial or anticancer agents. Upon binding, Q22 unfolds and elongates on both grapheme and MoS₂ surfaces, regardless of its initial conformation, with graphene showing slightly stronger effect. In contrast, initially collapsed Q46 remains mostly collapsed within our simulation time on both nanosheets even though they do provide some “stretching” to Q46 as well. Further analyses indicate that the hydrophobic nature of graphene/MoS₂ promotes the stretching of polyQ on nanosheets. However, there is strong competition with the intra-polyQ interactions (mainly internal hydrogen bonds) leading to the disparate folding/binding behaviors of Q22 and Q46. Our results present distinct Q-length specific behavior of the polyQ domain upon binding to two types of 2D-nanomaterials which holds clinical relevance for Huntington's disease.

Directional extraction and penetration of phosphorene nanosheets to cell membranes

Recently, phosphorene, a novel two-dimensional nanomaterial with a puckered surface morphology, was shown to exhibit cytotoxicity, but its underlying molecular mechanisms remain unknown. Herein, using large scale molecular dynamics simulations, we show that phosphorene nanosheets can penetrate into and extract large amounts of phospholipids from the cell membranes due to the strong dispersion interaction between phosphorene and lipid molecules, which would reduce cell viability. The extracted phospholipid molecules are aligned along the wrinkle direction of the phosphorene nanosheet because of its unique puckered structure. Our results also reveal that small phosphorene nanosheets penetrate into the cell membrane in a specific direction which is determined by the size and surface topography of phosphorene and the thickness of the membrane. These findings might shed light on understanding phosphorene's cytotoxicity and would be helpful for the future potential biomedical applications of phosphorene, such as biosensors and antibacterial agents.

The molecular mechanism of robust macrophage immune responses induced by PEGylated molybdenum disulfide

Molybdenum disulfide (MoS2), a representative hexagonal transition metal dichalcogenide (TMD), has been extensively exploited in biomedical applications due to its unique physicochemical properties and biocompatibility. However, the lack of adequate data regarding how MoS2 activates immunological responses of macrophages remains a key concern for its risk assessment. Here, we employ a combined theoretical and experimental approach to investigate the interactions of MoS2 and PEGylated MoS2 (MoS2-PEG) with macrophages. We first perform molecular dynamics simulations to examine the atomic-detailed interactions of MoS2 and MoS2-PEG nanoflakes with a realistic model of the macrophage membrane. We show that a small MoS2 nanoflake (edge length of 2.86 nm) is capable of penetrating the macrophage membrane independent of its concentration. We also demonstrate that when initiated with a corner point-on configuration, the surface-bound PEG chains of MoS2-PEG hinder its membrane insertion process, leading to a prolonged passage through the membrane. Moreover, when placed in a face-on arrangement initially, the MoS2-PEG exhibits a lower binding free energy than pristine MoS2 after its adsorption on the membrane surface. The PEG chains can even insert and get buried in the outer leaflet of the membrane, providing additional contact for membrane adsorption. Our flow cytometric experiments then show that the responses of macrophages to either MoS2-PEG or MoS2 are significantly higher than that of the control (no nanomaterial stimulus), with MoS2-PEG eliciting stronger cytokine secretion than the pristine MoS2. The characteristics of slower/prolonged membrane penetration and stronger membrane adsorption of MoS2-PEG compared to pristine MoS2 explain why it triggers more sustained stimulation and higher cytokine secretion in macrophages as observed in our experiments. Our findings reveal the underlying molecular mechanism of how MoS2-PEG influences the immune responses and suggest its potential applications in nanomedicine involving immune stimulation.

Modeling and Structural Characterization of the Sweet Taste Receptor Heterodimer

Sweet taste receptor, a heterodimer belonging to the class C G-protein coupled receptor (GPCR) family and composed of the T1R2 and T1R3 subunits, is responsible for the perception of natural sugars, sweet proteins, various d-amino acids, as well as artificial sweeteners. Despite the critical importance of the sweet receptor not only in mediating gustation but also in its role in the food industry, the architecture of the T1R2-T1R3 complex and the mechanism by which extracellular stimuli induce conformational changes that are propagated to the intracellular milieu, i.e., the signal transduction pathway, remain largely unknown. Here, we constructed and characterized a full-length structural model of the T1R2-T1R3 receptor, including both the transmembrane (TM) and extracellular (EC) domains of the heterodimer, using comparative modeling and extensive all-atom molecular dynamics simulations. Several heterodimer interfaces were first examined for the TM domain, and conformational changes occurring at the intracellular side and associated with the receptor's activation were characterized. From the analysis on the simulated data, putative allosteric binding sites for ligands, ions, and cholesterol were proposed. Also, insights into the protein interface of the TM domain upon activation are provided. The EC domain of the heterodimer, including both the Venus flytrap and cysteine-rich domains, was also investigated. Several important intersubunit interactions located at regions responsible for the receptor's proper function were observed, which resemble those recently identified in other class C GPCR members. Integration of the results from the TM and EC domains facilitates the generation of a full-length T1R2-T1R3 receptor. These findings along with the full-length structural model of the T1R2-T1R3 receptor provide a structural framework that may assist in understanding the mechanistic details associated with the receptor activation process for the sweet T1R2-T1R3 receptor as well as other members of the same family.

Different platinum crystal surfaces show very distinct protein denaturation capabilities

Different platinum (Pt) surfaces of nanocrystals usually exhibit significant distinctions with regard to various biological, physical, and chemical characteristics, such as bio-recognition, surface wetting, and catalytic activities. In this study, we report for the first time that two shape-controlled Pt nanocrystals with the most common low-index surfaces, Pt(100) and Pt(111), show very dissimilar protein denaturation capabilities based on all-atom molecular dynamics simulations employing the widely used model protein, villin headpiece (HP35). We demonstrate that HP35 is well preserved on the Pt(100) crystal surface, whereas it is severely disrupted on the Pt(111) crystal surface. This surprising difference originates from the distinct water behavior in the first solvation shell (FSS) of the two Pt crystal surfaces. Within the FSS of the Pt(100) crystal surface, water molecules form a very compact and stable monolayer through a highly uniform rhombic hydrogen-bond network. This water monolayer prefers the adsorption of acidic residues (such as Glu and Asp) and acts as a shield to prevent other residues from directly coming into contact with the metal surface. On the other hand, the hydrogen bond network in the water monolayer in the FSS of the Pt(111) crystal surface is very sparse and quite defective, which makes it more vulnerable to the penetration of various residues, particularly those with planar side chains such as Phe, Trp and Arg due to strong dispersion interactions, leading to subsequent protein unfolding. The binding free energy calculations for some key amino acids on the two different crystal surfaces further uncover the molecular origin behind their distinct protein denaturation capability. Our study reveals the vital importance of interfacial water in determining the structure of proteins when binding to different metal crystal surfaces. The discovered molecular mechanisms may be helpful for the future development of a bio-assisted programmable synthetic strategy of sophisticated Pt nanostructures for biomedical applications.

Lanosterol Disrupts the Aggregation of Amyloid-β Peptides

Lanosterol, an amphipathic molecule, was discovered only very recently to effectively hinder the aggregation of lens proteins and dissolve the extremely stable fibrillar aggregates in cataracts. Here, we combined computational and experimental approaches to study how lanosterol disrupts the aggregation of another important peptide, amyloid-β (Aβ) peptide, associated with the Alzheimer's Disease (AD). Molecular dynamics simulations using the core amyloidogenic segment (KLVFFA) of Aβ peptide revealed that lanosterol exhibits at least two types of inhibition mechanism on the self-assembly of Aβ peptides. First, lanosterol entangles with peptides and forms a hydrophobic core with residues Phe-19 and Phe-20 in particular. Second, it interferes with the steric zipper interaction at the β-sheet-β-sheet interface. These simulation data suggest that lanosterol induces the unfolding of the Aβ peptide and the separation of the β-sheet layers. This predicted inhibition effect of lanosterol was then confirmed by an in vitro ThT fluorescence assay and AFM imaging. The cell toxicity assay also showed that the treatment of lanosterol indeed mitigates the cytotoxicity of the Aβ peptide in PC-12 cells. Moreover, lanosterol shows a stronger suppression effect on Aβ peptides' aggregation than cholesterol because of its higher hydrophobicity. This result establishes a foundation for the development of lanosterol-based potential therapies for AD and other protein conformational diseases.

Parameterization of Molybdenum Disulfide Interacting with Water Using the Free Energy Perturbation Method

Water contact angles (WCA) are often used to parametrize force field parameters of novel 2D nanomaterials, such as molybdenum disulfide (MoS₂), which has emerged as a promising nanomaterial in many biomedical applications due to its unique and impressive properties. However, there is a wide range of water-MoS₂ contact angles in the literature depending on the aging process on the surface of a MoS₂ nanosheet and/or substrate material. In this study, we revisit and optimize existing parameters for the basal plane of MoS₂ with two popular water models, TIP3P and SPC/E, using the wide range of WCAs from various experiments. We develop and deploy the free energy perturbation method for parametrizing MoS₂ with experimentally determined WCAs for both fresh and aged surfaces. Energy decomposition analysis on the simulation trajectories reveals that MoS₂-water interaction is dominated by van der Waals interaction, which mainly comes from the top layer of MoS₂. We conclude that to describe both fresh and aged MoS₂ surfaces it is convenient to only adjust the Lennard-Jones parameter ε_S (the depth of the potential well of a sulfur atom), which displays a surprisingly linear correlation with WCAs.

Understanding the molecular interaction of human argonaute-2 and miR-20a complex: A molecular dynamics approach

Argonaute-2 (AGO2), a member of the Argonaute family, is the only member possessing catalytic and RNA silencing activity. In here, a molecular dynamics (MDs) simulation was performed using the crystal structure of human AGO2 protein complex with miR-20a. miR-20a is involved with various kind of biological process like heart and lung development, oncogenic process, etc. In precise, MD simulation was carried out with AGO2 protein complex with wild type, two mutant sites and four mutant sites in guided microRNA (miRNA). It has been noted that root-mean-square deviation (RMSD) of atomic positions of nucleic acid for wild type and two mutant sites guided miRNA has the same pattern of fluctuations, which stabilizes around 0.27 nm after 2 ns. Cα atom of AGO2 protein in the complex shows that this complex with wild type and two mutant site mutation duplex has a stable RMSD value after 20 ns, ranging between 0.14 and 0.21 nm. From the root-mean-square fluctuation (RMSF), we observed an increased pattern of fluctuations for the atoms of four mutant complex of AGO2-miR-20a complex. This increased RMSF of non-mutated nucleic acids is contributed by U-A bond breaking at the site of the nucleotide of U2 of guided miRNA, as observed from the duplex structure taken at different time steps of the simulation. Superimposed structure of the miRNA-mRNA duplex for the three complexes depicts that the three miRNA-mRNA duplexes are stable during the simulation. Current work demonstrates the possible correlations between the conformational changes of this AGO2-miR-20a duplex structure and the interactions of different atoms.

Short RNA regulators: the past, the present, the future, and implications for precision medicine and health disparities

We herein provide a brief review of the trajectory that the field of short RNA research followed in the last 25 years. We place emphasis on the unexpected discoveries and the ramifications of these discoveries for the field, as well as offer some thoughts about what the next 25 years may bring. Arguably, the uncovered dependence of different types of short RNAs on individual attributes such as a person's sex, population origin, race, and on tissue type, tissue state, and disease was most unexpected. This dependence has important ramifications in that it will provide a boost to our understanding of the molecular mechanisms of health disparities as well as pave the way for novel approaches to designing improved and personalized diagnostics and therapeutics.

Robust Antibacterial Activity of Tungsten Oxide (WO3-x) Nanodots

Antibacterial agents are an important tool in the prevention of bacterial infections. Inorganic materials are attractive due to their high stability under a variety of conditions compared to organic antibacterial agents. Herein tungsten oxide nanodots (WO_3-x), synthesized by a simple one-pot synthetic approach, were found to exhibit strong antibacterial capabilities. The analyses with colony-forming units (CFU) showed an excellent antibacterial activity of WO_3-x against both Gram-negative Escherichia coli (E. coli) and Gram-positive Staphylococcus aureus (S. aureus) strains. The scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images revealed clear damages to the bacterial cell membranes, which was further confirmed by molecular dynamics simulations. Additionally, exposure to simulated sunlight was found to further increase the germicidal activity of WO_3-x nanodots, a 30 min exposure to sunlight combined with 50 μg/mL WO_3-x nanodots showed a 70% decrease in E. coli viability compared to without exposure. Electron spin resonance spectroscopy (ESR) was used to elucidate the underlying mechanism of this photocatalytic activity through the generation of hydroxyl radical species. The cell counting kit-8 (CCK-8) and the live/dead assay were further employed to evaluate the cytotoxicity of WO_3-x nanodots on eukaryotic cells, which demonstrated their general biocompatibility. In summary, our results suggest WO_3-x nanodots have considerable potential in antibacterial applications, while also being biocompatible at large.

Mapping the Pax6 3' untranslated region microRNA regulatory landscape

BACKGROUND

PAX6 is a homeodomain transcription factor that acts in a highly dosage-sensitive manner to regulate the development and function of the eyes, nose, central nervous system, gut, and endocrine pancreas. Several individual microRNAs (miRNA) have been implicated in regulating PAX6 in different cellular contexts, but a more general view of how they contribute to the fine-tuning and homeostasis of PAX6 is poorly understood.

RESULTS

Here, a comprehensive analysis of the Pax6 3' untranslated region was performed to map potential miRNA recognition elements and served as a backdrop for miRNA expression profiling experiments to identify potential cell/tissue-specific miRNA codes. Pax6 3'UTR pull-down studies identified a cohort of miRNA interactors in pancreatic αTC1-6 cells that, based on the spacing of their recognition sites in the Pax6 3'UTR, revealed 3 clusters where cooperative miRNA regulation may occur. Some of these interacting miRNAs have been implicated in α cell function but have not previously been linked to Pax6 function and may therefore represent novel PAX6 regulators.

CONCLUSIONS

These findings reveal a regulatory landscape upon which miRNAs may participate in the developmental control, fine-tuning and/or homeostasis of PAX6 levels.

Mechanism by which DHA inhibits the aggregation of KLVFFA peptides: A molecular dynamics study

Docosahexaenoic acid (DHA) is one of the omega-3 polyunsaturated fatty acids, which has shown promising applications in lowering Aβ peptide neurotoxicity in vitro by preventing aggregation of Aβ peptides and relieving accumulation of Aβ fibrils. Unfortunately, the underlying molecular mechanisms of how DHA interferes with the aggregation of Aβ peptides remain largely enigmatic. Herein, aggregation behaviors of amyloid-β(Aβ)_16-21 peptides (KLVFFA) with or without the presence of a DHA molecule were comparatively studied using extensive all-atom molecular dynamics simulations. We found that DHA could effectively suppress the aggregation of KLVFFA peptides by redirecting peptides to unstructured oligomers. The highly hydrophobic and flexible nature of DHA made it randomly but tightly entangled with Leu-17, Phe-19, and Phe-20 residues to form unstructured but stable complexes. These lower-ordered unstructured oligomers could eventually pass through energy barriers to form ordered β-sheet structures through large conformational fluctuations. This study depicts a microscopic picture for understanding the role and mechanism of DHA in inhibition of aggregation of Aβ peptides, which is generally believed as one of the important pathogenic mechanisms of Alzheimer's disease.

Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study

Thermus thermophiles Argonaute (TtAgo) is a complex, which is consisted of 5′-phosphorylated guide DNA and a series of target DNA with catalytic activities at high temperatures. To understand why high temperatures are needed for the catalytic activities, three molecular dynamics simulations and binding free energy calculations at 310, 324, and 338K were performed for the TtAgo-DNA complex to explore the conformational changes between 16-mer guide DNA/15-mer target DNA and TtAgo at different temperatures. The simulation results indicate that a collapse of a small β-strand (residues 507–509) at 310 K caused Glu512 to move away from the catalytic residues Asp546 and Asp478, resulting in a decrease in catalytic activity, which was not observed in the simulations at 324 and 338 K. The nucleic acid binding channel became enlarged at 324 and 338K, thereby facilitating the DNA to slide in. Binding free energy calculations and hydrogen bond occupancy indicated that the interaction between TtAgo and the DNA was more stable at 324K and 338K than at 310 K. The DNA binding pocket residues Lys575 and Asn590 became less solvent accessible at 324 and 338K than at 310 K to influence hydrophilic interaction with DNA. Our simulation studies shed some light on the mechanism of TtAgo and explained why a high temperature was needed by TtAgo during gene editing of CRISPR.

Novel and Haplotype Specific MicroRNAs Encoded by the Major Histocompatibility Complex

The MHC is recognized for its importance in human health and disease. However, many disease-associated variants throughout the region remain of unknown significance, residing predominantly within non-coding regions of the MHC. The characterization of non-coding RNA transcripts throughout the MHC is thus central to understanding the genetic contribution of these variants. Therefore, we characterize novel miRNA transcripts throughout the MHC by performing deep RNA sequencing of two B lymphoblastoid cell lines with completely characterized MHC haplotypes. Our analysis identifies 89 novel miRNA transcripts, 48 of which undergo Dicer-dependent biogenesis and are loaded onto the Argonaute silencing complex. Several of the identified mature miRNA and pre-miRNA transcripts are unique to specific MHC haplotypes and overlap common SNPs. Furthermore, 43 of the 89 identified novel miRNA transcripts lie within linkage disequilibrium blocks that contain a disease-associated SNP. These disease associated SNPs are associated with 65 unique disease phenotypes, suggesting that these transcripts may play a role in the etiology of numerous diseases associated with the MHC. Additional in silico analysis reveals the potential for thousands of putative pre-miRNA encoding loci within the MHC that may be expressed by different cell types and at different developmental stages.

Exploring PAZ/3'-overhang interaction to improve siRNA specificity. A combined experimental and modeling study

The understanding of the dynamical and mechanistic aspects that lie behind siRNA-based gene regulation is a requisite to boost the performance of siRNA therapeutics. A systematic experimental and computational study on the 3'-overhang structural requirements for the design of more specific and potent siRNA molecules was carried out using nucleotide analogues differing in structural parameters, such as sugar constraint, lack of nucleobase, distance between the phosphodiester backbone and nucleobase, enantioselectivity, and steric hindrance. The results established a set of rules governing the siRNA-mediated silencing, indicating that the thermodynamic stability of the 5'-end is a crucial determinant for antisense-mediated silencing but is not sufficient to avoid sense-mediated silencing. Both theoretical and experimental approaches consistently evidence the existence of a direct connection between the PAZ/3'-overhang binding affinity and siRNA's potency and specificity. An overall description of the systems is thus achieved by atomistic simulations and free energy calculations that allow us to propose a robust and self-contained procedure for studying the factors implied in PAZ/3'-overhang siRNA interactions. A higher RNAi activity is associated with a moderate-to-strong PAZ/3'-overhang binding. Contrarily, lower binding energies compromise siRNA potency, increase specificity, and favor siRNA downregulation by Ago2-independent mechanisms. This work provides in-depth details for the design of powerful and safe synthetic nucleotide analogues for substitution at the 3'-overhang, enabling some of the intrinsic siRNA disadvantages to be overcome.

Differential Pd-nanocrystal facets demonstrate distinct antibacterial activity against Gram-positive and Gram-negative bacteria

Noble metal-based nanomaterials have shown promise as potential enzyme mimetics, but the facet effect and underlying molecular mechanisms are largely unknown. Herein, with a combined experimental and theoretical approach, we unveil that palladium (Pd) nanocrystals exhibit facet-dependent oxidase and peroxidase-like activities that endow them with excellent antibacterial properties via generation of reactive oxygen species. The antibacterial efficiency of Pd nanocrystals against Gram-positive bacteria is consistent with the extent of their enzyme-like activity, that is {100}-faceted Pd cubes with higher activities kill bacteria more effectively than {111}-faceted Pd octahedrons. Surprisingly, a reverse trend of antibacterial activity is observed against Gram-negative bacteria, with Pd octahedrons displaying stronger penetration into bacterial membranes than Pd nanocubes, thereby exerting higher antibacterial activity than the latter. Our findings provide a deeper understanding of facet-dependent enzyme-like activities and might advance the development of noble metal-based nanomaterials with both enhanced and targeted antibacterial activities.

Metal–organic framework as an efficient filter for the removal of heavy metal cations in water

Efficient removal of heavy metal cations by a zirconium phosphonate based metal–organic framework (MOF) filter.

Non-canonical targets destabilize microRNAs in human Argonautes

Although much is known about microRNA (miRNA) biogenesis and the mechanism by which miRNAs regulate their targets, little is known about the regulation of miRNA stability. Mature miRNAs are stabilized by binding to Argonaute (Ago) proteins, the core components of the RNA-induced silencing complex (RISC). Recent studies suggest that interactions between miRNAs and their highly complementary target RNAs promote release of miRNAs from Ago proteins, and this in turn can lead to destabilization of miRNAs. However, the physiological triggers of miRNA destabilization with molecular mechanisms remain largely unknown. Here, using an in vitro system that consists of a minimal human Ago2-RISC in HEK293T cell lysates, we sought to understand how miRNAs are destabilized by their targets. Strikingly, we showed that miRNA destabilization is dramatically enhanced by an interaction with seedless, non-canonical targets. We then showed that this process entails not only unloading of miRNAs from Ago, but also 3΄ end destabilization of miRNAs occurred within Ago. Furthermore, our mutation analysis indicates that conformational changes in the hinge region of the Ago PAZ domain are likely to be the main driving force of the miRNA destabilization. Our collective results suggest that non-canonical targets may provide a stability control mechanism in the regulation of miRNAs in humans.

Exploring the RNA-bound and RNA-free human Argonaute-2 by molecular dynamics simulation method

Argonaute 2 (Ago2) protein is the major vehicle of microRNAs (miRNAs)-guided gene repression and silencing processes. Although the crystal structure of human Ago2 (hAgo2) has recently been disclosed, the information of dynamically structural character of protein-RNA recognition is still lacking. Molecular dynamics simulations were used to systematically explore hAgo2 in the presence and absence of RNA duplex. Stable direct and water-mediated hydrogen bonds were observed between guide RNA backbone atoms and hAgo2, especially for nucleotides 2-7. In addition, water-mediated hydrogen bonds are indicated to be critical in the specific recognition between hAgo2 and the conserved adenine in position 1 of target RNA. The core domains (N, PAZ, MID, and PIWI) possess rigid body movements during the simulations. The motions of N-PAZ and PIWI-MID are negatively correlated with or without RNA binding and PAZ domain is identified as the most mobile domain in both systems. The reorientation of PAZ domain not only influences the binding of helix-7 and RNA duplex, the initial pairing process, but also the shape of N-PAZ cleft, where the supplemental base pairing occurs. It is speculated that PAZ domain could be a key regulator in hAgo2-mediated miRNA-induced gene regulation.

Graphene-Induced Pore Formation on Cell Membranes

Examining interactions between nanomaterials and cell membranes can expose underlying mechanisms of nanomaterial cytotoxicity and guide the design of safer nanomedical technologies. Recently, graphene has been shown to exhibit potential toxicity to cells; however, the molecular processes driving its lethal properties have yet to be fully characterized. We here demonstrate that graphene nanosheets (both pristine and oxidized) can produce holes (pores) in the membranes of A549 and Raw264.7 cells, substantially reducing cell viability. Electron micrographs offer clear evidence of pores created on cell membranes. Our molecular dynamics simulations reveal that multiple graphene nanosheets can cooperate to extract large numbers of phospholipids from the membrane bilayer. Strong dispersion interactions between graphene and lipid-tail carbons result in greatly depleted lipid density within confined regions of the membrane, ultimately leading to the formation of water-permeable pores. This cooperative lipid extraction mechanism for membrane perforation represents another distinct process that contributes to the molecular basis of graphene cytotoxicity.

Membrane destruction-mediated antibacterial activity of tungsten disulfide (WS2)

Tungsten disulfide (WS₂) demonstrates clear antibacterial activity through inducing mechanical damage to the bacteria membrane integrity.

Structural influence of proteins upon adsorption to MoS2 nanomaterials: comparison of MoS2 force field parameters

Molybdenum disulfide (MoS₂) has recently emerged as a promising nanomaterial in a wide range of applications due to its unique and impressive properties.

Understanding the core of RNA interference: The dynamic aspects of Argonaute-mediated processes

At the core of RNA interference, the Argonaute proteins (Ago) load and utilize small guide nucleic acids to silence mRNAs or cleave foreign nucleic acids in a sequence specific manner. In recent years, based on extensive structural studies of Ago and its interaction with the nucleic acids, considerable progress has been made to reveal the dynamic aspects of various Ago-mediated processes. Here we review these novel insights into the guide-strand loading, duplex unwinding, and effects of seed mismatch, with a focus on two representative Agos, the human Ago 2 (hAgo2) and the bacterial Thermus thermophilus Ago (TtAgo). In particular, comprehensive molecular simulation studies revealed that although sharing similar overall structures, the two Agos have vastly different conformational landscapes and guide-strand loading mechanisms because of the distinct rigidity of their L1-PAZ hinge. Given the central role of the PAZ motions in regulating the exposure of the nucleic acid binding channel, these findings exemplify the importance of protein motions in distinguishing the overlapping, yet distinct, mechanisms of Ago-mediated processes in different organisms.

Post-transcriptional Regulation of BRCA2 through Interactions with miR-19a and miR-19b

Breast cancer type 2, early onset susceptibility gene (BRCA2) is a major component of the homologous recombination DNA repair pathway. It acts as a tumor suppressor whose function is often lost in cancers. Patients with specific mutations in the BRCA2 gene often display discrete clinical, histopathological, and molecular features. However, a subset of sporadic cancers has wild type BRCA2 and display defects in the homology-directed repair pathway, which is the hallmark of ‘BRCAness.’ The mechanisms by which BRCAness arises are not well understood but post-transcriptional regulation of BRCA2 gene expression by microRNAs (miRNAs) may contribute to this phenotype. Here, we examine the post-transcriptional effects that some members of the six-miRNA cluster known as the miR-17/92 cluster have on the abundance of BRCA2’s messenger RNA (mRNA) and protein. We discuss two interactions involving the miR-19a and miR-19b members of the cluster and the 3′UTR of BRCA2’s mRNA. We investigated these miRNA:mRNA interactions in 15 cell lines derived from pancreatic, breast, colon, and kidney tissue. We show that over-expression of these two miRNAs results in a concomitant decrease of BRCA2’s mRNA and protein expression in a subset of the tested cell lines. Additionally, using luciferase reporter assays we identified direct interactions between miR-19a/miR-19b and a miRNA response element (MRE) in BRCA2’s 3′UTR. Our results suggest that BRCA2 is subject to a complex post-transcriptional regulatory program that has specific dependencies on the genetic and phenotypic background of cell types.

Robust Denaturation of Villin Headpiece by MoS2 Nanosheet: Potential Molecular Origin of the Nanotoxicity

MoS₂ nanosheet, a new two-dimensional transition metal dichalcogenides nanomaterial, has attracted significant attentions lately due to many potential promising biomedical applications. Meanwhile, there is also a growing concern on its biocompatibility, with little known on its interactions with various biomolecules such as proteins. In this study, we use all-atom molecular dynamics simulations to investigate the interaction of a MoS₂ nanosheet with Villin Headpiece (HP35), a model protein widely used in protein folding studies. We find that MoS₂ exhibits robust denaturing capability to HP35, with its secondary structures severely destroyed within hundreds of nanosecond simulations. Both aromatic and basic residues are critical for the protein anchoring onto MoS₂ surface, which then triggers the successive protein unfolding process. The main driving force behind the adsorption process is the dispersion interaction between protein and MoS₂ monolayer. Moreover, water molecules at the interface between some key hydrophobic residues (e.g. Trp-64) and MoS₂ surface also help to accelerate the process driven by nanoscale drying, which provides a strong hydrophobic force. These findings might have shed new light on the potential nanotoxicity of MoS₂ to proteins with atomic details, which should be helpful in guiding future biomedical applications of MoS₂ with its nanotoxicity mitigated.

A Flexible Domain-Domain Hinge Promotes an Induced-fit Dominant Mechanism for the Loading of Guide-DNA into Argonaute Protein in Thermus thermophilus

Argonaute proteins (Ago) are core components of the RNA Induced Silencing Complex (RISC) that load and utilize small guide nucleic acids to silence mRNAs or cleave foreign DNAs. Despite the essential role of Ago in gene regulation and defense against virus, the molecular mechanism of guide-strand loading into Ago remains unclear. We explore such a mechanism in the bacterium Thermus thermophilus Ago (TtAgo), via a computational approach combining molecular dynamics, bias-exchange metadynamics, and protein-DNA docking. We show that apo TtAgo adopts multiple closed states that are unable to accommodate guide-DNA. Conformations able to accommodate the guide are beyond the reach of thermal fluctuations from the closed states. These results suggest an induced-fit dominant mechanism for guide-strand loading in TtAgo, drastically different from the two-step mechanism for human Ago 2 (hAgo2) identified in our previous study. Such a difference between TtAgo and hAgo2 is found to mainly originate from the distinct rigidity of their L1-PAZ hinge. Further comparison among known Ago structures from various species indicates that the L1-PAZ hinge may be flexible in general for prokaryotic Ago's but rigid for eukaryotic Ago's.

A structural view of microRNA-target recognition

It is well established that the correct identification of the messenger RNA targeted by a given microRNA (miRNA) is a difficult problem, and that available methods all suffer from low specificity. We hypothesize that the correct identification of the pairing should take into account the effect of the Argonaute protein (AGO), an essential catalyst of the recognition process. Therefore, we developed a strategy named MiREN for building and scoring three-dimensional models of the ternary complex formed by AGO, a miRNA and 22 nt of a target mRNA that putatively interacts with it. We show here that MiREN can be used to assess the likelihood that an RNA molecule is the target of a given miRNA and that this approach is more accurate than other existing methods, usually based on sequence or sequence-related features. Our results also suggest that AGO plays a relevant role in the selection of the miRNA targets. Our method can represent an additional step for refining predictions made by faster but less accurate classical methods for the identification of miRNA targets.

Non-Specific Blocking of miR-17-5p Guide Strand in Triple Negative Breast Cancer Cells by Amplifying Passenger Strand Activity

Conventional wisdom holds that only one of the two strands in a micro ribonucleic acid (miRNA) precursor duplex is selected as the active miRNA guide strand. The complementary miRNA passenger strand, however, is thought to be inactive. High levels of the oncogenic miRNA (oncomiR) guide strand called miR-17-5p is overexpressed in triple negative breast cancer (TNBC) and can inhibit ribosomal translation of tumor suppressor gene mRNAs, such as programmed cell death 4 (PDCD4) or phosphatase and tensin homolog (PTEN). We hypothesized that knocking down the oncogenic microRNA (oncomiR) miR-17-5p might restore the expression levels of PDCD4 and PTEN tumor suppressor proteins, illustrating a route to oligonucleotide therapy of TNBC. Contrary to conventional wisdom, antisense knockdown of oncomiR miR-17-5p guide strand reduced PDCD4 and PTEN proteins by 1.8±0.3 fold in human TNBC cells instead of raising them. Bioinformatics analysis and folding energy calculations revealed that mRNA targets of miR-17-5p guide strand, such as PDCD4 and PTEN, could also be regulated by miR-17-3p passenger strand. Due to high sequence homology between the antisense molecules and miR-17-3p passenger strand, as well as the excess binding sites for the passenger strand on the 3’UTR of PDCD4 and PTEN mRNAs, introducing a miR-17-3p DNA-LNA mimic to knock down miR-17-5p reduced PDCD4 and PTEN protein expression instead of raising them. Our results imply that therapeutic antisense sequences against miRNAs should be designed to target the miRNA strand with the greatest number of putative binding sites in the target mRNAs, while minimizing affinity for the minor strand.

Dysregulated co-stimulatory molecule expression in a Sjögren's syndrome mouse model with potential implications by microRNA-146a

Sjögren's syndrome (SjS) is an autoimmune condition that primarily affects salivary and lacrimal glands, causing loss of secretion. We have previously shown that microRNA-146a (miR-146a) is over-expressed in the salivary glands and peripheral blood mononuclear cells (PBMC) of SjS-prone mice (C57BL/6.NOD-Aec1Aec2, B6DC) and in PBMC of SjS patients. The purpose of this research was to identify a target molecule of miR-146a and identify subpopulations of cells affected by altered miR-146a in the salivary glands of SjS-prone mice. In silico analyses identified costimulatory molecule CD80 as a potential target of miR-146a. Luciferase assay of the human CD80 3'untranslated region demonstrated miR-146a directly inhibited CD80 protein expression as indicated by reduced luciferase reporter expression and an examination of B6DC salivary glands revealed a reduction in CD80 protein. More interestingly, the specific reduction in CD80 protein was detected from the salivary gland epithelial cell population and in interstitial dendritic cells in the glands as well. The reduction in CD80 protein levels in salivary gland epithelial cells were negatively associated with elevated miR-146a expression. Therefore, this study provides the first indication that salivary gland epithelial cells may be critically involved in SjS progression by altering CD86:CD80 protein ratio in response to miR-146a upregulation.

Beyond the one-locus-one-miRNA paradigm: microRNA isoforms enable deeper insights into breast cancer heterogeneity

Here we describe our study of miRNA isoforms (isomiRs) in breast cancer (BRCA) and normal breast data sets from the Cancer Genome Atlas (TCGA) repository. We report that the full isomiR profiles, from both known and novel human-specific miRNA loci, are particularly rich in information and can distinguish tumor from normal tissue much better than the archetype miRNAs. IsomiR expression is also dependent on the patient's race, exemplified by miR-183-5p, several isomiRs of which are upregulated in triple negative BRCA in white but not black women. Additionally, we find that an isomiR's 5′ endpoint and length, but not the genomic origin, are key determinants of the regulation of its expression. Overexpression of distinct miR-183-5p isomiRs in MDA-MB-231 cells followed by microarray analysis revealed that each isomiR has a distinct impact on the cellular transcriptome. Parallel integrative analysis of mRNA expression from BRCA data sets of the TCGA repository demonstrated that isomiRs can distinguish between the luminal A and luminal B subtypes and explain in more depth the molecular differences between them than the archetype molecules. In conclusion, our findings provide evidence that post-transcriptional studies of BRCA will benefit from transcending the one-locus-one-miRNA paradigm and taking into account all isoforms from each miRNA locus as well as the patient's race.

Markov State Models Reveal a Two-Step Mechanism of miRNA Loading into the Human Argonaute Protein: Selective Binding followed by Structural Re-arrangement

Argonaute (Ago) proteins and microRNAs (miRNAs) are central components in RNA interference, which is a key cellular mechanism for sequence-specific gene silencing. Despite intensive studies, molecular mechanisms of how Ago recognizes miRNA remain largely elusive. In this study, we propose a two-step mechanism for this molecular recognition: selective binding followed by structural re-arrangement. Our model is based on the results of a combination of Markov State Models (MSMs), large-scale protein-RNA docking, and molecular dynamics (MD) simulations. Using MSMs, we identify an open state of apo human Ago-2 in fast equilibrium with partially open and closed states. Conformations in this open state are distinguished by their largely exposed binding grooves that can geometrically accommodate miRNA as indicated in our protein-RNA docking studies. miRNA may then selectively bind to these open conformations. Upon the initial binding, the complex may perform further structural re-arrangement as shown in our MD simulations and eventually reach the stable binary complex structure. Our results provide novel insights in Ago-miRNA recognition mechanisms and our methodology holds great potential to be widely applied in the studies of other important molecular recognition systems.

In RNA interference, Argonaute proteins and microRNAs together form the functional core that regulates the gene expression with high sequence specificity. Elucidating the detailed mechanism of molecular recognition between Argonaute proteins and microRNAs is thus important not only for the fundamental understanding of RNA interference, but also for the further development of microRNA-based therapeutic application. In this work, we propose a two-step model to understand the mechanism of microRNA loading into human Argonaute-2: selective binding followed by structural re-arrangement. Our model is based on the results from a combined approach of molecular dynamics simulations, Markov State Models and protein-RNA docking. In particular, we identify a metastable open state of apo hAgo2 in rapid equilibrium with other states. Some of conformations in this open state have largely exposed RNA binding groove that can accommodate microRNA. We further show that the initial Argonaute-microRNA binding complex undergoes structural re-arrangement to reach stable binary crystal structure. These results provide novel insights into the underlying mechanism of Argonaute-microRNA recognition. In addition, our method is readily applicable to the investigation of other complex molecular recognition events such as protein-protein interactions and protein-ligand binding.

Hepatitis B virus and microRNAs: Complex interactions affecting hepatitis B virus replication and hepatitis B virus-associated diseases

Chronic infection with the hepatitis B virus (HBV) is the leading risk factor for the development of hepatocellular carcinoma (HCC). With nearly 750000 deaths yearly, hepatocellular carcinoma is the second highest cause of cancer-related death in the world. Unfortunately, the molecular mechanisms that contribute to the development of HBV-associated HCC remain incompletely understood. Recently, microRNAs (miRNAs), a family of small non-coding RNAs that play a role primarily in post-transcriptional gene regulation, have been recognized as important regulators of cellular homeostasis, and altered regulation of miRNA expression has been suggested to play a significant role in virus-associated diseases and the development of many cancers. With this in mind, many groups have begun to investigate the relationship between miRNAs and HBV replication and HBV-associated disease. Multiple findings suggest that some miRNAs, such as miR-122, and miR-125 and miR-199 family members, are playing a role in HBV replication and HBV-associated disease, including the development of HBV-associated HCC. In this review, we discuss the current state of our understanding of the relationship between HBV and miRNAs, including how HBV affects cellular miRNAs, how these miRNAs impact HBV replication, and the relationship between HBV-mediated miRNA regulation and HCC development. We also address the impact of challenges in studying HBV, such as the lack of an effective model system for infectivity and a reliance on transformed cell lines, on our understanding of the relationship between HBV and miRNAs, and propose potential applications of miRNA-related techniques that could enhance our understanding of the role miRNAs play in HBV replication and HBV-associated disease, ultimately leading to new therapeutic options and improved patient outcomes.

"Off-Spotter": very fast and exhaustive enumeration of genomic lookalikes for designing CRISPR/Cas guide RNAs

BACKGROUND

CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated nucleases) is a powerful component of the prokaryotic immune system that has been adapted for targeted genetic engineering in higher organisms. A key element of CRISPR/Cas is the "guide" RNA (gRNA) that is ~20 nucleotides (nts) in length and designed to be complementary to the intended target site. An integral requirement of the CRISPR/Cas system is that the target site be followed by a protospacer adjacent motif (PAM). Care needs to be exercised during gRNA design to avoid unintended ("off-target") interactions.

RESULTS

We designed and implemented the Off-Spotter algorithm to assist with the design of optimal gRNAs. When presented with a candidate gRNA sequence and a PAM, Off-Spotter quickly and exhaustively identifies all genomic sites that satisfy the PAM constraint and are identical or nearly-identical to the provided gRNA. Off-Spotter achieves its extreme performance through purely algorithmic means and not through hardware accelerators such as graphical processing units (GPUs). Off-Spotter also allows the user to identify on-the-fly how many and which nucleotides of the gRNA comprise the "seed". Off-Spotter's output includes a histogram showing the number of potential off-targets as a function of the number of mismatches. The output also includes for each potential off-target the site's genomic location, a human genome browser hyperlink to the corresponding location, genomic annotation in the vicinity of the off-target, GC content, etc.

CONCLUSION

Off-Spotter is very fast and flexible and can help in the design of optimal gRNAs by providing several PAM choices, a run-time definition of the seed and of the allowed number of mismatches, and a flexible output interface that allows sorting of the results, optional viewing/hiding of columns, etc. A key element of Off-Spotter is that it does not have a rigid definition of the seed: instead, the user can declare both the seed's location and extent on-the-fly. We expect that this flexibility in combination with Off-Spotter's speed and richly annotated output will enable experimenters to interactively and quickly explore different scenarios and gRNA possibilities.

miR-CLIP capture of a miRNA targetome uncovers a lincRNA H19-miR-106a interaction

Identifying the interaction partners of noncoding RNAs is essential for elucidating their functions. We have developed an approach, termed microRNA crosslinking and immunoprecipitation (miR-CLIP), using pre-miRNAs modified with psoralen and biotin to capture their targets in cells. Photo-crosslinking and Argonaute 2 immunopurification followed by streptavidin affinity purification of probe-linked RNAs provided selectivity in the capture of targets, which were identified by deep sequencing. miR-CLIP with pre-miR-106a, a miR-17-5p family member, identified hundreds of putative targets in HeLa cells, many carrying conserved sequences complementary to the miRNA seed but also many that were not predicted computationally. miR-106a overexpression experiments confirmed that miR-CLIP captured functional targets, including H19, a long noncoding RNA that is expressed during skeletal muscle cell differentiation. We showed that miR-17-5p family members bind H19 in HeLa cells and myoblasts. During myoblast differentiation, levels of H19, miR-17-5p family members and mRNA targets changed in a manner suggesting that H19 acts as a 'sponge' for these miRNAs.

MiR-103a-3p targets the 5' UTR of GPRC5A in pancreatic cells

MicroRNAs (miRNAs) are short noncoding RNAs that regulate the expression of their targets in a sequence-dependent manner. For protein-coding transcripts, miRNAs regulate expression levels through binding sites in either the 3' untranslated region (3' UTR) or the amino acid coding sequence (CDS) of the targeted messenger RNA (mRNA). Currently, for the 5' untranslated region (5' UTR) of mRNAs, very few naturally occurring examples exist whereby the targeting miRNA down-regulates the expression of the corresponding mRNA in a seed-dependent manner. Here we describe and characterize two miR-103a-3p target sites in the 5' UTR of GPRC5A, a gene that acts as a tumor suppressor in some cancer contexts and as an ongocene in other cancer contexts. In particular, we show that the interaction of miR-103a-3p with each of these two 5' UTR targets reduces the expression levels of both GPRC5A mRNA and GPRC5A protein in one normal epithelial and two pancreatic cancer cell lines. By ectopically expressing "sponges" that contain instances of the wild-type 5' UTR targets we also show that we can reduce miR-103a-3p levels and increase GPRC5A mRNA and protein levels. These findings provide some first knowledge on the post-transcriptional regulation of this tumor suppressor/oncogene and present additional evidence for the participation of 5' UTRs in miRNA driven post-transcriptional regulatory control.

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets, including new miRNA target-identification software

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.

Argonaute CLIP-Seq reveals miRNA targetome diversity across tissue types

To date, analyses of individual targets have provided evidence of a miRNA targetome that extends beyond the boundaries of messenger RNAs (mRNAs) and can involve non-Watson-Crick base pairing in the miRNA seed region. Here we report our findings from analyzing 34 Argonaute HITS-CLIP datasets from several human and mouse cell types. Investigation of the architectural (i.e. bulge vs. contiguous pairs) and sequence (Watson-Crick vs. G:U pairs) preferences for human and mouse miRNAs revealed that many heteroduplexes are “non-canonical” i.e. their seed region comprises G:U and bulge combinations. The genomic distribution of miRNA targets differed distinctly across cell types but remained congruent across biological replicates of the same cell type. For some cell types intergenic and intronic targets were more frequent whereas in other cell types mRNA targets prevailed. The findings suggest an expanded model of miRNA targeting that is more frequent than the standard model currently in use. Lastly, our analyses of data from different cell types and laboratories revealed consistent Ago-loaded miRNA profiles across replicates whereas, unexpectedly, the Ago-loaded targets exhibited a much more dynamic behavior across biological replicates.

The miR-17/92 cluster: a comprehensive update on its genomics, genetics, functions and increasingly important and numerous roles in health and disease

The miR-17/92 cluster is among the best-studied microRNA clusters. Interest in the cluster and its members has been increasing steadily and the number of publications has grown exponentially since its discovery with more than 1000 articles published in 2012 alone. Originally found to be involved in tumorigenesis, research work in recent years has uncovered unexpected roles for its members in a wide variety of settings that include normal development, immune diseases, cardiovascular diseases, neurodegenerative diseases and aging. In light of its ever-increasing importance and ever-widening regulatory roles, we review here the latest body of knowledge on the cluster's involvement in health and disease as well as provide a novel perspective on the full spectrum of protein-coding and non-coding transcripts that are likely regulated by its members.

Reduced polyphenol oxidase gene expression and enzymatic browning in potato (Solanum tuberosum L.) with artificial microRNAs

BACKGROUND

Polyphenol oxidase (PPO), often encoded by a multi-gene family, causes oxidative browning, a significant problem in many food products. Low-browning potatoes were produced previously through suppression of PPO gene expression, but the contribution of individual PPO gene isoform to the oxidative browning process was unknown. Here we investigated the contributions of different PPO genes to total PPO protein activity, and the correlations between PPO protein level, PPO activity and tuber tissue browning potential by suppression of all previously characterized potato PPO genes, both individually and in combination using artificial microRNAs (amiRNAs) technology.

RESULTS

Survey of the potato genome database revealed 9 PPO-like gene models, named StuPPO1 to StuPPO9 in this report. StuPPO1, StuPPO2, StuPPO3 and StuPPO4 are allelic to the characterized POTP1/P2, POT32, POT33 and POT72, respectively. Fewer ESTs were found to support the transcriptions of StuPPO5 to StuPPO8. StuPPO9 related ESTs were expressed at significant higher levels in pathogen-infected potato tissues. A series of browning phenotypes were obtained by suppressing StuPPO1 to StuPPO4 genes alone and in combination. Down-regulation of one or several of the PPO genes did not usually cause up-regulation of the other PPO genes in the transgenic potato tubers, but resulted in reduced PPO protein levels. The different PPO genes did not contribute equally to the total PPO protein content in the tuber tissues, with StuPPO2 accounting for ~ 55% as the major contributor, followed by StuPPO1, ~ 25-30% and StuPPO3 and StuPPO4 together with less than 15%. Strongly positive correlations between PPO protein level, PPO activity and browning potential were demonstrated in our analysis. Low PPO activity and low-browning potatoes were produced by simultaneous down-regulation of StuPPO2 to StuPPO4, but the greatest reduction occurred when StuPPO1 to StuPPO4 were all suppressed.

CONCLUSION

StuPPO1 to StuPPO4 genes contributed to browning reactions in tuber tissues but their effect was not equal. Different PPO genes may be regulated independently reflecting their diversified functions. Our results show that amiRNAs can be used to suppress closely related members of highly conserved multi-gene family. This approach also suggests a new strategy for breeding low-browning crops using small DNA inserts.

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.

Stability of miRNA 5'terminal and seed regions is correlated with experimentally observed miRNA-mediated silencing efficacy

MicroRNAs (miRNAs) are key regulators of sequence-specific gene silencing. However, crucial factors that determine the efficacy of miRNA-mediated target gene silencing are poorly understood. Here we mathematized base-pairing stability and showed that miRNAs with an unstable 5′ terminal duplex and stable seed-target duplex exhibit strong silencing activity. The results are consistent with the previous findings that an RNA strand with unstable 5′ terminal in miRNA duplex easily loads onto the RNA-induced silencing complex (RISC), and miRNA recognizes target mRNAs with seed-complementary sequences to direct posttranscriptional repression. Our results suggested that both the unwinding and target recognition processes of miRNAs could be proficiently controlled by the thermodynamics of base-pairing in protein-free condition. Interestingly, such thermodynamic parameters might be evolutionarily well adapted to the body temperatures of various species.