Abstract
G protein Coupled Receptors (GPCRs) are the largest family of cell surface receptors in humans. Somatic mutations in GPCRs are implicated in cancer progression and metastasis, but mechanisms are poorly understood. Emerging evidence implicates perturbation of intra-receptor activation pathway motifs whereby extracellular signals are transmitted intracellularly. Recently, sufficiently sensitive methodology was described to calculate structural strain as a function of missense mutations in AlphaFold-predicted model structures, which was extensively validated on experimental and predicted structural datasets. When paired with Molecular Dynamics (MD) simulations, these tools provide a facile approach to screen mutations in silico. We applied this framework to calculate the structural and dynamic effects of cancer-associated mutations in the chemokine receptor CCR3, a Class A GPCR involved in cancer and autoimmune disorders. Residue-residue contact scoring refined effective strain results, highlighting significant remodeling of inter- and intra-motif contacts along the highly conserved GPCR activation pathway network. We then integrated AlphaFold-derived predicted Local Distance Difference Test scores with per-residue Root Mean Square Fluctuations and activation pathway Contact Analysis (CONAN) from coarse grain MD simulations to identify statistically significant changes in receptor dynamics upon mutation. Finally, analysis of negative control mutants suggests false positive results in AlphaFold pipelines should be considered but can be mitigated with stricter control of statistical analysis. Our results indicate selected mutants influence structural plasticity of CCR3 related to ligand interaction, activation, and G protein coupling, using a framework that could be applicable to a wide range of biochemically relevant protein targets following further validation.