CRISPR-Cas9 conformational activation as elucidated from enhanced molecular simulations

Multiple Parameter Replica Exchange Gaussian Accelerated Molecular Dynamics for Enhanced Sampling and Free Energy Calculation of Biomolecular Systems

This study introduces a novel method named multiple parameter replica exchange Gaussian accelerated molecular dynamics (MP-Rex-GaMD), building on the Gaussian accelerated molecular dynamics (GaMD) algorithm. GaMD enhances sampling and retrieves free energy information for biomolecular systems by adding a harmonic boost potential to smooth the potential energy surface without the need for predefined reaction coordinates. Our innovative approach advances the acceleration power and energetic reweighting accuracy of GaMD by incorporating a replica exchange algorithm that enables the exchange of multiple parameters, including the GaMD boost parameters of force constant and energy threshold, as well as temperature. Applying MP-Rex-GaMD to the three model systems of dialanine, chignolin, and HIV protease, we demonstrate its superior capability over conventional molecular dynamics and GaMD simulations in exploring protein conformations and effectively navigating various biomolecular states across energy barriers. MP-Rex-GaMD allows users to accurately map free energy landscapes through energetic reweighting, capturing the ensemble of biomolecular states from low-energy conformations to rare high-energy transitions within practical computational time scales.

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.

Activation of Polycystin-1 Signaling by Binding of Stalk-derived Peptide Agonists

Polycystin-1 (PC1) is the membrane protein product of the PKD1 gene whose mutation is responsible for 85% of the cases of autosomal dominant polycystic kidney disease (ADPKD). ADPKD is primarily characterized by the formation of renal cysts and potential kidney failure. PC1 is an atypical G protein-coupled receptor (GPCR) consisting of 11 transmembrane helices and an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal (NTF) and membrane-embedded C-terminal (CTF) fragments. Recently, signaling activation of the PC1 CTF was shown to be regulated by a stalk tethered agonist (TA), a distinct mechanism observed in the adhesion GPCR family. A novel allosteric activation pathway was elucidated for the PC1 CTF through a combination of Gaussian accelerated molecular dynamics (GaMD), mutagenesis and cellular signaling experiments. Here, we show that synthetic, soluble peptides with 7 to 21 residues derived from the stalk TA, in particular, peptides including the first 9 residues (p9), 17 residues (p17) and 21 residues (p21) exhibited the ability to re-activate signaling by a stalkless PC1 CTF mutant in cellular assays. To reveal molecular mechanisms of stalk peptide-mediated signaling activation, we have applied a novel Peptide GaMD (Pep-GaMD) algorithm to elucidate binding conformations of selected stalk peptide agonists p9, p17 and p21 to the stalkless PC1 CTF. The simulations revealed multiple specific binding regions of the stalk peptide agonists to the PC1 protein including an "intermediate" bound yet inactive state. Our Pep-GaMD simulation findings were consistent with the cellular assay experimental data. Binding of peptide agonists to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of the PC1 CTF signaling activation mechanism. Using sequence covariation analysis of PC1 homologs, we further showed that the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. Therefore, structural dynamic insights into the mechanisms of PC1 activation by stalk-derived peptide agonists have enabled an in-depth understanding of PC1 signaling. They will form a foundation for development of PC1 as a therapeutic target for the treatment of ADPKD.

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

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

Transient interactions modulate the affinity of NF-κB transcription factors for DNA

