ISPUB.com / IJGP/4/2/7920
  • Author/Editor Login
  • Registration
  • Facebook
  • Google Plus

ISPUB.com

Internet
Scientific
Publications

  • Home
  • Journals
  • Latest Articles
  • Disclaimers
  • Article Submissions
  • Contact
  • Help
  • The Internet Journal of Genomics and Proteomics
  • Volume 4
  • Number 2

Original Article

Prediction of three-dimensional structure of Chlamydia trachomatis heat shock protein 60: An immunodominant antigen.

R Jha, H Jha, S Biswas, A Mittal

Keywords

3-d model, heat shock protein 60, homology modeling, mgenthreader, ramachandran plot

Citation

R Jha, H Jha, S Biswas, A Mittal. Prediction of three-dimensional structure of Chlamydia trachomatis heat shock protein 60: An immunodominant antigen.. The Internet Journal of Genomics and Proteomics. 2008 Volume 4 Number 2.

Abstract


Chlamydia trachomatis infections are the most prevalent sexually transmitted infection in the world and are often associated with various fertility disorders. C. trachomatis heat shock protein 60 (cHSP60) or GroEL has been identified as an immunodominant antigen, being strong target of both humoral and cell mediated immune responses. A homology modeling was used for the prediction of three-dimensional (3-D) structure of cHSP60 by using the comparative modeling module of MODELLER. Since the MODELLER-aligned alignments led to asymmetric matching, the necessary sequence alignments were done using ClustalW and mGenTHREADER. The template used was the file with code 1kp8 from the Protein Data Bank (PDB), which gives the crystal structure of GroEL of Escherichia coli. The structure was energy minimized by MODELLER’s inherent minimization module. From the Ramachandran plot analysis, it was found that 91.8% of residues falls in the most favoured region. Further, this study could shed light on characterizing cHSP60 in wet laboratory.

 

Introduction

Chlamydia trachomatis infections are the most prevalent sexually transmitted disease recognized throughout the world. According to World Health Organization (WHO), annually 92 million C. trachomatis infections are detected worldwide, of which 43 million cases are reported from South East Asia (WHO, 2001). Up to 80% of women with genital chlamydial infection are asymptomatic and even upper genital tract infection is often clinically silent (Gaydos et al., 1998). Untreated chlamydial infection of lower genital tract can lead to pelvic inflammatory disease (PID) in 10-40% of affected women, which can result in infertility, ectopic pregnancy and chronic pelvic pain (Patric et al., 1997). It is also reported that urogenital infection with C. trachomatis may lead to development of an acute inflammatory arthritis (Inman et al., 2000).

Heat shock proteins (HSPs) are highly conserved proteins present in almost all prokaryotic and eukaryotic organisms. The 60-KDa HSPs have been extensively studied for their chaperone function in protein folding and their cooperation with other chaperones in cellular trafficking. More than 150 homologues of HSP60 sequence are currently available with pair-wise similarity at the amino-acid level extending from 40 to 100% (Brocchieri et al., 2000). In Chlamydiae, the 60-KDa (cHSP60) or GroEL is encoded by genes arranged on the bicistronic groESL operon (Morrison et al., 1990). In a panel of 116 recombinant C. trachomatis proteins, cHSP60 (Ct110) has recently been identified as an immunodominant antigen, being a strong target of both humoral and cell-mediated immune responses (Follman et al., 2008). In addition, cHSP60 is highly homologous and conserved throughout the evolution. It belongs to a group of chaperone proteins ranging in molecular weight from 15 to 110 KDa which are produced in response to infection, stress, inflammation and other cellular insults (Engel et al., 1990). cHSP60 exhibit greater than 80% homology between chlamydia spp., 60% identity with other bacteria and 50% homology with human HSP60 (hHSP60) and exhibits immuno-regulatory properties primarily by inducing proinflammatory response (Zugel et al., 1999). Recently, Habich et. al. reported distinct epitopes responsible for binding to innate immune cells mediating the pro-inflammatory effects of hHSP60 (Habich et al., 2007).

