Computational antisense oligo prediction with a neural network model

MOTIVATION: The expression of a gene can be selectively inhibited by antisense oligonucleotides (AOs) targeting the mRNA. However, if the target site in the mRNA is picked randomly, typically 20% or less of the AOs are effective inhibitors in vivo. The sequence properties that make an AO effective are not well understood, thus many AOs need to be tested to find good inhibitors, which is time consuming and costly. So far computational models have been based exclusively on RNA structure prediction or motif searches while ignoring information from other aspects of AO design into the model. RESULTS: We present a computational model for AO prediction based on a neural network approach using a broad range of input parameters. Collecting sequence and efficacy data from AO scanning experiments in the literature generated a database of 490 AO molecules. Using a set of derived parameters based on AO sequence properties we trained a neural network model. The best model, an ensemble of 10 networks, gave an overall correlation coefficient of 0.30 (p=10(-8)). This model can predict effective AOs (>50% inhibition of gene expression) with a success rate of 92%. Using these thresholds the model predicts on average 12 effective AOs per 1000 base pairs, making it a stringent yet practical method for AO prediction.     

Enzymatic activity of the Arabidopsis sulfurtransferase resides in the C-terminal domain but is boosted by the N-terminal domain and the linker peptide in the full-length enzyme

Sulfurtransferases/rhodaneses are a group of enzymes widely distributed in plants, animals, and bacteria that catalyze the transfer of sulfur from a donor molecule to a thiophilic acceptor substrate. Sulfurtransferases (STs) consist of two globular domains of nearly identical size and conformation connected by a short linker sequence. In plant STs this linker sequence is exceptionally longer than in sequences from other species. The Arabidopsis ST1 protein (AJ131404) contains five cysteine residues: one residue is universally conserved in all STs and considered to be catalytically essential; a second one, closely located in the primary sequence, is conserved only in sequences from eukaryotic species. Of the remaining three cysteine residues two are conserved in the so far known plant STs and one is unique to the Arabidopsis ST1. The aim of our study was to investigate the role of the two-domain structure, of the unique plant linker sequence and of each cysteine residue. The N- and C-terminal domains of the Arabidopsis ST1, the full-length protein with a shortened linker sequence and several point-mutated proteins were overexpressed in E. coli, purified and used for enzyme activity measurements. The C-terminal domain itself displayed ST activity which could be increased by adding the separately prepared N-terminal domain. The activity of an ST1 derivative with a shortened linker sequence was reduced by more than 60% of the wild-type activity, probably because of a drastically reduced protein stability. The replacement of each cysteine residue resulted in mutant forms which differed significantly in their stability, in the specific ST activities, and in their kinetic parameters which were determined for 3-mercaptopyruvate as well as thiosulfate as sulfur substrates: mutation of the putative active site cysteine (C332) essentially abolished activity; for C339 a crucial role at least for the turnover of thiosulfate could be identified.     

The domain structure of the dihydrolipoyl transacetylase component of the pyruvate dehydrogenase complex from Azotobacter vinelandii

Limited proteolysis with trypsin has been used to study the domain structure of the dihydrolipoyltransacetylase (E2) component of the pyruvate dehydrogenase complex of Azotobacter vinelandii. Two stable end products were obtained and identified as the N-terminal lipoyl domain and the C-terminal catalytic domain. By performing proteolysis of E2, which was covalently attached via its lipoyl groups to an activated thiol-Sepharose matrix, a separation was obtained between the catalytic domain and the covalently attached lipoyl domain. The latter was removed from the column after reduction of the S-S bond and purified by ultrafiltration. The lipoyl domain is monomeric with a mass of 32.6 kDa. It is an elongated structure with f/fo = 1.62. Circulair dichroic studies indicates little secondary structure. The catalytic domain is polymeric with S20.w = 17 S and mass = 530 kDa. It is a compact structure with f/fo = 1.24 and shows 40% of the secondary structure of E2. The cubic structure of the native E2 is retained by this fragment as observed by electron microscopy. Ultracentrifugation in 6 M guanidine hydrochloride in the presence of 2 mM dithiothreitol yields a mass of 15.8 kDa. An N-terminal sequence of 36 amino acids is homologous with residues 370-406 of Escherichia coli E2. The catalytic domain possesses the catalytic site, but in contrast to the E. coli subunit binding domain the pyruvate dehydrogenase (E1) and lipoamide dehydrogenase (E3) binding sites are lost during proteolysis. From comparison with the E. coli E2 sequence a model is presented in which the several functions, such as lipoyl domain, the E3 binding site, the catalytic site, the E2/E2 interaction sites, and the E1 binding site, are indicated.     

