Ensemble-Based Binding Free Energy Profiling and Network Analysis of the KRAS Interactions with DARPin Proteins Targeting Distinct Binding Sites: Revealing Molecular Determinants and Universal Architecture of Regulatory Hotspots and Allosteric Binding

TL;DR

This study explores how engineered DARPin proteins interact with the KRAS oncoprotein, revealing shared and unique molecular mechanisms that could aid in developing targeted cancer therapies.

Contribution

The study introduces a multi-scale approach combining simulations and network analysis to uncover universal and mechanism-specific hotspots in KRAS-DARPin interactions.

Findings

KRAS residues I36, Y40, M67, and H95 are consistently critical for binding stability across DARPin systems.

Allosteric binders K13 and K19 stabilize a KRAS pocket and rewire communication without full immobilization.

A unifying allosteric architecture is preserved across complexes, mirroring KRAS's intrinsic communication framework.

Abstract

KRAS is a pivotal oncoprotein that regulates cell proliferation and survival through interactions with downstream effectors such as RAF1. Despite significant advances in understanding KRAS biology, the structural and dynamic mechanisms of KRAS allostery remain poorly understood. In this study, we employ microsecond molecular dynamics simulations, mutational scanning, and binding free energy calculations together with dynamic network modeling to dissect how engineered DARPin proteins K27, K55, K13, and K19 engage KRAS through diverse molecular mechanisms ranging from effector mimicry to conformational restriction and allosteric modulation. Mutational scanning across all four DARPin systems identifies a core set of evolutionarily constrained residues that function as universal hotspots in KRAS recognition. KRAS residues I36, Y40, M67, and H95 consistently emerge as critical contributors to binding stability. Binding free energy computations show that, despite similar binding modes, K27 relies heavily on electrostatic contributions from major binding hotspots while K55 exploits a dense hydrophobic cluster enhancing its effector-mimetic signature. The allosteric binders K13 and K19, by contrast, stabilize a KRAS-specific pocket in the α3–loop–α4 motif, introducing new hinges and bottlenecks that rewire the communication architecture of KRAS without full immobilization. Network-based analysis reveals a strikingly consistent theme: despite their distinct mechanisms of recognition, all systems engage a unifying allosteric architecture that spans multiple functional motifs. This architecture is not only preserved across complexes but also mirrors the intrinsic communication framework of KRAS itself, where specific residues function as central hubs transmitting conformational changes across the protein. By integrating dynamic profiling, energetic mapping, and network modeling, our study provides a multi-scale mechanistic roadmap for targeting KRAS, revealing how engineered proteins can exploit both conserved motifs and isoform-specific features to enable precision modulation of KRAS signaling in oncogenic contexts.

Figures