Modeling the human intestinal Mucin (MUC2) C-terminal cystine knot dimer

Intestinal mucus, a viscous secretion that lines the mucosa, is believed to be a barrier to absorption of many therapeutic compounds and carriers, and is known to play an important physiological role in controlling pathogen invasion. Nevertheless, there is as yet no clear understanding of the barrier properties of mucus, such as the nature of the molecular interactions between drug molecules and mucus components as well as those that govern gel formation. Secretory mucins, large and complex glycoprotein molecules, are the principal determinants of the viscoelastic properties of intestinal mucus. Despite the important role that mucins play in controlling transport and in diseases such as cystic fibrosis, their structures remain poorly characterized. The major intestinal secretory mucin gene, MUC2, has been identified and fully sequenced. The present study was undertaken to determine a detailed structure of the cysteine-rich region within the C-terminal end of human intestinal mucin (MUC2) via homology modeling, and explore possible configurations of a dimer of this cysteine-rich region, which may play an important role in governing mucus gel formation. Based on sequence–structure alignments and three-dimensional modeling, a cystine knot tertiary structure homologous to that of human chorionic gonadotropin (HCG) is predicted at the C-terminus of MUC2. Dimers of this C-terminal cystine knot (CTCK) were modeled using sequence alignment based on HCG and TGF-beta, followed by molecular dynamics and simulated annealing. Results support the formation of a cystine knot dimer with a structure analogous to that of HCG.      

Comparative genomic analysis of the eight-membered ring cystine knot-containing bone morphogenetic protein antagonists

