Three-dimensional studies of pathogenic peptides from the c-terminal of Trypanosoma cruziribosomal P proteins and their interaction with a monoclonal antibody structural model
© Martín et al 2009
Received: 25 September 2008
Accepted: 27 May 2009
Published: 27 May 2009
The acidic C-terminal peptides from Trypanosoma cruzi ribosomal P proteins are the major target of the antibody response in patients suffering Chagas chronic heart disease. It has been proposed that the disease is triggered by the cross-reaction of these antibodies with the second extra cellular loop of the β1-adrenoreceptor, brought about by the molecular mimicry between the acidic C-terminal peptides and the receptor's loop. To improve the understanding of the structural basis of the autoimmune response against heart receptors, the 3-dimensional structure of the C-terminal peptides of Trypanosoma cruzi ribosomal proteins P0 (EDDDDDFGMGALF) and P2β (EEEDDDMGFGLFD) were solved using the Electrostaticaly Driven MonteCarlo method. Their structures were compared with the second extra-cellular loop of our homology model of human rhodopsin and the existing experimental NMR structures of the C-terminal peptides from human P0 (EESDDDMGFGLFD) and from Leishmania braziliensis P0 (EEADDDMGFGLFD). Docking of Trypanosoma cruzi peptides P0, P2β and human rhodopsin loop into our anti-P2β monoclonal antibody homology model allowed to explore their interactions.
The solution structure of peptides P0 and P2β can be briefly described as a bend. Although the global conformations of the peptides are not identical they shared a common region of four residues (3 to 6) that have a similar structure. The structural alignment of the five peptides also showed a surprising conformational similarity for the same residues. The antibody model and docking studies revealed a most remarkable feature in the active site, a positively charged, narrow and deep cavity where the acidic residues 3 to 6 were accommodated. These results suggest that the most important elements in the molecular peptide recognition by the antibody may be the shape of the loop and the presence of negative charges in positions 3–5 (P0, P2β) or a negative charge in position 4 (rhodopsin loop). This work describes clearly the interactions of the structural elements involved in the autoimmune mechanism of anti-P auto-antibodies cross-reaction and stimulation of the β1-adrenoreceptor and the visual pigment rhodopsin. Results from this study could lead eventually to the development of treatments to abolish receptor mediated symptoms in Chagas. PACS code: 87.15.-v
The Chagas disease or American tripanosomiasis is a parasitosis produced by the flagellated protozoa Trypanosoma cruzi. It is transmitted by the insect vector Triatoma infestans. According to a recent report of the Scientific Working Group on Chagas Disease the overall prevalence of human Trypanosoma cruzi infection is estimated at 16–18 million cases. Approximately 120 million people or 25% of Latin American inhabitants are at risk of contracting the infection.
People infected with Trypanosoma cruzi may suffer cardiac, gastrointestinal or neurological damage. Although disease manifestations vary widely from one area to another, it is estimated that 25–30% of the infected people will progress to irreversible cardiac, oesophageal and colonic pathologies, causing considerable morbidity and mortality .
Chagas's chronic heart disease (cChHD) is a myocarditis characterized by arrhythmia, tachycardia and heart failure. In the chronic phase of the disease, patients present antibodies against the acidic rich C-terminal of Trypanosoma cruzi ribosomal P proteins . It has been proposed that cChHD is an autoimmune disease produced by the existence of molecular mimicry between the C-terminal of ribosomal P proteins and the second extra cellular loop of the β1-adrenergic receptor (H26R) .
Ribosomal P proteins are a protein complex which forms the long and protruding region called the stalk in the large ribosome subunit. In eukaryotes the P family encompasses protein P0 (34 kD), P1 and P2 (~10 kD) . A third protein, P3, has been found in plants . The number of proteins in the P1 and P2 families varies among species. In mammals, insects and fungus, the complex is formed by P0 and two copies of P1 and P2 [6, 7]. Yeast and protozoa like Trypanosoma cruzi have five members: P0, P1α, P1β, P2α and P2β [8–10]. All of them have a conserved acidic motif in the C-terminal end and show high sequence identity with the last one hundred amino acids of P0 C-terminal.
P013 (EDDDDDFGMGALF) and R13 (EEEDDDMGFGLFD) are the peptides corresponding to the last 13 amino acids of Trypanosoma cruzi P0 and P2β respectively and a major target of the humoral immune response against the parasite .
The active immunization with P0 and P2β recombinant proteins from Trypanosoma cruzi has been shown to induce stimulation of the β1-adrenoreceptor in murine models. This activity has been correlated with ventricular arrhythmia, abnormal repolarizations and conduction defects . Moreover, it has been proved that monoclonal antibodies against H26R induce apoptotic effects in cardiomyocites . Another study has associated the dysfunction of the retina, suffered by chagasic patients, with a similar autoimmune mechanism in which anti-P auto-antibodies would cross-react with the visual pigment rhodopsin .
A preliminary study of the cross reaction between the anti-P2β monoclonal antibody (mAb17.2) and the β1-adrenergic receptor was done by Smulski et al. . They prove that an antibody induced by a Trypanosoma cruzi intracellular protein has the ability to bind a peptide which is part of the extracellular domain of a human cardiac receptor. Although such results include a paragraph about molecular modelling studies describing the manual docking of minimized structures of charged hepta peptides from P2β C-terminal (EDDDMGF) and the second extra cellular loop of the β1 adrenergic receptor (ESDEARRCYN) into an homology model of the antibody, they did not report detailed interaction studies.
The aim of this work is to contribute to a deeper understanding of the molecular basis of Chagas's chronic heart disease through the structural study of the C-terminal peptides and their interaction with the homology model of the anti-P2β monoclonal antibody(mAb17). The solution conformation structures of peptides P013 and R13 are the result of the global optimization of the all atom force-field ECEPP/3 plus implicit hydration models using the EDMC method (Electrostatically Driven Monte Carlo) [15, 16].
The EDMC structures were compared with the experimental NMR structures of human P0 H13, PDB code 1S4J (EESDDDMGFGLFD), Leishmania braziliensis P0 A13, PDB code 1S4H (EEADDDMGFGLFD) , and SEL (RYIPEGLQCSCGIDYYTLKPEVNNES), the second extra-cellular loop of our human rhodopsin homology model.
Conformational search to define the three dimensional structure of R13 and P013.
Homology modeling of rhodopsin and monoclonal antibody mAb17.2
Molecular docking. The antigen-antibody interaction was analyzed through the docking of the different peptides onto the target antibody mAb17.2
2.1. Conformational search
Electrostaticaly driven Monte Carlo (EDMC). The first step consisted in generating feasible target peptide conformations using the EDMC method. The method combines stochastic Monte Carlo algorithms and the optimization of the electrostatic interactions, whereby the molecule evolves toward conformations in which the charge distribution becomes energetically more favourable.
2.2. Evaluation of the Conformational Energy
An all-atom chain representation was used with the ECEPP3 force field [18–22]. Two alternative forms of the potential energy function were used to evaluate the total energy during the conformational search. One, in vacuum corresponding to the ECEPP/3 force field and the other adding a simple solvent interaction model of surface accessibility [23–25].
2.3. The Starting Point and conformational search
The amino and carboxyl termini of the target fragment of the P protein molecule were blocked by amino-COCH3 and carboxyl-NH2 groups, respectively. An ensemble of conformations was generated for each peptide starting from (I) canonical α-helix, assigning the values (-60°,-40°) to dihedral angles (Φ,Ψ) of each amino acid residue, while keeping all values of ω to 180°. (II) polyproline conformation, assigning the values (-80°,150°) to dihedral angles (Φ,Ψ) of each amino acid residue, while keeping all values of ω to 180° (III) canonical β-strand, assigning the values (-139°,135°) to dihedral angles (Φ,Ψ) of each amino acid residue, while keeping all values of ω to 180° (IV) Random conformation, the simulations were started from two different initial random conformations, for each of these random conformations, all backbone and side-chain dihedral angles were chosen randomly between 180.0° and -180.0°, with the exception of the dihedral angles ω of the peptide group, which were always chosen in the trans (180.0°) conformation. All backbone and side-chain dihedral angles were allowed to vary freely during the simulations.
2.4. pH conformational coupling
Many research groups [22, 26–31] have devoted a large amount of effort to develop fast and reliable algorithms provide a solution to the problem of ionization equilibrium. However, the incorporation of these algorithms into conformational searches of short peptides is computationally too time consuming. Here, we used an alternative approach in which the energy of the whole ensemble of conformations, generated in a conformational search, is re-evaluated by considering the problem of ionization equilibrium at a given pH explicitly . This methodology permitted us to compute the average degree of charge to each of the ionizable groups in the sequence. Allowing us to find out the importance and influence of the ionizable groups in the peptides, the influence of solvent polarization effects and electrostatic interactions in the discrimination of the lowest-energy conformation.
Where Eint(rp) is the internal conformational energy of the molecule in the absence of solvent, corresponding to the ECEPP/3 energy of the neutral molecule [18–22]; Fvib(rp) is the conformational entropy contribution calculated by the harmonic vibrational contribution with zero charge [33, 34]. Fcav(rp) is the free energy associated with the process of cavity creation when transferring the molecule from the gas phase into the aqueous solution; this term is proportional to the solvent accessible surface of the molecuele. Fsolv(rp) is the free energy associated with the polarization of the aqueous solution calculated by using the Multigrid Boundary Element method developed by Vorobjev and Scheraga . Finz(rp, pH) is the free energy associated with the change in the state of ionization of the ionisable groups due to the transfer of the molecule from the gas phase to the solvent at a fixed pH value. This energy is calculated by using the general multi-site titration formalism [26, 28]. This approximation computes the 2N ionization states, where N is the total number of ionisable groups. The present method provides an accurate calculation of the potential of mean forces between ionized groups of the protein and the pK shifts of the ionisable groups as a function of the peptide environment.
2.5. Homology Modeling
First, the structure of the mAb17.2 Fv domain was modelled with MODELLER , using as template the PDB code 1NBV (82% sequence identity). PDB code 1S5I was used to model the heavy chain CDR3 (TYWGQGTLVTVS, 100% sequence identity). BLAST  and T-Coffee  were used to search for the templates and perform the alignment, respectively.
Second, the model of human rhodopsin was generated in a similar way using the bovine rhodopsin structure as template (PDB 1U19, 93% sequence identity). The quality of the models was assessed by the WHAT_CHECK standalone program .
The docking of peptides R13, P013 and SEL into mAb17.2 was done using BIGGER (Biomolecular complex Generation with Global Evaluation and Ranking) implemented on Chemera 2.0 . BIGGER has a soft docking algorithm, an intermediate approach between flexible and rigid docking. The flexibility of the side chains was treated in a crude implicit way by allowing partial superposition of the side chains. Default options were used except for the angular step that was set at 12°. For R13, 5000 putative complexes were generated and then clustered using a 2.5 Å RMSD. P013 was treated in the same way. About 8000 complexes were generated with full length rhodopsin and then clustered using a 2.5 Å RMSD. An additional score based on the distance from the amino acids of the probe to the amino acids of the CDRs was generated. In order to discard unreal solutions, the best results were analysed visually. Three results were chosen for rhodopsin. For the R13 and P013 peptides the final sets were around 15 possible solutions. The final complexes were minimized using the Polak-Ribiere conjugate gradient method using the AMBER99 force field . The value of the potential energy, visual analysis and the comparison among all the different docking runs allowed the selection of just one ligand-antibody complex. They were submitted, for further refinement, to a simulated annealing protocol using the molecular dynamic software NAMD .
2.7. Docking Refinement
The complexes obtained from molecular docking were submitted to a regular molecular dynamics simulation with NAMD  of 7.5 ns, at 300 K, allowing the whole system to move. This simulation did not show significant changes with respect to unrefined docking results.
The next step was to refine the complexes obtained from the molecular docking, using simulated annealing through molecular dynamics simulation with the software NAMD . The parameters of the simulation were optimized, through shorts runs, allowing all system parameters to move.
The optimized protocol was as follows: the system was solvated using a sphere of TIP3  water model. Spherical boundary conditions were applied to avoid the evaporation of water molecules. The motion of the antibody was restricted, allowing only the CDRs residues to move. A langevin piston [43, 44] was implemented to keep the pressure at 1 atm. The regulation of the temperature was enforced by a langevin thermostat . An integration time step of 2 fs was used. The simulation started at 300 K and the whole system was heated up to 600 K by 50 K jumps every 10 ps. Then, the temperature was held at 600 K during 120 ps. Afterwards, the temperature was decreased by 50 K jumps. The first jump was done after 10 ps, the second one after 20 ps, then 40 ps, 60 ps and 100 ps. Then, the system was kept at 300 K during 160 ps. The whole process was iterated three times.
3.1. Peptide structure
The structure of P013 can be briefly described as a bend in which the five negatively charged residues (E1, D3, D4, D5, D6) are pointing towards the solvent. The uncharged D2 is interacting with M9 which is in a cluster of residues that form a kind of hydrophobic core (F7, G8, M9, G10, A11, L12 and F13).
3.2. Conformational comparison of P013, R13, H13, A13 and rhodopsin
Experimental data demonstrated that the affinity of these peptides towards monoclonal antibody 17.2 follows the order R13 > H13 ≈ P013 > A13 [14, 17]. Remarkably, their conformational differences, which are reflected in the RMSD values of the structural alignment, followed an identical order suggesting that there may be a correlation between the conformation of the bend and its affinity. See figure 3.
Another important outcome was finding out that SEL13, a portion of the second extra cellular loop of rhodopsin, has a matching shape in spite of not having sequence identity with these peptides. This fact strongly reinforces the results of studies that have associated the dysfunction of the retina in chagasic patients with an autoimmune mechanism in which anti-P auto-antibodies would cross-react with the visual pigment rhodopsin . Our findings also support one of the conclusions of a previous study by Soares et al.  about the importance of the conformational state in the peptide-antibody interaction.
3.3. Structure of antibody mAb17.2 Fv domain and Antibody-peptide docking analysis
The docking of the R13 peptide into the antibody showed that the curve formed by E3, D4, D5 and D6 binds very deeply into the cavity while the residues forming the hydrophobic core (residues 7 to 13) were associated with the first hypervariable loop of the light chain.
Other observed electrostatic interactions were of the hydrogen bond type. E2 backbone carbonyl is hydrogen bonded to R165 (light chain) and G101 (heavy chain). E3 hydroxyl is hydrogen bonded to Y215 hydroxyl (light chain) D4 is forming and hydrogen bond with L152 and G210 (light chain).
The docking of the P013 peptide showed that the bend residues (3 to 6) bind similarly but less deeply than R13. D3 and 5 are again neutralized by R50 and 52 (heavy chain). The main difference observed when comparing electrostatic interactions between these two peptides is the ion pair formed between D4, which is unprotonated in P013, and R165 (light chain). The remaining residues (7 to 13) are exposed to the solvent instead of interacting with the CDR1 of the light chain as in the case of R13.
Another difference observed between the interactions of the two peptides with the antibody is the variation of the solvent accessible surface between the free and the docked state. The change in the solvent accessible surface was 10% for P013 and 13% for R13 after refinement. This parameter is a measure of the entropy increase in the complex formation and is given by the water molecule exclusion in the interaction surface.
The third docking study was about rhodopsin (this complex was not refined due to the size of the system). It was found that the residues PEVN of the second extra cellular loop bind in a similar fashion as to residues 3 to 6 in the R13-P013 complexes. Rhodopsin E197 is occupying the same position as E3 in the R13 complex. It is stabilized by the formation of two ions pairs: the first with R50 (heavy chain), remaining in the same conformation as it had in the previous complexes, and the second with R52 (heavy chain), which has moved from the position it had in the R13 complex being now in a conformation that allows the interaction with E197.
The first part of the study concerning the calculation of peptide structures by theoretical methods and their comparison with existing NMR and crystal structures reveals a most important common feature. Namely the curve formed by residues 3 to 6 having a surprising conformational similarity. Remarkably, the affinity of these peptides towards the antibody, following the order R13 > H13 ≈ P013 > A13 according to previous experimental work of Smulski et al  and Soares et al , is reflected in the order of the conformational differences of the peptides with respect to R13, according to the present work.
The outcome from this study clearly describes the interactions of the structural elements involved in the autoimmune mechanism of anti-P auto-antibodies cross-reaction, the stimulation of the β1-adrenoreceptor and the visual pigment rhodopsin. It also gives us some qualitative ideas to judge the potential binding of the anti-P auto-antibodies to receptor loops in different tissues. Results from this study may eventually lead to the development of treatments to abolish receptor mediated symptoms of Chagas's disease and other autoimmune pathologies.
This research was supported by grants from Agencia Nacional de Promoción Científica y Tecnológica (FONCyT), PICTR 2003-00286, FONCyT-ANPCyT (PAE 22642/22672), the National Research Council of Argentina (CONICET), and from the Secretaría de Ciencia y Técnica de la Universidad Nacional de San Luis, Argentina. This research was conducted using the resources of: two Beowulf-type clusters located at (a) the Instituto de Matemática Aplicada San Luis (CONICET-UNSL); and (b) the Baker Laboratory of Chemistry and Chemical Biology, Cornell University.
- World Health Organization on behalf of the Special Programme for Research and Training in Tropical Diseases: Edited by: Lazdins-Helds JK, Guhl F. 2007, http://apps.who.int/tdr/svc/publications/tdr-research-publications/reporte-enfermedad-chagas, ,
- Aznar C, López Bergami P, Brandariz S, Mariette C, Liegeard P, de Deus Alves MC, Luna Barreiro E, Carrasco R, Lafon S, Kaplan D, Levitus G, Levin MJ, Hontebeyrie-Joskowicz M: FEMS Immunol Med Microbiol. 1995, 12: 231-238.View ArticleGoogle Scholar
- Ferrari I, Levin MJ, Wallukat G, Elies R, Lebesgue D, Chiale P, Elizari M, Rosenbaum M, Hoebeke J: J Exp Med. 1995, 182: 59-65.View ArticleGoogle Scholar
- Liljas A: Int Rev Cytol. 1991, 124: 103-36.View ArticleGoogle Scholar
- Bailey-Serres J, Vangala S, Szick K, Lee CH: Plant Physiol. 1997, 114 (4): 1293-305.View ArticleGoogle Scholar
- Lin A, Wittman-Leibold B, McNelly T, Wool IG: J Biol Chem. 1982, 257: 9189-9197.Google Scholar
- Prieto J, Candel E, Fernandez-Renart M, Coloma A: Biochim Biophys Acta. 1991, 115: 6-14.View ArticleGoogle Scholar
- Newton CH, Shimmin LC, Yee J, Dennis PP: J Bacteriol. 1990, 172: 579-588.Google Scholar
- Beltrame M, Bianchi ME: Mol Cell Biol. 1990, 10: 2341-2348.View ArticleGoogle Scholar
- Schijman AG, Vazquez MP, Ben Dov C, Ghio S, Lorenzi H, Levin MJ: Biochim Biophys Acta. 1995, 1264: 15-18.View ArticleGoogle Scholar
- Lopez Bergami P, Scaglione J, Levin MJ: FASEB J. 2001, 15: 2602-2612.View ArticleGoogle Scholar
- Staudt Y, Mobini R, Fu M, Jens S, Kühn P, Staudt A: European Journal of Pharmacology. 2003, 466: 1-6.View ArticleGoogle Scholar
- Matsumoto SC, Labovsky V, Roncoroni M, Guida MC, Gimenez L, Mitelman J, Gori H, Jurgelevicius R, Grillo A, Manfredi P, Levin MJ, Paveto C: FASEB J. 2006, 20 (3): 550-2.Google Scholar
- Smulski C, Labovsky V, Levy G, Hontebeyrie M, Hoebeke J, Levin MJ: FASEB. 2006, 20 (9): 1396-406.View ArticleGoogle Scholar
- Ripoll DR, Scheraga HA: Biopolymers. 1988, 27: 1283-1303.View ArticleGoogle Scholar
- Ripoll DR, Liwo A, Czaplewski C: TASK Quart. 1999, 3: 313-331.Google Scholar
- Soares MR, Bisch PM, Campos de Carvalho AC, Valente AP, Almeida FC: FEBS Lett. 2004, 560 (1–3): 134-140.View ArticleGoogle Scholar
- Momany FA, McGuire RF, Burgess AW, Scheraga HA: J Phys Chem. 1975, 79: 2361-2381.View ArticleGoogle Scholar
- Némethy G, Pottle MS, Scheraga HA: J Phys Chem. 1983, 87: 1883-1887.View ArticleGoogle Scholar
- Sippl MJ, Némethy G, Scheraga HA: J Phys Chem. 1984, 88: 6231-6233.View ArticleGoogle Scholar
- Némethy G, Gibson KD, Palmer KA, Yoon CN, Paterlini G, Zagari A, Rumsey S, Scheraga HA: J Phys Chem. 1992, 96: 6472-6484.View ArticleGoogle Scholar
- Ripoll DR, Vorobjev YN, Liwo A, Vila JA, Scheraga HA: J Mol Biol. 1996, 264 (4): 770-783.View ArticleGoogle Scholar
- Connolly ML: J Appl Cryst. 1983, 16: 548-555.View ArticleGoogle Scholar
- Connolly ML: Science. 1983, 221: 709-713.ADSView ArticleGoogle Scholar
- Connolly ML: Appl Cryst. 1985, 18: 499-505.View ArticleGoogle Scholar
- Yang AS, Honig B: J Mol Biol. 1993, 231: 459-474.View ArticleGoogle Scholar
- Gilson MK: Proteins: Struc Funct Genet. 1993, 15: 266-282.View ArticleGoogle Scholar
- Vorobjev YN, Scheraga HA, Hitz B, Honig B: J Phys Chem. 1994, 98: 10940-10948.View ArticleGoogle Scholar
- You TJ, Bashford D: Biophys J. 1995, 69: 1721-1733.View ArticleGoogle Scholar
- Beroza P, Fredkin DR, Okamura MY, Feher G: Biophys J. 1995, 68: 2233-2250.View ArticleGoogle Scholar
- Vorobjev YN, Scheraga HA: J Comput Chem. 1997, 18: 569-583.View ArticleGoogle Scholar
- Vila J, Williams RL, Vasquez M, Scheraga HA: Proteins: Structure, Function, and Genetics. 1991, 10: 199-218.View ArticleGoogle Scholar
- Gô N, Scheraga HA: J Chem Phys. 1969, 51: 4751-4767.ADSView ArticleGoogle Scholar
- Zimmerman SS, Pottle MS, Némethy G, Scheraga HA: Macromolecules. 1977, 10: 1-9.ADSView ArticleGoogle Scholar
- Sali A, Blundell TL: J Mol Biol. 1993, 234: 779-815.View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Nucleic Acids Res. 1997, 25 (17): 3389-3402.View ArticleGoogle Scholar
- Notredame C, Higgins D, Heringa J: J Mol Biol. 2000, 302: 205-217.View ArticleGoogle Scholar
- Hooft RW, Vriend G, Sander C, Abola EE: Nature. 1996, 381: 272-272.ADSView ArticleGoogle Scholar
- Palma PN, Krippahl L, Wampler JE, Moura JJ: Proteins. 2000, 39 (4): 372-384.View ArticleGoogle Scholar
- Wang JM, Cieplak P, Kollman PA: J Comput Chem. 2000, 21: 1049-1074.View ArticleGoogle Scholar
- Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kale L, Schulten K: Journal of Computational Chemistry. 2005, 26: 1781-1802.View ArticleGoogle Scholar
- Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML: J Chem Phys. 1983, 79: 926-935.ADSView ArticleGoogle Scholar
- Feller SE, Zhang Y, Pastor RW, Brooks BR: J Chem Phys. 1995, 103: 4613-4621.ADSView ArticleGoogle Scholar
- Quigley D, Probert MIJ: J Chem Phys. 2004, 120: 11432-11442.ADSView ArticleGoogle Scholar
- Kubo R, Toda M, Hashitsume N: Statistical Physics II: Nonequilibrium Statistical Mechanics. 1991, New York: Springer, 2MATHView ArticleGoogle Scholar
- Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA: Proc Natl Acad Sci USA. 2001, 98: 10037-10041.ADSView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.