Structural basis for Cas9 off-target activity.

The target DNA specificity of the CRISPR-associated genome editor nuclease Cas9 is determined by complementarity to a 20-nucleotide segment in its guide RNA. However, Cas9 can bind and cleave partially complementary off-target sequences, which raises safety concerns for its use in clinical applications. Here, we report crystallographic structures of Cas9 bound to bona fide off-target substrates, revealing that off-target binding is enabled by a range of noncanonical base-pairing interactions within the guide:off-target heteroduplex. Off-target substrates containing single-nucleotide deletions relative to the guide RNA are accommodated by base skipping or multiple noncanonical base pairs rather than RNA bulge formation. Finally, PAM-distal mismatches result in duplex unpairing and induce a conformational change in the Cas9 REC lobe that perturbs its conformational activation. Together, these insights provide a structural rationale for the off-target activity of Cas9 and contribute to the improved rational design of guide RNAs and off-target prediction algorithms.

Cas12a domain flexibility guides R-loop formation and forces RuvC resetting.

The specific nature of CRISPR-Cas12a makes it a desirable RNA-guided endonuclease for biotechnology and therapeutic applications. To understand how R-loop formation within the compact Cas12a enables target recognition and nuclease activation, we used cryo-electron microscopy to capture wild-type Acidaminococcus sp. Cas12a R-loop intermediates and DNA delivery into the RuvC active site. Stages of Cas12a R-loop formation-starting from a 5-bp seed-are marked by distinct REC domain arrangements. Dramatic domain flexibility limits contacts until nearly complete R-loop formation, when the non-target strand is pulled across the RuvC nuclease and coordinated domain docking promotes efficient cleavage. Next, substantial domain movements enable target strand repositioning into the RuvC active site. Between cleavage events, the RuvC lid conformationally resets to occlude the active site, requiring re-activation. These snapshots build a structural model depicting Cas12a DNA targeting that rationalizes observed specificity and highlights mechanistic comparisons to other class 2 effectors.

Unveiling the RNA-mediated allosteric activation discloses functional hotspots in CRISPR-Cas13a.

Cas13a is a recent addition to the CRISPR-Cas toolkit that exclusively targets RNA, which makes it a promising tool for RNA detection. It utilizes a CRISPR RNA (crRNA) to target RNA sequences and trigger a composite active site formed by two 'Higher Eukaryotes and Prokaryotes Nucleotide' (HEPN) domains, cleaving any solvent-exposed RNA. In this system, an intriguing form of allosteric communication controls the RNA cleavage activity, yet its molecular details are unknown. Here, multiple-microsecond molecular dynamics simulations are combined with graph theory to decipher this intricate activation mechanism. We show that the binding of a target RNA acts as an allosteric effector, by amplifying the communication signals over the dynamical noise through interactions of the crRNA at the buried HEPN1-2 interface. By introducing a novel Signal-to-Noise Ratio (SNR) of communication efficiency, we reveal critical allosteric residues-R377, N378, and R973-that rearrange their interactions upon target RNA binding. Alanine mutation of these residues is shown to select target RNA over an extended complementary sequence beyond guide-target duplex for RNA cleavage, establishing the functional significance of these hotspots. Collectively our findings offer a fundamental understanding of the Cas13a mechanism of action and pave new avenues for the development of highly selective RNA-based cleavage and detection tools.

New design strategies for ultra-specific CRISPR-Cas13a-based RNA detection with single-nucleotide mismatch sensitivity.

An increasingly pressing need for clinical diagnostics has required the development of novel nucleic acid-based detection technologies that are sensitive, fast, and inexpensive, and that can be deployed at point-of-care. Recently, the RNA-guided ribonuclease CRISPR-Cas13 has been successfully harnessed for such purposes. However, developing assays for detection of genetic variability, for example single-nucleotide polymorphisms, is still challenging and previously described design strategies are not always generalizable. Here, we expanded our characterization of LbuCas13a RNA-detection specificity by performing a combination of experimental RNA mismatch tolerance profiling, molecular dynamics simulations, protein, and crRNA engineering. We found certain positions in the crRNA-target-RNA duplex that are particularly sensitive to mismatches and establish the effect of RNA concentration in mismatch tolerance. Additionally, we determined that shortening the crRNA spacer or modifying the direct repeat of the crRNA leads to stricter specificities. Furthermore, we harnessed our understanding of LbuCas13a allosteric activation pathways through molecular dynamics and structure-guided engineering to develop novel Cas13a variants that display increased sensitivities to single-nucleotide mismatches. We deployed these Cas13a variants and crRNA design strategies to achieve superior discrimination of SARS-CoV-2 strains compared to wild-type LbuCas13a. Together, our work provides new design criteria and Cas13a variants to use in future easier-to-implement Cas13-based RNA detection applications.

