The history of leprosy and its interaction with man is one of suffering and misunderstanding. The newest research suggests that at least as early as 4000 B.C. individuals had been infected with M. leprae, while the first known written reference to the disease was found on Egyptian papyrus in about 1550 B.C. It is a chronic infectious disease afflicting between 10 and 15 million people, is caused by the obligate intracellular parasite Mycobacterium leprae. It causes damage to the skin and the peripheral nervous system. The disease develops slowly and results in skin lesions and deformities, most often affecting the cooler places on the body such as eyes, nose, earlobes, hands, feet, and testicles. The skin lesions and deformities can be very disfiguring and are the reason that infected individuals historically were considered outcasts in many cultures. Although human-to-human transmission is the primary source of infection, three other species can carry and (rarely) transfer M. leprae to humans: chimpanzees, mangabey monkeys, and nine-banded armadillos.
Leprosy is caused mainly by Mycobacterium leprae, a rod-shaped bacillus that is an obligate intracellular (only grows inside of certain human and animal cells) bacterium. The early signs and symptoms of leprosy are very subtle and occur slowly. Numbness and loss of temperature sensation are some of the first symptoms that patients experience. As the disease progresses, the sensations of touch, then pain, and eventually deep pressure are decreased or lost.

The protein selected for pharmacophore studies and structure based drug designing is serine hydroxymethyltransferase. Enzyme commission (E.C) number of the and the sequence length is 426aa. Gene name for the protein is glyA. The main function of serine hydroxymethyltransferase is the interconversion of serine and glycine.

Materials required for the further proceeding of the project work includes information collection regarding the protein sequence (SHMT) in FASTA format from the SwissProt database and the Protein Information Resource (PIR) database, relevant literature support from the PubMed database, metabolic pathway information of the protein from the KEGG pathway database and selection of ligand molecules from the DrugBank database in SMILES notation.

The entire project work was performed with the help of Accelrys Discovery Studio 2.5 and using various tools present within the software. The list of tools being used during the project work can be mentioned as such:

1. For sequence alignment the tool "BLAST search (NCBI server)" was used.
2. The tool "Align Sequence to Template" was being used for the alignment of the protein target sequence to the template sequence.
3. For modeling of the protein sequence (SHMT) with reference to the template sequence, the tool "Build Homology Models" was used.
4. The verification of the protein 3D structure was done by using the the tool "Verify Protein (Profiles-3D)".
5. For common pharmacophore studies of both the ligand and the receptor molecule, the tool "Common Feature Pharmacophore" was being used and for searching of ligands with common pharmacophore, the tool "Search 3D Database" was used.
6. Two tools namely "ADMET Descriptors" and "Toxicity Prediction (TOPKAT)" under the protocol "ADMET" was being used to check whether the ligand molecules satisfies all the required properties for receptor-ligand interaction such as "Aqueous solubility", "Blood brain barrier penetration(BBB)", "Cytochrome P450 (CYP450) 2D6 inhibition", "Hepatotoxicity", "Human intestinal absorption (HIA)", and "Plasma protein binding" and to check whether the ligand molecule is toxic under any condition before binding to the receptor molecule, if so we can easily sort out that specific ligand molecule before further proceeding.
7. The protein receptor molecule has to be optimized before binding to the ligand molecule. For doing so, the receptor molecule has to be stimulated for the generation of force field "CHARMm", and then specifying the active sites for binding of the ligand molecule by using the tool "Receptor-Ligand Interaction".
8. Finally, docking of the receptor to the ligand molecule is done by using the tool "Dock Ligands (LigandFit)" which is based upon the "Ligand Fit" algorithm.

The results of the project work can be summurised in the following way:
Sequence analysis of the target protein through the BLAST server returned with 41 hits, out of which, the sequence 3H7F which has a sequence length of 447aa, with 89% identity and 93% positive was selected as the template sequence for homology modeling of the unstructured protein serine hydroxymethyltransferase (SHMT). Three homology models of the protein SHMT was modeled, out of which, the model with a DOPE score of -46959.18 was selected for further studies. The DOPE score gives the structural stability of the protein structure. A low DOPE score indicates for a better stability of the 3D structure. For binding of the receptor to ligand molecules, the ligand molecules were scanned for pharmacophore regions common to both the receptor and the ligand. The ADMET results shows that only two amino acids falls out of the graph region. These two amino acid comes outside of the domain region of the protein. The domain region of the protein starts from 26 to 419, and that the two amino acids lies in the region of 119 and 430. Since all the amino acids comes under the domain region so no energy minimization was performed. Toxicity prediction of the ligand BTB03887 on male rat, female rat, male mouse, female mouse and ames mutagenicity shows no carcinogenicity. This shows that the ligand molecule can be considered for docking process.

Fig 1. Active site "Site1"
The above given figure shows the active site "Site 1" for binding of the ligand BTB03887 to the receptor molecule.
Finally, docking of the ligand and the receptor molecule was done by using the Dock Ligands (LigandFit) tool. A snapshot of the docked receptor and ligand molecule is shown.

