Evolutionary trace analysis of TGF-beta and related growth factors: implications for site-directed mutagenesis

Identification and functional characterization of multiple inositol polyphosphate phosphatase1 (Minpp1) isoform-2 in exosomes with potential to modulate tumor microenvironment

Inositol polyphosphates (InsPs) play key signaling roles in diverse cellular functions, including calcium homeostasis, cell survival and death. Multiple inositol polyphosphate phosphatase 1 (Minpp1) affects the cellular levels of InsPs and cell functions. The Minpp1 is an endoplasmic reticulum (ER) resident but localizes away from its cytosolic InsPs substrates. The current study examines the heterogeneity of Minpp1 and the potential physiologic impact of Minpp1 isoforms, distinct motifs, subcellular distribution, and enzymatic potential. The NCBI database was used to analyze the proteome diversity of Minpp1 using bioinformatics tools. The analysis revealed that translation of three different Minpp1 variants resulted in three isoforms of Minpp1 of varying molecular weights. A link between the minpp1 variant-2 gene and ER-stress, using real-time PCR, suggests a functional similarity between minpp1 variant-1 and variant-2. A detailed study on motifs revealed Minpp1 isoform-2 is the only other isoform, besides isoform-1, that carries a phosphatase motif for InsPs hydrolysis but no ER-retention signal. The confocal microscopy revealed that the Minpp1 isoform-1 predominantly localized near the nucleus with a GRP-78 ER marker, while Minpp1 isoform-2 was scattered more towards the cell periphery where it co-localizes with the plasma membrane-destined multivesicular bodies biomarker CD63. MCF-7 cells were used to establish that Minpp1 isoform-2 is secreted into exosomes. Brefeldin A treatment resulted in overexpression of the exosome-associated Minpp1 isoform-2, suggesting its secretion via an unconventional route involving endocytic-generated vesicles and a link to ER stress. Results further demonstrated that the exosome-associated Minpp1 isoform-2 was enzymatically active. Overall, the data support the possibility that an extracellular form of enzymatically active Minpp1 isoform-2 mitigates any anti-proliferative actions of extracellular InsPs, thereby also impacting the makeup of the tumor microenvironment.

In Silico Characterization of Growth Differentiation Factors as Inhibitors of TNF-Alpha and IL-6 in Immune-Mediated Inflammatory Disease Rheumatoid Arthritis

Tumor necrosis factor alpha (TNF-α) plays a critical role in the progression of inflammation and affects the cells of the synovial membrane. Another key factor in the progression of rheumatoid inflammation is interleukin-6 (IL-6). Both TNF-α and IL-6 promote the proliferation of synovial membrane cells thus stimulating the production of matrix metalloproteinases and other cytotoxins and leading towards bone erosion and destruction of the cartilage. Growth differentiation factor-11 (GDF11) and growth differentiation factor-8 (GDF8) which is also known as myostatin are members of the transforming growth factor-β family and could be used as antagonists to inflammatory responses which are associated with rheumatoid arthritis. In the current study, to elucidate the evolutionary relationships of GDF11 with its homologs from other closely related organisms, a comprehensive phylogenetic analysis was performed. From the phylogram, it was revealed that the clade of Primates that belong to superorder Euarchontoglires showed close evolutionary relationships with order Cetartiodactyla of the Laurasiatheria superorder. Fifty tetrapeptides were devised from conserved regions of GDF11 which served as ligands in protein-ligand docking against TNF-α and IL-6 followed by drug scanning and ADMET profiling of best selected ligands. The peptides SAGP showed strong interactions with IL-6, and peptides AFDP and AGPC showed strong interactions with TNF-α, and all three peptides fulfilled all the pharmacokinetic parameters which are important for bioavailability. The potential of GDF8 as an antagonist to TNF-α and IL-6 was also explored using a protein-protein docking approach. The binding patterns of GDF8 with TNF-α and IL-6 showed that GDF8 could be used as a potential inhibitor of TNF-α and IL-6 to treat rheumatoid arthritis.

Alterations in activin A-myostatin-follistatin system associate with disease activity in inflammatory myopathies

OBJECTIVES

The aim of this study was to investigate the systemic and skeletal muscle levels of atrophy-associated myokines in patients with idiopathic inflammatory myopathies (IIM) and their association with clinical characteristics of myositis.

METHODS

A total of 94 IIM patients and 162 healthy controls were recruited. Of those, 20 IIM patients and 28 healthy controls underwent a muscle biopsy. Circulating concentrations of myostatin, follistatin, activin A and TGF-β1 were assessed by ELISA. The expression of myokines and associated genes involved in the myostatin signalling pathway in muscle tissue was determined by real-time PCR.

RESULTS

We report decreased levels of circulating myostatin (median 1817 vs 2659 pg/ml; P = 0.003) and increased follistatin (1319 vs 1055 pg/ml; P = 0.028) in IIM compared with healthy controls. Activin A levels were also higher in IIM (414 vs 309 pg/ml; P = 0.0005) compared with controls. Myostatin was negatively correlated to muscle disease activity assessed by physician on visual analogue scale (MDA) (r = -0.289, P = 0.015) and positively to manual muscle testing of eight muscles (r = 0.366, P = 0.002). On the other hand, follistatin correlated positively with MDA (r = 0.235, P = 0.047). Gene expression analysis showed higher follistatin (P = 0.003) and myostatin inhibitor follistatin-like 3 protein (FSTL3) (P = 0.008) and lower expression of activin receptor type 1B (ALK4) (P = 0.034), signal transducer SMAD3 (P = 0.023) and atrophy marker atrogin-1 (P = 0.0009) in IIM muscle tissue compared with controls.

CONCLUSION

This study shows lower myostatin and higher follistatin levels in circulation and attenuated expression of myostatin pathway signalling components in skeletal muscle of patients with myositis, a newly emerging pattern of the activin A-myostatin-follistatin system in muscle wasting diseases.

A novel strategy for molecular interfaces optimization: The case of Ferritin-Transferrin receptor interaction

Protein-protein interactions regulate almost all cellular functions and rely on a fine tune of surface amino acids properties involved on both molecular partners. The disruption of a molecular association can be caused even by a single residue mutation, often leading to a pathological modification of a biochemical pathway. Therefore the evaluation of the effects of amino acid substitutions on binding, and the ad hoc design of protein-protein interfaces, is one of the biggest challenges in computational biology. Here, we present a novel strategy for computational mutation and optimization of protein-protein interfaces. Modeling the interaction surface properties using the Zernike polynomials, we describe the shape and electrostatics of binding sites with an ordered set of descriptors, making possible the evaluation of complementarity between interacting surfaces. With a Monte Carlo approach, we obtain protein mutants with controlled molecular complementarities. Applying this strategy to the relevant case of the interaction between Ferritin and Transferrin Receptor, we obtain a set of Ferritin mutants with increased or decreased complementarity. The extensive molecular dynamics validation of the method results confirms its efficacy, showing that this strategy represents a very promising approach in designing correct molecular interfaces.

Antibiotic Substrate Selectivity of Pseudomonas aeruginosa MexY and MexB Efflux Systems Is Determined by a Goldilocks Affinity

Resistance-nodulation-division (RND) efflux pumps are important contributors to bacterial antibiotic resistance. In this study, we combined evolutionary sequence analyses, computational structural modeling, and ligand docking to develop a framework that can explain the known antibiotic substrate selectivity differences between two Pseudomonas aeruginosa RND transporters, MexY and MexB. For efficient efflux, antibiotic substrates must possess a "Goldilocks affinity": binding strong enough to allow interaction with transporter but not so tight as to impede movement through the pump.

A Group of Olfactory Receptor Alleles that Encode Full Length Proteins are Down-Regulated as Olfactory Sensory Neurons Mature

