RESUMO
Restriction endonucleases protect bacterial cells against bacteriophage infection by cleaving the incoming foreign DNA into fragments. In presence of Mg2+ ions, EcoRV is able to cleave the DNA but not in presence of Ca2+ , although the protein binds to DNA in presence of both metal ions. We make an attempt to understand this difference using conformational thermodynamics. We calculate the changes in conformational free energy and entropy of conformational degrees of freedom, like DNA base pair steps and dihedral angles of protein residues in Mg2+ (A)-EcoRV-DNA complex compared to Ca2+ (S)-EcoRV-DNA complex using all-atom molecular dynamics (MD) trajectories of the complexes. We find that despite conformational stability and order in both complexes, the individual degrees of freedom behave differently in the presence of two different metal ions. The base pairs in cleavage region are highly disordered in Ca2+ (S)-EcoRV-DNA compared to Mg2+ (A)-EcoRV-DNA. One of the acidic residues ASP90, coordinating to the metal ion in the vicinity of the cleavage site, is conformationally destabilized and disordered, while basic residue LYS92 gets conformational stability and order in Ca2+ (S) bound complex than in Mg2+ (A) bound complex. The enhanced fluctuations hinder placement of the metal ion in the vicinity of the scissile phosphate of DNA. Similar loss of conformational stability and order in the cleavage region is observed by the replacement of the metal ion. Considering the placement of the metal ion near scissile phosphate as requirement for cleavage action, our results suggest that the changes in conformational stability and order of the base pair steps and the protein residues lead to cofactor sensitivity of the enzyme. Our method based on fluctuations of microscopic conformational variables can be applied to understand enzyme activities in other protein-DNA systems.
Assuntos
Clivagem do DNA , DNA/química , DNA/metabolismo , Desoxirribonucleases de Sítio Específico do Tipo II/metabolismo , Magnésio/metabolismo , Manganês/metabolismo , Sítios de Ligação , Catálise , Desoxirribonucleases de Sítio Específico do Tipo II/química , Desoxirribonucleases de Sítio Específico do Tipo II/genética , Magnésio/química , Manganês/química , Modelos Moleculares , Conformação Proteica , Especificidade por Substrato , TermodinâmicaRESUMO
DNA is an essential target for the treatment of various pathologies, especially cancer. Hence targeting DNA double helix for alteration of its function has been attempted by several ways. Drug-DNA intercalation, one such biophysical process, could not be studied extensively as this requires significant deformation of the receptor DNA. Here we report thorough theoretical investigation of intercalation process in daunomycin-DNA interaction, by performing molecular dynamics simulations of the drug-DNA complexes for various DNA sequences, followed by Free-energy analysis and density functional theory (DFT) based studies to understand the binding preference. The classical energy based analyses indicate that the drug prefers to bind to TC/GA sequence over others. The DFT based energies of supra-molecular complexes are always contaminated with basis set superposition error (BSSE), which can be corrected by counterpoise method. This method is quite effective for systems containing two molecular fragments but is not appropriate for studying interaction between two base pair fragments and the drug intercalated between them. We have adopted an extension of the counterpoise method for BSSE corrected interaction energy calculation. These interaction energies, along with the energy penalty due to un-stacking of the base pairs, also indicate TC/GA sequence is the most preferred sequence for binding.
Assuntos
Antibióticos Antineoplásicos/farmacologia , DNA/química , Daunorrubicina/farmacologia , Substâncias Intercalantes/farmacologia , Teoria da Densidade Funcional , Simulação de Dinâmica Molecular , Conformação de Ácido Nucleico/efeitos dos fármacos , TermodinâmicaRESUMO
Electrorheological (ER) fluids are known to exhibit enhanced viscous effects under an electric field stimulus. The present article reports the hitherto unreported phenomenon of greatly enhanced thermal conductivity in such electro-active colloidal dispersions in the presence of an externally applied electric field. Typical ER fluids are synthesized employing dielectric fluids and nanoparticles and experiments are performed employing an in-house designed setup. Greatly augmented thermal conductivity under a field's influence was observed. Enhanced thermal conduction along the fibril structures under the field effect is theorized as the crux of the mechanism. The formation of fibril structures has also been experimentally verified employing microscopy. Based on classical models for ER fluids, a mathematical formalism has been developed to predict the propensity of chain formation and statistically feasible chain dynamics at given Mason numbers. Further, a thermal resistance network model is employed to computationally predict the enhanced thermal conduction across the fibrillary colloid microstructure. Good agreement between the mathematical model and the experimental observations is achieved. The domineering role of thermal conductivity over relative permittivity has been shown by proposing a modified Hashin-Shtrikman (HS) formalism. The findings have implications towards better physical understanding and design of ER fluids from both 'smart' viscoelastic as well as thermally active materials points of view.
RESUMO
Denaturation of the DNA double helix inside the cell is essential for cellular processes such as replication and transcription for the growth of the cells. However, the growth of unwanted cells, which are responsible for cancerous kind of disease, is one of the biggest challenges of modern therapeutics. DNA cross-linking agents may kill cancer cells by damaging their DNA and stopping them from dividing. In the present study, we have carried out steered molecular dynamics simulations to study the effects of rupture and unzipping forces on the stability of dsDNA in the absence and presence of covalently bonded drugs. We have found that the stability of dsDNA increases strongly in the presence of covalently bonded drugs. The microscopic study of disruption of hydrogen-bonds associated with base-pairs of the dsDNA and the study of the variation of stacking overlap parameters gives evidence of symmetry during the rupture and asymmetry in the unzip event. The significance of the mechanism of force-induced melting study of the dsDNA in the absence and presence of antitumor drugs might have a biological relevance as it provides a pathway to open the double helix in a specific position and may help for the pharmaceutical design of drugs.Communicated by Ramaswamy H. Sarma.
Assuntos
Antineoplásicos , Simulação de Dinâmica Molecular , Desnaturação de Ácido Nucleico , DNA/química , Pareamento de Bases , Antineoplásicos/farmacologia , Antineoplásicos/química , Conformação de Ácido NucleicoRESUMO
Lamins maintain the shape and rigidity of the nucleus in the form of a proteinaceous scaffold underneath the inner nuclear membrane (INM) and provide anchorage to chromatin and other nuclear proteins. Mutations in the human LMNA gene encoding lamin A/C cause about 16 different diseases with distinct phenotypes collectively termed as laminopathies which affect primarily the muscle tissues as well as adipose tissues, neuromuscular junctions and multiple other organs in progeroid syndromes. Lamins contain several domains of which Ig-fold is one of the well characterized and structured domains that harbours many mutations leading to deleterious interactions with other nuclear proteins. In this work, we have elucidated the effects of 3 such mutations namely R453W, W498C and W498R on the dynamics and flexibility of the Ig-fold domain and the consequent effect on the assembly into lamina by live cell imaging, fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulations. From our simulation studies, we concluded that R453W exhibits the highest fluctuation at the residues 475 and 525 in the Ig fold domain compared to the wild type and other mutants. This resulted in pronounced random self-association which could be corroborated by lower diffusivity values obtained from FCS. This is the first report where such an alteration in the full length has been documented by gross changes in diffusional properties as a sequel to a mutation in the Ig fold domain.
RESUMO
Mycobacterium tuberculosis resistance to commercially available drugs is increasing day by day. To address this issue, various strategies were planned and are being implemented. However, there is a need for new drugs and rapid diagnostic methods. For this endeavour, in this paper, we present the synthesis of acetylene containing 2-(2-hydrazinyl) thiazole derivatives and in vitro evaluation against the H37Rv strain of Mycobacterium tuberculosis. Among the developed 26 acetylene containing 2-(2-hydrazinyl) thiazole derivatives, eight compounds inhibited the growth of Mycobacterium tuberculosis with MIC values ranging from 100 µg ml-1 to 50 µg ml-1. The parent acetylene containing thiosemicarbazones showed promising antimycobacterial activity by inhibiting up to 75% of the Mycobacterium at 50 µg ml-1. In addition, in silico studies were employed to understand the binding mode of all the novel acetylene-containing derivatives against the KasA protein of the Mycobacterium. Interestingly, the KasA protein interactions with the compounds were similar to the interactions of KasA protein with thiolactomycin and rifampicin. Cytotoxicity study results indicate that the compounds tested are non-toxic to human embryonic kidney cells.
RESUMO
Here, we study microscopic mechanism of complex formation between Ca2+-bound calmodulin (holoCaM) and Orai1 that regulates Ca2+-dependent inactivation process in eukaryotic cells. We compute conformational thermodynamic changes in holoCaM with respect to complex of Orai1 bound to C-terminal domain of holoCaM using histograms of dihedral angles of the proteins over trajectories from molecular dynamics simulations. Our analysis shows that the N-terminal domain residues L4, T5, Q41, N42, T44 and E67 of holoCaM get destabilized and disordered due to Orai1 binding to C-terminal domain of calmodulin affect the N-terminal domain residues. Among these residues, polar T44, having maximum destabilization and disorder via backbone fluctuations, shows the largest change in solvent exposure. This suggests that N-terminal domain is allosterically regulated via T44 by the binding of Orai1 to the C-terminal domain.
Assuntos
Calmodulina/química , Simulação de Dinâmica Molecular , Proteína ORAI1/química , Conformação Proteica , Termodinâmica , Regulação Alostérica , Calmodulina/metabolismo , Simulação de Acoplamento Molecular , Proteína ORAI1/metabolismo , Ligação Proteica , Domínios e Motivos de Interação entre Proteínas , Relação Estrutura-AtividadeRESUMO
Molecular understanding of bio-macromolecular binding is a challenging task due to large sizes of the molecules and presence of variety of interactions. Here, we study the molecular mechanism of calmodulin (CaM) binding to Orai1 that regulates Ca2+-dependent inactivation process in eukaryotic cells. Although experimental observations indicate that Orai1 binds to the C-terminal of Ca2+-loaded CaM, it is not decisive if N-domain of CaM interacts with Orai1. We address the issue of interaction of different domains of CaM with Orai1 using conformational thermodynamic changes, computed from histograms of dihedral angles over simulated trajectories of CaM, CaM-binding domain of Orai1 and complexes of CaM with Orai1. The changes for all residues of both C and N terminal domains of CaM upon Orai1 binding are compared. Our analysis shows that Orai1binds to both C-terminal and N-terminal domains of CaM, indicating 1:2 stoichiometry. The Orai1 binding to N-terminal domain of CaM is less stable than that to the C-terminal domain. The binding residues are primarily hydrophobic. These observations are in qualitative agreement to the experiments. The conformational thermodynamic changes thus provide a useful computational tool to provide atomic details of interactions in bio-macromolecular binding.
Assuntos
Calmodulina/química , Complexos Multiproteicos/química , Proteína ORAI1/química , Termodinâmica , Sequência de Aminoácidos/genética , Cálcio/metabolismo , Calmodulina/genética , Humanos , Interações Hidrofóbicas e Hidrofílicas , Simulação de Dinâmica Molecular , Complexos Multiproteicos/genética , Proteína ORAI1/genética , Ligação Proteica/genética , Conformação ProteicaRESUMO
Mutations in the different domains of A-type lamin proteins cause a diverse plethora of diseases collectively termed as laminopathies which can affect multiple organs. Ig fold is one such domain of lamin A which is implicated in numerous nuclear interactions wherein the mutations lead to different laminopathies. W514R is one such mutation in the Ig fold which leads to severe phenotypes in Skeletal Muscle Dystrophy (SMD) which is a class of laminopathies. In this report, we elucidated gross alterations in structure and dynamics at the level of individual amino acids. These studies indicate altered conformational features of residues in the close vicinity of W514. Imaging of mammalian cells transfected with the mutant have shown distinct perturbation of the nuclear meshwork with concomitant alteration in nuclear interactions as a result of increased oligomerization of Ig W514R. Hence, this novel approach of amalgamating theoretical and experimental procedures to predict the severity of a mutant in the context of laminopathies could be extended for numerous lamin A mutants.
Assuntos
Lamina Tipo A/química , Músculo Esquelético , Distrofias Musculares , Mutação de Sentido Incorreto , Dobramento de Proteína , Substituição de Aminoácidos , Humanos , Lamina Tipo A/genética , Domínios ProteicosRESUMO
Mycobacterium tuberculosis is an obligate pathogen of mammals and is responsible for more than two million deaths annually. The ability to acquire iron from the extracellular environment is a key determinant of pathogenicity in mycobacteria. M. tuberculosis acquires iron exclusively through the siderophores. Several lines of evidence suggest that siderophores have a critical role in bacterial growth and virulence. Hence, in the present study, we have used a combined ligand and structure-based drug design approach for identification of novel inhibitors against salicylate synthase MbtI, a unique and essential enzyme for the biosynthesis of siderophores in M. tuberculosis. We have generated the ligand based and structure based pharmacophores and validated exhaustively. From the validation results it was found that GH (Goodness of Hit) scores for the selected ligand based and structure based pharmacophore models were 0.89 and 0.97, respectively, which indicate that the quality of the pharmacophore models are acceptable as GH value is >0.7. The validated pharmacophores were used for screening the ZINC database. A total of 73 hits, obtained through various insilico screening techniques, were further enriched to 17 hits using docking studies. Molecular dynamics simulations were carried out to compare the binding mode and stability of complexes of MbtI bound with substrate, known inhibitors, and three top ranked hits. The results obtained in this study gave assurance about the identified hits as prospective inhibitors of MbtI.
Assuntos
Proteínas de Bactérias/antagonistas & inibidores , Inibidores Enzimáticos/análise , Inibidores Enzimáticos/química , Simulação de Dinâmica Molecular , Mycobacterium tuberculosis/enzimologia , Proteínas de Bactérias/química , Domínio Catalítico , Ácido Corísmico/química , Ácido Corísmico/metabolismo , Bases de Dados de Proteínas , Avaliação Pré-Clínica de Medicamentos , Inibidores Enzimáticos/farmacologia , Ligação de Hidrogênio , Ligantes , Simulação de Acoplamento Molecular , Mycobacterium tuberculosis/efeitos dos fármacos , Análise de Componente Principal , Reprodutibilidade dos Testes , Ácido Salicílico/química , Ácido Salicílico/metabolismo , Relação Estrutura-Atividade , Especificidade por Substrato/efeitos dos fármacosRESUMO
Mycobacterium tuberculosis, the pathogen responsible for tuberculosis, uses various strategies to survive in a variety of host lesions. The re-emergence of multi-drug-resistant strains of M. tuberculosis underlines the necessity to discover new molecules. Inhibitors of aryl acid adenylating enzyme, MbtA, involved in siderophore biosynthesis in M. tuberculosis, are being explored as potential anti tubercular agents. In this study, we have used 3D-QSAR models and shape based virtual screening to identify novel MbtA inhibitors. 3D-QSAR studies were carried out on nucleoside bisubstrate derivatives. Both Comparative Molecular Field Analysis (r(2) = .944 and r(2)(pred) = .938) and Comparative Molecular Similarity Indices Analysis (r(2) = .892 and r(2)(pred) = .842) models, developed using Gasteiger charges with all fields, predicted efficiently. A total of 13 hits were identified as novel prospective inhibitors for MbtA by utilizing an insilico workflow. Out of 13 hits, five top ranked hits were used for further molecular dynamics studies to gain more insights about the stability of the complexes.
Assuntos
Antituberculosos/química , Ligases/antagonistas & inibidores , Mycobacterium tuberculosis/enzimologia , Domínio Catalítico , Desenho de Fármacos , Ligação de Hidrogênio , Interações Hidrofóbicas e Hidrofílicas , Ligases/química , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Ligação Proteica , Relação Quantitativa Estrutura-AtividadeRESUMO
Multidrug resistance capacity of Mycobacterium tuberculosis demands urgent need for developing new antitubercular drugs. The present work is on M. tuberculosis-MbtA, an enzyme involved in the biosynthesis of siderophores, having a critical role in bacterial growth and virulence. The molecular models of both holo and apo forms of M. tuberculosis-MbtA have been constructed and validated. A docking study with a series of 42 5'-O-[N-(salicyl) sulfamoyl] adenosine derivatives, using GOLD software, revealed significant correlation (R(2) = 0.8611) between Goldscore and the reported binding affinity data. Further, binding energies of the docked poses were calculated and compared with the observed binding affinities (R(2) = 0.901). All-atom molecular dynamics simulation was performed for apo form, holo form without ligand and holo form with ligands. The holo form without ligand on molecular dynamics simulation for 20 ns converged to the apo form and the apo form upon induced fit docking of the natural substrate, 2,3-dihydroxybenzoic acid-adenylate, yielded the holo structure. The molecular dynamics simulation of the holo form with ligands across the time period of 20 ns provided with the insights into ligand-receptor interactions for inhibition of the enzyme. A thorough study involving interaction energy calculation between the ligands and the active site residues of MbtA model identified the key residues implicated in ligand binding. The holo model was capable to differentiate active compounds from decoys. In the absence of experimental structure of MbtA, the homology models together with the insights gained from this study will promote the rational design of potent and selective MbtA inhibitors as antitubercular therapeutics. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:33.
Assuntos
Antituberculosos/farmacologia , Ligases/química , Mycobacterium tuberculosis/enzimologia , Sítios de Ligação , Farmacorresistência Bacteriana Múltipla , Hidroxibenzoatos/química , Ligases/ultraestrutura , Testes de Sensibilidade Microbiana , Modelos Moleculares , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Mycobacterium tuberculosis/efeitos dos fármacos , Análise de Componente Principal , Ligação Proteica , Tuberculose/tratamento farmacológicoRESUMO
Chronic alcoholism is one of the most common causes of liver diseases worldwide. Nitric oxide (NO) has been proposed to have potential for clinical application against chronic hepatocellular injuries. However, mechanisms underlying hepatoprotective functions of NO in ethanol-induced apoptosis are largely unknown. Sprauge-Dawley rats were exposed to ethanol for 8 weeks. Half of the ethanol-fed animals received 14-deoxyandrographolide (14-DAG) treatment for the last 4 weeks of study. Preventive effect of 14-DAG against ethanol-induced hepatotoxicity involved constitutive nitric oxide synthase (cNOS) activation followed by up-regulation of γ-glutamylcysteine synthetase activity and reduced oxidative stress. Enhanced interaction of cNOS with caveolin-1 caused down-regulation of enzyme activity and led to depletion of NO in the hepatocytes of ethanol-fed animals. 14-DAG acted as activator of adenylate cyclase and modulated cyclic AMP (cAMP) mediated expression of caveolin-1 and calmodulin. This eventually favored activation of cNOS through inhibition of cNOS-caveolin-1 interaction. Our results suggest that, protective effect of 14-DAG against ethanol-induced hepatic injury is based on its ability to reduce oxidative stress through cNOS dependent improvement of redox status. 14-DAG mediated activation of adenylate cyclase-cAMP signaling leading to up-regulation of cNOS may provide a promising approach in the prevention of liver diseases during chronic alcoholism.
Assuntos
Adenilil Ciclases/metabolismo , Antioxidantes/uso terapêutico , Diterpenos/uso terapêutico , Hepatopatias Alcoólicas/prevenção & controle , Fígado/efeitos dos fármacos , Óxido Nítrico Sintase Tipo III/metabolismo , Sistemas do Segundo Mensageiro/efeitos dos fármacos , Adenilil Ciclases/química , Andrographis/química , Animais , Antioxidantes/isolamento & purificação , Antioxidantes/farmacologia , Apoptose/efeitos dos fármacos , Células Cultivadas , AMP Cíclico/agonistas , AMP Cíclico/metabolismo , Diterpenos/isolamento & purificação , Diterpenos/farmacologia , Feminino , Glutamato-Cisteína Ligase/química , Glutamato-Cisteína Ligase/metabolismo , Células Hep G2 , Humanos , Fígado/citologia , Fígado/metabolismo , Fígado/patologia , Hepatopatias Alcoólicas/metabolismo , Hepatopatias Alcoólicas/patologia , Óxido Nítrico/agonistas , Óxido Nítrico/metabolismo , Óxido Nítrico Sintase Tipo III/química , Oxirredução , Estresse Oxidativo/efeitos dos fármacos , Folhas de Planta/química , Ratos , Ratos Sprague-DawleyRESUMO
Aimed at identification and structural characterization of novel putative therapeutic targets in H. pylori, the etiological agent of numerous gastrointestinal diseases including peptic ulcer and gastric cancer, the present study comprised of three phases. First, through subtractive analysis of metabolic pathways of Helicobacter pylori HPAG1 and human, as documented in the KEGG database, 11 pathogen-specific pathways were identified. Next, all proteins involved in these pathogen-specific pathways were scrutinized in search of promising targets and the study yielded 25 candidate target proteins that are likely to be essential for the pathogen viability, but have no homolog in human. The lipopolysaccharide (LPS) biosynthesis pathway was found to be the largest contributor (nine proteins) to this list of candidate proteins. Considering the importance of LPS in H. pylori virulence, 3D structural models of three predicted target enzymes of this pathway, namely 2-dehydro-3-deoxy-phosphooctonate aldolase, UDP-3-O-[3-hydroxymyristoyl] N-acetylglucosamine deacetylase and Phosphoheptose isomerase, were then built up using the homology modeling approaches. Binding site analysis and docking of the known biological substrate PEP to 2-dehydro-3-deoxyphosphooctonate aldolase revealed the potential binding pocket present in the single monomeric form of the enzyme and identified 11 amino acid residues that might play the key roles in this protein-ligand interaction.
Assuntos
Aldeído Liases/química , Aldose-Cetose Isomerases/química , Amidoidrolases/química , Proteínas de Bactérias/química , Helicobacter pylori/química , Lipopolissacarídeos/biossíntese , Sequência de Aminoácidos , Aminoácidos/química , Sítios de Ligação , Simulação por Computador , Bases de Dados de Proteínas , Desenho de Fármacos , Helicobacter pylori/enzimologia , Humanos , Ligantes , Metaboloma , Modelos Moleculares , Dados de Sequência Molecular , Fosfoenolpiruvato/química , Ligação Proteica , Homologia Estrutural de Proteína , TermodinâmicaRESUMO
The re-emergence of tuberculosis infections, which are resistant to conventional drug therapy, has steadily risen in the last decade. Inhibitors of aryl acid adenylating enzyme known as MbtA, involved in siderophore biosynthesis in Mycobacterium tuberculosis, are being explored as potential antitubercular agents. The ability to identify fragments that interact with a biological target is a key step in fragment based drug design (FBDD). To expand the boundaries of quantitative structure activity relationship (QSAR) paradigm, we have proposed a Fragment Based QSAR methodology, referred here in as FB-QSAR, for deciphering the structural requirements of a series of nucleoside bisubstrate analogs for inhibition of MbtA, a key enzyme involved in siderophore biosynthetic pathway. For the development of FB-QSAR models, statistical techniques such as stepwise multiple linear regression (SMLR), genetic function approximation (GFA) and GFAspline were used. The predictive ability of the generated models was validated using different statistical metrics, and similarity-based coverage estimation was carried out to define applicability boundaries. To aid the creation of novel antituberculosis compounds, a bioisosteric database was enumerated using the combichem approach endorsed mining in a lead-like chemical space. The generated library was screened using an integrated in-silico approach and potential hits identified.
RESUMO
Dihydrofolate reductase (DHFR) has been used successfully as a drug target in the area of anti-bacterial, anti-cancer and anti-malarial therapy. It also acts as a drug target for Leishmaniasis. Inhibition of DHFR leads to cell death through lack of thymine (nucleotide metabolism). Although the crystal structures of Leishmania major and Trypanosoma cruzi DHFR-thymidylate synthase (TS) have been resolved, to date there is no three-dimensional (3D)-structural information on DHFR-TS of Leishmania donovani chagasi, which causes visceral leishmaniasis. Our aim in this study was to model the 3D structure of L. donovani chagasi DHFR-TS, and to investigate the structural requirements for its inhibition. In this paper we describe a highly refined homology model of L. donovani chagasi DHFR-TS based on available crystallographic structures by using the Homology module of Insight II. Structural refinement and minimization of the generated L. donovani chagasi DHFR-TS model employed the Discover 3 module of Insight II and molecular dynamic simulations. The model was further validated through use of the PROCHECK, Verify_3D, PROSA, PSQS and ERRAT programs, which confirm that the model is reliable. Superimposition of the model structure with the templates L. major A chain, L. major B chain And T. cruzi A chain showed root mean square deviations of 0.69 A, 0.71 A and 1.11 A, respectively. Docking analysis of the L. donovani chagasi DHFR-TS model with methotrexate enabled us to identify specific residues, viz. Val156, Val30, Lys95, Lys75 and Arg97, within the L. donovani chagasi DHFR-TS binding pocket, that play an important role in ligand or substrate binding. Docking studies clearly indicated that these five residues are important determinants for binding as they have strong hydrogen bonding interactions with the ligand.