Water and chloride as allosteric inhibitors in WNK kinase osmosensing.

Previous observations showed that chloride and osmotic stress regulate the autophosphorylation and activity of the kinase domains of WNK1 and WNK3. Further, prior crystallography on the asymmetric dimeric of the unphosphorylated WNK1 kinase domain (WNK1/S382A, WNK1/SA) revealed conserved waters in the active site. Here we show by crystallography that PEG400 applied to crystals of dimeric WNK1/SA grown in space group P1 induces de-dimerization with a change in space group to P2 1 . Both the conserved waters, referred to here as conserved water network 1 (CWN1) and the chloride binding site are disrupted by PEG400. CWN1 is surrounded and stabilized by a pan-WNK-conserved cluster of charged residues. Here we mutagenized these charges in WNK3 to probe the importance of the CWN1 to WNK regulation. Two mutations at E314 in the Activation Loop (WNK3/E314Q and WNK3/E314A) enhanced activity, consistent with the idea that the CWN1 is inhibitory. Mutations of other residues in the cluster had similar or less activity than wild-type. PEG400 activation of WNK3 was not significantly reduced in the point mutants tested. The crystallographic and assay data support a role for CWN1 and the charged cluster in stabilizing an inactive configuration of WNKs and suggest that water functions as an allosteric inhibitor of WNKs.

Hydrostatic Pressure Sensing by WNK kinases.

Previous study has demonstrated that the WNK kinases 1 and 3 are direct osmosensors consistent with their established role in cell-volume control. WNK kinases may also be regulated by hydrostatic pressure. Hydrostatic pressure applied to cells in culture with N2 gas or to Drosophila Malpighian tubules by centrifugation induces phosphorylation of downstream effectors of endogenous WNKs. In vitro, the autophosphorylation and activity of the unphosphorylated kinase domain of WNK3 (uWNK3) is enhanced to a lesser extent than in cells by 190 kPa applied with N2 gas. Hydrostatic pressure measurably alters the structure of uWNK3. Data from size exclusion chromatography in line with multi-angle light scattering (SEC-MALS), SEC alone at different back pressures, analytical ultracentrifugation (AUC), NMR, and chemical crosslinking indicate a change in oligomeric structure in the presence of hydrostatic pressure from a WNK3 dimer to a monomer. The effects on the structure are related to those seen with osmolytes. Potential mechanisms of hydrostatic pressure activation of uWNK3 and the relationships of pressure activation to WNK osmosensing are discussed.

Identification of a Class of WNK Isoform-Specific Inhibitors Through High-Throughput Screening.

Introduction: WNK [with no lysine (K)] kinases are serine/threonine kinases associated with familial hyperkalemic hypertension (FHHt). WNKs are therapeutic targets for blood pressure regulation, stroke and several cancers including triple negative breast cancer and glioblastoma. Here, we searched for and characterized novel WNK kinase inhibitors. Methods: We used a ~210,000-compound library in a high-throughput screen, re-acquisition and assay, commercial specificity screens and crystallography to identify WNK-isoform-selective inhibitors. Results: We identified five classes of compounds that inhibit the kinase activity of WNK1: quinoline compounds, halo-sulfones, cyclopropane-containing thiazoles, piperazine-containing compounds, and nitrophenol-derived compounds. The compounds are strongly pan-WNK selective, inhibiting all four WNK isoforms. A class of quinoline compounds was identified that further shows selectivity among the WNK isoforms, being more potent toward WNK3 than WNK1. The crystal structure of the quinoline-derived SW120619 bound to the kinase domain of WNK3 reveals active site binding, and comparison to the WNK1 structure reveals the potential origin of isoform specificity. Discussion: The newly discovered classes of compounds may be starting points for generating pharmacological tools and potential drugs treating hypertension and cancer.

Synthesis and Structural Characterization of Novel Trihalo-sulfone Inhibitors of WNK1.

With No lysine (K) [WNK] kinases are structurally unique serine/threonine protein kinases that have therapeutic potential for blood pressure regulation and cancer. A novel class of trihalo-sulfone compounds was identified by high-throughput screening. Trihalo-sulfone 1 emerged as an effective inhibitor of WNK1 with an IC50 value of 1.6 µM. Herein, we define chemical features necessary for inhibition of WNK1 using chemical synthesis and X-ray crystallography. Analogues that probed the role of specific functional groups to the inhibitory activity were synthesized. X-ray structures of trihalo-sulfone 1 and a second trihalo-sulfone 23 bound to WNK1 revealed active site binding to two of the three previously defined canonical inhibitor binding pockets as well as a novel binding site for the trihalo-sulfone moiety. The elucidation of these novel interaction sites may allow for the strategic design of even more selective and potent WNK inhibitors.

