Striking Plasticity of CRISPR-Cas9 and Key Role of Non-target DNA, as Revealed by Molecular Simulations

Conformational plasticity of SpyCas9 induced by AcrIIA4 and AcrIIA2: Insights from molecular dynamics simulation

CRISPR-Cas9 systems constitute bacterial adaptive immune systems that protect against phage infections. Bacteriophages encode anti-CRISPR proteins (Acrs) that mitigate the bacterial immune response. However, the structural basis for their inhibitory actions from a molecular perspective remains elusive. In this study, through microsecond atomistic molecular dynamics simulations, we demonstrated the remarkable flexibility of Streptococcus pyogenes Cas9 (SpyCas9) and its conformational adaptability during interactions with AcrIIA4 and AcrIIA2. Specifically, we demonstrated that the binding of AcrIIA4 and AcrIIA2 to SpyCas9 induces a conformational rearrangement that causes spatial separation between the nuclease and cleavage sites, thus making the endonuclease inactive. This separation disrupts the transmission of signals between the protospacer adjacent motif recognition and nuclease domains, thereby impeding the efficient processing of double-stranded DNA. The simulation also reveals that AcrIIA4 and AcrIIA2 cause different structural variations of SpyCas9. Our research illuminates the precise mechanisms underlying the suppression of SpyCas9 by AcrIIA4 and AcrIIA2, thus presenting new possibilities for controlling genome editing with higher accuracy.

Engineering salt-tolerant Cas12f1 variants for gene-editing applications

CRISPR has revolutionized the field of genome editing in life sciences by serving as a versatile and state-of-the-art tool. Cas12f1 is a small nuclease of the bacterial immunity CRISPR system with an ideal size for cellular delivery, in contrast to CRISPR-associated (Cas) proteins like Cas9 or Cas12. However, Cas12f1 works best at low salt concentrations. In this study, we find that the plasticity of certain Cas12f1 regions (K196-Y202 and I452-L515) is negatively affected by increased salt concentrations. On this line, key protein domains (REC1, WED, Nuc, lid) that are involved in the DNA-target recognition and the activation of the catalytic RuvC domain are in turn also affected. We suggest that salt concentration should be taken in to consideration for activity assessments of Cas engineered variants, especially if the mutations are on the protospacer adjacent motif interacting domain. The results can be exploited for the engineering of Cas variants and the assessment of their activity at varying salt concentrations. We propose that the K198Q mutation can restore at great degree the compromised plasticity and could potentially lead to salt-tolerant Cas12f1 variants. The methodology can be also employed for the study of biomolecules in terms of their salinity tolerance.Communicated by Ramaswamy H. Sarma.

Massively parallel analysis of single-molecule dynamics on next-generation sequencing chips

Single-molecule techniques are ideally poised to characterize complex dynamics but are typically limited to investigating a small number of different samples. However, a large sequence or chemical space often needs to be explored to derive a comprehensive understanding of complex biological processes. Here we describe multiplexed single-molecule characterization at the library scale (MUSCLE), a method that combines single-molecule fluorescence microscopy with next-generation sequencing to enable highly multiplexed observations of complex dynamics. We comprehensively profiled the sequence dependence of DNA hairpin properties and Cas9-induced target DNA unwinding-rewinding dynamics. The ability to explore a large sequence space for Cas9 allowed us to identify a number of target sequences with unexpected behaviors. We envision that MUSCLE will enable the mechanistic exploration of many fundamental biological processes.

Exploring the Thermostability of CRISPR-Cas12b using Molecular Dynamics Simulations

CRISPR (clustered regularly interspaced short palindromic repeat)-based diagnostics are at the forefront of rapid detection platforms of infectious diseases. The integration of reverse transcription-loop-mediated isothermal amplification (RT-LAMP) with CRISPR-Cas protein systems has led to the creation of advanced one-pot assays. The sensitivity of these assays has been bolstered by the utilization of a thermophilic Cas12 protein, BrCas12b, and its engineered variant, which exhibits enhanced thermal stability and allows for broader operation temperatures of the assay. Here, we perform all-atom molecular dynamics (MD) simulations on wild-type and mutant BrCas12b to reveal the mechanism of stabilization conferred by the mutation. High-temperature simulations reveal a small structural change along with greater flexibility in the PAM-interacting domain of the mutant BrCas12b, with marginal structural and flexibility changes in the other mutated domains. Comparative essential dynamics analysis between the wild-type and mutant BrCas12b at both ambient and elevated temperatures provides insights into the stabilizing effects of the mutations. Our findings not only offer a comprehensive insight into the dynamic alterations induced by mutations but reveal important motions in BrCas12b, important for the rational design of diagnostic and therapeutic platforms of Cas12 proteins.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