The family of olfactory receptors (ORs) subserves the sense of smell and includes both functional alleles and pseudogenes, the latter identified by mutations resulting in frame shift or premature truncation. During neuronal differentiation, nonfunctional ORs are expressed initially but then are switched out, and/or the olfactory sensory neurons (OSNs) expressing them die. We carried out a transcriptomic analysis of FACS-isolated cells from ΔSox2-eGFP, Neurog1-eGFP BAC and ΔOMP-eGFP strains of uninjured and olfactory bulbectomized transgenic mice that correspond to distinct stages in the progression from globose basal cell stem cells to fully mature OSNs. We analyzed the expression pattern of 1094 unique receptors across this progression and found that the vast majority were characterized by a typical and expected pattern of expression; i.e., levels of OR mRNA peaking in mature OSNs. However, 43 ORs, including several known pseudogenes, were different, such that mRNA expression declined in the mature OSNs relative to earlier stages. Protein and promoter sequence analysis of the atypical group did not uncover any obvious differences between them and more typical ORs. Nonetheless, the pattern of expression suggests that atypical ORs may be non-functional despite the lack of any obvious abnormality in the sequence analyses.

In silico analysis of class I adenylate-forming enzymes reveals family and group-specific conservations

Luciferases, aryl- and fatty-acyl CoA synthetases, and non-ribosomal peptide synthetase proteins belong to the class I adenylate-forming enzyme superfamily. The reaction catalyzed by the adenylate-forming enzymes is categorized by a two-step process of adenylation and thioesterification. Although all of these proteins perform a similar two-step process, each family may perform the process to yield completely different results. For example, luciferase proteins perform adenylation and oxidation to produce the green fluorescent light found in fireflies, while fatty-acyl CoA synthetases perform adenylation and thioesterification with coenzyme A to assist in metabolic processes involving fatty acids. This study aligned a total of 374 sequences belonging to the adenylate-forming superfamily. Analysis of the sequences revealed five fully conserved residues throughout all sequences, as well as 78 more residues conserved in at least 60% of sequences aligned. Conserved positions are involved in magnesium and AMP binding and maintaining enzyme structure. Also, ten conserved sequence motifs that included most of the conserved residues were identified. A phylogenetic tree was used to assign sequences into nine different groups. Finally, group entropy analysis identified novel conservations unique to each enzyme group. Common group-specific positions identified in multiple groups include positions critical to coordinating AMP and the CoA-bound product, a position that governs active site shape, and positions that help to maintain enzyme structure through hydrogen bonds and hydrophobic interactions. These positions could serve as excellent targets for future research.

Extended Impact of Pin1 Catalytic Loop Phosphorylation Revealed by S71E Phosphomimetic

Pin1 is a two-domain human protein that catalyzes the cis-trans isomerization of phospho-Ser/Thr-Pro (pS/T-P) motifs in numerous cell-cycle regulatory proteins. These pS/T-P motifs bind to Pin1's peptidyl-prolyl isomerase (PPIase) domain in a catalytic pocket, between an extended catalytic loop and the PPIase domain core. Previous studies showed that post-translational phosphorylation of S71 in the catalytic loop decreases substrate binding affinity and isomerase activity. To define the origins for these effects, we investigated a phosphomimetic Pin1 mutant, S71E-Pin1, using solution NMR. We find that S71E perturbs not only its host loop but also the nearby PPIase core. The perturbations identify a local network of hydrogen bonds and salt bridges that is more extended than previously thought, and includes interactions between the catalytic loop and the α2/α3 turn in the PPIase core. Explicit-solvent molecular dynamics simulations and phylogenetic analysis suggest that these interactions act as conserved "latches" between the loop and PPIase core that enhance binding of phosphorylated substrates, as they are absent in PPIases lacking pS/T-P specificity. Our results suggest that S71 is a hub residue within an electrostatic network primed for phosphorylation, and may illustrate a common mechanism of phosphorylation-mediated allostery.

Structural complexity and functional diversity of plant NADPH oxidases

Plant NADPH oxidases also known as respiratory burst oxidase homologs (Rbohs) are a family of membrane-bound enzymes that play diverse roles in the defense response and morphogenetic processes via regulated generation of reactive oxygen species. Rbohs are associated with a variety of functions, although the reason for this is not clear. To evaluate using bioinformatics, the possible mechanisms for the observed functional diversity within the plant kingdom, 127 Rboh protein sequences representing 26 plant species were analyzed. Multiple clusters were identified with gene duplications that were both dicot as well as monocot-specific. The N-terminal sequences were observed to be highly variable. The conserved cysteine (equivalent of Cys890) in C-terminal of AtRbohD suggested that the redox-based modification like S-nitrosylation may regulate the activity of other Rbohs. Three-dimensional models corresponding to the N-terminal domain for Rbohs from Arabidopsis thaliana and Oryza sativa were constructed and molecular dynamics studies were carried out to study the role of Ca2+ in the folding of Rboh proteins. Certain mutations indicated possibly affect the structure and function of the plant NADPH oxidases, thereby providing the rationale for further experimental validation.

The Rv3799-Rv3807 Gene Cluster in Mycobacterium Tuberculosis Genome Corresponds to the ‘Ancient Conserved Region’ in CMN Mycolyltransferases

We have identified based on gene cluster analysis that the genes between Rv3799–Rv3807 in M. tuberculosis have orthologs in Corynebacteria, Mycobacteria and Nocardia (CMN) genomes. Therefore, this gene cluster possibly corresponds to the ‘Ancient Conserved Region’ of CMN mycolyltransferases. The evolutionary trace analysis suggests that twelve amino acid residues; Leu39, Trp51, Pro71, Trp82, Trp97, Phe100, Gly124, Ser126, Asp192, Glu230, Gly260 and Trp264 are ‘absolutely conserved’. These amino acid residues constitute the active site and conserved hydrophobic tunnel in CMN mycolyltransferases.

In silico analysis of heme oxygenase structural homologues identifies group-specific conservations

Heme oxygenases (HO) catalyze the breakdown of heme, aiding the recycling of its components. Several other enzymes have homologous tertiary structures to HOs, while sharing little sequence homology. These homologues include thiaminases, the hydroxylase component of methane monooxygenases, and the R2 component of Class I ribonucleotide reductases (RNR). This study compared these structural homologues of HO, using a large number of protein sequences for each homologue. Alignment of a total of 472 sequences showed little sequence conservation, with no residues having conservation in more than 80% of aligned sequences and only five residues conserved in at least 60% of the sequences. Fourteen additional positions, most of which were critical for hydrophobic packing, displayed amino acid similarity of 60% or higher. Ten conserved sequence motifs were identified in HOs and RNRs. Phylogenetic analysis verified the existence of the four distinct groups of HO homologues, which were then analyzed by group entropy analysis to identify residues critical to the unique function of each enzyme. Other methods for determining functional residues were also performed. Several common index positions identified represent critical evolutionary changes that resulted in the unique function of each enzyme, suggesting potential targets for site‐directed mutagenesis. These positions included residues that coordinate ligands, form the active sites, and maintain enzyme structure.

Enzymes

Heme oxygenase (EC 1.14.14.18), methane monooxygenase (EC 1.14.13.25), ribonucleotide reductase (EC 1.17.4.1), thiaminase II (EC 3.5.99.2).

Transmission analysis of TGFB1 gene polymorphisms in non-syndromic cleft lip with or without cleft palate

OBJECTIVES

Transforming growth factor beta1 (TGF-β1) plays a significant role in craniofacial development. Previous linkage studies reported that the TGF-β1-locus at 19q13.1 harbour predisposing genes for non-syndromic oral clefts. In the present study case parents triads were evaluated to find the transmission effects of genetic variants in TGF- β1 towards non-syndromic cleft lip or palate (NSCL/P).

METHODS

Using allelic discrimination method148 families (case-parent triads) were assessed for single nucleotide polymorphisms (SNPs) in TGF-β1 gene. The SNPs were checked for mendelian errors and Hardy-Weinberg equilibrium (HWE). Transmission disequilibrium test and haplotype frequencies were estimated.

RESULTS

