Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of the influenza A virus

Viral proteins frequently mutate to evade host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present a deep analysis of the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveals how the mutations rewired interresidue communications, which underlie long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects, which are further tuned by epistasis. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, relative to 1918 NS1, enhancing its affinity to p85β. Complementary molecular dynamics simulations and graph theory-based network analysis for conformational dynamics on the submicrosecond timescales uncover how these mutations rewire the dynamic network, which underlies the allosteric epistasis. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provide insight into their role in the molecular evolution of NS1.

Dynamic sampling of a surveillance state enables DNA proofreading by Cas9

CRISPR-Cas9 has revolutionized genome engineering applications by programming its single-guide RNA, where high specificity is required. However, the precise molecular mechanism underscoring discrimination between on/off-target DNA sequences, relative to the guide RNA template, remains elusive. Here, using methyl-based NMR to study multiple holoenzymes assembled in vitro, we elucidate a discrete protein conformational state which enables recognition of DNA mismatches at the protospacer adjacent motif (PAM)-distal end. Our results delineate an allosteric pathway connecting a dynamic conformational switch at the REC3 domain, with the sampling of a catalytically competent state by the HNH domain. Our NMR data show that HiFi Cas9 (R691A) increases the fidelity of DNA recognition by stabilizing this "surveillance state" for mismatched substrates, shifting the Cas9 conformational equilibrium away from the active state. These results establish a paradigm of substrate recognition through an allosteric protein-based switch, providing unique insights into the molecular mechanism which governs Cas9 selectivity.

Flexibility in PAM recognition expands DNA targeting in xCas9

xCas9 is an evolved variant of the CRISPR-Cas9 genome editing system, engineered to improve specificity and reduce undesired off-target effects. How xCas9 expands the DNA targeting capability of Cas9 by recognising a series of alternative protospacer adjacent motif (PAM) sequences while ignoring others is unknown. Here, we elucidate the molecular mechanism underlying xCas9's expanded PAM recognition and provide critical insights for expanding DNA targeting. We demonstrate that while wild-type Cas9 enforces stringent guanine selection through the rigidity of its interacting arginine dyad, xCas9 introduces flexibility in R1335, enabling selective recognition of specific PAM sequences. This increased flexibility confers a pronounced entropic preference, which also improves recognition of the canonical TGG PAM. Furthermore, xCas9 enhances DNA binding to alternative PAM sequences during the early evolution cycles, while favouring binding to the canonical PAM in the final evolution cycle. This dual functionality highlights how xCas9 broadens PAM recognition and underscores the importance of fine-tuning the flexibility of the PAM-interacting cleft as a key strategy for expanding the DNA targeting potential of CRISPR-Cas systems. These findings deepen our understanding of DNA recognition in xCas9 and may apply to other CRISPR-Cas systems with similar PAM recognition requirements.

Flexibility in PAM Recognition Expands DNA Targeting in xCas9

xCas9 is an evolved variant of the CRISPR-Cas9 genome editing system, engineered to improve specificity and reduce undesired off-target effects. How xCas9 expands the DNA targeting capability of Cas9 by recognizing a series of alternative Protospacer Adjacent Motif (PAM) sequences while ignoring others is unknown. Here, we elucidate the molecular mechanism underlying xCas9's expanded PAM recognition and provide critical insights for expanding DNA targeting. We demonstrate that while wild-type Cas9 enforces stringent guanine selection through the rigidity of its interacting arginine dyad, xCas9 introduces flexibility in R1335, enabling selective recognition of specific PAM sequences. This increased flexibility confers a pronounced entropic preference, which also improves recognition of the canonical TGG PAM. Furthermore, xCas9 enhances DNA binding to alternative PAM sequences during the early evolution cycles, while favouring binding to the canonical PAM in the final evolution cycle. This dual functionality highlights how xCas9 broadens PAM recognition and underscores the importance of fine-tuning the flexibility of the PAM-interacting cleft as a key strategy for expanding the DNA targeting potential of CRISPR-Cas systems. These findings deepen our understanding of DNA recognition in xCas9 and may apply to other CRISPR-Cas systems with similar PAM recognition requirements.