AcrIIA28 is a metalloprotein that specifically inhibits targeted-DNA loading to SpyCas9 by binding to the REC3 domain

CRISPR-Cas systems serve as adaptive immune systems in bacteria and archaea, protecting against phages and other mobile genetic elements. However, phages and archaeal viruses have developed countermeasures, employing anti-CRISPR (Acr) proteins to counteract CRISPR-Cas systems. Despite the revolutionary impact of CRISPR-Cas systems on genome editing, concerns persist regarding potential off-target effects. Therefore, understanding the structural and molecular intricacies of diverse Acrs is crucial for elucidating the fundamental mechanisms governing CRISPR-Cas regulation. In this study, we present the structure of AcrIIA28 from Streptococcus phage Javan 128 and analyze its structural and functional features to comprehend the mechanisms involved in its inhibition of Cas9. Our current study reveals that AcrIIA28 is a metalloprotein that contains Zn2+ and abolishes the cleavage activity of Cas9 only from Streptococcus pyrogen (SpyCas9) by directly interacting with the REC3 domain of SpyCas9. Furthermore, we demonstrate that the AcrIIA28 interaction prevents the target DNA from being loaded onto Cas9. These findings indicate the molecular mechanisms underlying AcrIIA28-mediated Cas9 inhibition and provide valuable insights into the ongoing evolutionary battle between bacteria and phages.

Recent Updates of the CRISPR/Cas9 Genome Editing System: Novel Approaches to Regulate Its Spatiotemporal Control by Genetic and Physicochemical Strategies

The genome editing approach by clustered regularly interspaced short palindromic repeats (CRISPR)/associated protein 9 (CRISPR/Cas9) is a revolutionary advancement in genetic engineering. Owing to its simple design and powerful genome-editing capability, it offers a promising strategy for the treatment of different infectious, metabolic, and genetic diseases. The crystal structure of Streptococcus pyogenes Cas9 (SpCas9) in complex with sgRNA and its target DNA at 2.5 Å resolution reveals a groove accommodating sgRNA:DNA heteroduplex within a bilobate architecture with target recognition (REC) and nuclease (NUC) domains. The presence of a PAM is significantly required for target recognition, R-loop formation, and strand scission. Recently, the spatiotemporal control of CRISPR/Cas9 genome editing has been considerably improved by genetic, chemical, and physical regulatory strategies. The use of genetic modifiers anti-CRISPR proteins, cell-specific promoters, and histone acetyl transferases has uplifted the application of CRISPR/Cas9 as a future-generation genome editing tool. In addition, interventions by chemical control, small-molecule activators, oligonucleotide conjugates and bioresponsive delivery carriers have improved its application in other areas of biological fields. Furthermore, the intermediation of physical control by using heat-, light-, magnetism-, and ultrasound-responsive elements attached to this molecular tool has revolutionized genome editing further. These strategies significantly reduce CRISPR/Cas9's undesirable off-target effects. However, other undesirable effects still offer some challenges for comprehensive clinical translation using this genome-editing approach. In this review, we summarize recent advances in CRISPR/Cas9 structure, mechanistic action, and the role of small-molecule activators, inhibitors, promoters, and physical approaches. Finally, off-target measurement approaches, challenges, future prospects, and clinical applications are discussed.

Unity among the diverse RNA-guided CRISPR-Cas interference mechanisms

CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated) systems are adaptive immune systems that protect bacteria and archaea from invading mobile genetic elements (MGEs). The Cas protein-CRISPR RNA (crRNA) complex uses complementarity of the crRNA "guide" region to specifically recognize the invader genome. CRISPR effectors that perform targeted destruction of the foreign genome have emerged independently as multi-subunit protein complexes (Class 1 systems) and as single multi-domain proteins (Class 2). These different CRISPR-Cas systems can cleave RNA, DNA, and protein in an RNA-guided manner to eliminate the invader, and in some cases, they initiate programmed cell death/dormancy. The versatile mechanisms of the different CRISPR-Cas systems to target and destroy nucleic acids have been adapted to develop various programmable-RNA-guided tools and have revolutionized the development of fast, accurate, and accessible genomic applications. In this review, we present the structure and interference mechanisms of different CRISPR-Cas systems and an analysis of their unified features. The three types of Class 1 systems (I, III, and IV) have a conserved right-handed helical filamentous structure that provides a backbone for sequence-specific targeting while using unique proteins with distinct mechanisms to destroy the invader. Similarly, all three Class 2 types (II, V, and VI) have a bilobed architecture that binds the RNA-DNA/RNA hybrid and uses different nuclease domains to cleave invading MGEs. Additionally, we highlight the mechanistic similarities of CRISPR-Cas enzymes with other RNA-cleaving enzymes and briefly present the evolutionary routes of the different CRISPR-Cas systems.