The dimeric nuclear factor kappa B (NF-κB) transcription factors (TFs) regulate gene expression by binding to a variety of κB DNA elements with conserved G:C-rich flanking sequences enclosing a degenerate central region. Toward defining mechanistic principles of affinity regulated by degeneracy, we observed an unusual dependence of the affinity of RelA on the identity of the central base pair, which appears to be noncontacted in the complex crystal structures. The affinity of κB sites with A or T at the central position is ~10-fold higher than with G or C. The crystal structures of neither the complexes nor the free κB DNAs could explain the differences in affinity. Interestingly, differential dynamics of several residues were revealed in molecular dynamics simulation studies, where simulation replicates totaling 148 μs were performed on NF-κB:DNA complexes and free κB DNAs. Notably, Arg187 and Arg124 exhibited selectivity in transient interactions that orchestrated a complex interplay among several DNA-interacting residues in the central region. Binding and simulation studies with mutants supported these observations of transient interactions dictating specificity. In combination with published reports, this work provides insights into the nuanced mechanisms governing the discriminatory binding of NF-κB family TFs to κB DNA elements and sheds light on cancer pathogenesis of cRel, a close homolog of RelA.

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.

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.

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.

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.

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.

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

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.

Structure and dynamics of pteridine reductase 1: the key phenomena relevant to enzyme function and drug design

Pteridine reductase 1 (PTR1) is a folate and pterin pathway enzyme unique for pathogenic trypanosomatids. As a validated drug target, PTR1 has been the focus of recent research efforts aimed at finding more effective treatments against human parasitic diseases such as leishmaniasis or sleeping sickness. Previous PTR1-centered structural studies highlighted the enzyme characteristics, such as flexible regions around the active site, highly conserved structural waters, and species-specific differences in pocket properties and dynamics, which likely impacts the binding of natural substrates and inhibitors. Furthermore, several aspects of the PTR1 function, such as the substrate inhibition phenomenon and the level of ligand binding cooperativity in the enzyme homotetramer, likely related to the global enzyme dynamics, are poorly known at the molecular level. We postulate that future drug design efforts could greatly benefit from a better understanding of these phenomena through studying both the local and global PTR1 dynamics. This review highlights the key aspects of the PTR1 structure and dynamics relevant to structure-based drug design that could be effectively investigated by modeling approaches. Particular emphasis is given to the perspective of molecular dynamics, what has been accomplished in this area to date, and how modeling could impact the PTR1-targeted drug design in the future.

Application of CRISPR-Cas9 genome editing technology in various fields: A review

CRISPR-Cas9 has emerged as a revolutionary tool that enables precise and efficient modifications of the genetic material. This review provides a comprehensive overview of CRISPR-Cas9 technology and its applications in genome editing. We begin by describing the fundamental principles of CRISPR-Cas9 technology, explaining how the system utilizes a single guide RNA (sgRNA) to direct the Cas9 nuclease to specific DNA sequences in the genome, resulting in targeted double-stranded breaks. In this review, we provide in-depth explorations of CRISPR-Cas9 technology and its applications in agriculture, medicine, environmental sciences, fisheries, nanotechnology, bioinformatics, and biotechnology. We also highlight its potential, ongoing research, and the ethical considerations and controversies surrounding its use. This review might contribute to the understanding of CRISPR-Cas9 technology and its implications in various fields, paving the way for future developments and responsible applications of this transformative technology.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.

Molecular dynamics study reveals key disruptors of MEIG1-PACRG interaction

Interactions between the meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) protein are critical in the formation of mature sperm cells. Targeting either MEIG1 or PACRG protein could be a contraceptive strategy. The W50A and Y68A mutations on MEIG1 are known to interrupt the MEIG1-PACRG interactions resulting in defective sperm cells. However, the details about how the mutants disrupt the protein-protein binding are not clear. In this study, we reveal insights on MEIG1 and PACRG protein dynamics by applying Gaussian-accelerated molecular dynamics (GaMD) simulations and post-GaMD analysis. Our results show that the mutations destabilize the protein-protein interfacial interaction. The effect of the Y68A mutation is more significant than W50A as Y68 forms stronger polar interactions with PACRG. Because both human and mouse models demonstrate similar dynamic properties, the findings from mouse proteins can be applied to the human system. Moreover, we report a potential ligand binding pocket on the MEIG1 and PACRG interaction surface that could be a target for future drug design to inhibit the MEIG1-PACRG interaction. PACRG shows more qualified pockets along the protein-protein interface, implying that it is a better target than MEIG1. Our work provides a fundamental understanding of MEIG1 and PACRG protein dynamics, paving the way for drug discovery in male-based contraception.