Exploring CRISPR-Cas9 HNH-Domain-Catalyzed DNA Cleavage Using Accelerated Quantum Mechanical Molecular Mechanical Free Energy Simulation

The target DNA (tDNA) cleavage catalyzed by the CRISPR Cas9 enzyme is a critical step in the Cas9-based genome editing technologies. Previously, the tDNA cleavage from an active SpyCas9 enzyme conformation was modeled by Palermo and co-workers (Nierzwicki et al., Nat. Catal. 2022 5, 912) using ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulations, where the free energy barrier was found to be more favorable than that from a pseudoactive enzyme conformation. In this work, we performed ai-QM/MM simulations based on another catalytically active conformation (PDB 7Z4J) of the Cas9 HNH domain from cryo-electron microscopy experiments. For the wildtype enzyme, we acquired a free energy profile for the tDNA cleavage that is largely consistent with the previous report. Furthermore, we explored the role of the active-site K866 residue on the catalytic efficiency by modeling the K866A mutant and found that the K866A mutation increased the reaction free energy barrier, which is consistent with the experimentally observed reduction in the enzyme activity.

Graph theory approaches for molecular dynamics simulations

Graph theory, a branch of mathematics that focuses on the study of graphs (networks of nodes and edges), provides a robust framework for analysing the structural and functional properties of biomolecules. By leveraging molecular dynamics (MD) simulations, atoms or groups of atoms can be represented as nodes, while their dynamic interactions are depicted as edges. This network-based approach facilitates the characterization of properties such as connectivity, centrality, and modularity, which are essential for understanding the behaviour of molecular systems. This review details the application and development of graph theory-based models in studying biomolecular systems. We introduce key concepts in graph theory and demonstrate their practical applications, illustrating how innovative graph theory approaches can be employed to design biomolecular systems with enhanced functionality. Specifically, we explore the integration of graph theoretical methods with MD simulations to gain deeper insights into complex biological phenomena, such as allosteric regulation, conformational dynamics, and catalytic functions. Ultimately, graph theory has proven to be a powerful tool in the field of molecular dynamics, offering valuable insights into the structural properties, dynamics, and interactions of molecular systems. This review establishes a foundation for using graph theory in molecular design and engineering, highlighting its potential to transform the field and drive advancements in the understanding and manipulation of biomolecular systems.

CRISPR innovations in tissue engineering and gene editing

The CRISPR/Cas9 system is a powerful tool for genome editing, utilizing the Cas9 nuclease and programmable single guide RNA (sgRNA). However, the Cas9 nuclease activity can be disabled by mutation, resulting in catalytically deactivated Cas9 (dCas9). By combining the customizable sgRNA with dCas9, researchers can inhibit specific gene expression (CRISPR interference, CRISPRi) or activate the expression of a target gene (CRISPR activation, CRISPRa). In this review, we present the principles and recent advancements of these CRISPR technologies, as well as their delivery vectors. We also explore their applications in stem cell engineering and regenerative medicine, with a focus on in vitro stem cell fate manipulation and in vivo treatments. These include the prevention of retinal and muscular degeneration, neural regeneration, bone regeneration, cartilage tissue engineering, and the treatment of blood, skin, and liver diseases. Furthermore, we discuss the challenges of translating CRISPR technologies into regenerative medicine and provide future perspectives. Overall, this review highlights the potential of CRISPR in advancing regenerative medicine and offers insights into its application in various areas of research and therapy.

Probing Electrostatic Interactions in DNA-Bound CRISPR/Cas9 Complexes by Molecular Dynamics Simulations

Engineered protein mutations may be exploited to tune molecular interactions in the cellular environment. Here, we have explored the structural consequences of different Cas9 mutations in genome-editing CRISPR/Cas9 systems by means of Molecular Dynamics simulations. We have characterized mutation-induced structural changes and their implications for changes in protein-DNA, DNA-RNA, and DNA-DNA interactions. We present the analysis of multiple trajectories over the cumulative time scale of 7.7 μs, focusing on triple mutations that have been associated with enhancement of genome editing specificity, as well as control mutations. We find that the structural changes induced by the protein mutations are consistent with decreasing the strength of the interaction between Cas9 and the nontarget DNA strand. We discuss the implications of this finding for genome editing specificity.

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

