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

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.

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.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

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.

The Electronic Structure of Genome Editors from the First Principles

Genome editing based on the CRISPR-Cas9 system has paved new avenues for medicine, pharmaceutics, biotechnology, and beyond. This article reports the role of first-principles (ab-initio) molecular dynamics (MD) in the CRISPR-Cas9 revolution, achieving a profound understanding of the enzymatic function and offering valuable insights for enzyme engineering. We introduce the methodologies and explain the use of ab-initio MD simulations to characterize the two-metal dependent mechanism of DNA cleavage in the RuvC domain of the Cas9 enzyme, and how a second catalytic domain, HNH, cleaves the target DNA with the aid of a single metal ion. A detailed description of how ab-initio MD is combined with free-energy methods - i.e., thermodynamic integration and metadynamics - to break and form chemical bonds is given, explaining the use of these methods to determine the chemical landscape and establish the catalytic mechanism in CRISPR-Cas9. The critical role of classical methods is also discussed, explaining theory and application of constant pH MD simulations, used to accurately predict the catalytic residues' protonation states. Overall, first-principles methods are shown to unravel the electronic structure of the Cas9 enzyme, providing valuable insights that can serve for the design of genome editing tools with improved catalytic efficiency or controllable activity.

Full-Length Model of SaCas9-sgRNA-DNA Complex in Cleavage State

Staphylococcus aureus Cas9 (SaCas9) is a widely used genome editing tool. Understanding its molecular mechanisms of DNA cleavage could effectively guide the engineering optimization of this system. Here, we determined the first cryo-electron microscopy structure of the SaCas9-sgRNA-DNA ternary complex. This structure reveals that the HNH nuclease domain is tightly bound to the cleavage site of the target DNA strand, and is in close contact with the WED and REC domains. Moreover, it captures the complete structure of the sgRNA, including the previously unresolved stem-loop 2. Based on this structure, we build a full-length model for the ternary complex in cleavage state. This model enables identification of the residues for the interactions between the HNH domain and the WED and REC domains. Moreover, we found that the stem-loop 2 of the sgRNA tightly binds to the PI and RuvC domains and may also regulate the position shift of the RuvC domain. Further mutagenesis and molecular dynamics simulations supported the idea that the interactions of the HNH domain with the WED and REC domains play an important role in the DNA cleavage. Thus, this study provides new mechanistic insights into the DNA cleavage of SaCas9 and is also useful for guiding the future engineering of SaCas9-mediated gene editing systems.

Constructing next-generation CRISPR-Cas tools from structural blueprints

Clustered regularly interspaced short palindromic repeats - CRISPR-associated protein (CRISPR-Cas) systems are a critical component of the bacterial adaptive immune response. Since the discovery that they can be reengineered as programmable RNA-guided nucleases, there has been significant interest in using these systems to perform diverse and precise genetic manipulations. Here, we outline recent advances in the mechanistic understanding of CRISPR-Cas9, how these findings have been leveraged in the rational redesign of Cas9 variants with altered activities, and how these novel tools can be exploited for biotechnology and therapeutics. We also discuss the potential of the ubiquitous, yet often-overlooked, multisubunit CRISPR effector complexes for large-scale genomic deletions. Furthermore, we highlight how future structural studies will bolster these technologies.