Identification of key amino acid residues in AtUMAMIT29 for transport of glucosinolates

Glucosinolates are key defense compounds of plants in Brassicales order, and their accumulation in seeds is essential for the protection of the next generation. Recently, members of the Usually Multiple Amino acids Move In and Out Transporter (UMAMIT) family were shown to be essential for facilitating transport of seed-bound glucosinolates from site of synthesis within the reproductive organ to seeds. Here, we set out to identify amino acid residues responsible for glucosinolate transport activity of the main seed glucosinolate exporter UMAMIT29 in Arabidopsis thaliana. Based on a predicted model of UMAMIT29, we propose that the substrate transporting cavity consists of 51 residues, of which four are highly conserved residues across all the analyzed homologs of UMAMIT29. A comparison of the putative substrate binding site of homologs within the brassicaceous-specific, glucosinolate-transporting clade with the non-brassicaceous-specific, non-glucosinolate-transporting UMAMIT32 clade identified 11 differentially conserved sites. When each of the 11 residues of UMAMIT29 was individually mutated into the corresponding residue in UMAMIT32, five mutant variants (UMAMIT29#V27F, UMAMIT29#M86V, UMAMIT29#L109V, UMAMIT29#Q263S, and UMAMIT29#T267Y) reduced glucosinolate transport activity over 75% compared to wild-type UMAMIT29. This suggests that these residues are key for UMAMIT29-mediated glucosinolate transport activity and thus potential targets for blocking the transport of glucosinolates to the seeds.

Peptidoglycan recycling mediated by an ABC transporter in the plant pathogen Agrobacterium tumefaciens

During growth and division, the bacterial cell wall peptidoglycan (PG) is remodelled, resulting in the liberation of PG muropeptides which are typically reinternalized and recycled. Bacteria belonging to the Rhizobiales and Rhodobacterales orders of the Alphaproteobacteria lack the muropeptide transporter AmpG, despite having other key PG recycling enzymes. Here, we show that an alternative transporter, YejBEF-YepA, takes over this role in the Rhizobiales phytopathogen Agrobacterium tumefaciens. Muropeptide import by YejBEF-YepA governs expression of the β-lactamase AmpC in A. tumefaciens, contributing to β-lactam resistance. However, we show that the absence of YejBEF-YepA causes severe cell wall defects that go far beyond lowered AmpC activity. Thus, contrary to previously established Gram-negative models, PG recycling is vital for cell wall integrity in A. tumefaciens. YepA is widespread in the Rhizobiales and Rhodobacterales, suggesting that YejBEF-YepA-mediated PG recycling could represent an important but overlooked aspect of cell wall biology in these bacteria.

Evolutionary aspect of Miltefosine transporter proteins in Leishmania major

Transporter proteins, P-glycoprotein (P-gp) and P4ATPase-CDC50, are responsible for the transport of Miltefosine drug across cell membrane of a protozoan parasite Leishmania major. Mutations or change in activity of these proteins may lead to emergence of resistance in the parasite. Owing to the structural and functional importance of these transporter proteins, we have tried to decipher the evolutionary divergence of these Miltefosine transporter proteins across different forms of life including Protists, Fungi, Plants and Animals. We retrieved 96, 207, and 189 sequences of P-gp, P4ATPase and CDC50 proteins respectively, across diverse variety of organisms for the conserved analysis. Phylogenetic trees were constructed for these three transporter proteins based on Bayesian posterior probability inference. The evolutionary analysis concluded that these proteins remain highly conserved throughout the species diversity but still substantial differences in the proteins for host (Homo sapiens) and parasite (L. major) were observed which have led in targeting these Miltefosine transporter proteins in a parasite specific manner. The functional and structural components observed in terms of pattern resulting from the variability in the phylogenetic tree are outlined.

Identification of enzymes that have helminth-specific active sites and are required for Rhodoquinone-dependent metabolism as targets for new anthelmintics

Soil transmitted helminths (STHs) are major human pathogens that infect over a billion people. Resistance to current anthelmintics is rising and new drugs are needed. Here we combine multiple approaches to find druggable targets in the anaerobic metabolic pathways STHs need to survive in their mammalian host. These require rhodoquinone (RQ), an electron carrier used by STHs and not their hosts. We identified 25 genes predicted to act in RQ-dependent metabolism including sensing hypoxia and RQ synthesis and found 9 are required. Since all 9 have mammalian orthologues, we used comparative genomics and structural modeling to identify those with active sites that differ between host and parasite. Together, we found 4 genes that are required for RQ-dependent metabolism and have different active sites. Finding these high confidence targets can open up in silico screens to identify species selective inhibitors of these enzymes as new anthelmintics.

An Evolutionary Perspective of the Lipocalin Protein Family

The protein family of Lipocalins is ubiquitously present throughout the tree of life, with the exception of the phylum Archaea. Phylogenetic relationships of chordate Lipocalins have been proposed in the past based on protein sequence similarities, but their highly divergent primary structures and a shortage of experimental annotations in genome projects have precluded a well-supported hypothesis for their evolution. In this work we propose a novel topology for the phylogenetic tree of chordate Lipocalins, inferred from multiple amino acid sequence alignments. Sixteen jawed vertebrates with fair coverage by genomic sequencing were compared. The selected species span an evolutionary range of ∼400 million years, allowing for a balanced representation of all major vertebrate clades. A consensus phylogenetic tree is proposed following a comparison of sequence-based maximum-likelihood trees and protein structure dendrograms. This new phylogeny suggests an APOD-like common ancestor in early chordates, which gave rise, via whole-genome or tandem duplications, to the six Lipocalins currently present in fish (APOD, RBP4, PTGDS, AMBP, C8G, and APOM). Further gene duplications of APOM and PTGDS resulted in the altogether 15 Lipocalins found in contemporary mammals. Insights into the functional impact of relevant amino acid residues in early diverging Lipocalins are also discussed. These results should foster the experimental exploration of novel functions alongside the identification of new members of the Lipocalin family.

Structure, Function, and Thermal Adaptation of the Biotin Carboxylase Domain Dimer from Hydrogenobacter thermophilus 2-Oxoglutarate Carboxylase

2-Oxoglutarate carboxylase (OGC), a unique member of the biotin-dependent carboxylase family from the order Aquificales, captures dissolved CO₂ via the reductive tricarboxylic acid (rTCA) cycle. Structure and function studies of OGC may facilitate adaptation of the rTCA cycle to increase the level of carbon fixation for biofuel production. Here we compare the biotin carboxylase (BC) domain of Hydrogenobacter thermophilus OGC with the well-studied mesophilic homologues to identify features that may contribute to thermal stability and activity. We report three OGC BC X-ray structures, each bound to bicarbonate, ADP, or ADP-Mg²⁺, and propose that substrate binding at high temperatures is facilitated by interactions that stabilize the flexible subdomain B in a partially closed conformation. Kinetic measurements with varying ATP and biotin concentrations distinguish two temperature-dependent steps, consistent with biotin's rate-limiting role in organizing the active site. Transition state thermodynamic values derived from the Eyring equation indicate a larger positive ΔH^⧧ and a less negative ΔS^⧧ compared to those of a previously reported mesophilic homologue. These thermodynamic values are explained by partially rate limiting product release. Phylogenetic analysis of BC domains suggests that OGC diverged prior to Aquificales evolution. The phylogenetic tree identifies mis-annotations of the Aquificales BC sequences, including the Aquifex aeolicus pyruvate carboxylase structure. Notably, our structural data reveal that the OGC BC dimer comprises a "wet" dimerization interface that is dominated by hydrophilic interactions and structural water molecules common to all BC domains and likely facilitates the conformational changes associated with the catalytic cycle. Mutations in the dimerization domain demonstrate that dimerization contributes to thermal stability.

Molecular cloning of polymeric immunoglobulin receptor-like (pIgRL) in flounder (Paralichthys olivaceus) and its expression in response to immunization with inactivated Vibrio anguillarum

In the present work, the polymeric immunoglobulin receptor-like (pIgRL) from flounder (Paralichthys olivaceus) was firstly cloned and identified. The full length cDNA of flounder pIgRL was of 1393 bp including an open reading frame of 1053 bp, and the deduced pIgRL sequence encoded 350 amino acids, with a predicted molecular mass of 39 kDa. There were two immunoglobulin-like domains in flounder pIgRL. In healthy flounder, the transcriptional level of pIgRL was detected in different tissues by real-time PCR, showing the highest level in the skin and gills, and higher levels in the spleen and hindgut. After flounders were vaccinated with inactivated Vibrio anguillarum via intraperitoneal injection and immersion, the pIgRL mRNA level increased firstly and then declined in all tested tissues during 48 h, and the maximum expression levels in the gills, skin, spleen and hindgut in immersion group, or in the spleen, head kidney, skin and gills in injection group, were higher than in other tested tissues. In addition, recombinant protein of the extracellular region of flounder pIgRL was expressed in Escherichia coli BL21 (DE3), and rabbit anti-pIgRL polyclonal antibodies were prepared, which specifically reacted with the recombinant pIgRL, and a 39 kDa protein confirmed as natural pIgRL by liquid chromatography-mass spectrometry in skin mucus of flounder. Co-immunoprecipitation assay and western-blotting demonstrated that the pIgRL, together with IgM, could be immunoprecipitated by anti-pIgRL antibody in gut, skin and gill mucus of flounder, suggesting the existence of pIgRL-IgM complexes. These results indicated that the flounder pIgRL was probably involved in the mucosal IgM transportation and played important roles in mucosal immunity.

In silico identification and experimental validation of amino acid motifs required for the Rho-of-plants GTPase-mediated activation of receptor-like cytoplasmic kinases

Several amino acid motifs required for Rop-dependent activity were found to form a common surface on RLCKVI_A kinases. This indicates a unique mechanism for Rho-type GTPase-mediated kinase activation in plants. Rho-of-plants (Rop) G-proteins are implicated in the regulation of various cellular processes, including cell growth, cell polarity, hormonal and pathogen responses. Our knowledge about the signalling pathways downstream of Rops is continuously increasing. However, there are still substantial gaps in this knowledge. One reason for this is that these pathways are considerably different from those described for yeast and/or animal Rho-type GTPases. Among others, plants lack all Rho/Rac/Cdc42-activated kinase families. Only a small group of plant-specific receptor-like cytoplasmic kinases (RLCK VI_A) has been shown to exhibit Rop-binding-dependent in vitro activity. These kinases do not carry any known GTPase-binding motifs. Based on the sequence comparison of the Rop-activated RLCK VI_A and the closely related but constitutively active RLCK VI_B kinases, several distinguishing amino acid residues/motifs were identified. All but one of these were found to be required for the Rop-mediated regulation of the in vitro activity of two RLCK VI_A kinases. Structural modelling indicated that these motifs might form a common Rop-binding surface. Based on in silico data mining, kinases that have the identified Rop-binding motifs are present in Embryophyta but not in unicellular green algae. It can, therefore, be supposed that Rops recruited these plant-specific kinases for signalling at an early stage of land plant evolution.

Effect of the sequence data deluge on the performance of methods for detecting protein functional residues

BACKGROUND

The exponential accumulation of new sequences in public databases is expected to improve the performance of all the approaches for predicting protein structural and functional features. Nevertheless, this was never assessed or quantified for some widely used methodologies, such as those aimed at detecting functional sites and functional subfamilies in protein multiple sequence alignments. Using raw protein sequences as only input, these approaches can detect fully conserved positions, as well as those with a family-dependent conservation pattern. Both types of residues are routinely used as predictors of functional sites and, consequently, understanding how the sequence content of the databases affects them is relevant and timely.

RESULTS

In this work we evaluate how the growth and change with time in the content of sequence databases affect five sequence-based approaches for detecting functional sites and subfamilies. We do that by recreating historical versions of the multiple sequence alignments that would have been obtained in the past based on the database contents at different time points, covering a period of 20 years. Applying the methods to these historical alignments allows quantifying the temporal variation in their performance. Our results show that the number of families to which these methods can be applied sharply increases with time, while their ability to detect potentially functional residues remains almost constant.

CONCLUSIONS

These results are informative for the methods' developers and final users, and may have implications in the design of new sequencing initiatives.

A Bioinformatics Analysis Reveals a Group of MocR Bacterial Transcriptional Regulators Linked to a Family of Genes Coding for Membrane Proteins

The MocR bacterial transcriptional regulators are characterized by an N-terminal domain, 60 residues long on average, possessing the winged-helix-turn-helix (wHTH) architecture responsible for DNA recognition and binding, linked to a large C-terminal domain (350 residues on average) that is homologous to fold type-I pyridoxal 5′-phosphate (PLP) dependent enzymes like aspartate aminotransferase (AAT). These regulators are involved in the expression of genes taking part in several metabolic pathways directly or indirectly connected to PLP chemistry, many of which are still uncharacterized. A bioinformatics analysis is here reported that studied the features of a distinct group of MocR regulators predicted to be functionally linked to a family of homologous genes coding for integral membrane proteins of unknown function. This group occurs mainly in the Actinobacteria and Gammaproteobacteria phyla. An analysis of the multiple sequence alignments of their wHTH and AAT domains suggested the presence of specificity-determining positions (SDPs). Mapping of SDPs onto a homology model of the AAT domain hinted at possible structural/functional roles in effector recognition. Likewise, SDPs in wHTH domain suggested the basis of specificity of Transcription Factor Binding Site recognition. The results reported represent a framework for rational design of experiments and for bioinformatics analysis of other MocR subgroups.

Guanine nucleotide binding to the Bateman domain mediates the allosteric inhibition of eukaryotic IMP dehydrogenases

Inosine-5′-monophosphate dehydrogenase (IMPDH) plays key roles in purine nucleotide metabolism and cell proliferation. Although IMPDH is a widely studied therapeutic target, there is limited information about its physiological regulation. Using Ashbya gossypii as a model, we describe the molecular mechanism and the structural basis for the allosteric regulation of IMPDH by guanine nucleotides. We report that GTP and GDP bind to the regulatory Bateman domain, inducing octamers with compromised catalytic activity. Our data suggest that eukaryotic and prokaryotic IMPDHs might have developed different regulatory mechanisms, with GTP/GDP inhibiting only eukaryotic IMPDHs. Interestingly, mutations associated with human retinopathies map into the guanine nucleotide-binding sites including a previously undescribed non-canonical site and disrupt allosteric inhibition. Together, our results shed light on the mechanisms of the allosteric regulation of enzymes mediated by Bateman domains and provide a molecular basis for certain retinopathies, opening the door to new therapeutic approaches.

IMP dehydrogenase (IMPDH) plays essential roles in purine metabolism and cell proliferation. Here Buey et al. describe a guanine nucleotides regulated molecular mechanism for allosteric communication between the regulatory and catalytic domains of IMPDH.

Large-scale determination of previously unsolved protein structures using evolutionary information

The prediction of the structures of proteins without detectable sequence similarity to any protein of known structure remains an outstanding scientific challenge. Here we report significant progress in this area. We first describe de novo blind structure predictions of unprecendented accuracy we made for two proteins in large families in the recent CASP11 blind test of protein structure prediction methods by incorporating residue-residue co-evolution information in the Rosetta structure prediction program. We then describe the use of this method to generate structure models for 58 of the 121 large protein families in prokaryotes for which three-dimensional structures are not available. These models, which are posted online for public access, provide structural information for the over 400,000 proteins belonging to the 58 families and suggest hypotheses about mechanism for the subset for which the function is known, and hypotheses about function for the remainder.

Strigolactone promotes degradation of DWARF14, an α/β hydrolase essential for strigolactone signaling in Arabidopsis

Strigolactones (SLs) are phytohormones that play a central role in regulating shoot branching. SL perception and signaling involves the F-box protein MAX2 and the hydrolase DWARF14 (D14), proposed to act as an SL receptor. We used strong loss-of-function alleles of the Arabidopsis thaliana D14 gene to characterize D14 function from early axillary bud development through to lateral shoot outgrowth and demonstrated a role of this gene in the control of flowering time. Our data show that D14 distribution in vivo overlaps with that reported for MAX2 at both the tissue and subcellular levels, allowing physical interactions between these proteins. Our grafting studies indicate that neither D14 mRNA nor the protein move over a long range upwards in the plant. Like MAX2, D14 is required locally in the aerial part of the plant to suppress shoot branching. We also identified a mechanism of SL-induced, MAX2-dependent proteasome-mediated degradation of D14. This negative feedback loop would cause a substantial drop in SL perception, which would effectively limit SL signaling duration and intensity.

Molecular determinants of the ATP hydrolysis asymmetry of the CCT chaperonin complex

The eukaryotic cytosolic chaperonin CCT is a molecular machine involved in assisting the folding of proteins involved in important cellular processes. Like other chaperonins, CCT is formed by a double-ring structure but, unlike all of them, each ring is composed of eight different, albeit homologous subunits. This complexity has probably to do with the specificity in substrate interaction and with the mechanism of protein folding that takes place during the chaperonin functional cycle, but its detailed molecular basis remains unknown. We have analyzed the known proteomes in search of residues that are differentially conserved in the eight subunits, as predictors of functional specificity (specificity-determining positions; SDPs). We have found that most of these SDPs are located near the ATP binding site, and that they define four CCT clusters, corresponding to subunits CCT3, CCT6, CCT8 and CCT1/2/4/5/7. Our results point to a spatial organisation of the CCT subunits in two opposite areas of the ring and provide a molecular explanation for the previously described asymmetry in the hydrolysis of ATP.

Engineering the acyltransferase substrate specificity of assembly line polyketide synthases

Polyketide natural products act as a broad range of therapeutics, including antibiotics, immunosuppressants and anti-cancer agents. This therapeutic diversity stems from the structural diversity of these small molecules, many of which are produced in an assembly line manner by modular polyketide synthases. The acyltransferase (AT) domains of these megasynthases are responsible for selection and incorporation of simple monomeric building blocks, and are thus responsible for a large amount of the resulting polyketide structural diversity. The substrate specificity of these domains is often targeted for engineering in the generation of novel, therapeutically active natural products. This review outlines recent developments that can be used in the successful engineering of these domains, including AT sequence and structural data, mechanistic insights and the production of a diverse pool of extender units. It also provides an overview of previous AT domain engineering attempts, and concludes with proposed engineering approaches that take advantage of current knowledge. These approaches may lead to successful production of biologically active 'unnatural' natural products.

Getting personalized cancer genome analysis into the clinic: the challenges in bioinformatics

Progress in genomics has raised expectations in many fields, and particularly in personalized cancer research. The new technologies available make it possible to combine information about potential disease markers, altered function and accessible drug targets, which, coupled with pathological and medical information, will help produce more appropriate clinical decisions. The accessibility of such experimental techniques makes it all the more necessary to improve and adapt computational strategies to the new challenges. This review focuses on the critical issues associated with the standard pipeline, which includes: DNA sequencing analysis; analysis of mutations in coding regions; the study of genome rearrangements; extrapolating information on mutations to the functional and signaling level; and predicting the effects of therapies using mouse tumor models. We describe the possibilities, limitations and future challenges of current bioinformatics strategies for each of these issues. Furthermore, we emphasize the need for the collaboration between the bioinformaticians who implement the software and use the data resources, the computational biologists who develop the analytical methods, and the clinicians, the systems' end users and those ultimately responsible for taking medical decisions. Finally, the different steps in cancer genome analysis are illustrated through examples of applications in cancer genome analysis.