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

A handheld isothermal fluorescence detector for duplex visualization of aquatic pathogens via enhanced one-pot LAMP-PfAgo assay

The expansion of large-scale aquaculture has exacerbated the challenge of aquatic diseases, resulting in substantial economic losses annually. Currently, traditional laboratory-based diagnostic methods are time-consuming and costly, hindering on-site testing for individual farmers. We address this issue by developing a state-of-the-art handheld isothermal nucleic acid amplification device (WeD-1) capable of fluorescence tracking of reactions and integrating it with an enhanced one-pot Prokaryotic Argonaute based nucleic acid detection method, enabling duplex visual detection of aquatic pathogens. WeD-1 is portable, reusable, user-friendly, and cost-effective, offering real-time smartphone interaction and enabling real-time fluorescence observation during the reaction. The enhanced one-pot Loop-Mediated Isothermal Amplification (LAMP)-PfAgo method, incorporating paraffin-encapsulated lyophilized PfAgo protein, achieves precise target-specific cleavage, significantly enhancing multiplex nucleic acid detection. This innovation streamlines on-site testing, negating the need for specialized laboratory conditions while ensuring an aerosol-free system. With newly developed and highly sensitive LAMP primer sets, our compact WeD-1/LAMP-PfAgo nucleic acid rapid testing system exhibits remarkable sensitivity, readily detecting aquatic pathogens with naked eyes from rapidly prepared fish and shrimp samples within 40 min, even when the Ct values are as high as 34.

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.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.

Argonaute-integrated isothermal amplification for rapid, portable, multiplex detection of SARS-CoV-2 and influenza viruses

Isothermal amplification methods are a promising trend in virus detection because of their superiority in rapidity and sensitivity. However, the generation of false positives and limited multiplexity are major bottlenecks that must be addressed. In this study, we developed a multiplex Argonaute (Ago)-based nucleic acid detection system (MULAN) that integrates rapid isothermal amplification with the multiplex inclusiveness of a single Ago for simultaneous detection of multiple targets such as SARS-CoV-2 and influenza viruses. Owing to its high specificity, MULAN can distinguish targets at a single-base resolution for mutant genotyping. Moreover, MULAN also supports portable and visible devices with a limit of detection of five copies per reaction. Validated by SARS-CoV-2 pseudoviruses and clinical samples of influenza viruses, MULAN showed 100% agreement with quantitative reverse-transcription PCR. These results demonstrated that MULAN has great potential to facilitate reliable, easy, and quick point-of-care diagnosis for promoting the control of infectious diseases.

Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors

Bromodomain and extra-terminal domain (BET) subfamily is the most studied subfamily of bromodomain-containing proteins (BCPs) family which can modulate acetylation signal transduction and produce diverse physiological functions. Thus, the BET family can be treated as an alternative strategy for targeting androgen-receptor (AR)-driven cancers. In order to explore the effect of inhibitors binding to BRD4 (the most studied member of BET family), four 150 ns molecular dynamic simulations were performed (free BRD4, Cpd4-BRD4, Cpd9-BRD4 and Cpd19-BRD4). Docking studies showed that Cpd9 and Cpd19 were located at the active pocket, as well as Cpd4. Molecular dynamics (MD) simulations indicated that only Cpd19 binding to BRD4 can induce residue Trp81-Ala89 partly become α-helix during MD simulations. MM-GBSA calculations suggested that Cpd19 had the best binding effect with BRD4 followed by Cpd4 and Cpd9. Computational alanine scanning results indicated that mutations in Phe83 made the greatest effects in Cpd9-BRD4 and Cpd19-BRD4 complexes, showing that Phe83 may play crucial roles in Cpd9 and Cpd19 binding to BRD4. Our results can provide some useful clues for further BCPs family search.

Binding-induced functional-domain motions in the Argonaute characterized by adaptive advanced sampling

Argonaute proteins in combination with short microRNA (miRNAs) can target mRNA molecules for translation inhibition or degradation and play a key role in many regulatory processes. The miRNAs act as guide RNAs that associate with Argonaute and the complementary mRNA target region. The complex formation results in activation of Argonaute and specific cleavage of the target mRNA. Both the binding and activation processes involve essential domain rearrangements of functional importance. For the Thermus Thermophilus Argonaute (TtAgo) system guide-bound (binary) and guide/target-bound (ternary) complexes are known but how the binding of guide and target mediate domain movements is still not understood. We have studied the Argonaute domain motion in apo and guide/target bound states using Molecular Dynamics simulations and a Hamiltonian replica exchange (H-REMD) method that employs a specific biasing potential to accelerate domain motions. The H-REMD technique indicates sampling of a much broader distribution of domain arrangements both in the apo as well as binary and ternary complexes compared to regular MD simulations. In the apo state domain arrangements corresponding to more compact (closed) states are mainly sampled which undergo an opening upon guide and guide/target binding. Whereas only limited overlap in domain geometry between apo and bound states was found, a larger similarity in the domain distribution is observed for the simulations of binary and ternary complexes. Comparative simulations on ternary complexes with 15 or 16 base pairs (bp) formed between guide and target strands (instead of 14) resulted in dissociation of the 3’-guide strand from the PAZ domain and domain rearrangement. This agrees with the experimental observation that guide-target pairing beyond 14 bps is required for activation and gives a mechanistic explanation for the experimentally observed activation process.

