Isoliquiritigenin, a potent human monoamine oxidase inhibitor, modulates dopamine D1, D3, and vasopressin V1A receptors

Isoliquiritigenin (= 4,2′,4′-Trihydroxychalcone) (ILG) is a major constituent of the Glycyrrhizae Rhizoma that has significant neuroprotective functions. In the present study, we re-examined the potential of ILG to inhibit human monoamine oxidase (hMAO) in vitro and established its mechanism of inhibition through a kinetics study and molecular docking examination. ILG showed competitive inhibition of hMAO-A and mixed inhibition of hMAO-B with IC50 values of 0.68 and 0.33 µM, respectively, which varied slightly from the reported IC50 values. Since ILG has been reported to reduce dopaminergic neurodegeneration and psychostimulant-induced toxicity (both of which are related to dopamine and vasopressin receptors), we investigated the binding affinity and modulatory functions of ILG on dopamine and vasopressin receptors. ILG was explored as an antagonist of the D1 receptor and an agonist of the D3 and V1A receptors with good potency. An in silico docking investigation revealed that ILG can interact with active site residues at target receptors with low binding energies. These activities of ILG on hMAO and brain receptors suggest the potential role of the compound to ameliorate dopaminergic deficits, depression, anxiety, and associated symptoms in Parkinson’s disease and other neuronal disorders.


