Obtaining extremely large and accurate protein multiple sequence alignments from curated hierarchical alignments

Genome-wide identification and spatiotemporal expression analysis of cadherin superfamily members in echinoderms

BACKGROUND

Cadherins are calcium-dependent transmembrane cell-cell adhesion proteins that are essential for metazoan development. They consist of three subfamilies: classical cadherins, which bind catenin, protocadherins, which contain 6-7 calcium-binding repeat domains, and atypical cadherins. Their functions include forming adherens junctions, establishing planar cell polarity (PCP), and regulating cell shape, proliferation, and migration. Because they are basal deuterostomes, echinoderms provide important insights into bilaterian evolution, but their only well-characterized cadherin is G-cadherin, a classical cadherin that is expressed by many embryonic epithelia. We aimed to better characterize echinoderm cadherins by conducting phylogenetic analyses and examining the spatiotemporal expression patterns of cadherin-encoding genes during Strongylocentrotus purpuratus development.

RESULTS

Our phylogenetic analyses conducted on two echinoid, three asteroid, and one crinoid species identified ten echinoderm cadherins, including one deuterostome-specific ortholog, cadherin-23, and an echinoderm-specific atypical cadherin that possibly arose in an echinoid-asteroid ancestor. Catenin-binding domains in dachsous-2 orthologs were found to be a deuterostome-specific innovation that was selectively lost in mouse, while those in Fat4 orthologs appeared to be Ambulacraria-specific and were selectively lost in non-crinoid echinoderms. The identified suite of echinoderm cadherins lacks vertebrate-specific innovations but contains two proteins that are present in protostomes and absent from mouse. The spatiotemporal expression patterns of four embryonically expressed cadherins (fat atypical cadherins 1 and 4, dachsous-2, and protocadherin-9) were dynamic and mirrored the expression pattern of Frizzled 5/8, a non-canonical Wnt PCP pathway receptor protein essential for archenteron morphogenesis.

CONCLUSIONS

The echinoderm cadherin toolkit is more similar to that of an ancient bilaterian predating protostomes and deuterostomes than it is to the suite of cadherins found in extant vertebrates. However, it also appears that deuterostomes underwent several cadherin-related innovations. Based on their similar spatiotemporal expression patterns and orthologous relationships to PCP-related and tumor-suppressing proteins, we hypothesize that sea urchin cadherins may play a role in regulating the shape and growth of embryonic epithelia and organs. Future experiments will examine cadherin expression in non-echinoid echinoderms and explore the functions of cadherins during echinoderm development.

DeepHomo2.0: improved protein-protein contact prediction of homodimers by transformer-enhanced deep learning

Protein-protein interactions play an important role in many biological processes. However, although structure prediction for monomer proteins has achieved great progress with the advent of advanced deep learning algorithms like AlphaFold, the structure prediction for protein-protein complexes remains an open question. Taking advantage of the Transformer model of ESM-MSA, we have developed a deep learning-based model, named DeepHomo2.0, to predict protein-protein interactions of homodimeric complexes by leveraging the direct-coupling analysis (DCA) and Transformer features of sequences and the structure features of monomers. DeepHomo2.0 was extensively evaluated on diverse test sets and compared with eight state-of-the-art methods including protein language model-based, DCA-based and machine learning-based methods. It was shown that DeepHomo2.0 achieved a high precision of >70% with experimental monomer structures and >60% with predicted monomer structures for the top 10 predicted contacts on the test sets and outperformed the other eight methods. Moreover, even the version without using structure information, named DeepHomoSeq, still achieved a good precision of >55% for the top 10 predicted contacts. Integrating the predicted contacts into protein docking significantly improved the structure prediction of realistic Critical Assessment of Protein Structure Prediction homodimeric complexes. DeepHomo2.0 and DeepHomoSeq are available at http://huanglab.phys.hust.edu.cn/DeepHomo2/.

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ⁵⁴, and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ⁵⁴ binding loops. Based on the DnaB analysis, we hypothesize that DnaB 'grabs' ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

eCOMPASS: evaluative comparison of multiple protein alignments by statistical score

Motivation

Detecting subtle biologically relevant patterns in protein sequences often requires the construction of a large and accurate multiple sequence alignment (MSA). Methods for constructing MSAs are usually evaluated using benchmark alignments, which, however, typically contain very few sequences and are therefore inappropriate when dealing with large numbers of proteins.

Results

eCOMPASS addresses this problem using a statistical measure of relative alignment quality based on direct coupling analysis (DCA): to maintain protein structural integrity over evolutionary time, substitutions at one residue position typically result in compensating substitutions at other positions. eCOMPASS computes the statistical significance of the congruence between high scoring directly coupled pairs and 3D contacts in corresponding structures, which depends upon properly aligned homologous residues. We illustrate eCOMPASS using both simulated and real MSAs.

Availability and implementation

The eCOMPASS executable, C++ open source code and input data sets are available at https://www.igs.umaryland.edu/labs/neuwald/software/compass

Supplementary information

Supplementary data are available at Bioinformatics online.