Therapeutic and Toxic Effects of Valproic Acid Metabolites

Abstract

Valproic acid (VPA) and its salts are psychotropic drugs that are widely used in neurological diseases (epilepsy, neuropathic pain, migraine, etc.) and psychiatric disorders (schizophrenia, bipolar affective disorder, addiction diseases, etc.). In addition, the indications for the appointment of valproate have been expanding in recent years in connection with the study of new mechanisms of action of therapeutic and toxic metabolites of VPA in the human body. Thus, VPA is considered a component of disease-modifying therapy for multiple tumors, neurodegenerative diseases (Huntington’s disease, Parkinson’s disease, Duchenne progressive dystrophy, etc.), and human immunodeficiency syndrome. The metabolism of VPA is complex and continues to be studied. Known pathways of VPA metabolism include: β-oxidation in the tricarboxylic acid cycle (acetylation); oxidation with the participation of cytochrome P-450 isoenzymes (P-oxidation); and glucuronidation. The complex metabolism of VPA explains the diversity of its active and inactive metabolites, which have therapeutic, neutral, or toxic effects. It is known that some active metabolites of VPA may have a stronger clinical effect than VPA itself. These reasons explain the relevance of this narrative review, which summarizes the results of studies of blood (serum, plasma) and urinary metabolites of VPA from the standpoint of the pharmacogenomics and pharmacometabolomics. In addition, a new personalized approach to assessing the cumulative risk of developing VPA-induced adverse reactions is presented and ways for their correction are proposed depending on the patient’s pharmacogenetic profile and the level of therapeutic and toxic VPA metabolites in the human body fluids (blood, urine).

1. Introduction

Valproic acid (CH3CH2CH2)2CHCOOH 2-propylvaleric acid, VPA) is a fatty acid derivative originally synthesized in 1881 by Beverly S. Burton. For almost a century, it has been used as a popular organic solvent in industry and pharmaceuticals [1] (Figure 1).

In 1963, George Karraz made an unexpected discovery while researching the anticonvulsant effects of hellin, finding that all samples dissolved in VPA had the same degree of antiepileptic activity. VPA (Depakene) was first approved by the United States Food and Drug Administration in 1978 [2,3]. In subsequent years, valproates began to be used for many neurological diseases and mental disorders, including: epilepsy [2]; bipolar affective disorder [4,5]; migraine [6]; neuropathic pain [7]; schizophrenia [8]; addiction to alcohol consumption [9], as well as to reduce addictive behavior and detoxification after taking narcotic substances (for example, methamphetamine) [10,11]; compulsive sexual behavior [12]; human immunodeficiency virus infection [13,14]; and progressive Duchenne muscular dystrophy [15]. Moreover, VPA has additional indications as a component of disease-modifying therapy for massive tumors due to its oncostatic effect [16].

As a rule, for most of the above-mentioned disorders, the use of VPA is chronic. This requires maintaining an optimal balance between its effectiveness and safety. This balance depends on the genetically determined features of VPA metabolism and transport in the human body [17], the dose, and duration of VPA intake. From the standpoint of metabolomics, VPA metabolites are of great interest, since they can have both the expected therapeutic effect and predicted or unexpected toxic effects [18]. Thus, toxic metabolites of VPA can lead to an increase in the severity of neurological diseases and mental disorders [19], congenital malformations in the fetus [20], dysfunction of the neuroendocrine system [21], and impaired hematopoiesis [22]. The hepatotoxic effect of VPA metabolites has been well studied [23,24].

Some of the toxic metabolites of VPA lead to disturbances in replication and deoxyribonucleic acid (DNA) synthesis [25]. However, these same metabolites can be useful in the treatment of multiple tumors [26,27,28].

The process of formation of toxic metabolic products is called “toxification”, and biotransformation products with high toxicity are referred to as toxic metabolites [29]. In some cases, the toxic VPA metabolite is an unstable compound that undergoes further transformations in the human body. In this case, it is also called an intermediate or reactive (also known as active) metabolite. It is reactive VPA metabolites that are the substances that cause damage to biosystems at the molecular level [30].

A common property of almost all reactive metabolites is their electrically deficient state, i.e., high electrophilicity [31]. Reactive toxic metabolites of VPA can interact with electron rich (nucleophilic) molecules, damaging them [32]. These nucleophilic molecules include cell macromolecules, the structure of which includes large number of oxygen, nitrogen, and sulfur atoms. These are, first of all, proteins, DNA, and ribonucleic acid (RNA) [33,34]. Reactive VPA metabolites either attach to nucleophilic molecules, forming covalent bonds with them, or cause their oxidation. In both cases, the structure of macromolecules is disturbed, and therefore their functions are also impaired [35]. Bioactivation is not always accompanied by damage to the bio-substrate, since the processes of detoxification and reparation proceed simultaneously in the human body [36]. The intensity of these two processes may be sufficient to compensate for the damage associated with the formation of reactive toxic metabolites of VPA. However, when using high doses of VPA and/or chronic pharmacotherapy, protective mechanisms may fail, leading to the development of serious VPA-induced adverse drug reactions (ADRs) [37].

The main organ of VPA metabolism in the human body is the liver, mainly due to the diversity and high activity of biotransformation enzymes there [38]. Despite the dominant role of the liver in VPA metabolism, other organs also take part in this process. The kidneys and lungs contain enzymes of I and II phases of VPA metabolism [39,40]. In addition, some of these enzymes are also expressed in the brain [41]. The role of the kidneys is particularly great because this organ has a specific system for the capture and catabolism of VPA conjugation products formed in the liver [42]. The activity of other organs, such as the brain, blood, intestine, spleen, muscle tissue, and placenta, is much lower [43,44,45,46,47]. However, genetically determined expression level and functional activity of enzymes catalyzing VPA biotransformation processes in these human organs and tissues will play a key role in the development of VPA-induced ADRs in these organs and tissues with a systemic and local elevation of reactive toxic VPA metabolites. In the process of extrahepatic metabolism of VPA, compounds (molecules) can be formed, both similar to products of hepatic origin, and different from them. In some cases, extrahepatic metabolism may activate VPA metabolites formed in the liver.

The metabolism of VPA is complicated and continues to be studied. Known pathways of VPA metabolism include: β-oxidation in the tricarboxylic acid cycle (acetylation); oxidation with the participation of cytochrome P-450 isoenzymes (P-oxidation); and glucuronidation [19,48,49]. The complex metabolism of VPA explains the diversity of its active and inactive metabolites, which have therapeutic, neutral, or toxic effects. It is known that some active metabolites of VPA may have a stronger clinical effect than VPA itself.

These reasons explain the relevance of this narrative review, which summarizes the results of studies of the blood (serum, plasma) and urinary metabolites of VPA from the standpoint of pharmacogenetics and pharmacometabolomics. We have searched full-text publications in PubMed, Web of Science, Springer, Google Scholar, and e-Library databases. Keywords were as follows: valproic acid; blood metabolite; urinal metabolite; oxidation; glucuronidation; acetylation; pharmacometabolomics; pharmacogenomics; adverse drug reaction; risk factor; and personalized approach. In addition, earlier publications of historical interest are included in the narrative review.

2. Mechanism of Valproic Acid Action

The exact mechanisms of action of VPA (Figure 2) and its metabolites in neurological disease and psychiatric disorders are unknown, but there are several ways to explain this. VPA and its active therapeutic metabolites are known to inhibit succinic semialdehyde dehydrogenase [49]. This inhibition leads to an increase in succinic semialdehyde, which acts as a gamma-aminobutyric acid (GABA) transaminase inhibitor, decreasing GABA metabolism and increasing GABAergic neurotransmission. Because GABA is an inhibitory neurotransmitter, this increase results in increased inhibitory activity in the central nervous system. A possible secondary factor in the inhibition of neurons in the cerebral cortex is the direct suppression of the activity of voltage-dependent sodium channels and indirect suppression of the activity of cortical neurons through the effect on GABAergic neurotransmission [50,51]. VPA and its active (therapeutic) metabolites can affect the extracellular signal-related kinase (ERK) pathway [52]. These effects are dependent on mitogen-activated protein kinase (MEK) and result in ERK1/2 phosphorylation. This MEK-induced activation increases the expression of several downstream targets, including ELK-1, followed by an increase in c-fos, a growth cone-associated protein-43 that promotes neuronal plasticity; B-cell lymphoma/leukemia-2, which is an anti-apoptotic protein; and brain-derived neurotrophic factor, which is also involved in neuronal plasticity and growth [2]. The increase in neurogenesis and neurite growth under the influence of VPA and its active metabolites is explained by the effects of this pathway. An additional effect of the VPA-induced increase in brain-derived neurotrophic factor expression is an increase in the number of GABA receptors, which contributes to an increase in inhibitory GABAergic activity [53].

VPA has a non-competitive indirect inhibitory effect on myo-inositol-1-phosphate synthetase. This leads to a decrease in de novo synthesis of inositol monophosphatase and subsequent depletion of inositol [54]. VPA and its active metabolites influence also leads to the downregulation of protein kinase C (PKC)–α and -ε proteins, which are potentially associated with bipolar disorder, since PKC is not regulated in the frontal cortex in patients with bipolar disorder [55]. Inhibition of the PKC pathway may also contribute to the prevention of migraine attacks [56]. The level of the substrate PKC is also reduced by VPA and may contribute to changes in synapse remodeling by affecting the cytoskeleton [57].

VPA and its metabolites are also known to affect fatty acid metabolism. It is believed that less incorporation of fatty acid substrates into sterols and glycerol-lipids affects the fluidity of the neuronal cell membrane and leads to an increase in the action potential threshold, which may contribute to the antiepileptic effect of valproates [58]. It has been shown that VPA is a non-competitive direct inhibitor of microsomal long-chain fatty acyl-CoA synthetase in the brain. Inhibition of this enzyme reduces the availability of arachidonyl-CoA, a substrate for the production of inflammatory prostaglandins. It is thought that this may be one of the mechanisms behind VPA’s effectiveness in migraine prevention, as migraines are usually treated with non-steroidal anti-inflammatory drugs, which also suppress prostaglandin production [59].

Finally, VPA and its active metabolites act as direct inhibitors of histone deacetylase (HDAC) [51]. Hyperacetylation of lysine residues on histones promotes the unwinding of DNA and allows for an increase in gene transcription. Such genomic effects have a wide spectrum: the activity of the 461 gene is increased or suppressed under the influence of VPA and its active metabolites (both therapeutic and toxic). The relationship of these genomic effects to the therapeutic value of VPA is not fully understood; however, hyperacetylation of H3 and H4 correlates with a decrease in the severity of bipolar affective disorder [60]. Histone hyperacetylation in the BDNF (Brain-Derived Neurotrophic Factor) gene, which increases brain-derived neurotrophic factor expression, occurs after an epileptic seizure and is thought to be a neuroprotective mechanism that VPA may enhance or prolong. Hyperacetylation of H3 is associated with decreased activity of glyceraldehyde 3-phosphate dehydrogenase, a pro-apoptotic enzyme, further contributing to the neuroprotective effects of VPA [61].

3. Pharmacometabolomics and Pharmacogenomics of Valproic Acid

Numerous studies have shown that human cytochrome P450 (CYP) isoenzymes play a crucial role in VPA metabolism [62,63]. A key CYP-mediated branch of the VPA pathway produces the metabolites 4-hydroxy(OH)-VPA and 5-OH-VPA via hepatic cytochrome P450 isoenzymes CYP2C9, CYP2B6, and CYP2A6. The use of a selective inhibitor of the CYP2C9 isoenzyme, sulfafenazole, can dramatically reduce the formation of these three metabolites; however, sulfafenazole is most likely a catalyst for the formation of 4-ene-VPA rather than 4-OH-VPA or 5-OH-VPA [5]. Isoenzymes CYP2A6 and CYP2B6 together provide the formation of 20–25% of 4-ene-VPA, 4-OH-VPA, and 5-OH-VPA. In addition, the CYP2A6 isoenzyme mediates the oxidation of VPA to 3-OH-VPA, and the use of coumarin, a potent inhibitor of CYP2A6, significantly reduces the formation of 3-OH-VPA in human liver microsomes [5].

Analysis of VPA metabolism (1 nM) in vitro showed that the common (wild) single nucleotide variant (SNV) CYP2C9*1 (wild type) of the CYP2C9 gene, encoding a fully functional CYP2C9 isoenzyme, is associated with an active process of 4-hydroxylation and 5-hydroxylation of VPA by 75–80%. In contrast, the CYP2A6 isoenzyme contributes about 50% to the VPA 3-hydroxylation process. The role of CYP2A6 and CYP2B6 isoenzymes in the oxidative metabolism of VPA varies depending on the catalytic ability of these isoenzymes in human hepatic microsomes [64].