Enhancement of CRISPR/Cas12a trans-cleavage activity using hairpin DNA reporters.

The RNA programmed non-specific (trans) nuclease activity of CRISPR-Cas Type V and VI systems has opened a new era in the field of nucleic acid-based detection. Here, we report on the enhancement of trans-cleavage activity of Cas12a enzymes using hairpin DNA sequences as FRET-based reporters. We discover faster rate of trans-cleavage activity of Cas12a due to its improved affinity (Km) for hairpin DNA structures, and provide mechanistic insights of our findings through Molecular Dynamics simulations. Using hairpin DNA probes we significantly enhance FRET-based signal transduction compared to the widely used linear single stranded DNA reporters. Our signal transduction enables faster detection of clinically relevant double stranded DNA targets with improved sensitivity and specificity either in the presence or in the absence of an upstream pre-amplification step.

Substrate-independent activation pathways of the CRISPR-Cas9 HNH nuclease.

A hallmark of tightly regulated high-fidelity enzymes is that they become activated only after encountering cognate substrates, often by an induced-fit mechanism rather than conformational selection. Upon analysis of molecular dynamics trajectories, we recently discovered that the Cas9 HNH domain exists in three conformations: 1) Y836 (which is two residues away from the catalytic D839 and H840 residues) is hydrogen bonded to the D829 backbone amide, 2) Y836 is hydrogen bonded to the backbone amide of D861 (which is one residue away from the third catalytic residue N863), and 3) Y836 is not hydrogen bonded to either residue. Each of the three conformers differs from the active state of HNH. The conversion between the inactive and active states involves a local unfolding-refolding process that displaces the Cα and side chain of the catalytic N863 residue by ∼5 Å and ∼10 Å, respectively. In this study, we report the two largest principal components of coordinate variance of the HNH domain throughout molecular dynamics trajectories to establish the interconversion pathways of these conformations. We show that conformation 2 is an obligate step between conformations 1 and 3, which are not directly interconvertible without conformation 2. The loss of hydrogen bonding of the Y836 side chain in conformation 3 likely plays an essential role in activation during local unfolding-refolding of an α-helix containing the catalytic N863. Three single Lys-to-Ala mutants appear to eliminate this substrate-independent activation pathway of the wild-type HNH nuclease, thereby enhancing the fidelity of HNH cleavage.

Acetylcholine-Binding Protein Affinity Profiling of Neurotoxins in Snake Venoms with Parallel Toxin Identification.

Snakebite is considered a concerning issue and a neglected tropical disease. Three-finger toxins (3FTxs) in snake venoms primarily cause neurotoxic effects since they have high affinity for nicotinic acetylcholine receptors (nAChRs). Their small molecular size makes 3FTxs weakly immunogenic and therefore not appropriately targeted by current antivenoms. This study aims at presenting and applying an analytical method for investigating the therapeutic potential of the acetylcholine-binding protein (AChBP), an efficient nAChR mimic that can capture 3FTxs, for alternative treatment of elapid snakebites. In this analytical methodology, snake venom toxins were separated and characterised using high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS) and high-throughput venomics. By subsequent nanofractionation analytics, binding profiling of toxins to the AChBP was achieved with a post-column plate reader-based fluorescence-enhancement ligand displacement bioassay. The integrated method was established and applied to profiling venoms of six elapid snakes (Naja mossambica, Ophiophagus hannah, Dendroaspis polylepis, Naja kaouthia, Naja haje and Bungarus multicinctus). The methodology demonstrated that the AChBP is able to effectively bind long-chain 3FTxs with relatively high affinity, but has low or no binding affinity towards short-chain 3FTxs, and as such provides an efficient analytical platform to investigate binding affinity of 3FTxs to the AChBP and mutants thereof and to rapidly identify bound toxins.

GPU-Enhanced DFTB Metadynamics for Efficiently Predicting Free Energies of Biochemical Systems.