Functional study of rat 5-HT2A receptors using antisense oligonucleotides

We studied the effects in rats of a 6-day intracerebroventricular (i.c.v) infusion of four different end-capped phosphorothioate-modified antisense oligonucleotides (AOs), specifically targeting different regions of the 5-hydroxytryptamine2A (5-HT2A) receptor mRNA, on central 5-HT2A receptor expression and 5-HT2A receptor-mediated behaviours. Only one of the AOs (sequence 4), directed against the 5'-untranslated region (from + 557 to + 577), specifically affected central 5-HT2A receptor expression and receptor-mediated behaviour. This AO (sequence 4) reduced binding of the 5-HT2A agonist 1-(2,5-dimethoxy-4-[125I]iodophenyl)-2-aminopropane ([125I]DOI) up to 25% in cortical areas, as measured by quantitative autoradiography. Cortical binding of the antagonist [3H]ketanserin was not affected. As the specific AO treatment presumably affects the synthesis of new receptor, we hypothesize that this newly synthesized receptor represents the major part of the functionally active, G protein coupled receptor. A 5-day infusion of AO (sequence 4) resulted in profound inhibition of the head-twitch response (HTR) to 1-(2,5-dimethoxy-4-methylphenyl)-2-aminopropane (DOM). In contrast, treatment with vehicle, sense oligonucleotides (SOs) and other AOs (sequences 1, 2 and 3) caused an increased DOM-induced HTR as well as a spontaneous HTR. The latter was abolished by treatment with the 5-HT2 receptor antagonist, ritanserin. Systematic investigation of the surgical and infusion procedures revealed that the enhanced HTR already appeared following drilling of the skull. This wounding can probably damage the blood-brain barrier and cause a stress-induced increase in serotonergic transmission. AO (sequence 4) treatment also abolished the spontaneous HTR. AO (sequence 4) treatment allowed the identification of specific central 5-HT2A receptor-mediated behaviours in the complex serotonergic syndrome induced by tryptamine in rats. Only bilateral convulsions and body tremors were significantly inhibited. The backward locomotion, hunched back and Straub tail were not affected, nor was cyanosis, an index of vasoconstriction induced by peripheral 5-HT2A receptor activation. Labelling of central 5-HT2C receptors by [3H]mesulergine, and 5-HT2C receptor-mediated anxiety were not attenuated by AO or SO treatment. Rats treated with AO (sequence 4) showed increased locomotor activity and a strong reactivity towards touching. We hypothesize that the down-regulation of functional 5-HT2A receptors may shift the balance between various 5-HT receptor subtypes. Our analysis of the behavioural consequences of AO treatment and the use of different AOs and SOs has shown that specific receptor-mediated behaviour can be identified.     

Insilco analysis of functionally important residues in folate receptors

Lack of crystal structure data of folate binding proteins has left so many questions unanswered (for example, important residues in active site, binding domain, important amino acid residues involved in interactions between ligand and receptor). With sequence alignment and PROSITE motif identification, we attempted to answer evolutionarily significant residues that are of functional importance for ligand binding and that form catalytic sites. We have analyzed 46 different FRs and FBP sequences of various organisms obtained from Genbank. Multiple sequence alignment identified 44 highly conserved identical amino acid residues with 10 cysteine residues and 12 motifs including ECSPNLGPW (which might help in the structural stability of FR).     

VLDL-受体的配体结合结构域结构分析

极低密度脂蛋白受体(VLDL-R)的配体结合域具有8个富含半胱氨酸的配体结合重复序列(ligand-binding repeats,LBR),被认为是与配体结合的部位。该受体与含7个类似重复序列的低密度脂蛋白受体(LDL-R)的配体结合特性明显不同。为了明确VLDL-R中8个LBR在配体结合中的作用并探讨结合位点的结构,本研究采用计算机辅助蛋白质结构预测方法,在二级结构分析的基础上,通过同源建模方法预测受体N-端328个氨基酸的配体结合域空间结构,结果显示该区域呈现弧形口袋样结构,其中前3个LBR结构紧凑,呈棒状,负电荷相对集中,推测这一特征结构是配体结合的重要结构基础,结构分析同时表明在LBR5与LBR6之间连接区的弹性结构可以赋予结合位点一定的伸缩性,利于其与不同配体的结合。本研究预测结果首次提出了VLDL-R配体结构域结构及结合位点的结构特征,并与已有的实验结果一致。     