Machines on Genes through the Computational Microscope

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

The Electronic Structure of Genome Editors from the First Principles

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

Integrating comparative modeling and accelerated simulations reveals conformational and energetic basis of actomyosin force generation

Muscle contraction is performed by arrays of contractile proteins in the sarcomere. Serious heart diseases, such as cardiomyopathy, can often be results of mutations in myosin and actin. Direct characterization of how small changes in the myosin-actin complex impact its force production remains challenging. Molecular dynamics (MD) simulations, although capable of studying protein structure-function relationships, are limited owing to the slow timescale of the myosin cycle as well as a lack of various intermediate structures for the actomyosin complex. Here, employing comparative modeling and enhanced sampling MD simulations, we show how the human cardiac myosin generates force during the mechanochemical cycle. Initial conformational ensembles for different myosin-actin states are learned from multiple structural templates with Rosetta. This enables us to efficiently sample the energy landscape of the system using Gaussian accelerated MD. Key myosin loop residues, whose substitutions are related to cardiomyopathy, are identified to form stable or metastable interactions with the actin surface. We find that the actin-binding cleft closure is allosterically coupled to the myosin motor core transitions and ATP-hydrolysis product release from the active site. Furthermore, a gate between switch I and switch II is suggested to control phosphate release at the prepowerstroke state. Our approach demonstrates the ability to link sequence and structural information to motor functions.

Elucidating the molecular determinants in the process of gastrin C-terminal pentapeptide amide end activating cholecystokinin 2 receptor by Gaussian accelerated molecular dynamics simulations

Gastrin plays important role in stimulating the initiation and development of many gastrointestinal diseases through interacting with the cholecystokinin 2 receptor (CCK2R). The smallest bioactive unit of gastrin activating CCK2R is the C-terminal tetrapeptide capped with an indispensable amide end. Understanding the mechanism of this smallest bioactive unit interacting with CCK2R on a molecular basis could provide significant insights for designing CCK2R antagonists, which can be used to treat gastrin-related diseases. To this end, we performed extensive Gaussian accelerated molecular dynamics simulations to investigate the interaction between gastrin C-terminal pentapeptide capped with/without amide end and CCK2R. The amide cap influences the binding modes of the pentapeptide with CCK2R by weakening the electrostatic attractions between the C-terminus of the pentapeptide and basic residues near the extracellular domain in CCK2R. The C-terminus with the amide cap penetrates into the transmembrane domain of CCK2R while floating at the extracellular domain without the amide cap. Different binding modes induced different conformational dynamics of CCK2R. Residue pairs in CCK2R had stronger correlated motions when binding with the amidated pentapeptide. Key residues and interactions important for CCK2R binding with the amidated pentagastrin were also identified. Our results provide molecular insights into the determinants of the bioactive unit of gastrin activating CCK2R, which would be of great help for the design of CCK2R antagonists.

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.

Allosteric Boost by TAB1 on the TAK1 Kinase Favorably Sculpts the Thermodynamic Landscape of Activation

The intricate mechanisms of allosteric regulation in kinases are of general interest to the scientific community for potential therapeutic implications. However, the diversity among kinases and their regulatory routes requires a case-by-case study to widen the repertoire of known mechanisms. The present study achieves this by understanding TAK1 kinase activation by TAB1 as a model phenomenon for the first time. Despite the known capacity of TAK1 to switch between its inactive ("DFG-out") and active-like ("DFG-in") conformations, the questionable role of TAB1 in offering an energetic favor to this has been addressed here using sequential combination of enhanced sampling methods like targeted molecular dynamics (TMD) and Gaussian accelerated molecular dynamics (GaMD). It reveals how a minimal domain of TAB1 sufficiently acts like a "catalytic gear" by favorably sculpting TAK1's thermodynamic landscape (potential of mean force in 2D) that accelerates "in"-"out" conformational switching of the conserved DFG motif. Standard molecular dynamics simulations (∼5 μs) reveal that TAB1 fascinatingly exploits the "lever-like" αF helix of TAK1 kinase domain to remotely propel the DFG motif via subtle helical "unfolding-folding" modifications within the kinase activation loop. The presence of two charged residues on terminal poles of αF helix imparts it, with this unique "lever-like" utility, and this turns out to be one important signature of co-evolution between TAK1 and TAB1. The entire mechanism of TAB1's impact transduction, which is found to be analogous to the moves in the popular "Chinese checker" game, gives a clear proof of the "dynamics-driven allostery" concept in kinases. The findings further benchmark TAK1's known autophosphorylation capacity. A novel insight into kinase allostery is thus provided, which potentiates investigation of similar capacities in other kinases.

