Search

Transthyretin was discovered in the 1940s, named after its ability to bind thyroid hormones and retinol. In the genomic era, transthyretins were found to be part of a larger family with homologs of no obvious function, then called transthyretin-related proteins. Thus, it was proposed that the transthyretin gene could be the result of gene duplication of an ancestral of this newly identified homolog, later found out to be an enzyme involved in uric acid degradation, then named HIUase (5-hydroxy-isourate hydrolase). Here, we sought to re-enact the evolutionary history of this protein family by reconstructing, from a phylogeny inferred from 123 vertebrate sequences, three ancestors corresponding to key moments in their evolution-before duplication; the common transthyretin ancestor after gene duplication and the common ancestor of Eutheria transthyretins. Experimental and computational characterization showed the reconstructed ancestor before duplication was unable to bind thyroxine and likely presented the modern HIUase reaction mechanism, while the substitutions after duplication prevented that activity and were enough to provide stable thyroxine binding, as confirmed by calorimetry and x-ray diffraction. The Eutheria transthyretin ancestor was less prone to characterization, but limited data suggested thyroxine binding as expected. Sequence/structure analysis suggests an early ability to bind the Retinol Binding Protein. We solved the X-ray structures from the two first ancestors, the first at 1.46 resolution, the second at 1.55 resolution with well-defined electron density for thyroxine, providing a useful tool for the understanding of structural adaptation from enzyme to hormone distributor.

Protein 3D structure can be remarkably robust to the accumulation of mutations during evolution. On the other hand, sometimes a single amino acid substitution can be sufficient to generate dramatic and completely unpredictable structural consequences. In an attempt to rationally alter the preferences for the metal ion at the active site of a member of the Iron/Manganese superoxide dismutase family, two examples of the latter phenomenon were identified. Site directed mutants of SOD from Trichoderma reesei were generated and studied crystallographically together with the wild type enzyme. Despite being chosen for their potential impact on the redox potential of the metal, two of the mutations (D150G and G73A) in fact resulted in significant alterations to the protein quaternary structure. The D150G mutant presented alternative inter-subunit contacts leading to a loss of symmetry of the wild type tetramer, whereas the G73A mutation transformed the tetramer into an octamer despite not participating directly in any of the inter-subunit interfaces. We conclude that there is considerable intrinsic plasticity in the Fe/MnSOD fold that can be unpredictably affected by single amino acid substitutions. In much the same way as phenotypic defects at the organism level can reveal much about normal function, so too can such mutations teach us much about the subtleties of protein structure.

The Golgi complex is a membranous organelle located in the heart of the eukaryotic secretory pathway. A subfamily of the Golgi matrix proteins, called GRASPs, are key players in the stress-induced unconventional secretion, the Golgi dynamics during mitosis/apoptosis, and Golgi ribbon formation. The Golgi ribbon is vertebrate-specific and correlates with the appearance of two GRASP paralogs (GRASP55/GRASP65) and two coiled-coil Golgins (GM130/Golgin45), which interact with each otherin vivo. Although essential for the Golgi ribbon formation and the increase in Golgi structural complexity, the molecular details leading to their appearance only in this subphylum are still unknown. Moreover, despite the new functionalities supported by the GRASP paralogy, little is known about the structural and evolutionary differences between these paralogues. In this context, we used ancestor sequence reconstruction and several biophysical/biochemical approaches to assess the evolution of the GRASP structure, flexibility, and how they started anchoring their Golgin partners. Our data showed that the Golgins appeared in evolution and were anchored by the single GRASP ancestor beforegoraspgene duplication and divergence in Metazoans. After thegoraspdivergence, variations inside the GRASP binding pocket determined which paralogue would recruit each Golgin partner (GRASP55 with Golgin45 and GRASP65 with GM130). These interactions are responsible for the protein’s specific Golgi locations and the appearance of the Golgi ribbon. We also suggest that the capacity of GRASPs to form supramolecular structures is a long-standing feature, which likely affects GRASP’s participation as a trigger of the stress-induced secretory pathway.

Coenzyme Q (CoQ) is an essential molecule that consists of a highly substituted benzene ring attached to a polyprenyl tail anchored in the inner mitochondrial membrane. CoQ transfers electrons from NADH dehydrogenase and succinate dehydrogenase complexes toward ubiquinol-cytochrome c reductase, and that allows aerobic growth of cells. In Saccharomyces cerevisiae, the synthesis of CoQ depends on fourteen proteins Coq1p-Co11p, Yah1p, Arh1p, and Hfd1p. Some of these proteins are components of CoQ synthome. Using ab initio molecular modeling and site-directed mutagenesis, we identified the functional residues of the O-methyltransferase Coq3p, which depends on S-adenosylmethionine for catalysis and is necessary for two O-methylation steps required for CoQ maturation. Conserved residues as well as those that coevolved in the protein structure were found to have important roles in respiratory growth, CoQ biosynthesis, and also in the stability of CoQ synthome proteins. Finally, a multiple sequence alignment showed that S. cerevisiae Coq3p has a 45 amino acid residues insertion that is poorly conserved or absent in oleaginous yeast, cells that can store up to 20% of their dry weight as lipids. These results point to the Coq3p structural determinants of its biological and catalytic function and could contribute to the development of lipid-producing yeast for biotechnology.