Persistent chlamydial infection has been reported to cause chronic stimulation of host immune system against immunogenic antigen such as cHSP60 (Hessel et al., 2001). Numerous studies have shown a consistent immunopathogenic association between antibody response to the cHSP60 and the development of pelvic inflammatory disease (PID), ectopic pregnancy and tubal infertility (Tiitinen et al., 2006; Hartog et al., 2005; Clad et al., 2003). Also, presence of circulating antibodies to a conserved epitope of the cHSP60 is associated with a lower spontaneous conception rate, and increased likelihood of adverse pregnancy outcome in women treated by salpingectomy for the first episode of ectopic pregnancy (Sziller et al., 2008). We previously reported high seroprevalence of anti-cHSP60 in C. trachomatis infected women (Dutta et al., 2007) and that detection of anti-cHSP60 antibodies would help in the early prognosis of immunopathological sequelae in C. trachomatis infected women (Dutta et al., 2008). In addition, we evaluated cHSP60-specific proliferative responses and our study suggested the probable role of cHSPs in modulation of mucosal immune responses (Agrawal et al., 2007; Srivastava et al., 2008). Recently, we showed higher level of cHSP60 in infertile women than in fertile women in cervical epithelial cells from C. trachomatis infected site (Jha et al., 2009).

Till now 3-D structure is not available for cHSP60 which is important since homology modeling has become a useful tool to understand the particular protein. Although several crystal structures of GroEL from other organisms are available and two of them show GroEL in its 14 unit complex, one associated with ATPγS (Boisvert et al., 1996) and a second in an asymmetrical conformation bound to GroES and ADP (Xu et al., 1997). These authors have further reported seven identical GroEL monomers assembled in rings that dimerize in a 14-meric structure with an equatorial interface formed by the E Domains and two non-communicating central cavities. Each Domain E also provides a binding pocket for an ATP/ADP molecule and a Mg+ ion. Therefore, the present study was undertaken for predicting 3-D structure model of cHSP60, which is required for identifying major target epitopes for wet lab and also for designing new drugs.

Material and Methods

