COMER2: GPU-accelerated sensitive and specific homology searches

Gardnerella fibrinogen-binding protein as a candidate adherence factor

Bacterial vaginosis (BV), a form of vaginal dysbiosis, is associated with numerous adverse reproductive and obstetric outcomes. Gardnerella spp. are among the key bacteria identified in most BV cases. The formation of a polymicrobial Gardnerella-dominated biofilm on the vaginal epithelium is a characteristic diagnostic marker of BV. Gardnerella colonization and biofilm formation indicate a significant adhesion potential, the determinants of which remain unexplored. In this initial approach to identify Gardnerella adhesins, we analyzed the Cna protein located on the G. vaginalis ATCC 14018 cell surface as determined previously. Structure modeling of Cna (designated Grd Cna) revealed that the protein contains N2 and N3 domains with an immunoglobulin (IgG)-like fold, which shows structural homology to the corresponding domains in SdrD and UafA proteins of the microbial surface component recognizing adhesive matrix molecules (MSCRAMMs) family. A single B domain shares structural similarity with the corresponding domain of Sdr proteins. The R region is rich in PKD repeats, while the C-terminal contains a non-canonical LVNTG cell wall sorting motif. The cna gene was predominantly detected in G. vaginalis isolates but was absent in other commonly identified Gardnerella species isolates. The recombinant Grd Cna protein binds dose-dependently to human fibrinogen but does not interact with fibronectin or collagen types I, III, or IV. Cna-positive G. vaginalis cells adhered to immobilized fibrinogen; however, recombinant Cna did not inhibit this binding, suggesting that Cna may not be a major adhesin mediating G. vaginalis adherence to this ECM component.

GTalign: spatial index-driven protein structure alignment, superposition, and search

With protein databases growing rapidly due to advances in structural and computational biology, the ability to accurately align and rapidly search protein structures has become essential for biological research. In response to the challenge posed by vast protein structure repositories, GTalign offers an innovative solution to protein structure alignment and search-an algorithm that achieves optimal superposition at high speeds. Through the design and implementation of spatial structure indexing, GTalign parallelizes all stages of superposition search across residues and protein structure pairs, yielding rapid identification of optimal superpositions. Rigorous evaluation across diverse datasets reveals GTalign as the most accurate among structure aligners while presenting orders of magnitude in speedup at state-of-the-art accuracy. GTalign's high speed and accuracy make it useful for numerous applications, including functional inference, evolutionary analyses, protein design, and drug discovery, contributing to advancing understanding of protein structure and function.

Frequent nonhomologous replacement of replicative helicase loaders by viruses in Vibrionaceae

Several microbial genomes lack textbook-defined essential genes. If an essential gene is absent from a genome, then an evolutionarily independent gene of unknown function complements its function. Here, we identified frequent nonhomologous replacement of an essential component of DNA replication initiation, a replicative helicase loader gene, in Vibrionaceae. Our analysis of Vibrionaceae genomes revealed two genes with unknown function, named vdhL1 and vdhL2, that were substantially enriched in genomes without the known helicase-loader genes. These genes showed no sequence similarities to genes with known function but encoded proteins structurally similar with a viral helicase loader. Analyses of genomic syntenies and coevolution with helicase genes suggested that vdhL1/2 encodes a helicase loader. The in vitro assay showed that Vibrio harveyi VdhL1 and Vibrio ezurae VdhL2 promote the helicase activity of DnaB. Furthermore, molecular phylogenetics suggested that vdhL1/2 were derived from phages and replaced an intrinsic helicase loader gene of Vibrionaceae over 20 times. This high replacement frequency implies the host's advantage in acquiring a viral helicase loader gene.

Prediction of protein assemblies by structure sampling followed by interface-focused scoring

Proteins often function as part of permanent or transient multimeric complexes, and understanding function of these assemblies requires knowledge of their three-dimensional structures. While the ability of AlphaFold to predict structures of individual proteins with unprecedented accuracy has revolutionized structural biology, modeling structures of protein assemblies remains challenging. To address this challenge, we developed a protocol for predicting structures of protein complexes involving model sampling followed by scoring focused on the subunit-subunit interaction interface. In this protocol, we diversified AlphaFold models by varying construction and pairing of multiple sequence alignments as well as increasing the number of recycles. In cases when AlphaFold failed to assemble a full protein complex or produced unreliable results, additional diverse models were constructed by docking of monomers or subcomplexes. All the models were then scored using a newly developed method, VoroIF-jury, which relies only on structural information. Notably, VoroIF-jury is independent of AlphaFold self-assessment scores and therefore can be used to rank models originating from different structure prediction methods. We tested our protocol in CASP15 and obtained top results, significantly outperforming the standard AlphaFold-Multimer pipeline. Analysis of our results showed that the accuracy of our assembly models was capped mainly by structure sampling rather than model scoring. This observation suggests that better sampling, especially for the antibody-antigen complexes, may lead to further improvement. Our protocol is expected to be useful for modeling and/or scoring protein assemblies.

The COMER web server for protein analysis by homology

SUMMARY

Sequence homology is a basic concept in protein evolution, structure and function studies. However, there are not many different tools and services for homology searches being sensitive, accurate and fast at the same time. We present a new web server for protein analysis based on COMER2, a sequence alignment and homology search method that exhibits these characteristics. COMER2 has been upgraded since its last publication to improve its alignment quality and ease of use. We demonstrate how the user can benefit from using it by providing examples of extensive annotation of proteins of unknown function. Among the distinctive features of the web server is the user's ability to submit multiple queries with one click of a button. This and other features allow for transparently running homology searches-in a command-line, programmatic or graphical environment-across multiple databases with multiple queries. They also promote extensive simultaneous protein analysis at the sequence, structure and function levels.

AVAILABILITY AND IMPLEMENTATION

The COMER web server is available at https://bioinformatics.lt/comer.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein sequence-to-structure learning: Is this the end(-to-end revolution)?

The potential of deep learning has been recognized in the protein structure prediction community for some time, and became indisputable after CASP13. In CASP14, deep learning has boosted the field to unanticipated levels reaching near-experimental accuracy. This success comes from advances transferred from other machine learning areas, as well as methods specifically designed to deal with protein sequences and structures, and their abstractions. Novel emerging approaches include (i) geometric learning, that is, learning on representations such as graphs, three-dimensional (3D) Voronoi tessellations, and point clouds; (ii) pretrained protein language models leveraging attention; (iii) equivariant architectures preserving the symmetry of 3D space; (iv) use of large meta-genome databases; (v) combinations of protein representations; and (vi) finally truly end-to-end architectures, that is, differentiable models starting from a sequence and returning a 3D structure. Here, we provide an overview and our opinion of the novel deep learning approaches developed in the last 2 years and widely used in CASP14.