CEvADA is a database of amino acid coevolution networks aimed to detect specificity determinant and function related sites in protein families. The database was also designed to provide an easy access to protein coevolutionary constraints that can be incorporated in machine learning classification models, just as sequence annotation and structure prediction methods. The data can be accessed for the whole protein family and specific protein sequences. We also provide sequence search and a REST API for programmatic access in the database. The current version of the database contains data related to 6.301 conserved domains and 45 million protein sequences. CeVADA is free and can be accessed at http://bioinfo.icb.ufmg.br/cevada.

Transthyretin was discovered in the 1940s, named after its ability to bind thyroid hormones and retinol. In the genomic era, transthyretins were found to be part of a larger family with homologs of no obvious function, then called transthyretin-related proteins. Thus, it was proposed that the transthyretin gene could be the result of gene duplication of an ancestral of this newly identified homolog, later found out to be an enzyme involved in uric acid degradation, then named HIUase (5-hydroxy-isourate hydrolase). Here, we sought to re-enact the evolutionary history of this protein family by reconstructing, from a phylogeny inferred from 123 vertebrate sequences, three ancestors corresponding to key moments in their evolution-before duplication; the common transthyretin ancestor after gene duplication and the common ancestor of Eutheria transthyretins. Experimental and computational characterization showed the reconstructed ancestor before duplication was unable to bind thyroxine and likely presented the modern HIUase reaction mechanism, while the substitutions after duplication prevented that activity and were enough to provide stable thyroxine binding, as confirmed by calorimetry and x-ray diffraction. The Eutheria transthyretin ancestor was less prone to characterization, but limited data suggested thyroxine binding as expected. Sequence/structure analysis suggests an early ability to bind the Retinol Binding Protein. We solved the X-ray structures from the two first ancestors, the first at 1.46 resolution, the second at 1.55 resolution with well-defined electron density for thyroxine, providing a useful tool for the understanding of structural adaptation from enzyme to hormone distributor.

Summary CONAN is a web application developed to detect specificity determinants and function-related sites by amino acids co-variation networks analysis, emphasizing local coevolutionary constraints. The software allows the characterization of structurally and functionally relevant groups of residues and their relationship with subsets of sequences by automatic cross-referencing with GO terms, UniprotKb annotations and INTERPRO. Availability and implementation CONAN is free and open-source, being distributed in the terms of the GPLV3 license. The software is available as a web application and python script versions and can be accessed at http://bioinfo.icb.ufmg.br/conan. We also provide running instructions, the source code and a user guide.

Trypsin-like serine proteases are a group of homologous enzymes which exert multiple roles in both vertebrate and invertebrate organisms. Key properties of these enzymes include their activation from an inactive zymogen form to their active form by cleavage of residues in their N-terminus, the presence of a conserved catalytic triad of residues, and the existence of different patterns of substrate selectivity for residue cleavage between the various members of this protein family. In this article, we apply the decomposition of residue coevolution networks computational method to find sets of residues related to some of these key properties, especially to zymogen activation. Positive selection detection, normal modes analysis, and the calculation of thermal couplings between the bovine trypsinogen and bovine trypsin structures residues yielded further information for understanding the zymogen activation process and highlighted the importance of some of the coevolved set residues during these transitions.

Voltage-gated sodium channels play important roles in human physiology. However, their complexity hinders the understanding of their physiology and pathology at atomic level. We took advantage of the structural reports of similar channels obtained by cryo-EM (EeNav1.4, and NavPaS), and constructed models of human Nav1.4 channels at closed and open states. The open-state model is very similar to the recently published cryo-EM structure of hNav1.4. The comparison of both models shows shifts of the voltage sensors (VS) of DIII and DIV. The activated position of VS-DII in the closed model was demonstrated by Ts1 docking, thereby confirming the requirement that VS-DI, VS-DII and VS-DIII must be activated for the channel to open. The interactions observed with VS-DIII suggest a stepwise, yet fast, transition from resting to activated state. These models provide structural insights on the closed-open transition of the channel.

PFstats is a software developed for the extraction of useful information from protein multiple sequence alignments. By analyzing positional conservation and residue coevolution networks, the software allows the identification of structurally and functionally important residue groups and the discovery of probable functional subclasses. Furthermore, it contains tools for the identification of the possible biological significance of these findings. PFstats contains methods for maximizing the significance of alignments through filtering and weighting, residue conservation and coevolution analysis, automatic UniprotKb queries for residue-position annotation and many possible data visualization methods.