Retrieval of target sequence: The amino acid sequence of cHSP60 (AAC67701.1), hHSP60 (ABB01006.1) and E. coli HSP60 (ABF67773.1) was obtained from sequence database of NCBI (http://www.ncbi.nlm.nih.gov). It was ascertained that 3-D structure of cHSP60 is not available in Protein Data Bank (http://www.rcsb.org/pdb/home.do), hence the present exercise of developing the 3-D model of cHSP60 was ascertained. The cHSP60 is 543 amino acid in length and having molecular weight of 58 KDa.

Template searching: A suitable template for modeling of target protein was searched through mGenTHREADER (McGuffin et al., 2003), which is an online tool for searching similar sequences. From the homology searching high resolution of X-ray crystallography structure of the Escherichia coli [1kp8] was selected as template protein from protein data bank (Wang et al., 2003). Of the 14 chains, only chain A was used for modeling.

Alignment of Sequence: Amino acid sequence of target and template protein was done using MODELLER. Since the MODELLER-aligned alignments led to asymmetric matching, the necessary sequence alignments were done using ClustalW (http://www.ebi.ac.uk/Tools/clustalw2/) (Thompson et al., 1994).

3-D structure modeling: A rough 3-D model of cHSP60 was constructed from the sequence alignment between cHSP60 and template proteins using Modeller 8v0 (Marti-Renom et al., 2000). The structure was energy minimized by MODELLER’s inherent minimization module. In the last stage of modeling the structure was evaluated by the inspection of the Psi/Phi Ramachandran plot obtained from PROCHECK (http://biotech.ebi.ac.uk:8400/cgi-bin/sendquery) (Laskowsky et al., 1993).

Results and Discussion

The Chlamydia genome database (http://chlamydia-www.berkeley.edu) containing both nucleotide and amino acid sequences are available. The retrieved amino acid sequence of cHSP60 from Chlamydia genome database and NCBI showed high homology with sequence of E. coli GroEL. There are two other ortholog of cHSP60 protein copies in Chlamydiae species have diverged functionally after the gene duplication events (McNally et al., 2007) but cHSP60 was found more abundant and immunodominant than the other two. The transcription of cHSP60 gene in C. trachomatis has been shown to be immediately induced in infected cells exposed to higher temperature at 45ºC (Karunakaran et al., 2003), however, its specific function both under stress and normal conditions remains unclear. We therefore undertook bioinformatics analysis of 3-D structure that will help in understanding chlamydial disease pathogenesis.

Homology modeling of proteins has been successful in gaining better understanding of protein function and intermolecular interactions. Manual methods have been successful but, with the advent of genomic sequences database and the significant increase in data, more efficient methods are needed. Therefore, homology modeling was used to produce best model of unknown structure by comparing it with a known template from structurally highly similar molecule. In the multiple alignments of HSP60 proteins, the most conserved residue is glycine (G) (Fig. 3). Brocchieri et. al. earlier reported that sequence positions essential to protein function and structure are conserved over a broad evolutionary range. Freely variable regions of a protein family are generally non functional. It is well recognized that G residues can be a “hinge” for backbone conformations that are generally excluded for other side chains. Conserved G multiplets are also observed in other HSP families, like HSP70 (Karlin et al., 1998) or RecA (Brendel et al., 1997). The Ramachandran plot of all perfectly conserved G residues reveals that many of them belong to backbone conformations or to hinge connections relevant to structural reconfiguration (Fig. 1).

In our study based on the results obtained from mGenTHREADER and ClustalW, the chain-A of X-ray structure of the E. coli GroEL were selected as template which was 61% identical with cHSP60. The MODELLER was used for building the model and was further examined and validated using PROCHECK. Further, the total energy values of the predicted 3-D model were calculated as 91.8% of Ramachandran plot (Fig. 1) and the refined model was analyzed by different protein analysis programs including PROCHECK for the evaluation of the Ramachandran plot quality, and WHATIF (Vriend et al., 1990) for the calculation of packing quality. This structure for the corresponding coordinates was found to be satisfactory. The residue properties of cHSP60 indicates that 3-D structure possess 22 α-helix and 18 β-strand (Fig. 2). One of the epitopes (151SANNDAEIGNLI162) in the structure falls in the accessible region. Our earlier report suggests that detection of antibodies to this epitope would help in early prognosis of immunopathological sequelae in C. trachomatis infected women and thereby instituting appropriate therapy for controlling infection at any stage (Dutta et al., 2008).

Figure 1
Figure 1: (A) Predicted 3-D structure of cHSP60. (B) Ramachandran plot analysis of predicted 3-D structure of HSP60. Based on analysis of 118 structures of resolution with at least 2 Å and R factor no greater than 20%. A good quality model would be expected to have over 90% in the most favored regions. The Plot statistics are: Residues in most favored regions [A,B,L]- 438 (91.8%); residues in additional allowed regions [a,b,l.p] - 29 (6.1%); residues in generously allowed regions [-a,-b,-l,-p] - 4 (0.8%); residues in disallowed regions – 6 (1.3 %); number of non-glycine and non-proline residues- 477 (100.0%); number of end residues (excl. Gly and Pro) - 2; number of glycine residues (shown as triangles) - 45; number of proline residues - 17; total number of residues- 541.

Since, cHSPs are phylogenetically conserved and have high sequence homology, the immune response against epitopes within the region of cHSP60 may elicit an autoimmune response due to cross-reactivity to epitopes of the hHSP60. Also, the homology between bacterial and hHSP60 and the fact that hHSP60 may be increased locally at the site of infection may indicate a role for HSP60 autoimmunity in the immunopathogenicity of chlamydial infections (Jones et al., 1993). Further, Yi et. al. reported that seven of 13 B-cell peptide epitopes recognized by human anti-sera exhibited cross-reactive antibody binding to homologous peptide sequences in hHSP60. Such self-reactive B-cell immunity to HSP60 may contribute to chlamydial disease pathogenesis (Yi et al., 1993). Similarly, Domeika et al. examined the association between humoral immunity to unique and conserved epitopes of the cHSP60 and immunity to hHSP60. They have demonstrated that sensitization to hHSP60 is associated with humoral immune response to a conserved epitope of the cHSP60 in patients with PID and also there was serologic evidence of exposure to chlamydiae. This suggests that an autoimmune response to hHSP60 may arise as a consequence of a C.trachomatis upper genital tract infection in those women who develop sensitivity to cHSP60 epitopes that cross-react with epitopes of the hHSP60 (Domeika et al., 1998).

Figure 2
Figure 2: Amino-acid sequence of cHSP60: Secondary structural elements are indicated as α-Helix, β-strand and random coil. Blue (Buried) and white (accessible) regions in the sequence are shown as accessibility shading.

These observations were also supported by co-immunization with mouse HSP60 and cHSP60 inducing strong T and B-cell responses to self HSP60. In particular, T-cell proliferation in response to mouse HSP60 was dramatically increased upon co-immunization. The induction of strong autoimmune responses depended on the homology of amino acid sequences between mouse and cHSP60, since a strong response was not observed when a non-homologous but highly immunogenic protein such as ovalbumin was co-administered with mouse HSP60. The importance of sequence homology implicates cross reaction between the antigen recognition structures on T and/or B cells in the autoimmune response suggesting that molecular mimicry is the basis for autoimmunity induced by HSP60 (Lin et al., 1991). On the other hand, production of cHSP60 antibodies may also dominate over antibodies to another antigenic component of chlamydiae or high degree of amino acid sequence homology between cHSP60 and other microbial HSP60 that may incite an autoimmune reaction through molecular mimicry (Bavoil et al., 1990). Further, immune responses generated after cross reaction may lead to mediate inflammation in lower and upper genital tract, a key event in the cervicitis, PID, ectopic pregnancy and infertility.

Recently, Kumar et al. implicated HSP90 from Plasmodium falciparum as potential drug target against malaria. Using inhibitor specific to nucleotide binding domain of HSP90, they have shown potent growth inhibitory effect on development of malarial parasite in human erythrocyte (Kumar et al., 2007). Thus, knowledge of 3-D structural details of cHSP60 would also greatly advance the development of drugs targeting this molecule. Further, analysis is required to found the differences between hHSP60 and cHSP60 may make it possible to design inhibitors specific against Chlamydia.

Figure 3
Figure 3: Multiple sequence alignment of HSP60 amino-acid residues of and . The *(star) in the sequence represents identical residues. (C.tra- , E.coli- , H.sap- )

Acknowledgements

We acknowledge the Indian Council of Medical Research (ICMR) and University Grants Commission (UGC) for providing financial assistance to Rajneesh Jha, Hem Chandra Jha and Sumit Biswas in the form of a research fellowship. We also thank Mr. Harsh Vardhan, Mr. Rishein Gupta, Ms. Pragya Srivastava and Ms. Apurb Rashmi Bhengraj for his suggestions and critical reading of the manuscript.

Correspondence to

Aruna Mittal Deputy Director (Sr. Grade) Institute of Pathology- ICMR Safdarjung Hospital Campus Post Box no. 4909 New Delhi – 110 029, India Tel. no. - 091-011-26198 402-05 Fax no. - 091-011-26198 401 E-mail - amittal_cp@rediffmail.com

References

1. Global Prevalence and Incidence of selected curable Sexually Transmitted Infections-Overviews and estimates. (1st February, 2001) World Health Organization, Geneva. Available at: http://www.who.int/hiv/pub/sti/who_hiv_aids_2001.02.pdf.
2. Gaydos CA, Howell MR, Pare B, Clark KL, Ellis DA, Hendrix RM, Gaydos JC, McKee KT Jr, Quinn TC. Chlamydia trachomatis infections in female military recruits. N Engl J Med 1998;339(11):739–44.
3. Patrick DM. Chlamydia control: Components of an effective control strategy to reduce the incidence of Chlamydia trachomatis. Can J Human Sex 1997;6:143-9.
4. Inman RD, Whittum-Hudson JA, Schumacher HR, Hudson AP. Chlamydia-associated arthritis. Curr Opin Rheumatol 2000;12(4):254-62.
5. Brocchieri L, Karlin S. Conservation among HSP60 sequences in relation to structure, function, and evolution. Protein Science 2000;9(3):476–86.
6. Morrison RP, Su H, Lyng K, Yuan Y. The Chlamydia trachomatis hyp operon is homologous to the groE stress response operon of Escherichia coli. Infect Immun 1990;58(8):2701-5.
7. Follmann F, Olsen AW, Jensen KT, Hansen PR, Andersen P, Theisen M. Antigenic profiling of a Chlamydia trachomatis gene-expression library. J Infect Dis 2008;197(6):897-905.
8. Engel JN, Pollack J, Perara E and Gamen D. Heat shock response of murine Chlamydia trachomatis. J Bacteriol 1990;172(12):6959-72.
9. Zugel U, Kaufmann SHE. Role of Heat shock proteins in protection from and pathogenesis of infectious disease. Clin microbiol rev 1999;12(1):19-39.
10. Habich C, Burkart V. Heat shock protein 60: regulatory role on innate immune cells. Cell Mol Life Sci 2007;64(6):742-51.
11. Hessel T, Dhital, SP, Plank R, Dean D. Immune response to 60-kilodalton heat shock protein in tears from Nepali trachoma patients. Infect Immun 2001;69(8):4996-5000.
12. Tiitinen A, Surcel HM, Halttunen M, Birkelund S, Bloigu A, Christiansen G, Koskela P, Morrison SG, Morrison RP, Paavonen J. Chlamydia trachomatis and chlamydial heat shock protein 60-specific antibody and cell-mediated responses predict tubal factor infertility. Hum Reprod 2006;21(6):1533-8.
13. Hartog den J E, Land JA, Stassen FR, Kessels AG, and Bruggeman CA. Serological markers of persistent C. trachomatis infections in women with tubal factor subfertility. Hum Reprod 2005;20(4):986-90.
14. Clad A, Petersen EE, Dettlaff S. Antibodies to Chlamydia trachomatis heat shock protein 60 (cHSP60) and Chlamydia trachomatis major outer membrane protein (MOMP) in women with different tubal status. Clin Lab 2003;49(5-6):269-71.
15. Sziller I, Fedorcsa´k P, Csapo´ Zs, Szirmai K, Linhares IM, Papp Z, Witkin SS. Circulating Antibodies to a Conserved Epitope of the Chlamydia trachomatis 60-kDa Heat Shock Protein is Associated with Decreased Spontaneous Fertility Rate in Ectopic Pregnant Women Treated by Salpingectomy. Am J Reprod Immunol 2008;59(2):99–104.
16. Dutta R, Jha R, Gupta S, Gupta R, Salhan S, Mittal A. Seroprevalence of antibodies to conserved region of Chlamydia trachomatis heat shock proteins 60 and 10 in women attending Gynecology out patient department. Br J Biomed Sci 2007;64(2):78-83.
17. Dutta R, Jha R, Salhan S, Mittal A. Chlamydia trachomatis-specific heat shock proteins 60 antibodies can serve as prognostic marker in secondary infertile women. Infection 2008;36(4):374-8.
18. Agrawal T, Vats V, Salhan S, Mittal A. Mucosal and peripheral immune responses to chlamydial heat shock proteins in women infected with Chlamydia trachomatis. Clin Exp Immunol 2007;148(3):461-8.
19. Srivastava P, Jha R, Bas S, Salhan S, Mittal A. In infertile women, cells from Chlamydia trachomatis infected site release higher levels of Interferon-gamma, Interleukin-10 and Tumor Necrosis Factor-alpha upon heat shock protein stimulation than fertile women. Repro Biol Endocrinol 2008;6:20.
20. Jha R, Vardhan H, Bas S, Salhan S, Mittal A. In Infertile women, cervical epithelial cells from Chlamydia trachomatis infected site co-express higher level of chlamydial heat shock protein 60 and 10 than in fertile women. Gynecol Obstet Invest 2009. (In press)
21. Boisvert DC, Wang J, Otwinowski Z, HorwichAL, Sigler PB. The 2.4 Å crystal structure of the bacterial chaperonin GroELcomplexed with ATPγS. Nat Struct Biol 1996;3(2):170–7.
22. Xu Z, Horwich AL, Sigler PB. The crystal structure of the asymmetric GroEL-GroES-(ADP)7 chaperonin complex. Nature 1997;388(6644):741–50.
23. McGuffin LJ, Jones DT. Improvement of the GenTHREADER method for genomic fold recognition. Bioinformatics 2003;19(7):874-81.
24. Wang J, Boisvert DC. Structural Basis for GroEL-assisted Protein Folding from the Crystal Structure of (GroEL-KMgATP)14 at 2.0 A° Resolution. J Mol Biol 2003;327(4):843–55.
25. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994;22(22):4673-80.
26. Marti-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, Sali A. Comparative protein structure modeling of genes and genomics. Annu Rev Biophys Biomol Struct 2000;29:291-325.
27. Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst 1993;26(2):283-91.
28. McNally D, Fares AM. In silico identification of functional divergence between the multiple groEL gene paralogs in Chlamydiae. BMC Evolutionary Biology 2007;7:81.
29. Karunakaran KP, Noguchi Y, Read TD, Cherkasov A, Kwee J, Shen C, Nelson CC, Brunham RC. Molecular analysis of the multiple groEL proteins of Chlamydiae. Journal of Bacteriol 2003;185(6):1958–66.
30. Karlin S, Brocchieri L. Heat shock protein 70 family: Multiple sequence comparisons, function, and evolution. J Mol Evol 1998;47(5):565–77.
31. Brendel V, Brocchieri L, Sandler SJ, Clark AJ, Karlin S. Evolutionary comparisons of RecA-like proteins across all major kingdoms of living organisms. J Mol Evol 1997;44(5):528–41.
32. Vriend G. WHAT IF: A molecular modeling and drug design program. J Mol Graph 1990;8(1):52-6.
33. Jones DB, Coulson AWF, Duff GV. Sequence homologies between hsp60 and autoantigens. Immunol Today 1993;14(3):115–8.
34. Yi Y, Zhomg G, Brunham R. Continuous B-cell epitopes in Chlamydia trachomatis heat shock protein 60. Infect Immun 1993;61(3):1117-20.
35. Domeika M, Domeika K, Paavonen J, Mardh PA, Witkin SS. Humoral Immune Response to Conserved Epitopes of Chlamydia trachomatis and Human 60-kDa Heat-Shock Protein in Women with Pelvic Inflammatory Disease. The J Infect Dis 1998;177(3):714–9.
36. Lin RH, Mamula MJ, Hardin JA, Janeway CA. Induction of autoreactive B cells allows priming of autoreactive T cells. J Exp Med 1991;173(6):1433-9.
37. Bavoil P, Stephens RS, Falkow S. A soluble 60 kilodalton antigen of Chlamydia spp. Is a homologue of Escherichia coli GroEL. Mol Microbiol 1990;4(3):461–9.
38. Kumar R, Pavithra SR, Tatu U. Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria; J Biosci 2007;32(3):531–6.

Author Information

Rajneesh Jha, M.Sc.
Safdarjang Hospital Campus, Institute of Pathology- ICMR

Hem Chandra Jha, M.Sc.
Safdarjang Hospital Campus, Institute of Pathology- ICMR

Sumit Biswas, M.Sc.
Department of Biochemistry, Bose Institute

Aruna Mittal, Ph.D.
Safdarjang Hospital Campus, Institute of Pathology- ICMR

Download PDF

Your free access to ISPUB is funded by the following advertisements:

 

BACK TO TOP
  • Facebook
  • Google Plus

© 2013 Internet Scientific Publications, LLC. All rights reserved.    UBM Medica Network Privacy Policy