Isoliquiritigenin (= 4,2′,4′-Trihydroxychalcone) (ILG) is a major constituent of the Glycyrrhizae Rhizoma that has significant neuroprotective functions. In the present study, we re-examined the potential of ILG to inhibit human monoamine oxidase (hMAO) in vitro and established its mechanism of inhibition through a kinetics study and molecular docking examination. ILG showed competitive inhibition of hMAO-A and mixed inhibition of hMAO-B with IC
values of 0.68 and 0.33 µM, respectively, which varied slightly from the reported IC 50 values. Since ILG has been reported to reduce dopaminergic neurodegeneration and psychostimulant-induced toxicity (both of which are related to dopamine and vasopressin receptors), we investigated the binding affinity and modulatory functions of ILG on dopamine and vasopressin receptors. ILG was explored as an antagonist of the D 1 receptor and an agonist of the D 3 and V 1A receptors with good potency. An in silico docking investigation revealed that ILG can interact with active site residues at target receptors with low binding energies. These activities of ILG on hMAO and brain receptors suggest the potential role of the compound to ameliorate dopaminergic deficits, depression, anxiety, and associated symptoms in Parkinson's disease and other neuronal disorders.
Neurodegenerative diseases (NDDs) include a heterogeneous set of disorders characterized by progressive and irreversible damage to the neurons of the central nervous system (CNS) that can lead to functional and mental impairments 1 . Alzheimer's disease (AD) and Parkinson's disease (PD) are the common neurodegenerative disorders; each has a prevalence of approximately 1% of the elderly population in their late sixties and even more in the latter decades of life 2 . PD is a nervous system disorder characterized by progressive degeneration of nigrostriatal dopaminergic neurons in the midbrain, resulting in depletion of striatal dopamine (DA). Motor symptoms (such as tremor at rest, rigidity, bradykinesia, and postural instability) are the hallmarks of PD and often are accompanied by non-motor features, like cognitive dysfunction, depression, anxiety, and other socio-behavioral anomalies. Multiple mechanisms that provoke neurodegeneration have been proposed, such as complex molecular changes involving DA and non-dopamine neurons (such as cholinergic and gamma-aminobutyric acid-B [GABA]ergic neurons) and non-neuronal cells (such as microglia and astrocytes), perturbed cell functions like mitochondrial dysfunction resulting from excessive reactive oxygen species (ROS) or either underexpression or mutation of protective genes like peroxisome proliferator-activated receptor-gamma activator-1α (PGC-1α) and parkin-associated proteins (PARK6, PARK7, and PARK8). In addition, α-synuclein aggregation is a pathological marker of PD. Overall, the pathogenesis of PD represents a complex multifactorial process involving neurodegeneration and progression of the disease 3,4 .
For multifactorial diseases like PD that have intertwined pathophysiological routes, the "one gene, one target, one drug" paradigm might be inadequate to achieve the required therapeutic effect. Therefore, multi-target treatment strategies like combination therapy (also known as cocktail-drug therapy) that combines two or more www.nature.com/scientificreports/ drugs with an independent mechanism of action to obtain either an added or a synergistic effect are common. For instance, levodopa (l-Dopa) often is combined with a DA receptor agonist, a monoamine oxidase-B (MAO-B) inhibitor, or a catechol-O-methyltransferase inhibitor to allay motor symptoms in advanced PD. Such combination nevertheless holds drawbacks, such as polypharmacology-related adverse drug reactions due to varying pharmacokinetic and pharmacodynamics profiles of each drug and medication non-adherence and noncompliance. Thus, another multi-target strategy in which a single chemical entity can influence more than one target has piqued the interest of medicinal chemists. This approach can avoid adverse effects associated with combination drugs and attain a more predictable pharmacokinetic profile of a single drug compared to the multiple drugs administered in combination therapy 3,5 .
Monoamine oxidases (MAOs) are the principal flavoproteins responsible for the catalytic breakdown of monoamine neurotransmitters, and two isoforms of MAO (MAO-A and MAO-B) are present in human tissues. Serotonin and norepinephrine are selectively metabolized by MAO-A, whereas phenylethylamine and benzylamine are preferentially degraded by MAO-B. DA, tyramine, and tryptamine are deaminated by both forms of MAOs 15 . Thus, the inhibitors of MAO can be used as both prophylactic and therapeutic agents in neurodegenerative disorders, including AD, PD, anxiety, depression, and schizophrenia, where a loss of monoamines is observed [16][17][18][19] . Earlier studies have shown that ILG can inhibit MAOs, however, its inhibitory potential varies according to animal models. ILG displayed a weak inhibition of MAO from bovine serum (50% inhibition of MAO activity at concentrations > 200 µM) 20 . Whereas, ILG isolated from Sinofranchetia chinensis demonstrated significant inhibition of rat brain MAO-A and MAO-B (rMAO-A/B), with IC 50 values of 13.9 and 47.2 µM, respectively 21 . ILG inhibited human MAO-A (hMAO-A) and human MAO-B (hMAO-B), with IC 50 values 41.5 and 21.8 µM, respectively 22 .
A study by Zhuo et al. on the hMAO inhibitory activities of ILG derivatives found that ILG is selective toward hMAO-B, in contrast to higher selectivity toward rMAO-A. The synthesized compound C8 showed higher potency and selectivity for hMAO-B than ILG; the enzymatic activity of C8 against hMAO-B was about 16 times greater than that of ILG. However, molecular docking showed greater proximity of ILG toward the catalytic site active residues of hMAO-B than of C8 [22][23][24] . To examine the incongruent results between the enzymatic activity and the in silico docking study and to investigate the reproducibility of the inhibitory action of ILG on recombinant hMAOs, we re-evaluated the functions of ILG on hMAO-A and hMAO-B using a homogenous luminescent assay, enzyme kinetics, and molecular simulation studies. The present work includes a detailed study of inhibition kinetics and an in silico docking study focusing on the ILG-hMAO (A/B) complex, which has not been reported.
ILG was found to suppress cocaine-induced DA release in the rat brain by modulating the gamma-aminobutyric acid-B (GABA B ) receptor 25 and also, inhibit N-methyl-d-aspartate receptor-induced Ca 2+ influx 26 . Pretreatment with ILG controlled the hyperlocomotion caused by methamphetamine (METH) treatment in mice but did not cause a significant change in the monoamine levels in tissues of the cerebral cortex, striatum, nucleus acumen, thalamus, or hypothalamus 27 . METH is a psychoactive agent that increases the extracellular DA level either by suppressing DA reuptake or internalization of dopamine transporter (DAT) from the plasma www.nature.com/scientificreports/ membrane or stimulation of DA efflux, and results in rewarding and addictive effect 28 . Repeated administration of METH was demonstrated to affect the nigrostriatal DA pathway causing the neurodegeneration of dopaminergic terminals and subsequent reduction in tyrosine hydroxylase, DAT, and striatal DA levels, predisposing the patient to PD. Dopamine receptors (DARs) mediate various cognitive and behavioral functions associated with DA, and DA D 1 type receptor (D 1 R and D 5 R) antagonism was found to prevent amphetamine-induced striatal neurodegeneration 29,30 . In addition, Everett et al. demonstrated that V 1A R plays a substantial role to mediate the inhibitory action of oxytocin on METH-primed reinstatement and drug-seeking behaviors 31 . Thus, our present study aimed to assess the plausible role of ILG on G-protein coupled receptors (GPCRs), such as DARs (D 1 R, D 2 R, D 3 R, and D 4 R), and the vasopressin receptor (V 1A R) via radioligand binding and functional GPCR assays.

