PeptideMine - A webserver for the design of peptides for protein-peptide binding studies derived from protein-protein interactomes

Background: Signal transduction events often involve transient, yet specific, interactions between structurally conserved protein domains and polypeptide sequences in target proteins. The identification and validation of these associating domains is crucial to understand signal transduction pathways that modulate different cellular or developmental processes. Bioinformatics strategies to extract and integrate information from diverse sources have been shown to facilitate the experimental design to understand complex biological events. These methods, primarily based on information from high-throughput experiments, have also led to the identification of new connections thus providing hypothetical models for cellular events. Such models, in turn, provide a framework for directing experimental efforts for validating the predicted molecular rationale for complex cellular processes. In this context, it is envisaged that the rational design of peptides for protein-peptide binding studies could substantially facilitate the experimental strategies to evaluate a predicted interaction. This rational design procedure involves the integration of protein-protein interaction data, gene ontology, physico-chemical calculations, domain-domain interaction data and information on functional sites or critical residues. Results: Here we describe an integrated approach called ``PeptideMine'' for the identification of peptides based on specific functional patterns present in the sequence of an interacting protein. This approach based on sequence searches in the interacting sequence space has been developed into a webserver, which can be used for the identification and analysis of peptides, peptide homologues or functional patterns from the interacting sequence space of a protein. To further facilitate experimental validation, the PeptideMine webserver also provides a list of physico-chemical parameters corresponding to the peptide to determine the feasibility of using the peptide for in vitro biochemical or biophysical studies. Conclusions: The strategy described here involves the integration of data and tools to identify potential interacting partners for a protein and design criteria for peptides based on desired biochemical properties. Alongside the search for interacting protein sequences using three different search programs, the server also provides the biochemical characteristics of candidate peptides to prune peptide sequences based on features that are most suited for a given experiment. The PeptideMine server is available at the URL: http://caps.ncbs.res.in/peptidemine

The C-Terminal Domain of the MutL Homolog from Neisseria gonorrhoeae Forms an Inverted Homodimer

The mismatch repair (MMR) pathway serves to maintain the integrity of the genome by removing mispaired bases from the newly synthesized strand. In E. coli, MutS, MutL and MutH coordinate to discriminate the daughter strand through a mechanism involving lack of methylation on the new strand. This facilitates the creation of a nick by MutH in the daughter strand to initiate mismatch repair. Many bacteria and eukaryotes, including humans, do not possess a homolog of MutH. Although the exact strategy for strand discrimination in these organisms is yet to be ascertained, the required nicking endonuclease activity is resident in the C-terminal domain of MutL. This activity is dependent on the integrity of a conserved metal binding motif. Unlike their eukaryotic counterparts, MutL in bacteria like Neisseria exist in the form of a homodimer. Even though this homodimer would possess two active sites, it still acts a nicking endonuclease. Here, we present the crystal structure of the C-terminal domain (CTD) of the MutL homolog of Neisseria gonorrhoeae (NgoL) determined to a resolution of 2.4 A. The structure shows that the metal binding motif exists in a helical configuration and that four of the six conserved motifs in the MutL family, including the metal binding site, localize together to form a composite active site. NgoL-CTD exists in the form of an elongated inverted homodimer stabilized by a hydrophobic interface rich in leucines. The inverted arrangement places the two composite active sites in each subunit on opposite lateral sides of the homodimer. Such an arrangement raises the possibility that one of the active sites is occluded due to interaction of NgoL with other protein factors involved in MMR. The presentation of only one active site to substrate DNA will ensure that nicking of only one strand occurs to prevent inadvertent and deleterious double stranded cleavage.

A simplified method for the purification of riboflavin-binding protein from hen's egg yolk

A simple and efficient procedure for the purification of the riboflavin-binding protein from hen's egg yolk is described. This method involves the removal by exclusion of lipoproteins and subsequent fractionation of soluble yolk proteins held on a DEAE-cellulose column by a salt gradient which is followed by purification by gel filtration on Sephadex G-100. The protein thus isolated is homogeneous by various physicoehemical, immunological, and functional criteria.

Distinct Allostery Induced in the Cyclic GMP-binding, Cyclic GMP-specific Phosphodiesterase (PDE5) by Cyclic GMP, Sildenafil, and Metal Ions

