Identification of distant co-evolving residues in antigen 85C from Mycobacterium tuberculosis using statistical coupling analysis of the esterase family proteins

Veeky Baths, Utpal Roy

Department of Biological Sciences, BITS Pilani K. K Birla Goa Campus, GOA 403726, India.

INTRODUCTION

Communication between distant sites in proteins is fundamental to their function and often defines the biological role of a protein family. In signaling proteins,it represents information transfer -- the transmission of signals initiated at one functional surface to a distinct surface mediating downstream signaling. For example, ligand binding at an externally accessible site in G protein-coupled receptors (GPCRs)reliably triggers structural changes at distant cytoplasmic domains that mediate interaction with heterotrimeric G proteins[1,2].Studies in many other protein systems indicate that long-range interactions of amino acids also are important in binding (and catalytic)specificity. Substrate recognition in the chymotrypsin family of serine proteases[3,4], the tuning of antibody specificity through B-cell maturation[5]and the cooperativity of oxygen binding in hemoglobin[6]all depend not only on residues directly contacting the substrate, but also on distant residues located in supporting loops and other secondary structural elements.

Statistical coupling analysis (SCA)is a technique used to identify communication between distant sites in proteins. More specifically, it quantifies how much the amino acid distribution at some position i changes upon a perturbation of the amino acid distribution at another position j. The resulting statistical coupling energy indicates the degree of evolutionary dependence between the residues, with higher coupling energy corresponding to increased dependence. Previous application of SCA done on the PDZ domain family [(PDZ is an acronym combining the first letters of three proteins — post synaptic density protein(PSD95), Drosophila disc large tumor suppressor(DlgA), and zonula occludens-1 protein (zo-1)]. These PDZ domain family proteins were first discovered to share the domain[7]. Apart from sharing of domains,this protein family has also got a predicted set of energetically coupled positions for the binding site residue that contained unexpected long-range interactions.Further mutational analysis confirmed the prediction that statistical energy function is a good indicator for thermodynamic coupling in proteins[8]. Application of SCA to three structurally and functionally distinct protein families (GPCR, the chymotrypsin class of serine proteases and hemoglobin)revealed a simple architecture for amino acid interaction in protein families and link distant functional sites in the tertiary structure[9]. Further application to the S1A serine protease family indicated the presence of quasi-independent groups of correlated amino acids termed as"Protein Sectors". Each of these sectors was found to be physically connected in the tertiary structure, had a distinct functional role and constituted an independent mode of sequence divergence in the protein family[10].

The antigen 85 (ag85)complex, composed of three proteins (ag85 A, B and C), is a major component of the Mycobacterium (M.)tuberculosis cell wall. Each protein possesses a mycolyltransferase activity required for the biogenesis of trehalose dimycolate (cord factor), a dominant structure necessary for maintaining cell wall integrity[11,12]. The protein sequence of ag85 C, which is a member of the esterase family proteins, seems to be responsible for the high affinity of Mycobacterium to fibronectin[12]. Since ag85 proteins are important for cell wall biosynthesis, it could be a target for a novel antitubercular drug. With the aim of identifying residue positions, which could be responsible for the function of ag85 C and hence act as potential drug targets for M. tuberculosis, we applied SCA to the esterase protein family.

MATERIALS AND METHODS

Methods

Datasets

The putative esterase family (pfam ID: PF00756 alignment (of 4945 sequences)was downloaded (updated till April 2010)from the pfam database (http://pfam.sanger.ac.uk/)[10]. Sequence redundancies were removed from the database at the threshold value of 99% using JalView[14](846 out of 4945 sequences left). The dataset was realigned using ClustalX and no further manual adjustments were done[15]. The dataset was further refined by the removal of all columns with more than 20% gaps (using MATLAB)(149 columns left). This prevented any trivial over-representation of gaps in the alignment and ensured that the calculations were only made at largely non-gapped sequences positions. SCA was performed using the code provided by Ranganathan[10]. Clustering analysis was performed in the form of hierarchical clustering analysis by construction of dendrograms of the correlation matrix using MATLAB. Coupled residues were mapped on to the structure of the secreted form of ag85 C and the crystal structure of the secreted form of ag85 C (PDB ID: 1DQZ)was available on the PDB database for structural analysis.

RESULTS

Fig. 1 shows the positional conservation of all the 149 residue positions (left after removal of redundancy and gap reduction on the original downloaded MSA. We call the new alignment: "truncated alignment")in the truncated alignment of the esterase family. It also shows the distribution of Di(a)values of all the positions as a measure of conservation.

Histogram of pairwise correlation values of the 149 residue positions as shown in Fig. 1 indicates that most of the values (nearly 99%)lie close to the mean value of the histogram and only a few pairs show significant correlation and only about 1% of the total pairs show significant correlation (values greater than eXo+2.36σ).This indicates that only a few pairs of residues correlate with each other with high correlation values. This essentially explains that in putative esterase protein family only a few residue pairs show interaction and communicate with each other even from a long range.According to the CASP guidelines (http://www.predictioncenter.llnl-.gov/CASP9), a "long-range" interresidue contact is defined as two residues that are separated by at least nine residue positions in the linear sequence in which Cβ-Cβ distance is ≤ 8 Å[16]. Fig. 2 shows the correlation values.

Fig. 1 Sequence conservation in the esterase family multiple sequence alignment. The degree of sequence conservation is plotted along the protein sequence with high values indicating high conservation. The plot corresponds to the diagonal elements.

Fig. 2 Histogram of pairwise correlation values. The figure suggests that only a few pairs correlate with a significant correlation where as most of the correlation values (99%)lie around the mean value of 0.040978.

Fig. 3 Heat map of the correlation matrix. The red pixels represent high correlation values and as we scale down the blue pixels represent low correlation values.

For a correlation matrix M, the value M(i,j)would mean the correlation value between the ith and jth residue in the multiple sequence alignment (MSA). The heat map in Fig. 3 is a representation of this correlation matrix for the esterase family alignments where most of the pixels are blue to indicate that most of the correlation values are significantly low. However, there are few red spots, highlighting the significant correlations. The diagonal values represent auto-correlation values, which is a representation of the conservation of that particular position. In order to identify the pairs of residues with high correlation values, we performed hierarchical clustering. This is a method of cluster analysis, which seeks to build a hierarchy of clusters.Hierarchical clustering enables grouping of the data over a variety of scales by creating a cluster tree or dendrogram. The tree is not a single set of clusters,but rather a multilevel hierarchy, where clusters at one level are joined as clusters to the next level[17].Fig. 4 represents the result obtained by performing hierarchical clustering on the pair-wise correlation.Values obtained from the SCA of the esterase family alignments as can be seen in Fig. 4 where a cluster of 7 residues (in green)was identified. Positions 22, 47,67, 92, 94, 121 and 146 from the truncated alignment formed the cluster of closely related residue position,which correlated with a high correlation value. This essentially means that these residue positions are statistically correlated to one other with high correlation value. Our hypothesis is that these are significant to the esterase family protein functions.

Fig. 4 Dendogram representing clustering analysis of pairwise correlation values. The cluster represented in green is the cluster with the most significantly coupled residue positions. X-axis represents the residue position in the truncated alignment and Y-axis represents the clustering distance between any two-residue positions. The arrow indicates a cluster of seven residues.

Esterase family alignments contain the sequence of ag85 C secreted by M. tuberculosis. We mapped these residue positions from the truncated alignment onto the crystal structure of the secreted form of ag85 C manually. Fig. 5 represents the mapping of seven residues. There are a total of 8-282 residues, the apo form of ag85 C (PDB accession 1DQZ)

As can be seen in Fig. 5, the residues that have been mapped are Asp 38, Val 66, Phe 98, Leu 123, Met 125, Trp 186 and Ser 215. A close look at the structural positions of these residues highlights some interesting and important information about these residues and ag85 C. Three of the seven residues (shown in red in Fig. 5)seem to be physically interacting in the form of electrostatic interactions, and further Asp38 and Met125, Asp38 and Leu123 also qualify for the CASP definition of long-range residues[16]. Fig. 6 shows other esterase.

Fig. 5 Mapping of the 7 residue positions highlighted by SCA onto the crystal structure of antigen 85 C(PDB ID: 1DQZ). The residues shown are Asp38, Val66,Phe98, Leu123, Met125, Trp186 and Ser215. The pairs in red Asp38, Leu123 and Asp38, Met125 seem to show physical interaction through allostery.

Fig. 6 Structures of the Mycobacterium tuberculosis 30 kDa major secretory protein (Antigen 85B). A mycolyl transferase (PDB ID: 1FOP), another esterase showing similar pocket.

DISCUSSION

In this work, we have presented the SCA of MSA of putative esterase family of proteins, which is known to contain several seemingly unrelated proteins. Three(Asp38, Leu123 and Met125)out of seven residues as identified by SCA seem to interact physically as seen from the mapping on the crystal structure of ag85 C(Fig. 5). The pair Arg 38 and Met 125 has been previously reported to be important for the function of ag85 C. Arg38 is one of the residues, which form a pocket with a high negative electrostatic potential. Ag85 C when complexed with a covalent inhibitor implicates residues Leu40 (close to Arg38)and Met125 as components of oxyanion hole[10]. Mutating such a pair may result in the non-functioning of ag85 C.The functioning of ag85 C is crucial for the survival of M. tuberculosis in host environment[18]. The results of this investigation suggest that Arg38 and Met125 can be potential targets for antituburcular drugs, as allosteric changes in the structure of these two residues may lead to the non-functioning of ag85 C and hence inactivation of the organism. Active site of ag85 C illustrates the binding mode of the substrate and extends the knowledge concerning specific protein/substrate interactions. This is the first step toward designing potent inhibitors to the mycolyltransferase activity of the ag85 enzymes[13].

[1]Gether U. Uncovering molecular mechanisms involved in activation of G protein coupled receptors. Endocr Rev 2000;21,90-113.

[2]Menon ST, Han M, Sakmar TP. Rhodospin: structural basis of molecular physiology, Physiol Rev 2001;81:1659-88.

[3]Hedstrom L, Szilagyi L, Rutter WJ. Converting trypsin to chymotrypsin: the role of surface loops. Science 1993;255:1249-53.

[4]Hedstrom L. Trypsin: a case study in the structural determinants of enzyme specificity. Biol Chem 1996;377:465-70.

[5]Patten PA, Gray NS, Yang PL, Marks CB, Wedemayer GY, Jay Boniface J. The immunological evolution of catalysis. Science 1996;271:1086-91.

[6]Pertuz MF. Stereochemistry of cooperative effects in haemoglobin. Nature 1970;228:726-39.

[7]Ponting CP, Phillips C, Davies KE, Blake DJ (1997).PDZ domains: targeting signaling molecules to submembranous sites. Bioessays 1997;(6):469-79.

[8]Lockless SW, Ranganathan R. Evolutionarily conserved pathways of energetic connectivity in protein families.Science 1999;286:295-9.

[9]Suel GM, Lockless SW, Wall MA, Ranganathan R.Evolutionarily conserved networks of residues mediate allosteric communication in proteins. Nat Struct Biol 2003;10:59-69.

[10]Halabi N, Rivoire O, Leibler S, Ranganathan R. Protein sectors: Evolutionary units of three-dimensional structure. Cell 2009;138:774-86.

[11]Ronning DR, Klabunde T, Besra GS, Vissa VD, Belisle JT, Sacchettini JC. Crystal structure of the secreted form of antigen 85C reveals potential drug targets for mycobacterial drugs and vaccines. Nat Struct Biol 2000;7:141-6.

[12]Belisle JT, Vissa VD, Sievert T, Takayama K, Brennan PJ, Besra GS. Role of major antigen of Mycobacterium tuberculosis in cell wall biogenesis. Science 1997;276:1420-2.

[13]Ronning DR, Vissa V, Belisle JT, Sacchettini JC. Mycobacterium tuberculosis antigen 85A and 85C structures confirm binding orientations and conserved substrate specificity. J Biol Chem 2004;279:771-7.

[14]Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, et al. The Pfam protein families database. Nucleic Acids Research Database Issue 2010;38:D211-22.

[15]Waterhouse AM, Procter JB, Martin DMA, Clamp M,Barton GJ. Jalview Version 2 - a multiple sequence alignment editor and analysis workbench. Bioinformatics 2009;25:1189-91.

[16]Larkin MA, Blackshields G, Brown NP, Chenna R, Mc-Gettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007;23:2947-8.

[17]Anthony A, Fodor Richard W, Aldrich. Influence of conservation on calculations of amino acid covariance in multiple sequence alignments. Proteins 2004;56:211-21.

[18]Ward Joe H. Hierarchical grouping to optimize an objective function. J Am Stat Assoc 1963;58:236-44.