Why Does the E1219V Mutation Expand T-Rich PAM Recognition in Cas9 from Streptococcus pyogenes?

Popular RNA-guided DNA endonuclease Cas9 from Streptococcus pyogenes (SpCas9) recognizes the canonical 5'-NGG-3' protospacer adjacent motif (PAM) and triggers double-stranded DNA cleavage activity. Mutations in SpCas9 were demonstrated to expand the PAM readability and hold promise for therapeutic and genome editing applications. However, the energetics of the PAM recognition and its relation to the atomic structure remain unknown. Using the X-ray structure (precatalytic SpCas9:sgRNA:dsDNA) as a template, we calculated the change in the PAM binding affinity in response to SpCas9 mutations using computer simulations. The E1219V mutation in SpCas9 fine-tunes the water accessibility in the PAM binding pocket and promotes new interactions in the SpCas9:noncanonical T-rich PAM, thus weakening the PAM stringency. The nucleotide-specific interaction of two arginine residues (i.e., R1333 and R1335 of SpCas9) ensured stringent 5'-NGG-3' PAM recognition. R1335A substitution (SpCas9R1335A) completely disrupts the direct interaction between SpCas9 and PAM sequences (canonical or noncanonical), accounting for the loss of editing activity. Interestingly, the double mutant (SpCas9R1335A,E1219V) boosts DNA binding affinity by favoring protein:PAM electrostatic contact in a desolvated pocket. The underlying thermodynamics explain the varied DNA cleavage activity of SpCas9 variants. A direct link between the energetics, structures, and activity is highlighted, which can aid in the rational design of improved SpCas9-based genome editing tools.

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.

Recent Advances in Understanding the Molecular Mechanisms of Multidrug Resistance and Novel Approaches of CRISPR/Cas9-Based Genome-Editing to Combat This Health Emergency

The rapid spread of multidrug resistance (MDR), due to abusive use of antibiotics has led to global health emergency, causing substantial morbidity and mortality. Bacteria attain MDR by different means such as antibiotic modification/degradation, target protection/modification/bypass, and enhanced efflux mechanisms. The classical approaches of counteracting MDR bacteria are expensive and time-consuming, thus, it is highly significant to understand the molecular mechanisms of this resistance to curb the problem from core level. The revolutionary approach of clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated sequence 9 (CRISPR/Cas9), considered as a next-generation genome-editing tool presents an innovative opportunity to precisely target and edit bacterial genome to alter their MDR strategy. Different bacteria possessing antibiotic resistance genes such as mecA, ermB, ramR, tetA, mqrB and blaKPC that have been targeted by CRISPR/Cas9 to re-sensitize these pathogens against antibiotics, such as methicillin, erythromycin, tigecycline, colistin and carbapenem, respectively. The CRISPR/Cas9 from S. pyogenes is the most widely studied genome-editing tool, consisting of a Cas9 DNA endonuclease associated with tracrRNA and crRNA, which can be systematically coupled as sgRNA. The targeting strategies of CRISPR/Cas9 to bacterial cells is mediated through phage, plasmids, vesicles and nanoparticles. However, the targeting approaches of this genome-editing tool to specific bacteria is a challenging task and still remains at a very preliminary stage due to numerous obstacles awaiting to be solved. This review elaborates some recent updates about the molecular mechanisms of antibiotic resistance and the innovative role of CRISPR/Cas9 system in modulating these resistance mechanisms. Furthermore, the delivery approaches of this genome-editing system in bacterial cells are discussed. In addition, some challenges and future prospects are also described.

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.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Differential Divalent Metal Binding by SpyCas9's RuvC Active Site Contributes to Nonspecific DNA Cleavage

