Lucas Bleicher

Using Amino Acid Correlation and Community Detection Algorithms to Identify Functional Determinants in Protein Families

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/α-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.

Novel arsenic-transforming bacteria and the diversity of their arsenic-related genes and enzymes arising from arsenic-polluted freshwater sediment

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.

Sequence, structure and function relationships in flaviviruses as assessed by evolutive aspects of its conserved non-structural protein domains

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.

Coq7p relevant residues for protein activity and stability

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.

Her2p molecular modeling, mutant analysis and intramitochondrial localization.

Bacterial GatCAB amidotransferases are responsible for the transamidation of mischarged glutamyl-tRNA(Gln) into glutaminyl-tRNA(Gln). Mitochondria matrix also has a multienzymatic complex necessary for the transamidation of glutamyl-tRNA(Gln). Gtf1p, Her2p and Pet112p are the constituents of mitochondrial GatFAB amidotransferase complex. Her2p is subunit A of GatFAB complex, while Gtf1p is subunit F, a connector protein between Pet112p (subunit B) and Her2p. Here we evaluate through molecular modeling and amino acid correlation analysis the HER2 protein family. Localization studies indicated that Her2p is predominantly localized in the mitochondrial outer membrane, but it is also located in the mitochondrial matrix where together with Pet112p and Gtf1p constitutes the GatFAB complex. Finally, HER2 random mutagenesis unveiled important residues that provide thermo stability for the complex and are differently suppressed by overexpression of GTF1 or PET112. For instance, her2/ts11 mutant showed its fermentative growth impaired, and poorly rescued by GTF1 indicating that Her2p unknown function in the mitochondria outer membrane affects cell viability.

Residue correlation networks in nuclear receptors reflect functional specialization and the formation of the nematode-specific P-box

BackgroundNuclear receptors (NRs) are transcription factors which bind small hormones, whose evolutionary history and the presence of different functional surfaces makes them an interesting target for a correlation based analysis.ResultsCorrelation analysis of ligand binding domains shows that correlated residue subsets arise from the differences between functional sites in different nuclear receptor subfamilies. For the DNA binding domain, particularly, the analysis shows that the main source of correlation comes from residues that regulate hormone response element specificity, and one of the conserved residue sub-sets arises due to the presence of an unusual sequence for the DNA binding motif known as P-box in nematodes, suggesting the existence of different DBD-DNA specificities in nuclear receptors.ConclusionsWe conclude that DNA specificity and functional surface specialization has independently driven nuclear receptor evolution, and suggest possible binding modes for the class of divergent nuclear receptors in nematodes.

Conservation analysis and decomposition of residue correlation networks in the phospholipase A2 superfamily (PLA2s): Insights into the structure-function relationships of snake venom toxins

&NA; 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. HighlightsOur results serve as a guide to structurally and functionally relevant residues.14 highly conserved positions and 3 sets of coevolved residues are disclosed.Phylogenetic analysis reveals class‐specific features and divergent group of PLA2s.

Thyroxine-binding globulin deficiency due to a novel SERPINA7 mutation: Clinical characterization, analysis of X-chromosome inactivation pattern and protein structural modeling

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.

Mitochondrial ribosome bL34 mutants present diminished translation of cytochrome c oxidase subunits: MRPL34 codes for a component of subunit 54S of mitoribosome

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.

A new method bridging graph theory and residue co-evolutionary networks for specificity determinant positions detection

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.