Constructing next-generation CRISPR-Cas tools from structural blueprints

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

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.

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 Mg²⁺, unveiling its cardinal structural role. This activated Mg²⁺-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.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Thermodynamic Swings: How Ideal Complex of Cas9-RNA/DNA Forms

Most processes of the recognition and formation of specific complexes in living systems begin with collisions in solutions or quasi-solutions. Then, the thermodynamic regulation of complex formation and fine tuning of complexes come into play. Precise regulation is very important in all cellular processes, including genome editing using the CRISPR-Cas9 tool. The Cas9 endonuclease is an essential component of the CRISPR-Cas-based genome editing systems. The attainment of high-specificity and -efficiency Cas9 during targeted DNA cleavage is the main problem that limits the practical application of the CRISPR-Cas9 system. In this study, we analyzed the thermodynamics of interaction of a complex's components of Cas9-RNA/DNA through experimental and computer simulation methods. We found that there is a small energetic preference during Cas9-RNA/DNA formation from the Cas9-RNA and DNA/DNA duplex. The small difference in binding energy is relevant for biological interactions and could be part of the sequence-specific recognition of double-stranded DNA by the CRISPR-Cas9 system.

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.

Implementing accelerated dynamics to unravel the effects of high-fidelity Cas9 mutants on target DNA and guide RNA hybrid stability

The clustered regularly interspersed short palindromic repeats (CRISPR) and its associated nuclease (Cas9) offers a unique and easily reprogrammable system for editing eukaryotic genomes. Cas9 is guided to the target by an RNA strand, and precise edits are created by introducing double-stranded breaks. However, nuclease activity of Cas9 is also triggered at other sites other than the target sit, which is a major limitation for various applications. Cas9 variants have been designed to improve the efficacy of the tool by introducing certain mutations. However, the on-target activity of such Cas9 variants is often seen as compromised. Hence, understanding the sub-molecular differences in the variants is essential to elucidate the factors that contribute to efficiency. The study reveals distortions in the PAM-distal regions of the nucleic hybrids as well as changes in the interactions between the Cas9 variants and RNA-DNA hybrid, contributing to the explanation for differences in on-target activity.Communicated by Ramaswamy H. Sarma.

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.

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 and Computational Strategies to Increase the Efficiency of CRISPR-Based Techniques

The prokaryote-derived Clustered Regularly Interspaced Palindromic Repeats (CRISPR)/Cas mediated gene editing tools have revolutionized our ability to precisely manipulate specific genome sequences in plants and animals. The simplicity, precision, affordability, and robustness of this technology have allowed a myriad of genomes from a diverse group of plant species to be successfully edited. Even though CRISPR/Cas, base editing, and prime editing technologies have been rapidly adopted and implemented in plants, their editing efficiency rate and specificity varies greatly. In this review, we provide a critical overview of the recent advances in CRISPR/Cas9-derived technologies and their implications on enhancing editing efficiency. We highlight the major efforts of engineering Cas9, Cas12a, Cas12b, and Cas12f proteins aiming to improve their efficiencies. We also provide a perspective on the global future of agriculturally based products using DNA-free CRISPR/Cas techniques. The improvement of CRISPR-based technologies efficiency will enable the implementation of genome editing tools in a variety of crop plants, as well as accelerate progress in basic research and molecular breeding.

