schizophrenia ( sz ) affects approximately 1% of the world s population , having similar prevalence throughout diverse ethnic groups.1 sz is a chronic , mental , highly heritable and common disease characterized by neurophysiology impairment and neuropsychological abnormalities.13 the susceptibility of sz is influenced by environment factors , polygenic components , and their interactions.4 molecular mechanisms activating sz remains intangible . sz gene identification is particularly exigent and challenging due to the fact that the disorder may include limited sz diagnostic accuracy as phenotypic definition and numerous entities have not been yet defined . additionally , the lack of conclusive linkage from scans of the genome could also be due to the existence of several sz vulnerability genes that are difficult to replicate and detect.5 the variations in d - amino acid oxidase activator ( daoa ) ( 13q34 ) gene initially were linked with sz.6 furthermore , this gene has been linked with other phenotypes and psychiatric disorders such as severe bipolar disorder,7 bipolar disorder,8 and major depression.9 
 daoa genetic variations have contributed to different central nervous system ( cns ) disorders linked with dysfunction of glutamatergic signaling.6,10,11 the longest open reading frame of daoa ( g72 ) is predicted to encode a putative 153 amino acid protein isolated from testis , spinal cord , caudate nucleus , and amygdale libraries.6 daoa ( g72 ) , a primate - specific protein , has been characterized as dao activity modulator.6 
 daoa expression in transgenic mice induced schizophrenic - related behavioral phenotypes.12,13 sz patients show daoa overexpression in the dorsolateral prefrontal cortex when compared with healthy controls.14 sz susceptibility genes have been identified in genetic studies,1518 but genetic interactions among sz genes and their interplay with clinical subtypes and neurobiological abnormalities remains obscure . the product of dao is an enzyme that degrades d - serine amino acid which acts as a co - agonist at the glycine site of the n - methyl - d - aspartic acid ( nmda ) receptors.19 the product of daoa ( g72 ) activates the dao enzyme.6 the biological function of daoa and dao are involved in the hypothesized hypofunction of nmda receptor complex as the prospective pathogenesis of sz ( figure 1).20 the neurotransmission of nmda is the dominant molecular mechanism for memory , cognition , and synaptic plasticity . numerous psychiatric and neurological disorders are linked with dysfunction of nmda receptor  mediated neurotransmission.21 hyperactivity and overexpression of brain dao have been associated with sz.22,23 dao hyperfunction resulting from the decreased level of d - serine and nmda receptor hypofunction has been suggested to be an underlying dysregulation in this disorder.24,25 in the cerebrospinal fluid26 and serum of schizophrenic patients , the level of d - serine is reduced . in addition to its role in cognition and psychosis , dao may be involved in amyotrophic lateral sclerosis , resulting in d - serine accumulation , motoneuron death , and excitotoxicity.28 dao activity manipulation could be an approach to regulate the level of d - serine29 ( figure 1 ) . daoa has been identified in many different parts of the cns , including spinal cord , nucleus caudatus , and amygdala.26 it has also been implicated in the regulation of dendritic branching and mitochondrial functions.30 daoa was identified to activate dao , while recent studies showed that daoa modulates human dao function as a negative effector.26,31 dao catalyzes the oxidative deamination of d - amino acids such as d - serine to a keto - acids and d-3,4-dihydroxyphenylalanine ( d - dopa).32,33 l - dopa enters in the basic biosynthetic pathway to homovanillic acid ( hva ) and dopamine . dopamine--hydroxylase converts dopamine into noradrenaline , which is degraded into 3-methoxy-4-hydroxyphenylglycol ( mhpg ) after entering its metabolic pathway . kinetic data showed that the maximal velocity for the oxidative deamination of d - serine is much lower than for d - dopa.19 d - serine is an allosteric modulator of nmda - type glutamate receptor ( nmdar),19 having a d - serine modulatory site . d - serine occupation at this site is required for glutamate to stimulate the flow of cations.34,35 study on the interactions between noradrenaline36 and glutamate propose that daoa may be linked with noradrenaline via glutamatergic mechanism . a bidirectional interaction between the dopamine system and nmdar the activation of nmdar leads to enhance the recruitment of the dopamine d1 receptor ( drd1 ) to the plasma membrane.37,38 there is also a direct protein  protein coupling between nmdar an drd1.36,39,40 it has also been suggested that nmdar and drd1 , early after their biosynthesis , form heteromeric complexes , which are transported as performed units to the plasma membrane.41 nmda antagonists lead to the increase of dopamine neuron firing rates in mid brain.42 a direct association between dopamine release and d - serine has been shown , as high doses of d - serine attenuated amphetamine - induced dopamine release.31 the present work focused on the ligand analysis to reveal novel inhibitor against sz . experimental validation of daoa isoforms by using x - ray crystallography and nmr ( nuclear magnetic resonance ) are not available yet . so as to solve structural insights , the isoforms of daoa three - dimensional structures were predicted by applying template crystal structures from protein data bank ( pdb ) . the inclusive in silico analysis may provide an evidence for a reliable framework that could assist in the designing and development of lead molecules that could be potential drug targets against sz . the significant role of daoa on dao regulation has been suggested;6 the specific role daoa plays in regulating the activity of dao in schizophrenic patients is still controversial . sacchi et al31 demonstrated that a 1:2 molar ratio of daoa to dao is optimal for the inhibition of dao . on the other hand , chumakov et al6 revealed that daoa enhances the activity of dao , resulting in a reduction of d - serine concentration . daoa is elevated in the serum of nonmediated and mediated patients of sz,6 despite the fact that daoa and daoa mrna are hardly detectable in the brain.31,43 as a result , dysfunction of nmda and its related sz pathological symptoms may result from daoa overexpression that changes the activity of dao . sehgal et al44 explored the protein interaction between daoa and dao by the exercise of molecular docking . four isoforms of daoa ( isoform 1 , isoform 2 , isoform 3 , isoform 4 ) were selected from the uniprot protein database . conserved amino acid residues are found in the selected four daoa isoforms suggesting that the conserved region is important for daoa folding and function . a daoa transcribing two transcripts led to the production of four isoforms collectively , namely , daoa-153 , daoa-126 , daoa-82 , and daoa-125 . while daoa-125 has the accession number a2t115 , the rest of the three isoforms were represented by the single accession number p59103 in uniprot . in the current study , structure prediction , sequence analysis , phylogenetic analysis , docking studies , and comparative proteomics analysis were performed on hp core - i-5 workstation . the amino acid sequences of three known isoforms of daoa and one another variant of daoa were retrieved for comparative modeling . the amino acid sequences of daoa isoforms were retrieved from uniprot and subjected to a protein protein blast for the identification of suitable template structure against pdb.45 the automated protein modeling program modeller 9v1046 was employed to predict the three - dimensional structures of daoa by satisfying spatial restraints . five evaluation tools were used : rampage,47 procheck,48 anolea,49 and errat50 were applied to assess the predicted models , and the structures were further evaluated by molprobity.51 finally , the poor ramachandran outliers and rotamers were removed by utilizing wincoot52 tool . sehgal et al44 utilized a ligand daoa : o , o - diethylthiophosphoryl ( z)-2-(2-aminothiazol-4-yl)-2-trityloxyiminoacetate ( c28h28n3o5ps2 ) to analyze the binding sites of daoa ( accession number a2t115 ) . the ligand ( c28h28n3o5ps2 ) was utilized for protein  ligand docking analysis using the autodock tools.53 the number of h - bond donors , rotatable bonds , and h - bond acceptors were obtained by employing pubchem,54 mcule , and molinspiration . mcule and cresset were employed to estimate the mutagenesis of novel molecules , and no mutagenic risks were detected . the aim of the docking analysis was to identify the relative binding specificities and binding pattern of protein against the specific ligand ( c28h28n3o5ps2 ) . the geometrical optimization of three - dimensional structures of proteins and ligand were performed using chemdraw ultra and vega zz . admet refers to the absorption , distribution , metabolism , excretion , and toxicity properties of compound and was analyzed by utilizing admetsar online server . p - distance ( difference between two amino acid sequences on the basis of their amino acid positions ) , complete - deletion , and uncorrected options were utilized for analysis . national center for biotechnology information ( ncbi ) , ensembl , ucsc genome browsers , and biological databases were utilized to analyze all the mentioned species in these databases . to analyze the daoa ancestral relationship , molecular evolutionary genetic algorithm ( mega 5 ) tool was employed . a considered bootstrap value of 1,000 replications was applied and using neighbor - joining ( nj ) approach to construct the phylogenetic tree . ensembl , ncbi , and ucsc genome browser were utilized to search for evidence of daoa paralogs and orthologs . the protein sequences were run through ensembl blat and ncbi blast against all the species included in ensembl and primates . the objective of our study was based on the relation of daoa isoforms with sz and its bioinformatics analysis for analyzing the ligand . the optimally aligned templates against all the daoa - selected proteins with maximum identity , and query coverage are mentioned in table 1 . the suitable , scrutinized templates were utilized for three - dimensional structure predictions of daoa isoforms . consequently , selected templates gave better evaluation results . evaluation tools showed the reliability and efficacy of daoa three - dimensional predicted structures ( figure 2 ) . favored and outlier regions of predicted models were detected in ramachandran plot and mentioned in table 1 . no amino acid residue was observed in outlier region , and all the amino acids of the daoa-82 were observed in favored and allowed region . in daoa-125 , only one residue was observed in outlier region . poor rotamers and outliers were corrected to refine the predicted models of daoa . however , the docking analysis of ligand revealed variations in their binding energies . initially , the docking analysis was performed with ten poses , and 150 runs were saved , out of which the best poses with lowest binding energy were chosen from daoa isoforms . our findings indicated that the ligand ( daoa : o ) efficiently binds to daoa isoforms ( table 2).44 ligand docking analysis with daoa isoforms were analyzed on the basis of their binding energy values and drug properties ( table 3 ) . the ligand compound may be considered as potential antischizophrenic agents . among all the selected isoforms , autodock tool was utilized for docking analysis , and top complexes of daoa isoforms having lowest binding energies were selected for further analysis . perhaps , the stability of the ligand may be due to the allowing of more conserved binding affinities . it was revealed that the daoa isoforms  docked complexes with analyzed ligand molecule showed reliable results . it was observed that the selected ligand binds at the conserved region of daoa isoforms and revealed the binding pocket . it is also possible that the combination of observed binding residues leads to the lower binding energy achieved in the current study . it was observed that ser-100 of daoa-153 , ser-29 of daoa-82 , and ser-72 of daoa-125 ( figure 4 ) having good binding interactions with the same docked ligand . the residues having different positions due to variation in the size of isoforms are indicated with matching superscript letters in table 2 . the conserved region has same binding residues but have different positions due to different size of isoforms ( figure 5 ) . in an effort to comprehend the better interactions observed between ligand and amino acid residues in the active site of protein , a plot of amino acids  ligand interactions were generated employing ligplot and chimera 1.6 as shown in figure 6 . this emerging inclination has immense importance in decreasing the time required and amplifying the selected compound with better biological activity and minimal side effects for specific disease targets . the chemical structures are evaluated to be an efficacious drug - like compound subjected to lipinski s rule of five55 and for their oral bioavailability . admet refers to the absorption , distribution , metabolism , excretion , and toxicity properties of compound and was analyzed by utilizing admetsar online server . eleven mathematical models ( ames toxicity , blood  brain barrier penetration , fish toxicity , aqueous solubility [ logs ] , honey bee toxicity , carcinogens , acute oral toxicity , caco2 permeability , cytochrome p450 2d6 inhibition , rat acute toxicity , and human intestinal absorption ) were employed to predict the admet properties of selected compound . various toxicities are predicted often and are used in drug designing ( table 3 ) . these toxicities help in evaluating intermediates , metabolites , and pollutants , along with adjusting the dose range for animal assays . the aqueous solubility ( at 25c ) prediction showed that the selected compound is soluble in water . lipophilicity measurement ( logp ) solubility is the ratio of compound in octanol compared to its solubility in water . it was observed that the compound showed less logp value and follows lipinki s rule of five , thus showing better oral bioavailability . the excretion process that eliminates the compound from the body depends on logp.56 the compound must be absorbed by human intestine . the intestinal human absorption results showed that the compound could easily be absorbed by the intestine . the selected molecule was found to be a noninhibitor of cytochrome p450 2d6 , which indicated that the analyzed molecule may be well metabolized in phase i metabolism . the cytochrome p450 2d6 is one of the significant enzymes involved in drug metabolism.57 carcinogenicity and toxicity risk assessment were analyzed for the selected molecule , which showed it to be noncarcinogenic . the results of docking analysis suggested that an ideal ligand must be the one that comes under the designed parameters for satisfying its drug likeness and has a less docking score . coinciding with the statement and satisfying the least docking energy , it is suggested that c28h28n3o5ps2 is a potential drug molecule for sz treatment that works by targeting daoa isoforms . our docking results revealed that binding residues mentioned in table 2 are crucial , and mutational studies of these binding residues could be highly effective in further studies . it stands to reason that the selected ligand molecule in this study has the propensity to be a good candidate for the treatment of sz by targeting daoa . the proposed and used strategy led to simplification of the process of novel drug designing as biological investigation and drug analysis can be carried out on candidates with good results and without the need to waste money and time on poor activity compounds . protein data sets for a huge range of vertebrate and invertebrate genomes are presently available for analyzing the phylogenetic history of candidate genes of sz . the phylogenetic nj tree presented in current analysis reveals many interesting characteristics of the candidate gene in vertebrates ( figure 7 ) . there are approximately 80 genomes of vertebrates that are completely sequenced , and sequences are embedded in the ensembl genome browser . the blat alignment tool of ensembl was employed for sequence alignment of daoa canonical protein sequence against sequences of all the species present in ensembl . it was observed that no orthologs of daoa were found in any other species through ncbi and ensembl . daoa may be absent in all other sequenced species except hominoids and great apes . in silico analysis revealed that the gene was only found in humans , chimpanzees , gorillas , orangutans , and crab - eating macaques . the phylogenetic tree was reconciled with primate phylogeny.58 it is hypothesized that daoa is a primate - specific gene according to the analysis of current available biological sequences data in biological databases ( ncbi , ensembl , ucsc ) . it is also hypothesized that may be this gene evolved approximately 35 million years ago before the divergence of the new - world monkey from the old - world monkey clade . the old - world monkeys and new - world monkeys shared derived features , and the primate s brain is more developed than those of other species . although the significance of dao regulation in the function of neurophysiology has been revealed , the authenticated exact mechanism remains vague . the molecular characterization of daoa showed that it was an endogenous modulator of the activity of dao . this study reveals the binding interactions of daoa isoforms with ligand . by employing computational and in silico approaches , it is demonstrated that daoa isoforms have conserved regions ( figure 8) and conserved c - terminal region . this work also identified common binding residues that interact with selected ligand and normalize the expression of daoa . the bioinformatics analysis indicated that binding residues within c - terminal of daoa are essential for ligand , instead of n - terminal . the manipulation in the expression and activity of daoa can be a fruitful option for the treatment of sz . the results of the docking analysis , by selected ligand , could be used for designing novel chemical compounds and synthetic peptides to reduce the overexpression of daoa that leads to sz . the level of d - serine plays a crucial role in neurotransmission of glutamate nmda . the level of d - serine must be maintained to avoid over- or underexpression of nmda receptor function within a relatively narrow range for avoiding cns function impairment . in schizophrenic patients , the d - serine level in serum is 17.7% lower than normal controls , and the ratio of d - serine in serum to total serine is 28.4% lower . in conclusion , the selected ligand molecule is efficacious in the treatment of sz for targeting daoa . the insilico analysis of daoa has higher probability and efficacy on the basis of binding energy and other used parameters . it is also hypothesized that daoa is a primate - specific gene , and only exists in great apes . the future analysis of this gene will be helpful in understanding the behavior of homo sapiens . further studies and synthesis of novel compounds considering these findings can expect similar response rates and can cure the effects of sz .