The Vibrio cholerae Cascade-TniQ complex unveiled a new paradigm in biology, demonstrating that CRISPR-associated proteins can direct DNA transposition. Despite the tremendous potential of "knocking in" genes at desired sites, the mechanisms underlying DNA binding and transposition remain elusive. In this system, a conformational change of the Cas8 protein is essential for DNA binding, yet how it occurs is unclear. Here, structural modeling and free energy simulations reconstruct the Cas8 helical bundle and reveal an open-to-close conformational transition at key steps of the complex's function. We show that when Cascade-TniQ binds RNA, the Cas8 bundle changes conformation mediated by the interaction with the Cas7.1 protein. This interaction alleviates unfavorable contacts and synchronizes Cas8's shift with neighboring subunits, lowering the barrier for the transition to the open state, a critical requirement for DNA binding. As DNA fully pairs with RNA, the open state becomes increasingly accessible, favoring interactions with DNA and aiding the formation of an R-loop. These outcomes provide the first dynamic representation of a critical conformational change in one of the largest CRISPR systems and illustrate its role at critical steps of the Cascade-TniQ biophysical function, advancing our understanding of nucleic acid binding and transposition mechanisms.

Deciphering Internal Regulatory Patterns within the p53 Core Tetramer: Insights from Community Network Analysis

Gene therapy is one of the most effective strategies for cancer treatment. The p53 protein, commonly known as the "guardian of the genome", plays a critical role in gene activation and tumor suppression. Tetramerization of the p53 core domain is an essential allosteric process that supports its suppression functions. This letter presents a framework to analyze the structure, function, and dynamic connectivity of the p53 tetramer, using community network analysis based on all-atom molecular dynamics simulations. The communities within the p53 monomer exhibit distinct functional roles, while interactions between molecules establish a symmetrical network structure. We identified direct evidence of single, double, and multiple pathway regulations within the p53 tetramer and crucial residue pairs involved in these connections. Our study provides a comprehensive framework to understand the community network of the p53 tetramer, offering new insights into the stable formation of the p53 core tetramer.

A parameterized two-domain thermodynamic model explains diverse mutational effects on protein allostery

New experimental findings continue to challenge our understanding of protein allostery. Recent deep mutational scanning study showed that allosteric hotspots in the tetracycline repressor (TetR) and its homologous transcriptional factors are broadly distributed rather than spanning well-defined structural pathways as often assumed. Moreover, hotspot mutation-induced allostery loss was rescued by distributed additional mutations in a degenerate fashion. Here, we develop a two-domain thermodynamic model for TetR, which readily rationalizes these intriguing observations. The model accurately captures the in vivo activities of various mutants with changes in physically transparent parameters, allowing the data-based quantification of mutational effects using statistical inference. Our analysis reveals the intrinsic connection of intra- and inter-domain properties for allosteric regulation and illustrate epistatic interactions that are consistent with structural features of the protein. The insights gained from this study into the nature of two-domain allostery are expected to have broader implications for other multi-domain allosteric proteins.

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus

Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.

Secondary Conformational Checkpoint in CRISPR-Cas9

A specific checkpoint between target DNA binding and cleavage primarily governs the precision of Cas9 gene editing. Although various CRISPR-Cas9 variants have been developed to improve DNA cleavage accuracy, we still lack a comprehensive understanding of how they work at the molecular level. Herein, we have focused on studying the late-stage conformational transitions of Cas9 and an evolved Cas9 mutant (evoCas9) that start from the precleavage state. Our submilliseconds of dynamic simulations reveal that the presence of base mismatches leads the HNH nuclease domain of Cas9 to alter its principal functional modes of motion, thereby impairing its conformational activation. This observation suggests the existence of a secondary conformational checkpoint that fine-tunes the final DNA cleavage activation. Remarkably, evoCas9 is prone to deviating from the normal activation pathway with base mismatches. This is characterized by a noticeable shift in the positioning of the HNH domain and a significantly perturbed allosteric communication network within the enzyme. Therefore, the mutations evolved in evoCas9 also reinforce the secondary checkpoint in addition to the previously identified primary checkpoint, collectively ensuring this variant's high gene-editing accuracy. This mechanism should also apply to other Cas9-guide RNA variants with enhanced fidelity.