Results
Inhibitory activity of ILG on recombinant human monoamine oxidases. ILG exhibited a strong inhibitory effect against hMAO-A and hMAO-B, as indicated by the low IC 50 and K i values (Table 1). ILG inhibited hMAO-A with an IC 50 0.68 ± 0.03 µM in a competitive manner and had a K i value of 0.16 ± 0.00 µM. Whereas, it showed mixed inhibition of hMAO-B, with an IC 50 0.33 ± 0.02 µM, and K ic and K iu values of 0.094 ± 0.00 and 0.71 ± 0.01 µM, respectively. The modes of enzyme inhibition were determined by enzyme kinetics analyses and were represented graphically by Lineweaver-Burk, Dixon, and secondary plots (Fig. 2). The results indicated that ILG can suppress the activity of hMAO-B more efficiently than that of hMAO-A provided that a lower concentration is required to achieve the same level of activity (IC 50 values); as K ic < K iu , the compound has a greater affinity to inhibit unbound enzymes compared to substrate-bound enzymes.
Computational interaction study of ILG with hMAOs. The overall docking results including the binding energy and the interacting residues of hMAO-A and hMAO-B when ILG binds to these enzymes are www.nature.com/scientificreports/ shown in Table 2 and Fig. 3. To optimize and verify the docking simulation study, the reference ligand, harmine, was docked to hMAO-A (2BXR), while C17 was docked to hMAO-B (2V60). From the best pose obtained for the ligand-enzyme complex, the binding modes of ILG at active sites of the hMAO isoenzymes were visualized. ILG interacted with the catalytic site residues of hMAO-A with a binding energy of − 7.44 kcal/mol, which was lower than that required for the standard ligand harmine (− 6.46 kcal/mol). The docking pose revealed that the 2′,4′-dihydroxyphenyl moiety of ILG occupied the central catalytic site of the hMAO-A and interacted with Tyr407 via a π-π T-shaped hydrophobic bond, and with Tyr444 and flavin adenine dinucleotide (FAD) via H-bonds. A carbonyl group at C1 was aligned toward FAD, which enabled its interaction with FAD through an H-bond. The aromatic ring of the 4″-hydrophenyl moiety engaged in hydrophobic interactions with active site residues Ile335, Leu337, and Phe208 through π-σ and π-π stack bonds. A H-bond between 4″-OH and Thr336 was also formed in the ILG-hMAO-A complex.
Unlike hMAO-A, the orthosteric site structure of hMAO-B entails two cavities: a substrate cavity lined by numerous aromatic and aliphatic amino acid residues and an entrance cavity that lies adjacent to the substrate cavity. An aromatic cage formed by Tyr398 and Tyr435 along with FAD represents the recognition site for catalysis, while Tyr 326, Ile199, and Leu171 act as the key gating residues between the entrance and substrate binding Table 2. Binding energy and interacting residues during the inhibition of hMAO-A and hMAO-B by isoliquiritigenin and reference ligands. a The estimated binding energy (kcal/mol), which signifies the binding affinity of a ligand to the active site of hMAO-A and hMAO-B. b The interacting amino acid residues in the ligand-enzyme complex were determined using the AutoDock 4.2 program. c Harmine, reference ligand;

Ligands
Binding energy a www.nature.com/scientificreports/ cavity 23 . Molecular analysis revealed that ILG could bind with the residues of the substrate-binding cavity and the hydrophobic cavity between the substrate cavity and entrance cavity of the enzyme with a binding energy of − 8.69 kcal/mol. The 4″-hydroxyphenyl moiety of the compound interacted with Tyr398 and FAD, while another aromatic ring (2′,4′-dihydroxyphenyl moiety) was aligned to interact with gating residues Leu171 and Ile199. Pi-sulfur and H-bond interactions were formed with Cys172, a residue in the hydrophobic cavity, via 2′,4′-dihydroxyphenyl moiety and 1-CO group, respectively. In addition to catalytic interactions, ILG formed bonds with allosteric site residues of hMAO-B, with − 8.03 kcal/mol of binding energy (Fig. 3c,f). Polar H-bond associations with Pro104, His115, Asn116, and Asp123 and non-polar interactions with Val106, Trp119, and Arg120 were observed at the allosteric site.