Saccharomyces cerevisiae mitoribosomes are specialized in the translation of a few number of highly hydrophobic membrane proteins, components of the oxidative phosphorylation system. Mitochondrial characteristics, such as the membrane system and its redox state driven mitoribosomes evolution through great diversion from their bacterial and cytosolic counterparts. Therefore, mitoribosome presents a considerable number of mitochondrial-specific proteins, as well as new protein extensions. In this work we characterize temperature sensitive mutants of the subunit bL34 present in the 54S large subunit. Although bL34 has bacterial homologs, in yeast it has a long 65 aminoacids mitochondrial N-terminal addressing sequence, here we demonstrate that it can be replaced by the mitochondrial addressing sequence of Neurospora crassa ATP9 gene. The bL34 temperature sensitive mutants present lowered translation of mitochondrial COX1 and COX3, which resulted in reduced cytochrome c oxidase activity and respiratory growth deficiency. The sedimentation properties of bL34 in sucrose gradients suggest that similarly to its bacterial homolog, bL34 is also a later participant in the process of mitoribosome biogenesis.

Flaviviruses are responsible for serious diseases such as dengue, yellow fever, and zika fever. Their genomes encode a polyprotein which, after cleavage, results in three structural and seven non-structural proteins. Homologous proteins can be studied by conservation and coevolution analysis as detected in multiple sequence alignments, usually reporting positions which are strictly necessary for the structure and/or function of all members in a protein family or which are involved in a specific sub-class feature requiring the coevolution of residue sets. This study provides a complete conservation and coevolution analysis on all flaviviruses non-structural proteins, with results mapped on all well-annotated available sequences. A literature review on the residues found in the analysis enabled us to compile available information on their roles and distribution among different flaviviruses. Also, we provide the mapping of conserved and coevolved residues for all sequences currently in SwissProt as a supplementary material, so that particularities in different viruses can be easily analyzed.

Pb27 antigen is an interesting alternative to immunological diagnosis of Paracoccidioidomycosis (PCM) and has demonstrated to be protective in experimental PCM. Its tertiary structure and possible function remained unknown till now. To study Pb27 at the atomic level, the recombinant protein was expressed in Escherichia coli BL21(DE3), purified, and its three-dimensional structure was solved by X-ray crystallography. Based on this structure, we performed a residue correlation analysis and in silico ligand search assays to address a possible biological function to Pb27. We identified Pb27 as a member of the extensive nucleotidyltransferase superfamily. The protein has an αβαβαβ topology with two domains (N- and C-terminal domains) and adopts a monomeric form as its biological unit in solution. Structural comparisons with similar members of the superfamily clearly indicate Pb27 C-terminal domain is singular and may play an important role in its biological function. Bioinformatics analysis suggested that Pb27 might bind to ATP and CTP. This suggestion is corroborated by the fact that a magnesium cation is coordinated by two aspartic acid residues present at the active site (between N- and C-terminal domains), as evidenced by X-ray diffraction data. Besides, NMR assays (1H-15N HSQC spectra) confirmed the binding of CTP to Pb27, demonstrating for the first time an interaction between a nucleotide and this protein. Moreover, we evaluated the reactivity of sera from patients with Paracoccidioides brasiliensis infection against the recombinant form of Pb27 and showed that it was recognized by sera from infected and treated patients. Predicted B and T cell epitopes were synthesized and further evaluated against sera of PCM patients, providing information of the most reactive peptides in Pb27 primary structure which interact with specific Pb27 antibodies.

Motivation Computational studies of molecular evolution are usually performed from a multiple alignment of homologous sequences, on which sequences resulting from a common ancestor are aligned so that equivalent residues are placed in the same position. Residues frequency patterns of a full alignment or from a subset of its sequences can be highly useful for suggesting positions under selection. Most methods mapping co-evolving or specificity determinant sites are focused on positions, however, they do not consider the case for residues that are specificity determinants in one subclass, but variable in others. In addition, many methods are impractical for very large alignments, such as those obtained from Pfam, or require a priori information of the subclasses to be analyzed. Results In this paper we apply the complex networks theory, widely used to analyze co-affiliation systems in the social and ecological contexts, to map groups of functional related residues. This methodology was initially evaluated in simulated environments and then applied to four different protein families datasets, in which several specificity determinant sets and functional motifs were successfully detected. Availability and implementation The algorithms and datasets used in the development of this project are available on http://www.biocomp.icb.ufmg.br/biocomp/software-and-databases/networkstats/. Supplementary information Supplementary data are available at Bioinformatics online.

Objective Thyroxine-binding globulin (TBG) is the major human thyroid hormone transport protein, encoded by the SERPINA7 gene (Xq22.2). We aim to investigate the molecular basis of partial TBG deficiency (TBG-PD) in a female, by evaluating the X-chromosome inactivation pattern as well as the mutant protein structural modeling. Design and methods Sequencing of the coding region of the SERPINA7 gene was performed in a female with a TBG-PD phenotype and her first-degree relatives. The proband presented with low serum levels of total T3 (TT3) and total T4 (TT4), serum TSH level of 5.4 μUI/mL (normal range, 0.35–5.5), and serum TBG level of 5.5 mg/L (normal range, 13.6–27.2). X-chromosome inactivation pattern was evaluated by methylation analysis of the androgen receptor gene (Xq11.2). Structural analysis of the SERPIN family was performed using Pymol and Areaimol, and PFSTATS for conservation analysis and family-wide investigation of equivalent positions in human homologs. Modeller was used for point mutation structural modeling. Results A novel missense SERPINA7 mutation (p.R35W; c.163C > T) was found in heterozygosity in the proband, and in hemizygosity in her affected siblings. The proband X-chromosome inactivation ratio was 20:80. The substitution of an arginine by a tryptophan is predicted to disrupt the protein surface and main electrostatic interactions. Tryptophans are extremely rare (0.1%) in this position. Conclusions We report a new SERPINA7 variant associated with TBG-PD in three siblings. We named this variant TBG-Brasilia. The X-chromosome inactivation pattern may have accounted for the rare phenotypic expression in a female. The hydrophobic nature of the mutant is predicted to create an apolar patch at the surface, which results in protein aggregation and/or misfolding, potentially responsible for thyroid hormone transport defect.