Engineering CRISPR/Cas9 therapeutics for cancer precision medicine

The discovery of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) technology has revolutionized field of cancer treatment. This review explores usage of CRISPR/Cas9 for editing and investigating genes involved in human carcinogenesis. It provides insights into the development of CRISPR as a genetic tool. Also, it explores recent developments and tools available in designing CRISPR/Cas9 systems for targeting oncogenic genes for cancer treatment. Further, we delve into an overview of cancer biology, highlighting key genetic alterations and signaling pathways whose deletion prevents malignancies. This fundamental knowledge enables a deeper understanding of how CRISPR/Cas9 can be tailored to address specific genetic aberrations and offer personalized therapeutic approaches. In this review, we showcase studies and preclinical trials that show the utility of CRISPR/Cas9 in disrupting oncogenic targets, modulating tumor microenvironment and increasing the efficiency of available anti treatments. It also provides insight into the use of CRISPR high throughput screens for cancer biomarker identifications and CRISPR based screening for drug discovery. In conclusion, this review offers an overview of exciting developments in engineering CRISPR/Cas9 therapeutics for cancer treatment and highlights the transformative potential of CRISPR for innovation and effective cancer treatments.

High-fidelity, hyper-accurate, and evolved mutants rewire atomic-level communication in CRISPR-Cas9

The high-fidelity (HF1), hyper-accurate (Hypa), and evolved (Evo) variants of the CRISPR-associated protein 9 (Cas9) endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in recognition subdomain 3 (Rec3) mediate specificity in these variants are poorly understood. Here, solution nuclear magnetic resonance and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic His-Asn-His (HNH) domain. The variants remodel the communication sourcing from the Rec3 α helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

A parametrized two-domain thermodynamic model explains diverse mutational effects on protein allostery

New experimental findings continue to challenge our understanding of protein allostery. Recent deep mutational scanning study showed that allosteric hotspots in the tetracycline repressor (TetR) and its homologous transcriptional factors are broadly distributed rather than spanning well-defined structural pathways as often assumed. Moreover, hotspot mutation-induced allostery loss was rescued by distributed additional mutations in a degenerate fashion. Here, we develop a two-domain thermodynamic model for TetR, which readily rationalizes these intriguing observations. The model accurately captures the in vivo activities of various mutants with changes in physically transparent parameters, allowing the data-based quantification of mutational effects using statistical inference. Our analysis reveals the intrinsic connection of intra- and inter-domain properties for allosteric regulation and illustrate epistatic interactions that are consistent with structural features of the protein. The insights gained from this study into the nature of two-domain allostery are expected to have broader implications for other multidomain allosteric proteins.

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.

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.

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.

High-Fidelity, Hyper-Accurate, and Evolved Mutants Rewire Atomic Level Communication in CRISPR-Cas9

The Cas9-HF1, HypaCas9, and evoCas9 variants of the Cas9 endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in the Rec3 domain mediate specificity in these variants are poorly understood. Here, solution NMR and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and peculiar dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic HNH domain. The variants remodel the communication sourcing from the Rec3 α-helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Examining the origins of observed terahertz modes from an optically pumped atomistic model protein in aqueous solution

The microscopic origins of terahertz (THz) vibrational modes in biological systems are an active and open area of current research. Recent experiments [Phys Rev X. 8, 031061 (2018)] have revealed the presence of a pronounced mode at ∼0.3 THz in fluorophore-decorated bovine serum albumin (BSA) protein in aqueous solution under nonequilibrium conditions induced by optical pumping. This result was heuristically interpreted as a collective elastic fluctuation originating from the activation of a low-frequency phonon mode. In this work, we show that the sub-THz spectroscopic response emerges in a statistically significant manner (> 2 σ ) from such collective behavior, illustrating how photoexcitation can alter specific THz vibrational modes. We revisit the theoretical analysis with proof-of-concept molecular dynamics that introduce optical excitations into the simulations. Using information theory techniques, we show that these excitations can give rise to a multiscale response involving two optically excited chromophores (tryptophans), other amino acids in the protein, ions, and water. Our results motivate new experiments and fully nonequilibrium simulations to probe these phenomena, as well as the refinement of atomistic models of Fröhlich condensates that are fundamentally determined by nonlinear interactions in biology.

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.