To protect against mobile genetic elements (MGEs), some bacteria and archaea have clustered regularly interspaced short palindromic repeats-CRISPR associated (CRISPR-Cas) adaptive immune systems. CRISPR RNAs (crRNAs) bound to Cas nucleases hybridize to MGEs based on sequence complementarity to guide the nucleases to cleave the MGEs. This programmable DNA cleavage has been harnessed for gene editing. Safety concerns include off-target and guide RNA (gRNA)-free DNA cleavages, both of which are observed in the Cas nuclease commonly used for gene editing, Streptococcus pyogenes Cas9 (SpyCas9). We developed a SpyCas9 variant (SpyCas9H982A) devoid of gRNA-free DNA cleavage activity that is more selective for on-target cleavage. The H982A substitution in the metal-dependent RuvC active site reduces Mn2+-dependent gRNA-free DNA cleavage by ∼167-fold. Mechanistic molecular dynamics analysis shows that Mn2+, but not Mg2+, produces a gRNA-free DNA cleavage competent state that is disrupted by the H982A substitution. Our study demonstrates the feasibility of modulating cation:protein interactions to engineer safer gene editing tools.

Genome-wide CRISPR off-target prediction and optimization using RNA-DNA interaction fingerprints

The powerful CRISPR genome editing system is hindered by its off-target effects, and existing computational tools achieved limited performance in genome-wide off-target prediction due to the lack of deep understanding of the CRISPR molecular mechanism. In this study, we propose to incorporate molecular dynamics (MD) simulations in the computational analysis of CRISPR system, and present CRISOT, an integrated tool suite containing four related modules, i.e., CRISOT-FP, CRISOT-Score, CRISOT-Spec, CRISORT-Opti for RNA-DNA molecular interaction fingerprint generation, genome-wide CRISPR off-target prediction, sgRNA specificity evaluation and sgRNA optimization of Cas9 system respectively. Our comprehensive computational and experimental tests reveal that CRISOT outperforms existing tools with extensive in silico validations and proof-of-concept experimental validations. In addition, CRISOT shows potential in accurately predicting off-target effects of the base editors and prime editors, indicating that the derived RNA-DNA molecular interaction fingerprint captures the underlying mechanisms of RNA-DNA interaction among distinct CRISPR systems. Collectively, CRISOT provides an efficient and generalizable framework for genome-wide CRISPR off-target prediction, evaluation and sgRNA optimization for improved targeting specificity in CRISPR genome editing.

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.

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.

Dynamics of Target DNA Binding and Cleavage by Staphylococcus aureus Cas9 as Revealed by High-Speed Atomic Force Microscopy

Programmable DNA binding and cleavage by CRISPR-Cas9 has revolutionized the life sciences. However, the off-target cleavage observed in DNA sequences with some homology to the target still represents a major limitation for a more widespread use of Cas9 in biology and medicine. For this reason, complete understanding of the dynamics of DNA binding, interrogation and cleavage by Cas9 is crucial to improve the efficiency of genome editing. Here, we use high-speed atomic force microscopy (HS-AFM) to investigate Staphylococcus aureus Cas9 (SaCas9) and its dynamics of DNA binding and cleavage. Upon binding to single-guide RNA (sgRNA), SaCas9 forms a close bilobed structure that transiently and flexibly adopts also an open configuration. The SaCas9-mediated DNA cleavage is characterized by release of cleaved DNA and immediate dissociation, confirming that SaCas9 operates as a multiple turnover endonuclease. According to present knowledge, the process of searching for target DNA is mainly governed by three-dimensional diffusion. Independent HS-AFM experiments show a potential long-range attractive interaction between SaCas9-sgRNA and its target DNA. The interaction precedes the formation of the stable ternary complex and is observed exclusively in the vicinity of the protospacer-adjacent motif (PAM), up to distances of several nanometers. The direct visualization of the process by sequential topographic images suggests that SaCas9-sgRNA binds to the target sequence first, while the following binding of the PAM is accompanied by local DNA bending and formation of the stable complex. Collectively, our HS-AFM data reveal a potential and unexpected behavior of SaCas9 during the search for DNA targets.

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.

The Many (Inter)faces of Anti-CRISPRs: Modulation of CRISPR-Cas Structure and Dynamics by Mechanistically Diverse Inhibitors