The TGF-β1 SNPs showed very low minor allele frequencies (MAFs) and observed heterozygosity (Hobs). The transmission disequilibrium test (TDT) and parent-of-origin likelihood ratio tests (PO-LRT) were not significant for any of the SNPs tested. Strong linkage disequilibrium (r2 = 0.722) was found between rs1800469 and rs1800470 SNPs. Haplotype analysis ignoring parent of origin showed strong evidence of excess transmission but it is not significant (p-value = 0.293).

CONCLUSION

Transmission of minor alleles were not observed from either parent indicating that the TGF-β1 gene polymorphisms by themselves do not confer risk for non-syndromic oral clefts but, rather, modify the stability and the activation process of TGF-β1. As the number of families included in the study are less, results must be considered still preliminary and require replication using more families.

Structural and evolutionary analyses reveal determinants of DNA binding specificities of nucleoid-associated proteins HU and IHF

Nucleoid-associated proteins (NAPs) are chromosome-organizing factors, which affect the transcriptional landscape of a bacterial cell. HU is an NAP, which binds to DNA with a broad specificity while homologous IHF (Integration Host Factor), binds DNA with moderately higher specificity. Specificity and differential binding affinity of HU/IHF proteins towards their target binding sites play a crucial role in their regulatory dynamics. Decades of biochemical and genomic studies have been carried out for HU and IHF like proteins. Yet, questions related to their DNA binding specificity, and differential ability to bend DNA thus affecting the binding site length remained unanswered. In addition, the problem has not been investigated from an evolutionary perspective. Our phylogenetic analysis revealed three major clades belonging to HU, IHFα and IHFβ like proteins with reference to E. coli. We carried out a comparative analysis of three-dimensional structures of HU/IHF proteins to gain insight into the structural basis of clade division. The present study revealed three major features which contribute to differential DNA binding specificity of HU/IHF proteins, (I) conformational restriction of DNA binding residues due to salt-bridge formation, (II) the enrichment of alanine in the DNA binding site increasing conformational space of flexible side chains in its vicinity and (III) nature of DNA binding residue (Arg to Lys bias in different clades) which interacts differentially to DNA bases. We observed an extended electropositive surface at the DNA draping site for IHF clade proteins compared to HU, which stabilizes the DNA bend. Differences in the dimer stabilization strategies between HU and IHF were also observed. Our analysis reveals a comprehensive evolutionary picture, which rationalizes the origin of multi-specificity of HU/IHF proteins using sequence and structure-based determinants, which could also be applied to understand differences in binding specificities of other nucleic acid binding proteins.

Evolution of the Twist Subfamily Vertebrate Proteins: Discovery of a Signature Motif and Origin of the Twist1 Glycine-Rich Motifs in the Amino-Terminus Disordered Domain

Twist proteins belong to the basic helix-loop-helix (bHLH) family of multifunctional transcriptional factors. These factors are known to use domains other than the common bHLH in protein-protein interactions. There has been much work characterizing the bHLH domain and the C-terminus in protein-protein interactions but despite a few attempts more focus is needed at the N-terminus. Since the region of highest diversity in Twist proteins is the N-terminus, we analyzed the conservation of this region in different vertebrate Twist proteins and study the sequence differences between Twist1 and Twist2 with emphasis on the glycine-rich regions found in Twist1. We found a highly conserved sequence motif in all Twist1 (SSSPVSPADDSLSNSEEE) and Twist2 (SSSPVSPVDSLGTSEEE) mammalian species with unknown function. Through sequence comparison we demonstrate that the Twist protein family ancestor was “Twist2-like” and the two glycine-rich regions found in Twist1 sequences were acquired late in evolution, apparently not at the same time. The second glycine-rich region started developing first in the fish vertebrate group, while the first glycine region arose afterwards within the reptiles. Disordered domain and secondary structure predictions showed that the amino acid sequence and disorder feature found at the N-terminus is highly evolutionary conserved and could be a functional site that interacts with other proteins. Detailed examination of the glycine-rich regions in the N-terminus of Twist1 demonstrate that the first region is completely aliphatic while the second region contains some polar residues that could be subject to post-translational modification. Phylogenetic and sequence space analysis showed that the Twist1 subfamily is the result of a gene duplication during Twist2 vertebrate fish evolution, and has undergone more evolutionary drift than Twist2. We identified a new signature motif that is characteristic of each Twist paralog and identified important residues within this motif that can be used to distinguish between these two paralogs, which will help reduce Twist1 and Twist2 sequence annotation errors in public databases.

Computational analysis on conformational dynamics of bone morphogenetic protein-2 (BMP-2)

BMP-2 is widely used for bone regeneration because of its ability to induce osteoblast differentiation and proliferation. The pharmaceutical application of BMP-2 as bone implant makes the studies on stability and conformational dynamics very relevant as proteins are functional only in their native three-dimensional state. Knowing the factors affecting BMP-2 structure becomes essential for designing bone implants activated by BMP-2. In order to explore the influence of temperature and hydration on protein conformation, we have performed the molecular dynamics (MD) simulations at the time scale of 100 ns with two different force fields. We have examined the dynamic behaviour of BMP-2 monomer and dimer in aqueous medium as well as in vacuum at four different temperatures (300, 350, 400 and 450 K). MD simulation of BMP-2 monomer and dimer in water and vacuum environments shows the major contribution of water in structure stabilization. Temperature of the system affects the secondary structure differently in case of monomer and dimer simulation and the dynamics also depends on the environment viz. vacuum and aqueous. Vacuum simulations show very early loss of the major secondary structure content. On the other hand, BMP-2 monomer and dimer in aqueous environment show the unfolding of α-helix with increasing temperature. This unfolded α-helix is converted into β-sheet at 400 K in monomer of BMP-2. Contrary to this, we did not observe β-sheet formation in dimer BMP-2 even at 450 K indicating that monomers are more aggregation prone entity as compared to dimers of BMP-2.

Animal Mitochondrial DNA Replication

Recent advances in the field of mitochondrial DNA (mtDNA) replication highlight the diversity of both the mechanisms utilized and the structural and functional organization of the proteins at mtDNA replication fork, despite the relative simplicity of the animal mtDNA genome. DNA polymerase γ, mtDNA helicase and mitochondrial single-stranded DNA-binding protein-the key replisome proteins, have evolved distinct structural features and biochemical properties. These appear to be correlated with mtDNA genomic features in different metazoan taxa and with their modes of DNA replication, although substantial integrative research is warranted to establish firmly these links. To date, several modes of mtDNA replication have been described for animals: rolling circle, theta, strand-displacement, and RITOLS/bootlace. Resolution of a continuing controversy relevant to mtDNA replication in mammals/vertebrates will have a direct impact on the mechanistic interpretation of mtDNA-related human diseases. Here we review these subjects, integrating earlier and recent data to provide a perspective on the major challenges for future research.

Crystal structure of human GDF11

Members of the TGF-β family of proteins are believed to play critical roles in cellular signaling processes such as those involved in muscle differentiation. The extent to which individual family members have been characterized and linked to biological function varies greatly. The role of myostatin, also known as growth differentiation factor 8 (GDF8), as an inhibitor of muscle differentiation is well understood through genetic linkages. In contrast, the role of growth differentiation factor 11 (GDF11) is much less well understood. In humans, the mature forms of GDF11 and myostatin are over 94% identical. In order to understand the role that the small differences in sequence may play in the differential signaling of these molecules, the crystal structure of GDF11 was determined to a resolution of 1.50 Å. A comparison of the GDF11 structure with those of other family members reveals that the canonical TGF-β domain fold is conserved. A detailed structural comparison of GDF11 and myostatin shows that several of the differences between these proteins are likely to be localized at interfaces that are critical for the interaction with downstream receptors and inhibitors.

Antiplasmodial Activity Is an Ancient and Conserved Feature of Tick Defensins