hMAO-B (2V60)
Receptor binding profile of ILG. The binding affinity of ILG to the human dopaminergic D 1 , D 2 , D 3 , and D 4 receptors was measured in vitro by displacement of antagonist radioligands [ 3 H]SCH23390 and [ 3 H]methylspiperone from the recombinant Chinese hamster ovary (CHO)-D 1 R, human embryonic kidney (HEK)-D 2L R, CHO-D 3 R, and CHO-D 4 R cell membranes. The V 1A R binding property was evaluated by the agonist radioligand binding assay using [ 3 H] arginine vasopressin (AVP) on the V 1A R-transfected CHO cells and validated using AVP as a reference. Screening the binding property of ILG at 100 μM revealed significant inhibition of control-specific binding (> 50%) for D 1 R, D 3 R, and V 1A R and moderate inhibition of D 4 R with 44.8% inhibition of control-specific binding (Table 3). High radioligand displacement was observed on the CHO-D 1 R membrane (96.0%), whereas no binding affinity was seen on the HEK-D 2L R membrane.
Modulatory action of ILG on D 1 , D 3 , and V 1A receptors. The preliminary in vitro binding screening using 100 μM ILG indicated a high binding affinity for the D 1 R, D 3 R, and V 1A R receptors. Since binding assays are limited to determining the affinity and selectivity of ligands to receptors, we conducted cellular functional assays to reveal agonist or antagonist behavior of ILG on D 1 R, D 3 R, and V 1A R. Functional GPCR assays showed that ILG is a D 1 R antagonist and a D 3 R and V 1A R agonist. As shown in Table 4, 100 µM ILG inhibited the D 1 R control agonist response by 91.9 ± 3.1% and stimulated D 3 R and V 1A R by 92.7 ± 0.5% and 73.3 ± 4.0% of the control agonist response, respectively. Figure 4 shows the dose-response curve obtained as the percent inhibition of the control agonist response on D 1 R and the percent of the control agonist response on D 3 R and V 1A R. The concentration required to produce 50% of agonist response (EC 50 ) or antagonist response (IC 50 ) of ILG along with those of the reference drugs on the tested receptors are presented in Table 4.    Table 5. ILG was located within the ligandbinding site of a human model of D 1 R (hD 1 R) and interacted with residues Asn292, Asp103, and Ile154 through H-bonds to three hydroxyl groups at C2′, C4′, and C4″ of ILG. Two phenyl rings of the chalcone moiety anchored with conserved residues Ile104, Phe288, and Ser198 of transmembranes (TM) 3, 6, and 5, respectively, via strong hydrophobic connections. Molecular docking of ILG with D 3 R revealed that ILG can bind to the orthosteric binding site (OBS) of D 3 R with a lower binding energy (− 7.24 kcal/mol) than that required for the reference agonist DA (− 5.47 kcal/mol). Polar interaction of 4′-OH was found with Asp110 of TM3, which was similar to the case of DA but without the formation of a salt bridge. Major interactions were observed with the amino acid residues of TMs 3 and 5, such as Val107 (via π-alkyl bonding to the 2′,4′-dihydroxyphenyl ring) and Val189 (via H-bond to 4′-OH and π-alkyl bonding to the 4″-hydroxyphenyl ring). An electrostatic π-cation linkage with the key OBS residue His348 was formed with the 4″-hydroxyphenyl moiety. Moreover, ILG extends toward the second extracellular loop (ECL2) and binds to Cys181 via a H-bond and to Ile183 via a non-polar bond.
Similarly, ILG was predicted to be buried within helices 1, 2, and 3 of V 1A R with a binding energy of − 7.03 kcal/mol, which is lower than that required for AVP binding (− 5.98 kcal/mol). In silico docking displayed five hydrophilic connections of ILG with the residues of V 1A R that were assisted by the 2′,4′,4″-hydroxy functional groups and a carbonyl group at C1. An electrostatic interaction of the 4″-hydroxyphenyl moiety with a nonconserved hydrophilic residue (Asp112 of TM2) was observed. Hydrophobic interactions of the 2′,4′-dihydroxy phenyl ring of ILG occurred with the Val127 and Lys128 residues of TM3. Of all molecular interactions, Glu54, Asp112, and Lys128 were the common binding residues compared with the reference agonist AVP.

Prediction of the drug-likeness and ADMET profile of ILG. Pharmacokinetic parameters and toxic-
ity of ILG were predicted using SwissADMET and pkCSM applications. SwissADMET showed drug-likeness and lead-likeness of ILG with a lipophilicity of 2.37 (log Po/w) and a solubility of − 3.06 (soluble). A pkCSM application indicated high intestinal absorption (HIA) of ILG (> 90%) along with its probability to cross the blood brain barrier (BBB) and reach the CNS. It also predicted ILG as a safe drug-like molecule with no toxicity to hepatocytes or upon Ames testing (Table 6).

