Evolution of substrate specificity in a recipient's enzyme following horizontal gene transfer

Despite the prominent role of horizontal gene transfer (HGT) in shaping bacterial metabolism, little is known about the impact of HGT on the evolution of enzyme function. Specifically, what is the influence of a recently acquired gene on the function of an existing gene? For example, certain members of the genus Corynebacterium have horizontally acquired a whole L-tryptophan biosynthetic operon, whereas in certain closely related actinobacteria, for example, Mycobacterium, the trpF gene is missing. In Mycobacterium, the function of the trpF gene is performed by a dual-substrate (βα)8 phosphoribosyl isomerase (priA gene) also involved in L-histidine (hisA gene) biosynthesis. We investigated the effect of a HGT-acquired TrpF enzyme upon PriA’s substrate specificity in Corynebacterium through comparative genomics and phylogenetic reconstructions. After comprehensive in vivo and enzyme kinetic analyses of selected PriA homologs, a novel (βα)8 isomerase subfamily with a specialized function in L-histidine biosynthesis, termed subHisA, was confirmed. X-ray crystallography was used to reveal active-site mutations in subHisA important for narrowing of substrate specificity, which when mutated to the naturally occurring amino acid in PriA led to gain of function. Moreover, in silico molecular dynamic analyses demonstrated that the narrowing of substrate specificity of subHisA is concomitant with loss of ancestral protein conformational states. Our results show the importance of HGT in shaping enzyme evolution and metabolism

GeMMA: functional subfamily classification within superfamilies of predicted protein structural domains

GeMMA (Genome Modelling and Model Annotation) is a new approach to automatic functional subfamily classification within families and superfamilies of protein sequences. A major advantage of GeMMA is its ability to subclassify very large and diverse superfamilies with tens of thousands of members, without the need for an initial multiple sequence alignment. Its performance is shown to be comparable to the established high-performance method SCI-PHY. GeMMA follows an agglomerative clustering protocol that uses existing software for sensitive and accurate multiple sequence alignment and profile–profile comparison. The produced subfamilies are shown to be equivalent in quality whether whole protein sequences are used or just the sequences of component predicted structural domains. A faster, heuristic version of GeMMA that also uses distributed computing is shown to maintain the performance levels of the original implementation. The use of GeMMA to increase the functional annotation coverage of functionally diverse Pfam families is demonstrated. It is further shown how GeMMA clusters can help to predict the impact of experimentally determining a protein domain structure on comparative protein modelling coverage, in the context of structural genomics

Non-monotonic variation with salt concentration of the second virial coefficient in protein solutions

The osmotic virial coefficient $B_2$ of globular protein solutions is calculated as a function of added salt concentration at fixed pH by computer simulations of the ``primitive model''. The salt and counter-ions as well as a discrete charge pattern on the protein surface are explicitly incorporated. For parameters roughly corresponding to lysozyme, we find that $B_2$ first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength. Our studies demonstrate that the existence of a discrete charge pattern on the protein surface profoundly influences the effective interactions and that non-linear Poisson Boltzmann and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory fail for large ionic strength. The observed non-monotonicity of $B_2$ is compared to experiments. Implications for protein crystallization are discussed.Comment: 43 pages, including 17 figure

MACSIMS : multiple alignment of complete sequences information management system

BACKGROUND: In the post-genomic era, systems-level studies are being performed that seek to explain complex biological systems by integrating diverse resources from fields such as genomics, proteomics or transcriptomics. New information management systems are now needed for the collection, validation and analysis of the vast amount of heterogeneous data available. Multiple alignments of complete sequences provide an ideal environment for the integration of this information in the context of the protein family. RESULTS: MACSIMS is a multiple alignment-based information management program that combines the advantages of both knowledge-based and ab initio sequence analysis methods. Structural and functional information is retrieved automatically from the public databases. In the multiple alignment, homologous regions are identified and the retrieved data is evaluated and propagated from known to unknown sequences with these reliable regions. In a large-scale evaluation, the specificity of the propagated sequence features is estimated to be >99%, i.e. very few false positive predictions are made. MACSIMS is then used to characterise mutations in a test set of 100 proteins that are known to be involved in human genetic diseases. The number of sequence features associated with these proteins was increased by 60%, compared to the features available in the public databases. An XML format output file allows automatic parsing of the MACSIM results, while a graphical display using the JalView program allows manual analysis. CONCLUSION: MACSIMS is a new information management system that incorporates detailed analyses of protein families at the structural, functional and evolutionary levels. MACSIMS thus provides a unique environment that facilitates knowledge extraction and the presentation of the most pertinent information to the biologist. A web server and the source code are available at

Evolutionary origins of the estrogen signaling system : insights from amphioxus