Ancestral sequence reconstruction has been widely used to test evolution-based hypotheses. The genome of the European tick vector, Ixodes ricinus, encodes for defensin peptides with diverse antimicrobial activities against distantly related pathogens. These pathogens include fungi, Gram-negative, and Gram-positive bacteria, i.e., a wide antimicrobial spectrum. Ticks do not transmit these pathogens, suggesting that these defensins may act against a wide range of microbes encountered by ticks during blood feeding or off-host periods. As demonstrated here, these I. ricinus defensins are also effective against the apicomplexan parasite Plasmodium falciparum. To study the general evolution of antimicrobial activity in tick defensins, the ancestral amino acid sequence of chelicerate defensins, which existed approximately 444 million years ago, was reconstructed using publicly available scorpion and tick defensin sequences (named Scorpions-Ticks Defensins Ancestor, STiDA). The activity of STiDA was tested against P. falciparum and the same Gram-negative and Gram-positive bacteria that were used for the I. ricinus defensins. While some extant tick defensins exhibit a wide antimicrobial spectrum, the ancestral defensin showed moderate activity against one of the tested microbes, P. falciparum. This study suggests that amino acid variability and defensin family expansion increased the antimicrobial spectrum of ancestral tick defensins.

UET: a database of evolutionarily-predicted functional determinants of protein sequences that cluster as functional sites in protein structures

The structure and function of proteins underlie most aspects of biology and their mutational perturbations often cause disease. To identify the molecular determinants of function as well as targets for drugs, it is central to characterize the important residues and how they cluster to form functional sites. The Evolutionary Trace (ET) achieves this by ranking the functional and structural importance of the protein sequence positions. ET uses evolutionary distances to estimate functional distances and correlates genotype variations with those in the fitness phenotype. Thus, ET ranks are worse for sequence positions that vary among evolutionarily closer homologs but better for positions that vary mostly among distant homologs. This approach identifies functional determinants, predicts function, guides the mutational redesign of functional and allosteric specificity, and interprets the action of coding sequence variations in proteins, people and populations. Now, the UET database offers pre-computed ET analyses for the protein structure databank, and on-the-fly analysis of any protein sequence. A web interface retrieves ET rankings of sequence positions and maps results to a structure to identify functionally important regions. This UET database integrates several ways of viewing the results on the protein sequence or structure and can be found at http://mammoth.bcm.tmc.edu/uet/.

Enamel matrix proteins exhibit growth factor activity: A review of evidence at the cellular and molecular levels

Enamel matrix derivative (EMD) is a commercially available protein extract, mainly comprising amelogenins. A number of other polypeptides have been identified in EMD, mostly growth factors, which promote cementogenesis and osteogenesis during the regeneration processes through the regulation of cell proliferation, differentiation and activity; however, not all of their functions are clear. Enamel extracts have been proposed to have numerous activities such as bone morphogenetic protein- and transforming growth factor β (TGF-β)-like activity, and activities similar to those of insulin-like growth factor, fibroblast growth factor, platelet-derived growth factor, vascular endothelial growth factor and epidermal growth factor. These activities have been observed at the molecular and cellular levels and in numerous animal models. Furthermore, it has been suggested that EMD contains an unidentified biologically active factor that acts in combination with TGF-β1, and several studies have reported functional similarities between growth factors and TGF-β in cellular processes. The effects of enamel extracts on the cell cycle and biology are summarized and discussed in this review.

DNA repair. PAXX, a paralog of XRCC4 and XLF, interacts with Ku to promote DNA double-strand break repair

XRCC4 and XLF are two structurally related proteins that function in DNA double-strand break (DSB) repair. Here, we identify human PAXX (PAralog of XRCC4 and XLF, also called C9orf142) as a new XRCC4 superfamily member and show that its crystal structure resembles that of XRCC4. PAXX interacts directly with the DSB-repair protein Ku and is recruited to DNA-damage sites in cells. Using RNA interference and CRISPR-Cas9 to generate PAXX(-/-) cells, we demonstrate that PAXX functions with XRCC4 and XLF to mediate DSB repair and cell survival in response to DSB-inducing agents. Finally, we reveal that PAXX promotes Ku-dependent DNA ligation in vitro and assembly of core nonhomologous end-joining (NHEJ) factors on damaged chromatin in cells. These findings identify PAXX as a new component of the NHEJ machinery.

Polymorphisms at amino acid residues 141 and 154 influence conformational variation in ovine PrP

Polymorphisms in ovine PrP at amino acid residues 141 and 154 are associated with susceptibility to ovine prion disease: Leu141Arg154 with classical scrapie and Phe141Arg154 and Leu141His154 with atypical scrapie. Classical scrapie is naturally transmissible between sheep, whereas this may not be the case with atypical scrapie. Critical amino acid residues will determine the range or stability of structural changes within the ovine prion protein or its functional interaction with potential cofactors, during conversion of PrPC to PrPSc in these different forms of scrapie disease. Here we computationally identified that regions of ovine PrP, including those near amino acid residues 141 and 154, displayed more conservation than expected based on local structural environment. Molecular dynamics simulations showed these conserved regions of ovine PrP displayed genotypic differences in conformational repertoire and amino acid side-chain interactions. Significantly, Leu141Arg154 PrP adopted an extended beta sheet arrangement in the N-terminal palindromic region more frequently than the Phe141Arg154 and Leu141His154 variants. We supported these computational observations experimentally using circular dichroism spectroscopy and immunobiochemical studies on ovine recombinant PrP. Collectively, our observations show amino acid residues 141 and 154 influence secondary structure and conformational change in ovine PrP that may correlate with different forms of scrapie.

Cube - an online tool for comparison and contrasting of protein sequences

When comparing sequences of similar proteins, two kinds of questions can be asked, and the related two kinds of inference made. First, one may ask to what degree they are similar, and then, how they differ. In the first case one may tentatively conclude that the conserved elements common to all sequences are of central and common importance to the protein's function. In the latter case the regions of specialization may be discriminative of the function or binding partners across subfamilies of related proteins. Experimental efforts - mutagenesis or pharmacological intervention - can then be pointed in either direction, depending on the context of the study.

Cube simplifies this process for users that already have their favorite sets of sequences, and helps them collate the information by visualization of the conservation and specialization scores on the sequence and on the structure, and by spreadsheet tabulation. All information can be visualized on the spot, or downloaded for reference and later inspection.

Server homepage: http://eopsf.org/cube

Recent advances in functional region prediction by using structural and evolutionary information - Remaining problems and future extensions

Structural genomics projects have solved many new structures with unknown functions. One strategy to investigate the function of a structure is to computationally find the functionally important residues or regions on it. Therefore, the development of functional region prediction methods has become an important research subject. An effective approach is to use a method employing structural and evolutionary information, such as the evolutionary trace (ET) method. ET ranks the residues of a protein structure by calculating the scores for relative evolutionary importance, and locates functionally important sites by identifying spatial clusters of highly ranked residues. After ET was developed, numerous ET-like methods were subsequently reported, and many of them are in practical use, although they require certain conditions. In this mini review, we first introduce the remaining problems and the recent improvements in the methods using structural and evolutionary information. We then summarize the recent developments of the methods. Finally, we conclude by describing possible extensions of the evolution- and structure-based methods.

Importin alpha subtypes determine differential transcription factor localization in embryonic stem cells maintenance

We recently demonstrated that the expression of the importin α subtype is switched from α2 to α1 during neural differentiation in mouse embryonic stem cells (ESCs) and that this switching has a major impact on cell differentiation. In this study, we report a cell-fate determination mechanism in which importin α2 negatively regulates the nuclear import of certain transcription factors to maintain ESC properties. The nuclear import of Oct6 and Brn2 was inhibited via the formation of a transport-incompetent complex of the cargo bound to a nuclear localization signal binding site in importin α2. Unless this dominant-negative effect was downregulated upon ESC differentiation, inappropriate cell death was induced. We propose that although certain transcription factors are necessary for differentiation in ESCs, these factors are retained in the cytoplasm by importin α2, thereby preventing transcription factor activity in the nucleus until the cells undergo differentiation.

