Structure-based Multiple Sequence Alignment Research Articles

SH2db, an information system for the SH2 domain.

SH2 domains are key mediators of phosphotyrosine-based signalling, and therapeutic targets for diverse, mostly oncological, disease indications. They have a highly conserved structure with a central beta sheet that divides the binding surface of the protein into two main pockets, responsible for phosphotyrosine binding (pY pocket) and substrate specificity (pY+3 pocket). In recent years, structural databases have proven to be invaluable resources for the drug discovery community, as they contain highly relevant and up-to-date information on important protein classes. Here, we present SH2db, a comprehensive structural database and webserver for SH2 domain structures. To organize these protein structures efficiently, we introduce (i) a generic residue numbering scheme to enhance the comparability of different SH2 domains, (ii) a structure-based multiple sequence alignment of all 120 human wild-type SH2 domain sequences and their PDB and AlphaFold structures. The aligned sequences and structures can be searched, browsed and downloaded from the online interface of SH2db (http://sh2db.ttk.hu), with functions to conveniently prepare multiple structures into a Pymol session, and to export simple charts on the contents of the database. Our hope is that SH2db can assist researchers in their day-to-day work by becoming a one-stop shop for SH2 domain related research.

Multidimensional Phylogenetic Metrics Identify Class I Aminoacyl-tRNA Synthetase Evolutionary Mosaicity and Inter-Modular Coupling.

The role of aminoacyl-tRNA synthetases (aaRS) in the emergence and evolution of genetic coding poses challenging questions concerning their provenance. We seek evidence about their ancestry from curated structure-based multiple sequence alignments of a structurally invariant “scaffold” shared by all 10 canonical Class I aaRS. Three uncorrelated phylogenetic metrics—mutation frequency, its uniformity, and row-by-row cladistic congruence—imply that the Class I scaffold is a mosaic assembled from successive genetic sources. Metrics for different modules vary in accordance with their presumed functionality. Sequences derived from the ATP– and amino acid– binding sites exhibit specific two-way coupling to those derived from Connecting Peptide 1, a third module whose metrics suggest later acquisition. The data help validate: (i) experimental fragmentations of the canonical Class I structure into three partitions that retain catalytic activities in proportion to their length; and (ii) evidence that the ancestral Class I aaRS gene also encoded a Class II ancestor in frame on the opposite strand. A 46-residue Class I “protozyme” roots the Class I tree prior to the adaptive radiation of the Rossmann dinucleotide binding fold that refined substrate discrimination. Such rooting implies near simultaneous emergence of genetic coding and the origin of the proteome, resolving a conundrum posed by previous inferences that Class I aaRS evolved after the genetic code had been implemented in an RNA world. Further, pinpointing discontinuous enhancements of aaRS fidelity establishes a timeline for the growth of coding from a binary amino acid alphabet.

Structural Insights into the Unique Activation Mechanisms of a Non-classical Calpain and Its Disease-Causing Variants

SUMMARYIncreased calpain activity is linked to neuroinflammation including a heritable retinal disease caused by hyper-activating mutations in the calcium-activated calpain-5 (CAPN5) protease. Although structures for classical calpains are known, the structure of CAPN5, a non-classical calpain, remains undetermined. Here we report the 2.8 Å crystal structure of the human CAPN5 protease core (CAPN5-PC). Compared to classical calpains, CAPN5-PC requires high calcium concentrations for maximal activity. Structure-based phylogenetic analysis and multiple sequence alignment reveal that CAPN5-PC contains three elongated flexible loops compared to its classical counterparts. The presence of a disease-causing mutation (c.799G>A, p.Gly267Ser) on the unique PC2L2 loop reveals a function in this region for regulating enzymatic activity. This mechanism could be transferred to distant calpains, using synthetic calpain hybrids, suggesting an evolutionary mechanism for fine-tuning calpain function by modifying flexible loops. Further, the open (inactive) conformation of CAPN5-PC provides structural insight into CAPN5-specific residues that can guide inhibitor design.

A Structurally-Validated Multiple Sequence Alignment of 497 Human Protein Kinase Domains

Studies on the structures and functions of individual kinases have been used to understand the biological properties of other kinases that do not yet have experimental structures. The key factor in accurate inference by homology is an accurate sequence alignment. We present a parsimonious, structure-based multiple sequence alignment (MSA) of 497 human protein kinase domains excluding atypical kinases. The alignment is arranged in 17 blocks of conserved regions and unaligned blocks in between that contain insertions of varying lengths present in only a subset of kinases. The aligned blocks contain well-conserved elements of secondary structure and well-known functional motifs, such as the DFG and HRD motifs. From pairwise, all-against-all alignment of 272 human kinase structures, we estimate the accuracy of our MSA to be 97%. The remaining inaccuracy comes from a few structures with shifted elements of secondary structure, and from the boundaries of aligned and unaligned regions, where compromises need to be made to encompass the majority of kinases. A new phylogeny of the protein kinase domains in the human genome based on our alignment indicates that ten kinases previously labeled as “OTHER” can be confidently placed into the CAMK group. These kinases comprise the Aurora kinases, Polo kinases, and calcium/calmodulin-dependent kinase kinases.

RStrucFam: a web server to associate structure and cognate RNA for RNA-binding proteins from sequence information.

BackgroundRNA-binding proteins (RBPs) interact with their cognate RNA(s) to form large biomolecular assemblies. They are versatile in their functionality and are involved in a myriad of processes inside the cell. RBPs with similar structural features and common biological functions are grouped together into families and superfamilies. It will be useful to obtain an early understanding and association of RNA-binding property of sequences of gene products. Here, we report a web server, RStrucFam, to predict the structure, type of cognate RNA(s) and function(s) of proteins, where possible, from mere sequence information.ResultsThe web server employs Hidden Markov Model scan (hmmscan) to enable association to a back-end database of structural and sequence families. The database (HMMRBP) comprises of 437 HMMs of RBP families of known structure that have been generated using structure-based sequence alignments and 746 sequence-centric RBP family HMMs. The input protein sequence is associated with structural or sequence domain families, if structure or sequence signatures exist. In case of association of the protein with a family of known structures, output features like, multiple structure-based sequence alignment (MSSA) of the query with all others members of that family is provided. Further, cognate RNA partner(s) for that protein, Gene Ontology (GO) annotations, if any and a homology model of the protein can be obtained. The users can also browse through the database for details pertaining to each family, protein or RNA and their related information based on keyword search or RNA motif search.ConclusionsRStrucFam is a web server that exploits structurally conserved features of RBPs, derived from known family members and imprinted in mathematical profiles, to predict putative RBPs from sequence information. Proteins that fail to associate with such structure-centric families are further queried against the sequence-centric RBP family HMMs in the HMMRBP database. Further, all other essential information pertaining to an RBP, like overall function annotations, are provided. The web server can be accessed at the following link: http://caps.ncbs.res.in/rstrucfam.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1289-x) contains supplementary material, which is available to authorized users.

In Silico Structure and Sequence Analysis of Bacterial Porins and Specific Diffusion Channels for Hydrophilic Molecules: Conservation, Multimericity and Multifunctionality.

Diffusion channels are involved in the selective uptake of nutrients and form the largest outer membrane protein (OMP) family in Gram-negative bacteria. Differences in pore size and amino acid composition contribute to the specificity. Structure-based multiple sequence alignments shed light on the structure-function relations for all eight subclasses. Entropy-variability analysis results are correlated to known structural and functional aspects, such as structural integrity, multimericity, specificity and biological niche adaptation. The high mutation rate in their surface-exposed loops is likely an important mechanism for host immune system evasion. Multiple sequence alignments for each subclass revealed conserved residue positions that are involved in substrate recognition and specificity. An analysis of monomeric protein channels revealed particular sequence patterns of amino acids that were observed in other classes at multimeric interfaces. This adds to the emerging evidence that all members of the family exist in a multimeric state. Our findings are important for understanding the role of members of this family in a wide range of bacterial processes, including bacterial food uptake, survival and adaptation mechanisms.

Structural Insights into a Novel Interkingdom Signaling Circuit by Cartography of the Ligand-Binding Sites of the Homologous Quorum Sensing LuxR-Family

Recent studies have identified a novel interkingdom signaling circuit, via plant signaling molecules, and a bacterial sub-family of LuxR proteins, bridging eukaryotes and prokaryotes. Indeed pivotal plant-bacteria interactions are regulated by the so called Plant Associated Bacteria (PAB) LuxR solo regulators that, although closely related to the quorum sensing (QS) LuxR family, do not bind or respond to canonical quorum sensing N-acyl homoserine lactones (AHLs), but only to specific host plant signal molecules. The large body of structural data available for several members of the QS LuxR family complexed with different classes of ligands (AHLs and other compounds), has been exploited to dissect the cartography of their regulatory domains through structure-based multiple sequence alignments, structural superimposition and a comparative analysis of the contact residues involved in ligand binding. In the absence of experimentally determined structures of members of the PAB LuxR solos subfamily, an homology model of its prototype OryR is presented, aiming to elucidate the architecture of its ligand-binding site. The obtained model, in combination with the cartography of the regulatory domains of the homologous QS LuxRs, provides novel insights into the 3D structure of its ligand-binding site and unveils the probable molecular determinants responsible for differences in selectivity towards specific host plant signal molecules, rather than to canonical QS compounds.

Analyses of the radiation of birnaviruses from diverse host phyla and of their evolutionary affinities with other double-stranded RNA and positive strand RNA viruses using robust structure-based multiple sequence alignments and advanced phylogenetic methods

BackgroundBirnaviruses form a distinct family of double-stranded RNA viruses infecting animals as different as vertebrates, mollusks, insects and rotifers. With such a wide host range, they constitute a good model for studying the adaptation to the host. Additionally, several lines of evidence link birnaviruses to positive strand RNA viruses and suggest that phylogenetic analyses may provide clues about transition.ResultsWe characterized the genome of a birnavirus from the rotifer Branchionus plicalitis. We used X-ray structures of RNA-dependent RNA polymerases and capsid proteins to obtain multiple structure alignments that allowed us to obtain reliable multiple sequence alignments and we employed “advanced” phylogenetic methods to study the evolutionary relationships between some positive strand and double-stranded RNA viruses. We showed that the rotifer birnavirus genome exhibited an organization remarkably similar to other birnaviruses. As this host was phylogenetically very distant from the other known species targeted by birnaviruses, we revisited the evolutionary pathways within the Birnaviridae family using phylogenetic reconstruction methods. We also applied a number of phylogenetic approaches based on structurally conserved domains/regions of the capsid and RNA-dependent RNA polymerase proteins to study the evolutionary relationships between birnaviruses, other double-stranded RNA viruses and positive strand RNA viruses.ConclusionsWe show that there is a good correlation between the phylogeny of the birnaviruses and that of their hosts at the phylum level using the RNA-dependent RNA polymerase (genomic segment B) on the one hand and a concatenation of the capsid protein, protease and ribonucleoprotein (genomic segment A) on the other hand. This correlation tends to vanish within phyla. The use of advanced phylogenetic methods and robust structure-based multiple sequence alignments allowed us to obtain a more accurate picture (in terms of probability of the tree topologies) of the evolutionary affinities between double-stranded RNA and positive strand RNA viruses. In particular, we were able to show that there exists a good statistical support for the claims that dsRNA viruses are not monophyletic and that viruses with permuted RdRps belong to a common evolution lineage as previously proposed by other groups. We also propose a tree topology with a good statistical support describing the evolutionary relationships between the Picornaviridae, Caliciviridae, Flaviviridae families and a group including the Alphatetraviridae, Nodaviridae, Permutotretraviridae, Birnaviridae, and Cystoviridae families.

Ribokinase family evolution and the role of conserved residues at the active site of the PfkB subfamily representative, Pfk-2 from Escherichia coli

Phosphofructokinase-2 (Pfk-2) belongs to the ribokinase family and catalyzes the ATP-dependent phosphorylation of fructose-6-phosphate, showing allosteric inhibition by a second ATP molecule. Several structures have been deposited on the PDB for this family of enzymes. A structure-based multiple sequence alignment of a non-redundant set of these proteins was used to infer phylogenetic relationships between family members with different specificities and to dissect between globally conserved positions and those common to phosphosugar kinases. We propose that phosphosugar kinases appeared early in the evolution of the ribokinase family. Also, we identified two conserved sequence motifs: the TR motif, not described previously, present in phosphosugar kinases but not in other members of the ribokinase family, and the globally conserved GXGD motif. Site-directed mutagenesis of R90 and D256 present in these motifs, indicate that R90 participates in the binding of the phosphorylated substrate and that D256 is involved in the phosphoryl transfer mechanism.

Relationships between functional subclasses and information contained in active-site and ligand-binding residues in diverse superfamilies

To investigate the relationships between functional subclasses and sequence and structural information contained in the active-site and ligand-binding residues (LBRs), we performed a detailed analysis of seven diverse enzyme superfamilies: aldolase class I, TIM-barrel glycosidases, alpha/beta-hydrolases, P-loop containing nucleotide triphosphate hydrolases, collagenase, Zn peptidases, and glutamine phosphoribosylpyrophosphate, subunit 1, domain 1. These homologous superfamilies, as defined in CATH, were selected from the enzyme catalytic-mechanism database. We defined active-site and LBRs based solely on the literature information and complex structures in the Protein Data Bank. From a structure-based multiple sequence alignment for each CATH homologous superfamily, we extracted subsequences consisting of the aligned positions that were used as an active-site or a ligand-binding site by at least one sequence. Using both the subsequences and full-length alignments, we performed cluster analysis with three sequence distance measures. We showed that the cluster analysis using the subsequences was able to detect functional subclasses more accurately than the clustering using the full-length alignments. The subsequences determined by only the literature information and complex structures, thus, had sufficient information to detect the functional subclasses. Detailed examination of the clustering results provided new insights into the mechanism of functional diversification for these superfamilies.

Structural basis for the substrate specificity of PepA from Streptococcus pneumoniae, a dodecameric tetrahedral protease

Regulated cytosolic proteolysis is one of the key cellular processes ensuring proper functioning of a cell. M42 family proteases show a broad spectrum of substrate specificities, but the structural basis for such diversity of the substrate specificities is lagging behind biochemical data. Here we report the crystal structure of PepA from Streptococcus pneumoniae, a glutamyl aminopeptidase belonging to M42 family (SpPepA). We found that Arg-257 in the substrate binding pocket is strategically positioned so that Arg-257 can make electrostatic interactions with the acidic residue of a substrate at its N-terminus. Structural comparison of the substrate binding pocket of the M42 family proteases, along with the structure-based multiple sequence alignment, argues that the appropriate electrostatic interactions contribute to the selective substrate specificity of SpPepA.

Solution Structure and Phylogenetics of Prod1, a Member of the Three-Finger Protein Superfamily Implicated in Salamander Limb Regeneration

BackgroundFollowing the amputation of a limb, newts and salamanders have the capability to regenerate the lost tissues via a complex process that takes place at the site of injury. Initially these cells undergo dedifferentiation to a state competent to regenerate the missing limb structures. Crucially, dedifferentiated cells have memory of their level of origin along the proximodistal (PD) axis of the limb, a property known as positional identity. Notophthalmus viridescens Prod1 is a cell-surface molecule of the three-finger protein (TFP) superfamily involved in the specification of newt limb PD identity. The TFP superfamily is a highly diverse group of metazoan proteins that includes snake venom toxins, mammalian transmembrane receptors and miscellaneous signaling molecules.Methodology/Principal FindingsWith the aim of identifying potential orthologs of Prod1, we have solved its 3D structure and compared it to other known TFPs using phylogenetic techniques. The analysis shows that TFP 3D structures group in different categories according to function. Prod1 clusters with other cell surface protein TFP domains including the complement regulator CD59 and the C-terminal domain of urokinase-type plasminogen activator. To infer orthology, a structure-based multiple sequence alignment of representative TFP family members was built and analyzed by phylogenetic methods. Prod1 has been proposed to be the salamander CD59 but our analysis fails to support this association. Prod1 is not a good match for any of the TFP families present in mammals and this result was further supported by the identification of the putative orthologs of both CD59 and N. viridescens Prod1 in sequence data for the salamander Ambystoma tigrinum.Conclusions/SignificanceThe available data suggest that Prod1, and thereby its role in encoding PD identity, is restricted to salamanders. The lack of comparable limb-regenerative capability in other adult vertebrates could be correlated with the absence of the Prod1 gene.

Correlated Mutation Analysis on the Catalytic Domains of Serine/Threonine Protein Kinases

BackgroundProtein kinases (PKs) have emerged as the largest family of signaling proteins in eukaryotic cells and are involved in every aspect of cellular regulation. Great progresses have been made in understanding the mechanisms of PKs phosphorylating their substrates, but the detailed mechanisms, by which PKs ensure their substrate specificity with their structurally conserved catalytic domains, still have not been adequately understood. Correlated mutation analysis based on large sets of diverse sequence data may provide new insights into this question.Methodology/Principal FindingsStatistical coupling, residue correlation and mutual information analyses along with clustering were applied to analyze the structure-based multiple sequence alignment of the catalytic domains of the Ser/Thr PK family. Two clusters of highly coupled sites were identified. Mapping these positions onto the 3D structure of PK catalytic domain showed that these two groups of positions form two physically close networks. We named these two networks as θ-shaped and γ-shaped networks, respectively.Conclusions/SignificanceThe θ-shaped network links the active site cleft and the substrate binding regions, and might participate in PKs recognizing and interacting with their substrates. The γ-shaped network is mainly situated in one side of substrate binding regions, linking the activation loop and the substrate binding regions. It might play a role in supporting the activation loop and substrate binding regions before catalysis, and participate in product releasing after phosphoryl transfer. Our results exhibit significant correlations with experimental observations, and can be used as a guide to further experimental and theoretical studies on the mechanisms of PKs interacting with their substrates.

Infernal 1.0: inference of RNA alignments

Abstract Summary: infernal builds consensus RNA secondary structure profiles called covariance models (CMs), and uses them to search nucleic acid sequence databases for homologous RNAs, or to create new sequence- and structure-based multiple sequence alignments. Availability: Source code, documentation and benchmark downloadable from http://infernal.janelia.org. infernal is freely licensed under the GNU GPLv3 and should be portable to any POSIX-compliant operating system, including Linux and Mac OS/X. Contact: nawrockie,kolbed,eddys@janelia.hhmi.org

Crystal structure of an alkaline serine protease from Nesterenkonia sp. defines a novel family of secreted bacterial proteases

Crystal structure of an alkaline serine protease from Nesterenkonia sp. defines a novel family of secreted bacterial proteases Na Yang, Jie Nan, Erik Brostromer, Rajni Hatti-Kaul, and Xiao-Dong Su* 1 Shenzhen Graduate School of Peking University, Shenzhen 518055, China 2 National Laboratory of Protein Engineering and Plant Genetic Engineering, Peking University, Beijing 100871, China 3 Department of Biotechnology, Center for Chemistry and Chemical Engineering, Lund University, S-221 00 Lund, Sweden

Algorithms for multiple protein structure alignment and structure-derived multiple sequence alignment.

Primary amino acid content and the geometry of the folded protein 3D structure are major parameters of protein function. During the course of evolution the protein 3D structure is more preserved than its primary sequence. Thus, analysis of protein structures is expected to lead to a deep insight into protein function. Recognition of a structural core common to a set of protein structures serves as a basic tool for the studies of protein evolution and classification, analysis of similar structural motifs and functional binding sites, and for homology modeling and threading. In this chapter, we discuss several biologically related computational aspects of the multiple structure alignment and propose a method that provides solutions to these problems. Finally, we address the problem of structure-based multiple sequence alignment and propose an optimization method that unifies primary sequence and 3D structure information.

Structure-guided identification of a new catalytic motif of oligosaccharyltransferase

Asn-glycosylation is widespread not only in eukaryotes but also in archaea and some eubacteria. Oligosaccharyltransferase (OST) catalyzes the co-translational transfer of an oligosaccharide from a lipid donor to an asparagine residue in nascent polypeptide chains. Here, we report that a thermophilic archaeon, Pyrococcus furiosus OST is composed of the STT3 protein alone, and catalyzes the transfer of a heptasaccharide, containing one hexouronate and two pentose residues, onto peptides in an Asn-X-Thr/Ser-motif-dependent manner. We also determined the 2.7-A resolution crystal structure of the C-terminal soluble domain of Pyrococcus STT3. The structure-based multiple sequence alignment revealed a new motif, DxxK, which is adjacent to the well-conserved WWDYG motif in the tertiary structure. The mutagenesis of the DK motif residues in yeast STT3 revealed the essential role of the motif in the catalytic activity. The function of this motif may be related to the binding of the pyrophosphate group of lipid-linked oligosaccharide donors through a transiently bound cation. Our structure provides the first structural insights into the formation of the oligosaccharide-asparagine bond.

Definition by Functional and Structural Analysis of Two Malonyl-CoA Sites in Carnitine Palmitoyltransferase 1A

Carnitine palmitoyltransferase 1 (CPT1) catalyzes the conversion of palmitoyl-CoA to palmitoylcarnitine in the presence of l-carnitine, thus facilitating the entry of fatty acids to mitochondria, in a process that is physiologically inhibited by malonyl-CoA. To examine the mechanism of CPT1 liver isoform (CPT1A) inhibition by malonyl-CoA, we constructed an in silico model of both its NH2- and COOH-terminal domains. Two malonyl-CoA binding sites were found. One of these, the "CoA site" or "A site," is involved in the interactions between NH2- and COOH-terminal domains and shares the acyl-CoA hemitunnel. The other, the "opposite-to-CoA site" or "O site," is on the opposite side of the enzyme, in the catalytic channel. The two sites share the carnitine-binding locus. To prevent the interaction between NH2- and COOH-terminal regions, we produced CPT1A E26K and K561E mutants. A double mutant E26K/K561E (swap), which was expected to conserve the interaction, was also produced. Inhibition assays showed a 12-fold decrease in the sensitivity (IC50) toward malonyl-CoA for CPT1A E26K and K561E single mutants, whereas swap mutant reverts to wild-type IC50 value. We conclude that structural interaction between both domains is critical for enzyme sensitivity to malonyl-CoA inhibition at the "A site." The location of the "O site" for malonyl-CoA binding was supported by inhibition assays of expressed R243T mutant. The model is also sustained by kinetic experiments that indicated linear mixed type malonyl-CoA inhibition for carnitine. Malonyl-CoA alters the affinity of carnitine, and there appears to be an exponential inverse relation between carnitine Km and malonyl-CoA IC50.

RILM: a web-based resource to aid comparative and functional analysis of the insulin and IGF-1 receptor family.

The metazoan receptors for insulin (INSR), insulin-like growth factor 1 (IGF1R), and other insulin-like molecules are transmembrane tyrosine kinases involved in the regulation of cell size, cell proliferation, development, signaling of nutritional and environmental conditions, and aging. Historically, mutations in the human insulin receptor have been studied because such changes often lead to severe insulin resistance. More recently, amino acid sequence alterations in the insulin receptor-like receptors of Drosophila melanogaster and Caenorhabditis elegans, as well as in the mouse insulin receptor have been the focus of attention. These modifications can have profound effects on growth, body size, metabolism, and aging. To integrate the many findings on insulin/IGF1 receptor structure and function across species we have created "Receptors for Insulin and Insulin-like Molecules" (RILM), a curated computer-based resource that displays residue-by-residue information on sequence homology, three-dimensional structure, structure/function annotation, and documented mutations. The resource includes data obtained from sequence and structure analysis tools, primary database resources, and published reports. The information is integrated via a structure-based multiple sequence alignment of diverse members of the family. RILM was designed to provide easy access to multiple data types that could prove useful in the analysis of the effect of mutations on protein structure and ligand binding within this receptor family. RILM is available at www.biochem.ucl.ac.uk/RILM.

The Solution Structure of Escherichia coli Wzb Reveals a Novel Substrate Recognition Mechanism of Prokaryotic Low Molecular Weight Protein-tyrosine Phosphatases

Low molecular weight protein-tyrosine phosphatases (LMW-PTPs) are small enzymes that ubiquitously exist in various organisms and play important roles in many biological processes. In Escherichia coli, the LMW-PTP Wzb dephosphorylates the autokinase Wzc, and the Wzc/Wzb pair regulates colanic acid production. However, the substrate recognition mechanism of Wzb is still poorly understood thus far. To elucidate the molecular basis of the catalytic mechanism, we have determined the solution structure of Wzb at high resolution by NMR spectroscopy. The Wzb structure highly resembles that of the typical LMW-PTP fold, suggesting that Wzb may adopt a similar catalytic mechanism with other LMW-PTPs. Nevertheless, in comparison with eukaryotic LMW-PTPs, the absence of an aromatic amino acid at the bottom of the active site significantly alters the molecular surface and implicates Wzb may adopt a novel substrate recognition mechanism. Furthermore, a structure-based multiple sequence alignment suggests that a class of the prokaryotic LMW-PTPs may share a similar substrate recognition mechanism with Wzb. The current studies provide the structural basis for rational drug design against the pathogenic bacteria.