Purification and partial physical-chemical characterization of a new bovine trypsin proteoform (zeta-trypsin).

Recent advancements in enzyme research have unveiled a new proteoform of bovine trypsin, expanding our understanding of this well-characterized enzyme. While generally similar to other trypsins, this novel proteoform comprises three polypeptide chains, marking a significant difference in activity, kinetic properties, and conformational stability. Compared with the already known bovine trypsin proteoforms, the results showed a lower: activity, kcat and kcat.KM-1 and protein 'foldedness' ratio for the new proteoform. Molecular autolysis, a common feature in trypsin and chymotrypsin, has been explored through comparative physical chemistry properties with other proteoforms. This new proteoform of trypsin not only enriches the existing enzyme repertoire but also promises to shed light on the intricate physiological pathway for enzyme inactivation. Our results suggest that the new trypsin proteoform is one of the likely final pathways for enzyme inactivation in a physiological environment. This discovery opens up new avenues for further research into the functional implications of this new trypsin proteoform.

Design of a multi-epitope vaccine (vme-VAC/MST-1) against cholera and vibriosis based on reverse vaccinology and immunoinformatics approaches.

Vibriosis and cholera are serious diseases distributed worldwide and caused by six marine bacteria of the Vibrio genus. Thousands of deaths occur each year due to these illnesses, necessitating the development of new preventive measures. Presently, the existing cholera vaccine demonstrates an effectiveness of approximately 60%. Here we describe a new multi-epitope vaccine, 'vme-VAC/MST-1' based on vaccine targets identified by reverse vaccinology and epitopes predicted by immunoinformatics, two currently effective tools for predicting new vaccines for bacterial pathogens. The vaccine was designed to combat vibriosis and cholera by incorporating epitopes predicted for CTL, HTL, and B cells. These epitopes were identified from six vaccine targets revealed through subtractive genomics, combined with reverse vaccinology, and were further filtered using immunoinformatics approaches based on their predicted immunogenicity. To construct the vaccine, 28 epitopes (24 CTL/B and 4 HTL/B) were linked to the sequence of the cholera toxin B subunit adjuvant. In silico analyses indicate that the resulting immunogen is stable, soluble, non-toxic, and non-allergenic. Furthermore, it exhibits no homology to the host and demonstrates a strong capacity to elicit innate, B-cell, and T-cell immune responses. Our analysis suggests that it is likely to elicit immune reactions mediated through the TLR5 pathway, as evidenced by the molecular docking of the vaccine with the receptor, which revealed high affinity and a favorable reaction. Thus, vme-VAC/MST-1 is predicted to be a safe and effective solution against pathogenic Vibrio spp. However, further experimental analyses are required to measure the vaccine's effects In vivo.Communicated by Ramaswamy H. Sarma.

Prediction of enzymatic function with high efficiency and a reduced number of features using genetic algorithm.

The post-genomic era has raised a growing demand for efficient procedures to identify protein functions, which can be accomplished by applying machine learning to the characteristics set extracted from the protein. This approach is feature-based and has been the focus of several works in bioinformatics. In this work, we investigated the characteristics of proteins, representing the primary, secondary, tertiary, and quaternary structures of the protein, that improve the model's quality by applying dimensionality reduction techniques and using the Support Vector Machine classifier for predicting the enzymes' classes. During the investigation, two approaches were evaluated: feature extraction/transformation, which was performed using the statistical technique Factor Analysis, and feature selection methods. For feature selection, we proposed an approach based on a genetic algorithm to face the optimization conflict between the simplicity and reliability of an ideal representation of the characteristics of the enzymes and also compared and employed other methods for this purpose. The best result was accomplished using a feature subset generated by our implementation of a multi-objective genetic algorithm enriched with features that this work identified as relevant to represent the enzymes. This subset representation reduced the dataset by about 87% and reached 85.78% of F-measure performance, improving the overall quality of the model classification. In addition, we verified in this work a subset addressed with only 28 features out of a total of 424 that reached a performance above 80% of F-measure for four of the six evaluated classes, showing that satisfactory classification performance can be achieved with a reduced number of enzymes's characteristics. The datasets and implementations are openly available.

Structural and biochemical characterization of Leptospira interrogans Lsa45 reveals a penicillin-binding protein with esterase activity.

Leptospirosis is a bacterial disease that affects humans and animals and is caused by Leptospira. The recommended treatment for leptospirosis is antibiotic therapy, which should be given early in the course of the disease. Despite the use of these antibiotics, their role during the course of the disease is still not completely clear because of the lack of effective clinical trials, particularly for severe cases of the disease. Here, we present the characterization of L. interrogans Lsa45 protein by gel filtration, protein crystallography, SAXS, fluorescence and enzymatic assays. The oligomeric studies revealed that Lsa45 is monomeric in solution. The crystal structure of Lsa45 revealed the presence of two subdomains: a large α/ß subdomain and a small α-helical subdomain. The large subdomain contains the amino acids Ser122, Lys125, and Tyr217, which correspond to the catalytic triad that is essential for ß-lactamase or serine hydrolase activity in similar enzymes. Additionally, we also confirmed the bifunctional promiscuity of Lsa45, in hydrolyzing both the 4-nitrophenyl acetate (p-NPA) and nitrocefin ß-lactam antibiotic. Therefore, this study provides novel insights into the structure and function of enzymes from L. interrogans, which furthers our understanding of this bacterium and the development of new therapies for the prevention and treatment of leptospirosis.

