Bi-Domain Protein Tyrosine Phosphatases Reveal an Evolutionary Adaptation to Optimize Signal Transduction

Significance: The bi-domain protein tyrosine phosphatases (PTPs) exemplify functional evolution in signaling proteins for optimal spatiotemporal signal transduction. Bi-domain PTPs are products of gene duplication. The catalytic activity, however, is often localized to one PTP domain. The inactive PTP domain adopts multiple functional roles. These include modulation of catalytic activity, substrate specificity, and stability of the bi-domain enzyme. In some cases, the inactive PTP domain is a receptor for redox stimuli. Since multiple bi-domain PTPs are concurrently active in related cellular pathways, a stringent regulatory mechanism and selective cross-talk is essential to ensure fidelity in signal transduction. Recent Advances: The inactive PTP domain is an activator for the catalytic PTP domain in some cases, whereas it reduces catalytic activity in other bi-domain PTPs. The relative orientation of the two domains provides a conformational rationale for this regulatory mechanism. Recent structural and biochemical data reveal that these PTP domains participate in substrate recruitment. The inactive PTP domain has also been demonstrated to undergo substantial conformational rearrangement and oligomerization under oxidative stress. Critical Issues and Future Directions: The role of the inactive PTP domain in coupling environmental stimuli with catalytic activity needs to be further examined. Another aspect that merits attention is the role of this domain in substrate recruitment. These aspects have been poorly characterized in vivo. These lacunae currently restrict our understanding of neo-functionalization of the inactive PTP domain in the bi-domain enzyme. It appears likely that more data from these research themes could form the basis for understanding the fidelity in intracellular signal transduction.

Thermodynamic characterization of the reversible, two-state unfolding of maltose binding protein, a large two-domain protein

The folding and stability of maltose binding protein (MBP) have been investigated as a function of pH and temperature by intrinsic tryptophan fluorescence, far- and near-UV circular dichroism, and high-sensitivity differential scanning calorimetric measurements. MBP is a monomeric, two-domain protein containing 370 amino acids. The protein is stable in the pH range of 4-10.5 at 25 degrees C. The protein exhibits reversible, two-state, thermal and guanidine hydrochloride-mediated denaturation at neutral pH. The thermostability of MBP is maximal at pH 6, with a Tm of 64.9 degrees C and a deltaHm of 259.7 kcal mol(-1). The linear dependence of deltaHm on Tm was used to estimate a value of deltaCp of 7.9 kcal mol(-1) K(-1) or 21.3 cal (mol of residue)(-1) K(-1). These values are higher than the corresponding deltaCp's for most globular proteins studied to date. However, the extrapolated values of deltaH and deltaS (per mole of residue) at 110 degrees C are similar to those of other globular proteins. These data have been used to show that the temperature at which a protein undergoes cold denaturation depends primarily on the deltaCp (per mol of residue) and that this temperature increases with an increase in deltaCp. The predicted decrease in stability of MBP at low temperatures was experimentally confirmed by carrying out denaturant-mediated unfolding studies at neutral pH at 2 and 28 degrees C.

Modulation of Catalytic Activity in Multi-Domain Protein Tyrosine Phosphatases

Signaling mechanisms involving protein tyrosine phosphatases govern several cellular and developmental processes. These enzymes are regulated by several mechanisms which include variation in the catalytic turnover rate based on redox stimuli, subcellular localization or protein-protein interactions. In the case of Receptor Protein Tyrosine Phosphatases (RPTPs) containing two PTP domains, phosphatase activity is localized in their membrane-proximal (D1) domains, while the membrane-distal (D2) domain is believed to play a modulatory role. Here we report our analysis of the influence of the D2 domain on the catalytic activity and substrate specificity of the D1 domain using two Drosophila melanogaster RPTPs as a model system. Biochemical studies reveal contrasting roles for the D2 domain of Drosophila Leukocyte antigen Related (DLAR) and Protein Tyrosine Phosphatase on Drosophila chromosome band 99A (PTP99A). While D2 lowers the catalytic activity of the D1 domain in DLAR, the D2 domain of PTP99A leads to an increase in the catalytic activity of its D1 domain. Substrate specificity, on the other hand, is cumulative, whereby the individual specificities of the D1 and D2 domains contribute to the substrate specificity of these two-domain enzymes. Molecular dynamics simulations on structural models of DLAR and PTP99A reveal a conformational rationale for the experimental observations. These studies reveal that concerted structural changes mediate inter-domain communication resulting in either inhibitory or activating effects of the membrane distal PTP domain on the catalytic activity of the membrane proximal PTP domain.