Discussion
Of the multiple kinds of phytochemicals, flavonoids have been known to possess significant pharmacological action against different NDDs, including AD and PD 32 . ILG is a plant-derived chalcone that has shown promising effects against NDDs on multiple investigations 9 . Though earlier records of the MAO inhibition potential www.nature.com/scientificreports/ of ILG are available, we re-evaluated the human recombinant MAO inhibition potential of ILG obtained from Glycyrrhizae radix together with the investigation of enzyme kinetic parameters and the computational docking study to establish its mechanism of action.
In the in-vitro recombinant hMAO inhibition assay, ILG exhibited a remarkable inhibition of both hMAO-A and hMAO-B, with respective IC 50 values 0.68 ± 0.03 and 0.33 ± 0.02 µM. This activity is significantly higher than the formerly reported activity by Zhou et al. at 61 times more potent against hMAO-A and 66 times more potent against hMAO-B. Nonetheless, in both of these studies, we found an approximately twofold higher selectivity of ILG for hMAO-B than for hMAO-A. Mohamed and colleagues had earlier investigated hMAO inhibition of ILG isolated from Colvillea racemosa and found selective hMAO-B inhibition with an IC 50 value of 0.51 ± 0.01 µM, while that for hMAO-A inhibition was 22.66 ± 1.84 µM 33 . The variations in the potency of ILG for hMAO inhibition might have occurred due to different experimental conditions and procedures. To validate our results, we used the standard hMAO-A inhibitor clorgyline⋅HCl and the hMAO-B inhibitor l-deprenyl⋅HCl.
Although reports on hMAO inhibition by ILG are available, a detailed study on the mechanism of hMAO inhibition is lacking. Thus, our study incorporated an enzyme kinetics study of ILG for inhibition of hMAO isoenzymes and demonstrated concordant inhibition constant (K i ) values with respect to the IC 50 values for each enzyme (Table 1). In a previous study, ILG showed non-competitive inhibition of rMAO-A and mixed inhibition of rMAO-B 21 . The present study found competitive inhibition of hMAO-A by ILG along with a similar mode of action against hMAO-B as in rMAO-B. Though rat and human MAO-A have a 92% similar sequence identity, crystallographic analysis has revealed that a significant structural difference between hMAO-A and rMAO-A exists in the conformation of residues 210-216, which are vital for the structure of the hMAO-A active site. Unlike hMAO-A, rMAO-A is a dimer and has a smaller substrate cavity. A recent study of hMAO inhibition found liquiritigenin (LG) to have very potent action against hMAO isoenzymes compared to its activity against rMAO isoenzymes 34 . In addition, screening of coumarin and 5H-indeno[1,2-c]pyridazin-5one derivatives for inhibition of MAOs of human and rat showed greater effectiveness of compounds against hMAO-B than against rMAO-B and no correlation between human and rat pIC 50 values 24,35 . These findings support that results drawn from experimentation in one animal model cannot consistently be extrapolated to humans.
Our docking study showed good proximity of ILG toward the isoalloxazine of the FAD cofactor in hMAO-A and hMAO-B. ILG bonded with key residues of the hMAO isoenzymes through multiple polar and non-polar interactions that resulted in strong inhibition. In hMAO-A, the 2′,4′-dihydroxyphenyl ring was engaged in interaction with the major active site residues FAD, Tyr444, and Tyr407, while in the case of hMAO-B, the 4-hydroxyphenyl moiety was responsible for forming strong interactions with FAD and a catalytic residue Tyr398. LG, a structurally interrelated flavanone form of ILG, was reported to have strong inhibition of hMAO-A and hMAO-B, with IC 50 values of 0.27 ± 0.04 and 0.098 ± 0.00 µM, respectively. Although these IC 50 values suggest LG to be more potent than ILG, both compounds exhibited a greater affinity for hMAO-B as the binding energy for the hMAO-B was lower than that for hMAO-A 34 . 3-Deoxysappanchalcone, a 2′-methoxy-4′,4″-dihydroxychalcone, is related structurally to ILG and demonstrates a potent hMAO inhibition potential similar to that of ILG (IC 50 values: 10.11 and 0.68 µM for hMAO-A and hMAO-B, respectively) with a higher affinity for hMAO-B (− 8.90 kcal/mol) 36 . Other natural chalcones, such as 4-hydroxyderricin, 2,2′-dihydroxy-4′,6′-dimethoxychalcone, and broussochalcone, also have notable hMAO inhibition potential 36 . In vitro and in vivo studies have shown that ILG can protect dopaminergic neurons from METH 37 , 6-OHDA 38 , and glutamate-induced neurotoxicity 39 and alleviate cognitive impairment 12,14 , anxiety, and locomotor sensitization 40 . Since the modulation of DA and vasopressin receptors has been found to affect cognition, anxiety, and neuronal survival, we evaluated the binding affinity of ILG to the D 1 R, D 2 R, D 3 R, D 4 R, and V 1A R receptors to determine whether ILG can target these receptors and identified significant binding to D 1 R, D 3 R, and V 1A R. The functional GPCR assay showed ILG to have D 1 R antagonist and D 3 R and V 1A R agonist properties. A ligand-D 1 R interaction study revealed that ILG can bind with the active site residues Asp103, Ile104, and Phe288 similar to the reference antagonist SCH23390. It has been reported that interactions with Asp103, Asn292, and Ser198 residues are vital for strong binding affinity of the ligand to the D 1 R 41 . The molecular docking pose of D 3 R with ILG indicates that ILG extends between helices 3, 5, and 6 to the ECL2. Some ECL2 residues (182-185) form a part of the ligandbinding pocket of D 3 R, and interactions with ECL2 residues have been acknowledged for their selectivity of ligands between D 2 and D 3 receptors 42 . Thus, the binding of ILG with Cys181 and Ile183 of ELC2 might be responsible for the higher affinity of ILG for D 3 R compared to D 2 R. Cotte et al. reported that aromatic residues in V 1A R play an insignificant role in the agonist's binding affinity; however, the conserved hydrophilic residues Glu108, Lys128, and Gln185 contribute significantly to the V 1A R agonist AVP binding 43 . The H-bonding with Gln108 and non-polar linkage with Lys128, along with other known agonist binding residues (such as Glu54, Cys203, Asp112, and Val127) 44 , were detected by computational docking of ILG with a 3D model of human V 1A R.
For treatment of PD, l-Dopa is the gold standard therapy; however, prolonged use of l-Dopa causes dyskinesia and dystonia. DAR agonists stimulate the postsynaptic DARs and enhance the function of the DA system. Thus, they are considered the first choice in de nova patients to delay the onset of l-Dopa therapy 45 . Another pharmacological principle used to alleviate DA deficiency in PD is the use of reversible and irreversible hMAO-B inhibitors, such as selegiline, rasagiline, and safinamide. Clinical trials have shown improvement of motor   Table 5. Binding energy and interacting residues of isoliquiritigenin and reference ligands at the active sites of human dopamine and vasopressin receptors (hD 1 R, hD 3 R, and hV 1A R). a Reference ligands for hD 1 R, hD 3 R, and hV 1A R. www.nature.com/scientificreports/ To treat the motor and non-motor complications of PD that can arise from fluctuations in l-Dopa concentration in the plasma and the pulsatile stimulation of DARs, DAR agonists can be used as an add-on therapy to l-Dopa 48 . The D 2 /D 3 receptor agonist pramipexole was found to improve depression in PD, anhedonia, and motor deficits 49 . D 1 R activation has shown improvement in hippocampal neurogenesis and anxiolytic and antidepressant-like effects in a rat model of PD, whereas the D 1 R antagonist SCH23390 was found to inhibit stimulant-induced DAR supersensitivity 50 . Addictive psychostimulants like cocaine, amphetamines, and METH increase D 1 R activation, leading to the amplification of c-fos and fosB expression in the striatum. These elevations in c-fos and fosB along with the METH-induced increases in c-Jun, JunB, JunD, and Fra2 expression were all inhibited by pretreatment with the D 1 R antagonist SCH23390 51 . Although the expression of D 1 Rs is reduced in PD, increased stimulation of D 1 R signaling pathways has been linked to l-Dopa-induced dyskinesia 52 . The activation of D 1 Rs by METH is known to cause neurotoxicity of nigrostriatal dopaminergic neurons, and the deletion of D 1 Rs or D 1 R antagonism considerably protected against METH-induced neurotoxic effects 53 . Thus, the neuroprotective effects of ILG, which had previously been shown to impede METH-induced neurotoxicity 37 and hyperlocomotion without alterations in monoamine levels 27 , can be attributed to its D 1 R antagonist behavior and MAO-inhibition potential.
Former investigations in neuroscience have reported that D 3 R agonists (such as 7-OH-DPAT) can restore nigrostriatal integrity by inducing neurogenesis and therefore improve locomotor function in PD models 54 . Piribedil, another D 3 R agonist, was found to be effective in enhancing learning and memory in an animal model of cerebral ischemia-reperfusion 55 . Vasopressin V 1A R agonist activity has been associated with cognitive function and social recognition, while V 1A R antagonists have an anxiolytic effect. Studies conducted on V 1A R knockout mice have identified abnormalities in social interaction and recognition on social behavior tests and reduced anxiety levels in elevated plus-maze and marble-burying behavior tests 56,57 . Therefore, the modulatory effects of ILG on V 1A R might be related to reduced anxiety-like behavior, which was observed following ILG treatment in nicotine-withdrawn mice 40 .
The computational prediction of pharmacokinetic parameters such as drug-likeness, solubility, absorption, and CNS permeability for ILG shows favorable characteristics. In silico prediction also reveals no toxicity of ILG on AMES test and hepatotoxicity assessment. Experimentally, ILG has been found to be hepatoprotective 58 and induces Nfr2-dependent detoxification genes 59 . Licorice is reported to be used as an adjuvant agent to enhance the efficacy of other drugs, however, herb-drug interactions have also been expected with its prolonged use. ILG, as one of the phytoconstituents, has been found to inhibit cytochrome P450 enzymes such as CYP1A2, CYP2C19, CYP2C9, and CYP3A4. These enzymes are responsible for major drug metabolism, and inactivation of these enzymes can reduce drug metabolism, thus, leading to increased plasma concentration and the risks of side effects 60 . Treatment with higher doses (36 and 100 µM) of ILG inhibited follicular growth and steroidogenesis owing to the perturbation in the expression of key steroidogenesis regulators such as CYP17A1, CYP19A1, and HSD17B1 61 .
Overall, this study reinvestigated the hMAO-A and hMAO-B inhibitory activity of ILG and established the enzyme inhibition mode through kinetic assays and computational docking. Our study recognized that ILG is a potent hMAO-inhibitor showing competitive inhibition against hMAO-A and a mixed mode of inhibition against hMAO-B. ILG was observed to have a significant binding affinity for the D 1 , D 3 , and V 1A receptors. Probing the modulatory activity of ILG on DA (D 1 and D 3 ) and vasopressin (V 1A ) receptors revealed that ILG has an antagonist effect on the D 1 receptor and agonist effects on the D 3 and V 1A receptors. The multi-target nature of ILG together with the prediction of suitable pharmacokinetics and the toxicity profile renders this compound as a potential flavonoid for the management of PD and its related neurological symptoms. Future experimental studies on ILG's pharmacological and toxicological properties and adverse effects will be important to confirm its therapeutic efficacy and benefits.  Table 6. Prediction of the pharmacokinetic parameters and toxicity of isoliquiritigenin. a Octanol-water partition coefficient. b LogS scale: Insoluble < − 10 < Poorly < − 6 < Moderately < − 4 < Soluble < − 2 < Very < 0 < Highly. c Human intestinal absorption: If < 30%, poorly absorbed. d Log BB > 0.3 considered to readily cross the blood brain barrier, while log BB < − 1 considered to be poorly distributed to the brain. e Log PS values > − 2 are considered to penetrate the CNS, while log PS values < − 3 are considered unable to penetrate the CNS. www.nature.com/scientificreports/ In vitro hMAO-A and -B inhibitory assay and enzyme kinetics. Experimental conditions and procedures for this experiment were as described 62 . Briefly, we added 12.5 µl of the test compound or l-deprenyl to an aliquot of 12.5 µl of beetle luciferin derivative substrate (initial concentrations of 160 µM and 16 µM for hMAO-A and hMAO-B, respectively) in each well of the plate. A 25-µl enzyme solution was added to the test samples to initiate the reaction. After an hour of incubation of the mixture at 25 °C, 50 µl reconstituted luciferin detection reagent was added to every well to stop the hMAO reaction, and an additional 20 min of incubation at 25 °C was performed. Then, the luminescence reading was taken on a FilterMax F5 Multi-Mode microplate reader (Molecular Devices, LLC., CA, USA). The enzyme inhibition kinetics were analyzed using varying concentrations of hMAO substrate (40-160 μM for hMAO-A and 4-16 μM for hMAO-B) and ILG (0-1 μM), as shown in Fig. 2. Inhibition constants (K i ) for each enzyme inhibition were obtained from secondary plots analyzed using SigmaPlot 12.0 TM software (SPCC, Inc., Chicago IL, USA).