Cerebral creatine deficiency disorders - A clinical, genetic and follow up study from India.

INTRODUCTION: Cerebral creatine deficiency syndromes (CCDS) are a group of potentially treatable neurometabolic disorders. The clinical, genetic profile and follow up outcome of Indian CCDS patients is presented. MATERIALS AND METHODS: This was a retrospective cohort of CCDS patients seen over six-years. Diagnosis was based either on low creatine peak on proton magnetic resonance spectroscopy (MRS) and/or genetic evaluation. RESULTS: Thirteen patients were eligible [8 creatine transporter deficiency (CTD), 4 guanidinoacetate methyltransferase (GAMT) deficiency and 1 could not be classified]. The mean (±SD) age at diagnosis was 7.2(±5.0) years. Clinical manifestations included intellectual disability (ID) with significant expressive speech delay in all. Most had significant behavior issues (8/13) and/or autism (8/13). All had history of convulsive seizures (11/13 had epilepsy; 2 patients only had febrile seizures) and 2/13 had movement disorder. Constipation was the commonest non-neurological manifestation (5/13 patients). Cranial MRI was normal in all CTD patients but showed globus pallidus hyperintensity in all four with GAMT deficiency. MRS performed in 11/13 patients, revealed abnormally low creatine peak. A causative genetic variant (novel mutation in nine) was identified in 12 patients. Three GAMT deficiency and one CTD patient reported neurodevelopmental improvement and good seizure control after creatine supplementation. CONCLUSION: Intellectual disability, disproportionate speech delay, autism, and epilepsy, were common in our CCDS patients. A normal structural neuroimaging with easily controlled febrile and/or afebrile seizures differentiated CTD from GAMT deficiency patients who had abnormal neuroimaging and often difficult to control epilepsy and movement disorder.

3-M syndrome - a primordial short stature disorder with novel CUL7 mutation in two Indian patients.

OBJECTIVE: To evaluate the cause of short stature in children. CASE PRESENTATION: Two children with suspected skeletal dysplasia and short stature were evaluated. CONCLUSIONS: The 3-M syndrome is a primordial growth disorder manifesting severe postnatal growth restriction, skeletal anomalies and prominent fleshy heels. The 3-M syndrome is a genetically heterogeneous disorder and the phenotype is similar. This is a rare autosomal recessive disorder with normal intellect. Two affected children have been identified by whole-exome sequencing. One patient harboured a compound heterozygous variant and the other was a homozygous missense variant. The genetic diagnosis helped in counselling the families and facilitated prenatal diagnosis in one (case 1) family.

Molecular insights into the role of genetic determinants of congenital hypothyroidism.

In our previous studies, we have demonstrated the association of certain variants of the thyroid-stimulating hormone receptor (TSHR), thyroid peroxidase (TPO), and thyroglobulin (TG) genes with congenital hypothyroidism. Herein, we explored the mechanistic basis for this association using different in silico tools. The mRNA 3'-untranslated region (3'-UTR) plays key roles in gene expression at the post-transcriptional level. In TSHR variants (rs2268477, rs7144481, and rs17630128), the binding affinity of microRNAs (miRs) (hsa-miR-154-5p, hsa-miR-376a-2-5p, hsa-miR-3935, hsa-miR-4280, and hsa-miR-6858-3p) to the 3'-UTR is disrupted, affecting post-transcriptional gene regulation. TPO and TG are the two key proteins necessary for the biosynthesis of thyroid hormones in the presence of iodide and H2O2. Reduced stability of these proteins leads to aberrant biosynthesis of thyroid hormones. Compared to the wild-type TPO protein, the p.S398T variant was found to exhibit less stability and significant rearrangements of intra-atomic bonds affecting the stoichiometry and substrate binding (binding energies, ΔG of wild-type vs. mutant: ‒15 vs. ‒13.8 kcal/mol; and dissociation constant, Kd of wild-type vs. mutant: 7.2E-12 vs. 7.0E-11 M). The missense mutations p.G653D and p.R1999W on the TG protein showed altered ΔG (0.24 kcal/mol and 0.79 kcal/mol, respectively). In conclusion, an in silico analysis of TSHR genetic variants in the 3'-UTR showed that they alter the binding affinities of different miRs. The TPO protein structure and mutant protein complex (p.S398T) are less stable, with potentially deleterious effects. A structural and energy analysis showed that TG mutations (p.G653D and p.R1999W) reduce the stability of the TG protein and affect its structure-functional relationship.

