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.

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.

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.

Making it Rain: Cloud-Based Molecular Simulations for Everyone.

We present a user-friendly front-end for running molecular dynamics (MD) simulations using the OpenMM toolkit on the Google Colab framework. Our goals are (1) to highlight the usage of a cloud-computing scheme for educational purposes for a hands-on approach when learning MD simulations and (2) to exemplify how low-income research groups can perform MD simulations in the microsecond time scale. We hope this work facilitates teaching and learning of molecular simulation throughout the community.

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.

The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains.

Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc3-Man9-GlcNAc2-PP-Dolichol), bacterial (Glc1-GalNAc5-Bac1-PP-Undecaprenol), and archaeal (Glc1-Man1-Gal1-Man1-Glc1-Gal1-Glc1-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.

In silico Investigation of the PglB Active Site Reveals Transient Catalytic States and Octahedral Metal Ion Coordination.

The last step of the bacterial N-glycosylation pathway involves PglB, an oligosaccharyltransferase, which is responsible for the en bloc transfer of a fully assembled oligosaccharide chain to a protein possessing the extended motif D/E-X-N-X-S/T. Recently, this molecule had its full structure elucidated, enabling the description of its domains and the proposition of a catalytic mechanism. By employing molecular dynamics simulations, we were able to evaluate structural aspects of PglB, suggesting prevalent motions that may bring insights into the mechanism of the glycosylated peptide detachment. Additionally, we identified transient states at the catalytic site, in which the previously described carboxamide twisting mechanism was observed. Aided by quantum mechanics calculations for each different conformational states of the catalytic site, we determined the presence of an octahedral metal coordination, along with the presence of one water molecule at the catalytic site.

Natural Plant Alkaloid (Emetine) Inhibits HIV-1 Replication by Interfering with Reverse Transcriptase Activity.

Ipecac alkaloids are secondary metabolites produced in the medicinal plant Psychotria ipecacuanha. Emetine is the main alkaloid of ipecac and one of the active compounds in syrup of Ipecac with emetic property. Here we evaluated emetine's potential as an antiviral agent against Human Immunodeficiency Virus. We performed in vitro Reverse Transcriptase (RT) Assay and Natural Endogenous Reverse Transcriptase Activity Assay (NERT) to evaluate HIV RT inhibition. Emetine molecular docking on HIV-1 RT was also analyzed. Phenotypic assays were performed in non-lymphocytic and in Peripheral Blood Mononuclear Cells (PBMC) with HIV-1 wild-type and HIV-harboring RT-resistant mutation to Nucleoside Reverse Transcriptase Inhibitors (M184V). Our results showed that HIV-1 RT was blocked in the presence of emetine in both models: in vitro reactions with isolated HIV-1 RT and intravirion, measured by NERT. Emetine revealed a strong potential of inhibiting HIV-1 replication in both cellular models, reaching 80% of reduction in HIV-1 infection, with low cytotoxic effect. Emetine also blocked HIV-1 infection of RT M184V mutant. These results suggest that emetine is able to penetrate in intact HIV particles, and bind and block reverse transcription reaction, suggesting that it can be used as anti-HIV microbicide. Taken together, our findings provide additional pharmacological information on the potential therapeutic effects of emetine.

Conformational characterization of ipomotaosides and their recognition by COX-1 and 2.

The aerial parts of Ipomoea batatas are described herein to produce four new resin glycosides, designated as ipomotaosides A, B, C, and D. Ipomotaoside A was found to present inhibitory activity on both cyclooxygenases. However, the conformational elucidation of these molecules may be difficult due to their high flexibility. In this context, the current work presents a conformational characterization of ipomotaosides A-D in aqueous and nonaqueous solvents. The employed protocol includes metadynamics evaluation and unrestrained molecular dynamics simulations (MD). The obtained data provided structural models for the ipomotaosides in good agreement with previous ROESY distances measured in pyridine. Accordingly, the most abundant conformation of ipomotaoside A in solution was employed in flexible docking studies, providing a structural basis for the compound's inhibition of COX enzymes. The so-obtained complex supports resin glycosides' role as original scaffolds for future studies, aiming at structural optimization and development of potential new anti-inflammatory agents.

Unveiling Cas8 Dynamics and Regulation within a transposon-encoded Cascade-TniQ Complex.

Cascade is a class 1, type 1 CRISPR-Cas system with a variety of roles in prokaryote defense, specifically against DNA-based viruses. The Vibrio Cholerae transposon, Tn6677, encodes a variant of the type 1F Cascade known as type 1F-3. This Cascade variant complexes with a homodimer of the transposition protein TniQ and leverages the sequence specificity of Cascade to direct the integration activity of the heteromeric transposase tnsA/B, resulting in site-specific transposition of Tn6677. We desire to uncover the molecular details behind R Loop formation of 'Cascade-TniQ.' Due to the lack of a complete model of Cascade-TniQ available at atom-level resolution, we first build a complete model using AlphaFold V2.1. We then simulate this model via classical molecular dynamics and umbrella sampling to study an important regulatory component within Cascade-TniQ, known as the Cas8 'bundle.' Particularly, we show that this alpha helical bundle experiences a free energy barrier to its large-scale translatory motions and relative free energies of its states primarily dependent on a loop within a Cas7 subunit in Cascade-TniQ. Further, we comment on additional structural and dynamical regulatory points of Cascade-TniQ during R Loop formation, such as Cascade-TniQ backbone rigidity, and the potential role TniQ plays in regulating bundle dynamics. In summary, our outcomes provide the first all-atom dynamic representation of one of the largest CRISPR systems, with information that can contribute to understanding the mechanism of nucleic acid binding and, eventually, to transposase recruitment itself. Such information may prove informative to advance genome engineering efforts.

An alpha-helical lid guides the target DNA toward catalysis in CRISPR-Cas12a.

CRISPR-Cas12a is a powerful RNA-guided genome-editing system that generates double-strand DNA breaks using its single RuvC nuclease domain by a sequential mechanism in which initial cleavage of the non-target strand is followed by target strand cleavage. How the spatially distant DNA target strand traverses toward the RuvC catalytic core is presently not understood. Here, continuous tens of microsecond-long molecular dynamics and free-energy simulations reveal that an α-helical lid, located within the RuvC domain, plays a pivotal role in the traversal of the DNA target strand by anchoring the crRNA:target strand duplex and guiding the target strand toward the RuvC core, as also corroborated by DNA cleavage experiments. In this mechanism, the REC2 domain pushes the crRNA:target strand duplex toward the core of the enzyme, while the Nuc domain aids the bending and accommodation of the target strand within the RuvC core by bending inward. Understanding of this critical process underlying Cas12a activity will enrich fundamental knowledge and facilitate further engineering strategies for genome editing.

Machines on Genes through the Computational Microscope.

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

RNA-mediated Allosteric Activation 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. The protein uses a CRISPR RNA (crRNA) to target RNA sequences, which are cleaved by a composite active site formed by two 'Higher Eukaryotes and Prokaryotes Nucleotide' (HEPN) catalytic domains. In this system, an intriguing form of allosteric communication controls RNA cleavage activity, yet its molecular details are unknown. Here, multiple-microsecond molecular dynamics simulations are combined with graph theory and RNA cleavage assays to decipher this activation mechanism. We show that the binding of a target RNA acts as an allosteric effector of the spatially distant HEPN catalytic cleft, by amplifying the allosteric signals over the dynamical noise, that passes through the buried HEPN interface. Critical residues in this region - N378, R973, and R377 - rearrange their interactions upon target RNA binding, and alter allosteric signalling. Alanine mutation of these residues is experimentally shown to select target RNA over an extended complementary sequence beyond guide-target duplex, for RNA cleavage. Altogether, our findings offer a fundamental understanding of the Cas13a mechanism of action and pave new avenues for the development of more selective RNA-based cleavage and detection tools.

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

The 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 new Cas13a variants for easier-to-implement Cas13-based diagnostics. KEY POINTS: Certain positions in the Cas13a crRNA-target-RNA duplex are particularly sensitive to mismatches.Understanding Cas13a's allosteric activation pathway allowed us to develop novel high-fidelity Cas13a variants.These Cas13a variants and crRNA design strategies achieve superior discrimination of SARS-CoV-2 strains.

Dynamics and mechanisms of CRISPR-Cas9 through the lens of computational methods.

The clustered regularly interspaced short palindromic repeat (CRISPR) genome-editing revolution established the beginning of a new era in life sciences. Here, we review the role of state-of-the-art computations in the CRISPR-Cas9 revolution, from the early refinement of cryo-EM data to enhanced simulations of large-scale conformational transitions. Molecular simulations reported a mechanism for RNA binding and the formation of a catalytically competent Cas9 enzyme, in agreement with subsequent structural studies. Inspired by single-molecule experiments, molecular dynamics offered a rationale for the onset of off-target effects, while graph theory unveiled the allosteric regulation. Finally, the use of a mixed quantum-classical approach established the catalytic mechanism of DNA cleavage. Overall, molecular simulations have been instrumental in understanding the dynamics and mechanism of CRISPR-Cas9, contributing to understanding function, catalysis, allostery, and specificity.