Druglikeness
Radioligand binding assays. The binding assays were carried out using validated methods and standard operating procedures used by Eurofins Cerep (catalog items 0044, 1405, 0048, 0049, and 0159). For the human D 1 R, D 3 R, and D 4 R binding assays, membrane homogenates (8-80 µg) of the respective receptor-expressed CHO cells were suspended in a buffer containing 50 mM Tris-HCl, 5 mM KCl, 5 mM MgCl 2 , 1.5 mM CaCl 2 /120 mM NaCl, and 5 mM EDTA and incubated for 60 min at 22 °C with appropriate radioligand [ 3 H]SCH23390 or [ 3 H] methylspiperone either in the presence or absence of the test compound. The human D 2L R binding experiments were run on membrane preparations prepared from HEK cells. The plasma membrane homogenates (40 µg) of D 2L R-expressing HEK cells were incubated with 0.3 nM [ 3 H] methylspiperone for 60 min at 22 °C after being suspended in a binding buffer containing 50 mM Tris-HCl, 5 mM MgCl 2 , 1 mM EDTA, 1 UI/mL ADA, 1 µg/ml leupeptin, 1 µM pepstatin, and 10 µg/mL trypsin inhibitor. Human V 1A R binding was assessed using membrane homogenates of transfected CHO cells suspended in a buffer containing 5 mM Tris-HCl, 5 mM MgCl2, and 0.1% BSA and incubated with [ 3 H]AVP at 22 °C for 60 min. 1 µM SCH23390 (for D 1 R), 10 µM (+)-butaclamol (for D 2L R, D 3 R and D 4 R), and 1 µM AVP (for V 1A R) were used to define non-specific binding.
Following incubation, the binding reaction within the samples was terminated by vacuum filtration through 0.3% polyethyleneimine-treated glass fiber filters (GF/B, Packard). The filters were washed several times with an ice-cold wash buffer (50 mM Tris-HCl) using a 96-sample cell harvester (Unifilter, Packard) and dried. Scintillation cocktail (Microscint O, Packard) was added to the dried filters, and their radioactivity was determined using a scintillation counter (Topcount, Packard). The radioligand binding results were expressed as the percent inhibition of control specific binding given by the following equation: 100 -[(measured specific binding/control specific binding) × 100].
To verify the results, the standard compounds were tested at different concentrations in each of the binding assays to obtain competition curves from which their IC 50 values were determined. The binding affinity of ILG to the receptors was screened at 100 μM.
In vitro functional GPCR assay. Cellular and nuclear receptor functional assays were performed at Eurofins Cerep using human recombinant CHO cells transfected with the GPCR genes of interest (D 1 R, D 3 R, and V 1A R). The in-house assay protocol for the assay was as described in our earlier report 63 . The results of the functional assay were based on the measurement of effects on the cAMP level and calcium ion mobilization.