GRAPHICAL ABSTRACT

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.

MDiGest: A Python package for describing allostery from molecular dynamics simulations

Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.

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.

Exploring and Learning the Universe of Protein Allostery Using Artificial Intelligence Augmented Biophysical and Computational Approaches

Allosteric mechanisms are commonly employed regulatory tools used by proteins to orchestrate complex biochemical processes and control communications in cells. The quantitative understanding and characterization of allosteric molecular events are among major challenges in modern biology and require integration of innovative computational experimental approaches to obtain atomistic-level knowledge of the allosteric states, interactions, and dynamic conformational landscapes. The growing body of computational and experimental studies empowered by emerging artificial intelligence (AI) technologies has opened up new paradigms for exploring and learning the universe of protein allostery from first principles. In this review we analyze recent developments in high-throughput deep mutational scanning of allosteric protein functions; applications and latest adaptations of Alpha-fold structural prediction methods for studies of protein dynamics and allostery; new frontiers in integrating machine learning and enhanced sampling techniques for characterization of allostery; and recent advances in structural biology approaches for studies of allosteric systems. We also highlight recent computational and experimental studies of the SARS-CoV-2 spike (S) proteins revealing an important and often hidden role of allosteric regulation driving functional conformational changes, binding interactions with the host receptor, and mutational escape mechanisms of S proteins which are critical for viral infection. We conclude with a summary and outlook of future directions suggesting that AI-augmented biophysical and computer simulation approaches are beginning to transform studies of protein allostery toward systematic characterization of allosteric landscapes, hidden allosteric states, and mechanisms which may bring about a new revolution in molecular biology and drug discovery.

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.

Insights into the Mechanism of CRISPR/Cas9-Based Genome Editing from Molecular Dynamics Simulations

The CRISPR/Cas9 system is a popular genome-editing tool with immense therapeutic potential. It is a simple two-component system (Cas9 protein and RNA) that recognizes the DNA sequence on the basis of RNA:DNA complementarity, and the Cas9 protein catalyzes the double-stranded break in the DNA. In the past decade, near-atomic resolution structures at various stages of the CRISPR/Cas9 DNA editing pathway have been reported along with numerous experimental and computational studies. Such studies have boosted knowledge of the genome-editing mechanism. Despite such advancements, the application of CRISPR/Cas9 in therapeutics is still limited, primarily due to off-target effects. Several studies aim at engineering high-fidelity Cas9 to minimize the off-target effects. Molecular Dynamics (MD) simulations have been an excellent complement to the experimental studies for investigating the mechanism of CRISPR/Cas9 editing in terms of structure, thermodynamics, and kinetics. MD-based studies have uncovered several important molecular aspects of Cas9, such as nucleotide binding, catalytic mechanism, and off-target effects. In this Review, the contribution of MD simulation to understand the CRISPR/Cas9 mechanism has been discussed, preceded by an overview of the history, mechanism, and structural aspects of the CRISPR/Cas9 system. These studies are important for the rational design of highly specific Cas9 and will also be extremely promising for achieving more accurate genome editing in the future.

Twisting and swiveling domain motions in Cas9 to recognize target DNA duplexes, make double-strand breaks, and release cleaved duplexes