Metadynamics calculations of large chemical systems with ab initio methods are computationally prohibitive due to the extensive sampling required to simulate the large degrees of freedom in these systems. To address this computational bottleneck, we utilized a GPU-enhanced density functional tight binding (DFTB) approach on a massively parallelized cloud computing platform to efficiently calculate the thermodynamics and metadynamics of biochemical systems. To first validate our approach, we calculated the free-energy surfaces of alanine dipeptide and showed that our GPU-enhanced DFTB calculations qualitatively agree with computationally-intensive hybrid DFT benchmarks, whereas classical force fields give significant errors. Most importantly, we show that our GPU-accelerated DFTB calculations are significantly faster than previous approaches by up to two orders of magnitude. To further extend our GPU-enhanced DFTB approach, we also carried out a 10 ns metadynamics simulation of remdesivir, which is prohibitively out of reach for routine DFT-based metadynamics calculations. We find that the free-energy surfaces of remdesivir obtained from DFTB and classical force fields differ significantly, where the latter overestimates the internal energy contribution of high free-energy states. Taken together, our benchmark tests, analyses, and extensions to large biochemical systems highlight the use of GPU-enhanced DFTB simulations for efficiently predicting the free-energy surfaces/thermodynamics of large biochemical systems.

Structural Basis for Reduced Dynamics of Three Engineered HNH Endonuclease Lys-to-Ala Mutants for the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Associated 9 (CRISPR/Cas9) Enzyme.

Many bacteria possess type-II immunity against invading phages or plasmids known as the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) system to detect and degrade the foreign DNA sequences. The Cas9 protein has two endonucleases responsible for double-strand breaks (the HNH domain for cleaving the target strand of DNA duplexes and RuvC domain for the nontarget strand, respectively) and a single-guide RNA-binding domain where the RNA and target DNA strands are base-paired. Three engineered single Lys-to-Ala HNH mutants (K810A, K848A, and K855A) exhibit an enhanced substrate specificity for cleavage of the target DNA strand. We report in this study that in the wild-type (wt) enzyme, D835, Y836, and D837 within the Y836-containing loop (comprising E827-D837) adjacent to the catalytic site have uncharacterizable broadened 1H15N nuclear magnetic resonance (NMR) features, whereas remaining residues in the loop have different extents of broadened NMR spectra. We find that this loop in the wt enzyme exhibits three distinct conformations over the duration of the molecular dynamics simulations, whereas the three Lys-to-Ala mutants retain only one conformation. The versatility of multiple alternate conformations of this loop in the wt enzyme could help to recruit noncognate DNA substrates into the HNH active site for cleavage, thereby reducing its substrate specificity relative to the three mutants. Our study provides further experimental and computational evidence that Lys-to-Ala substitutions reduce dynamics of proteins and thus increase their stability.

Structural and dynamic insights into the HNH nuclease of divergent Cas9 species.

CRISPR-Cas9 is a widely used biochemical tool with applications in molecular biology and precision medicine. The RNA-guided Cas9 protein uses its HNH endonuclease domain to cleave the DNA strand complementary to its endogenous guide RNA. In this study, novel constructs of HNH from two divergent organisms, G. stearothermophilus (GeoHNH) and S. pyogenes (SpHNH) were engineered from their respective full-length Cas9 proteins. Despite low sequence similarity, the X-ray crystal structures of these constructs reveal that the core of HNH surrounding the active site is conserved. Structure prediction of the full-length GeoCas9 protein using Phyre2 and AlphaFold2 also showed that the crystallographic construct of GeoHNH represents the structure of the domain within the full-length GeoCas9 protein. However, significant differences are observed in the solution dynamics of structurally conserved regions of GeoHNH and SpHNH, the latter of which was shown to use such molecular motions to propagate the DNA cleavage signal. Indeed, molecular simulations show that the intradomain signaling pathways, which drive SpHNH function, are non-specific and poorly formed in GeoHNH. Taken together, these outcomes suggest mechanistic differences between mesophilic and thermophilic Cas9 species.

Disruption of electrostatic contacts in the HNH nuclease from a thermophilic Cas9 rewires allosteric motions and enhances high-temperature DNA cleavage.