Application of CRISPR-Cas9 System to Study Biological Barriers to Drug Delivery

In recent years, sequence-specific clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) systems have been widely used in genome editing of various cell types and organisms. The most developed and broadly used CRISPR-Cas system, CRISPR-Cas9, has benefited from the proof-of-principle studies for a better understanding of the function of genes associated with drug absorption and disposition. Genome-scale CRISPR-Cas9 knockout (KO) screen study also facilitates the identification of novel genes in which loss alters drug permeability across biological membranes and thus modulates the efficacy and safety of drugs. Compared with conventional heterogeneous expression models or other genome editing technologies, CRISPR-Cas9 gene manipulation techniques possess significant advantages, including ease of design, cost-effectiveness, greater on-target DNA cleavage activity and multiplexing capabilities, which makes it possible to study the interactions between membrane proteins and drugs more accurately and efficiently. However, many mechanistic questions and challenges regarding CRISPR-Cas9 gene editing are yet to be addressed, ranging from off-target effects to large-scale genetic alterations. In this review, an overview of the mechanisms of CRISPR-Cas9 in mammalian genome editing will be introduced, as well as the application of CRISPR-Cas9 in studying the barriers to drug delivery.

Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor

Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1-4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.

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.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Enhanced sampling without borders: on global biasing functions and how to reweight them

Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

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.

A Review: Computational Approaches to Design sgRNA of CRISPR-Cas9

:

Clustered regularly interspaced short palindromic repeats along with CRISPR-associated protein mechanisms preserve the memory of previous experiences with DNA invaders, in particular spacers that are embedded in CRISPR arrays between coordinate repeats. There has been a fast progression in the comprehension of this immune system and its implementations; however, there are numerous points of view that anticipate explanations to make the field an energetic research zone. The efficiency of CRISPR-Cas depends upon well-considered single guide RNA; for this purpose, many bioinformatics methods and tools are created to support the design of greatly active and precise single guide RNA. Insilico single guide RNA architecture is a crucial point for effective gene editing by means of the CRISPR technique. Persistent attempts have been made to improve in-silico single guide RNA formulation having great on-target effectiveness and decreased off-target effects. This review offers a summary of the CRISPR computational tools to help different researchers pick a specific tool for their work according to pros and cons, along with new thoughts to make new computational tools to overcome all existing limitations.

Mechanism of RNA recognition by a Musashi RNA-binding protein

The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.

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

Markov State Models and Molecular Dynamics Simulations Provide Understanding of the Nucleotide-Dependent Dimerization-Based Activation of LRRK2 ROC Domain

Mutations in leucine-rich repeat kinase 2 (LRRK2) are recognized as the most frequent cause of Parkinson's disease (PD). As a multidomain ROCO protein, LRRK2 is characterized by the presence of both a Ras-of-complex (ROC) GTPase domain and a kinase domain connected through the C-terminal of an ROC domain (COR). The bienzymatic ROC-COR-kinase catalytic triad indicated the potential role of GTPase domain in regulating kinase activity. However, as a functional GTPase, the detailed intrinsic regulation of the ROC activation cycle remains poorly understood. Here, combining extensive molecular dynamics simulations and Markov state models, we disclosed the dynamic structural rearrangement of ROC's homodimer during nucleotide turnover. Our study revealed the coupling between dimerization extent and nucleotide-binding state, indicating a nucleotide-dependent dimerization-based activation scheme adopted by ROC GTPase. Furthermore, inspired by the well-known R1441C/G/H PD-relevant mutations within the ROC domain, we illuminated the potential allosteric molecular mechanism for its pathogenetic effects through enabling faster interconversion between inactive and active states, thus trapping ROC in a prolonged activated state, while the implicated allostery could provide further guidance for identification of regulatory allosteric pockets on the ROC complex. Our investigations illuminated the thermodynamics and kinetics of ROC homodimer during nucleotide-dependent activation for the first time and provided guidance for further exploiting ROC as therapeutic targets for controlling LRRK2 functionality in PD treatment.