The CRISPR-associated protein 9 (Cas9) has been engineered as a precise gene editing tool to make double-strand breaks. CRISPR-associated protein 9 binds the folded guide RNA (gRNA) that serves as a binding scaffold to guide it to the target DNA duplex via a RecA-like strand-displacement mechanism but without ATP binding or hydrolysis. The target search begins with the protospacer adjacent motif or PAM-interacting domain, recognizing it at the major groove of the duplex and melting its downstream duplex where an RNA-DNA heteroduplex is formed at nanomolar affinity. The rate-limiting step is the formation of an R-loop structure where the HNH domain inserts between the target heteroduplex and the displaced non-target DNA strand. Once the R-loop structure is formed, the non-target strand is rapidly cleaved by RuvC and ejected from the active site. This event is immediately followed by cleavage of the target DNA strand by the HNH domain and product release. Within CRISPR-associated protein 9, the HNH domain is inserted into the RuvC domain near the RuvC active site via two linker loops that provide allosteric communication between the two active sites. Due to the high flexibility of these loops and active sites, biophysical techniques have been instrumental in characterizing the dynamics and mechanism of the CRISPR-associated protein 9 nucleases, aiding structural studies in the visualization of the complete active sites and relevant linker structures. Here, we review biochemical, structural, and biophysical studies on the underlying mechanism with emphasis on how CRISPR-associated protein 9 selects the target DNA duplex and rejects non-target sequences.

Probing conformational landscapes and mechanisms of allosteric communication in the functional states of the ABL kinase domain using multiscale simulations and network-based mutational profiling of allosteric residue potentials

In the current study, multiscale simulation approaches and dynamic network methods are employed to examine the dynamic and energetic details of conformational landscapes and allosteric interactions in the ABL kinase domain that determine the kinase functions. Using a plethora of synergistic computational approaches, we elucidate how conformational transitions between the active and inactive ABL states can employ allosteric regulatory switches to modulate intramolecular communication networks between the ATP site, the substrate binding region, and the allosteric binding pocket. A perturbation-based network approach that implements mutational profiling of allosteric residue propensities and communications in the ABL states is proposed. Consistent with biophysical experiments, the results reveal functionally significant shifts of the allosteric interaction networks in which preferential communication paths between the ATP binding site and substrate regions in the active ABL state become suppressed in the closed inactive ABL form, which in turn features favorable allosteric coupling between the ATP site and the allosteric binding pocket. By integrating the results of atomistic simulations with dimensionality reduction methods and Markov state models, we analyze the mechanistic role of macrostates and characterize kinetic transitions between the ABL conformational states. Using network-based mutational scanning of allosteric residue propensities, this study provides a comprehensive computational analysis of long-range communications in the ABL kinase domain and identifies conserved regulatory hotspots that modulate kinase activity and allosteric crosstalk between the allosteric pocket, ATP binding site, and substrate binding regions.

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.

Structural biology of CRISPR-Cas immunity and genome editing enzymes

CRISPR-Cas systems provide resistance against foreign mobile genetic elements and have a wide range of genome editing and biotechnological applications. In this Review, we examine recent advances in understanding the molecular structures and mechanisms of enzymes comprising bacterial RNA-guided CRISPR-Cas immune systems and deployed for wide-ranging genome editing applications. We explore the adaptive and interference aspects of CRISPR-Cas function as well as open questions about the molecular mechanisms responsible for genome targeting. These structural insights reflect close evolutionary links between CRISPR-Cas systems and mobile genetic elements, including the origins and evolution of CRISPR-Cas systems from DNA transposons, retrotransposons and toxin-antitoxin modules. We discuss how the evolution and structural diversity of CRISPR-Cas systems explain their functional complexity and utility as genome editing tools.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

Deep mutational scanning and machine learning reveal structural and molecular rules governing allosteric hotspots in homologous proteins

A fundamental question in protein science is where allosteric hotspots - residues critical for allosteric signaling - are located, and what properties differentiate them. We carried out deep mutational scanning (DMS) of four homologous bacterial allosteric transcription factors (aTFs) to identify hotspots and built a machine learning model with this data to glean the structural and molecular properties of allosteric hotspots. We found hotspots to be distributed protein-wide rather than being restricted to 'pathways' linking allosteric and active sites as is commonly assumed. Despite structural homology, the location of hotspots was not superimposable across the aTFs. However, common signatures emerged when comparing hotspots coincident with long-range interactions, suggesting that the allosteric mechanism is conserved among the homologs despite differences in molecular details. Machine learning with our large DMS datasets revealed global structural and dynamic properties to be a strong predictor of whether a residue is a hotspot than local and physicochemical properties. Furthermore, a model trained on one protein can predict hotspots in a homolog. In summary, the overall allosteric mechanism is embedded in the structural fold of the aTF family, but the finer, molecular details are sequence-specific.