Comparative Genomic and Sequence Analysis Provides Insight into the Molecular Functionality of NOD1 and NOD2

Amino acids with functional or key structural roles display higher degrees of conservation through evolution. The comparative analysis of protein sequences from multiple species and/or between homologous proteins can be highly informative in the identification of key structural and functional residues. Residues which in turn provide insight into the molecular mechanisms of protein function. We have explored the genomic and amino acid conservation of the prototypic innate immune genes NOD1 and NOD2. NOD1 orthologs were found in all vertebrate species analyzed, whilst NOD2 was absent from the genomes of avian, reptilian and amphibian species. Evolutionary trace analysis was used to identify highly conserved regions of NOD1 and NOD2 across multiple species. Consistent with the known functions of NOD1 and NOD2 highly conserved patches were identified that matched the Walker A and B motifs and provided interaction surfaces for the adaptor protein RIP2. Other patches of high conservation reflect key structural functions as predicted by homology models. In addition, the pattern of residue conservation within the leucine-rich repeat (LRR) region of NOD1 and NOD2 is indicative of a conserved mechanism of ligand recognition involving the concave surface of the LRRs.

The glaucoma-associated olfactomedin domain of myocilin is a novel calcium binding protein

Myocilin is a protein found in the trabecular meshwork extracellular matrix tissue of the eye that plays a role in regulating intraocular pressure. Both wild-type and certain myocilin variants containing mutations in the olfactomedin (OLF) domain are linked to the optic neuropathy glaucoma. Because calcium ions are important biological cofactors that play numerous roles in extracellular matrix proteins, we examined the calcium binding properties of the myocilin OLF domain (myoc-OLF). Our study reveals an unprecedented high affinity calcium binding site within myoc-OLF. The calcium ion remains bound to wild-type OLF at neutral and acidic pH. A glaucoma-causing OLF variant, myoc-OLF(D380A), is calcium-depleted. Key differences in secondary and tertiary structure between myoc-OLF(D380A) and wild-type myoc-OLF, as well as limited access to chelators, indicate that the calcium binding site is largely buried in the interior of the protein. Analysis of six conserved aspartate or glutamate residues and an additional 18 disease-causing variants revealed two other candidate residues that may be involved in calcium coordination. Our finding expands our knowledge of calcium binding in extracellular matrix proteins; provides new clues into domain structure, function, and pathogenesis for myocilin; and offers insights into highly conserved, biomedically relevant OLF domains.

The immunoglobulin domain of the sodium channel β3 subunit contains a surface-localized disulfide bond that is required for homophilic binding

The β subunits of voltage-gated sodium (Na_v) channels possess an extracellular immunoglobulin (Ig) domain that is related to the L1 family of cell-adhesion molecules (CAMs). Here we show that in HEK293 cells, secretion of the free Ig domain of the β3 subunit is reduced significantly when it is coexpressed with the full-length β3 and β1 subunits but not with the β2 subunit. Using immunoprecipitation, we show that the β3 subunit can mediate trans homophilic-binding via its Ig domain and that the β3-Ig domain can associate heterophilically with the β1 subunit. Evolutionary tracing analysis and structural modeling identified a cluster of surface-localized amino acids fully conserved between the Ig domains of all known β3 and β1 sequences. A notable feature of this conserved surface cluster is the presence of two adjacent cysteine residues that previously we have suggested may form a disulfide bond. We now confirm the presence of the disulfide bond in β3 using mass spectrometry, and we show that its integrity is essential for the association of the full-length, membrane-anchored β3 subunit with itself. However, selective reduction of this surface disulfide bond did not inhibit homophilic binding of the purified β3-Ig domain in free solution. Hence, the disulfide bond itself is unlikely to be part of the homophilic binding site. Rather, we suggest that its integrity ensures the Ig domain of the membrane-tethered β3 subunit adopts the correct orientation for productive association to occur in vivo.—Yereddi, N. R., Cusdin, F. S., Namadurai, S., Packman, L. C., Monie, T. P., Slavny, P., Clare, J. C., Powell, A. J., Jackson, A. P. The immunoglobulin domain of the sodium channel β3 subunit contains a surface-localized disulfide bond that is required for homophilic binding.

Structural and functional analysis of the NLRP4 pyrin domain

NLRP4 is a member of the nucleotide-binding and leucine-rich repeat receptor (NLR) family of cytosolic receptors and a member of an inflammation signaling cascade. Here, we present the crystal structure of the NLRP4 pyrin domain (PYD) at 2.3 Å resolution. The NLRP4 PYD is a member of the death domain (DD) superfamily and adopts a DD fold consisting of six α-helices tightly packed around a hydrophobic core, with a highly charged surface that is typical of PYDs. Importantly, however, we identified several differences between the NLRP4 PYD crystal structure and other PYD structures that are significant enough to affect NLRP4 function and its interactions with binding partners. Notably, the length of helix α3 and the α2−α3 connecting loop in the NLRP4 PYD are unique among PYDs. The apoptosis-associated speck-like protein containing a CARD (ASC) is an adaptor protein whose interactions with a number of distinct PYDs are believed to be critical for activation of the inflammatory response. Here, we use co-immunoprecipitation, yeast two-hybrid, and nuclear magnetic resonance chemical shift perturbation analysis to demonstrate that, despite being important for activation of the inflammatory response and sharing several similarities with other known ASC-interacting PYDs (i.e., ASC2), NLRP4 does not interact with the adaptor protein ASC. Thus, we propose that the factors governing homotypic PYD interactions are more complex than the currently accepted model, which states that complementary charged surfaces are the main determinants of PYD–PYD interaction specificity.

Functional region prediction with a set of appropriate homologous sequences--an index for sequence selection by integrating structure and sequence information with spatial statistics

Background

The detection of conserved residue clusters on a protein structure is one of the effective strategies for the prediction of functional protein regions. Various methods, such as Evolutionary Trace, have been developed based on this strategy. In such approaches, the conserved residues are identified through comparisons of homologous amino acid sequences. Therefore, the selection of homologous sequences is a critical step. It is empirically known that a certain degree of sequence divergence in the set of homologous sequences is required for the identification of conserved residues. However, the development of a method to select homologous sequences appropriate for the identification of conserved residues has not been sufficiently addressed. An objective and general method to select appropriate homologous sequences is desired for the efficient prediction of functional regions.

Results

We have developed a novel index to select the sequences appropriate for the identification of conserved residues, and implemented the index within our method to predict the functional regions of a protein. The implementation of the index improved the performance of the functional region prediction. The index represents the degree of conserved residue clustering on the tertiary structure of the protein. For this purpose, the structure and sequence information were integrated within the index by the application of spatial statistics. Spatial statistics is a field of statistics in which not only the attributes but also the geometrical coordinates of the data are considered simultaneously. Higher degrees of clustering generate larger index scores. We adopted the set of homologous sequences with the highest index score, under the assumption that the best prediction accuracy is obtained when the degree of clustering is the maximum. The set of sequences selected by the index led to higher functional region prediction performance than the sets of sequences selected by other sequence-based methods.

Conclusions

Appropriate homologous sequences are selected automatically and objectively by the index. Such sequence selection improved the performance of functional region prediction. As far as we know, this is the first approach in which spatial statistics have been applied to protein analyses. Such integration of structure and sequence information would be useful for other bioinformatics problems.

Characterization of follistatin-type domains and their contribution to myostatin and activin A antagonism