The discovery of protein inhibitors of CRISPR-Cas systems, called anti-CRISPRs (Acrs), has enabled the development of highly controllable and precise CRISPR-Cas tools. Anti-CRISPRs share very little structural or sequential resemblance to each other or to other proteins, which raises intriguing questions regarding their modes of action. Many structure-function studies have shed light on the mechanism(s) of Acrs, which can act as orthosteric or allosteric inhibitors of CRISPR-Cas machinery, as well as enzymes that irreversibly modify CRISPR-Cas components. Only recently has the breadth of diversity of Acr structures and functions come to light, and this remains a rapidly evolving field. Here, we draw attention to a plethora of Acr mechanisms, with particular focus on how their action toward Cas proteins modulates conformation, dynamic (allosteric) signaling, nucleic acid binding, and cleavage ability.

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.

Controlling CRISPR-Cas9 by guide RNA engineering

The clustered regularly interspaced short palindromic repeats (CRISPR) system is a product of million years of evolution by microbes to fight against invading genetic materials. Around 10 years ago, scientists started to repurpose the CRISPR as genetic tools by molecular engineering approaches. The guide RNA provides a versatile and unique platform for the innovation to improve and expand the application of CRISPR-Cas9 system. In this review, we will first introduce the basic sequence and structure of guide RNA and its role during the function of CRISPR-Cas9. We will then summarize recent progress on the development of various guide RNA engineering strategies. These strategies have been dedicated to improve the performance of CRISPR-Cas9, to achieve precise spatiotemporal control of CRISPR-Cas9, and to broaden the application of CRISPR-Cas9. Finally, we will briefly discuss the uniqueness and advantage of guide RNA-engineering based systems versus those with engineered Cas9 proteins and speculate potential future directions in guide RNA engineering. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Methods > RNA Nanotechnology Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Mechanistic insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase

Calcium (Ca²⁺) signaling plays an important role in the regulation of many cellular functions. Ca²⁺-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca²⁺/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca²⁺ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca²⁺-bound and-unbound states to reveal the importance of Ca²⁺. MD simulations revealed that removal of Ca²⁺ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca²⁺-unbound state. Structural analysis further revealed that Ca²⁺ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca²⁺ in physiological processes.

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.

Frankenstein Cas9: engineering improved gene editing systems

The discovery of CRISPR-Cas9 and its widespread use has revolutionised and propelled research in biological sciences. Although the ability to target Cas9's nuclease activity to specific sites via an easily designed guide RNA (gRNA) has made it an adaptable gene editing system, it has many characteristics that could be improved for use in biotechnology. Cas9 exhibits significant off-target activity and low on-target nuclease activity in certain contexts. Scientists have undertaken ambitious protein engineering campaigns to bypass these limitations, producing several promising variants of Cas9. Cas9 variants with improved and alternative activities provide exciting new tools to expand the scope and fidelity of future CRISPR applications.

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.

Decrypting the mechanistic basis of CRISPR/Cas9 protein

CRISPR/Cas system, a newly but extensively investigated genome-editing method, harbors practical solutions for various genetic problems. It relies on short guide RNAs (gRNAs) to recruit the Cas9 protein, a DNA cleaving enzyme, to its genomic target DNAs. The Cas9 enzyme exhibits some unique properties, like the ability to differentiate self vs. non-self - DNA strands using the base-pairing potential of crRNA, i.e., only CRISPR DNA is entirely complementary to the CRISPR repeat sequences at the crRNA whereas the presence of mismatches in the upstream region of the spacer permit CRISPR interference which is inhibited in case of CRISPR-DNA, allosteric regulation in its domains, and domain reorientation on sgRNA binding. Several groups have contributed their efforts in understanding the functioning of the CRISPR/Cas system, but even then, there is a lot more to explore in this area. The structural and sequence-based understanding of the whole CRISPR-associated bacterial ortholog family landscape is still ambiguous. A better understanding of the underlying energetics of the CRISPR/Cas9 system should reveal critical parameters to design better CRISPR/Cas9s.

Comparative Structural and Dynamics Study of Free and gRNA-bound FnCas9 and SpCas9 Proteins

•

gRNA binding caused Cas9 to open up, exposing hydrophobic residues.

•

Concerted domain movement was observed in both SpCas9 and FnCas9.

•

gRNA binding led to decrease in helicity and increase in structural transitions.

•

High binding affinity of FnCas9 with gRNA is mostly due to electrostatic interaction.

•

Arginine-helix has shown a pivotal role in the binding of gRNA in FnCas9 protein.