The activity of many proteins orchestrating different biological processes is regulated by allostery, where ligand binding at one site alters the function of another site. Allosteric changes can be brought about by either a change in the dynamics of a protein, or alteration in its mean structure. We have investigated the mechanisms of allostery induced by chemically distinct ligands in the cGMP-binding, cGMP-specific phosphodiesterase, PDE5. PDE5 is the target for catalytic site inhibitors, such as sildenafil, that are used for the treatment of erectile dysfunction and pulmonary hypertension. PDE5 is a multidomain protein and contains two N-terminal cGMP-specific phosphodiesterase, bacterial adenylyl cyclase, FhLA transcriptional regulator (GAF) domains, and a C-terminal catalytic domain. Cyclic GMP binding to the GAFa domain and sildenafil binding to the catalytic domain result in conformational changes, which to date have been studied either with individual domains or with purified enzyme. Employing intramolecular bioluminescence resonance energy transfer, which can monitor conformational changes both in vitro and in intact cells, we show that binding of cGMP and sildenafil to PDE5 results in distinct conformations of the protein. Metal ions bound to the catalytic site also allosterically modulated cGMP- and sildenafil-induced conformational changes. The sildenafil-induced conformational change was temperature-sensitive, whereas cGMP-induced conformational change was independent of temperature. This indicates that different allosteric ligands can regulate the conformation of a multidomain protein by distinct mechanisms. Importantly, this novel PDE5 sensor has general physiological and clinical relevance because it allows the identification of regulators that can modulate PDE5 conformation in vivo.

Crystal structure of Salmonella typhimurium 2-methylcitrate synthase: Insights on domain movement and substrate specificity

2-Methylcitric acid (2-MCA) cycle is one of the well studied pathways for the utilization of propionate as a source of carbon and energy in bacteria such as Salmonella typhimurium and Escherichia coli. 2-Methylcitrate synthase (2-MCS) catalyzes the conversion of oxaloacetate and propionyl-CoA to 2-methylcitrate and CoA in the second step of 2-MCA cycle. Here, we report the X-ray crystal structure of S. typhimurium 2-MCS (StPrpC) at 2.4 A resolution and its functional characterization. StPrpC was found to utilize propionyl-CoA more efficiently than acetyl-CoA or butyryl-CoA. The polypeptide fold and the catalytic residues of StPrpC are conserved in citrate synthases (CSs) suggesting similarities in their functional mechanisms. In the triclinic P1 cell, StPrpC molecules were organized as decamers composed of five identical dimer units. In solution, StPrpC was in a dimeric form at low concentrations and was converted to larger oligomers at higher concentrations. CSs are usually dimeric proteins. In Gram-negative bacteria, a hexameric form, believed to be important for regulation of activity by NADH, is also observed. Structural comparisons with hexameric E. coil CS suggested that the key residues involved in NADH binding are not conserved in StPrpC. Structural comparison with the ligand free and bound states of CSs showed that StPrpC is in a nearly closed conformation despite the absence of bound ligands. It was found that the Tyr197 and Leu324 of StPrpC are structurally equivalent to the ligand binding residues His and Val, respectively, of CSs. These substitutions might determine the specificities for acyl-CoAs of these enzymes. (C) 2010 Elsevier Inc. All rights reserved.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

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.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

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.

Inhibition of Mycobacterium tuberculosis RNA Polymerase by Binding of a Gre Factor Homo log to the Secondary Channel

Because of its essential nature, each step of transcription, viz., initiation, elongation, and termination, is subjected to elaborate regulation. A number of transcription factors modulate the rates of transcription at these different steps, and several inhibitors shut down the process. Many modulators, including small molecules and proteinaceous inhibitors, bind the RNA polymerase (RNAP) secondary channel to control transcription. We describe here the first small protein inhibitor of transcription in Mycobacterium tuberculosis. Rv3788 is a homolog of the Gre factors that binds near the secondary channel of RNAP to inhibit transcription. The factor also affected the action of guanosine pentaphosphate (pppGpp) on transcription and abrogated Gre action, indicating its function in the modulation of the catalytic center of RNAP. Although it has a Gre factor-like domain organization with the conserved acidic residues in the N terminus and retains interaction with RNAP, the factor did not show any transcript cleavage stimulatory activity. Unlike Rv3788, another Gre homolog from Mycobacterium smegmatis, MSMEG_6292 did not exhibit transcription-inhibitory activities, hinting at the importance of the former in influencing the lifestyle of M. tuberculosis.

C2 domain is responsible for targeting rice phosphoinositide specific phospholipase C