Follistatin (FST)-type proteins are important antagonists of some members of the large TGF-β family of cytokines. These include myostatin, an important negative regulator of muscle growth, and the closely related activin A, which is involved in many physiological functions, including maintenance of a normal reproductive axis. FST-type proteins, including FST and FST-like 3 (FSTL3), differentially inhibit various TGF-β family ligands by binding each ligand with two FST-type molecules. In this study, we sought to examine features that are important for ligand antagonism by FST-type proteins. Previous work has shown that a modified construct consisting of the FST N-terminal domain (ND) followed by two repeating follistatin domains (FSD), herein called FST ND-FSD1-FSD1, exhibits strong specificity for myostatin over activin A. Using cell-based assays, we show that FST ND-FSD1-FSD1 is unique in its specificity for myostatin as compared with similar constructs containing domains from FSTL3 and that the ND is critical to its activity. Furthermore, we demonstrate that FSD3 of FST provides affinity to ligand inhibition and confers resistance to perturbations in the ND and FSD2, likely through the interaction of FSD3 of one FST molecule with the ND of the other FST molecule. Additionally, our data suggest that this contact provides cooperativity to ligand antagonism. Cross-linking studies show that this interaction also potentiates formation of 1:2 ligand-FST complexes, whereas lack of FSD3 allows formation of 1:1 complexes. Altogether, these studies support that domain differences generate FST-type molecules that are each uniquely suited ligand antagonists.

Analyzing effects of naturally occurring missense mutations

Single-point mutation in genome, for example, single-nucleotide polymorphism (SNP) or rare genetic mutation, is the change of a single nucleotide for another in the genome sequence. Some of them will produce an amino acid substitution in the corresponding protein sequence (missense mutations); others will not. This paper focuses on genetic mutations resulting in a change in the amino acid sequence of the corresponding protein and how to assess their effects on protein wild-type characteristics. The existing methods and approaches for predicting the effects of mutation on protein stability, structure, and dynamics are outlined and discussed with respect to their underlying principles. Available resources, either as stand-alone applications or webservers, are pointed out as well. It is emphasized that understanding the molecular mechanisms behind these effects due to these missense mutations is of critical importance for detecting disease-causing mutations. The paper provides several examples of the application of 3D structure-based methods to model the effects of protein stability and protein-protein interactions caused by missense mutations as well.

Structure of myostatin·follistatin-like 3: N-terminal domains of follistatin-type molecules exhibit alternate modes of binding

TGF-β family ligands are involved in a variety of critical physiological processes. For instance, the TGF-β ligand myostatin is a staunch negative regulator of muscle growth and a therapeutic target for muscle-wasting disorders. Therefore, it is important to understand the molecular mechanisms of TGF-β family regulation. One form of regulation is through inhibition by extracellular antagonists such as the follistatin (Fst)-type proteins. Myostatin is tightly controlled by Fst-like 3 (Fstl3), which is the only Fst-type molecule that has been identified in the serum bound to myostatin. Here, we present the crystal structure of myostatin in complex with Fstl3. The structure reveals that the N-terminal domain (ND) of Fstl3 interacts uniquely with myostatin as compared with activin A, because it utilizes different surfaces on the ligand. This results in conformational differences in the ND of Fstl3 that alter its position in the type I receptor-binding site of the ligand. We also show that single point mutations in the ND of Fstl3 are detrimental to ligand binding, whereas corresponding mutations in Fst have little effect. Overall, we have shown that the NDs of Fst-type molecules exhibit distinctive modes of ligand binding, which may affect overall affinity of ligand·Fst-type protein complexes.

Energy deprivation alters in a leptin- and cortisol-independent manner circulating levels of activin A and follistatin but not myostatin in healthy males

CONTEXT

Activin A, myostatin, and follistatin have recently emerged as important regulatory molecules of reproduction and the musculoskeletal system. Little is known, however, about their day/night patterns of secretion and their physiological regulation by energy availability.

OBJECTIVE

The objective of the study was to explore day/night patterns of secretion and assess whether energy deprivation alters circulating levels of activin A, myostatin, follistatin, and cortisol and to examine whether leptin may mediate this effect. DESIGN, SETTING AND PATIENTS, AND INTERVENTIONS: Seven healthy lean men (aged 23.2 ± 3.7 yr, body mass index 23.6 ± 1.7 kg/m(2)) were studied for 72 h under three different conditions: on their baseline/isocaloric diet and in a complete fasting state with administration of either placebo or metreleptin. The two fasting studies were randomized and double blinded. Blood samples were obtained every 15 min from 0800 h on d 3 until 0800 h on d 4 and pooled hourly.

MAIN OUTCOME MEASURES

Serum concentrations of activin A, myostatin, follistatin, cortisol, and leptin were measured.

RESULTS

In contrast to cortisol, we demonstrated no day/night pattern of activin A, myostatin, and follistatin secretion. Activin A concentrations decreased significantly in response to energy deprivation (P < 0.01). Follistatin and cortisol concentrations increased significantly (P < 0.01 and P < 0.01, respectively). Myostatin remained unaffected (P = 0.40). Leptin administration reversed cortisol response (P < 0.01) but failed to alter activin A, follistatin, or myostatin concentrations.

CONCLUSIONS

Unlike cortisol, there is no day/night variation in the concentrations of activin A, myostatin, and follistatin in healthy young males. Although energy deprivation-induced cortisol changes are leptin mediated, the changes in follistatin and activin A concentrations occur through a leptin-independent pathway.

Evolutionary traces decode molecular mechanism behind fast pace of myosin XI

Background

Cytoplasmic class XI myosins are the fastest processive motors known. This class functions in high-velocity cytoplasmic streaming in various plant cells from algae to angiosperms. The velocities at which they process are ten times faster than its closest class V homologues.

Results

To provide sequence determinants and structural rationale for the molecular mechanism of this fast pace myosin, we have compared the sequences from myosin class V and XI through Evolutionary Trace (ET) analysis. The current study identifies class-specific residues of myosin XI spread over the actin binding site, ATP binding site and light chain binding neck region. Sequences for ET analysis were accumulated from six plant genomes, using literature based text search and sequence searches, followed by triple validation viz. CDD search, string-based searches and phylogenetic clustering. We have identified nine myosin XI genes in sorghum and seven in grape by sequence searches. Both the plants possess one gene product each belonging to myosin type VIII as well. During this process, we have re-defined the gene boundaries for three sorghum myosin XI genes using fgenesh program.

Conclusion

Molecular modelling and subsequent analysis of putative interactions involving these class-specific residues suggest a structural basis for the molecular mechanism behind high velocity of plant myosin XI. We propose a model of a more flexible switch I region that contributes to faster ADP release leading to high velocity movement of the algal myosin XI.

Determinants, discriminants, conserved residues--a heuristic approach to detection of functional divergence in protein families

In this work, belonging to the field of comparative analysis of protein sequences, we focus on detection of functional specialization on the residue level. As the input, we take a set of sequences divided into groups of orthologues, each group known to be responsible for a different function. This provides two independent pieces of information: within group conservation and overlap in amino acid type across groups. We build our discussion around the set of scoring functions that keep the two separated and the source of the signal easy to trace back to its source.We propose a heuristic description of functional divergence that includes residue type exchangeability, both in the conservation and in the overlap measure, and does not make any assumptions on the rate of evolution in the groups other than the one under consideration. Residue types acceptable at a certain position within an orthologous group are described as a distribution which evolves in time, starting from a single ancestral type, and is subject to constraints that can be inferred only indirectly. To estimate the strength of the constraints, we compare the observed degrees of conservation and overlap with those expected in the hypothetical case of a freely evolving distribution.Our description matches the experiment well, but we also conclude that any attempt to capture the evolutionary behavior of specificity determining residues in terms of a scalar function will be tentative, because no single model can cover the variety of evolutionary behavior such residues exhibit. Especially, models expecting the same type of evolutionary behavior across functionally divergent groups tend to miss a portion of information otherwise retrievable by the conservation and overlap measures they use.

The pyrokinin/ pheromone biosynthesis-activating neuropeptide (PBAN) family of peptides and their receptors in Insecta: evolutionary trace indicates potential receptor ligand-binding domains