Phospholipases A2 (PLA2s) comprise a superfamily of glycerophospholipids hydrolyzing enzymes present in many organisms in nature, whose catalytic activity was majorly unveiled by analysis of snake venoms. The latter have pharmaceutical and biotechnological interests and can be divided into different functional sub-classes. Our goal was to identify important residues and their relation to the functional and class-specific characteristics in the PLA2s family with special emphasis on snake venom PLA2s (svPLA2s). We identified such residues by conservation analysis and decomposition of residue coevolution networks (DRCN), annotated the results based on the available literature on PLA2s, structural analysis and molecular dynamics simulations, and related the results to the phylogenetic distribution of these proteins. A filtered alignment of PLA2s revealed 14 highly conserved positions and 3 sets of coevolved residues, which were annotated according to their structural or functional role. These residues are mostly involved in ligand binding and catalysis, calcium-binding, the formation of disulfide bridges and a hydrophobic cluster close to the binding site. An independent validation of the inference of structure-function relationships from our co-evolution analysis on the svPLA2s family was obtained by the analysis of the pattern of selection acting on the Viperidae and Elapidae lineages. Additionally, a molecular dynamics simulation on the Lys49 PLA2 from Agkistrodon contortrix laticinctus was carried out to further investigate the correlation of the Lys49-Glu69 pair. Our results suggest this configuration can result in a novel conformation where the binding cavity collapses due to the approximation of two loops caused by a strong salt bridge between Glu69 and Arg34. Finally, phylogenetic analysis indicated a correlation between the presence of residues in the coevolved sets found in this analysis and the clade localization. The results provide a guide for important positions in the family of PLA2s, and potential new objects of investigation.

Bacteria are essential in arsenic cycling. However, few studies have addressed 16S rRNA and arsenic-related functional gene diversity in long-term arsenic-contaminated tropical sediment. Here, using culture-based, metagenomic and computational approaches, we describe the diversity of bacteria, genes and enzymes involved in AsIII and AsV transformation in freshwater sediment and in anaerobic AsIII- and AsV-enrichment cultures (ECs). The taxonomic profile reveals significant differences among the communities. Arcobacter, Dechloromonas, Sedimentibacter and Clostridium thermopalmarium were exclusively found in ECs, whereas Anaerobacillus was restricted to AsV-EC. Novel taxa that are both AsV-reducers and AsIII-oxidizers were identified: Dechloromonas, Acidovorax facilis, A. delafieldii, Aquabacterium, Shewanella, C. thermopalmarium and Macellibacteroides fermentans. Phylogenic discrepancies were revealed among the aioA, arsC and arrA genes and those of other species, indicating horizontal gene transfer. ArsC and AioA have sets of amino acids that can be used to assess their functional and structural integrity and familial subgroups. The positions required for AsV reduction are conserved, suggesting strong selective pressure for maintaining the functionality of ArsC. Altogether, these findings highlight the role of freshwater sediment bacteria in arsenic mobility, and the untapped diversity of dissimilatory arsenate-reducing and arsenate-resistant bacteria, which might contribute to arsenic toxicity in aquatic environments.

In Saccharomyces cerevisiae mitochondrial dysfunction induces retrograde signaling, a pathway of communication from mitochondria to the nucleus that promotes a metabolic remodeling to ensure sufficient biosynthetic precursors for replication. Rtg2p is a positive modulator of this pathway that is also required for cellular longevity. This protein belongs to the ASKHA superfamily, and contains a putative N-terminal ATP-binding domain, but there is no detailed structural and functional map of the residues in this domain that accounts for their contribution to retrograde signaling and aging. Here we use Decomposition of Residue Correlation Networks and site-directed mutagenesis to identify Rtg2p structural determinants of retrograde signaling and longevity. We found that most of the residues involved in retrograde signaling surround the ATP-binding loops, and that Rtg2p N-terminus is divided in three regions whose mutants have different aging phenotypes. We also identified E137, D158 and S163 as possible residues involved in stabilization of ATP at the active site. The mutants shown here may be used to map other Rtg2p activities that crosstalk to other pathways of the cell related to genomic stability and aging.

