"Virtual patch clamp analysis" for predicting the functional significance of pathogenic variants in sodium channels.

OBJECTIVE: Opening of voltage-gated sodium channels is crucial for neuronal depolarization. Proper channel opening and influx of Na+ through the ion pore, is dependent upon binding of Na+ ion to a specific amino-acid motif (DEKA) within the pore. In this study we used molecular dynamic simulations, an advanced bioinformatic tool, to research the dysfunction caused by pathogenic variants in SCN1a, SCN2a and SCN8a genes. METHOD: Molecular dynamic simulations were performed in six patients: three patients with Dravet syndrome (p.Gly177Ala,p.Ser259Arg and p.Met1267Ile, SCN1a), two patients with early onset drug resistant epilepsy(p.Ala263Val, SCN2a and p.Ile251Arg, SCN8a), and a patient with autism (p.Thr155Ala, SCN2a). After predicting the 3D-structure of mutated proteins by homology modeling, time dependent molecular dynamic simulations were performed, using the Schrödinger algorithm. The opening of the sodium channel, including the detachment of the sodium ion to the DEKA motif and pore diameter were assessed. Results were compared to the existent patch clamp analysis in four patients, and consistency with clinical phenotype was noted. RESULTS: The Na+ ion remained attached to DEKA filter longer when compared to wild type in the p.Gly177Ala, p.Ser259Arg,SCN1a, and p.Thr155Ala, SCN2a variants, consistent with loss-of-function. In contrast, it detached quicker from DEKA than wild type in the p.Ala263Val,SCN2a variant, consistent with gain-of-function. In the p.Met1267Ile,SCN1a variant, detachment from DEKA was quicker, but pore diameter decreased, suggesting partial loss-of-function. In the p.Leu251Arg,SCN8a variant, the pore remained opened longer when compared to wild type, consistent with a gain-of-function. The molecular dynamic simulation results were consistent with the existing patch-clamp analysis studies, as well as the clinical phenotype. SIGNIFICANCE: Molecular dynamic simulation can be useful in predicting pathogenicity of variants and the disease phenotype, and selecting targeted treatment based on channel dysfunction. Further development of these bioinformatic tools may lead to "virtual patch-clamp analysis".

Creatine transporter deficiency: Novel mutations and functional studies.

X-linked cerebral creatine deficiency (MIM 300036) is caused by deficiency of the creatine transporter encoded by the SLC6A8 gene. Here we report three patients with this condition from Israel. These unrelated patients were evaluated for global developmental delays and language apraxia. Borderline microcephaly was noted in one of them. Diagnosis was prompted by brain magnetic resonance imaging and spectroscopy which revealed normal white matter distribution, but absence of the creatine peak in all three patients. Biochemical testing indicated normal plasma levels of creatine and guanidinoacetate, but an increased urine creatine/creatinine ratio. The diagnosis was confirmed by demonstrating absent ([14])C-creatine transport in fibroblasts. Molecular studies indicated that the first patient is hemizygous for a single nucleotide change substituting a single amino acid (c.619 C > T, p.R207W). Expression studies in HeLa cells confirmed the causative role of the R207W substitution. The second patient had a three base pair deletion in the SLC6A8 gene (c.1222_1224delTTC, p.F408del) as well as a single base change (c.1254 + 1G > A) at a splicing site in the intron-exon junction of exon 8, the latter occurring de novo. The third patient, had a three base pair deletion (c.1006_1008delAAC, p.N336del) previously reported in other patients with creatine transporter deficiency. These three patients are the first reported cases of creatine transporter deficiency in Israel.