The pyrokinin/pheromone biosynthesis-activating neuropeptide (PBAN) family of G-protein-coupled receptors and their ligands have been identified in various insects. Physiological functions of pyrokinin peptides include muscle contraction, whereas PBAN regulates, among other functions, pheromone production in moths which indicates the pleiotropic nature of these peptides. Based on the alignment of annotated genomic sequences, the pyrokinin/PBAN family of receptors have similarity with the corresponding structures of the capa or periviscerokinin receptors of insects and the neuromedin U receptors of vertebrates. In our study, evolutionary trace (ET) analysis on the insect receptor sequences was conducted to predict the putative ligand recognition and binding sites. The ET analysis of four class-specific receptors indicated several amino acid residues that are conserved in the transmembrane domains. The receptor extracellular domains exhibit several class-specific amino acid residues, which could indicate putative domains for activation of these receptors by ligand recognition and binding.

Molecular modeling and in silico characterization of Mycobacterium tuberculosis TlyA: possible misannotation of this tubercle bacilli-hemolysin

Background

The TlyA protein has a controversial function as a virulence factor in Mycobacterium tuberculosis (M. tuberculosis). At present, its dual activity as hemolysin and RNA methyltransferase in M. tuberculosis has been indirectly proposed based on in vitro results. There is no evidence however for TlyA relevance in the survival of tubercle bacilli inside host cells or whether both activities are functionally linked. A thorough analysis of structure prediction for this mycobacterial protein in this study shows the need for reevaluating TlyA's function in virulence.

Results

Bioinformatics analysis of TlyA identified a ribosomal protein binding domain (S4 domain), located between residues 5 and 68 as well as an FtsJ-like methyltranferase domain encompassing residues 62 and 247, all of which have been previously described in translation machinery-associated proteins. Subcellular localization prediction showed that TlyA lacks a signal peptide and its hydrophobicity profile showed no evidence of transmembrane helices. These findings suggested that it may not be attached to the membrane, which is consistent with a cytoplasmic localization. Three-dimensional modeling of TlyA showed a consensus structure, having a common core formed by a six-stranded β-sheet between two α-helix layers, which is consistent with an RNA methyltransferase structure. Phylogenetic analyses showed high conservation of the tlyA gene among Mycobacterium species. Additionally, the nucleotide substitution rates suggested purifying selection during tlyA gene evolution and the absence of a common ancestor between TlyA proteins and bacterial pore-forming proteins.

Conclusion

Altogether, our manual in silico curation suggested that TlyA is involved in ribosomal biogenesis and that there is a functional annotation error regarding this protein family in several microbial and plant genomes, including the M. tuberculosis genome.

Feature amplified voting algorithm for functional analysis of protein superfamily

BACKGROUND

Identifying the regions associated with protein function is a singularly important task in the post-genomic era. Biological studies often identify functional enzyme residues by amino acid sequences, particularly when related structural information is unavailable. In some cases of protein superfamilies, functional residues are difficult to detect by current alignment tools or evolutionary strategies when phylogenetic relationships do not parallel their protein functions. The solution proposed in this study is Feature Amplified Voting Algorithm with Three-profile alignment (FAVAT). The core concept of FAVAT is to reveal the desired features of a target enzyme or protein by voting on three different property groups aligned by three-profile alignment method. Functional residues of a target protein can then be retrieved by FAVAT analysis. In this study, the amidohydrolase superfamily was an interesting case for verifying the proposed approach because it contains divergent enzymes and proteins.

RESULTS

The FAVAT was used to identify critical residues of mammalian imidase, a member of the amidohydrolase superfamily. Members of this superfamily were first classified by their functional properties and sources of original organisms. After FAVAT analysis, candidate residues were identified and compared to a bacterial hydantoinase in which the crystal structure (1GKQ) has been fully elucidated. One modified lysine, three histidines and one aspartate were found to participate in the coordination of metal ions in the active site. The FAVAT analysis also redressed the misrecognition of metal coordinator Asp57 by the multiple sequence alignment (MSA) method. Several other amino acid residues known to be related to the function or structure of mammalian imidase were also identified.

CONCLUSIONS

The FAVAT is shown to predict functionally important amino acids in amidohydrolase superfamily. This strategy effectively identifies functionally important residues by analyzing the discrepancy between the sequence and functional properties of related proteins in a superfamily, and it should be applicable to other protein families.

Crystal structures of mutant forms of the yeast F1 ATPase reveal two modes of uncoupling

The mitochondrial ATP synthase couples the flow of protons with the phosphorylation of ADP. A class of mutations, the mitochondrial genome integrity (mgi) mutations, has been shown to uncouple this process in the yeast mitochondrial ATP synthase. Four mutant forms of the yeast F(1) ATPase with mgi mutations were crystallized; the structures were solved and analyzed. The analysis identifies two mechanisms of structural uncoupling: one in which the empty catalytic site is altered and in doing so, apparently disrupts substrate (phosphate) binding, and a second where the steric hindrance predicted between γLeu83 and β(DP) residues, Leu-391 and Glu-395, located in Catch 2 region, is reduced allowing rotation of the γ-subunit with less impedance. Overall, the structures provide key insights into the critical interactions in the yeast ATP synthase involved in the coupling process.

Rapid evolution of a novel signalling mechanism by concerted duplication and divergence of a BMP ligand and its extracellular modulators

Gene duplication and divergence is widely considered to be a fundamental mechanism for generating evolutionary novelties. The Bone Morphogenetic Proteins (BMPs) are a diverse family of signalling molecules found in all metazoan genomes that have evolved by duplication and divergence from a small number of ancestral types. In the fruit fly Drosophila, there are three BMPs: Decapentaplegic (Dpp) and Glass bottom boat (Gbb), which are the orthologues of vertebrate BMP2/4 and BMP5/6/7/8, respectively, and Screw (Scw), which, at the sequence level, is equally divergent from Dpp and Gbb. It has recently been shown that Scw has arisen from a duplication of Gbb in the lineage leading to higher Diptera. We show that since this duplication event, Gbb has maintained the ancestral BMP5/6/7/8 functionality while Scw has rapidly diverged. The evolution of Scw was accompanied by duplication and divergence of a suite of extracellular regulators that continue to diverge together in the higher Diptera. In addition, Scw has become restricted in its receptor specificity: Gbb proteins can signal through the Type I receptors Thick veins (Tkv) and Saxophone (Sax), while Scw signals through Sax. Thus, in a relatively short span of evolutionary time, the duplication event that gave rise to Scw produced not only a novel ligand but also a novel signalling mode that is functionally distinct from the ancestral Gbb mode. Our results demonstrate the plasticity of the BMP pathway not only in evolving new family members and new functions but also new signalling modes by redeploying key regulators in the pathway.

Coevolution within a transcriptional network by compensatory trans and cis mutations

Transcriptional networks have been shown to evolve very rapidly, prompting questions as to how such changes arise and are tolerated. Recent comparisons of transcriptional networks across species have implicated variations in the cis-acting DNA sequences near genes as the main cause of divergence. What is less clear is how these changes interact with trans-acting changes occurring elsewhere in the genetic circuit. Here, we report the discovery of a system of compensatory trans and cis mutations in the yeast AP-1 transcriptional network that allows for conserved transcriptional regulation despite continued genetic change. We pinpoint a single species, the fungal pathogen Candida glabrata, in which a trans mutation has occurred very recently in a single AP-1 family member, distinguishing it from its Saccharomyces ortholog. Comparison of chromatin immunoprecipitation profiles between Candida and Saccharomyces shows that, despite their different DNA-binding domains, the AP-1 orthologs regulate a conserved block of genes. This conservation is enabled by concomitant changes in the cis-regulatory motifs upstream of each gene. Thus, both trans and cis mutations have perturbed the yeast AP-1 regulatory system in such a way as to compensate for one another. This demonstrates an example of "coevolution" between a DNA-binding transcription factor and its cis-regulatory site, reminiscent of the coevolution of protein binding partners.

Evolution of a signaling nexus constrained by protein interfaces and conformational States