Wastewater treatment plants (WWTPs) harbor bacteria and antimicrobial resistance genes, favoring gene exchange events and resistance dissemination. Here, a culture-based and metagenomic survey of qnrA, qnrB, qnrS, and aac(6′)-Ib genes from raw sewage (RS) and activated sludge (AS) of a full-scale municipal WWTP was performed. A total of 96 bacterial isolates were recovered from nalidixic acid-enrichment cultures. Bacteria harboring the aac(6′)-Ib gene predominated in RS, whereas qnrS-positive isolates were specific to AS. Novel qnrS- and aac(6′)-Ib-cr positive species were identified: Morganella morganii, Providencia rettgeri, and Pseudomonas guangdongensis (qnrS), and Alcaligenes faecalis and P. rettgeri (aac(6′)-Ib-cr). Analysis of qnrS and aac(6′)-Ib sequences from isolates and clone libraries suggested that the diversity of qnrS is wider than that of aac(6′)-Ib. A large number of amino acid mutations were observed in the QnrS and AAC(6′)-Ib proteins at previously undetected positions, whose structural implications are not clear. An accumulation of mutations at the C72, Q73, L74, A75 and M76 positions of QnrS, and D181 of AAC(6′)-Ib might be important for resistance. These findings add significant information on bacteria harboring qnrS and aac(6′)-Ib genes, and the presence of novel mutations that may eventually emerge in clinical isolates.

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 Å resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting β-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

Human transthyretin (TTR) is a homotetrameric protein involved in several amyloidoses. Zn2+ enhances TTR aggregation in vitro, and is a component of ex vivo TTR amyloid fibrils. We report the first crystal structure of human TTR in complex with Zn2+ at pH 4.6–7.5. All four structures reveal three tetra-coordinated Zn2+-binding sites (ZBS 1–3) per monomer, plus a fourth site (ZBS 4) involving amino acid residues from a symmetry-related tetramer that is not visible in solution by NMR. Zn2+ binding perturbs loop E-α-helix-loop F, the region involved in holo-retinol-binding protein (holo-RBP) recognition, mainly at acidic pH; TTR affinity for holo-RBP decreases ∼5-fold in the presence of Zn2+. Interestingly, this same region is disrupted in the crystal structure of the amyloidogenic intermediate of TTR formed at acidic pH in the absence of Zn2+. HNCO and HNCA experiments performed in solution at pH 7.5 revealed that upon Zn2+ binding, although the α-helix persists, there are perturbations in the resonances of the residues that flank this region, suggesting an increase in structural flexibility. While stability of the monomer of TTR decreases in the presence of Zn2+, which is consistent with the tertiary structural perturbation provoked by Zn2+ binding, tetramer stability is only marginally affected by Zn2+. These data highlight structural and functional roles of Zn2+ in TTR-related amyloidoses, as well as in holo-RBP recognition and vitamin A homeostasis.

Transthyretin (TTR) is a tetrameric β-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation ( Green et al., 2005 ) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

Superoxide dismutases (SODs) are a crucial class of enzymes in the combat against intracellular free radical damage. They eliminate superoxide radicals by converting them into hydrogen peroxide and oxygen. In spite of their very different life cycles and infection strategies, the human parasites Plasmodium falciparum, Trypanosoma cruzi and Trypanosoma brucei are known to be sensitive to oxidative stress. Thus the parasite Fe-SODs have become attractive targets for novel drug development. Here we report the crystal structures of FeSODs from the trypanosomes T. brucei at 2.0 A and T. cruzi at 1.9 A resolution, and that from P. falciparum at a higher resolution (2.0 A) to that previously reported. The homodimeric enzymes are compared to the related human MnSOD with particular attention to structural aspects which are relevant for drug design. Although the structures possess a very similar overall fold, differences between the enzymes at the entrance to the channel which leads to the active site could be identified. These lead to a slightly broader and more positively charged cavity in the parasite enzymes. Furthermore, a statistical coupling analysis (SCA) for the whole Fe/MnSOD family reveals different patterns of residue coupling for Mn and Fe SODs, as well as for the dimeric and tetrameric states. In both cases, the statistically coupled residues lie adjacent to the conserved core surrounding the metal center and may be expected to be responsible for its fine tuning, leading to metal ion specificity. Proteins 2009. © 2009 Wiley-Liss, Inc.

The orphan receptor nerve growth factor-induced B (NGFI-B) is a member of the nuclear receptor's subfamily 4A (Nr4a). NGFI-B was shown to be capable of binding both as a monomer to an extended half-site containing a single AAAGGTCA motif and also as a homodimer to a widely separated everted repeat, as opposed to a large number of nuclear receptors that recognize and bind specific DNA sequences predominantly as homo- and/or heterodimers. To unveil the structural organization of NGFI-B in solution, we determined the quaternary structure of the NGFI-B LBD by a combination of ab initio procedures from small-angle X-ray scattering (SAXS) data and hydrogen–deuterium exchange followed by mass spectrometry. Here we report that the protein forms dimers in solution with a radius of gyration of 2.9 nm and maximum dimension of 9.0 nm. We also show that the NGFI-B LBD dimer is V-shaped, with the opening angle significantly larger than that of classical dimer's exemplified by estrogen receptor (ER) or retinoid X receptor (RXR). Surprisingly, NGFI-B dimers formation does not occur via the classical nuclear receptor dimerization interface exemplified by ER and RXR, but instead, involves an extended surface area composed of the loop between helices 3 and 4 and C-terminal fraction of the helix 3. Remarkably, the NGFI-B dimer interface is similar to the dimerization interface earlier revealed for glucocorticoid nuclear receptor (GR), which might be relevant to the recognition of cognate DNA response elements by NGFI-B and to antagonism of NGFI-B–dependent transcription exercised by GR in cells.