Post-transcriptional gene silencing is an important process to regulate protein synthesis in eukaryotes and prokaryotes. The Argonaute proteins as part of the RNA-induced-silencing-complex (RISC) form a central element of the process by silencing of a target messenger RNA (mRNA) via degradation or repression of translation. The Argonaute protein binds initially a short RNA that acts as a guide to promote binding of a complementary target mRNA. The complex formation can lead to activation of Argonaute and specific cleavage of the target mRNA. The whole process involves domain rearrangements that are not fully understood. We applied an advanced Molecular Dynamics sampling technique to specifically accelerate domain motions of the Thermus Thermophilus Argonaute (TtAgo) system in apo, guide bound and guide/target bound states. The simulations indicate only limited overlap of domain arrangements in apo and bound states and identified domain opening motions necessary for guide and target binding. The study also offers an explanation why a minimum of 15 or 16 base pairs between guide and target strands are necessary for Argonaute activation.

Expression and Functional Analysis of the Argonaute Protein of Thermus thermophilus (TtAgo) in E. coli BL21(DE3)

The prokaryotic Argonaute proteins (pAgos) have been reported to cleave or interfere with DNA targets in a guide-dependent or independent manner. It is often difficult to characterize pAgos in vivo due to the extreme environments favored by their hosts. In the present study, we expressed functional Thermus thermophilus pAgo (TtAgo) in E. coli BL21 (DE3) cells at 37 °C. Initial attempts to express TtAgo in BL21(DE3) cells at 37 °C failed. This was not because of TtAgo mediated general toxicity to the host cells, but instead because of TtAgo-induced loss of its expression plasmid. We employed this discovery to establish a screening system for isolating loss-of-function mutants of TtAgo. The E. colifabI gene was used to help select for full-length TtAgo loss of function mutants, as overexpression of fabI renders the cell to be resistant to the triclosan. We isolated and characterized eight mutations in TtAgo that abrogated function. The ability of TtAgo to induce loss of its expression vector in vivo at 37 °C is an unreported function that is mechanistically different from its reported in vitro activity. These results shed light on the mechanisms by which TtAgo functions as a defense against foreign DNA invasion.

How Different Substitution Positions of F, Cl Atoms in Benzene Ring of 5-Methylpyrimidine Pyridine Derivatives Affect the Inhibition Ability of EGFRL858R/T790M/C797S Inhibitors: A Molecular Dynamics Simulation Study

Lung cancer is the most frequent cause of cancer-related deaths worldwide, and mutations in the kinase domain of the epidermal growth factor receptor (EGFR) are a common cause of non-small-cell lung cancers, which is a major subtype of lung cancers. Recently, a series of 5-methylpyrimidine-pyridinone derivatives have been designed and synthesized as novel selective inhibitors of EGFR and EGFR mutants. However, the binding-based inhibition mechanism has not yet been determined. In this study, we carried out molecular dynamic simulations and free-energy calculations for EGFR derivatives to fill this gap. Based on the investigation, the three factors that influence the inhibitory effect of inhibitors are as follows: (1) The substitution site of the Cl atom is the main factor influencing the activity through steric effect; (2) The secondary factors are repulsion between the F atom (present in the inhibitor) and Glu762, and the blocking effect of Lys745 on the phenyl ring of the inhibitor. (3) The two factors function synergistically to influence the inhibitory capacity of the inhibitor. The theoretical results of this study can provide further insights that will aid the design of oncogenic EGFR inhibitors with high selectivity.

Gene Silencing Mechanisms Revealed by Dynamics of Guide, Target, and Duplex Binding to Argonaute

Argonaute (Ago) protein plays a central role in silencing gene expression by binding a "guide" strand to the base-pair with a complementary mRNA and degrading the mRNA. The current understanding of how Ago-guide and Ago-guide-mRNA complexes assemble is based mainly on static crystal structures; the associated kinetic pathways remain unknown/unclear. By simulating the successive binding of guide/target strand to Thermus thermophilus Ago (TtAgo) and computing the respective free energy landscapes, we directly visualize how TtAgo silencing complexes form and function. We show that the guide binding rate depends on its initial loading position onto TtAgo. Subsequent target recognition beyond the scissile 10-11 nucleotides must overcome a substantial energy barrier for TtAgo's nucleotide-binding groove to expand widely. This work reveals novel roles for the core TtAgo domains and shows how the kinetic barriers that must be overcome for critical structural changes to occur lead to target repression/cleavage.