Deciphering the Allosteric Inhibition Mechanism of CXCL12 by Gypsogenin: A Molecular Dynamics and Free Energy Landscape Study

The CXCL12/CXCR4 signaling axis plays a pivotal role in driving cancer metastasis, making it a high-priority therapeutic target. However, due to its flat, highly flexible, and solvent-exposed interaction interface, CXCL12 has historically been classified as an elusive or undruggable ligand. Herewith, we challenged this paradigm by investigating the antagonistic potential and molecular mechanism of Gypsogenin, a bioactive natural pentacyclic triterpenoid, against human CXCL12 using comprehensive 100 ns all-atom molecular dynamics (MD) simulations. Interface overlap analysis revealed that Gypsogenin establishes a remarkable 57% structural overlap with the native CXCR4-binding footprint, primarily anchoring to the critical Lys24 hotspot with a high interaction frequency (40.5%). This interaction effectively masks the electropositive patch required to recruit the sulfotyrosine (sY7/sY12) residues of the receptor. Furthermore, Dynamic Cross-Correlation Map (DCCM) analysis demonstrated that Gypsogenin binding restricts the chemokine's structural plasticity, locking the core functional domain (residues 10–50) into a highly rigid, synchronized motion that prevents induced fit transitions. Essential dynamics coupled with Free Energy Landscape (FEL) mapping confirmed that the complex is thermodynamically trapped within a single, deep global energy minimum valley, ensuring a severe kinetic barrier against dissociation. Taken together, our findings provide the first detailed structural and thermodynamic evidence of Gypsogenin as a powerful natural modulator capable of directly neutralizing CXCL12, offering a new blueprint for ligand-targeted anti-metastatic drug discovery.