A Análise de Acoplamentos Estatísticos é uma técnica computacional capaz de identificar resíduos importantes para a estrutura e função de proteínas em uma família por meio da quantificação de conservação posicional, correlação entre posições e identificação de grupos de resíduos correlacionados entre si. Neste trabalho, a análise de acoplamentos estatísticos foi implementada e aplicada ao estudo das proteínas tirosina fosfatases. Em conjunto com as proteínas tirosina quinases (PTKs), que adicionam um grupo fosforil a um resíduo de tirosina em uma proteína, as proteínas tirosina fosfatases (PTPs), que o removem, são responsáveis por diversos processos de sinalização celular. Elas são um caso de evolução convergente, onde um subgrupo (as proteínas tirosina fosfatases de baixo peso molecular) não apresenta homologia às chamadas PTPs \"clássicas\", capazes de defosforilar apenas resíduos de tirosina, e às fosfatases de especifidicade dupla, capazes de defosforilar também resíduos de serina e treonina, além de substratos não-protéicos. Em comum, as três subfamílias apresentam apenas o motivo CX5R, característico para todas as PTPs. Através do estudo das três subfamílias utilizando a análise de acoplamentos estatísticos, foi possível obter uma descrição detalhada de suas características conservadas e correlacionadas, relacionando-as ao conhecimento acumulado sobre proteínas tirosina fosfatases e a questões em aberto como a regulação por dimerização, a especificidade e mutações relacionadas a patologias. Foi possível também apresentar um método capaz de distinguir proteínas tirosina fosfatases de baixo peso molecular das arsenato redutases, derivadas das primeiras por evolução divergente. Adicionalmente, a técnica foi aplicada ao estudo das hexoquinases, às superóxido dismutases e às peroxidases. A tese descreve também estudos desenvolvidos pelo autor na área de cristalografia de proteínas a determinação das estruturas da Transtirretina humana em complexo com genisteína, da holo-Hexoquinase PI de S. cerevisae, do complexo IL-22/IL-22R1 e da Laminarinase de R. marinus. The statistical coupling analysis is a computational technique which can identify important residues for the structure and function of proteins in a family by quantifying positional conservation, correlation between positions and identifying groups of self-correlating residues. Its implementation in this research group was applied to the study of the protein tyrosine phosphatases. Together with the protein tyrosine kinases (PTKs), which add a phosphoryl group to a tyrosine residue in proteins, the protein tyrosine phosphatases (PTPs), which remove it, are responsible for a variety of cell signaling processes. They are a case of convergent evolution, since one subgroup (the low molecular weight protein tyrosine phosphatases) are not homologous to the classical phosphatases, which can only dephosphorilate tyrosine residues, and the dual-specificity phosphatases, which can also dephosphorilate serine and threonine residues, and also non-proteinaceous substrates. All three sub-families have, in common, the CX5R motif, a characteristic of all PTPs. By applying the statistical coupling analysis to the study of the three sub-families, it was possible to obtain a detailed depiction of their conserved and correlated characteristics, relating them to the accumulated knowledge on protein tyrosine phosphatases and open questions such as protein regulation by dimerization, specificity and disease-related mutations. It was also possible to present a method to distinguish between low molecular weight phosphatases and arsenate reductases, which are derived by the former by divergent evolution. In addition, the technique was applied to the study of hexokinases, superoxide dismutases and peroxidases. The thesis also describe studies developed by the author in the field of protein crystallography the structure determination of human transthyretin in complex with genistein, holo-hexokinase PI from S. cerevisae, the IL-22/IL-22R1 complex and the laminarinase from R. marinus.

Coenzyme Q (Q) is an isoprenylated benzoquinone electron carrier required for electronic transport in the mitochondrial respiratory chain, shuttling electrons from complexes I and II to complex III. Q synthesis requires proteins termed Coq (Coq1-Coq11). Coq7p is part of the multimeric complex involved in Q synthesis catalyzing the hydroxylation of demethoxy-Q6 (DMQ6), the last monooxygenase step in Q synthesis with a catalytic center containing a carboxylate-bridged di-iron at the active site of the enzyme. Here we indicate a group of Coq7p residues that modulate protein activity: D53, R57, V111 and S114. R57, V111 and S114 are very conserved residues; V111 and S114 are present in separated communities of amino acid correlation analysis. The coq7 double mutant V111G/S114A and S114E show respiratory deficiency at non permissive temperature, DMQ6 accumulation and lower content of Q6. Therefore we conclude that phosphomimetic S114E inhibit Coq7p activity, and propose that S114 phosphorylation is required to move a non-structured loop of 25 amino acids between helix 2 and 3, and that affects the di-iron coordination in Coq7p catalytic center.