Allosteric signaling within multidomain proteins is a driver of communication between spatially distant functional sites. Understanding the mechanism of allosteric coupling in large multidomain proteins is the most promising route to achieving spatial and temporal control of the system. The recent explosion of CRISPR-Cas9 applications in molecular biology and medicine has created a need to understand how the atomic level protein dynamics of Cas9, which are the driving force of its allosteric crosstalk, influence its biophysical characteristics. In this study, we used a synergistic approach of nuclear magnetic resonance (NMR) and computation to pinpoint an allosteric hotspot in the HNH domain of the thermostable GeoCas9. We show that mutation of K597 to alanine disrupts a salt-bridge network, which in turn alters the structure, the timescale of allosteric motions, and the thermostability of the GeoHNH domain. This homologous lysine-to-alanine mutation in the extensively studied mesophilic S. pyogenes Cas9 similarly alters the dynamics of the SpHNH domain. We have previously demonstrated that the alteration of allostery via mutations is a source for the specificity enhancement of SpCas9 (eSpCas9). Hence, this may also be true in GeoCas9.

All-atom simulations disentangle the functional dynamics underlying gene maturation in the intron lariat spliceosome.

The spliceosome (SPL) is a majestic macromolecular machinery composed of five small nuclear RNAs and hundreds of proteins. SPL removes noncoding introns from precursor messenger RNAs (pre-mRNAs) and ligates coding exons, giving rise to functional mRNAs. Building on the first SPL structure solved at near-atomic-level resolution, here we elucidate the functional dynamics of the intron lariat spliceosome (ILS) complex through multi-microsecond-long molecular-dynamics simulations of ∼1,000,000 atoms models. The ILS essential dynamics unveils (i) the leading role of the Spp42 protein, which heads the gene maturation by tuning the motions of distinct SPL components, and (ii) the critical participation of the Cwf19 protein in displacing the intron lariat/U2 branch helix. These findings provide unprecedented details on the SPL functional dynamics, thus contributing to move a step forward toward a thorough understanding of eukaryotic pre-mRNA splicing.

Folding-upon-Repair DNA Nanoswitches for Monitoring the Activity of DNA Repair Enzymes.

We present a new class of DNA-based nanoswitches that, upon enzymatic repair, could undergo a conformational change mechanism leading to a change in fluorescent signal. Such folding-upon-repair DNA nanoswitches are synthetic DNA sequences containing O6 -methyl-guanine (O6 -MeG) nucleobases and labelled with a fluorophore/quencher optical pair. The nanoswitches are rationally designed so that only upon enzymatic demethylation of the O6 -MeG nucleobases they can form stable intramolecular Hoogsteen interactions and fold into an optically active triplex DNA structure. We have first characterized the folding mechanism induced by the enzymatic repair activity through fluorescent experiments and Molecular Dynamics simulations. We then demonstrated that the folding-upon-repair DNA nanoswitches are suitable and specific substrates for different methyltransferase enzymes including the human homologue (hMGMT) and they allow the screening of novel potential methyltransferase inhibitors.

Key role of the REC lobe during CRISPR-Cas9 activation by 'sensing', 'regulating', and 'locking' the catalytic HNH domain.

Understanding the conformational dynamics of CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 is of the utmost importance for improving its genome editing capability. Here, molecular dynamics simulations performed using Anton-2 - a specialized supercomputer capturing micro-to-millisecond biophysical events in real time and at atomic-level resolution - reveal the activation process of the endonuclease Cas9 toward DNA cleavage. Over the unbiased simulation, we observe that the spontaneous approach of the catalytic domain HNH to the DNA cleavage site is accompanied by a remarkable structural remodeling of the recognition (REC) lobe, which exerts a key role for DNA cleavage. Specifically, the significant conformational changes and the collective conformational dynamics of the REC lobe indicate a mechanism by which the REC1-3 regions 'sense' nucleic acids, 'regulate' the HNH conformational transition, and ultimately 'lock' the HNH domain at the cleavage site, contributing to its catalytic competence. By integrating additional independent simulations and existing experimental data, we provide a solid validation of the activated HNH conformation, which had been so far poorly characterized, and we deliver a comprehensive understanding of the role of REC1-3 in the activation process. Considering the importance of the REC lobe in the specificity of Cas9, this study poses the basis for fully understanding how the REC components control the cleavage of off-target sequences, laying the foundation for future engineering efforts toward improved genome editing.

Mapping the Structural and Dynamic Determinants of pH-Sensitive Heparin Binding to Granulocyte Macrophage Colony Stimulating Factor.

