A Computational Model for Nme1Cas9 HNH Activation Driven by Dynamic Interface Engineering at Residues S593 and W596

Nme1Cas9 is an encouraging genome-editing tool with high fidelity and compactness, but its applications are limited by poor catalytic efficiency compared with SpyCas9. Understanding the dynamic activation mechanism of the HNH nuclease domain is the key to breaking the kinetic bottleneck. Here, we integrated Steered Molecular Dynamics (SMD) with the Traveling-Salesman-based automated Path Searching (TAPS) algorithm to reconstruct the atomic-level activation landscape of the L1-HNH module. The simulations suggest a complex “Lifting-Rearrangement-Sliding” pathway, revealing the critical role of a “Backbone Sliding” conformation; in this step, the HNH domain rotates across the R-loop surface. A thermodynamic analysis using free energy decomposition by MM/PBSA indicates that the intrinsic instability of the wild-type HNH/R-loop interface constitutes the predominant energetic barrier. Hyperactive variants (S593Q/W596K and S593Q/W596R) can overcome this barrier by substantially increasing binding affinity to the R-loop through a “Geometry–Electrostatics Synergism”: S593Q improves interfacial proximity, whereas W596K/R acts as an “Electrostatic Anchor.” The results of unbiased MD simulations demonstrate that strengthened interfacial interactions effectively promote spontaneous conformational drift toward the activated state. This computational study proposes a novel in silico model for “Dynamic Interface Engineering” in which reinforcing transient interfacial contacts during conformational sliding can be an effective strategy in developing high-efficiency CRISPR-Cas effectors.