The introduction of CRISPR/Cas9 based gene editing has greatly accelerated therapeutic genome editing. However, the off-target DNA cleavage by CRISPR/Cas9 protein hampers its clinical translation, hindering its widespread use as a programmable genome editing tool. Although Cas9 variants with better mismatch discrimination have been developed, they have significantly lower rates of on-target DNA cleavage. Here, we have compared the dynamics of a more specific naturally occurring Cas9 from Francisella novicida (FnCas9) to the most widely used, SpCas9 protein. Long-scale atomistic MD simulation of free and gRNA bound forms of both the Cas9 proteins was performed, and their domain rearrangements and binding affinity with gRNA were compared to decipher the possible reason behind the enhanced specificity of FnCas9 protein. The greater binding affinity with gRNA, high domain electrostatics, and more volatility of FnCas9 than SpCas9 may explain its increased specificity and lower tolerance for mismatches.

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.

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 (K_m) 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.

Engineering and Design of Programmable Genome Editors

Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.

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.

Probing the Dynamics of Streptococcus pyogenes Cas9 Endonuclease Bound to the sgRNA Complex Using Hydrogen-Deuterium Exchange Mass Spectrometry

The Cas9 endonuclease is an essential component of the CRISPR-Cas-based genome editing tools. The attainment of high specificity and efficiency of Cas9 during targetted DNA cleavage is the main problem that limits the clinical application of the CRISPR-Cas9 system. A deep understanding of the Cas9 mechanism and its structural-functional relationships is required to develop strategies for precise gene editing. Here, we present the first attempt to describe the solution structure of Cas9 from S. pyogenes using hydrogen-deuterium exchange mass spectrometry (HDX-MS) coupled to molecular dynamics simulations. HDX data revealed multiple protein regions with deuterium uptake levels varying from low to high. By analysing the difference in relative deuterium uptake by apoCas9 and its complex with sgRNA, we identified peptides involved in the complex formation and possible changes in the protein conformation. The REC3 domain was shown to undergo the most prominent conformational change upon enzyme-RNA interactions. Detection of the HDX in two forms of the enzyme provided detailed information about changes in the Cas9 structure induced by sgRNA binding and quantified the extent of the changes. The study demonstrates the practical utility of HDX-MS for the elucidation of mechanistic aspects of Cas9 functioning.

Deciphering the Mechanism of Gilteritinib Overcoming Lorlatinib Resistance to the Double Mutant I1171N/F1174I in Anaplastic Lymphoma Kinase

Anaplastic lymphoma kinase (ALK) is validated as a therapeutic molecular target in multiple malignancies, such as non-small cell lung cancer (NSCLC). However, the feasibility of targeted therapies exerted by ALK inhibitors is inevitably hindered owing to drug resistance. The emergence of clinically acquired drug mutations has become a major challenge to targeted therapies and personalized medicines. Thus, elucidating the mechanism of resistance to ALK inhibitors is helpful for providing new therapeutic strategies for the design of next-generation drug. Here, we used molecular docking and multiple molecular dynamics simulations combined with correlated and energetical analyses to explore the mechanism of how gilteritinib overcomes lorlatinib resistance to the double mutant ALK I1171N/F1174I. We found that the conformational dynamics of the ALK kinase domain was reduced by the double mutations I1171N/F1174I. Moreover, energetical and structural analyses implied that the double mutations largely disturbed the conserved hydrogen bonding interactions from the hinge residues Glu1197 and Met1199 in the lorlatinib-bound state, whereas they had no discernible adverse impact on the binding affinity and stability of gilteritinib-bound state. These discrepancies created the capacity of the double mutant ALK I1171N/F1174I to confer drug resistance to lorlatinib. Our result anticipates to provide a mechanistic insight into the mechanism of drug resistance induced by ALK I1171N/F1174I that are resistant to lorlatinib treatment in NSCLC.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

Coordinated Actions of Cas9 HNH and RuvC Nuclease Domains Are Regulated by the Bridge Helix and the Target DNA Sequence

