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

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.

Structures of Neisseria meningitidis Cas9 Complexes in Catalytically Poised and Anti-CRISPR-Inhibited States

High-resolution Cas9 structures have yet to reveal catalytic conformations due to HNH nuclease domain positioning away from the cleavage site. Nme1Cas9 and Nme2Cas9 are compact nucleases for in vivo genome editing. Here, we report structures of meningococcal Cas9 homologs in complex with sgRNA, dsDNA, or the AcrIIC3 anti-CRISPR protein. DNA-bound structures represent an early step of target recognition, a later HNH pre-catalytic state, the HNH catalytic state, and a cleaved-target-DNA-bound state. In the HNH catalytic state of Nme1Cas9, the active site is seen poised at the scissile phosphodiester linkage of the target strand, providing a high-resolution view of the active conformation. The HNH active conformation activates the RuvC domain. Our structures explain how Nme1Cas9 and Nme2Cas9 read distinct PAM sequences and how AcrIIC3 inhibits Nme1Cas9 activity. These structures provide insights into Cas9 domain rearrangements, guide-target engagement, cleavage mechanism, and anti-CRISPR inhibition, facilitating the optimization of these genome-editing platforms.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.

Disclosing the Impact of Carcinogenic SF3b Mutations on Pre-mRNA Recognition Via All-Atom Simulations

The spliceosome accurately promotes precursor messenger-RNA splicing by recognizing specific noncoding intronic tracts including the branch point sequence (BPS) and the 3'-splice-site (3'SS). Mutations of Hsh155 (yeast)/SF3B1 (human), which is a protein of the SF3b factor involved in BPS recognition and induces altered BPS binding and 3'SS selection, lead to mis-spliced mRNA transcripts. Although these mutations recur in hematologic malignancies, the mechanism by which they change gene expression remains unclear. In this study, multi-microsecond-long molecular-dynamics simulations of eighth distinct ∼700,000 atom models of the spliceosome Bact complex, and gene sequencing of SF3B1, disclose that these carcinogenic isoforms destabilize intron binding and/or affect the functional dynamics of Hsh155/SF3B1 only when binding non-consensus BPSs, as opposed to the non-pathogenic variants newly annotated here. This pinpoints a cross-talk between the distal Hsh155 mutation and BPS recognition sites. Our outcomes unprecedentedly contribute to elucidating the principles of pre-mRNA recognition, which provides critical insights on the mechanism underlying constitutive/alternative/aberrant splicing.

Molecular basis for the PAM expansion and fidelity enhancement of an evolved Cas9 nuclease

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems have been harnessed as powerful genome editing tools in diverse organisms. However, the off-target effects and the protospacer adjacent motif (PAM) compatibility restrict the therapeutic applications of these systems. Recently, a Streptococcus pyogenes Cas9 (SpCas9) variant, xCas9, was evolved to possess both broad PAM compatibility and high DNA fidelity. Through determination of multiple xCas9 structures, which are all in complex with single-guide RNA (sgRNA) and double-stranded DNA containing different PAM sequences (TGG, CGG, TGA, and TGC), we decipher the molecular mechanisms of the PAM expansion and fidelity enhancement of xCas9. xCas9 follows a unique two-mode PAM recognition mechanism. For non-NGG PAM recognition, xCas9 triggers a notable structural rearrangement in the DNA recognition domains and a rotation in the key PAM-interacting residue R1335; such mechanism has not been observed in the wild-type (WT) SpCas9. For NGG PAM recognition, xCas9 applies a strategy similar to WT SpCas9. Moreover, biochemical and cell-based genome editing experiments pinpointed the critical roles of the E1219V mutation for PAM expansion and the R324L, S409I, and M694I mutations for fidelity enhancement. The molecular-level characterizations of the xCas9 nuclease provide critical insights into the mechanisms of the PAM expansion and fidelity enhancement of xCas9 and could further facilitate the engineering of SpCas9 and other Cas9 orthologs.

Structural and functional insights into the bona fide catalytic state of Streptococcus pyogenes Cas9 HNH nuclease domain

The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (SpyCas9), along with a programmable single-guide RNA (sgRNA), has been exploited as a significant genome-editing tool. Despite the recent advances in determining the SpyCas9 structures and DNA cleavage mechanism, the cleavage-competent conformation of the catalytic HNH nuclease domain of SpyCas9 remains largely elusive and debatable. By integrating computational and experimental approaches, we unveiled and validated the activated Cas9-sgRNA-DNA ternary complex in which the HNH domain is neatly poised for cleaving the target DNA strand. In this catalysis model, the HNH employs the catalytic triad of D839-H840-N863 for cleavage catalysis, rather than previously implicated D839-H840-D861, D837-D839-H840, or D839-H840-D861-N863. Our study contributes critical information to defining the catalytic conformation of the HNH domain and advances the knowledge about the conformational activation underlying Cas9-mediated DNA cleavage.

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.

Deciphering Off-Target Effects in CRISPR-Cas9 through Accelerated Molecular Dynamics