The thyroid hormone receptor (TR) D-domain links the ligand-binding domain (LBD, EF-domain) to the DNA-binding domain (DBD, C-domain), but its structure, and even its existence as a functional unit, are controversial. The D domain is poorly conserved throughout the nuclear receptor family and was originally proposed to comprise an unfolded hinge that facilitates rotation between the LBD and the DBD. Previous TR LBD structures, however, have indicated that the true unstructured region is three to six amino acid residues long and that the D-domain N terminus folds into a short amphipathic alpha-helix (H0) contiguous with the DBD and that the C terminus of the D-domain comprises H1 and H2 of the LBD. Here, we solve structures of TR-LBDs in different crystal forms and show that the N terminus of the TRalpha D-domain can adopt two structures; it can either fold into an amphipathic helix that resembles TRbeta H0 or form an unstructured loop. H0 formation requires contacts with the AF-2 coactivator-binding groove of the neighboring TR LBD, which binds H0 sequences that resemble coactivator LXXLL motifs. Structural analysis of a liganded TR LBD with small angle X-ray scattering (SAXS) suggests that AF-2/H0 interactions mediate dimerization of this protein in solution. We propose that the TR D-domain has the potential to form functionally important extensions of the DBD and LBD or unfold to permit TRs to adapt to different DNA response elements. We also show that mutations of the D domain LXXLL-like motif indeed selectively inhibit TR interactions with an inverted palindromic response element (F2) in vitro and TR activity at this response element in cell-based transfection experiments.

We present a program that uses an optimization algorithm to fit rocking curves of ion-implanted semiconductors. This is an inverse problem that cannot be solved by simple methods. However, using recursion formulae for rocking curve calculations and a model of ion distribution after implantation, it is possible to fit experimental data with a general-purpose optimization method. In our case, we use a modified version of the genetic algorithm, which has been shown to be a good technique for this problem. The program also calculates rocking curves for a given ion profile, such as those generated by ion implantation simulation programs.

Made available in DSpace on 2013-08-28T14:09:29Z (GMT). No. of bitstreams: 1 WOS000298666200001.pdf: 402150 bytes, checksum: 3723ce8567931f389f45074a49611911 (MD5) Made available in DSpace on 2013-09-30T18:37:20Z (GMT). No. of bitstreams: 2 WOS000298666200001.pdf: 402150 bytes, checksum: 3723ce8567931f389f45074a49611911 (MD5) WOS000298666200001.pdf.txt: 56655 bytes, checksum: 10d0e581654ecca834117e332f3f69ca (MD5) Previous issue date: 2011-12-20 Submitted by Vitor Silverio Rodrigues (vitorsrodrigues@reitoria.unesp.br) on 2014-05-20T13:49:40Z No. of bitstreams: 2 WOS000298666200001.pdf: 402150 bytes, checksum: 3723ce8567931f389f45074a49611911 (MD5) WOS000298666200001.pdf.txt: 56655 bytes, checksum: 10d0e581654ecca834117e332f3f69ca (MD5) Made available in DSpace on 2014-05-20T13:49:40Z (GMT). No. of bitstreams: 2 WOS000298666200001.pdf: 402150 bytes, checksum: 3723ce8567931f389f45074a49611911 (MD5) WOS000298666200001.pdf.txt: 56655 bytes, checksum: 10d0e581654ecca834117e332f3f69ca (MD5) Previous issue date: 2011-12-20 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Correlated mutation analysis has a long history of interesting applications, mostly in the detection of contact pairs in protein structures. Based on previous observations that, if properly assessed, amino acid correlation data can also provide insights about functional sub-classes in a protein family, we provide a complete framework devoted to this purpose. An amino acid specific correlation measure is proposed, which can be used to build networks summarizing all correlation and anti-correlation patterns in a protein family. These networks can be submitted to community structure detection algorithms, resulting in subsets of correlated amino acids which can be further assessed by specific parameters and procedures that provide insight into the relationship between different communities, the individual importance of community members and the adherence of a given amino acid sequence to a given community. By applying this framework to three protein families with contrasting characteristics (the Fe/Mn-superoxide dismutases, the peroxidase-catalase family and the C-type lysozyme/alpha-lactalbumin family), we show how our method and the proposed parameters and procedures are related to biological characteristics observed in these protein families, highlighting their potential use in protein characterization and gene annotation. Univ Fed Minas Gerais, Inst Ciencias Biol, Dept Bioquim & Imunol, Belo Horizonte, MG, Brazil Univ Estadual Paulista, Dept Fis & Biofis, Botucatu, SP, Brazil Univ São Paulo, Inst Fis Sao Carlos, Dept Fis & Informat, Sao Carlos, SP, Brazil Univ Estadual Paulista, Dept Fis & Biofis, Botucatu, SP, Brazil FAPESP: 08/58734-1 FAPESP: 98/14138-2