Osmosensing by WNK Kinases.

With No Lysine (K) WNK kinases regulate electro-neutral cotransporters that are controlled by osmotic stress and chloride. We showed previously that autophosphorylation of WNK1 is inhibited by chloride, raising the possibility that WNKs are activated by osmotic stress. Here we demonstrate that unphosphorylated WNK isoforms 3 and 1 autophosphorylate in response to osmotic pressure in vitro, applied with the crowding agent polyethylene glycol (PEG)400 or osmolyte ethylene glycol (EG), and that this activation is opposed by chloride. Small angle x-ray scattering of WNK3 in the presence and absence of PEG400, static light scattering in EG, and crystallography of WNK1 were used to understand the mechanism. Osmosensing in WNK3 and WNK1 appears to occur through a conformational equilibrium between an inactive, unphosphorylated, chloride-binding dimer and an autophosphorylation-competent monomer. An improved structure of the inactive kinase domain of WNK1, and a comparison with the structure of a monophosphorylated form of WNK1, suggests that large cavities, greater hydration, and specific bound water may participate in the osmosensing mechanism. Our prior work showed that osmolytes have effects on the structure of phosphorylated WNK1, suggestive of multiple stages of osmotic regulation in WNKs.

WNKs are potassium-sensitive kinases.

With no lysine (K) (WNK) kinases regulate epithelial ion transport in the kidney to maintain homeostasis of electrolyte concentrations and blood pressure. Chloride ion directly binds WNK kinases to inhibit autophosphorylation and activation. Changes in extracellular potassium are thought to regulate WNKs through changes in intracellular chloride. Prior studies demonstrate that in some distal nephron epithelial cells, intracellular potassium changes with chronic low- or high-potassium diet. We, therefore, investigated whether potassium regulates WNK activity independent of chloride. We found decreased activity of Drosophila WNK and mammalian WNK3 and WNK4 in fly Malpighian (renal) tubules bathed in high extracellular potassium, even when intracellular chloride was kept constant at either ∼13 mM or 26 mM. High extracellular potassium also inhibited chloride-insensitive mutants of WNK3 and WNK4. High extracellular rubidium was also inhibitory and increased tubule rubidium. The Na+/K+-ATPase inhibitor, ouabain, which is expected to lower intracellular potassium, increased tubule Drosophila WNK activity. In vitro, potassium increased the melting temperature of Drosophila WNK, WNK1, and WNK3 kinase domains, indicating ion binding to the kinase. Potassium inhibited in vitro autophosphorylation of Drosophila WNK and WNK3, and also inhibited WNK3 and WNK4 phosphorylation of their substrate, Ste20-related proline/alanine-rich kinase (SPAK). The greatest sensitivity of WNK4 to potassium occurred in the range of 80-180 mM, encompassing physiological intracellular potassium concentrations. Together, these data indicate chloride-independent potassium inhibition of Drosophila and mammalian WNK kinases through direct effects of potassium ion on the kinase.

Newborn screening and single nucleotide variation profiling of TSHR, TPO, TG and DUOX2 candidate genes for congenital hypothyroidism.

High prevalence of congenital hypothyroidism (CH) among Indian newborns prompted us to establish population-specific reference ranges of TSH and to explore the contribution of the common genetic variants in TSHR, TPO, TG and DUOX2 genes towards CH. A total of 1144 newborns (593 males and 551 females) were screened for CH. SNV profiling (n = 22) spanning three candidate genes, i.e. TSHR, TPO and TG was carried out in confirmed CH cases (n = 45). In screen negative cases (n = 700), ten TSHR variants were explored to establish association with CH. No mutation found in DUOX2. The 2.5th to 97.5th percentiles of TSH in these newborns were 0.5 to 12.2 mU/L. In newborns with optimal birth weight, the cut-off was 10 mU/L. Lower or higher birth weight resulted in slightly higher TSH. Two TSHR variants, i.e. rs7144481 and rs17630128 were associated with agenesis, hypoplasia and goiter. The rs2268477 was associated with agenesis and hypoplasia. The rs1991517, rs2075176 and rs2241119 were associated with agenesis only. The rs7144481, rs17630128, rs1991517 and rs2268477 were associated with 2.17, 4.62, 2.91 and 2.29-fold increased risk for CH, respectively. Among the TPO variants, rs867983 and rs2175977 were associated with agenesis and goiter, respectively. Among the TG variants, rs2076740 showed association with agenesis and goiter. Two rare variants i.e. TPO g.IVS14-19 G>C and TG c.1262 C>T were observed in CH cases. No genetic variant identified in the two exons of DUOX2. To conclude, the current study established Indian population-specific normative values for TSH and demonstrates specific genotype-phenotype correlations among three candidate genes.

