Investigation of RNA-binding protein NOVA1 in silico: Comparison of the modern human V197 with the archaic I197 variant present in Neanderthals

By comparing the high-coverage archaic genome sequences to those of modern humans, specific genetic differences have been identified. For example, a human-specific substitution has been found in neuro-oncological ventral antigen 1 (NOVA1) - an RNA-binding protein that regulates the alternative splicing of neuronal pre-mRNA. The amino acid substitution results in an isoleucine-to-valine change at position 197 in NOVA1 (archaic: I197, modern human: V197). Previous studies have utilised gene editing technology to compare the archaic and modern human forms of NOVA1 in cortical organoids, however, the structural and molecular details require further investigation. Using an in silico approach, the modern human (WT) and archaic (V197I) structures of NOVA1 were generated. Moreover, the structure of NOVA1 containing a glycine-to-valine substitution at position 68 (G68V), which occurs at the RNA-binding interface, was examined for comparison. Protein-RNA docking was subsequently performed to model the interaction of NOVA1 variants with RNA and the complexes were evaluated further using classical molecular dynamics (MD) simulations. Based on the MM-PBSA analysis, the binding free energies were similar between the WT (-956.8 ± 32.6 kcal/mol), V197I (-975.4 ± 65.6 kcal/mol), and G68V (-946.7 ± 34.3 kcal/mol) complexes. The findings highlight the binding and stability of protein-RNA complexes with only modest structural changes observed in the archaic and G68V variants compared to the WT NOVA1 protein. Further clarification is required to enhance our understanding of the impact of NOVA1 mutations on alternative splicing and disease development. In particular, delineating the effect of multiple mutations in the NOVA1 gene is of importance.