Principles of target DNA cleavage and the role of Mg2+ in the catalysis of CRISPR-Cas9

At the core of the CRISPR-Cas9 genome-editing technology, the endonuclease Cas9 introduces site-specific breaks in DNA. However, precise mechanistic information to ameliorating Cas9 function is still missing. Here, multi-microsecond molecular dynamics, free-energy and multiscale simulations are combined with solution NMR and DNA cleavage experiments to resolve the catalytic mechanism of target DNA cleavage. We show that the conformation of an active HNH nuclease is tightly dependent on the catalytic Mg2+, unveiling its cardinal structural role. This activated Mg2+-bound HNH is consistently described through molecular simulations, solution NMR and DNA cleavage assays, revealing also that the protonation state of the catalytic H840 is strongly affected by active site mutations. Finally, ab-initio QM(DFT)/MM simulations and metadynamics establish the catalytic mechanism, showing that the catalysis is activated by H840 and completed by K866, rationalising DNA cleavage experiments. This information is critical to enhance the enzymatic function of CRISPR-Cas9 toward improved genome-editing.

Deciphering the QR Code of the CRISPR-Cas9 System: Synergy between Gln768 (Q) and Arg976 (R)

Markov state models (MSMs) and machine learning (ML) algorithms can extrapolate the long-time-scale behavior of large biomolecules from molecular dynamics (MD) trajectories. In this study, an MD-MSM-ML scheme has been applied to probe the large endonuclease (Cas9) in the bacterial adaptive immunity CRISPR-Cas9 system. CRISPR has become a programmable and state-of-the-art powerful genome editing tool that has already revolutionized life sciences. CRISPR-Cas9 is programmed to process specific DNA sequences in the genome. However, human/biomedical applications are compromised by off-target DNA damage. Characterization of Cas9 at the structural and biophysical levels is a prerequisite for the development of efficient and high-fidelity Cas9 variants. The Cas9 wild type and two variants (R63A-R66A-R70A, R69A-R71A-R74A-R78A) are studied herein. The configurational space of Cas9 is provided with a focus on the conformations of the side chains of two residues (Gln768 and Arg976). A model for the synergy between those two residues is proposed. The results are discussed within the context of experimental literature. The results and methodology can be exploited for the study of large biomolecules in general and for the engineering of more efficient and safer Cas9 variants for applications.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

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.

Molecular Dynamics Simulations Establish the Molecular Basis for the Broad Allostery Hotspot Distributions in the Tetracycline Repressor

It is imperative to identify the network of residues essential to the allosteric coupling for the purpose of rationally engineering allostery in proteins. Deep mutational scanning analysis has emerged as a function-centric approach for identifying such allostery hotspots in a comprehensive and unbiased fashion, leading to observations that challenge our understanding of allostery at the molecular level. Specifically, a recent deep mutational scanning study of the tetracycline repressor (TetR) revealed an unexpectedly broad distribution of allostery hotspots throughout the protein structure. Using extensive molecular dynamics simulations (up to 50 μs) and free energy computations, we establish the molecular and energetic basis for the strong anticooperativity between the ligand and DNA binding sites. The computed free energy landscapes in different ligation states illustrate that allostery in TetR is well described by a conformational selection model, in which the apo state samples a broad set of conformations, and specific ones are selectively stabilized by either ligand or DNA binding. By examining a range of structural and dynamic properties of residues at both local and global scales, we observe that various analyses capture different subsets of experimentally identified hotspots, suggesting that these residues modulate allostery in distinct ways. These results motivate the development of a thermodynamic model that qualitatively explains the broad distribution of hotspot residues and their distinct features in molecular dynamics simulations. The multifaceted strategy that we establish here for hotspot evaluations and our insights into their mechanistic contributions are useful for modulating protein allostery in mechanistic and engineering studies.

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.