The rs1991517 polymorphism is a genetic risk factor for congenital hypothyroidism.

The objective of the current study is to explore the association of thyroid-stimulating hormone receptor (TSHR) rs1991517 polymorphism (c.2337 C > G, p.D727E) with congenital hypothyroidism (CH) through a case-control study followed by a meta-analysis. The case-control study was based on 45 CH subjects and 700 healthy controls. Meta-analysis comprised of seven published studies and our current findings (1044 CH cases and 1649 healthy controls). The allele contrast model showed that the presence of G- allele increased CH risk by 45% (OR: 1.45, 95% CI 1.20-1.76) and 41% (OR: 1.41, 95% CI 1.03-1.93) in fixed effect and random effect models, respectively. The GG- genotype is associated with 2.3-fold (95% CI 1.32-3.99) increased risk for CH in the fixed-effect model. I 2 (0.58) and Cochran's Q test (Q: 16.72, p = 0.02) revealed evidence of heterogeneity in the association. No publication bias was observed by Egger's test (p = 0.70). Sensitivity analysis revealed that even after excluding any study this polymorphism is associated with risk for CH. The rs1991517 mutation alters the binding affinity to cAMP (ΔG of 727D vs.727E: - 7.27 vs. - 7.34 kcal/mol). In conclusion, rs1991517 is a genetic risk factor for CH and exerts its impact by altering cAMP-mediated signal transduction.

Mitochondrial carbonic anhydrase VA deficiency in three Indian infants manifesting early metabolic crisis.

BACKGROUND: Hyperammonemia and hyperlactatemia in neonates and young children with non-specific biochemical markers poses a diagnostic challenge. An accurate diagnosis is essential for effective management. CASE REPORTS: We present three infants from unrelated families, one with infantile and two with neonatal hyperammonemic encephalopathy, hypoglycaemia, and hyperlactatemia. The underlying cause was confirmed following whole exome sequencing as biochemical markers were not conclusive of a definite diagnosis. CONCLUSION: The combination of hyperammonemic encephalopathy, hyperlactatemia and hypoglycemia in neonates and infants should prompt physicians to suspect Carbonic anhydrase VA deficiency. Majority of these children can have a favourable long-term outcome with symptomatic treatment.

A Phosphorylated Intermediate in the Activation of WNK Kinases.

WNK kinases autoactivate by autophosphorylation. Crystallography of the kinase domain of WNK1 phosphorylated on the primary activating site (pWNK1) in the presence of AMP-PNP reveals a well-ordered but inactive configuration. This new pWNK1 structure features specific and unique interactions of the phosphoserine, less hydration, and smaller cavities compared with those of unphosphorylated WNK1 (uWNK1). Because WNKs are activated by osmotic stress in cells, we addressed whether the structure was influenced directly by osmotic pressure. pWNK1 crystals formed in PEG3350 were soaked in the osmolyte sucrose. Suc-WNK1 crystals maintained X-ray diffraction, but the lattice constants and pWNK1 structure changed. Differences were found in the activation loop and helix C, common switch loci in kinase activation. On the basis of these structural changes, we tested for effects on in vitro activity of two WNKs, pWNK1 and pWNK3. The osmolyte PEG400 enhanced ATPase activity. Our data suggest multistage activation of WNKs.

In silico analysis of the structural and functional implications of SLC19A1 R27H polymorphism.