Conformational Basis for Substrate Recruitment in Protein Tyrosine Phosphatase 10D

The coordinated activity of protein tyrosine phosphatases (PTPs) is crucial for the initiation, modulation, and termination of diverse cellular processes. The catalytic activity of this protein depends on a nucleophilic cysteine at the active site that mediates the hydrolysis of the incoming phosphotyrosine substrate. While the role of conserved residues in the catalytic mechanism of PTPs has been extensively examined, the diversity in the mechanisms of substrate recognition and modulation of catalytic activity suggests that other, less conserved sequence and structural features could contribute to this process. Here we describe the crystal structures of Drosophila melanogaster PTP10D in the apo form as well as in a complex with a substrate peptide and an inhibitor. These studies reveal the role of aromatic ring stacking interactions at the boundary of the active site of PTPs in mediating substrate recruitment. We note that phenylalanine 76, of the so-called KNRY loop, is crucial for orienting the phosphotyrosine residue toward the nucleophilic cysteine. Mutation of phenylalanine 76 to leucine results in a 60-fold decrease in the catalytic efficiency of the enzyme. Fluorescence measurements with a competitive inhibitor, p-nitrocatechol sulfate, suggest that Phe76 also influences the formation of the enzyme-substrate intermediate. The structural and biochemical data for PTP10D thus highlight the role of relatively less conserved residues in PTP domains in both substrate recruitment and modulation of reaction kinetics.

Inter-domain interactions influence the stability and catalytic activity of the bi-domain protein tyrosine phosphatase PTP99A

The two protein tyrosine phosphatase (PTP) domains in bi-domain PTPs share high sequence and structural similarity. However, only one of the two PIP domains is catalytically active. Here we describe biochemical studies on the two tandem PTP domains of the bi-domain PTP, PTP99A. Phosphatase activity, monitored using small molecule as well as peptide substrates, revealed that the inactive (D2) domain activates the catalytic (D1) domain. Thermodynamic measurements suggest that the inactive D2 domain stabilizes the bi-domain (D1-D2) protein. The mechanism by which the D2 domain activates and stabilizes the bi-domain protein is governed by few interactions at the inter-domain interface. In particular, mutating Lys990 at the interface attenuates inter-domain communication. This residue is located at a structurally equivalent location to the so-called allosteric site of the canonical single domain PIP, PTP1B. These observations suggest functional optimization in bi-domain PTPs whereby the inactive PTP domain modulates the catalytic activity of the bi-domain enzyme. (C) 2012 Elsevier B.V. All rights reserved.

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

Role of a PAS sensor domain in the Mycobacterium tuberculosis transcription regulator Rv1364c

The Mycobacterium tuberculosis transcriptional regulator Rv1364c regulates the activity of the stress response sigma factor sigma(F). This multi-domain protein has several components: a signaling PAS domain and an effector segment comprising of a phosphatase, a kinase and an anti-anti-sigma factor domain. Based on Small Angle X-ray Scattering (SAXS) data, Rv1364c was recently shown to be a homo-dimer and adopt an elongated conformation in solution. The PAS domain could not be modeled into the structural envelope due to poor sequence similarity with known PAS proteins. The crystal structure of the PAS domain described here provides a structural basis for the dimerization of Rv1364c. It thus appears likely that the PAS domain regulates the anti-sigma activity of Rv1364c by oligomerization. A structural comparison with other characterized PAS domains reveal several sequence and conformational features that could facilitate ligand binding - a feature which suggests that the function of Rv1364c could potentially be governed by specific cellular signals or metabolic cues. (C) 2010 Elsevier Inc. All rights reserved.

Crystal structures, dynamics and functional implications of molybdenum-cofactor biosynthesis protein MogA from two thermophilic organisms

Molybdenum-cofactor (Moco) biosynthesis is an evolutionarily conserved pathway in almost all kingdoms of life, including humans. Two proteins, MogA and MoeA, catalyze the last step of this pathway in bacteria, whereas a single two-domain protein carries out catalysis in eukaryotes. Here, three crystal structures of the Moco-biosynthesis protein MogA from the two thermophilic organisms Thermus thermophilus (TtMogA; 1.64 angstrom resolution, space group P2(1)) and Aquifex aeolicus (AaMogA; 1.70 angstrom resolution, space group P2(1) and 1.90 angstrom resolution, space group P1) have been determined. The functional roles and the residues involved in oligomerization of the protein molecules have been identified based on a comparative analysis of these structures with those of homologous proteins. Furthermore, functional roles have been proposed for the N- and C-terminal residues. In addition, a possible protein-protein complex of MogA and MoeA has been proposed and the residues involved in protein-protein interactions are discussed. Several invariant water molecules and those present at the subunit interfaces have been identified and their possible structural and/or functional roles are described in brief. In addition, molecular-dynamics and docking studies with several small molecules (including the substrate and the product) have been carried out in order to estimate their binding affinities towards AaMogA and TtMogA. The results obtained are further compared with those obtained for homologous eukaryotic proteins.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.