Diffraction data to 3.2 Å from crystals of the 3.6 MDa erythrocruorin from a Brazilian earthworm represent the highest resolution reported to date for similar complexes. An unambiguous molecular replacement solution shows the particle to belong to the type I class., Glossoscolex paulistus is a free-living earthworm encountered in south-east Brazil. Its oxygen transport requirements are undertaken by a giant extracellular haemoglobin, or erythrocruorin (HbGp), which has an approximate molecular mass of 3.6 MDa and, by analogy with its homologue from Lumbricus terrestris (HbLt), is believed to be composed of a total of 180 polypeptide chains. In the present work the full 3.6 MDa particle in its cyanomet state was purified and crystallized using sodium citrate or PEG8000 as precipitant. The crystals contain one-quarter of the full particle in the asymmetric unit of the I222 cell and have parameters of a = 270.8 Å, b = 320.3 Å and c = 332.4 Å. Diffraction data were collected to 3.15 Å using synchrotron radiation on beamline X29A at the Brookhaven National Laboratory and represent the highest resolution data described to date for similar erythrocruorins. The structure was solved by molecular replacement using a search model corresponding to one-twelfth of its homologue from HbLt. This revealed that HbGp belongs to the type I class of erythrocruorins and provided an interpretable initial electron density map in which many features including the haem groups and disulfide bonds could be identified.

Coq10p is a protein required for coenzyme Q function, but its specific role is still unknown. It is a member of the START domain superfamily that contains a hydrophobic tunnel implicated in the binding of lipophilic molecules. We used site-directed mutagenesis, statistical coupling analysis and molecular modeling to probe structural determinants in the Coq10p putative tunnel. Four point mutations were generated (coq10-K50E, coq10-L96S, coq10-E105K and coq10-K162D) and their biochemical properties analysed, as well as structural consequences. Our results show that all mutations impaired Coq10p function and together with molecular modeling indicate an important role for the Coq10p putative tunnel.

Royal palm tree peroxidase (RPTP) is a very stable enzyme in regards to acidity, temperature, H(2)O(2), and organic solvents. Thus, RPTP is a promising candidate for developing H(2)O(2)-sensitive biosensors for diverse applications in industry and analytical chemistry. RPTP belongs to the family of class III secretory plant peroxidases, which include horseradish peroxidase isozyme C, soybean and peanut peroxidases. Here we report the X-ray structure of native RPTP isolated from royal palm tree (Roystonea regia) refined to a resolution of 1.85A. RPTP has the same overall folding pattern of the plant peroxidase superfamily, and it contains one heme group and two calcium-binding sites in similar locations. The three-dimensional structure of RPTP was solved for a hydroperoxide complex state, and it revealed a bound 2-(N-morpholino) ethanesulfonic acid molecule (MES) positioned at a putative substrate-binding secondary site. Nine N-glycosylation sites are clearly defined in the RPTP electron-density maps, revealing for the first time conformations of the glycan chains of this highly glycosylated enzyme. Furthermore, statistical coupling analysis (SCA) of the plant peroxidase superfamily was performed. This sequence-based method identified a set of evolutionarily conserved sites that mapped to regions surrounding the heme prosthetic group. The SCA matrix also predicted a set of energetically coupled residues that are involved in the maintenance of the structural folding of plant peroxidases. The combination of crystallographic data and SCA analysis provides information about the key structural elements that could contribute to explaining the unique stability of RPTP.

Interleukin-22 (IL-22) plays an important role in the regulation of immune and inflammatory responses in mammals. The IL-22 binding protein (IL-22BP), a soluble receptor that specifically binds IL-22, prevents the IL-22/interleukin-22 receptor 1 (IL-22R1)/interleukin-10 receptor 2 (IL-10R2) complex assembly and blocks IL-22 biological activity. Here we present the crystal structure of the IL-22/IL-22BP complex at 2.75Å resolution. The structure reveals IL-22BP residues critical for IL-22 binding, which were confirmed by site-directed mutagenesis and functional studies. Comparison of IL-22/IL-22BP and IL-22/IL-22R1 crystal structures shows that both receptors display an overlapping IL-22 binding surface, which is consistent with the inhibitory role played by IL-22 binding protein.Structured summaryMINT-7010533: IL-22 BP (uniprotkb:Q969J5) and IL-22 (uniprotkb:Q9GZX6) bind (MI:0407) by X-ray crystallography (MI:0114)

Interleukin-22 (IL-22) is a member of the interleukin-10 cytokine family, which is involved in anti-microbial defenses, tissue damage protection and repair, and acute phase responses. Its signaling mechanism involves the sequential binding of IL-22 to interleukin-22 receptor 1 (IL-22R1), and of this dimer to interleukin-10 receptor 2 (IL-10R2) extracellular domain. We report a 1.9Å crystal structure of the IL-22/IL-22R1 complex, revealing crucial interacting residues at the IL-22/IL-22R1 interface. Functional importance of key residues was confirmed by site-directed mutagenesis and functional studies. Based on the X-ray structure of the binary complex, we discuss a molecular basis of the IL-22/IL-22R1 recognition by IL-10R2.Structured summaryMINT-6693956:Il22 (uniprotkb:Q9GZX6) and IL22R1 (uniprotkb:Q8N6P7) bind (MI:0407) by X-ray crystallography (MI:0114)