In view of the documented association of solute carrier family 19 member 1 (SLC19A1) G80A (R27H) polymorphism with the risk for different types of cancers and systemic lupus erythematosus (SLE), we have reanalysed the case-control study on breast cancer to ascertain the conditions in which this polymorphic variant exerts the risk of breast cancer. Association statistics have revealed that this polymorphism exerts the risk for breast cancer under the conditions of low folate intake, and in the absence of well-documented protective polymorphism in cytosolic serine hydroxymethyltransferase. To substantiate this observation, we have developed a homology model of SLC19A1 using glycerol-3-phosphate transporter (d1pw4a) as a template where 73% of the residues were modelled at 90% confidence while 162 residues were modelled ab initio. The wild and mutant proteins shared same topology in S3, S5, S6, S7, S11 and S12 transmembrane domains. The topology varied at S1 (28-43 residue vs 28-44 residue), S2 (66-87 residue vs 69-87 residue), S4 (117-140 residue vs 117-139 residue), S8 (305-325 residue vs 305-324 residue), S9 (336-356 residue vs 336-355residue), and S10 (361-386 residue vs 361-385 residue) transmembrane domains between wild versus mutant proteins. S2 domain is shortened by three amino acid residues in the mutant while in other domains the difference corresponds to one amino acid residue. The 3DLigandSite analysis revealed that the metallic-ligand-binding sites at 273Trp, 277Asn, 379Leu, 439Phe and 442Leu are although unaffected, there is a loss of active sites corresponding to nonmetallic ligand binding. Tetrahydrofolate and methotrexate have lesser affinity towards the mutant protein than the wild protein. To conclude, the R27H polymorphism affects the secondary and tertiary structures of SLC19A1 with the significant loss in ligand-binding sites.

Intracellular Chloride and Scaffold Protein Mo25 Cooperatively Regulate Transepithelial Ion Transport through WNK Signaling in the Malpighian Tubule.

Background With No Lysine kinase (WNK) signaling regulates mammalian renal epithelial ion transport to maintain electrolyte and BP homeostasis. Our previous studies showed a conserved role for WNK in the regulation of transepithelial ion transport in the Drosophila Malpighian tubule.Methods Using in vitro assays and transgenic Drosophila lines, we examined two potential WNK regulators, chloride ion and the scaffold protein mouse protein 25 (Mo25), in the stimulation of transepithelial ion flux.ResultsIn vitro, autophosphorylation of purified Drosophila WNK decreased as chloride concentration increased. In conditions in which tubule intracellular chloride concentration decreased from 30 to 15 mM as measured using a transgenic sensor, Drosophila WNK activity acutely increased. Drosophila WNK activity in tubules also increased or decreased when bath potassium concentration decreased or increased, respectively. However, a mutation that reduces chloride sensitivity of Drosophila WNK failed to alter transepithelial ion transport in 30 mM chloride. We, therefore, examined a role for Mo25. In in vitro kinase assays, Drosophila Mo25 enhanced the activity of the Drosophila WNK downstream kinase Fray, the fly homolog of mammalian Ste20-related proline/alanine-rich kinase (SPAK), and oxidative stress-responsive 1 protein (OSR1). Knockdown of Drosophila Mo25 in the Malpighian tubule decreased transepithelial ion flux under stimulated but not basal conditions. Finally, whereas overexpression of wild-type Drosophila WNK, with or without Drosophila Mo25, did not affect transepithelial ion transport, Drosophila Mo25 overexpressed with chloride-insensitive Drosophila WNK increased ion flux.Conclusions Cooperative interactions between chloride and Mo25 regulate WNK signaling in a transporting renal epithelium.

Mutational Analysis of Androgen Receptor Gene in Two Families with Androgen Insensitivity.

BACKGROUND: Androgen insensitivity syndrome (AIS) is a rare X-linked disorder due to mutations in the androgen receptor (AR) gene causing end-organ resistance to the androgenic hormone. SUBJECTS AND METHODS: Genetic studies were carried out in two families by karyotype and targeted exome sequencing of the AR gene. RESULTS: Two novel missense mutations were identified, p.L822P and p.P392S, in two families with complete androgen insensitivity (CAIS) and partial androgen insensitivity (PAIS), respectively. Both had 46, XY karyotype. The mother was a heterozygous carrier in PAIS and negative in CAIS. These two were novel mutations, reported for the first time, in the AR gene. In silico analysis predicted that both mutations were damaging. We reviewed the various reported Indian mutations in the AR gene. CONCLUSION: AR gene mutations cause a wide spectrum of disorders from CAIS to male infertility or primary amenorrhea. Early diagnosis is essential for gender assignment and further management, family counseling, and prenatal diagnosis.

Targeted exome sequencing for the identification of complementation groups in methylmalonic aciduria: A south Indian experience.