Granulocyte macrophage colony stimulating factor (GMCSF) is an immunomodulatory cytokine that is harnessed as a therapeutic. GMCSF is known to interact with other clinically important molecules, such as heparin, suggesting that endogenous and administered GMCSF has the potential to modulate orthogonal treatment outcomes. Thus, molecular level characterization of GMCSF and its interactions with biologically active compounds is critical to understanding these mechanisms and predicting clinical consequences. Here, we dissect the biophysical factors that facilitate the GMCSF-heparin interaction, previously shown to be pH-dependent, using nuclear magnetic resonance spectroscopy, surface plasmon resonance, and molecular dynamics simulations. We find that the affinity of GMCSF for heparin increases not only with a transition to acidic pH but also with an increase in heparin chain length. Changes in local flexibility, including a disruption of the N-terminal helix at acidic pH, also accompany the binding of heparin to GMCSF. We use molecular dynamics simulations to propose a mechanism in which a positive binding pocket that is not fully solvent accessible at neutral pH becomes more accessible at acidic pH, facilitating the binding of heparin to the protein.

Allosteric Motions of the CRISPR-Cas9 HNH Nuclease Probed by NMR and Molecular Dynamics.

CRISPR-Cas9 is a widely employed genome-editing tool with functionality reliant on the ability of the Cas9 endonuclease to introduce site-specific breaks in double-stranded DNA. In this system, an intriguing allosteric communication has been suggested to control its DNA cleavage activity through flexibility of the catalytic HNH domain. Here, solution NMR experiments and a novel Gaussian-accelerated molecular dynamics (GaMD) simulation method are used to capture the structural and dynamic determinants of allosteric signaling within the HNH domain. We reveal the existence of a millisecond time scale dynamic pathway that spans HNH from the region interfacing the adjacent RuvC nuclease and propagates up to the DNA recognition lobe in full-length CRISPR-Cas9. These findings reveal a potential route of signal transduction within the CRISPR-Cas9 HNH nuclease, advancing our understanding of the allosteric pathway of activation. Further, considering the role of allosteric signaling in the specificity of CRISPR-Cas9, this work poses the mechanistic basis for novel engineering efforts aimed at improving its genome-editing capability.

Structure and Dynamics of the Acyl Chains in the Membrane Trafficking and Enzymatic Processing of Lipids.

The regulatory chemical mechanisms of lipid trafficking and degradation are involved in many pathophysiological processes, being implicated in severe pain, inflammation, and cancer. In addition, the processing of lipids is also relevant for industrial and environmental applications. However, there is poor understanding of the chemical features that control lipid membrane trafficking and allow lipid-degrading enzymes to efficiently select and hydrolyze specific fatty acids from a complex cellular milieu of bioactive lipids. This is particularly true for lipid acyl chains, which have diverse structures that can critically affect the many complex reactions needed to elongate, desaturate, or transport fatty acids. Building upon our own contributions in this field, we will discuss how molecular simulations, integrated with experimental evidence, have revealed that the structure and dynamics of the lipid tail are actively involved in modulating membrane trafficking at cellular organelles, and enzymatic reactions at cell membranes. Further evidence comes from recent crystal structures of lipid receptors and remodeling enzymes. Taken together, these recent works have identified those structural features of the lipid acyl chain that are crucial for the regioselectivity and stereospecificity of essential desaturation reactions. In this context, we will first illustrate how atomistic and coarse-grained simulations have elucidated the structure-function relationships between the chemical composition of the lipid's acyl chains and the molecular properties of lipid bilayers. Particular emphasis will be given to the prominent chemical role of the number of double carbon-carbon bonds along the lipid acyl chain, that is, discriminating between saturated, monounsaturated, and polyunsaturated lipids. Different levels of saturation in fatty acid molecules dramatically influence the biophysical properties of lipid assemblies and their interaction with proteins. We will then discuss the processing of lipids by membrane-bound enzymes. Our focus will be on lipids such as anandamide and 2-arachidonoylglycerol. These are the main molecules that act as neurotransmitters in the endocannabinoid system. Specifically, recent findings indicate a crucial interplay between the level of saturation of the lipid tail, its energetically and sterically favored conformations, and the hydrophobic accessory cavities in lipid-degrading enzymes, which help form catalytically active conformations of the selected substrate. This Account will emphasize how the specific chemical structure of acyl chains affects the molecular mechanisms for modulating membrane trafficking and selective hydrolysis. The results examined here show that, by using molecular simulations to investigate lipid plasticity and substrate flexibility, researchers can enrich their interpretation of experimental results about the structure-function relationships of lipids. This could positively impact chemical and biological studies in the field and ultimately support protein engineering studies and structure-based drug discovery to target lipid-processing enzymes.

Molecular Dynamics Reveals a DNA-Induced Dynamic Switch Triggering Activation of CRISPR-Cas12a.

