Ligand distances as key predictors of pathogenicity and function in NMDA receptors

Genetic variants in the genes GRIN1, GRIN2A, GRIN2B, and GRIN2D, which encode subunits of the N-methyl-D-aspartate receptor (NMDAR), have been associated with severe and heterogeneous neurologic and neurodevelopmental disorders, including early onset epilepsy, developmental and epileptic encephalopathy, intellectual disability, and autism spectrum disorders. Missense variants in these genes can result in gain or loss of the NMDAR function, requiring opposite therapeutic treatments. Computational methods that predict pathogenicity and molecular functional effects of missense variants are therefore crucial for therapeutic applications. We assembled 223 missense variants from patients, 631 control variants from the general population, and 160 missense variants characterized by electrophysiological readouts that show whether they can enhance or reduce the function of the receptor. This includes new functional data from 33 variants reported here, for the first time. By mapping these variants onto the NMDAR protein structures, we found that pathogenic/benign variants and variants that increase/decrease the channel function were distributed unevenly on the protein structure, with spatial proximity to ligands bound to the agonist and antagonist binding sites being a key predictive feature for both variant pathogenicity and molecular functional consequences. Leveraging distances from ligands, we developed two machine-learning based predictors for NMDA variants: a pathogenicity predictor which outperforms currently available predictors and the first molecular function (increase/decrease) predictor. Our findings can have direct application to patient care by improving diagnostic yield for genetic neurodevelopmental disorders and by guiding personalized treatment informed by the knowledge of the molecular disease mechanism.