OBJECTIVES: In view of high incidence of methylmalonic aciduria (MMA) among South Indians, we have performed clinical, biochemical and molecular genetic evaluation of fifteen patients. DESIGN AND METHODS: Targeted exome sequencing was performed for a panel of MMA causing genes i.e. MUT, ABCD4, ACSF3, CD320, LMBRD1, MCEE, MMAA, MMAB, MMACHC, MMADHC. RESULTS: Methylmalonyl-CoA mutase (MUT), MMAB and MMAA genetic variants were found to contribute towards 40%, 33.3% and 6.6% etiology, respectively. Early onset of the disease (during the neonatal period) and presence of MUT and MMAB genetic variants was shown to be associated with higher mortality. The patients with MMAA variants had a milder disease. Among the identified mutations, 66% were already known. Three novel mutations, i.e.MUTp.Ala376Serfs, MMAB p.Glu112* and MMAA p.Tyr24* were identified. We also report three novel variants with predicted pathogenicity, MMAA intron 3 c.562+1_562+2insT, p.Ala668Pro in exon 12 of one of the alleles of the MUT gene and c.519+1G>A in intron 6 of one of the alleles in MMAB gene. We performed prenatal diagnosis in five of these families. CONCLUSIONS: MMA among South Indian patients is genetically heterogeneous, caused by different complementation groups. Both B12-responsive and non-responsive patients were diagnosed. In biochemically diagnosed patients, targeted exome sequencing is cost effective to identify different MMA causing mutations and facilitate genetic counseling.

Clinical and Mutational Characterizations of Ten Indian Patients with Beta-Ketothiolase Deficiency.

Beta-ketothiolase deficiency (mitochondrial acetoacetyl-CoA thiolase (T2) deficiency) is an inherited disease of isoleucine catabolism and ketone body utilization caused by ACAT1 mutations. We identified ten Indian patients who manifested with ketoacidotic episodes of variable severity. The patients showed increased urinary excretion of isoleucine-catabolic intermediates: 2-methyl-3-hydroxybutyrate, 2-methylacetoacetate, and tiglylglycine. Six patients had a favorable outcome, one died, and three developed neurodevelopmental sequela. Mutational analysis revealed a common (p.Met193Arg) and four novel (p.Ile323Thr, p.Ala215Asn, c.1012_1015dup, and c.730+1G>A) ACAT1 mutations. Transient expression analyses of wild-type and mutant cDNA were performed at 30, 37, and 40°C. A p.Ile323Thr mutant T2 was detected with relative enzyme activity and protein amount of 20% and 25%, respectively, compared with wild type at 37°C; it was more prevalent at 30°C but ablated at 40°C. These findings showed that p.Ile323Thr had a significant residual T2 activity with temperature-sensitive instability. Neither residual enzymatic activity nor mutant T2 protein was identified in p.Met193Arg, p.Ala215Asn, and c.1012_1015dup mutations using supernatants; however, these mutant T2 proteins were detected in insoluble pellets by immunoblot analysis. Expression analyses confirmed pathogenicity of these mutations. T2 deficiency has a likely high incidence in India and p.Met193Arg may be a common mutation in the Indian population.

Neuro-fuzzy model of homocysteine metabolism.

In view of well-documented association of hyperhomocysteinaemia with a wide spectrum of diseases and higher incidence of vitamin deficiencies in Indians, we proposed a mathematical model to forecast the role of demographic and genetic variables in influencing homocysteinemetabolism and investigated the influence of life style modulations in controlling homocysteine levels. Total plasma homocysteine levels were measured in fasting samples using reverse phase HPLC. Multiple linear regression (MLR) and neuro-fuzzy models were developed. The MLR model explained 64% variability in homocysteine, while the neurofuzzy model showed higher accuracy in predicting homocysteine with a mean absolute error of 0.00002 µmol/L. Methylene tetrahydrofolate reductase (MTHFR) C677T, 5-methyltetrahydrofolate homocysteine methyltransferase (MTR) A2756G and 5- methyltetrahydrofolate homocysteine methyltransferase reductase (MTRR) A66G were shown to be positively associatiated with homocysteine, while nonvegetarian diet, serine hydroxymethyltransferase 1 (SHMT1) C1420T and TYMS 5'-UTR 28 bp tandem repeat exhibited negative association with homocysteine. The protective role of SHMT1 C1420T was attributed to more H-bonding interactions in the mutant modelled compared to the wild type, as shown through in silico analysis. To conclude, polymorphisms in genes regulating remethylation of homocysteine strongly influence homocysteine levels. The restoration of one-carbon homeostasis by SHMT1 C1420T or increased flux of folate towards remethylation due to TYMS 5'-UTR 28 bp tandem repeat or nonvegetarian diet can lower homocysteine levels.