Architecture of the U5 small nuclear RNA.

We have used comparative sequence analysis and deletion analysis to examine the secondary structure of the U5 small nuclear RNA (snRNA), an essential component of the pre-mRNA splicing apparatus. The secondary structure of Saccharomyces cerevisiae U5 snRNA was studied in detail, while sequences from six other fungal species were included in the phylogenetic analysis. Our results indicate that fungal U5 snRNAs, like their counterparts from other taxa, can be folded into a secondary structure characterized by a highly conserved stem-loop (stem-loop 1) that is flanked by a moderately conserved internal loop (internal loop 1). In addition, several of the fungal U5 snRNAs include a novel stem-loop structure (ca. 30 nucleotides) that is adjacent to stem-loop 1. By deletion analysis of the S. cerevisiae snRNA, we have demonstrated that the minimal U5 snRNA that can complement the lethal phenotype of a U5 gene disruption consists of (i) stem-loop 1, (ii) internal loop 1, (iii) a stem-closing internal loop 1, and (iv) the conserved Sm protein binding site. Remarkably, all essential, U5-specific primary sequence elements are encoded by a 39-nucleotide domain consisting of stem-loop 1 and internal loop 1. This domain must, therefore, contain all U5-specific sequences that are essential for splicing activity, including binding sites for U5-specific proteins.     

Structural and functional conservation profiles of novel cathepsin L-like proteins identified in the Drosophila melanogaster genome

Cathepsin L is a cysteine protease which degrades connective tissue proteins including collagen, elastin, and fibronectin. In this study, five well-characterized cathepsin L proteins from different arthropods were used as query sequences for the Drosophila genome database. The search yielded 10 cathepsin L-like sequences, of which eight putatively represent novel cathepsin L-like proteins. To understand the phylogenetic relationship among these cathepsin L-like proteins, a phylogenetic tree was constructed based on their sequences. In addition, models of the tertiary structures of cathepsin L were constructed using homology modeling methods and subjected to molecular dynamics simulations to obtain reasonable structure to understand its dynamical behavior. Our findings demonstrate that all of the potential Drosophila cathepsin L-like proteins contain at least one cathepsin propeptide inhibitor domain. Multiple sequence alignment and homology models clearly highlight the conservation of active site residues, disulfide bonds, and amino acid residues critical for inhibitor binding. Furthermore, comparative modeling indicates that the sequence/structure/function profiles and active site architectures are conserved.     

The dihydrolipoyltransacetylase component of the pyruvate dehydrogenase complex from Azotobacter vinelandii. Molecular cloning and sequence analysis

The gene encoding the dihydrolipoyltransacetylase component (E2) of the pyruvate dehydrogenase complex from Azotobacter vinelandii has been cloned in Escherichia coli. A plasmid containing a 2.8-kbp insert of A. vinelandii chromosomal DNA was obtained and its nucleotide sequence determined. The gene comprises 1911 base pairs, 637 codons excluding the initiation codon GUG and stop codon UGA. It is preceded by the gene encoding the pyruvate dehydrogenase component (E1) of pyruvate dehydrogenase complex and by an intercistronic region of 11 base pairs containing a good ribosome binding site. The gene is followed downstream by a strong terminating sequence. The relative molecular mass (64913), amino acid composition and N-terminal sequence are in good agreement with information obtained from studies on the purified enzyme. Approximately the first half of the gene codes for the lipoyl domain. Three very homologous sequences are present, which are translated in three almost identical units, alternated with non-homologous regions which are very rich in alanyl and prolyl residues. The N-terminus of the catalytic domain is sited at residue 381. Between the lipoyl domain and the catalytic domain, a region of about 50 residues is found containing many charged amino acid residues. This region is characterized as a hinge region and is involved in the binding of the pyruvate dehydrogenase and lipoamide dehydrogenase components. The homology with the dihydrolipoyltransacetylase from E. coli is high: 50% amino acid residues are identical.     

Molecular characterization of the amplified aldehyde oxidase from insecticide resistant Culex quinquefasciatus.

