Search

Coarse-grained models have provided researchers with greatly improved computational efficiency in modeling structures and dynamics of biomacromolecules, but, to be practically useful, they need fast and accurate conversion methods back to the all-atom representation. Reconstruction of atomic details may also be required in the case of some experimental methods, like electron microscopy, which may provide Cα-only structures. In this contribution, we present a new method for recovery of all backbone atom positions from just the Cα coordinates. Our approach, called deepBBQ, uses a deep convolutional neural network to predict a single internal coordinate per peptide plate, based on Cα trace geometric features, and then proceeds to recalculate the cartesian coordinates based on the assumption that the peptide plate atoms lie in the same plane. Extensive comparison with similar programs shows that our solution is accurate and cost-efficient. The deepBBQ program is available as part of the open-source bioinformatics toolkit Bioshell and is free for download and the documentation is available online.

Ferredoxins, iron-sulfur (Fe-S) cluster proteins, play a key role in oxidoreduction reactions. To date, evolutionary analysis of these proteins across the domains of life have been confined to observing the abundance of Fe-S cluster types (2Fe-2S, 3Fe-4S, 4Fe-4S, 7Fe-8S (3Fe-4s and 4Fe-4S) and 2[4Fe-4S]) and the diversity of ferredoxins within these cluster types was not studied. To address this research gap, here we propose a subtype classification and nomenclature for ferredoxins based on the characteristic spacing between the cysteine amino acids of the Fe-S binding motif as a subtype signature to assess the diversity of ferredoxins across the living organisms. To test this hypothesis, comparative analysis of ferredoxins between bacterial groups, Alphaproteobacteria and Firmicutes and ferredoxins collected from species of different domains of life that are reported in the literature has been carried out. Ferredoxins were found to be highly diverse within their types. Large numbers of alphaproteobacterial species ferredoxin subtypes were found in Firmicutes species and the same ferredoxin subtypes across the species of Bacteria, Archaea, and Eukarya, suggesting shared common ancestral origin of ferredoxins between Archaea and Bacteria and lateral gene transfer of ferredoxins from prokaryotes (Archaea/Bacteria) to eukaryotes. This study opened new vistas for further analysis of diversity of ferredoxins in living organisms.

The assignment of secondary structure elements in protein conformations is necessary to interpret a protein model that has been established by computational methods. The process essentially involves labeling the amino acid residues with H (Helix), E (Strand), or C (Coil, also known as Loop). When particular atoms are absent from an input protein structure, the procedure becomes more complicated, especially when only the alpha carbon locations are known. Various techniques have been tested and applied to this problem during the last forty years. The application of machine learning techniques is the most recent trend. This contribution presents the HECA classifier, which uses neural networks to assign protein secondary structure types. The technique exclusively employs Cα coordinates. The Keras (TensorFlow) library was used to implement and train the neural network model. The BioShell toolkit was used to calculate the neural network input features from raw coordinates. The study’s findings show that neural network-based methods may be successfully used to take on structure assignment challenges when only Cα trace is available. Thanks to the careful selection of input features, our approach’s accuracy (above 97%) exceeded that of the existing methods.

Cytochrome P450 monooxygenases (CYPs/P450s) are heme thiolate proteins present in species across the biological kingdoms. By virtue of their broad substrate promiscuity and regio- and stereo-selectivity, these enzymes enhance or attribute diversity to secondary metabolites. Actinomycetes species are well-known producers of secondary metabolites, especially Salinispora species. Despite the importance of P450s, a comprehensive comparative analysis of P450s and their role in secondary metabolism in Salinispora species is not reported. We therefore analyzed P450s in 126 strains from three different species Salinispora arenicola, S. pacifica, and S. tropica. The study revealed the presence of 2643 P450s that can be grouped into 45 families and 103 subfamilies. CYP107 and CYP125 families are conserved, and CYP105 and CYP107 families are bloomed (a P450 family with many members) across Salinispora species. Analysis of P450s that are part of secondary metabolite biosynthetic gene clusters (smBGCs) revealed Salinispora species have an unprecedented number of P450s (1236 P450s-47%) part of smBGCs compared to other bacterial species belonging to the genera Streptomyces (23%) and Mycobacterium (11%), phyla Cyanobacteria (8%) and Firmicutes (18%) and the classes Alphaproteobacteria (2%) and Gammaproteobacteria (18%). A peculiar characteristic of up to six P450s in smBGCs was observed in Salinispora species. Future characterization Salinispora species P450s and their smBGCs have the potential for discovering novel secondary metabolites.

The impact of lifestyle on shaping the genome content of an organism is a well-known phenomenon and cytochrome P450 enzymes (CYPs/P450s), heme-thiolate proteins that are ubiquitously present in organisms, are no exception. Recent studies focusing on a few bacterial species such as Streptomyces, Mycobacterium, Cyanobacteria and Firmicutes revealed that the impact of lifestyle affected the P450 repertoire in these species. However, this phenomenon needs to be understood in other bacterial species. We therefore performed genome data mining, annotation, phylogenetic analysis of P450s and their role in secondary metabolism in the bacterial class Gammaproteobacteria. Genome-wide data mining for P450s in 1261 Gammaproteobacterial species belonging to 161 genera revealed that only 169 species belonging to 41 genera have P450s. A total of 277 P450s found in 169 species grouped into 84 P450 families and 105 P450 subfamilies, where 38 new P450 families were found. Only 18% of P450s were found to be involved in secondary metabolism in Gammaproteobacterial species, as observed in Firmicutes as well. The pathogenic or commensal lifestyle of Gammaproteobacterial species influences them to such an extent that they have the lowest number of P450s compared to other bacterial species, indicating the impact of lifestyle on shaping the P450 repertoire. This study is the first report on comprehensive analysis of P450s in Gammaproteobacteria.

BioShell is an open-source package for processing biological data, particularly focused on structural applications. The package provides parsers, data structures and algorithms for handling and analyzing macromolecular sequences, structures and sequence profiles. The most frequently used routines are accessible by a set of easy-to-use command line utilities for a Linux environment. The full functionality of the package assumes knowledge of C++ or Python to assemble an application using this software library. Since the last publication that announced the version 2.0, the package has been greatly expanded and rewritten in C++ standard 11 (C++11) to improve its modularity and efficiency. A new testing platform has been implemented to continuously test the correctness and integrity of the package. More than two hundred test programs have been published to provide simple examples that can be used as templates. This makes BioShell an easy to use library that greatly speeds up development of bioinformatics applications and web services without compromising computational efficiency.