Author Posting. © The Author(s), 2011. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Journal of Steroid Biochemistry and Molecular Biology 127 (2011): 176–188, doi:10.1016/j.jsbmb.2011.03.022.Classically, the estrogen signaling system has two core components: cytochrome P450 aromatase (CYP19), the enzyme complex that catalyzes the rate limiting step in estrogen biosynthesis; and estrogen receptors (ERs), ligand activated transcription factors that interact with the regulatory region of target genes to mediate the biological effects of estrogen. While the importance of estrogens for regulation of reproduction, development and physiology has been well-documented in gnathostome vertebrates, the evolutionary origins of estrogen as a hormone are still unclear. As invertebrates within the phylum Chordata, cephalochordates (e.g. the amphioxus of the genus Branchiostoma) are among the closest invertebrate relatives of the vertebrates and can provide critical insight into the evolution of vertebrate-specific molecules and pathways. To address this question, this paper briefly reviews relevant earlier studies that help to illuminate the history of the aromatase and ER genes, with a particular emphasis on insights from amphioxus and other invertebrates. We then present new analyses of amphioxus aromatase and ER sequence and function, including an in silico model of the amphioxus aromatase protein, and CYP19 gene analysis. CYP19 shares a conserved gene structure with vertebrates (9 coding exons) and moderate sequence conservation (40% amino acid identity with human CYP19). Modeling of the amphioxus aromatase substrate binding site and simulated docking of androstenedione in comparison to the human aromatase shows that the substrate binding site is conserved and predicts that androstenedione could be a substrate for amphioxus CYP19. The amphioxus ER is structurally similar to vertebrate ERs, but differs in sequence and key residues of the ligand binding domain. Consistent with results from other laboratories, amphioxus ER did not bind radiolabeled estradiol, nor did it modulate gene expression on an estrogen-responsive element (ERE) in the presence 59 of estradiol, 4-hydroxytamoxifen, diethylstilbestrol, bisphenol A or genistein. Interestingly, it has been shown that a related gene, the amphioxus “steroid receptor” (SR), can be activated by estrogens and that amphioxus ER can repress this activation. CYP19, ER and SR are all primarily expressed in gonadal tissue, suggesting an ancient paracrine/autocrinesignaling role, but it is not yet known how their expression is regulated and, if estrogen is actually synthesized in amphioxus, whether it has a role in mediating any biological effects . Functional studies are clearly needed to link emerging bioinformatics and in vitro molecular biology results with organismal physiology to develop an understanding of the evolution of estrogen signaling.Supported by grants from the NIEHS P42 ES07381 (GVC, SV) and EPA (STAR-RD831301) (GVC), a Ruth L Kirschstein National Research Service Award (AT, F32 ES013092-01), an NIH traineeship (SS, SG), a NATO Fellowship (AN) and the Boston University Undergraduate Research Program (LC)

A Comparative Structural Bioinformatics Analysis of the Insulin Receptor Family Ectodomain Based on Phylogenetic Information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals

Habitat partitioning and vulnerability of sharks in the Great Barrier Reef Marine Park

Sharks present a critical conservation challenge, but little is known about their spatial distribution and vulnerability, particularly in complex seascapes such as Australia's Great Barrier Reef Marine Park (GBRMP). We review (1) the distribution of shark species among the primary habitats of the GBRMP (coral reefs, inshore/shelf, pelagic and deep-water habitats) (2) the relative exploitation of each species by fisheries, and (3) how current catch rates interact with their vulnerability and trophic index. Excluding rays and chimaeras, we identify a total of 82 shark species in the GBRMP. We find that shark research in the GBRMP has yielded little quantitative information on most species. Reef sharks are largely site-fidelic, but can move large distances and some regularly use non-reef habitats. Inshore and shelf sharks use coastal habitats either exclusively or during specific times in their life cycle (e.g. as nurseries). Virtually nothing is known about the distribution and habitat use of the GBRMP's pelagic and deep-water sharks. At least 46 species (53.5 %) are caught in one or more fisheries, but stock assessments are lacking for most. At least 17 of the sharks caught are considered highly vulnerable to exploitation. We argue that users of shark resources should be responsible for demonstrating that a fishery is sustainable before exploitation is allowed to commence or continue. This fundamental change in management principle will safeguard against stock collapses that have characterised many shark fisheries

FuSe: a tool to move RNA-Seq analyses from chromosomal/gene loci to functional grouping of mRNA transcripts

Typical RNA sequencing (RNA-Seq) analyses are performed either at the gene level by summing all reads from the same locus, assuming that all transcripts from a gene make a protein or at the transcript level, assuming that each transcript displays unique function. However, these assumptions are flawed, as a gene can code for different types of transcripts and different transcripts are capable of synthesizing similar, different or no protein. As a consequence, functional changes are not well illustrated by either gene or transcript analyses. We propose to improve RNA-Seq analyses by grouping the transcripts based on their similar functions. We developed FuSe to predict functional similarities using the primary and secondary structure of proteins. To estimate the likelihood of proteins with similar functions, FuSe computes two confidence scores: knowledge (KS) and discovery (DS) for protein pairs. Overlapping protein pairs exhibiting high confidence are grouped to form 'similar function protein groups' and expression is calculated for each functional group. The impact of using FuSe is demonstrated on in vitro cells exposed to paracetamol, which highlight genes responsible for cell adhesion and glycogen regulation which were earlier shown to be not differentially expressed with traditional analysis methods