CRISPR-Cas12a is a genome-editing system, recently also harnessed for nucleic acid detection, which is promising for the diagnosis of the SARS-CoV-2 coronavirus through the DETECTR technology. Here, a collective ensemble of multimicrosecond molecular dynamics characterizes the key dynamic determinants allowing nucleic acid processing in CRISPR-Cas12a. We show that DNA binding induces a switch in the conformational dynamics of Cas12a, which results in the activation of the peripheral REC2 and Nuc domains to enable cleavage of nucleic acids. The simulations reveal that large-amplitude motions of the Nuc domain could favor the conformational activation of the system toward DNA cleavages. In this process, the REC lobe plays a critical role. Accordingly, the joint dynamics of REC and Nuc shows the tendency to prime the conformational transition of the DNA target strand toward the catalytic site. Most notably, the highly coupled dynamics of the REC2 region and Nuc domain suggests that REC2 could act as a regulator of the Nuc function, similar to what was observed previously for the HNH domain in the CRISPR-associated nuclease Cas9. These mutual domain dynamics could be critical for the nonspecific binding of DNA and thereby for the underlying mechanistic functioning of the DETECTR technology. Considering that REC is a key determinant in the system's specificity, our findings provide a rational basis for future biophysical studies aimed at characterizing its function in CRISPR-Cas12a. Overall, our outcomes advance our mechanistic understanding of CRISPR-Cas12a and provide grounds for novel engineering efforts to improve genome editing and viral detection.

CRISPR-Cas9 conformational activation as elucidated from enhanced molecular simulations.

CRISPR-Cas9 has become a facile genome editing technology, yet the structural and mechanistic features underlying its function are unclear. Here, we perform extensive molecular simulations in an enhanced sampling regime, using a Gaussian-accelerated molecular dynamics (GaMD) methodology, which probes displacements over hundreds of microseconds to milliseconds, to reveal the conformational dynamics of the endonuclease Cas9 during its activation toward catalysis. We disclose the conformational transition of Cas9 from its apo form to the RNA-bound form, suggesting a mechanism for RNA recruitment in which the domain relocations cause the formation of a positively charged cavity for nucleic acid binding. GaMD also reveals the conformation of a catalytically competent Cas9, which is prone for catalysis and whose experimental characterization is still limited. We show that, upon DNA binding, the conformational dynamics of the HNH domain triggers the formation of the active state, explaining how the HNH domain exerts a conformational control domain over DNA cleavage [Sternberg SH et al. (2015) Nature, 527, 110-113]. These results provide atomic-level information on the molecular mechanism of CRISPR-Cas9 that will inspire future experimental investigations aimed at fully clarifying the biophysics of this unique genome editing machinery and at developing new tools for nucleic acid manipulation based on CRISPR-Cas9.

Understanding the mechanistic basis of non-coding RNA through molecular dynamics simulations.

Noncoding RNA (ncRNA) has a key role in regulating gene expression, mediating fundamental processes and diseases via a variety of yet unknown mechanisms. Here, we review recent applications of conventional and enhanced Molecular Dynamics (MD) simulations methods to address the mechanistic function of large biomolecular systems that are tightly involved in the ncRNA function and that are of key importance in life sciences. This compendium focuses of three biomolecular systems, namely the CRISPR-Cas9 genome editing machinery, group II intron ribozyme and the ribonucleoprotein complex of the spliceosome, which edit and process ncRNA. We show how the application of a novel accelerated MD simulations method has been key in disclosing the conformational transitions underlying RNA binding in the CRISPR-Cas9 complex, suggesting a mechanism for RNA recruitment and clarifying the conformational changes required for attaining genome editing. As well, we discuss the use of mixed quantum-classical MD simulations in deciphering the catalytic mechanism of RNA splicing as operated by group II intron ribozyme, one of the largest ncRNA structures crystallized so far. Finally, we debate the future challenges and opportunities in the field, discussing the recent application of MD simulations for unraveling the functional biophysics of the spliceosome, a multi-mega Dalton complex of proteins and small nuclear RNAs that performs RNA splicing in humans. This showcase of applications highlights the current talent of MD simulations to dissect atomic-level details of complex biomolecular systems instrumental for the design of finely engineered genome editing machines. As well, this review aims at inspiring future investigations of several other ncRNA regulatory systems, such as micro and small interfering RNAs, which achieve their function and specificity using RNA-based recognition and targeting strategies.