The carriage of low-functioning (or non-functional) and high-functioning SNVs in the CYP2A6, CYP2B6, and CYP2C9 genes may partly explain the interindividual variability in VPA metabolism in the human body. Thus, patients with heterozygous or homozygous CYP2A6*4/*4 genotypes (CYP2A6 gene deletion), known as poor metabolizers (PM), had higher plasma concentrations of VPA than extensive metabolizers (EM), homozygous carriers of the fully functional allele. Other SNVs may have a different effect on plasma concentrations of VPA and its metabolites, for example: carriers of the highly functional CYP2B6*4 allele (rs2279343; 785A>G), known as ultra-rapid metabolizers (URM), had higher activity of the CYP2B6 isoenzyme compared to carriers of the extensive variant (wild type), although protein expression was slightly reduced in patients with the URM phenotype. On the contrary, patients with the PM phenotype, who are carriers of non-functional SNVs of the CYP2B6 gene (rs3745274, 516G>T; rs2279343, 785A>G), had a higher concentration of VPA in the blood plasma compared to patients with the EM phenotype [65].

Patients, heterozygous and homozygous carriers of the CYP2C9*3 (rs1057910; 1075A>C) and CYP2C9*2 (rs1799853; 430C>T) alleles, had an increased concentration of toxic VPA metabolites in plasma compared to homozygous carriers of the CYP2C9*1 allele [66]. It has been shown that elderly female patients require a 30–50% lower dose of VPA to achieve therapeutic concentrations of the drug compared to young men [67]. Studying the CYP2C9 profile in pediatric patients before starting treatment allowed us to optimize the dose of VPA and minimize exposure to suboptimal or toxic concentrations, thereby reducing the incidence of ADRs such as hyperammonemia. Moreover, carriers of CYP2C9*3 allelic variant, CYP2A6*1/*4 and CYP2A6*4/*4 genotypes, had elevated hepatotoxic levels of 4-en-VPA and/or 2,4-diene-VPA compared to homozygous carriers of the wild-type allele CYP2A6*1/*1 [68]. Low-functional SNVs in the CYP2C9 and CYP2A6 genes are risk factors for an increase in the plasma level of toxic VPA metabolites and an increase in the risk of hepatotoxicity up to 7.50 and 5.13 times, respectively. Therefore, the study of low-functioning SNVs in the CYP2C9 and CYP2A6 genes before treatment can predict or prevent serious VPA-induced ADRs [5].

The individual variability of alleles in the CYP2A6, CYP2B6, and CYP2C9 genes can explain the significant variability in the pharmacokinetics of VPA and the formation of its therapeutic and toxic metabolites in different people. For example, non-functional SNVs of CYP2C9*2 and CYP2C9*3 are associated with decreased VPA metabolism, so patients who are homozygous for CYP2C9*2 or CYP2C9*3 or who are heterozygous (CYP2C9*2/*3) have a PM phenotype and show decreased P-oxidation of VPA in liver microsomes [69]. Patients with homozygous CYP2C9*2/*2 and CYP2C9*3/*3 genotype had reduced oxidative biotransformation of VPA in liver microsomes compared to subjects with compound heterozygous CYP2C9*2/*3 genotype and catalyzed the formation of toxic metabolites 4-ene-VPA, 4-OH-VPA, and 5-OH-VPA [37]. Therefore, pharmacogenetic profiling of the CYP2C9 gene in patients with neurological diseases and psychiatric disorders may help optimize VPA dosage, and prevent ADRs and iatrogenic worsening of the disease.

In addition, CYP2C19*2 (rs4244285; 681G>A) and CYP2C19*3 (rs4986893; 636G>A) are associated with changes in the volume of distribution and plasma concentrations of VPA and its metabolites in patients with epilepsy. Thus, carriers of a highly functional CYP2C19*2 variant require higher doses of VPA to achieve target therapeutic plasma concentrations (450 μg/mL), indicating that the CYP2C19 isoenzyme is also involved in the metabolic pathway of VPA. Carriers of the CYP2C9*13 variant (rs72558187) did not require dose adjustment, as there was no correlation between the CYP2C9*13 variant and plasma VPA concentration [5]. The mean concentration/dose ratios (CDRs) of VPA were significantly higher in patients with the heterozygous CYP2C19*1/*2 genotype or compound heterozygous CYP2C19*2/*3 genotype than in patients with the homozygous CYP2C19*1/*1 genotype. The dose of VPA needs to be reduced for patients with intermediate metabolizer (IM) and PM phenotypes to achieve average blood levels of VPA in patients with the EM phenotype (CYP2C19*1/*1 genotype) [5].

A study in freshly extracted rat hepatocytes showed that VPA-induced oxidative stress and mitochondrial dysfunction preceded VPA hepatotoxicity. In addition, toxic metabolites of VPA induced hepatotoxicity by causing leakage of lysosomal membranes, as well as the formation of reactive oxygen species (ROS) as a result of metabolic activation of the CYP2E1 isoenzyme. It has been shown that the CYP2E1 isoenzyme is effective for the production of ROS and is one of the most powerful inducers of oxidative stress in cells. VPA-induced ROS formation was protected by CYP2E1 inhibitors (1-phenylimidazole and 4-methylpyrazole). In this regard, the role of the SNVs of the CYP2E1 gene encoding this isoenzyme needs to be studied as a potentially important genetic biomarker of metabolic disorders induced by chronic valproate therapy [47].

Although CYPs play a minor role in the VPA metabolic pathway, they are important in assessing the risk of elevated plasma levels of toxic VPA metabolites, especially in patients with impaired uridine glucuronosyltransferase system (UGTs) isoenzyme activity. Due to conflicting results on the effect of CYP isoenzymes SNVs on the pharmacokinetics of VPA, larger cohorts are needed to test these results and explore new candidate genes [17].

VPA metabolism is also affected by the SNVs in the gene encoding medium chain acyl-CoA synthetase, family member 2A (ACSM2A). Carriers of the low-functional and non-functional SNVs of the ACSM2A gene have higher levels of alanine aminotransferase and aspartate aminotransferase compared to carriers of the full-functional allele (wild-type) [70], which is associated with toxic damage to hepatocytes and can lead to an increase in plasma toxic metabolites of VPA.

However, glucuronide conjugation is the predominant pathway for the metabolism and excretion of VPA [48]. Approximately 20–70% of VPA is excreted in the urine as glucuronide conjugates. Research on VPA glucuronide conjugation has focused on several genes of the UGT1A and UGT2B families, including the UGT1A3, UGT1A4, UGT1A8, UGT1A9, UGT1A10, and UGT2B15 genes, as well as the best-studied UGT1A6 and UGT2B7 genes. It has been shown that low-functional and non-functional SNVs in the UGT family genes can affect the dosage and concentration of VPA in the blood plasma and can be directly related to an increase in the blood level of its toxic metabolites and the development of ADRs. A study of recombinant enzymes and human liver microsomes about the effect of three SNVs (19T>G, 541A>G, and 552A>C) in the human UGT1A6 gene demonstrated that, compared to full-functioning SNVs (UGT1A6*1), low-functional SNVs UGT1A6*2 were associated with a 2-fold increase in VPA glucuronidation activity [71].

A study of VPA monotherapy and stable control of epileptic seizures in patients with epilepsy showed that carriers of SNVs UGT1A6*3 (rs6759892; 19T>G), UGT1A6*5 (rs2070959; 541A>G) and UGT1A6*9 (rs1105879; 552A>C) had a reduction in blood levels of VPA therapeutic metabolites. These patients required higher doses of VPA and had a lower concentration-to-dose ratio. Carriers of highly functional alleles of the UGT1A6 gene (19T>G, 552A>C, and 541A>G) demonstrated increased activity of the UGT1A6 enzyme compared to homozygous carriers of the wild-type allele, and carriers of the 552A>C SNV in the UGT1A6 gene had a longer half-life of toxic metabolites of VPA and a lower clearance rate, which often led to VPA-induced ADRs such as ataxia, liver damage, metabolic syndrome, tremors, hallucinations, pancreatitis, and weight gain. However, the effect of these SNVs on VPA metabolism and the development of ADRs is difficult to determine, as alterations in VPA metabolism have been described either when VPA is given alone or in combination with carbamazepine. For example, in carriers of UGT1A3*5 treated with VPA as monotherapy, a decrease in the concentration of active therapeutic metabolites of VPA in blood plasma has been described. Therefore, homozygous and heterozygous carriers of the UGT1A3*5 allele required a higher dose of VPA to reach the therapeutic range in plasma (as measured by therapeutic drug monitoring) in the range of 50–100 μg/mL [72].

The role of the UGT2B7 gene in VPA metabolism has also been extensively studied. Several SNVs in this gene play an important role in the clearance of toxic VPA metabolites. However, studies on the association of a few SNVs, including UGT2B7*2 (802C > T; rs7439366), UGT2B7*3 (211G > T; rs12233719) or UGT2B7*4 (1192G > A; rs145725367), with plasma levels of therapeutic and toxic VPA metabolites are controversial. However, carriers of the homozygous and heterozygous genotypes (TT and CT, respectively) UGT2B7*2 had lower plasma concentrations of VPA compared to carriers of the homozygous CC genotype, suggesting that carriers of the TT and CT genotypes may require higher doses of VPA [5]. Furthermore, patient age has been shown to be positively correlated with adjusted serum VPA concentrations in patients with pediatric epilepsy. Thus, carriers of UGT1A3*5 require a higher dose of VPA to achieve the reference (therapeutic) corridor in the plasma in the range of 50–100 μg/mL [73]. In addition, SNVs in the UGT1A3 gene affect plasma levels of VPA therapeutic metabolites [74].

Drug transporter proteins play a critical role in the pharmacokinetics of VPA. Overexpression of efflux drug transporter proteins can also be regulated by the pregnane X receptor (PXR) [75]. Among them, the most studied VPA transporter is P-glycoprotein (P-gp) or multidrug resistance protein. P-gp is a volatile efflux pump that eliminates the action of several active metabolites’ VPA and is a product of the ATP-binding cassette subfamily b member 1 (ABCB1) gene, also known as the multidrug resistance gene 1 (MDR1) [76]. The ABCB1 gene is one of the most important genetic biomarkers affecting drug transport across the brain–blood barrier. Like ABCB1, ABCC2 gene expression is higher in brain endothelial cells of patients with VPA resistance [77]. It is known that SNVs in the ABCB1 gene can directly affect the uptake of VPA by the brain and also the excretion (efflux) of VPA from the brain into the blood. Homozygous TT genotypes of G2677 and C3435 polymorphisms, as well as TT, CTT, and TTT haplotypes of the ABCB1 gene, are significantly associated with the development of therapeutic resistance to VPA. At the same time, these more often had the ABCB1 3435 CC and CT genotypes. Such patients had reduced plasma levels of VPA therapeutic metabolites compared to the homozygous TT genotype [78,79]. In patients with temporal lobe epilepsy, an association was found between therapeutic resistance to VPA and the CGC haplotype, as well as with the CC or CT genotypes of two SNVs (rs1045642 and rs1128503) of the ABCB1 gene [80]. The results of a meta-analysis [81] demonstrated that the SNV G2677T>A of the ABCB1 gene is a biomarker for an increased risk of developing therapeutic resistance to valproate.

Five key enzymes are involved in the urea cycle, including carbamoyl phosphate synthetase 1 (SPS1), ornithine transcarbamoylase, argininosuccinate synthase, argininosuccinate lyase, and arginase 1. Another N-acetylglutamate synthase (NAGS) expressed in mitochondria is also important for the functioning of the urea cycle as it provides the necessary mutant N-acetylglutamic acid activator. A study in rat liver mitochondria showed that the valproil-CoA metabolite inhibits NAGS activity, leading to ammonia accumulation. Thus, CPS1 is the first limiting enzyme in the urea cycle, which accelerates the conversion of ammonium to carbamoyl phosphate in the liver. Dysregulation of CPS1 function may result from genetic changes or epigenetic regulation in hepatocellular carcinoma. Therefore, it can be assumed that the absence of CPS1 or a decrease in its activity may be largely associated with VPA-induced hyperammonemia. The 4217C > A polymorphism of the CPS1 gene leads to the amino acid conversion of threonine to asparagine, which correlates with the low functional activity of the CPS1 enzyme [82]. This SNV has been shown to be a risk factor for VPA-induced hyperammonemia [24,83] and the development of hyperammoniemic encephalopathy.

In addition, significant impairments in mitochondrial oxidation of toxic VPA metabolites are observed in patients with SNVs of the DNA polymerase gamma, catalytic subunit (POLG) gene encoding polymerase γ. The POLG enzyme is defined as a mitochondrial DNA polymerase that is associated with various metabolic disorders and monogenic hereditary diseases such as Alpers–Huttenlocher syndrome, which is associated with an increased risk of fatal hepatotoxicity of VPA metabolites [84,85]. Approximately one third of patients with Alpers–Huttenlocher syndrome develop liver failure within 3 months of starting VPA medications. Low-functional SNVs in the POLG gene occur in up to 0.5% of the population [86]. In addition, chronic administration of VPA drugs leads to significant overexpression of the POLG gene and an increase in mitochondrial biogenesis by altering the expression of several mitochondrial genes. Heterozygous mutations p.Q1236H and p.E1143G in the POLG1 gene have been associated with VPA-induced liver failure. Therefore, pharmacogenetic testing (PGx) of SNVs and mutations in the POLG gene is recommended in all patients with suspected mitochondrial disease prior to initiation of VPA treatment. In addition, PGx is useful in all patients with VPA-induced ADRs to rule out their association with the carriage of low-functioning (or non-functional) SNVs in the POLG1 gene. This approach will minimize the risk of liver failure in patients receiving valproate [85,87].

A significant increase in the level of γ-glutamyltransferase (GGT) was shown in patients with heterozygous (GSTM1/GSTM1^-) and homozygous (GSTM1^-/GSTT1^-) genotypes in patients treated with VPA. The increase in GGT enzyme activity in VPA-treated patients may be due to hepatic glutathione (GSH) elimination and has also been presented as an early biomarker of oxidative stress. However, the clinical relationship between the level of GGT and the SNVs in the GST gene and the formation of toxic metabolites of VPA remains unknown. At the same time, an association was found between the homozygous Val/Val genotype in the SOD2 (Superoxide Dismutase 2) gene with an increase in the level of alanine aminotransferase (OR = 3,5; p-value = 0.056) caused by hepatotoxic metabolites of VPA [88].

In addition to hepatotoxicity, VPA metabolites also cause teratogenicity by downregulating IGF2R (Insulin-like growth factor 2 receptor), RGS4 (Regulator of G protein signaling 4), COL6A3 (Collagen, type VI, alpha 3), EDNRB (Endothelin B receptor), and KLF6 (Krueppel-like factor 6) genes, as demonstrated in VPA-induced neural tube defects in a chick embryo prevalence model [89]. The 116C/G variant (rs226957) in the promoter region of the XBP1 (X-Box Binding Protein 1) gene is associated with increased response to active VPA metabolites; in patients carrying the minor G allele, the expected therapeutic response to VPA increased [90].

It has also been shown that SNV rs1137101 of the LEPR (Leptin Receptor) gene, rs1800497 of the ANKK1 (Ankyrin Repeat and Kinase Domain Containing 1) gene, and rs10789038 of the PRKAA2 (Protein Kinase AMP-Activated Catalytic Subunit Alpha 2) gene are associated with VPA-induced weight gain. The administration of valproate to obese patients results in an increase in serum leptin levels compared to non-obese patients [91], but there is no answer yet as to which VPA metabolites induce this ADR. However, it is known that VPA and its active toxic metabolites can reduce the expression of the up-regulated gene 4/up-regulator of cell proliferation URG4/URGCP (URG4/URGCP) and CCND1 (Cyclin D1) genes and suppress the proliferation of neuroblastoma cells is a human-derived line cell used in scientific research SHY5Y [27,92,93].

Thus, pharmacogenomics and pharmacometabolomics can be the most important sources for assessing interindividual differences in the pharmacokinetics of VPA and its metabolites. Optimization of valproate dosage based on the results of pharmacogenetic and pharmacometabolic testing plays an important role in psychopharmacotherapy, so the study of genetic biomarkers that affect the pharmacokinetics of VPA can help develop a new personalized strategy for prescribing valproates to patients with neurological diseases and mental disorders by achieving a balance of efficacy and safety [17].

4. Therapeutic Metabolites of Valproic Acid

A metabolite is any substance formed during metabolism (digestion or other chemical processes in the body) [94]. In addition, the term “metabolite” refers to substances that remain (circulate in the blood, accumulate in organs and tissues, excreted in the urine) after the drug is broken down (metabolized) by the body [95]. Research into therapeutic metabolites of VPA has been ongoing for many years. Such studies have intensified over the past decade due to the rapid development of a new direction in personalized medicine—pharmacometabolomics [96]. The number of studied VPA metabolites is close to fifty [97], but their clinical role (Figure 3) in the expected therapeutic response to valproate has not yet been sufficiently studied for many of them.

The most studied therapeutic metabolites of VPA in the blood and urine are presented in Table 1 and Table 2.

4.1. Valproate Acid Glucuronide

Valproate acid glucuronide (VPAG) is a glucuronidation product of VPA. It is known that VPAG is the main metabolite of VPA excreted in the urine (30–50% of the dose). With the age of the patient, the level of this metabolite in the blood does not change, despite the known decrease in liver weight [97].

4.2. 2-N-Propyl-2-Pentenoic Acid

2-n-propyl-2-pentenoic acid is also known as 2,4-diene VPA. This VPA metabolite is formed in the body by glucuronidation. This therapeutic metabolite of VPA exhibits a high antiepileptic potential similar to VPA, but it has a low teratogenic potential [102], which is of undoubted clinical interest.

4.3. Valproyl-CoA

Valproyl-CoA is a therapeutic metabolite of VPA, which is formed by the metabolism of VPA as a result of β-oxidation (acetylation) in the Krebs cycle [100]. Previously, this metabolite was found in liver cells (hepatocytes), but recent work has shown that valproyl-CoA can accumulate in brain neurons as a result of normal fatty acid turnover processes. Valproil-CoA is a therapeutic metabolite and increases the antiepileptic activity of VPA by stimulating the activity of Na+, K+-ATPase at low concentrations of adenosine triphosphate in the brain [100].

4.4. 4-Ene-Valproic Acid

It is known that 4-ene-valproic acid (4-ene-VPA) belongs to the class of organic compounds known as methyl-branched fatty acids. These are fatty acids with an acyl chain having a methyl branch. Usually, they are saturated and contain only one or several methyl groups [101]. However, other branches than the methyl branch may also be present. In addition, 4-ene-VPA is a very hydrophobic molecule; it is practically insoluble (in water) but refers to relatively neutral VPA metabolites, although this opinion is debatable. Several studies have demonstrated that 4-ene-VPA has an antiepileptic effect and also inhibits HDACs, which is a promising direction in the treatment of cancer, autoimmune disease, and acquired immunodeficiency syndrome [97]. 4-ene-VPA circulates in the blood and is excreted in the urine, which facilitates the study of its level in biological fluids of the human body.

4.5. 2-Ene-Valproic Acid

2-ene-valproic acid (2-ene-VPA) is a product of P-oxidation metabolism in hepatocytes. This therapeutic metabolite of VPA belongs to a class of organic compounds known as methyl-branched fatty acids. In addition, it inhibits HDACs, which is a promising direction in the treatment of cancer, autoimmune disease, and acquired immunodeficiency syndrome [101].

5. Role of Therapeutic Metabolites of Valproic Acid

Thus, VPA is a derivative of fatty acids and HDACs of classes I and IIa [103], which allows VPA and its active metabolites to regulate the expression of various genes [103]. In addition, pathophysiological processes include both excitotoxicity of VPA metabolites and a decrease in gene transcription due to a decrease in the level of histone acetylation, which is accompanied by a loss of GABAergic (gamma-aminobutyric acid) neurons in the striatum as a pathological hallmark of Huntington’s disease [13,104]. VPA is known to exert neuroprotective effects through these primary targets. Excitotoxicity is a mechanism of cell death caused by hyperactivation of excitatory amino acid receptors, which increases the ion permeability of the cell membrane and leads to an overload of intracellular calcium [105]. Notably, VPA and its active metabolites also strongly induce 70 kilodalton heat shock proteins through epigenetic mechanisms, reducing infarct volume and improving functional recovery in rodents subjected to cerebral ischemia [106]. HDAC1 has been implicated in reducing the human immunodeficiency viruses’ activity in infected cells. The inhibitory effect of VPA metabolites on this protein makes it a good candidate for acquired immunodeficiency syndrome therapy [107].

Neuroprotective effects of VPA therapeutic metabolites have been reported in a wide range of acute central nervous system injury models, including post-traumatic inflammation caused by spinal cord injury [108]. The results of the studies demonstrated that VPA and its therapeutic metabolites promote the phenotypic shift of microglia from the M1 phenotype to the M2 phenotype and inhibited microglial activation, thereby reducing the level of pro-inflammatory factors caused by spinal cord injury. The increase in acetylation of signal transducer and activator of transcription 1/ nuclear factor-κB (STAT1/NF-κB pathway) during VPA treatment was probably caused by HDAC3 translocation to the nucleus and its activity. In addition, VPA and its therapeutic metabolites have been shown to downregulate HDAC3 expression and activity and increase acetylation of STAT1 as well as NF-κB p65 after spinal cord injury. The acetylation status of NF-kB p65 and the complex with nuclear factor kappa-light-chain-enhancer of activated B cells NF-κB (NF-κB) p65 and STAT1 suppressed the transcriptional activity of NF-kB p65 and attenuated the central inflammatory response mediated by microglia and led to the neuroprotective effect of VPA therapeutic metabolites [109].

As mentioned earlier, VPA has been described as an HDAC, which makes it possible to consider not only the well-known therapeutic metabolites of VPA, but also the toxic metabolites of VPA, as promising potential therapeutic agents for the treatment of various forms of cancer and leukemia [110]. Chromatin is formed from DNA packaged in nucleosome structures of 146 base pairs of DNA wound around a histone octamer (two copies of each histone: H2A, H2B, H3, and H4) held in place by histone H1. In terms of transcription, the condensed form of chromatin (heterochromatin) is inactive. In contrast, the decondensed form (euchromatin) corresponds to the active form. The transition between euchromatin and heterochromatin depends on two families of proteins: histone acetyltransferases and HDACs. It has been established that histone acetylation leads to relaxation of the nucleosome structure, releasing DNA and allowing activation of transcription. Inhibition of HDAC promotes the formation of decondensed chromatin, thereby promoting gene expression [13]. The toxic metabolites of VPA are among the most promising molecules because, unlike other drugs that target the expression of a molecule or family of molecules, they target chromatin through associated proteins: HDAC; deoxyribonucleic acid methyltransferase; heterochromatin protein 1; and smooth muscle cells. Thus, toxic VPA metabolites affect the expression of many proteins and therefore can be considered as therapeutic and applicable to a wide range of neurodegenerative diseases and autoimmune diseases, as well as multiple tumors and leukemia, where these VPA metabolites repress many antioncogenes during carcinogenesis [13].

Parkinson’s disease is caused by the degeneration of nigrostriatal dopaminergic neurons. No direct effect of VPA on Parkinson’s disease has been shown, but in vitro studies have shown a positive effect of active VPA metabolites in animal models that mimic Parkinson’s disease at various levels. VPA and its active metabolites prevent rotenone-induced apoptosis, protect neurons from the dopaminergic toxin 1-methyl-phenylpyridinium, stimulate the release of neutrophils from neuroglial cells, protect cultured dopaminergic neurons from degeneration caused by excessive microglial activation by stimulating microglial apoptosis, and protect neurons, increasing the expression of α-synuclein and preventing its monoubiquitination and nuclear translocation [13]. Experiments on an animal model of Parkinson’s disease (rodents) have shown that selective changes in α-synuclein induced by rotenone (a decrease in native protein and an increase in monoubiquitination in substantia nigra and striatum) can be reduced with long-term administration of VPA [104]. However, there is still no convincing answer to the question about which of the active metabolites of VPA are the most promising molecules for the treatment of Parkinson’s disease and other neurodegenerative diseases, including Huntington’s disease and Alzheimer’s disease.

The role of active VPA metabolites in the treatment of primary muscular dystrophy (eg., Duchenne muscular dystrophy) has been demonstrated in vitro and in vivo in mdx/utrn^-/- knockout mice. It has been shown that VPA therapeutic metabolites are able to induce the Akt/mTOR/p70S6K pathway through the induction of phosphatidylinositol-3-OH kinase, which leads to a decrease in collagen content and fibrosis in skeletal muscles, a decrease in hind limb contractures, an increase in the integrity of the myocyte sarcolemma, a decrease in the number CD8-positive inflammatory cells, and an increase in the level of activated Akt in the skeletal muscles of experimental animals [13].

Additionally, therapeutic metabolites of VPA have the potential to delay the development of atrial remodeling in patients from the risk group. Therapeutic metabolites of VPA have been shown to attenuate many components of atrial remodeling that are present in transgenic mice, animal models, and humans. However, VPA therapeutic metabolites significantly reduced atrial dilatation, cardiomyocyte enlargement, atrial fibrosis, and myocyte ultrastructure disorganization. In addition, it was demonstrated that these metabolites significantly reduced the occurrence of atrial thrombi and reduced the severity of negative changes in the action potential of cardiomyocytes [111]. An increase in H4 histone acetylation in the atria of VPA-treated transgenic mice confirmed the efficacy of HDAC inhibition in vivo. Cardiomyocyte-specific genetic inactivation of HDAC2 in transgenic mice attenuated ultrastructural myocyte disorganization comparable to VPA. Finally, in transgenic mice, VPA therapeutic metabolites restrain dysregulation of proteins that are involved in multiple pathways associated with atrial fibrillation, such as oxidative phosphorylation or RhoA (Ras homolog gene family, member A) signaling, and disease functions such as cardiac fibrosis and muscle cells apoptosis [111].

6. Toxic Metabolites of Valproic Acid

The toxicity of VPA metabolites refers to the level of damage these compounds can cause to the human body. The toxic effects of VPA metabolites are dose-dependent and may affect the entire system, such as the central nervous system, or a specific organ, such as the liver.

Toxic metabolites of VPA can form in the human body as intermediates, byproducts, or end products of its metabolism. There are many examples of how hepatocytes organize metabolic processes in space and time to prevent the accumulation of these toxic VPA metabolites or reduce their level in organs and tissues in which biological processes that are most sensitive to toxic VPA metabolites occur [13] (Figure 4).

The most studied toxic metabolites in the blood (serum and plasma) and urine are presented in Table 3 and Table 4.

6.1. 3-Hydroxyvalproic Acid

3-Hydroxyvalproic acid (3-OH-VPA) belongs to a class of organic compounds known as hydroxy fatty acids [97]. These fatty acids do have a hydroxyl group in their chain. This VPA metabolite is formed in the body, mainly in hepatocytes, as a result of P-oxidation involving two isoenzymes: CYP3A5 and CYP2A6. This VPA metabolite has a pronounced hepatotoxic effect [97], which explains the importance of PGx (Appendix A) for identifying patients at risk.

6.2. 4-Hydroxyvalproic Acid

4-Hydroxyvalproic acid (4-OH-VPA) also belongs to the class of organic compounds known as hydroxy fatty acids. These are fatty acids that have a hydroxyl group in their chain. This toxic metabolite is the result of P-oxidation of VPA in hepatocytes (primarily) by four isoenzymes: CYP3A5; CYP2A6; CYP2C9; and CYP2B6. It is known that 4-OH-VPA has a pronounced hepatotoxic effect [97], which explains the importance of PGx (Appendix A) for identifying patients at risk.

6.3. 5-Hydroxyvalproic Acid

5-Hydroxyvalproic acid (5-OH-VPA) belongs to a class of organic compounds known as hydroxy fatty acids. These are fatty acids that have a hydroxyl group in their chain. This toxic metabolite of VPA is formed in the human body (mainly in hepatocytes) by P-oxidation involving three cytochrome P450 isoenzymes: CYP2A6; CYP2C9; and CYP2B6 [97]. 5-OH-VPA has been shown to have a hepatotoxic effect, which also explains the importance of PGx (Appendix A) in identifying patients at risk.

6.4. 3-Oxovalproic Acid

3-Oxovalproic acid belongs to a class of organic compounds known as short-chain keto acids and their derivatives [97]. These keto acids do have an alkyl chain containing fewer than six carbon atoms. This toxic metabolite of VPA has a suspected toxic effect [97], and studies of its toxic potential are ongoing.

6.5. Valproyl-CoenzymeA

Valproil-CoA is the main metabolite of VPA produced in mitochondria. In addition, valproyl-CoA affects the activity of succinate-CoA ligase, especially the ATP-specific enzyme [115]. Inhibition of activated succinate-CoA ligase can affect nucleoside-diphosphate kinase activity, causing an imbalance of nucleotides in mitochondria and, ultimately, depletion of mitochondrial DNA [115]. This may explain the hepatotoxic potential of this VPA metabolite and the development of liver failure in patients with a deficiency in the mitochondrial DNA replicase system, such as gamma polymerase [115].

6.6. 2-N-Propyl-4-Oxopentanoic Acid

2-N-propyl-4-oxopentanoic acid belongs to a class of organic compounds known as gamma-keto acids and their derivatives. These are organic compounds containing an aldehyde substituted by a keto group on the C4 carbon atom [114]. It has been previously shown that 4-keto-VPA causes abnormal and irregular growth of the neural tube neuroepithelium, microvesicular steatosis, and reduces ammonia excretion through inhibition of carbamoyl phosphate [114], which significantly increases the risk of VPA-induced ADRs, including teratogenesis.

6.7. 4-Ene-Valproic Acid

4-Ene-valproic acid (4-ene-VPA) belongs to a class of organic compounds known as methyl-branched fatty acids [97]. These are fatty acids with an acyl chain having a methyl branch. They are usually saturated and contain only one or a few methyl groups. However, branches other than the methyl branch may also be present. It is formed in the body as a result of P-oxidation of valproates in the liver with the participation of isoenzymes CYP2A6, CYP2C9, and CYP2B. As previously stated, 4-ene-VPA is a highly hydrophobic molecule that is considered a toxic metabolite of VPA in patients with epilepsy and psychiatric disorders; this metabolite has a hepatotoxic effect, causing hyperammonemia and Reye’s syndrome in rats [97,114]. At the same time, 4-ene-VPA is being actively studied as a potential therapeutic metabolite in multiple tumors, autoimmune diseases, and AIDS [114].

6.8. 2-Propyl-2,4-Pentadienoic Acid

2-Propyl-2,4-pentadienoic acid (also known as 2,4-diene-VPA) is metabolized in the human body by glucuronidation [97]. This toxic metabolite of VPA binds reversibly to the α-subunit of the trifunctional protein, inhibiting β-oxidation. As a result, 2,4-diene VPA reduces levels of acetyl-CoA and glutamate, which are required for the synthesis of N-acetylglucosamine 6-phosphate deacetylase. A decrease in the level of N-acetylglucosamine 6-phosphate deacetylase, a cofactor for the CPSI enzyme, leads to impaired urea production and the subsequent development of VPA-associated hyperammonemia and hyperammonemic encephalopathy, cerebral edema [68]. This VPA metabolite is formed in the body by P-oxidation with the help of CYP2A6 and CYP2C9 isoenzymes [114].

6.9. 2-Ene-Valproic Acid

2-Ene-valproic acid (2-ene-VPA) is a product of VPA metabolism in the liver by β-oxidation. This toxic metabolite of VPA belongs to a class of organic compounds known as methyl-branched fatty acids [97]. It is known that 2-ene-VPA is a neurotoxic metabolite of VPA [114]; therefore, an increase in the level of this metabolite in the blood and in the central nervous system can lead to VPA-induced neurological ADRs (diplopia, ataxia, seizures) and mental ADRs (cognitive disorders, behavioral disorders).

6.10. Valproylcarnitine

Valproyl-carnitine belongs to a class of organic compounds known as acylcarnitines [116]. These are organic compounds containing a fatty acid, the carboxylic acid of which is attached to carnitine through an ester bond. VPA reduces the concentration of L-carnitine in the body through the excretion of their acyl-carnitines, pivaloyl-carnitine, and valproyl-carnitine in the urine [113], which causes the development of such ADRs as: lipid myopathy; hypoglycemia; fatty liver disease; and hyperammonemia.

7. Role of Toxic Metabolites of Valproic Acid

As is known, the toxicity of VPA and its metabolites refers to how toxic or harmful these substances can be and which organ(s) or tissue(s) are most sensitive to their toxic effects (Figure 4). This occurs when there are too many toxic VPA metabolites in the blood as a result of a primary (genetically determined) disturbance of drug metabolism, which leads to the ADRs of toxic metabolites on the human body. In addition, an increase in the level of toxic metabolites of VPA in biological fluids, organs, and tissues can occur when too high a dose of VPA is administered, when there are secondary VPA metabolism violations in the liver (for example, in elderly patients or in patients with liver failure), or when there is a violation of toxic VPA metabolites excretion from the bloodstream through the kidneys as a result of their primary or secondary breakdown.

One of the reasons for the development of ADRs of toxic VPA metabolites is the induced carnitine deficiency due to the excretion of valproyl-carnitine by the urine [117] and, as a consequence, impaired mitochondrial oxidation [118]. In addition, some VPA metabolites are toxic to isolated liver mitochondria with an effect similar to VPA accumulation [119]. As a result, hepatic microvesicular steatosis develops, as well as hyperammonemia, anion gap acidosis, metabolic encephalopathy, and hepatotoxicity [24]. Toxic VPA metabolites also induce teratogenicity by downregulating the IGF2R, RGS4, COL6A3, EDNRB, and KLF6 genes, which are related to neural tube defects [90] in children whose mothers took valproate during pregnancy [120].

In addition, the co-administration of other drugs can lead to an increase in the level of toxic metabolites of VPA in the blood. For example, acetylsalicylic acid increases plasma levels of VPA and its toxic metabolites by competing for plasma protein binding [121,122]. This is demonstrated in a study where the authors profiled VPA and VPA metabolites in urine before (day 1) and after (day 2) administration of antipyretic doses of acetylsalicylic acid in seven subjects with stable serum VPA levels according to therapeutic drug monitoring. Of the 13 metabolites analyzed by gas-liquid chromatography / mass spectrometry, levels of (E)-2-ene-VPA and 3-keto-VPA decreased significantly on day 2, while levels of VPA (glucuronide) and 4-ene-VPA went up significantly. VPA β-oxidation pathway activity, including (E)-2-ene-VPA, 3-OH-VPA and 3-keto-VPA, decreased from 24.5% ± 10.3% of total metabolites excreted on day 1, up to 8.3% ± 4.2% on day 2. Thus, a slowdown in VPA metabolism was shown with the additional appointment of acetylsalicylic acid by 66% of the original. At the same time, the content of VPA glucuronide increased from 50.5 ± 12.6% on day 1 to 65.5 ± 14% of the total amount of excreted metabolites on day 2, which is 30% more than the initial value. The ratio of VPA glucuronide “day 2 / day 1” was statistically significantly correlated with the ratio of the average free fraction of VPA “day 2/day 1” (r = 0.9424; p-value = 0.005). Inhibition of VPA β-oxidation by salicylate was sufficient to counterbalance the increased elimination of VPA in the form of its conjugates and explain why the total clearance of VPA after administration of this salicylate remains unchanged, even if the free fraction of VPA increases. Metabolic profiles indicate that salicylate probably inhibits VPA β-oxidation by reducing the formation of valproil-CoA [123]. However, such a toxic effect of the combination of acetylsalicylic acid and VPA may have a synergistic killing effect on the proliferation and apoptosis of hepatocellular cancer cells [124]. These studies open new perspectives and re-evaluate the clinical role of toxic VPA metabolites.

Also, toxic metabolites of VPA can cause drug-induced parkinsonism. This VPA-induced neurological ADR occurs equally in males and females, but is more common in the elderly (usually over 55 years of age) and is not dose-dependent [125]. According to previous studies, in all patients with VPA-induced parkinsonism, serum levels of VPA were in the therapeutic range (50–100 µg/mL), although the level of neurotoxic metabolites of VPA was not investigated. Chronic maintenance doses of VPA appear to result in increased levels of its toxic metabolites in the central nervous system, which is the main cause of this ADR. This hypothesis is supported by the results of studies demonstrating that the symptoms of VPA-induced parkinsonism usually decrease within a few weeks and completely disappear within a few months after stopping valproate. However, if extrapyramidal symptoms do not improve, then it is likely that the toxic metabolites of VPA have unlocked the hidden potential (genetic predisposition) for the development of Parkinson’s disease. Such patients may benefit from the addition of levodopa therapy [126], although the mechanism of how toxic VPA metabolites cause the development of parkinsonism remains unknown [127].

8. Discussion

The results of the technical and scientific advances of the 21st century demonstrate the astonishing advances that have made it possible to detect, identify, and quantify a large number of VPA metabolites in biological fluids and/or tissues of the human body. In addition, similar technical and scientific advances have occurred in DNA sequencing. These advances explain the importance of moving beyond the study of individual VPA metabolites and individual genetic polymorphisms to the study of probably hundreds of VPA metabolites and even more genomic variants in a single cell or subject. In recent decades, it has become possible to combine and integrate large datasets from various “-omics” methods to improve our understanding of the molecular basis of VPA-induced ADRs risks and/or drug response phenotypes to VPA and its therapeutic and toxic metabolites.

Undoubtedly, the new concept that metabolomic profiles (metabolotypes) can help predict the phenotypes of the expected therapeutic response to VPA and its therapeutic and toxic metabolites [128] seems to be relevant and important from a clinical point of view. Obviously, the variation in metabotype among patients taking valproate reflects variations due to both environmental factors (eg., a dose of valproate or duration of therapy) and genetic factors. As a result, pharmacometabolomics may provide the practicing neurologist, psychiatrist, and other professionals (eg., pharmacologist, oncologist, or hematologist) with a better understanding of the overall molecular variability that contributes to individual differences in therapeutic response to valproate than pharmacogenomics alone [129].

Thus, the pharmacogenomics and pharmacometabolomics of VPA are closely interconnected. Known genetic biomarkers of changes in VPA metabolism are summarized by us and presented in Appendix A. VPA metabolites can be divided into three groups: therapeutic, neutral, and toxic (Figure 5). However, such an assessment should be considered conditional, since the effect of VPA metabolites will depend on the goal of the therapeutic strategy for a particular disease in a particular patient (individual). This is important to remember because the same VPA metabolite can be considered both therapeutic (e.g., in massive tumors and neurodegenerative diseases) and toxic (e.g., in congenital malformations of the fetus in women taking VPA during pregnancy).

However, despite the fact that the application of an integrated (pharmacogenomics + pharmacometabolomics) approach represents significant progress in the last decade, most researchers and clinicians are still focused on using only one -omics method in real clinical practice using VPA without evaluating the potential contribution of various -omics with potential new ideas that could provide such a combined approach [128]. For example, pharmacogenetics focusing on individual genes and their influence on the expected drug response to valproate and other drugs, as well as the assessment of individual (personalized) risk of VPA-induced ADRs (Figure 6), has turned into pharmacogenomics using genome-wide studies of clinical drug phenotypes [130].

At the same time, the addition of data from other studies, i.e., the further evolution of pharmacogenomics into “pharmacoomics”, including the pharmacogenomics and pharmacometabolomics of VPA (Figure 7), is now possible and may be required in real clinical practice, especially for studies of drug therapy with VPA drugs for complex diseases with many different underlying issues at the basis of pathophysiological mechanisms, which include, first of all, neurological diseases and mental disorders [131].

The pharmacogenomics and pharmacometabolomics of VPA can be particularly useful when used together (Figure 8).

The genome is consistent with various cell and tissue types in humans, although epigenomics and the regulation of genomic expression are highly dependent on cells and tissues, and since pharmacometabolomics is an “-omic” type of data that appears to be most closely associated with clinical phenotypes [128,132], the integration of pharmacogenomics and pharmacometabolomics may be particularly useful in achieving an optimal balance between the efficacy and safety of VPA. The use of genomic methods, such as genome-wide association studies, has been very successful [132], especially in situations where the phenotype of VPA drug metabolism in the liver in a particular patient is clearly defined (for example, P-oxidation involving cytochrome P450 isoenzymes) [37]. However, the use of pharmacogenomic methods alone was less successful when applied to complex and less defined phenotypes (for example, the «P-oxidation + glucuronidation + acetylation») of VPA metabolism in patients with common and rare neurological diseases and mental disorders. It is for this reason that the use of our proposed personalized research strategy, which starts with VPA pharmacogenomics and allows genomic data to “guide” the VPA pharmacometabolomics study, can be particularly useful in selected clinical situations.

For example, with a genetically determined increase in the level of toxic metabolites associated with impaired mitochondrial oxidation, in which VPA inhibits carnitine biosynthesis by reducing the concentration of alpha-ketoglutarate and can contribute to the development of carnitine deficiency, L-carnitine preparations are the drug of choice for the correction of VPA-induced ADRs. It is known that L-carnitine is an amino acid derivative that is an important cofactor in the beta-oxidation of fatty acids [117], including VPA. L-carnitine is increasingly recommended for the treatment of hyperammonemia caused by VPA. It is a corrective drug for VPA-induced hyperammonemia and other ADRs of this drug and acts as an acceptor of toxic metabolites such as valproil-CoA. Valproil-carnitine ester is subsequently excreted in the urine. As a result of the normalization of the ratio of acyl-CoA/CoA-OH in mitochondria, the reverse inhibition of β-oxidation occurs, which, in turn, leads to the formation of acetyl-CoA. The formation of acetyl-CoA allows the formation of NAGA and the subsequent detoxification of ammonia to urea, thereby correcting hyperammonemia. L-carnitine supplementation is well tolerated, and reports of possible ADRs are rare [38,133].

With an increase in the level of toxic metabolites caused by a violation of P-oxidation of VPA, in order to correct VPA-induced ADRs, it is recommended to prescribe ursodeoxycholic acid (UDCA) drugs. It is known that UDCA is a primary bile acid that is easily synthesized in the human body [134]. In addition, UDCA has a direct cytoprotective effect due to the stabilization of the hepatocyte membrane and improvement of mitochondrial oxidative phosphorylation [135]. In addition, UDCA prevents the transition of the mitochondrial membrane permeability and promotes a decrease in the VPA-induced high serum levels of alanine aminotransferase and aspartate aminotransferase [136]. An immunomodulatory effect was found in UDCA; its mechanism is explained by a decrease in the expression of HLA class I proteins in hepatocytes and the anti-apoptotic effect of UDCA due to a decrease in the production of reactive oxygen species [137].

With an increase in the level of metabolites associated with impaired glucuronidation [138] of VPA drugs, phenobarbital is a well-established corrector of VPA-induced ADRs [139]. Phenobarbital has been used for many years for detoxification in patients suffering from various forms of sedative dependence, as well as intoxication with psychotropic drugs. There is also a reliable clinical test using pentobarbital to assess the tolerance of psychotropic drugs and thus determine the dosing regimen to replace phenobarbital [140].

Thus, the presented personalized approach to the cumulative risk assessment for the development and correction of ADRs associated with an increase in the level of toxic VPA metabolites in the blood and urine allows clinicians to adapt previously known classical methods of treatment using valproate to individual patients from the standpoint of pharmacogenomics and pharmacometabolomics. It should be recognized that practicing neurologists and psychiatrists often lack adequate objective data to truly personalize psychopharmacotherapy [141]. In the future, a new era of personalized medicine may enable the effective use of a combination of pharmacogenomics and pharmacometabolomics to identify major molecular pathways and genetic variations in candidate genes encoding key VPA biotransformation enzymes. This may help us explain the individual phenotypes of patients with neurological disease and psychiatric disorders and ensure that an optimal balance is achieved between the efficacy and safety of valproates.

So, the term “pharmacometabolomics-informed pharmacogenomics”, first proposed 10 years ago (in 2012) [142], represents one promising strategy at the intersection of pharmacometabolomics and pharmacogenomics and gives us the key to understanding the mechanisms of formation, action, and utilization of VPA metabolites. The alternative approach proposed by us, “pharmacogenetics-informed pharmacometabolomics”, makes it possible to individually predict the risk of elevated levels of toxic VPA metabolites based on lifelong unchanged PGx results that do not change throughout a person’s life and optimize the choice and scope of metabolic studies (in particular, gas-liquid chromatography/mass spectrometry). This is an important aspect of the movement towards 4P medicine (prediction, prevention, personalization, and participation) [143]. In addition, it is promising and important to compare the results of pharmacogenomic and pharmacometabolic studies of therapeutic and toxic metabolites of VPA in the future in order to identify genetic variants that contribute to metabolomic variation but are not dependent on drug response to valproate.

9. Limitations

Within the limits of restrictions, we would like to note that the search was carried out only in PubMed, Web of Science, Springer, Google Scholar, and e-Library databases. We have focused on the well-known and potential role of therapeutic and toxic metabolites in the efficacy and safety of VPA in this narrative review, but we have not analyzed clinical randomized trials of the efficacy of valproates, including VPA and its salts. It is possible that the level of therapeutic and toxic metabolites of VPA in human biological fluids was also studied in some of these trails. Undoubtedly, it would be interesting to prepare a systematic review in the future.

10. Conclusions

Valproates are psychotropic drugs that are widely used for neurological diseases and mental disorders. In addition, the indications for the appointment of valproate have been expanding in recent years in connection with the study of the mechanism of action of therapeutic and toxic metabolites of VPA in the human body. The most studied pathways of VPA metabolism include: β-oxidation in the tricarboxylic acid cycle (acetylation); oxidation involving cytochrome P-450 isoenzymes in the liver (P-oxidation); and glucuronidation. The complex metabolism of VPA explains the diversity of its active and inactive metabolites, which have therapeutic, neutral, or toxic effects. It is known that some active metabolites of VPA may have a stronger clinical effect than VPA itself. In future studies of VPA metabolites, it is important to use pharmacometabolomics together with pharmacogenomics to improve understanding of the biochemical pathways for the therapeutic and toxic metabolites formation, which may play an important role in individual variations in the phenotypes of the expected drug response to VPA, as well as in individual variations in the metabotypes associated with the risk of developing VPA-induced ADRs. An illustration of the importance of this approach is the proposed personalized approach to the cumulative assessment of the risk of developing VPA-induced ADRs and the choice of tactics for their correction, depending on the pharmacogenetic profile of the patient and the level of therapeutic and toxic VPA metabolites in human body fluids (blood, urine), and possibly saliva.

Author Contributions

Conceptualization, N.A.S. and R.F.N.; methodology, N.A.S. and K.S.A.; software, V.V.G. and A.K.K.; validation, E.N.B., A.K.K., V.V.G. and E.A.D.; formal analysis, M.M.P., A.R.A. and G.A.S.; investigation, V.V.G. and M.A.-Z.; resources, V.V.G.; data curation, E.N.B. and A.K.K.; writing—original draft preparation, V.V.G., E.A.D. and A.K.K.; writing—review and editing, N.A.S. and R.F.N.; visualization, V.V.G. and N.A.S.; supervision, N.A.S. and R.F.N.; project administration, N.A.S. and R.F.N. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. 3D scheme of valproic acid molecule.

Figure 2. Mechanism of action of valproic acid. Note: GOT—glutamate-oxaloacetate transaminase; GABA—gamma-aminobutyric acid; ABAT—4-aminobutyrate aminotransferase; ALDHA1—aldehyde dehydrogenase; SUCL—succinate-CoA ligase; OGDH—oxoglutarate dehydrogenase; HDAC—histone deacetylase.

Figure 3. The most common clinical therapeutic effects of valproic acid metabolites. Notes: (1) blockade of voltage-dependent sodium, calcium, and potassium channels; (2) neuroprotection; (3) inhibition of signaling pathways in cancer cells; (4) modulation of neurogenesis; (5) angioprotection.

Figure 4. The most common clinical effects of toxic metabolites of valproic acid. Notes: (1) hepatotoxicity; (2) teratogenicity; (3) dermatitis, alopecia; (4) thrombocytopenia; (5) depression of the central nervous system.

Figure 5. Clinical classification of valproic acid metabolites: green sign–therapeutic metabolites; yellow sign–neutral metabolites; red sign–toxic metabolites.

Figure 6. Pharmacogenomic approach to assessing the risk of development- and correction-adverse reactions induced by valproic acid metabolites.

Figure 7. Pharmacometabolomic approach to assessing the risk of development and correction of adverse reactions induced by valproic acid metabolites.

Figure 8. «The pharmacogenetics-informed pharmacometabolomics» approach to assessing the risk of development and correction of adverse drug reactions induced by valproic acid metabolites.

Table 1

Serum and Plasma Therapeutic Metabolites of Valproic Acid.

Metabolite	HMDB Number	ClinicalEffect	References
VPA-glucuronide	0000901	TherapeuticNeutral	[18,98]
2-n-propyl-2-pentenoic acid	0013902	Therapeutic	[99]
Valproyl-CoA	0013115	Therapeutic	[100]
4-ene-VPA	0013897	Therapeutic	[101]
2-ene-VPA	0013902	Therapeutic	[101]

Note: CoA–coenzyme A; HMDB–The Human Metabolome Database [https://hmdb.ca/ (accessed on 12 November 2022)]; VPA–valproic acid.

Table 2

Urinary Therapeutic Metabolites of Valproic Acid.

Metabolite	HMDB Number	ClinicalEffect	References
VPA-glucuronide	0000901	TherapeuticNeutral	[18,98]
2-n-propyl-2-pentenoic acid	0013902	Therapeutic	[99]
4-ene-VPA	0013897	Therapeutic	[100]
2-ene-VPA	0013902	Therapeutic	[101]

Note: CoA–coenzyme A; HMDB–The Human Metabolome Database [https://hmdb.ca/ (accessed on 12 November 2022)]; VPA–valproic acid.

Table 3

Serum and Plasma Toxic Metabolites of Valproic Acid.

Metabolite	HMDB Number	ClinicalEffect	References
3-hydroxy valproic acid	0013899	Toxic	[95]
4-hydroxy valproic acid	0013900	Toxic	[112]
5-hydroxy valproic acid	0013898	Toxic	[95]
Valproyl Coenzyme A	0013115	Toxic	[95]
Valproylcarnitine	0259757	Toxic	[113]
2-n-propyl-2-pentenoic acid	0013902	Toxic	[114]
2-n-propyl-4-oxopentanoic acid	0060683	Toxic	[114]
2-Propyl-2,4-pentadienoic acid	0060682	Toxic	[23]
3-oxovalproic acid	0060685	Toxic	[95]
4-ene-valproic acid	0013897	Toxic	[95]
2-ene-valproic acid	0013902	Toxic	[114]
2-ene-valproic acid Coenzyme A	0060714	Toxic	[38]
3-ene-valproic acid Coenzyme A	0060740	Toxic	[95]
2-diene-valproic acid glucuronide		Toxic	[23]

Note: HMDB–The Human Metabolome Database [https://hmdb.ca/ (accessed on 12 November 2022)].

Table 4

Urinary Toxic Metabolites of Valproic Acid.

Metabolite	HMDB Number	ClinicalEffect	References
3-hydroxy valproic acid	0013899	Toxic	[95]
4-hydroxy valproic acid	0013900	Toxic	[112]
5-hydroxy valproic acid	0013898	Toxic	[95]
3-oxovalproic acid	0060685	Toxic	[95]
2-ene-valproic acid	0013902	Toxic	[114]
2-ene-valproic acid Coenzyme A	0060714	Toxic	[38]
3-ene-valproic acid Coenzyme A	0060740	Toxic	[95]
2-propyl-2,4-pentadienoic acid	0060682	Toxic	[23]
2-n-propyl-4-oxopentanoic acid	0060683	Toxic	[114]

Note: HMDB–The Human Metabolome Database [https://hmdb.ca/ (accessed on 12 November 2022)].

Appendix A

Table A1

Genetic Biomarkers of Valproic Acid Metabolism Retardation.

Gene(OMIM Number)	ChromosomalLocation	Protein	SNV (RSID)	References
			Oxidation
CYP2A6(122720)	19q13.2	Isoenzyme2A6 ofcytochrome P450	*2 (rs1801272)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*5 (rs5031017)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.The heterozygous GT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			rs (111033610)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			rs28399447	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*6 (rs4986891)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*7 (rs5031016)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*9A (s28399433)	The homozygous GG genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*11 (rs28399447)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*17 (rs28399454)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*18 (rs1809810)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*19 (rs1809810 )	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*23 (rs56256500)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*24A (rs143731390)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*24A (rs72549435)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*26 (rs59552350)	The homozygous GG genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*26 (rs4986891)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*27 (rs2839944)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*27 (rs28399445)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*35 (rs143731390)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*41 (rs140471703 7)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
CYP2C9(601130)	10q23.33	Isoenzyme2A6 ofcytochrome P450	*2 (rs1799853)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*3 (rs1057910)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*4 (rs56165452)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*5 (rs28371686)	The homozygous GG genotype– a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*8 (rs7900194)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*11 (rs28371685)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*13 (rs72558187)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*15 (rs72558190)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
CYP2B6(123930)	19q13.2	Isoenzyme2A6 ofcytochrome P450	*5A (rs3211371)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*8 (rs12721655)	The homozygous GG genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*18 (rs28399499)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*27 (rs36079186)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*28 (rs34097093)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
CYP3A4(124010)	7q22.1	Isoenzyme2A6 ofcytochrome P450	*3 (rs4986910)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*17 (rs4987161)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*18 (rs28371759)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*20 (rs67666821)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
CYP2C19(124020)	10q23.33	Isoenzyme2A6 ofcytochrome P450	*2 (rs4244285)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*4 (rs28399504)	The homozygous GG genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*5 (rs56337013)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*8 (rs41291556)	The homozygous CC genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
CYP2D6(124030)	22q13.2	Isoenzyme2A6 ofcytochrome P450	*4 (rs3892097)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*4F (rs3892097)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*4G (rs3892097)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*4H (rs3892097)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*8 (rs5030865)	The homozygous AA genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*10 (rs1065852)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			*17 (rs28371706)	The homozygous TT genotype—a biomarker for a high risk of impaired VPA oxidation and increased levels of toxic metabolites.
			Glucuronization
UGT1A1(191740)	2q37.1	UDP-glucuronosyltransferase	rs111033539	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033540	The homozygous GG genotype (Crigler Najjar syndrome, Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033541	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A3(606428)	2q37.1	UDP-glucuronosyltransferase	rs111033541	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs10929302	A biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033539	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
191740			rs111033540	The homozygous GG genotype (Crigler-Najjar syndrome, Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs11891311	Allele G—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites (встрeчaeтся чащe у азиатoв).
			rs281865418	The homozygous GG; AA genotypes (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs28934877	The homozygous GG; CC genotypes (Crigler-Najjar syndrome Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs3064744	The homozygous TA; TA genotypes—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs35003977	The homozygous GG genotype (Hyperbilirubinemia Gilbert’s syndrome; Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs34993780	The homozygous AA; CC; GG genotypes (Lucey-Driscoll syndrome; Crigler-Najjar syndrome; Hyperbilirubinemia)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs35350960	The homozygous AA; TT genotypes—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs3755319	The homozygous GG genotype (Lucey-Driscoll syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs386576623	The homozygous GG genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs4124874	The homozygous CC genotype (Gilbert syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs4148323	The homozygous AA genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776761	The homozygous AA genotype (Gilbert syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776763	The homozygous TT genotype (Gilbert syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776764	The homozygous CC genotype (Gilbert syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776765	The homozygous TT genotype (Crigler-Najjar syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776766	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587784535	The homozygous TT genotype (Hyperbilirubinemia) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs62625011	The homozygous AA genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs6742078	The homozygous TT genotype (16% bilirubin levels increased risk of gallstones) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551341	The homozygous AA genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551348	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551349	The homozygous TT genotype (Crigler Najjar syndrome Gilbert’s syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551350	The homozygous TT genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551351	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551353	The homozygous TT genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs755218546	The homozygous AA genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs886039771	The homozygous AA genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A4(606429)	2q37.1	UDP-glucuronosyltransferase	rs111033539	The homozygous TT genotype (Crigler-Najjar syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033540	The homozygous AA genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033541	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs2011425	The homozygous GG; AA genotypes—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs281865418	The homozygous GG; AA genotypes (Crigler-Najjar syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs28934877	The homozygous GG; CC genotypes (Crigler-Najjar syndrome; Gilbert’s syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs34993780	The homozygous GG; AA; CC genotypes (Lucey-Driscoll syndrome Crigler-Najjar syndrome Hyperbilirubinemia) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs35003977	The homozygous GG genotype (Hyperbilirubinemia Gilbert’s syndrome; Crigler-Najjar syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs3755319	The homozygous GG genotype (Lucey-Driscoll syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs386576623	The homozygous GG genotype (Crigler-Najjar syndrome Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs4124874	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776761	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776763	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776764	The homozygous CC genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776765	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776766	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587784535	The homozygous TT genotype (Hyperbilirubinemia) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs62625011	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs6742078	The homozygous TT genotype (+16% bilirubin levels increased risk of gallstones) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551341	The homozygous AA genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551348	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551349	The homozygous TT genotype (Crigler Najjar syndrome Gilbert’s syndrome) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551350	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551351	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551353	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs755218546	The homozygous AA genotype —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs797046090	Genotypes CGTA;CGTA (Hyperbilirubinemia) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs886039771	The homozygous AA genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A6(606431)	2q37.1	UDP-glucuronosyltransferase	rs111033539	The homozygous TT genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033540	The homozygous GG genotype (Crigler-Najjar syndrome Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033541	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs2011425	The homozygous AA; GG genotypes—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs281865418	The homozygous AA; GG genotypes (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs28934877	The homozygous CC; GG genotypes (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs34993780	The homozygous AA; GG genotypes (Lucey-Driscoll syndrome; Crigler-Najjar syndrome Hyperbilirubinemia)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs35003977	The homozygous GG genotype (Hyperbilirubinemia Gilbert’s syndrome; Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs3755319	The homozygous GG genotype (Lucey-Driscoll syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs386576623	The homozygous GG genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs4124874	The homozygous GG genotype (Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776761	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776763	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776764	The homozygous CC genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776765	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776766	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587784535	The homozygous TT genotype (Hyperbilirubinemia)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs62625011	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs6742078	The homozygous TT genotype (+16% bilirubin levels increased risk of gallstones) —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551341	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs72551348	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A8(606433)	2q37.1	UDP-glucuronosyltransferase	rs111033539	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033540	The homozygous GG genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033541	The homozygous GG genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs2011425	The homozygous GG; AA genotypes—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs281865418	The homozygous GG; AA genotypes (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs28934877	The homozygous GG; CC genotypes (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs34993780	The homozygous AA; GG; CC genotypes (Lucey-Driscoll syndrome; Crigler-Najjar syndrome Hyperbilirubinemia)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs35003977	The homozygous GG genotype (Hyperbilirubinemia Gilbert’s syndrome; Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs3755319	The homozygous GG genotype (Lucey-Driscoll syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs386576623	The homozygous GG genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs4124874	The homozygous CC genotype (Gilbert syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs587776761	The homozygous AA genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A9(606434)	2q37.1	UDP-glucuronosyltransferase	rs111033539	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033540	The homozygous GG genotype (Crigler-Najjar syndrome; Gilbert’s syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs111033541	The homozygous GG genotype (Crigler-Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT1A10(606435)	2q37.1	UDP-glucuronosyltransferase	rs72551350	The homozygous TT genotype (Crigler Najjar syndrome)—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs17863783	The homozygous TT genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT2B7(600068)	4q13.2	UDP-glucuronosyltransferase	rs12233719	The homozygous TT genotype —a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs7438135	The homozygous GG genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs7439366	The homozygous TT genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			rs7668258	The homozygous TT genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
UGT2B15(600069)	4q13.2	UDP-glucuronosyltransferase	rs1902023	The homozygous GG genotype—a biomarker for a high risk of impaired VPA glucuronization and increased levels of toxic metabolites.
			Acetylation
ACADSB(600301)	10q26.13	Acyl-CoA Dehydrogenase	rs58639322	The homozygous TT genotype —a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs142095937	The homozygous GG genotype —a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs137852649	The homozygous TT genotype—a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs760791287	The homozygous TT genotype—a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs147936696	The homozygous AA genotype —a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs387906409	The homozygous AA genotype—a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
ACSM1(614357)	16p12.3	Acyl-CoA Synthetase 1	rs163234	The homozygous AA genotype —a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs151222	The homozygous CC genotype —a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.
			rs433598	The homozygous CC genotype—a biomarker for a high risk of impaired VPA acetylation and increased levels of toxic metabolites.

Note: CYP—Cytochrome; UGT—Uridine diphosphate glucosyltransferase; ACADSB—Acyl-Coenzyme A dehydrogenase short/branched chain; ACSM1- Acyl-Coenzyme A synthetase medium chain family member 1; CoA—Coenzyme A; UDP- Glucuronosyltransferase.

References

1. Woolley, B.; Mills, J. Versatile Valproic Acid. Issues Ment. Health Nurs.; 2022; 43, pp. 1072-1074. [DOI: https://dx.doi.org/10.1080/01612840.2022.2122431] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36136610]

2. Romoli, M.; Mazzocchetti, P.; D’Alonzo, R.; Siliquini, S.; Rinaldi, V.E.; Verrotti, A.; Calabresi, P.; Costa, C. Valproic acid and epilepsy: From molecular mechanisms to clinical evidences. Curr. Neuropharmacol.; 2019; 17, pp. 926-946. [DOI: https://dx.doi.org/10.2174/1570159X17666181227165722] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30592252]

3. Peterson, G.M.; Naunton, M. Valproate: A simple chemical with so much to offer. J. Clin. Pharm. Ther.; 2005; 30, pp. 417-421. [DOI: https://dx.doi.org/10.1111/j.1365-2710.2005.00671.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16164485]

4. Macritchie, K.A.; Geddes, J.R.; Scott, J.; Haslam, D.R.; Goodwin, G.M. Valproic acid, valproate and divalproex in the maintenance treatment of bipolar disorder. Cochrane Database Syst. Rev.; 2001; 3, CD003196. [DOI: https://dx.doi.org/10.1002/14651858.CD003196]

5. Iannaccone, T.; Sellitto, C.; Manzo, V.; Colucci, F.; Giudice, V.; Stefanelli, B.; Iuliano, A.; Corrivetti, G.; Filippelli, A. Pharmacogenetics of Carbamazepine and Valproate: Focus on Polymorphisms of Drug Metabolizing Enzymes and Transporters. Pharmaceuticals; 2021; 14, 204. [DOI: https://dx.doi.org/10.3390/ph14030204]

6. Linde, M.; Mulleners, W.M.; Chronicle, E.P.; McCrory, D.C. Valproate (valproic acid or sodium valproate or a combination of the two) for the prophylaxis of episodic migraine in adults. Cochrane Database Syst. Rev.; 2013; 6, CD010611. [DOI: https://dx.doi.org/10.1002/14651858.CD010611]

7. Guo, A.; Li, J.; Luo, L.; Chen, C.; Lu, Q.; Ke, J.; Feng, X. Valproic acid mitigates spinal nerve ligation-induced neuropathic pain in rats by modulating microglial function and inhibiting neuroinflammatory response. Int. Immunopharmacol.; 2021; 92, 107332. [DOI: https://dx.doi.org/10.1016/j.intimp.2020.107332]

8. Wang, Y.; Xia, J.; Helfer, B.; Li, C.; Leucht, S. Valproate for schizophrenia. Cochrane Database Syst. Rev.; 2016; 11. [DOI: https://dx.doi.org/10.1002/14651858.CD004028.pub4]

9. Lambie, D.G.; Johnson, R.H.; Vijayasenan, M.E.; Whiteside, E.A. Sodium valproate in the treatment of the alcohol withdrawal syndrome. Aust. N. Z. J. Psychiatry; 1980; 14, pp. 213-215. [DOI: https://dx.doi.org/10.3109/00048678009159381]

10. Armeanu, S.; Bitzer, M.; Lauer, U.M.; Venturelli, S.; Pathil, A.; Krusch, M.; Kaiser, S.; Jobst, J.; Smirnow, I.; Wagner, A. et al. Natural killer cell-mediated lysis of hepatoma cells via specific induction of NKG2D ligands by the histone deacetylase inhibitor sodium valproate. Cancer Res.; 2005; 65, pp. 6321-6329. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-04-4252]

11. Kristensen, Ø.; Lølandsmo, T.; Isaksen, A.; Vederhus, J.K.; Clausen, T. Treatment of polydrug-using opiate dependents during withdrawal: Towards a standardisation of treatment. BMC Psychiatry; 2006; 6, 54. [DOI: https://dx.doi.org/10.1186/1471-244X-6-54]

12. Kraus, S.W.; Voon, V.; Potenza, M.N. Should compulsive sexual behavior be considered an addiction?. Addiction; 2016; 111, pp. 2097-2106. [DOI: https://dx.doi.org/10.1111/add.13297]

13. Chateauvieux, S.; Morceau, F.; Dicato, M.; Diederich, M. Molecular and therapeutic potential and toxicity of valproic acid. J. Biomed. Biotechnol.; 2010; 2010, 479364. [DOI: https://dx.doi.org/10.1155/2010/479364]

14. Lehrman, G.; Hogue, I.B.; Palmer, S.; Jennings, C.; Spina, C.A.; Wiegand, A.; Landay, A.L.; Coombs, R.W.; Richman, D.D.; Mellors, J.W. et al. Depletion of latent HIV-1 infection in vivo: A proof-of-concept study. Lancet; 2005; 366, pp. 549-555. [DOI: https://dx.doi.org/10.1016/S0140-6736(05)67098-5]

15. Gurpur, P.B.; Liu, J.; Burkin, D.J.; Kaufman, S.J. Valproic acid activates the PI3K/Akt/mTOR pathway in muscle and ameliorates pathology in a mouse model of Duchenne muscular dystrophy. Am. J. Pathol.; 2009; 174, pp. 999-1008. [DOI: https://dx.doi.org/10.2353/ajpath.2009.080537]

16. Duenas-Gonzalez, A.; Candelaria, M.; Perez-Plascencia, C.; Perez-Cardenas, E.; de la Cruz-Hernandez, E.; Herrera, L.A. Valproic acid as epigenetic cancer drug: Preclinical, clinical and transcriptional effects on solid tumors. Cancer Treat. Rev.; 2008; 34, pp. 206-222. [DOI: https://dx.doi.org/10.1016/j.ctrv.2007.11.003]

17. Zhu, M.M.; Li, H.L.; Shi, L.H.; Chen, X.P.; Luo, J.; Zhang, Z.L. The pharmacogenomics of valproic acid. J. Hum. Genet.; 2017; 62, pp. 1009-1014. [DOI: https://dx.doi.org/10.1038/jhg.2017.91]

18. Perucca, E. Pharmacological and therapeutic properties of valproate: A summary after 35 years of clinical experience. CNS Drugs.; 2002; 16, pp. 695-714. [DOI: https://dx.doi.org/10.2165/00023210-200216100-00004]

19. Silva, M.F.; Aires, C.C.; Luis, P.B.; Ruiter, J.P.; IJlst, L.; Duran, M.; Wanders, R.J.; de Almeida, I.T. Valproic acid metabolism and its effects on mitochondrial fatty acid oxidation: A review. J. Inherit. Metab. Dis.; 2008; 31, pp. 205-216. [DOI: https://dx.doi.org/10.1007/s10545-008-0841-x]

20. Benjamin, R.J.; Piorczynski, T.В.; Hansen, J.M. 285—VPA Inhibits P19 Neural Differentiation through Redox Dysregulation and Oxidative Stress. Free. Radic. Biol. Med.; 2017; 112, 190. [DOI: https://dx.doi.org/10.1016/j.freeradbiomed.2017.10.298]

21. Sepahi, S.; Riahi-Zanjani, B.; Ghorani-Azam, A. Effect of valproic acid on metabolic status and endocrine system in pediatric patients with epilepsy: Systematic literature review. Rev. Clin. Med.; 2017; 4, pp. 7-13. [DOI: https://dx.doi.org/10.22038/rcm.2016.6490]

22. Chateauvieux, S.; Eifes, S.; Morceau, F.; Grigorakaki, C.; Schnekenburger, M.; Henry, E.; Dicato, M.; Diederich, M. Valproic acid perturbs hematopoietic homeostasis by inhibition of erythroid differentiation and activation of the myelo-monocytic pathway. Biochem. Pharmacol.; 2011; 81, pp. 498-509. [DOI: https://dx.doi.org/10.1016/j.bcp.2010.11.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21108933]

23. Tang, W.; Abbott, F.S. Bioactivation of a toxic metabolite of valproic acid, (E)-2-propyl-2,4-pentadienoic acid, via glucuronidation. LC/MS/MS characterization of the GSH-glucuronide diconjugates. Chem. Res. Toxicol.; 1996; 9, pp. 517-526. [DOI: https://dx.doi.org/10.1021/tx950120y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8839057]

24. Kudin, A.P.; Mawasi, H.; Eisenkraft, A.; Elger, C.E.; Bialer, M.; Kunz, W.S. Mitochondrial Liver Toxicity of Valproic Acid and Its Acid Derivatives Is Related to Inhibition of α-Lipoamide Dehydrogenase. Int. J. Mol. Sci.; 2017; 18, 1912. [DOI: https://dx.doi.org/10.3390/ijms18091912] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28878165]

25. Lan, M.J.; Yuan, P.; Chen, G.; Manji, H.K. Neuronal peroxisome proliferator-activated receptor gamma signaling: Regulation by mood-stabilizer valproate. J. Mol. Neurosci.; 2008; 35, pp. 225-234. [DOI: https://dx.doi.org/10.1007/s12031-008-9056-8]

26. Togi, S.; Kamitani, S.; Kawakami, S.; Ikeda, O.; Muromoto, R.; Nanbo, A.; Matsuda, T. HDAC3 influences phosphorylation of STAT3 at serine 727 by interacting with PP2A. Biochem. Biophys. Res. Commun.; 2009; 379, pp. 616-620. [DOI: https://dx.doi.org/10.1016/j.bbrc.2008.12.132]

27. Dodurga, Y.; Gundogdu, G.; Tekin, V.; Koc, T.; Satiroglu-Tufan, N.L.; Bagci, G.; Kucukatay, V. Valproic acid inhibits the proliferation of SHSY5Y neuroblastoma cancer cells by downregulating URG4/URGCP and CCND1 gene expression. Mol. Biol. Rep.; 2014; 41, pp. 4595-5599. [DOI: https://dx.doi.org/10.1007/s11033-014-3330-3]

28. Li, Z.; Yang, L.; Zhang, S.; Song, J.; Sun, H.; Shan, C.; Wang, D.; Liu, S. Valproic acid Suppresses Breast Cancer Cell Growth Through Triggering Pyruvate Kinase M2 Isoform Mediated Warburg Effect. Cell Transplant.; 2021; 30, 9636897211027524. [DOI: https://dx.doi.org/10.1177/09636897211027524]

29. Ebmeyer, J.; Braeuning, A.; Glatt, H.; These, A.; Hessel-Pras, S.; Lampen, A. Human CYP3A4-mediated toxification of the pyrrolizidine alkaloid lasiocarpine. Food Chem. Toxicol.; 2019; 130, pp. 79-88. [DOI: https://dx.doi.org/10.1016/j.fct.2019.05.019]

30. Pal, R.; Singh, K.; Khan, S.A.; Chawla, P.; Kumar, B.; Akhtar, M.J. Reactive metabolites of the anticonvulsant drugs and approaches to minimize the adverse drug reaction. Eur. J. Med. Chem.; 2021; 226, 113890. [DOI: https://dx.doi.org/10.1016/j.ejmech.2021.113890]

31. Surendradoss, J.; Chang, T.K.; Abbott, F.S. Evaluation of in situ generated valproyl 1-O-β-acyl glucuronide in valproic acid toxicity in sandwich-cultured rat hepatocytes. Drug Metab. Dispos.; 2014; 42, pp. 1834-1842. [DOI: https://dx.doi.org/10.1124/dmd.114.059352]

32. Liu, X.; Lv, H.; Guo, Y.; Teka, T.; Wang, X.; Huang, Y.; Han, L.; Pan, G. Structure-Based Reactivity Profiles of Reactive Metabolites with Glutathione. Chem. Res. Toxicol.; 2020; 33, pp. 1579-1593. [DOI: https://dx.doi.org/10.1021/acs.chemrestox.0c00081]

33. Jovanović-Stević, S.; Radisavljević, S.; Scheurer, A.; Ćoćić, D.; Šmit, B.; Petković, M.; Živanović, M.N.; Virijević, K.; Petrović, B. Bis(triazinyl)pyridine complexes of Pt(II) and Pd(II): Studies of the nucleophilic substitution reactions, DNA/HSA interactions, molecular docking and biological activity. J. Biol. Inorg. Chem.; 2021; 26, pp. 625-637. [DOI: https://dx.doi.org/10.1007/s00775-021-01879-3]

34. Motorin, Y.; Helm, M. RNA nucleotide methylation. Wiley Interdiscip. Rev. RNA; 2011; 2, pp. 611-631. [DOI: https://dx.doi.org/10.1002/wrna.79]

35. Leone, A.; Nie, A.; Parker, J.B.; Sawant, S.; Piechta, L.A.; Kelley, M.F.; Kao, L.M.; Jim Proctor, S.; Verheyen, G.; Johnson, M.D. et al. Oxidative stress/reactive metabolite gene expression signature in rat liver detects idiosyncratic hepatotoxicants. Toxicol. Appl. Pharmacol.; 2014; 275, pp. 189-197. [DOI: https://dx.doi.org/10.1016/j.taap.2014.01.017]

36. Iyanagi, T. Molecular mechanism of phase I and phase II drug-metabolizing enzymes: Implications for detoxification. Int. Rev. Cytol.; 2007; 260, pp. 35-112. [DOI: https://dx.doi.org/10.1016/S0074-769660002-8]

37. Shnaider, N.A.; Dmitrenko, D.V. Chronic valproic acid intoxication in epileptology: Diagnosis and treatment. Neurol. Neuropsychiatry Psychosom.; 2016; 8, pp. 94-99. [DOI: https://dx.doi.org/10.14412/2074-2711-2016-2-94-99]

38. Guo, H.L.; Jing, X.; Sun, J.Y.; Hu, Y.H.; Xu, Z.J.; Ni, M.M.; Chen, F.; Lu, X.P.; Qiu, J.C.; Wang, T. Valproic Acid and the Liver Injury in Patients with Epilepsy: An Update. Curr. Pharm. Des.; 2019; 25, pp. 343-351. [DOI: https://dx.doi.org/10.2174/1381612825666190329145428]

39. Aly, M.I.; Abdel-Latif, A.A. Studies on distribution and metabolism of valproate in rat brain, liver, and kidney. Neurochem. Res.; 1980; 5, pp. 1231-1242. [DOI: https://dx.doi.org/10.1007/BF00964959]

40. Gao, Y.; Jiang, D.; Wang, C.; An, G.; Zhu, L.; Cui, C. Comprehensive Analysis of Metabolic Changes in Male Mice Exposed to Sodium Valproate Based on GC-MS Analysis. Drug Des. Devel. Ther.; 2022; 16, pp. 1915-1930. [DOI: https://dx.doi.org/10.2147/DDDT.S357530]

41. Ferguson, C.S.; Tyndale, R.F. Cytochrome P450 enzymes in the brain: Emerging evidence of biological significance. Trends Pharmacol. Sci.; 2011; 32, pp. 708-714. [DOI: https://dx.doi.org/10.1016/j.tips.2011.08.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21975165]

42. Knights, M.; Finlay, E. The effects of sodium valproate on the renal function of children with epilepsy. Pediatr. Nephrol.; 2014; 29, pp. 1131-1138. [DOI: https://dx.doi.org/10.1007/s00467-013-2512-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23720013]

43. Darius, J.; Meyer, F.P.; Sabel, B.A.; Schroeder, U. Influence of nanoparticles on the brain-to-serum distribution and the metabolism of valproic acid in mice. J. Pharm. Pharmacol.; 2000; 52, pp. 1043-1047. [DOI: https://dx.doi.org/10.1211/0022357001774958] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11045883]

44. Mancl, E.E.; Gidal, B.E. The effect of carbapenem antibiotics on plasma concentrations of valproic acid. Ann. Pharmacother.; 2009; 43, pp. 2082-2087. [DOI: https://dx.doi.org/10.1345/aph.1M296] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19934386]

45. Graf, R.; Gossrau, R.; Merker, H.J.; Schwabe, R.; Stahlmann, R.; Nau, H. Enzyme cytochemistry combined with electron microscopy, pharmacokinetics, and clinical chemistry for the evaluation of the effects of steady-state valproic acid concentrations on the mouse. Histochemistry; 1985; 83, pp. 347-358. [DOI: https://dx.doi.org/10.1007/BF00684382] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3934114]

46. Melegh, B.; Pap, M.; Morava, E.; Molnar, D.; Dani, M.; Kurucz, J. Carnitine-dependent changes of metabolic fuel consumption during long-term treatment with valproic acid. J. Pediatr.; 1994; 125, pp. 317-321. [DOI: https://dx.doi.org/10.1016/S0022-3476(94)70218-7]

47. Ma, L.; Wang, Y.; Chen, X.; Zhao, L.; Guo, Y. Involvement of CYP2E1-ROS-CD36/DGAT2 axis in the pathogenesis of VPA-induced hepatic steatosis in vivo and in vitro. Toxicology; 2020; 445, 152585. [DOI: https://dx.doi.org/10.1016/j.tox.2020.152585]

48. Zaccara, G.; Lattanzi, S. A review of pharmacokinetic drug interactions between antimicrobial and antiseizure medications in children. Epileptic Disord.; 2021; 23, pp. 229-256. [DOI: https://dx.doi.org/10.1684/epd.2021.1261]

49. Vanadia, E.; Gibson, K.M.; Pearl, P.L.; Trapolino, E.; Mangano, S.; Vanadia, F. Therapeutic efficacy of magnesium valproate in succinic semialdehyde dehydrogenase deficiency. JIMD Rep.; 2013; 8, pp. 133-137. [DOI: https://dx.doi.org/10.1007/8904_2012_170]

50. Muralidharan, A.; Rahman, J.; Banerjee, D.; Mohammed, A.R.H.; Malik, B.H. Parkinsonism: A Rare Adverse Effect of Valproic Acid. Cureus; 2020; 12, e8782. [DOI: https://dx.doi.org/10.7759/cureus.8782]

51. Lunke, S.; Maxwell, S.; Khurana, I.; KN, H.; Okabe, J.; Al-Hasani, K.; El-Osta, A. Epigenetic evidence of an Ac/Dc axis by VPA and SAHA. Clin. Epigenetics; 2021; 13, 58. [DOI: https://dx.doi.org/10.1186/s13148-021-01050-4]

52. Sakakibara, Y.; Kojima, A.; Asai, Y.; Nadai, M.; Katoh, M. Changes in uridine 5’-diphospho-glucuronosyltransferase 1A6 expression by histone deacetylase inhibitor valproic acid. Biopharm. Drug Dispos.; 2022; 43, pp. 175-182. [DOI: https://dx.doi.org/10.1002/bdd.2328]

53. Yamada, M.K.; Nakanishi, K.; Ohba, S.; Nakamura, T.; Ikegaya, Y.; Nishiyama, N.; Matsuki, N. Brain-derived neurotrophic factor promotes the maturation of GABAergic mechanisms in cultured hippocampal neurons. J. Neurosci.; 2002; 22, pp. 7580-7585. [DOI: https://dx.doi.org/10.1523/JNEUROSCI.22-17-07580.2002]

54. Shaltiel, G.; Mark, S.; Kofman, O.; Belmaker, R.H.; Agam, G. Effect of valproate derivatives on human brain myo-inositol-1-phosphate (MIP) synthase activity and amphetamine-induced rearing. Pharmacol. Rep.; 2007; 59, pp. 402-407.

55. Watterson, J.M.; Watson, D.G.; Meyer, E.M.; Lenox, R.H. A role for protein kinase C and its substrates in the action of valproic acid in the brain: Implications for neural plasticity. Brain Res.; 2002; 934, pp. 69-80. [DOI: https://dx.doi.org/10.1016/S0006-8993(02)02362-4]

56. Parikh, S.K.; Silberstein, S.D. Current Status of Antiepileptic Drugs as Preventive Migraine Therapy. Curr. Treat. Options Neurol.; 2019; 21, 16. [DOI: https://dx.doi.org/10.1007/s11940-019-0558-1]

57. Traetta, M.E.; Codagnone, M.G.; Uccelli, N.A.; Ramos, A.J.; Zárate, S.; Reinés, A. Hippocampal neurons isolated from rats subjected to the valproic acid model mimic in vivo synaptic pattern: Evidence of neuronal priming during early development in autism spectrum disorders. Mol. Autism.; 2021; 12, 23. [DOI: https://dx.doi.org/10.1186/s13229-021-00428-8]

58. Rosenberg, G. The mechanisms of action of valproate in neuropsychiatric disorders: Can we see the forest for the trees?. Cell Mol. Life Sci.; 2007; 64, pp. 2090-2103. [DOI: https://dx.doi.org/10.1007/s00018-007-7079-x]

59. Becker, W.J. Acute Migraine Treatment in Adults. Headache; 2015; 55, pp. 778-793. [DOI: https://dx.doi.org/10.1111/head.12550]

60. Sharma, R.P.; Rosen, C.; Kartan, S.; Guidotti, A.; Costa, E.; Grayson, D.R.; Chase, K. Valproic acid and chromatin remodeling in schizophrenia and bipolar disorder: Preliminary results from a clinical population. Schizophr. Res.; 2006; 88, pp. 227-231. [DOI: https://dx.doi.org/10.1016/j.schres.2006.07.015]

61. Bahna, S.G.; Niles, L.P. Epigenetic regulation of melatonin receptors in neuropsychiatric disorders. Br. J. Pharmacol.; 2018; 175, pp. 3209-3219. [DOI: https://dx.doi.org/10.1111/bph.14058] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28967098]

62. Yoon, H.Y.; Ahn, M.H.; Yee, J.; Lee, N.; Han, J.M.; Gwak, H.S. Influence of CYP2C9 and CYP2A6 on plasma concentrations of valproic acid: A meta-analysis. Eur. J. Clin. Pharmacol.; 2020; 76, pp. 1053-1058. [DOI: https://dx.doi.org/10.1007/s00228-020-02872-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32385545]

63. Nasyrova, R.F.; Sivakova, N.А.; Lipatova, L.V.; Ivashchenko, D.V.; Sosina, K.А.; Drokov, А.P.; Shnayder, N.А. Biological Markers Of The Antiepileptic Drugs Efficacy And Safety: Pharmacogenetics And Pharmacokinetics. Sib. Med. Rev.; 2017; 1, pp. 18-25. [DOI: https://dx.doi.org/10.20333/2500136-2017-1-17-25]

64. Kiang, T.K.; Ho, P.C.; Anari, M.R.; Tong, V.; Abbott, F.S.; Chang, T.K. Contribution of CYP2C9, CYP2A6, and CYP2B6 to valproic acid metabolism in hepatic microsomes from individuals with the CYP2C9*1/*1 genotype. Toxicol. Sci.; 2006; 94, pp. 261-271. [DOI: https://dx.doi.org/10.1093/toxsci/kfl096] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16945988]

65. Tan, L.; Yu, J.T.; Sun, Y.P.; Ou, J.R.; Song, J.H.; Yu, Y. The influence of cytochrome oxidase CYP2A6, CYP2B6, and CYP2C9 polymorphisms on the plasma concentrations of valproic acid in epileptic patients. Clin. Neurol. Neurosurg.; 2010; 112, pp. 320-323. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20089352][DOI: https://dx.doi.org/10.1016/j.clineuro.2010.01.002]

66. Song, C.; Li, X.; Mao, P.; Song, W.; Liu, L.; Zhang, Y. Impact of CYP2C19 and CYP2C9 gene polymorphisms on sodium valproate plasma concentration in patients with epilepsy. Eur. J. Hosp. Pharm.; 2022; 29, pp. 198-201. [DOI: https://dx.doi.org/10.1136/ejhpharm-2020-002367]

67. Chen, Y.B.; Liang, C.S.; Wang, L.J.; Hung, K.C.; Carvalho, A.F.; Solmi, M.; Vieta, E.; Tseng, P.T.; Lin, P.Y.; Tu, Y.K. et al. Comparative effectiveness of valproic acid in different serum concentrations for maintenance treatment of bipolar disorder: A retrospective cohort study using target trial emulation framework. EClinicalMedicine; 2022; 54, 101678. [DOI: https://dx.doi.org/10.1016/j.eclinm.2022.101678]

68. Zhao, M.; Zhang, T.; Li, G.; Qiu, F.; Sun, Y.; Zhao, L. Associations of CYP2C9 and CYP2A6 Polymorphisms with the Concentrations of Valproate and its Hepatotoxin Metabolites and Valproate-Induced Hepatotoxicity. Basic Clin. Pharmacol. Toxicol.; 2017; 121, pp. 138-143. [DOI: https://dx.doi.org/10.1111/bcpt.12776]

69. Shnayder, N.А.; Dmitrenko, D.V.; Bochanova, E.N. Pharmacogenetics. valproic acid. Clinical Psychopharmacogenetics: Monograph. Ad.; Nasyrova, R.F.; Nesnanov, N.G. Publisher House DEAN: St. Petersburg, Russia, 2019; pp. 265-278.

70. Wang, C.; Wang, P.; Yang, L.P.; Pan, J.; Yang, X.; Ma, H.Y. Association of CYP2C9, CYP2A6, ACSM2A, and CPT1A gene polymorphisms with adverse effects of valproic acid in Chinese patients with epilepsy. Epilepsy Res.; 2017; 132, pp. 64-69. [DOI: https://dx.doi.org/10.1016/j.eplepsyres.2017.02.015]

71. Guo, Y.; Hu, C.; He, X.; Qiu, F.; Zhao, L. Effects of UGT1A6, UGT2B7, and CYP2C9 genotypes on plasma concentrations of valproic acid in Chinese children with epilepsy. Drug Metab. Pharmacokinet.; 2012; 27, pp. 536-542. [DOI: https://dx.doi.org/10.2133/dmpk.DMPK-11-NT-144]

72. Hung, C.C.; Ho, J.L.; Chang, W.L.; Tai, J.J.; Hsieh, T.J.; Hsieh, Y.W.; Hsieh, Y.W.; Liou, H.H. Association of genetic variants in six candidate genes with valproic acid therapy optimization. Pharmacogenomics; 2011; 12, pp. 1107-1117. [DOI: https://dx.doi.org/10.2217/pgs.11.64]

73. Marahatta, A.; Bhandary, B.; Jeong, S.K.; Kim, H.R.; Chae, H.J. Soybean greatly reduces valproic acid plasma concentrations: A food-drug interaction study. Sci. Rep.; 2014; 4, 4362. [DOI: https://dx.doi.org/10.1038/srep04362]

74. Chu, X.M.; Zhang, L.F.; Wang, G.J.; Zhang, S.N.; Zhou, J.H.; Hao, H.P. Influence of UDP-glucuronosyltransferase polymorphisms on valproic acid pharmacokinetics in Chinese epilepsy patients. Eur. J. Clin. Pharmacol.; 2012; 68, pp. 1395-1401. [DOI: https://dx.doi.org/10.1007/s00228-012-1277-7]

75. Stępień, K.M.; Tomaszewski, M.; Tomaszewska, J.; Czuczwar, S.J. The multidrug transporter P-glycoprotein in pharmacoresistance to antiepileptic drugs. Pharmacol. Rep.; 2012; 64, pp. 1011-1019. [DOI: https://dx.doi.org/10.1016/S1734-1140(12)70900-3]

76. Hung, C.C.; Chen, C.C.; Lin, C.J.; Liou, H.H. Functional evaluation of polymorphisms in the human ABCB1 gene and the impact on clinical responses of antiepileptic drugs. Pharmacogenet. Genomics.; 2008; 18, pp. 390-402. [DOI: https://dx.doi.org/10.1097/FPC.0b013e3282f85e36]

77. Löscher, W.; Potschka, H. Blood-brain barrier active efflux transporters: ATP-binding cassette gene family. J. Am. Soc. Exp. Neurother.; 2005; 2, pp. 86-98. [DOI: https://dx.doi.org/10.1602/neurorx.2.1.86]

78. Sabin, O.; Pop, R.; Trifa, A.; Buzoianu, A.D. The influence of CYP2C9, CYP2C19 and ABCB1 polymorphisms on the plasma concentrations of valproic acid in epileptic patients. HVM Bioflux.; 2016; 8, pp. 29-33.

79. Siddiqui, A.; Kerb, R.; Weale, M.E.; Brinkmann, U.; Smith, A.; Goldstein, D.B.; Wood, N.W.; Sisodiya, S.M. Association of multidrug resistance in epilepsy with a polymorphism in the drug-transporter gene ABCB1. N. Engl. J. Med.; 2003; 348, pp. 1442-1448. [DOI: https://dx.doi.org/10.1056/NEJMoa021986]

80. Ajmi, M.; Boujaafar, S.; Zouari, N.; Amor, D.; Nasr, A.; Rejeb, N.B.; Amor, S.B.; Omezzine, A.; Benammou, S.; Bouslama, A. Association between ABCB1 polymorphisms and response to first-generation antiepileptic drugs in a Tunisian epileptic population. Int. J. Neurosci.; 2018; 128, pp. 705-714. [DOI: https://dx.doi.org/10.1080/00207454.2017.1412964]

81. Yu, L.; Liao, W.P.; Yi, Y.H.; Qiu, G. ABCB1 G2677T/A polymorphism is associated with the risk of drug-resistant epilepsy in Asians. Epilepsy Res.; 2015; 115, pp. 100-108. [DOI: https://dx.doi.org/10.1016/j.eplepsyres.2015.05.015]

82. Yagi, M.; Nakamura, T.; Okizuka, Y.; Oyazato, Y.; Kawasaki, Y.; Tsuneishi, S.; Sakaeda, T.; Matsuo, M.; Okumura, K.; Okamura, N. Effect of CPS14217C>A genotype on valproic-acid-induced hyperammonemia. Pediatr. Int.; 2010; 52, pp. 744-748. [DOI: https://dx.doi.org/10.1111/j.1442-200X.2010.03157.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20456087]

83. Chen, L.; Tian, Q.; Zhang, M.; Chen, D.; Gao, X.; Yang, H.; Li, H.; Li, C.; Wen, J.; Li, Y. et al. CPS1 T1405N polymorphism, HDL cholesterol, homocysteine and renal function are risk factors of VPA induced hyperammonemia among epilepsy patients. Epilepsy Res.; 2019; 154, pp. 139-143. [DOI: https://dx.doi.org/10.1016/j.eplepsyres.2019.05.010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31151073]

84. Hudson, G.; Chinnery, P.F. Mitochondrial DNA polymerase-gamma and human disease. Hum. Mol. Genet.; 2006; 15, pp. R244-R252. [DOI: https://dx.doi.org/10.1093/hmg/ddl233] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16987890]

85. Shilkina, O.S.; Shnayder, N.A.; Artyukhovm, I.P.; Moskaleva, P.V.; Panina, Y.S. Problems of differential diagnosis of myoclon us-epilepsy associated with the mutation of the POLG gene and juvenile myoclonic epilepsy: A clinical case. Russ. J. Child Neurol.; 2018; 13, pp. 57-63. [DOI: https://dx.doi.org/10.17650/2073-8803-2018-13-1-57-63]

86. Stewart, J.D.; Horvath, R.; Baruffini, E.; Ferrero, I.; Bulst, S.; Watkins, P.B.; Fontana, R.J.; Day, C.P.; Chinnery, P.F. Polymerase γ gene POLG determines the risk of sodium valproate-induced liver toxicity. Hepatology; 2010; 52, pp. 1791-1796. [DOI: https://dx.doi.org/10.1002/hep.23891]

87. Yang, K.F.; Meng, L.Y.; Zhang, Y.F.; Yang, Y.C.; Cheng, H.Y.; Jiang, Z.H.; Zhang, H.; Chen, Y.C. POLG Mutations Are Probably Rare in the Han Chinese Population. Chin. Med. Sci. J.; 2020; 35, pp. 350-356. [DOI: https://dx.doi.org/10.24920/003697]

88. Saruwatari, J.; Deguchi, M.; Yoshimori, Y.; Noai, M.; Yoshida, S.; Ogusu, N.; Oniki, K.; Yoshida, S.; Yasui-Furukori, N.; Kaneko, S. et al. Superoxide dismutase 2 Val16Ala polymorphism is a risk factor for the valproic acid-related elevation of serum aminotransferases. Epilepsy Res.; 2012; 99, pp. 183-186. [DOI: https://dx.doi.org/10.1016/j.eplepsyres.2011.10.033]

89. Hsieh, C.L.; Chen, K.C.; Ding, C.Y.; Tsai, W.J.; Wu, J.F.; Peng, C.C. Valproic acid substantially downregulated genes folr1, IGF2R, RGS2, COL6A3, EDNRB, KLF6, and pax-3, N-acetylcysteine alleviated most of the induced gene alterations in chicken embryo model. Rom. J. Morphol. Embryol.; 2013; 54, pp. 993-1004.

90. Kim, B.; Kim, C.Y.; Lee, M.J.; Joo, Y.H. Preliminary evidence on the association between XBP1-116C/G polymorphism and response to prophylactic treatment with valproate in bipolar disorders. Psychiatry Res.; 2009; 168, pp. 209-212. [DOI: https://dx.doi.org/10.1016/j.psychres.2008.05.010]

91. Li, H.; Wang, X.; Zhou, Y.; Ni, G.; Su, Q.; Chen, Z.; Chen, Z.; Li, J.; Chen, X.; Hou, X. et al. Association of LEPR and ANKK1 Gene Polymorphisms with Weight Gain in Epilepsy Patients Receiving Valproic Acid. Int. J. Neuropsychopharmacol.; 2015; 18, pyv021. [DOI: https://dx.doi.org/10.1093/ijnp/pyv021]

92. Drokov, A.P.; Lipatova, L.V.; Shnaĭder, N.A.; Nasyrova, R.F. Pharmacogenetic markers of metabolic disorders in the treatment with valproic acid. Zhurnal Nevrol. Psikhiatrii Im. S.S. Korsakova; 2018; 118, pp. 82-89. [DOI: https://dx.doi.org/10.17116/jnevro201811810282]

93. Song, A.; Cho, G.-W.; Vijayakumar, K.A.; Moon, C.; Ang, M.J.; Kim, J.; Park, I.; Jang, C.H. Neuroprotective Effect of Valproic Acid on Salicylate-Induced Tinnitus. Int. J. Mol. Sci.; 2021; 23, 23. [DOI: https://dx.doi.org/10.3390/ijms23010023]

94. Siemes, H.; Nau, H.; Schultze, K.; Wittfoht, W.; Drews, E.; Penzien, J.; Seidel, U. Valproate (VPA) metabolites in various clinical conditions of probable VPA-associated hepatotoxicity. Epilepsia; 1993; 34, pp. 332-346. [DOI: https://dx.doi.org/10.1111/j.1528-1157.1993.tb02419.x]

95. Ghodke-Puranik, Y.; Thorn, C.F.; Lamba, J.K.; Leeder, J.S.; Song, W.; Birnbaum, A.K.; Altman, R.B.; Klein, T.E. Valproic acid pathway: Pharmacokinetics and pharmacodynamics. Pharm. Genom.; 2013; 23, pp. 236-241. [DOI: https://dx.doi.org/10.1097/FPC.0b013e32835ea0b2]

96. Beger, R.D.; Schmidt, M.A.; Kaddurah-Daouk, R. Current Concepts in Pharmacometabolomics, Biomarker Discovery, and Precision Medicine. Metabolites; 2020; 10, 129. [DOI: https://dx.doi.org/10.3390/metabo10040129]

97. Wishart, D.S.; Guo, A.; Oler, E.; Wang, F.; Anjum, A.; Peters, H.; Dizon, R.; Sayeeda, Z.; Tian, S.; Lee, B.L. et al. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res.; 2022; 50, pp. D622-D631. [DOI: https://dx.doi.org/10.1093/nar/gkab1062]

98. Suzuki, E.; Yamamura, N.; Ogura, Y.; Nakai, D.; Kubota, K.; Kobayashi, N.; Miura, S.; Okazaki, O. Identification of valproic acid glucuronide hydrolase as a key enzyme for the interaction of valproic acid with carbapenem antibiotics. Drug Metab. Dispos.; 2010; 38, pp. 1538-1544. [DOI: https://dx.doi.org/10.1124/dmd.110.032938]

99. Silva, M.F.; Ruiter, J.P.; IJlst, L.; Jakobs, C.; Duran, M.; de Almeida, I.T.; Wanders, R.J. Differential effect of valproate and its Delta2- and Delta4-unsaturated metabolites, on the beta-oxidation rate of long-chain and medium-chain fatty acids. Chem. Biol. Interact.; 2001; 137, pp. 203-212. [DOI: https://dx.doi.org/10.1016/S0009-2797(01)00234-4]

100. Friel, P. Valproyl CoA: An active metabolite of valproate?. Med. Hypotheses; 1990; 31, pp. 31-32. [DOI: https://dx.doi.org/10.1016/0306-9877(90)90049-K]

101. Eyal, S.; Yagen, B.; Shimshoni, J.; Bialer, M. Histone deacetylases inhibition and tumor cells cytotoxicity by CNS-active VPA constitutional isomers and derivatives. Biochem. Pharmacol.; 2005; 69, pp. 1501-1508. [DOI: https://dx.doi.org/10.1016/j.bcp.2005.02.012]

102. Ehlers, K.; Stürje, H.; Merker, H.J.; Nau, H. The valproic acid metabolite E-2-n-propyl-2-pentenoic acid does not induce spina bifida in the mouse. Dev. Pharmacol. Ther.; 1992; 19, pp. 196-204. [DOI: https://dx.doi.org/10.1159/000457485] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/1343622]

103. Kim, T.; Song, S.; Park, Y.; Kang, S.; Seo, H. HDAC Inhibition by Valproic Acid Induces Neuroprotection and Improvement of PD-like Behaviors in LRRK2 R1441G Transgenic Mice. Exp. Neurobiol.; 2019; 28, pp. 504-515. [DOI: https://dx.doi.org/10.5607/en.2019.28.4.504] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31495079]

104. Zádori, D.; Geisz, A.; Vámos, E.; Vécsei, L.; Klivényi, P. Valproate ameliorates the survival and the motor performance in a transgenic mouse model of Huntington’s disease. Pharmacol. Biochem. Behav.; 2009; 94, pp. 148-153. [DOI: https://dx.doi.org/10.1016/j.pbb.2009.08.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19698736]

105. Binvignat, O.; Olloquequi, J. Excitotoxicity as a Target Against Neurodegenerative Processes. Curr. Pharm. Des.; 2020; 26, pp. 1251-1262. [DOI: https://dx.doi.org/10.2174/1381612826666200113162641]

106. Scheuing, L.; Chiu, C.T.; Liao, H.M.; Linares, G.R.; Chuang, D.M. Preclinical and clinical investigations of mood stabilizers for Huntington’s disease: What have we learned?. Int. J. Biol. Sci.; 2014; 10, pp. 1024-1038. [DOI: https://dx.doi.org/10.7150/ijbs.9898] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25285035]

107. Soria-Castro, R.; Schcolnik-Cabrera, A.; Rodríguez-López, G.; Campillo-Navarro, M.; Puebla-Osorio, N.; Estrada-Parra, S.; Estrada-García, I.; Chacón-Salinas, R.; Chávez-Blanco, A.D. Exploring the Drug Repurposing Versatility of Valproic Acid as a Multifunctional Regulator of Innate and Adaptive Immune Cells. J. Immunol. Res.; 2019; 2019, 9678098. [DOI: https://dx.doi.org/10.1155/2019/9678098] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31001564]

108. Reis, K.P.; Sperling, L.E.; Teixeira, C.; Sommer, L.; Colombo, M.; Koester, L.S.; Pranke, P. VPA/PLGA microfibers produced by coaxial electrospinning for the treatment of central nervous system injury. Braz. J. Med. Biol. Res.; 2020; 53, e8993. [DOI: https://dx.doi.org/10.1590/1414-431x20208993] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32294700]

109. Chen, S.; Ye, J.; Chen, X.; Shi, J.; Wu, W.; Lin, W.; Lin, W.; Li, Y.; Fu, H.; Li, S. Valproic acid attenuates traumatic spinal cord injury-induced inflammation via STAT1 and NF-κB pathway dependent of HDAC3. J. Neuroinflammation.; 2018; 15, 150. [DOI: https://dx.doi.org/10.1186/s12974-018-1193-6]

110. Heers, H.; Stanislaw, J.; Harrelson, J.; Lee, M.W. Valproic acid as an adjunctive therapeutic agent for the treatment of breast cancer. Eur. J. Pharmacol.; 2018; 835, pp. 61-74. [DOI: https://dx.doi.org/10.1016/j.ejphar.2018.07.057]

111. Scholz, B.; Schulte, J.S.; Hamer, S.; Himmler, K.; Pluteanu, F.; Seidl, M.D.; Stein, J.; Wardelmann, E.; Hammer, E.; Völker, U. et al. HDAC (Histone Deacetylase) Inhibitor Valproic Acid Attenuates Atrial Remodeling and Delays the Onset of Atrial Fibrillation in Mice. Circ. Arrhythm. Electrophysiol.; 2019; 12, e007071. [DOI: https://dx.doi.org/10.1161/CIRCEP.118.007071]

112. Abbott, F.S.; Kassam, J.; Orr, J.M.; Farrell, K. The effect of aspirin on valproic acid metabolism. Clin. Pharmacol. Ther.; 1986; 40, pp. 94-100. [DOI: https://dx.doi.org/10.1038/clpt.1986.144]

113. Ohnishi, S.; Okamura, N.; Sakamoto, S.; Hasegawa, H.; Norikura, R.; Kanaoka, E.; Takahashi, K.; Horie, K.; Sakamoto, K.; Baba, T. Role of Na+/L-carnitine transporter (OCTN2) in renal handling of pivaloylcarnitine and valproylcarnitine formed during pivalic acid-containing prodrugs and valproic acid treatment. Drug Metab. Pharmacokinet.; 2008; 23, pp. 293-303. [DOI: https://dx.doi.org/10.2133/dmpk.23.293]

114. Sztajnkrycer, M.D. Valproic acid toxicity: Overview and management. J. Toxicol. Clin. Toxicol.; 2002; 40, pp. 789-801. [DOI: https://dx.doi.org/10.1081/CLT-120014645]

115. Luís, P.B.; Ruiter, J.; IJlst, L.; de Almeida, I.T.; Duran, M.; Wanders, R.J.; Silva, M.F. Valproyl-CoA inhibits the activity of ATP- and GTP-dependent succinate:CoA ligases. J. Inherit. Metab. Dis.; 2014; 37, pp. 353-357. [DOI: https://dx.doi.org/10.1007/s10545-013-9657-4]

116. PubChem. Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information; 2004-. PubChem Compound Summary for CID 125849, Valproylcarnitine. Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Valproylcarnitine (accessed on 21 December 2022).

117. Lheureux, P.E.; Hantson, P. Carnitine in the treatment of valproic acid-induced toxicity. Clin. Toxicol.; 2009; 47, pp. 101-111. [DOI: https://dx.doi.org/10.1080/15563650902752376]

118. Millington, D.S.; Bohan, T.P.; Roe, C.R.; Yergey, A.L.; Liberato, D.J. Valproylcarnitine: A novel drug metabolite identified by fast atom bombardment and thermospray liquid chromatography-mass spectrometry. Clin. Chim. Acta.; 1985; 145, pp. 69-76. [DOI: https://dx.doi.org/10.1016/0009-8981(85)90020-8]

119. Torres, S.; Baulies, A.; Insausti-Urkia, N.; Alarcón-Vila, C.; Fucho, R.; Solsona-Vilarrasa, E.; Núñez, S.; Robles, D.; Ribas, V.; Wakefield, L. et al. Endoplasmic Reticulum Stress-Induced Upregulation of STARD1 Promotes Acetaminophen-Induced Acute Liver Failure. Gastroenterology; 2019; 157, pp. 552-568. [DOI: https://dx.doi.org/10.1053/j.gastro.2019.04.023]

120. Dmitrenko, D.V.; Shnaider, N.A.; Strotskaya, I.G.; Kichkaylo, A.S.; Zobova, S.N. Mechanisms of valproate-induced teratogenesis. Neurol. Neuropsychiatry Psychosom.; 2017; 9, pp. 89-96. [DOI: https://dx.doi.org/10.14412/2074-2711-2017-1S-89-96]

121. Doré, M.; San Juan, A.E.; Frenette, A.J.; Williamson, D. Clinical Importance of Monitoring Unbound Valproic Acid Concentration in Patients with Hypoalbuminemia. Pharmacotherapy; 2017; 37, pp. 900-907. [DOI: https://dx.doi.org/10.1002/phar.1965]

122. Patel, A.R.; Nagalli, S. Valproate Toxicity; StatPearls Publishing: Treasure Island, FL, USA, 2022; Available online: https://www.ncbi.nlm.nih.gov/books/NBK560898/ (accessed on 27 November 2022).

123. Becker, C.M.; Harris, R.A. Influence of valproic acid on hepatic carbohydrate and lipid metabolism. Arch. Biochem. Biophys.; 1983; 223, pp. 381-392. [DOI: https://dx.doi.org/10.1016/0003-9861(83)90602-1]

124. Li, X.; Zhu, Y.; He, H.; Lou, L.; Ye, W.; Chen, Y.; Wang, J. Synergistically killing activity of aspirin and histone deacetylase inhibitor valproic acid (VPA) on hepatocellular cancer cells. Biochem. Biophys. Res. Commun.; 2013; 436, pp. 259-264. [DOI: https://dx.doi.org/10.1016/j.bbrc.2013.05.088] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23726914]

125. Vaiman, E.E.; Shnaider, N.A.; Nesnanov, N.G.; Nasyrova, R.F. Drug-induced Parkinsonism. Soc. Clin. Psychiatry; 2021; 13, pp. 96-103.

126. Carriere, C.H.; Kang, N.H.; Niles, L.P. Neuroprotection by valproic acid in an intrastriatal rotenone model of Parkinson’s disease. Neuroscience; 2014; 267, pp. 114-121. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2014.02.028] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24613722]

127. Mahmoud, F.; Tampi, R.R. Valproic acid-induced parkinsonism in the elderly: A comprehensive review of the literature. Am. J. Geriatr. Pharmacother.; 2011; 9, pp. 405-412. [DOI: https://dx.doi.org/10.1016/j.amjopharm.2011.09.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21993183]

128. Neavin, D.; Kaddurah-Daouk, R.; Weinshilboum, R. Pharmacometabolomics informs Pharmacogenomics. Metabolomics; 2016; 12, 121. [DOI: https://dx.doi.org/10.1007/s11306-016-1066-x]

129. Ellero-Simatos, S.; Lewis, J.P.; Georgiades, A.; Yerges-Armstrong, L.M.; Beitelshees, A.L.; Horenstein, R.B.; Dane, A.; Harms, A.C.; Ramaker, R.; Vreeken, R.J. et al. Pharmacometabolomics reveals that serotonin is implicated in aspirin response variability. CPT Pharmacomet. Syst. Pharmacol.; 2014; 3, e125. [DOI: https://dx.doi.org/10.1038/psp.2014.22]

130. Wang, L.; McLeod, H.L.; Weinshilboum, R.M. Genomics and drug response. N. Engl. J. Med.; 2011; 364, pp. 1144-1153. [DOI: https://dx.doi.org/10.1056/NEJMra1010600]

131. Nasyrova, R.F.; Neznanov, N.G. Clinical Psychopharmacogenetics; Publisher DEAN: St. Petersburg, Russia, 2019; 405.

132. Kaddurah-Daouk, R.; Weinshilboum, R.M. Pharmacometabolomics Research Network. Pharmacometabolomics: Implications for clinical pharmacology and systems pharmacology. Clin. Pharmacol. Ther.; 2014; 95, pp. 154-167. [DOI: https://dx.doi.org/10.1038/clpt.2013.217]

133. Okumura, A.; Takagi, M.; Numoto, S.; Iwayama, H.; Azuma, Y.; Kurahashi, H. Effects of l-carnitine supplementation in patients with childhood-onset epilepsy prescribed valproate. Epilepsy Behav.; 2021; 122, 108220. [DOI: https://dx.doi.org/10.1016/j.yebeh.2021.108220]

134. Poupon, R. Ursodeoxycholic acid and bile-acid mimetics as therapeutic agents for cholestatic liver diseases: An overview of their mechanisms of action. Clin. Res. Hepatol. Gastroenterol.; 2012; 36, (Suppl. 1), pp. S3-S12. [DOI: https://dx.doi.org/10.1016/S2210-7401(12)70015-3]

135. Goossens, J.F.; Bailly, C. Ursodeoxycholic acid and cancer: From chemoprevention to chemotherapy. Pharmacol. Ther.; 2019; 203, 107396. [DOI: https://dx.doi.org/10.1016/j.pharmthera.2019.107396]

136. Yoon, S.; Lee, H.; Ji, S.C.; Yoon, S.H.; Cho, J.Y.; Chung, J.Y. Pharmacokinetics and Pharmacodynamics of Ursodeoxycholic Acid in an Overweight Population with Abnormal Liver Function. Clin. Pharmacol. Drug Dev.; 2021; 10, pp. 68-77. [DOI: https://dx.doi.org/10.1002/cpdd.790]

137. Asgarshirazi, M.; Shariat, M.; Dalili, H.; Keihanidoost, Z. Ursodeoxycholic Acid Can Improve Liver Transaminase Quantities in Children with Anticonvulsant Drugs Hepatotoxicity: A Pilot Study. Acta Med. Iran.; 2015; 53, pp. 351-355.

138. Plummer, S.; Beaumont, B.; Elcombe, M.; Wallace, S.; Wright, J.; Mcinnes, E.F.; Currie, R.A.; Cowie, D. Species differences in phenobarbital-mediated UGT gene induction in rat and human liver microtissues. Toxicol. Rep.; 2020; 8, pp. 155-161. [DOI: https://dx.doi.org/10.1016/j.toxrep.2020.12.019]

139. Tutty, M.A.; Movia, D.; Prina-Mello, A. Three-dimensional (3D) liver cell models—A tool for bridging the gap between animal studies and clinical trials when screening liver accumulation and toxicity of nanobiomaterials. Drug Deliv. Transl. Res.; 2022; 12, pp. 2048-2074. [DOI: https://dx.doi.org/10.1007/s13346-022-01147-0]

140. Ravi, N.V.; Maany, I.; Burke, W.M.; Dhopesh, V.; Woody, G.E. Detoxification with phenobarbital of alprazolam-dependent polysubstance abusers. J. Subst. Abuse Treat.; 1990; 7, pp. 55-58. [DOI: https://dx.doi.org/10.1016/0740-5472(90)90036-P]

141. Bochanova, E.N.; Snaider, N.A.; Dmitrenco, D.V.; Artukhov, I.P.; Gusev, S.D.; Yurjieva, E.A.; Shilkina, O.S. Process of Personalized Prescription of Valproic Acid as the Main Element of the Management of Epilepsy. Int. J. Biomed.; 2018; 8, pp. 26-32. [DOI: https://dx.doi.org/10.21103/Article8(1)_OA3]

142. Ji, Y.; Hebbring, S.; Zhu, H.; Jenkins, G.D.; Biernacka, J.; Snyder, K.; Drews, M.; Fiehn, O.; Zeng, Z.; Schaid, D. et al. Glycine and a glycine dehydrogenase (GLDC) SNP as citalopram/escitalopram response biomarkers in depression: Pharmacometabolomics-informed pharmacogenomics. Clin. Pharmacol. Ther.; 2011; 89, pp. 97-104. [DOI: https://dx.doi.org/10.1038/clpt.2010.250]

143. Sychev, D.А.; Shuev, G.N.; Torbenkov, Е.S.; Adriyanova, М.А. Personalized medicine: The view of a clinical pharmacologist. Cons. Med.; 2017; 19, pp. 61-68. [DOI: https://dx.doi.org/10.26442/2075-1753_2017.19.1.61-68]