Naview: A d3.js Based JavaScript Library for Drawing and Annotating Voltage-Gated Sodium Channels Membrane Diagrams.

Voltage-gated sodium channels (Nav) are membrane proteins essential to initiating and propagating action potential in neurons and other excitable cells. For a given organism there are often multiple, specialized sodium channels found in different tissues, whose mutations can cause deleterious effects observed in numerous diseases. Consequently, there is high medical and pharmacological interest in these proteins. Scientific literature often uses membrane diagrams to depict important patterns in these channels including the six transmembrane segments (S1-S6) present in four different homologous domains (D1-D4), the S4 voltage sensors, the pore-lining residue segments and the ion selectivity filter residues, glycosylation and phosphorylation residues, toxin binding sites and the inactivation loop, among others. Most of these diagrams are illustrated either digitally or by hand and programs specifically dedicated to the interactive and data-friendly generation of such visualizations are scarce or non-existing. This paper describes Naview, an open-source javascript visualization compatible with modern web browsers for the dynamic drawing and annotation of voltage-gated sodium channels membrane diagrams based on the D3.js library. By using a graphical user interface and combining user-defined annotations with optional UniProt code as inputs, Naview allows the creation and customization of membrane diagrams. In this interface, a user can also map and display important sodium channel properties, residues, regions and their relationships through symbols, colors, and edge connections. Such features can facilitate data exploration and provide fast, high-quality publication-ready graphics for this highly active area of research.

Unexpected plasticity of the quaternary structure of iron-manganese superoxide dismutases.

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.

Characterization of an Allosteric Pocket in Zika Virus NS2B-NS3 Protease.

The NS2B-NS3 protease from Zika virus (ZIKV NS2B-NS3pro) cleaves the viral polyprotein, being essential for its replication and a therapeutic target. Inhibitors that target the active site of ZIKV NS2B-NS3pro have been developed, but they tend to have unfavorable pharmacokinetic properties due to their highly positive charge. Thus, the characterization of allosteric sites in this protease provides new strategies for inhibitor development. Here, we characterized a new allosteric pocket in ZIKV NS2B-NS3pro, analogous to the one previously described for the dengue virus protease. Molecular dynamics simulations indicate the presence of cavities around the residue Ala125, sampling protein conformations in which they are connected to the active site. This link between the residue Ala125 and the active site residues was reinforced by correlation network analysis. To experimentally verify the existence of this allosteric mechanism, we expressed and purified the Ala125Cys mutant of ZIKV NS2B-NS3pro and demonstrated that this variant is inhibited by the thiol-containing chemical probes 5,5'-dithiobis-(2-nitrobenzoic acid) and aldrithiol, which do not affect the activity of the wild-type protein. Inhibition of the mutant protein is reversed by the addition of strong reducing agents, supporting the involvement of Cys125 in covalent bond formation and enzyme inhibition. Together, our results provide experimental evidence for an allosteric pocket in ZIKV NS2B-NS3pro, in the region around Ala125, and computational insights on the structural connection between this region and the enzyme active site.

Resurrecting Golgi proteins to grasp Golgi ribbon formation and self-association under stress.

The Golgi complex is an essential organelle of the eukaryotic exocytic pathway. A subfamily of Golgi matrix proteins, called GRASPs, is central in stress-induced unconventional secretion, Golgi dynamics during mitosis/apoptosis, and Golgi ribbon formation. The Golgi ribbon is vertebrate-specific and correlates with the appearance of two GRASP paralogues and two Golgins (GM130/Golgin45), which form specific GRASP-Golgin pairs. The molecular details of their appearance only in Metazoans are unknown. Moreover, despite new functionalities supported by GRASP paralogy, little is known about their structural and evolutionary differences. Here, we used ancestor sequence reconstruction and biophysical/biochemical approaches to assess the evolution of GRASPs structure/dynamics, fibrillation, and how they started anchoring their Golgin partners. Our data showed that a GRASP ancestor anchored Golgins before gorasp gene duplication in Metazoans. After gene duplication, variations within the GRASP binding pocket determined which paralogue would recruit which Golgin. These interactions are responsible for their specific Golgi location and Golgi ribbon appearance. We also suggest that GRASPs have a long-standing capacity to form supramolecular structures, affecting their participation in stress-induced processes.

Coq3p relevant residues for protein activity and stability.

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.

Reenacting the Birth of a Function: Functional Divergence of HIUases and Transthyretins as Inferred by Evolutionary and Biophysical Studies.

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.

Glutantßase: a database for improving the rational design of glucose-tolerant ß-glucosidases.

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the ß-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, ß-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of ß-glucosidase structures, called Glutantßase. Our database includes 3842 GH1 ß-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantßase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant ß-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant ß-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient ß-glucosidases. We hope that Glutantßase can help improve second-generation biofuel production. Glutantßase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

Coevolved Positions Represent Key Functional Properties in the Trypsin-Like Serine Proteases Protein Family.

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.

Structural and immunological characterization of a new nucleotidyltransferase-like antigen from Paracoccidioides brasiliensis.

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.

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.

Closed- and open-state models of human skeletal muscle sodium channel.

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.

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.

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

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.

PFstats: A Network-Based Open Tool for Protein Family Analysis.

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.

Mitochondrial ribosome bL34 mutants present diminished translation of cytochrome c oxidase subunits.

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.

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.