CRISPR-Cas9 is the state-of-the-art technology for editing and manipulating nucleic acids. However, the occurrence of off-target mutations can limit its applicability. Here, all-atom enhanced molecular dynamics (MD) simulations-using Gaussian accelerated MD (GaMD)-are used to decipher the mechanism of off-target binding at the molecular level. GaMD reveals that base pair mismatches in the target DNA at distal sites with respect to the protospacer adjacent motif (PAM) can induce an extended opening of the RNA:DNA heteroduplex, which leads to newly formed interactions between the unwound DNA and the L2 loop of the catalytic HNH domain. These conserved interactions constitute a "lock" effectively decreasing the conformational freedom of the HNH domain and hampering its activation for cleavage. Remarkably, depending on their positions at PAM distal sites, DNA mismatches responsible for off-target cleavages are unable to "lock" the HNH domain, thereby leading to the unselective cleavage of DNA sequences. In consistency with the available experimental data, the ability to "lock" the catalytic HNH domain in an inactive "conformational checkpoint" is shown to be a key determinant in the onset of off-target effects. This mechanistic rationale contributes in clarifying a long lasting open issue in the CRISPR-Cas9 function and poses the foundation for designing novel and more specific Cas9 variants, which could be obtained by magnifying the "locking" interactions between HNH and the target DNA in the presence of any incorrect off-target sequence, thus preventing undesired cleavages.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

The invisible dance of CRISPR-Cas9. Simulations unveil the molecular side of the gene-editing revolution

Since the discovery of the DNA double helix, the main molecular repository of genetic information, scientists have been struggling to find ways to efficiently manipulate genes. The ability to mark, modify, or regulate specific sequences of DNA in a controlled fashion is of key importance because of the ways that gene editing could be used to improve human life. For example, genetic therapies are being developed to permanently cure cancer and other life-threatening diseases.

Structure and Dynamics of the CRISPR-Cas9 Catalytic Complex

CRISPR-Cas9 is a bacterial immune system with exciting applications for genome editing. In spite of extensive experimental characterization, the active site chemistry of the RuvC domain-which performs DNA cleavages-has remained elusive. Its knowledge is key for structure-based engineering aimed at improving DNA cleavages. Here, we deliver an in-depth characterization by using quantum-classical (QM/MM) molecular dynamics (MD) simulations and a Gaussian accelerated MD method, coupled with bioinformatics analysis. We disclose a two-metal aided architecture in the RuvC active site, which is poised to operate DNA cleavages, in analogy with other DNA/RNA processing enzymes. The conformational dynamics of the RuvC domain further reveals that an "arginine finger" stably contacts the scissile phosphate, with the function of stabilizing the active complex. Remarkably, the formation of a catalytically competent state of the RuvC domain is only observed upon the conformational activation of the other nuclease domain of CRISPR-Cas9-i.e., the HNH domain-such allowing concerted cleavages of double stranded DNA. This structure is in agreement with the available experimental data and remarkably differs from previous models based on classical mechanics, demonstrating also that only quantum mechanical simulations can accurately describe the metal-aided active site in CRISPR-Cas9. This fully catalytic structure-in which both the HNH and RuvC domains are prone to perform DNA cleavages-constitutes a stepping-stone for understanding DNA cleavage and specificity. It calls for novel experimental verifications and offers the structural foundations for engineering efforts aimed at improving the genome editing capability of CRISPR-Cas9.

Exploring the Catalytic Mechanism of Cas9 Using Information Inferred from Endonuclease VII

Elucidating the nature of the gene editing mechanism of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) is an important task in view of the role of this breakthrough to the advancement of human medicine. In particular, it is crucial to understand the catalytic mechanism of Cas9 (one of the CRISPR associated proteins) and its role in confirming accurate editing. Thus, we focus in this work on an attempt to analyze the catalytic mechanism of Cas9. Considering the absence of detailed structural information on the active form of Cas9, we use an empirical valence bond (EVB) which is calibrated on the closely related mechanism of T4 endonuclease VII. The calibrated EVB is then used in studying the reaction of Cas9, while trying several structural models. It is found that the catalytic activation requires a large conformational change, where K848 or other positively charged group moves from a relatively large distance toward the scissile phosphate. This conformational change leads to the change in position of the Mg2+ ion and to a major reduction in the activation barrier for the catalytic reaction. Our finding provides an important clue on the nature of the catalytic activation of CAS9 and thus should help in elucidating a key aspect of the gene editing process. For example, the approach used here should be effective in exploring the nature of off target activation and its relationship to the energetics of the unwinding process. This strategy may offer ways to improve the selectivity of Cas9.

Eigenvector centrality for characterization of protein allosteric pathways

Determining the principal energy-transfer pathways responsible for allosteric communication in biomolecules remains challenging, partially due to the intrinsic complexity of the systems and the lack of effective characterization methods. In this work, we introduce the eigenvector centrality metric based on mutual information to elucidate allosteric mechanisms that regulate enzymatic activity. Moreover, we propose a strategy to characterize the range of correlations that underlie the allosteric processes. We use the V-type allosteric enzyme imidazole glycerol phosphate synthase (IGPS) to test the proposed methodology. The eigenvector centrality method identifies key amino acid residues of IGPS with high susceptibility to effector binding. The findings are validated by solution NMR measurements yielding important biological insights, including direct experimental evidence for interdomain motion, the central role played by helix h[Formula: see text], and the short-range nature of correlations responsible for the allosteric mechanism. Beyond insights on IGPS allosteric pathways and the nature of residues that could be targeted by therapeutic drugs or site-directed mutagenesis, the reported findings demonstrate the eigenvector centrality analysis as a general cost-effective methodology to gain fundamental understanding of allosteric mechanisms at the molecular level.