TGF-beta family proteins with a cystine knot motif serve as ligands for diverse families of plasma membrane receptors. Bone morphogenetic protein (BMP) antagonists represent a subgroup of these proteins, some of which bind BMPs and antagonize their actions during development and morphogenesis. Availability of completed genome sequences from diverse organisms allows bioinformatic analysis of the evolution of BMP antagonists and facilitates their classification. Using a regular expression algorithm (http://BioRegEx.stanford.edu), an exhaustive search of the human genome identified all cystine knot-containing BMP antagonists. Based on the size of the cystine ring, these proteins were divided into three subfamilies: CAN (eight-membered ring), twisted gastrulation (nine-membered ring), as well as chordin and noggin (10-membered ring). The CAN family can be divided further into four subgroups based on a conserved arrangement of additional cysteine residues-gremlin and PRDC, cerberus and coco, and DAN, together with USAG-1 and sclerostin. We searched for orthologs of human BMP antagonists in the genomes of model organisms and analyzed their phylogenetic relationship. New human paralogs were identified together with the verification of orthologous relationships of known genes. We also discuss the physiological roles of the CAN subfamily of BMP antagonists and the associated genetic defects. Based on the known three-dimensional structure of key cystine knot proteins, we postulated disulfide bondings for eight-membered ring BMP antagonists to predict their potential folding and dimerization.     

Aplysia temptin - the 'glue' in the water-borne attractin pheromone complex

Temptin, a component of the complex of water-borne protein pheromones that stimulate attraction and mating behavior in the marine mollusk Aplysia, has sequence homology to the epidermal growth factor (EGF)-like domains of higher organisms that mediate protein-cell surface contact during fertilization and blood coagulation. In this work, recombinant temptin for structural and functional studies was produced in Escherichia coli using a cold shock promoter and purified by RP-HPLC. CD spectra confirmed a predominantly beta-sheet structure. Two disulfide bonds were determined via limited proteolysis and MS. One internal disulfide (Cys57-Cys77) was predicted from initial alignments with class I EGF-like domains; the second, between Cys18 and Cys103, could protect temptin against proteolysis in seawater and stabilize its interacting surface. A three-dimensional model of temptin was prepared with our MPACK suite, based on the Ca(2+)-binding, EGF-like domain of the extracellular matrix protein fibrillin. Two temptin residues, Trp52 and Trp79, which align with cysteine residues conserved in fibrillins, lie adjacent to and could stabilize the disulfide bonds and a proposed metal-binding loop. The water-borne pheromone attractin in egg cordon eluates is complexed with other proteins. Docking results with our model and the NMR structure of attractin suggest that one face of temptin interacts with the pheromone, perhaps controlling its access to the cellular receptors. Gel shifts confirmed that temptin complexes with wild-type attractin. These results indicate that temptin, analogous to the role of fibrillin in controlling transforming growth factor-beta concentration, modulates pheromone signaling by direct binding to attractin.     

Bone morphogenetic proteins and their receptors in the eye

The human genome encodes at least 42 different members of the transforming growth factor-beta superfamily of growth factors. Bone morphogenetic proteins (BMPs) are the largest subfamily of proteins within the transforming growth factor-beta superfamily and are involved in numerous cellular functions including development, morphogenesis, cell proliferation, apoptosis, and extracellular matrix synthesis. This article first reviews BMPs and BMP receptors, BMP signaling pathways, and mechanisms controlling BMP signaling. Second, we review BMP and BMP receptor expression during embryonic ocular development/ differentiation and in adult ocular tissues. Lastly, future research directions with respect to BMP, BMP receptors, and ocular tissues are suggested.     

Identification of the gene encoding bursicon, an insect neuropeptide responsible for cuticle sclerotization and wing spreading

To accommodate growth, insects must periodically replace their exoskeletons. After shedding the old cuticle, the new soft cuticle must sclerotize. Sclerotization has long been known to be controlled by the neuropeptide hormone bursicon, but its large size of 30 kDa has frustrated attempts to determine its sequence and structure. Using partial sequences obtained from purified cockroach bursicon, we identified the Drosophila melanogaster gene CG13419 as a candidate bursicon gene. CG13419 encodes a peptide with a predicted final molecular weight of 15 kDa, which likely functions as a dimer. This predicted bursicon protein belongs to the cystine knot family, which includes vertebrate transforming growth factor-beta (TGF-beta) and glycoprotein hormones. Point mutations in the bursicon gene cause defects in cuticle sclerotization and wing expansion behavior. Bioassays show that these mutants have decreased bursicon bioactivity. In situ hybridization and immunocytochemistry revealed that bursicon is co-expressed with crustacean cardioactive peptide (CCAP). Transgenic flies that lack CCAP neurons also lacked bursicon bioactivity. Our results indicate that CG13419 encodes bursicon, the last of the classic set of insect developmental hormones. It is the first member of the cystine knot family to have a defined function in invertebrates. Mutants show that the spectrum of bursicon actions is broader than formerly demonstrated.     

Molecular architecture and biorecognition processes of the cystine knot protein superfamily: part I. The glycoprotein hormones

In this review article, the reader is introduced to recent advances in our knowledge on a subset of the cystine knot superfamily of homo- and hetero-dimeric proteins, from the perspective of the endocrine glycoprotein hormone family of proteins: follitropin (FSH), Iutropin (LH), thyrotropin. (TSH) and chorionic gonadotropin (CG). Subsequent papers will address the structure-function behaviour of other members of this increasingly significant family of proteins, including various members of the transforming growth factor-beta (TGF-beta) family of proteins, the activins, inhibins, bone morphogenic growth factor, platelet derived growth factor-beta, nerve growth factor and more than 35 other proteins with similar topological features. In the present review article, specific emphasis has been placed on advances with the glycoprotein hormones (GPHs) that have facilitated greater insight into their physiological functions, molecular structures and most importantly the basis of the molecular recognition events that lead to the formation of hetero-dimeric structures as well as their specific and selective recognition by their corresponding receptors and antibodies. Thus, this review article focuses on the structural motifs involved in receptor recognition and the current techniques available to identify these regions, including the role of immunological methodology, peptide fragment design and synthesis and mutagenesis to delineate their structure-function relationships and molecular recognition behaviour.     

Surface proteome analysis identifies platelet derived growth factor receptor-alpha as a critical mediator of transforming growth factor-beta-induced collagen secretion

Fibroblasts are extracellular matrix-producing cells in the lung. Fibroblast activation by transforming growth factor-beta leads to myofibroblast-differentiation and increased extracellular matrix deposition, a hallmark of pulmonary fibrosis. While fibroblast function with respect to migration, invasion, and extracellular matrix deposition has been well-explored, little is known about the surface proteome of lung fibroblasts in general and its specific response to fibrogenic growth factors, in particular transforming growth factor-beta.We thus performed a cell-surface proteome analysis of primary human lung fibroblasts in presence/absence of transforming growth factor-beta, followed by characterization of our findings using FACS analysis, Western blot, and siRNA-mediated knockdown experiments.We identified 213 surface proteins significantly regulated by transforming growth factor-beta, platelet derived growth factor receptor-alpha being one of the top down-regulated proteins. Transforming growth factor beta-induced downregulation of platelet derived growth factor receptor-alpha induced upregulation of platelet derived growth factor receptor-beta expression and phosphorylation of Akt, a downstream target of platelet derived growth factor signaling. Importantly, collagen type V expression and secretion was strongly increased after forced knockdown of platelet derived growth factor receptor-alpha, an effect that was potentiated by transforming growth factor-beta. We therefore show previously underappreciated cross-talk of transforming growth factor-beta and platelet derived growth factor signaling in human lung fibroblasts, resulting in increased extracellular matrix deposition in a platelet derived growth factor receptor-alpha dependent manner. These findings are of particular importance for the treatment of lung fibrosis patients with high pulmonary transforming growth factor-beta activity.     

Iterative sequence/secondary structure search for protein homologs: comparison with amino acid sequence alignments and application to fold recognition in genome databases

MOTIVATION: Sequence alignment techniques have been developed into extremely powerful tools for identifying the folding families and function of proteins in newly sequenced genomes. For a sufficiently low sequence identity it is necessary to incorporate additional structural information to positively detect homologous proteins. We have carried out an extensive analysis of the effectiveness of incorporating secondary structure information directly into the alignments for fold recognition and identification of distant protein homologs. A secondary structure similarity matrix based on a database of three-dimensionally aligned proteins was first constructed. An iterative application of dynamic programming was used which incorporates linear combinations of amino acid and secondary structure sequence similarity scores. Initially, only primary sequence information is used. Subsequently contributions from secondary structure are phased in and new homologous proteins are positively identified if their scores are consistent with the predetermined error rate. RESULTS: We used the SCOP40 database, where only PDB sequences that have 40% homology or less are included, to calibrate homology detection by the combined amino acid and secondary structure sequence alignments. Combining predicted secondary structure with sequence information results in a 8-15% increase in homology detection within SCOP40 relative to the pairwise alignments using only amino acid sequence data at an error rate of 0.01 errors per query; a 35% increase is observed when the actual secondary structure sequences are used. Incorporating predicted secondary structure information in the analysis of six small genomes yields an improvement in the homology detection of approximately 20% over SSEARCH pairwise alignments, but no improvement in the total number of homologs detected over PSI-BLAST, at an error rate of 0.01 errors per query. However, because the pairwise alignments based on combinations of amino acid and secondary structure similarity are different from those produced by PSI-BLAST and the error rates can be calibrated, it is possible to combine the results of both searches. An additional 25% relative improvement in the number of genes identified at an error rate of 0.01 is observed when the data is pooled in this way. Similarly for the SCOP40 dataset, PSI-BLAST detected 15% of all possible homologs, whereas the pooled results increased the total number of homologs detected to 19%. These results are compared with recent reports of homology detection using sequence profiling methods. AVAILABILITY: Secondary structure alignment homepage at http://lutece.rutgers.edu/ssas CONTACT: anders@rutchem.rutgers.edu; ronlevy@lutece.rutgers.edu Supplementary Information: Genome sequence/structure alignment results at http://lutece.rutgers.edu/ss_fold_predictions.     

The Arabidopsis ESCRT protein�Cprotein interaction network

In yeast, endosomal sorting of monoubiquitylated transmembrane proteins is performed by a subset of the 19 "class E vacuolar protein sorting" proteins. The core machinery consists of 11 proteins that are organised in three complexes termed ESCRT I-III (endosomal sorting complex required for transport I-III) and is conserved in eukaryotic cells. While the pathway is well understood in yeast and animals, the plant ESCRT system is largely unexplored. At least one sequence homolog for each ESCRT component can be found in the Arabidopsis genome. Generally, sequence conservation between yeast/animals and the Arabidopsis proteins is low. To understand details about participating proteins and complex organization we have performed a systematic pairwise yeast two hybrid analysis of all Arabidopsis proteins showing homology to the ESCRT core machinery. Positive interactions were validated using bimolecular fluorescence complementation. In our experiments, most putative ESCRT components exhibited interactions with other ESCRT components that could be shown to occur on endosomes suggesting that despite their low homology to their yeast and animal counterparts they represent functional components of the plant ESCRT pathway.     

Knots in rings. The circular knotted protein Momordica cochinchinensis trypsin inhibitor-II folds via a stable two-disulfide intermediate

The aim of this work was to elucidate the oxidative folding mechanism of the macrocyclic cystine knot protein MCoTI-II. We aimed to investigate how the six-cysteine residues distributed on the circular backbone of the reduced unfolded peptide recognize their correct partner and join up to form a complex cystine-knotted topology. To answer this question, we studied the oxidative folding of the naturally occurring peptide using a range of spectroscopic methods. For both oxidative folding and reductive unfolding, the same disulfide intermediate species was prevalent and was characterized to be a native-like two-disulfide intermediate in which the Cys1-Cys18 disulfide bond was absent. Overall, the folding pathway of this head-to-tail cyclized protein was found to be similar to that of linear cystine knot proteins from the squash family of trypsin inhibitors. However, the pathway differs in an important way from that of the cyclotide kalata B1, in that the equivalent two-disulfide intermediate in that case is not a direct precursor of the native protein. The size of the embedded ring within the cystine knot motif appears to play a crucial role in the folding pathway. Larger rings contribute to the independence of disulfides and favor an on-pathway native-like intermediate that has a smaller energy barrier to cross to form the native fold. The fact that macrocyclic proteins are readily able to fold to a complex knotted structure in vitro in the absence of chaperones makes them suitable as protein engineering scaffolds that have remarkable stability.     

Analysis on sliding helices and strands in protein structural comparisons: a case study with protein kinases

Protein structural alignments are generally considered as 'golden standard' for the alignment at the level of amino acid residues. In this study we have compared the quality of pairwise and multiple structural alignments of about 5900 homologous proteins from 718 families of known 3-D structures. We observe shifts in the alignment of regular secondary structural elements (helices and strands) between pairwise and multiple structural alignments. The differences between pairwise and multiple structural alignments within helical and beta-strand regions often correspond to 4 and 2 residue positions respectively. Such shifts correspond approximately to "one turn" of these regular secondary structures. We have performed manual analysis explicitly on the family of protein kinases. We note shifts of one or two turns in helix-helix alignments obtained using pairwise and multiple structural alignments. Investigations on the quality of the equivalent helix-helix, strand-strand pairs in terms of their residue side-chain accessibilities have been made. Our results indicate that the quality of the pairwise alignments is comparable to that of the multiple structural alignments and, in fact, is often better. We propose that pairwise alignment of protein structures should also be used in formulation of methods for structure prediction and evolutionary analysis.     

Accuracy of structure-based sequence alignment of automatic methods

Background  

Accurate sequence alignments are essential for homology searches and for building three-dimensional structural models of proteins. Since structure is better conserved than sequence, structure alignments have been used to guide sequence alignments and are commonly used as the gold standard for sequence alignment evaluation. Nonetheless, as far as we know, there is no report of a systematic evaluation of pairwise structure alignment programs in terms of the sequence alignment accuracy.     

Structure alignment based on coding of local geometric measures

Background  

A structure alignment method based on a local geometric property is presented and its performance is tested in pairwise and multiple structure alignments. In this approach, the writhing number, a quantity originating from integral formulas of Vassiliev knot invariants, is used as a local geometric measure. This measure is used in a sliding window to calculate the local writhe down the length of the protein chain. By encoding the distribution of writhing numbers across all the structures in the protein databank (PDB), protein geometries are represented in a 20-letter alphabet. This encoding transforms the structure alignment problem into a sequence alignment problem and allows the well-established algorithms of sequence alignment to be employed. Such geometric alignments offer distinct advantages over structural alignments in Cartesian coordinates as it better handles structural subtleties associated with slight twists and bends that distort one structure relative to another.     

Sorting nexin 6, a novel SNX, interacts with the transforming growth factor-beta family of receptor serine-threonine kinases

Sorting nexins (SNX) comprise a family of proteins with homology to several yeast proteins, including Vps5p and Mvp1p, that are required for the sorting of proteins to the yeast vacuole. Human SNX1, -2, and -4 have been proposed to play a role in receptor trafficking and have been shown to bind to several receptor tyrosine kinases, including receptors for epidermal growth factor, platelet-derived growth factor, and insulin as well as the long form of the leptin receptor, a glycoprotein 130-associated receptor. We now describe a novel member of this family, SNX6, which interacts with members of the transforming growth factor-beta family of receptor serine-threonine kinases. These receptors belong to two classes: type II receptors that bind ligand, and type I receptors that are subsequently recruited to transduce the signal. Of the type II receptors, SNX6 was found to interact strongly with ActRIIB and more moderately with wild type and kinase-defective mutants of TbetaRII. Of the type I receptors, SNX6 was found to interact only with inactivated TbetaRI. SNXs 1-4 also interacted with the transforming growth factor-beta receptor family, showing different receptor preferences. Conversely, SNX6 behaved similarly to the other SNX proteins in its interactions with receptor tyrosine kinases. Strong heteromeric interactions were also seen among SNX1, -2, -4, and -6, suggesting the formation in vivo of oligomeric complexes. These findings are the first evidence for the association of the SNX family of molecules with receptor serine-threonine kinases.     

Differential in vitro phenotype pattern, transforming growth factor-β1 activity and mRNA expression of transforming growth factor-β1 in Apert osteoblasts

The phenotype of Apert osteoblasts differs from that of normal osteoblasts in the accumulation of macromolecules in the extracellular matrix. Apert osteoblasts increase type I collagen, fibronectin and glycosaminoglycans secretion compared with normal osteoblasts. Because the extracellular matrix macromolecule accumulation is greatly modulated by transforming growth factor-beta(1), we examined the ability of normal and Apert osteoblasts to secrete transforming growth factor-beta(1) by CCL-64 assay and to produce transforming growth factor-beta(1 )by analysis of the mRNA expression of transforming growth factor-beta(1). Northern blot analysis revealed an increased amount of transforming growth factor-beta(1) mRNA expression in Apert osteoblasts compared with normal ones. Moreover, the level of the active transforming growth factor-beta(1) isoform was higher in Apert than in normal media. In pathologic cells, the increase in transforming growth factor-beta(1) gene expression was associated with a parallel increase in the factor secreted into the medium. The level of transforming growth factor-beta(1) was decreased by the addition of basic fibroblast growth factor. Transforming growth factor-beta(1) is controlled temporally and spatially during skeletal tissue development and produces complex stimulatory and inhibitory changes in osteoblast functions. We hypothesise that in vitro differences between normal and Apert osteoblasts may be correlated to different transforming growth factor-beta(1) cascade patterns, probably due to an altered balance between transforming growth factor-beta(1) and basic fibroblast growth factor.