Primary structural information including the complete nucleotide sequence of the first insect aldehyde oxidase (AO) was obtained from the common house mosquito Culex quinquefasciatus (Say) through cloning and sequencing of both genomic DNA and cDNA. The deduced amino-acid sequence encodes a 150-kDa protein of 1266 amino-acid residues, which is consistent with the expected monomeric subunit size of AO. The Culex AO sequence contains a molybdopterin cofactor binding domain and two iron-sulfur centres. A comparison of the partial sequences of AO from insecticide resistant and susceptible strains of C. quinquefasciatus shows two distinct alleles of this enzyme, one of which is amplified in the insecticide resistant strain on a 30-kb DNA amplicon alongside two resistance-associated esterases. The amplified AO gene results in elevated AO activity in all life stages, but activity is highest in 3rd instar larvae. The elevated enzyme can be seen as a separate band on polyacrylamide gel electrophoresis. The role of AO in xenobiotic oxidation in mammals and the partial inhibition of elevated AO activity by a range of insecticides in Culex, suggest that this AO may play a role in insecticide resistance.     

The cyanide-resistant alternative oxidases from the fungi Pichia stipitis and Neurospora crassa are monomeric and lack regulatory features of the plant enzyme

Both plant and fungal mitochondria have cyanide-resistant alternative oxidases that use reductant from the mitochondrial ubiquinone pool to reduce oxygen to water in a reaction that conserves no energy for ATP synthesis. The dimeric plant alternative oxidase is relatively inactive when its subunits are linked by a disulfide bond. When this bond is reduced, the enzyme can then be stimulated by its activators, alpha-keto acids. A Cys in the N-terminal section of the protein is responsible for both of these features. We examined the alternative oxidases in mitochondria isolated from two fungi Neurospora crassa and Pichia stipitis for dimeric structure, ability to form an intermolecular disulfide, and sensitivity to alpha-keto acids. Neither of the two fungal alternative oxidases could be covalently linked by diamide, which induces disulfide bond formation between nearby Cys residues, nor could they be cross-linked by a Lys-specific reagent or glutaraldehyde at concentrations which cross-link the plant alternative oxidase dimer completely. Alternative oxidase activity in fungal mitochondria was not stimulated by the alpha-keto acids pyruvate and glyoxylate. Pyruvate did stimulate activity when succinate was the respiratory substrate, but this was not a direct effect on the alternative oxidase. In contrast, added GMP was a strong activator of fungal alternative oxidase activity. Analysis of plant and fungal alternative oxidase protein sequences revealed a unique domain of about 40 amino acids surrounding the regulatory Cys in the plant sequences that is not present in the fungal sequences. This domain may be where dimerization of the plant enzymes occurs. In contrast to plant enzymes, the fungal alternative oxidases studied here are monomeric and their activities are independent of alpha-keto acids.     

Possible structure and function of the extra C-terminal sequence of pyruvate kinase from Bacillus stearothermophilus

The pyruvate kinases from Genus Bacillus and a few other bacteria have an extra C-terminal sequence with a phosphoenolpyruvate binding motif composed of about 110 amino acids. To elucidate the possible structure and function of this sequence, the enzyme lacking the sequence was prepared and characterized. The N-terminal sequences of the peptides, which were found only in the lysylendopeptidase digest of the wild enzyme and not in that of the truncated enzyme, were determined. All the determined sequences were found in the extra C-terminal sequence deduced from the DNA sequence. The truncated enzyme showed decreased affinity for phosphoenolpyruvate and the allosteric effector ribose 5-phosphate, and had a reduced thermostability. Other properties, such as tetrameric structure, specific activity, and allosteric characteristics were unchanged. A comparison of the CD spectra of the truncated enzyme and the recombinant enzyme indicated that the structure of the C-terminal sequence should be rich in beta-sheet. These findings suggest that the sequence actually exists and that it may form a steady domain interacting with the A-domain and C-domain, which are the catalytic domain and allosteric effector binding domain, respectively.     

Multiplicity of carbohydrate-binding sites in beta-prism fold lectins: occurrence and possible evolutionary implications