Synergy between the N-terminal and C-terminal domains of Mycobacterium tuberculosis HupB is essential for high-affinity binding, DNA supercoiling and inhibition of RecA-promoted strand exchange

The occurrence of DNA architectural proteins containing two functional domains derived from two different architectural proteins is an interesting emerging research theme in the field of nucleoid structure and function. Mycobacterium tuberculosis HupB, unlike Escherichia coli HU, is a two-domain protein that, in the N-terminal region, shows broad sequence homology with bacterial HU. The long C-terminal extension, on the other hand, contains seven PAKK/KAAK motifs, which are characteristic of the histone H1/H5 family of proteins. In this article, we describe several aspects of HupB function, in comparison with its truncated derivatives lacking either the C-terminus or N-terminus. We found that HupB binds a variety of DNA repair and replication intermediates with K(d) values in the nanomolar range. By contrast, the N-terminal fragment of M. tuberculosis HupB (HupB(MtbN)) showed diminished DNA-binding activity, with K(d) values in the micromolar range, and the C-terminal domain was completely devoid of DNA-binding activity. Unlike HupB(MtbN), HupB was able to constrain DNA in negative supercoils and introduce negative superhelical turns into relaxed DNA. Similarly, HupB exerted a robust inhibitory effect on DNA strand exchange promoted by cognate and noncognate RecA proteins, whereas HupB(MtbN), even at a 50-fold molar excess, had no inhibitory effect. Considered together, these results suggest that synergy between the N-terminal and C-terminal domains of HupB is essential for its DNA-binding ability, and to modulate the topological features of DNA, which has implications for processes such as DNA compaction, gene regulation, homologous recombination, and DNA repair.

Structural plasticity in Mycobacterium tuberculosis uracil-DNA glycosylase (MtUng) and its functional implications

17 independent crystal structures of family I uracil-DNA glycosylase from Mycobacterium tuberculosis (MtUng) and its complexes with uracil and its derivatives, distributed among five distinct crystal forms, have been determined. Thermodynamic parameters of binding in the complexes have been measured using isothermal titration calorimetry. The two-domain protein exhibits open and closed conformations, suggesting that the closure of the domain on DNA binding involves conformational selection. Segmental mobility in the enzyme molecule is confined to a 32-residue stretch which plays a major role in DNA binding. Uracil and its derivatives can bind to the protein in two possible orientations. Only one of them is possible when there is a bulky substituent at the 50 position. The crystal structures of the complexes provide a reasonable rationale for the observed thermodynamic parameters. In addition to providing fresh insights into the structure, plasticity and interactions of the protein molecule, the results of the present investigation provide a platform for structure-based inhibitor design.

Evolution of domain combinations in protein kinases and its implications for functional diversity

Protein kinases phosphorylating Ser/Thr/Tyr residues in several cellular proteins exert tight control over their biological functions. They constitute the largest protein family in most eukaryotic species. Protein kinases classified based on sequence similarity in their catalytic domains, cluster into subfamilies, which share gross functional properties. Many protein kinases are associated or tethered covalently to domains that serve as adapter or regulatory modules,naiding substrate recruitment, specificity, and also serve as scaffolds. Hence the modular organisation of the protein kinases serves as guidelines to their functional and molecular properties. Analysis of genomic repertoires of protein kinases in eukaryotes have revealed wide spectrum of domain organisation across various subfamilies of kinases. Occurrence of organism-specific novel domain combinations suggests functional diversity achieved by protein kinases in order to regulate variety of biological processes. In addition, domain architecture of protein kinases revealed existence of hybrid protein kinase subfamilies and their emerging roles in the signaling of eukaryotic organisms. In this review we discuss the repertoire of non-kinase domains tethered to multi-domain kinases in the metazoans. Similarities and differences in the domain architectures of protein kinases in these organisms indicate conserved and unique features that are critical to functional specialization. (C) 2009 Elsevier Ltd. All rights reserved.