Fig 2. Result for the docked molecule
The above figure shows the distance between between that are bonded. Distance between the atoms H34 and O is 1.760. Distance between the atoms H31 and O is 2.325. The overall dock score gives a value of 51.443.

 B. Myrvang, D. S. Ridley, S. S. Fröland, and Y. K. Song. "Immune responsiveness to Mycobacterium leprae and other mycobacterial antigens throughout the clinical and histopathological spectrum of leprosy". Advanced Bacterial Genetics (Cold Spring Harbor Laboratory, New York, 1980).
 Bullock WE, Jr, Fasal P. Studies of immune mechanisms in leprosy. "The role of cellular and humoral factors in impairment of the in vitro immune response". J Immunol. 106(4):888-899.
 Estrada-Parra S. 1972 Feb. "Immunochemical identification of a defined antigen of Mycobacterium leprae. Infect Immun". 5(2):258-259.
 Gaugas JM. 1968 Dec 21. "Enhancing effect of antilymphocytic globulin on human leprosy infection in thymectomized mice". Nature. 220(5173):1246-1248
 Godal T, Myklestad B, Samuel DR, Myrvang B. 1971 Dec. "Characterization of the cellular immune defect in lepromatous leprosy: a specific lack of circulating Mycobacterium leprae-reactive lymphocytes". Clin Exp Immunol. 9(6):821-831
 Jondal M, Holm G, Wigzell H. 1972 Aug 1. "Surface markers on human T and B lymphocytes. I. A large population of lymphocytes forming nonimmune rosettes with sheep red blood cells". J Exp Med. 136(2):207-215.
 Lane FC, Unanue ER. 1972 May 1. "Requirement of thymus (T) lymphocytes for resistance to listeriosis". J Exp Med. 135(5):1104-1112.
 Lerner RA, Glassock RJ, Dixon FJ. 1967 Dec 1. "The role of anti-glomerular basement membrane antibody in the pathogenesis of human glomerulonephritis". J Exp Med. 126(6):989-1004.
 Mackaness GB. 1969 May 1. "The influence of immunologically committed lymphoid cells on macrophage activity in vivo". J Exp Med. 129(5):973-992.
 Milde EJ, Tönder O. 1968 Aug. "Demonstration of rheumatoid factor in tissue by mixed agglutination with tissue sections". Arthritis Rheum. 11(4):537-545.
 Misch EA, Berrington WR, Vary JC Jr, Hawn TR. "Leprosy and the human genome". Department of Medicine, University of Washington School of Medicine, Seattle, Washington 98195, USA.
 Monot M, Honoré N, Garnier T, Zidane N, Sherafi D, Paniz-Mondolfi A, Matsuoka M, Taylor GM, Donoghue HD, Bouwman A, Mays S, Watson C, Lockwood D, Cole ST. 2010 Apr 42. "Comparative genomic and phylogeographic analysis of Mycobacterium leprae". Institut Pasteur, Paris, France, Nat Genet. (4):361. Khamispour, Ali .
 Rees RJ, Waters MF, Weddell AG, Palmer E. 1967 Aug 5. "Experimental lepromatous leprosy". Nature. 215(5101):599-602.
 Rees RJ, Weddell AG, Palmer E, Pearson JM. 1969 Jul 26. "Human leprosy in normal mice". Br Med J. 3(5664):216-217.
 Richard A. Young, Vijay Mehra, Douglas Sweetser,Thomas Buchanan, Josephine Clark-Curtiss, Ronald W. Davis & Barry R. Bloom. "Genes for the major protein antigens of the leprosy parasite Mycobacterium leprae". Department of Microbiology and Immunology, and Irvington House Institute for Medical Research, Albert Einstein College of Medicine, Bronx, New York 10461, USA.
 Ridley DS, Jopling WH. 1969 Jul. "Classification of leprosy according to immunity". A five-group system. Int J Lepr Other Mycobact Dis. 1966 Jul-Sep;34(3):255-273.
 S N Cho, D L Yanagihara, S W Hunter, R H Gelber and P J Brennan. "Serological specificity of phenolic glycolipid I from Mycobacterium leprae and use in serodiagnosis of leprosy". The Lancet, Volume 314, Issue 8150, Pages 994-996.
 Sandrine Roussel, Bénédicte Rognon, Coralie Barrera, Gabriel Reboux, Karine Salamin, Frédéric Grenouillet, Isabelle Thaon, Jean-Charles Dalphin, Isabelle Tillie-Leblond, Manfredo Quadroni, Michel Monod and Laurence. "Millon Immuno-reactive proteins from Mycobacterium immunogenum useful for serodiagnosis of metalworking fluid hypersensitivity pneumonitis". J. Biol. Chem. 255, 1885-1890.

About Author / Additional Info:
An optimistic author from India with lots of hopes that science will revolutionise the world