Exploring the Dynamic Interplay of Deleterious Variants on the RAF1–RAP1A Binding in Cancer: Conformational Analysis, Binding Free Energy, and Essential Dynamics
submitted on 2025-10-08, 10:23 and posted on 2025-10-08, 10:26authored byAbbas Khan, Syed Shujait Ali, Muhammad Ammar Zahid, Shahenda Salah Abdelsalam, Noorah Albekairi, Raed M. Al‐Zoubi, Mohanad Shkoor, Dong‐Qing Wei, Abdelali Agouni
<p dir="ltr">The RAF1–RAP1A interaction activates the MAPK/ERK pathway which is very crucial in the carcinogenesis process. This protein complex influences tumor formation, proliferation, and metastasis. Understanding aberrant interactions driven by clinical mutations is vital for targeted therapies. Hence, the current study focuses on the screening of clinically reported substitutions in the <i>RAF1</i> and <i>RAP1A</i> genes using predictive algorithms integrated with all‐atoms simulation, essential dynamics, and binding free energy methods. Survival analysis results revealed a strong association between <i>RAF1</i> and <i>RAP1A</i> expression levels and diminished survival rates in cancer patients across different cancer types. Integrated machine learning algorithms showed that among the 134 mutations reported for these 2 proteins, only 13 and 35 were classified as deleterious mutations in <i>RAF1</i> and <i>RAP1P</i>, respectively. Moreover, one mutation in RAF1 reported elevated levels of binding between <i>RAF1</i> and <i>RAP1P</i> while in <i>RAP1A</i>, 7 mutations were reported to increase the binding affinity. The high‐binding mutations, P34Q and V60F, were subjected to protein–protein coupling which confirmed the increase in the binding affinity. Wild‐type and mutant RAF1–RAP1P bound complexes were subjected to molecular simulation investigation, revealing enhanced structural stability, increased compactness, and stabilized residue fluctuations of the mutant systems in contrast to the wild‐type. In addition, hydrogen bonding analysis revealed a variation in the binding paradigm which further underscores the impact of these substitutions on the coupling of <i>RAF1</i> and <i>RAP1A</i>. Principal component analysis (PCA) and free energy landscape (FEL) evaluation further determined dynamical variations in the wild‐type and mutant complexes. Finally, the Gibbs free energy for each complex was estimated and found to be −71.94 ± 0.38 kcal/mol for the wild‐type, −95.57 ± 0.37 kcal/mol for the V60F, and −85.76 ± 0.72 kcal/mol for P34Q complex. These findings confirm the effect of these variants on increasing the binding affinity of RAF1 to RAP1P. These mutations can therefore be targeted for cancer therapy to modulate the activity of the MAPK/ERK signaling pathway.</p><h2>Other Information</h2><p dir="ltr">Published in: Proteins: Structure, Function, and Bioinformatics<br>License: <a href="http://creativecommons.org/licenses/by/4.0/" target="_blank">http://creativecommons.org/licenses/by/4.0/</a><br>See article on publisher's website: <a href="https://dx.doi.org/10.1002/prot.26759" target="_blank">https://dx.doi.org/10.1002/prot.26759</a></p>
Funding
Open Access funding provided by the Qatar National Library.
Qatar National Research Fund (ARG01- 0601-230451).