Phosphoinositide-specific phospholipase C (PLC) is involved in Ca2+ mediated signalling events that lead to altered cellular status. Using various sequence-analysis methods, we identified two conserved motifs in known PLC sequences. The identified motifs are located in the C2 domain of plant PLCs and are not found in any other protein. These motifs are specifically found in the Ca2+ binding loops and form adjoining beta strands. Further, we identified certain conserved residues that are highly distinct from corresponding residues of animal PLCs. The motifs reported here could be used to annotate plant-specific phospholipase C sequences. Furthermore, we demonstrated that the C2 domain alone is capable of targeting PLC to the membrane in response to a Ca2+ signal. We also showed that the binding event results from a change in the hydrophobicity of the C2 domain upon Ca2+ binding. Bioinformatic analyses revealed that all PLCs from Arabidopsis and rice lack a transmembrane domain, myristoylation and GPI-anchor protein modifications. Our bioinformatic study indicates that plant PLCs are located in the cytoplasm, the nucleus and the mitochondria. Our results suggest that there are no distinct isoforms of plant PLCs, as have been proposed to exist in the soluble and membrane associated fractions. The same isoform could potentially be present in both subcellular fractions, depending on the calcium level of the cytosol. Overall, these data suggest that the C2 domain of PLC plays a vital role in calcium signalling.

The Hinge Region of Human Thyroid-Stimulating Hormone (TSH) Receptor Operates as a Tunable Switch between Hormone Binding and Receptor Activation

The mechanism by which the hinge regions of glycoprotein hormone receptors couple hormone binding to activation of downstream effecters is not clearly understood. In the present study, agonistic (311.62) and antagonistic (311.87) monoclonal antibodies (MAbs) directed against the TSH receptor extracellular domain were used to elucidate role of the hinge region in receptor activation. MAb 311.62 which identifies the LRR/Cb-2 junction (aa 265-275), increased the affinity of TSHR for the hormone while concomitantly decreasing its efficacy, whereas MAb 311.87 recognizing LRR 7-9 (aa 201-259) acted as a non-competitive inhibitor of Thyroid stimulating hormone (TSH) binding. Binding of MAbs was sensitive to the conformational changes caused by the activating and inactivating mutations and exhibited differential effects on hormone binding and response of these mutants. By studying the effects of these MAbs on truncation and chimeric mutants of thyroid stimulating hormone receptor (TSHR), this study confirms the tethered inverse agonistic role played by the hinge region and maps the interactions between TSHR hinge region and exoloops responsible for maintenance of the receptor in its basal state. Mechanistic studies on the antibody-receptor interactions suggest that MAb 311.87 is an allosteric insurmountable antagonist and inhibits initiation of the hormone induced conformational changes in the hinge region, whereas MAb 311.62 acts as a partial agonist that recognizes a conformational epitope critical for coupling of hormone binding to receptor activation. The hinge region, probably in close proximity with the alpha-subunit in the hormone-receptor complex, acts as a tunable switch between hormone binding and receptor activation.

The antibodies against the computationally designed mimic of the glycoprotein hormone receptor transmembrane domain provide insights into receptor activation and suppress the constitutively activated receptor Mutants

The exoloops of glycoprotein hormone receptors (GpHRs) transduce the signal generated by the ligand-ectodomain interactions to the transmembrane helices either through direct hormonal contact and/or by modulating the interdomain interactions between the hinge region (HinR) and the transmembrane domain (TMD). The ligand-induced conformational alterations in the HinRs and the interhelical loops of luteinizing hormone receptor/follicle stimulating hormone receptor/thyroid stimulating hormone receptor were mapped using exoloop-specific antibodies generated against a mini-TMD protein designed to mimic the native exoloop conformations that were created by joining the thyroid stimulating hormone receptor exoloops constrained through helical tethers and library-derived linkers. The antibody against the mini-TMD specifically recognized all three GpHRs and inhibited the basal and hormone-stimulated cAMP production without affecting hormone binding. Interestingly, binding of the antibody to all three receptors was abolished by prior incubation of the receptors with the respective hormones, suggesting that the exoloops are buried in the hormone-receptor complexes. The antibody also suppressed the high basal activities of gain-of-function mutations in the HinRs, exoloops, and TMDs such as those involved in precocious puberty and thyroid toxic adenomas. Using the antibody and point/deletion/chimeric receptor mutants, we demonstrate that changes in the HinR-exoloop interactions play an important role in receptor activation. Computational analysis suggests that the mini-TMD antibodies act by conformationally locking the transmembrane helices by means of restraining the exoloops and the juxta-membrane regions. Using GpHRs as a model, we describe a novel computational approach of generating soluble TMD mimics that can be used to explain the role of exoloops during receptor activation and their interplay with TMDs.

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background: Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism. Results: Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 angstrom resolution) and citrate-bound (Form-II, 1.90 angstrom resolution) forms. The enzyme showed broad substrate specificity with k(cat)/K-m in the order of acetate > propionate > formate. Further, the K-m for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5'-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic beta beta beta alpha beta alpha beta alpha topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230-300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes. Conclusions: The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.