Measurement of cAMP level.
Stably transfected CHO cells containing the cDNA of human D 1 and D 3 receptors were suspended in a medium containing an HBSS buffer that had been supplemented with 20 mM HEPES buffer and 500 µM IBMX. This cell suspension was distributed into assay plates at a density of 5 × 10 3 cells/well and incubated at 25 °C (for the D 1 R transfected cell suspension) or 37 °C (for the D 3 R transfected suspension) for 10-30 min both with and without ILG or standard. Then, the D 2 -labeled cAMP conjugate was dispensed into the cell plate followed by the addition of europium cryptate-labeled anti-cAMP antibody. A lysis buffer was added to each well, and the resulting mixture was incubated for 1 h at 25 °C. The homogeneous timeresolved fluorescence (HTRF) reading was recorded using a PerkinElmer Envision microplate reader (Waltham, MA, USA) at an excitation intensity of 337 nm and emission intensities of 620 and 665 nm. The cAMP level was calculated as the ratio of the signal measured at 665 nm to that measured at 620 nm. The final results were illustrated as the percentage of the control agonist response and as the percentage inhibition of the control agonist response.

Measurement of intracellular calcium levels.
Calcium ion influx was determined fluorimetrically to establish the functional effect of ILG on V 1A R. Transfected CHO cells were suspended in an HBSS/20-mM HEPES buffer and distributed into microplate wells at a density of 1 × 10 5 cells/well. A fluorescent probe (Fluo8 Direct, Invitrogen, Carlsbad, CA, USA) was mixed with probenecid in HBSS/20 M HEPES (pH: 7.4) discretely and applied to each well, and left to equilibrate with the cells at 37 °C for 60 min. Thereafter, the assay plate was kept in a microplate reader (CellLux, PerkinElmer, Waltham, MA, USA), and ILG, the reference agonist, or HBSS buffer (control) was added to the plate; finally, the fluorescence was measured. The agonist effect was determined as the percentage of the control response to 1 µM AVP, while the antagonist behavior was calculated as the percentage inhibition of the control response to 10 nm AVP.
Homology modeling. The sequences for hD 1 R and hV 1A R were retrieved from the UniProt database with respective IDs of P21728 (DRD1_HUMAN) and P37288 (V1AR_HUMAN). Due to the high sequence similar-