CRISPR-Cas systems are RNA-guided nucleases that provide adaptive immune protection in bacteria and archaea against intruding genomic materials. Cas9, a type-II CRISPR effector protein, is widely used for gene editing applications since a single guide RNA can direct Cas9 to cleave specific genomic targets. The conformational changes associated with RNA/DNA binding are being modulated to develop Cas9 variants with reduced off-target cleavage. Previously, we showed that proline substitutions in the arginine-rich bridge helix (BH) of Streptococcus pyogenes Cas9 (SpyCas9-L64P-K65P, SpyCas92Pro) improve target DNA cleavage selectivity. In this study, we establish that kinetic analysis of the cleavage of supercoiled plasmid substrates provides a facile means to analyze the use of two parallel routes for DNA linearization by SpyCas9: (i) nicking by HNH followed by RuvC cleavage (the TS (target strand) pathway) and (ii) nicking by RuvC followed by HNH cleavage (the NTS (nontarget strand) pathway). BH substitutions and DNA mismatches alter the individual rate constants, resulting in changes in the relative use of the two pathways and the production of nicked and linear species within a given pathway. The results reveal coordinated actions between HNH and RuvC to linearize DNA, which is modulated by the integrity of the BH and the position of the mismatch in the substrate, with each condition producing distinct conformational energy landscapes as observed by molecular dynamics simulations. Overall, our results indicate that BH interactions with RNA/DNA enable target DNA discrimination through the differential use of the parallel sequential pathways driven by HNH/RuvC coordination.

Atomic-scale insights into allosteric inhibition and evolutional rescue mechanism of Streptococcus thermophilus Cas9 by the anti-CRISPR protein AcrIIA6

CRISPR-Cas systems are prokaryotic adaptive immunity against invading phages and plasmids. Phages have evolved diverse protein inhibitors of CRISPR-Cas systems, called anti-CRISPR (Acr) proteins, to neutralize this CRISPR machinery. In response, bacteria have co-evolved Cas variants to escape phage's anti-CRISPR strategies, called anti-anti-CRISPR systems. Here we explore the anti-CRISPR allosteric inhibition and anti-anti-CRISPR rescue mechanisms between Streptococcus thermophilus Cas9 (St1Cas9) and the anti-CRISPR protein AcrIIA6 at the atomic level, by generating mutants of key residues in St1Cas9. Extensive unbiased molecular dynamics simulations show that the functional motions of St1Cas9 in the presence of AcrIIA6 differ substantially from those of St1Cas9 alone. AcrIIA6 binding triggers a shift of St1Cas9 conformational ensemble towards a less catalytically competent state; this state significantly compromises protospacer adjacent motif (PAM) recognition and nuclease activity by altering interdependently conformational dynamics and allosteric signals among nuclease domains, PAM-interacting (PI) regions, and AcrIIA6 binding motifs. Via in vitro DNA cleavage assays, we further elucidate the rescue mechanism of efficiently escaping AcrIIA6 inhibition harboring St1Cas9 triple mutations (G993K/K1008M/K1010E) in the PI domain and identify the evolutionary landscape of such mutational escape within species. Our results provide mechanistic insights into Acr proteins as natural brakes for the CRISPR-Cas systems and a promising potential for the design of allosteric Acr peptidomimetics.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn²⁺) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn²⁺ and Mg²⁺. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg²⁺, Ca²⁺) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg²⁺- and Ca²⁺-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg²⁺ binding relative to Ca²⁺ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn²⁺ is the most preferred relative to both Mg²⁺ and Ca²⁺ in the wild-type and mutant Cas1. Substitution of Mg²⁺ by Ca²⁺ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg²⁺ versus Ca²⁺) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

All-Atom Simulations Reveal a Key Interaction Network in the HLA-E/NKG2A/CD94 Immune Complex Fine-Tuned by the Nonameric Peptide

Natural killer (NK) cells, an important part of the innate immune system, can clear a wide variety of pathological challenges, including tumor, senescent, and virally infected cells. They express various activating and inhibitory receptors on their surface, and the balance of interactions between them and specific ligands displayed on the surface of target cells is critical for NK cell cytolytic function and target cell protection. The CD94/NKG2A heterodimer is one of the inhibitory receptors that interacts with its trimeric ligand consisting of HLA-E, β2m, and a nonameric peptide. Here, multi-microsecond-long all-atom molecular dynamics simulations of eight immune complexes elucidate the subtleties of receptor (NKG2A/CD94)-ligand (HLA-E/β2m/peptide) molecular recognition that mediate the NK cell protection from a geometric and energetic perspective. We identify key differences in the interactions between the receptor and ligand complexes, which are via an entangled network of hydrogen bonds fine-tuned by the ligand-specific nonameric peptide. We further reveal that the receptor protein NKG2A regulates the NK cell activity, while its CD94 partner forms the majority of the energetically important interactions with the ligand. This knowledge rationalizes the atomistic details of the fundamental NK cell protection mechanism and may enable a variety of opportunities in rational-based drug discovery for diverse pathologies including viral infections and cancer and elimination of senescent cells associated with potential treatment of many age-related diseases.