The beta-prism II fold lectins of known structure, all from monocots, invariably have three carbohydrate-binding sites in each subunit/domain. Until recently, beta-prism I fold lectins of known structure were all from dicots and they exhibited one carbohydrate-binding site per subunit/domain. However, the recently determined structure of the beta-prism fold I lectin from banana, a monocot, has two very similar carbohydrate-binding sites. This prompted a detailed analysis of all the sequences appropriate for two-lectin folds and which carry one or more relevant carbohydrate-binding motifs. The very recent observation of a beta-prism I fold lectin, griffthsin, with three binding sites in each domain further confirmed the need for such an analysis. The analysis demonstrates substantial diversity in the number of binding sites unrelated to the taxonomical position of the plant source. However, the number of binding sites and the symmetry within the sequence exhibit reasonable correlation. The distribution of the two families of beta-prism fold lectins among plants and the number of binding sites in them, appear to suggest that both of them arose through successive gene duplication, fusion and divergent evolution of the same primitive carbohydrate-binding motif involving a Greek key. Analysis with sequences in individual Greek keys as independent units lends further support to this conclusion.It would seem that the preponderance of three carbohydrate-binding sites per domain in monocot lectins, particularly those with the beta-prism II fold, is related to the role of plant lectins in defence.     

Cloning, in silico characterization and interaction of cysteine protease and cystatin for establishing their role in early blight disease in tomato

Cysteine protease (CP) and Cysteine protease inhibitor (CPI) or cystatin constitute a critical point in programmed cell death (PCD), a basic biological phenomenon which takes place in the plants, when they are exposed to varying biotic and abiotic stresses. In the present study we isolated and cloned cDNAs encoding cysteine protease and cystatin from early blight infected tomato plants. Using computational biology tools the sequence-structure-function relationships for the tomato cystatin and cysteine protease were elucidated. Interaction between the cystatin and cysteine protease of host and pathogen is higher as compared to interaction shown by cystatin and cysteine protease within the host. The interaction energy of (a)tomato cystatin—tomato cysteine protease, (b)tomato cystatin—fungal cysteine protease and (c)tomato cysteine protease—fungal cystatin are ?319.33 Kcal/mol, ?504.71 Kcal/mol and ?373.731 Kcal/mol respectively. Comparative protein sequence analysis with different plant cystatins and cysteine protease were also done with the sequences of cystatin and cysteine protease isolated from tomato. Structures for all the cystatin and cysteine protease were modeled along with their interactions with fungal cystatin and cysteine protease in order to explore the structural variability and its manifestation at the functional level. This helped to relate the already known functions of these proteins with their sequences as well as the predicted structures. This also served to better understand the CP-CPI interaction operational in developing this protein family and its implication in plant defense during fungal pathogenesis in tomato plants.     

Solution structure of the PABC domain from wheat poly (A)-binding protein: an insight into RNA metabolic and translational control in plants

In animals, the PABC domain from poly (A)-binding protein recruits proteins containing a specific interacting motif (PAM-2) to the mRNP complex. These proteins include Paip1, Paip2, and eukaryotic release factor 3 (eRF3), all of which regulate PABP function in translation. The following reports the solution structure of PABC from Triticum avestium (wheat) poly (A)-binding protein determined by NMR spectroscopy. Wheat PABC (wPABC) is an alpha-helical protein domain, which displays a fold highly similar to the human PABC domain and contains a PAM-2 peptide binding site. Through a bioinformatics search, several plant proteins containing a PAM-2 site were identified including the early response to dehydration protein (ERD-15), which was previously shown to regulate PABP-dependent translation. The plant PAM-2 proteins contain a variety of conserved sequences including a PABP-interacting 1 motif (PAM-1), RNA binding domains, an SMR endonuclease domain, and a poly (A)-nuclease regulatory domain, all of which suggest a function in either translation or mRNA metabolism. The proteins identified are well conserved throughout plant species but have no sequence homologues in metazoans. We show that wPABC binds to the plant PAM-2 motif with high affinity through a conserved mechanism. Overall, our results suggest that plant species have evolved a distinct regulatory mechanism involving novel PABP binding partners.     

Trypsin-like proteins of the fungi as possible markers of pathogenicity