schizophrenia ( sz ) , a chronic mental and heritable disorder characterized by neurophysiological impairment and neuropsychological abnormalities , is strongly associated with d - amino acid oxidase activator ( daoa , g72 ) . 
 research studies emphasized that overexpression of daoa may be responsible for improper functioning of neurotransmitters , resulting in neurological disorders like sz . in the present study , 
 a hybrid approach of comparative modeling and molecular docking followed by inhibitor identification and structure modeling was employed . 
 screening was performed by two - dimensional similarity search against selected inhibitor , keeping in view the physiochemical properties of the inhibitor . here , we report an inhibitor compound which showed maximum binding affinity against four selected isoforms of daoa . 
 docking studies revealed that glu-53 , thr-54 , lys-58 , val-85 , ser-86 , tyr-87 , leu-88 , glu-90 , leu-95 , val-98 , ser-100 , glu-112 , tyr-116 , lys-120 , asp-121 , and arg-122 are critical residues for receptor  ligand interaction . the c - terminal of selected isoforms is conserved , and binding was observed on the conserved region of isoforms . 
 we propose that selected inhibitor might be more potent on the basis of binding energy values . 
 further analysis of this inhibitor through site - directed mutagenesis could be helpful for exploring the details of ligand - binding pockets . 
 overall , the findings of this study may be helpful in designing novel therapeutic targets to cure sz .