Nonameric Peptide Orchestrates Signal Transduction in the Activating HLA-E/NKG2C/CD94 Immune Complex as Revealed by All-Atom Simulations

The innate immune system's natural killer (NK) cells exert their cytolytic function against a variety of pathological challenges, including tumors and virally infected cells. Their activation depends on net signaling mediated via inhibitory and activating receptors that interact with specific ligands displayed on the surfaces of target cells. The CD94/NKG2C heterodimer is one of the NK activating receptors and performs its function by interacting with the trimeric ligand comprised of the HLA-E/β2m/nonameric peptide complex. Here, simulations of the all-atom multi-microsecond molecular dynamics in five immune complexes provide atomistic insights into the receptor-ligand molecular recognition, as well as the molecular events that facilitate the NK cell activation. We identify NKG2C, the HLA-Eα2 domain, and the nonameric peptide as the key elements involved in the molecular machinery of signal transduction via an intertwined hydrogen bond network. Overall, the study addresses the complex intricacies that are necessary to understand the mechanisms of the innate immune system.

Homology-based repair induced by CRISPR-Cas nucleases in mammalian embryo genome editing

Recent advances in genome editing, especially CRISPR-Cas nucleases, have revolutionized both laboratory research and clinical therapeutics. CRISPR-Cas nucleases, together with the DNA damage repair pathway in cells, enable both genetic diversification by classical non-homologous end joining (c-NHEJ) and precise genome modification by homology-based repair (HBR). Genome editing in zygotes is a convenient way to edit the germline, paving the way for animal disease model generation, as well as human embryo genome editing therapy for some life-threatening and incurable diseases. HBR efficiency is highly dependent on the DNA donor that is utilized as a repair template. Here, we review recent progress in improving CRISPR-Cas nuclease-induced HBR in mammalian embryos by designing a suitable DNA donor. Moreover, we want to provide a guide for producing animal disease models and correcting genetic mutations through CRISPR-Cas nuclease-induced HBR in mammalian embryos. Finally, we discuss recent developments in precise genome-modification technology based on the CRISPR-Cas system.

Establishing the allosteric mechanism in CRISPR-Cas9

Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.

CRISPR/Cas Genome Editing in Filamentous Fungi

The review describes the CRISPR/CAS system and its adaptation for the genome editing in filamentous fungi commonly used for production of enzyme complexes, enzymes, secondary metabolites, and other compounds used in industrial biotechnology and agriculture. In the second part of this review, examples of the CRISPR/CAS technology application for improving properties of the industrial strains of fungi from the Trichoderma, Aspergillus, Penicillium, and other genera are presented. Particular attention is given to the efficiency of genome editing, as well as system optimization for specific industrial producers.

Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules

Single-particle cryogenic electron microscopy (cryo-EM) has revolutionized the field of the structural biology, providing an access to the atomic resolution structures of large biomolecular complexes in their near-native environment. Today's cryo-EM maps can frequently reach the atomic-level resolution, while often containing a range of resolutions, with conformationally variable regions obtained at 6 Å or worse. Low resolution density maps obtained for protein flexible domains, as well as the ensemble of coexisting conformational states arising from cryo-EM, poses new challenges and opportunities for Molecular Dynamics (MD) simulations. With the ability to describe the biomolecular dynamics at the atomic level, MD can extend the capabilities of cryo-EM, capturing the conformational variability and predicting biologically relevant short-lived conformational states. Here, we report about the state-of-the-art MD procedures that are currently used to refine, reconstruct and interpret cryo-EM maps. We show the capability of MD to predict short-lived conformational states, finding remarkable confirmation by cryo-EM structures subsequently solved. This has been the case of the CRISPR-Cas9 genome editing machinery, whose catalytically active structure has been predicted through both long-time scale MD and enhanced sampling techniques 2 years earlier than cryo-EM. In summary, this contribution remarks the ability of MD to complement cryo-EM, describing conformational landscapes and relating structural transitions to function, ultimately discerning relevant short-lived conformational states and providing mechanistic knowledge of biological function.