Heterotrimeric G proteins act as the physical nexus between numerous receptors that respond to extracellular signals and proteins that drive the cytoplasmic response. The Gα subunit of the G protein, in particular, is highly constrained due to its many interactions with proteins that control or react to its conformational state. Various organisms contain differing sets of Gα-interacting proteins, clearly indicating that shifts in sequence and associated Gα functionality were acquired over time. These numerous interactions constrained much of Gα evolution; yet Gα has diversified, through poorly understood processes, into several functionally specialized classes, each with a unique set of interacting proteins. Applying a synthetic sequence-based approach to mammalian Gα subunits, we established a set of seventy-five evolutionarily important class-distinctive residues, sites where a single Gα class is differentiated from the three other classes. We tested the hypothesis that shifts at these sites are important for class-specific functionality. Importantly, we mapped known and well-studied class-specific functionalities from all four mammalian classes to sixteen of our class-distinctive sites, validating the hypothesis. Our results show how unique functionality can evolve through the recruitment of residues that were ancestrally functional. We also studied acquisition of functionalities by following these evolutionarily important sites in non-mammalian organisms. Our results suggest that many class-distinctive sites were established early on in eukaryotic diversification and were critical for the establishment of new Gα classes, whereas others arose in punctuated bursts throughout metazoan evolution. These Gα class-distinctive residues are rational targets for future structural and functional studies.

Role of Hsp70 ATPase domain intrinsic dynamics and sequence evolution in enabling its functional interactions with NEFs

Catalysis of ADP-ATP exchange by nucleotide exchange factors (NEFs) is central to the activity of Hsp70 molecular chaperones. Yet, the mechanism of interaction of this family of chaperones with NEFs is not well understood in the context of the sequence evolution and structural dynamics of Hsp70 ATPase domains. We studied the interactions of Hsp70 ATPase domains with four different NEFs on the basis of the evolutionary trace and co-evolution of the ATPase domain sequence, combined with elastic network modeling of the collective dynamics of the complexes. Our study reveals a subtle balance between the intrinsic (to the ATPase domain) and specific (to interactions with NEFs) mechanisms shared by the four complexes. Two classes of key residues are distinguished in the Hsp70 ATPase domain: (i) highly conserved residues, involved in nucleotide binding, which mediate, via a global hinge-bending, the ATPase domain opening irrespective of NEF binding, and (ii) not-conserved but co-evolved and highly mobile residues, engaged in specific interactions with NEFs (e.g., N57, R258, R262, E283, D285). The observed interplay between these respective intrinsic (pre-existing, structure-encoded) and specific (co-evolved, sequence-dependent) interactions provides us with insights into the allosteric dynamics and functional evolution of the modular Hsp70 ATPase domain.

The heat shock protein 70 (Hsp70) serves as a housekeeper in the cell, assisting in the correct folding, trafficking, and degradation of many proteins. The ATPase domain is the control unit of this molecular machine and its efficient functioning requires interactions with co-chaperones, including, in particular, the nucleotide exchange factors (NEFs). We examined the molecular motions of the ATPase domain in both NEF-bound and -unbound forms. We found that the NEF-binding surface enjoys large global movements prior to NEF binding, which presumably facilitates NEF recognition and binding. NEF binding stabilizes the ATPase domain in an open form and thereby facilitates the nucleotide exchange step of the chaperone cycle. A series of highly correlated amino acids were distinguished at the NEF-binding sites of the Hsp70 ATPase domain, which highlights the adaptability of the ATPase domain, both structurally and sequentially, to recognize NEFs. In contrast, the nucleotide-binding residues are tightly held near a global hinge center and are highly conserved. The contrasting properties of these two groups of residues point to an evolutionarily optimized balance between conserved/constrained and co-evolved/mobile amino acids, which enables the functional interactions of the modular Hps70 ATPase domains with NEFs.

Bone morphogenetic protein and growth differentiation factor cytokine families and their protein antagonists

The BMPs (bone morphogenetic proteins) and the GDFs (growth and differentiation factors) together form a single family of cystine-knot cytokines, sharing the characteristic fold of the TGFbeta (transforming growth factor-beta) superfamily. Besides the ability to induce bone formation, which gave the BMPs their name, the BMP/GDFs display morphogenetic activities in the development of a wide range of tissues. BMP/GDF homo- and hetero-dimers interact with combinations of type I and type II receptor dimers to produce multiple possible signalling complexes, leading to the activation of one of two competing sets of SMAD transcription factors. BMP/GDFs have highly specific and localized functions. These are regulated in a number of ways, including the developmental restriction of BMP/GDF expression and through the secretion of several specific BMP antagonist proteins that bind with high affinity to the cytokines. Curiously, a number of these antagonists are also members of the TGF-beta superfamily. Finally a number of both the BMP/GDFs and their antagonists interact with the heparan sulphate side chains of cell-surface and extracellular-matrix proteoglycans.

In silico point mutation and evolutionary trace analysis applied to nicotinic acetylcholine receptors in deciphering ligand-binding surfaces

The nicotinic acetylcholine receptors (nAChRs) are members of the Cys-loop superfamily and contain ligand gated ion channels (LGIC). These receptors are located mostly in the central nervous system (CNS) and peripheral nervous system (PNS). nAChRs reside at pre-synaptic regions to mediate acetylcholine neurotransmission and in the post synaptic membrane to propagate nerve impulses through neurons via acetylcholine. Malfunction of this neurotransmitter receptor is believed to cause various neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease and schizophrenia, and nAChRs are thus important drug targets. In the present work, starting from an earlier model of pentameric alpha7nAChR, a considerable effort has been taken to investigate interaction with ligands by performing docking studies with a diverse array of agonists and antagonists. Analysis of these docking complexes reveals identification of possible ligand-interacting residues. Some of these residues, e.g. Ser34, Gln55, Ser146, and Tyr166, which are evolutionarily conserved, were specifically subjected to virtual mutations based on their amino acid properties and found to be highly sensitive in the presence of antagonists by docking. Further, the study was extended using evolutionary trace analysis, revealing conserved and class-specific residues close to the putative ligand-binding site, further supporting the results of docking experiments.

Deciphering chemotaxis pathways using cross species comparisons

Background

Chemotaxis is the process by which motile bacteria sense their chemical environment and move towards more favourable conditions. Escherichia coli utilises a single sensory pathway, but little is known about signalling pathways in species with more complex systems.

Results

To investigate whether chemotaxis pathways in other bacteria follow the E. coli paradigm, we analysed 206 species encoding at least 1 homologue of each of the 5 core chemotaxis proteins (CheA, CheB, CheR, CheW and CheY). 61 species encode more than one of all of these 5 proteins, suggesting they have multiple chemotaxis pathways. Operon information is not available for most bacteria, so we developed a novel statistical approach to cluster che genes into putative operons. Using operon-based models, we reconstructed putative chemotaxis pathways for all 206 species. We show that cheA-cheW and cheR-cheB have strong preferences to occur in the same operon as two-gene blocks, which may reflect a functional requirement for co-transcription. However, other che genes, most notably cheY, are more dispersed on the genome. Comparison of our operons with shuffled equivalents demonstrates that specific patterns of genomic location may be a determining factor for the observed in vivo chemotaxis pathways.

We then examined the chemotaxis pathways of Rhodobacter sphaeroides. Here, the PpfA protein is known to be critical for correct partitioning of proteins in the cytoplasmically-localised pathway. We found ppfA in che operons of many species, suggesting that partitioning of cytoplasmic Che protein clusters is common. We also examined the apparently non-typical chemotaxis components, CheA3, CheA4 and CheY6. We found that though variants of CheA proteins are rare, the CheY6 variant may be a common type of CheY, with a significantly disordered C-terminal region which may be functionally significant.

Conclusions

We find that many bacterial species potentially have multiple chemotaxis pathways, with grouping of che genes into operons likely to be a major factor in keeping signalling pathways distinct. Gene order is highly conserved with cheA-cheW and cheR-cheB blocks, perhaps reflecting functional linkage. CheY behaves differently to other Che proteins, both in its genomic location and its putative protein interactions, which should be considered when modelling chemotaxis pathways.