Sequences of peptidases with conserved motifs around the active site residues that are characteristic of trypsins (similar to trypsin peptidases, STP) were obtained from publicly-available fungal genomes and related databases. Among the 75 fungal genomes, 29 species of parasitic Ascomycota contained genes encoding STP and their homologs. Searches of non-redundant protein sequences, patented protein sequences, and expressed sequence tags resulted in another 18 STP sequences in 10 fungal species from Ascomycota, Basidiomycota, and Zygomycota. A comparison of fungi species containing STP sequences revealed that almost all are pathogens of plants, animals or fungi. A comparison of the primary structure of homologous proteins, including the residues responsible for substrate binding and specificity of the enzyme, revealed three groups of homologous sequences, all presumably from S1 family: trypsin-like peptidases, chymotrypsin-like peptidases and serine peptidases with unknown substrate specificity. Homologs that are presumably functionally inactive were predicted in all groups. The results in general support the hypothesis that the expression of trypsin-like peptidases in fungi represents a marker of fungal phytopathogenicity. A phylogenetic tree was constructed using peptidase and homolog amino acid sequences, demonstrating that all have noticeable differences and almost immediately deviate from the common root. Therefore, we conclude that the changes that occurred in STP of pathogenic fungi in the course of evolution represent specific adaptations to proteins of their respective hosts, and mutations in peptidase genes are important components of life-style changes and taxonomic divergence.     

Characterization of a cellulosome dockerin domain from the anaerobic fungus Piromyces equi

The recycling of photosynthetically fixed carbon in plant cell walls is a key microbial process. In anaerobes, the degradation is carried out by a high molecular weight multifunctional complex termed the cellulosome. This consists of a number of independent enzyme components, each of which contains a conserved dockerin domain, which functions to bind the enzyme to a cohesin domain within the protein scaffoldin protein. Here we describe the first three-dimensional structure of a fungal dockerin, the N-terminal dockerin of Cel45A from the anaerobic fungus Piromyces equi. The structure contains a novel fold of 42 residues. The ligand binding site consists of residues Trp 35, Tyr 8 and Asp 23, which are conserved in all fungal dockerins. The binding site is on the opposite side of the N- and C-termini of the molecule, implying that tandem dockerin domains, seen in the majority of anaerobic fungal plant cell wall degrading enzymes, could present multiple simultaneous binding sites and, therefore, permit tailoring of binding to catalytic demands.     

Crystal structure of pyruvate kinase from Geobacillus stearothermophilus

The pyruvate kinase (PK) from a moderate thermophile, Geobacillus stearothermophilus, is an allosteric enzyme activated by AMP and ribose 5-phosphate but not fructose 1, 6-bisphosphate (FBP), which is a common activator of PKs. It has an extra C-terminal sequence (ECTS), which contains a highly conserved phosphoenolpyruvate (PEP) binding motif, but its function and structure remain unclear. To elucidate the structural characteristics of the effector-binding site and the ECTS, the crystal structure of the C9S/C268S mutant of the enzyme was determined at 2.4 A resolution. The crystal belonged to space group P6(2)22, with unit cell parameters a, b = 145.97 A, c = 118.03 A. The enzyme was a homotetramer and its overall domain structure was similar to the previously solved structures except that the ECTS formed a new domain (C' domain). The structure of the C' domain closely resembled that of the PEP binding domain of maize pyruvate phosphate dikinase. A sulphate ion was found in a pocket in the effector-binding C domain. This site corresponds to the 6-phosphate group-binding site in yeast PK bound FBP and seems to be the effector-binding site. Through comparison of the structure of the putative effector-binding site to that of the FBP binding site of the yeast enzyme, the structural basis of the effector specificity of the G. stearothermophilus PK is discussed.     

Study on Folate Binding Domain of Dihydrofolate Reductase in Different Plant species and Human beings

Data base (NCBI and TIGR) searches are made to retrieve protein sequences of different plant species namely Medicago truncatula, Pisum sativum, Ricinus communis, Arabidopsis thaliana, Vitis vinifera, Glycine max, Daucus carota, Oryza sativa Japonica Group, Arabidopsis lyrata subsp. lyrata, Brachypodium distachyon, Oryza sativa Indica Group, Zea mays and careful alignment of derived sequences shows 95% or higher identity. Similarly, DHFR sequence of human being is also retrieved from NCBI. A phylogenetic tree is constructed from different plant and human DHFR domain using the Neighbour – Joining method in MEGA 5.05. Conservation score is performed by using PARALINE. Result suggests that folate binding domain of dihydrofolare reductase is conserved (score 8.06) and excepting some minor variations the basic structure of the domain in both plant species and human being is rather similar. Human DHFR domain contains PEKN sequence near active site, though proline is common for all the selected organisms but the other sequences are different in plants. The plant domain is always associated with TS (Thymidylate synthase). Plant based system is predicted to be an effective model for assessment of MTX (Methotrexate) and other antifolate drugs.