Protein Sequence Similarity Research Articles

Two novel DnaJ chaperone proteins CG5001 and P58IPK regulate the pathogenicity of Huntington’s disease related aggregates

Huntington’s disease (HD) is a rare neurodegenerative disease caused due to aggregation of Huntingtin (HTT) protein. This study involves the cloning of 40 DnaJ chaperones from Drosophila, and overexpressing them in yeasts and fly models of HD. Accordingly, DnaJ chaperones were catalogued as enhancers or suppressors based on their growth phenotypes and aggregation properties. 2 of the chaperones that came up as targets were CG5001 and P58IPK. Protein aggregation and slow growth phenotype was rescued in yeasts, S2 cells, and Drosophila transgenic lines of HTT103Q with these overexpressed chaperones. Since DnaJ chaperones have protein sequence similarity across species, they can be used as possible tools to combat the effects of neurodegenerative diseases.

Enhancing protein-ligand binding affinity prediction through sequential fusion of graph and convolutional neural networks.

Predicting protein-ligand binding affinity is a crucial and challenging task in structure-based drug discovery. With the accumulation of complex structures and binding affinity data, various machine-learning scoring functions, particularly those based on deep learning, have been developed for this task, exhibiting superiority over their traditional counterparts. A fusion model sequentially connecting a graph neural network (GNN) and a convolutional neural network (CNN) to predict protein-ligand binding affinity is proposed in this work. In this model, the intermediate outputs of the GNN layers, as supplementary descriptors of atomic chemical environments at different levels, are concatenated with the input features of CNN. The model demonstrates a noticeable improvement in performance on CASF-2016 benchmark compared to its constituent CNN models. The generalization ability of the model is evaluated by setting a series of thresholds for ligand extended-connectivity fingerprint similarity or protein sequence similarity between the training and test sets. Masking experiment reveals that model can capture key interaction regions. Furthermore, the fusion model is applied to a virtual screening task for a novel target, PI5P4Kα. The fusion strategy significantly improves the ability of the constituent CNN model to identify active compounds. This work offers a novel approach to enhancing the accuracy of deep learning models in predicting binding affinity through fusion strategies.

P-type calcium ATPases play important roles in biotic and abiotic stress signaling.

Knowledge of Ca2+-ATPases is imperative for improving crop quality/ food security, highly threatened due to global warming. Ca2+-ATPases modulates calcium, essential for stress signaling and modulating growth, development, and immune activities. Calcium is considered a versatile secondary messenger and essential for short- and long-term responses to biotic and abiotic stresses in plants. Coordinated transport activities from both calcium influx and efflux channels are required to generate cellular calcium signals. Various extracellular stimuli cause an induction in cytosolic calcium levels. To cope with such stresses, it is important to maintain intracellular Ca2+ levels. Plants need to evolve efficient efflux mechanisms to maintain Ca2+ ion homeostasis. Plant Ca2+-ATPases are members of the P-type ATPase superfamily and localized in the plasma membrane and endoplasmic reticulum (ER). They are required for various cellular processes, including plant growth, development, calcium signaling, and even retorts to environmental stress. These ATPases play an essential role in Ca2+ homeostasis and are actively involved in Ca2+ transport. Plant Ca2+-ATPases are categorized into two major classes: type IIA and type IIB. Although these two classes of ATPases share similarities in protein sequence, they differ in their structure, cellular localization, and sensitivity to inhibitors. Due to the emerging role of Ca2+-ATPase in abiotic and biotic plant stress, members of this family may help promote agricultural improvement under stress conditions. This review provides a comprehensive overview of P-type Ca2+-ATPase, and their role in Ca2+ transport, stress signaling, and cellular homeostasis focusing on their classification, evolution, ion specificities, and catalytic mechanisms. It also describes the main aspects of the role of Ca2+-ATPase in transducing signals during plant biotic and abiotic stress responses and its role in plant development and physiology.

Protein features fusion using attributed network embedding for predicting protein-protein interaction

BackgroundProtein-protein interactions (PPIs) hold significant importance in biology, with precise PPI prediction as a pivotal factor in comprehending cellular processes and facilitating drug design. However, experimental determination of PPIs is laborious, time-consuming, and often constrained by technical limitations.MethodsWe introduce a new node representation method based on initial information fusion, called FFANE, which amalgamates PPI networks and protein sequence data to enhance the precision of PPIs’ prediction. A Gaussian kernel similarity matrix is initially established by leveraging protein structural resemblances. Concurrently, protein sequence similarities are gauged using the Levenshtein distance, enabling the capture of diverse protein attributes. Subsequently, to construct an initial information matrix, these two feature matrices are merged by employing weighted fusion to achieve an organic amalgamation of structural and sequence details. To gain a more profound understanding of the amalgamated features, a Stacked Autoencoder (SAE) is employed for encoding learning, thereby yielding more representative feature representations. Ultimately, classification models are trained to predict PPIs by using the well-learned fusion feature.ResultsWhen employing 5-fold cross-validation experiments on SVM, our proposed method achieved average accuracies of 94.28%, 97.69%, and 84.05% in terms of Saccharomyces cerevisiae, Homo sapiens, and Helicobacter pylori datasets, respectively.ConclusionExperimental findings across various authentic datasets validate the efficacy and superiority of this fusion feature representation approach, underscoring its potential value in bioinformatics.

Molecular Detection and Characterization of Rickettsia Species in Ixodid Ticks from Selected Regions of Namibia.

Rickettsial pathogens are among the emerging and re-emerging vector-borne zoonoses of public health importance. Reports indicate human exposure to Rickettsial pathogens in Namibia through serological surveys, but there is a lack of data on infection rates in tick vectors, hindering the assessment of the relative risk to humans. Our study sought to screen Ixodid ticks collected from livestock for the presence of Rickettsia species in order to determine infection rates in ticks and to determine the Rickettsia species circulating in the country. We collected and pooled Hyalomma and Rhipicephalus ticks from two adjacent regions of Namibia (Khomas and Otjozondjupa) and observed an overall minimum Rickettsia infection rate of 8.6% (26/304), with an estimated overall pooled prevalence of 9.94% (95% CI: 6.5-14.3). There were no statistically significant differences in the estimated pooled prevalence between the two regions or tick genera. Based on the nucleotide sequence similarity and phylogenetic analysis of the outer membrane protein A (n = 9) and citrate synthase (n = 12) genes, BLAST analysis revealed similarity between Rickettsia africae (n = 2) and Rickettsia aeschlimannii (n = 11), with sequence identities ranging from 98.46 to 100%. Our initial study in Namibia indicates that both zoonotic R. africae and R. aeschlimannii are in circulation in the country, with R. aeschlimannii being the predominant species.

Functional expression, purification, biochemical and biophysical characterizations, and molecular dynamics simulation of a histidine acid phosphatase from Saccharomyces cerevisiae.

A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75min and 55s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, β-Turns, and β-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.

An evolutionary timeline of the oxytocin signaling pathway.

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Fusarium sacchari hypovirus 1, a Member of Hypoviridae with Virulence Attenuation Capacity in Phytopathogenic Fusarium Species.

In a survey of mycoviruses in Fusarium species that cause sugarcane Pokkah boeng disease, twelve Fusarium strains from three Fusarium species (F. sacchari, F. andiyazi, and F. solani) were found to contain Fusarium sacchari hypovirus 1 (FsHV1), which we reported previously. The genomes of these variants range from 13,966 to 13,983 nucleotides, with 98.6% to 99.9% nucleotide sequence identity and 98.70% to 99.9% protein sequence similarity. Phylogenetic analysis placed these FsHV1 variants within the Alphahypovirus cluster of Hypoviridae. Intriguingly, no clear correlation was found between the geographic origin and host specificity of these viral variants. Additionally, six out of the twelve variants displayed segmental deletions of 1.5 to 1.8 kilobases, suggesting the existence of defective viral dsRNA. The presence of defective viral dsRNA led to a two-thirds reduction in the dsRNA of the wild-type viral genome, yet a tenfold increase in the total viral dsRNA content. To standardize virulence across natural strains, all FsHV1 strains were transferred into a single, virus-free Fusarium recipient strain, FZ06-VF, via mycelial fusion. Strains of Fusarium carrying FsHV1 exhibited suppressed pigment synthesis, diminished microspore production, and a marked decrease in virulence. Inoculation tests revealed varying capacities among different FsHV1 variants to modulate fungal virulence, with the strain harboring the FsHV1-FSA1 showing the lowest virulence, with a disease severity index (DSI) of 3.33, and the FsHV1-FS1 the highest (DSI = 17.66). The identification of highly virulent FsHV1 variants holds promise for the development of biocontrol agents for Pokkah boeng management.

PPI3D: a web server for searching, analyzing and modeling protein-protein, protein-peptide and protein-nucleic acid interactions.

Structure-resolved protein interactions with other proteins, peptides and nucleic acids are key for understanding molecular mechanisms. The PPI3D web server enables researchers to query preprocessed and clustered structural data, analyze the results and make homology-based inferences for protein interactions. PPI3D offers three interaction exploration modes: (i) all interactions for proteins homologous to the query, (ii) interactions between two proteins or their homologsand (iii) interactions within a specific PDB entry. The server allows interactive analysis of the identified interactions in both summarized and detailed manner. This includes protein annotations, structures, the interface residues and the corresponding contact surface areas. In addition, users can make inferences about residues at the interaction interface for the query protein(s) from the sequence alignments and homology models. The weekly updated PPI3D database includes all the interaction interfaces and binding sites from PDB, clustered based on both protein sequence and structural similarity, yielding non-redundant datasets without loss of alternative interaction modes. Consequently, the PPI3D users avoid being flooded with redundant information, a typical situation for intensely studied proteins. Furthermore, PPI3D provides a possibility to download user-defined sets of interaction interfaces and analyze them locally. The PPI3D web server is available at https://bioinformatics.lt/ppi3d.

The transcription factor RhMYB17 regulates the homeotic transformation of floral organs in rose (Rosa hybrida) under cold stress.

Low temperatures affect flower development in rose (Rosa hybrida), increasing petaloid stamen number and reducing normal stamen number. We identified the low-temperature-responsive R2R3-MYB transcription factor RhMYB17, which is homologous to Arabidopsis MYB17 by similarity of protein sequences. RhMYB17 was up-regulated at low temperatures, and RhMYB17 transcripts accumulated in floral buds. Transient silencing of RhMYB17 by virus-induced gene silencing decreased petaloid stamen number and increased normal stamen number. According to the ABCDE model of floral organ identity, class A genes APETALA 1 (AP1) and AP2 contribute to sepal and petal formation. Transcription factor binding analysis identified RhMYB17 binding sites in the promoters of rose APETALA 2 (RhAP2) and APETALA 2-LIKE (RhAP2L). Yeast one-hybrid assays, dual-luciferase reporter assays, and electrophoretic mobility shift assays confirmed that RhMYB17 directly binds to the promoters of RhAP2 and RhAP2L, thereby activating their expression. RNA sequencing further demonstrated that RhMYB17 plays a pivotal role in regulating the expression of class A genes, and indirectly influences the expression of the class C gene. This study reveals a novel mechanism for the homeotic transformation of floral organs in response to low temperatures.

Effect of moisture deficiency on the apple stress-associated protein expression level .

In plants, stress-associated proteins (SAP) are involved in response to adverse biotic and abiotic factors. In order to study the drought effect on changes in the SAP expression level in apple trees and their involvement in stress response, an expression level of fourteen genes belonging to this family was analyzed using qPCR. Their participation in a stress response to drought in the genome of the MM-106 apple tree rootstock was shown. At that, the genes demonstrated a tendency to increase their expression by the fourth hour of drought exposure followed by its decrease by the twenty-fourth hour.A comparison of the data obtained during the described experiment with the data obtained when modeling drought under milder conditions showed that under more severe drought conditions, an expression of a greater number of gene encoding stress-associated proteins increases, and a more pronounced response occurs. Evaluation of the effect of an identity degree of nucleotide sequences of gene encoding stress-associated proteins on the similarity of their expression profiles under drought exposure did not reveal any direct relationship between the primary structure of genes and the nature of their expression.In the structure of stress-associated proteins encoded by the genes most strongly response to drought in Malus, Solanum lycopersicum, Gossypium hirsutum, and Cucumis sativus, the A20-AN1 type domain occurs more frequently than other types of zinc finger domains. Analysis of the similarity of stress-associated protein sequences and the expression levels of genes encoding them under drought conditions in different plant species (apple, cotton, tomato, and cucumber) suggests that in each species, the evolution of adaptation mechanisms, including SAP encoding genes, occurred independently and, most likely, after separation of those species from their common ancestor. The results obtained will make it possible to clarify the role of individual genes, encoding SAP, in the formation of an apple stress response to the drought effect.

Transcriptome sequencing and screening of genes related to the MADS-box gene family in Clematis courtoisii.

The MADS-box gene family controls plant flowering and floral organ development; therefore, it is particularly important in ornamental plants. To investigate the genes associated with the MADS-box family in Clematis courtoisii, we performed full-length transcriptome sequencing on C. courtoisii using the PacBio Sequel third-generation sequencing platform, as no reference genome data was available. A total of 12.38 Gb of data, containing 9,476,585 subreads and 50,439 Unigenes were obtained. According to functional annotation, a total of 37,923 Unigenes (75.18% of the total) were assigned with functional annotations, and 50 Unigenes were identified as MADS-box related genes. Subsequently, we employed hmmerscan to perform protein sequence similarity search for the translated Unigene sequences and successfully identified 19 Unigenes associated with the MADS-box gene family, including MIKC*(1) and MIKCC (18) genes. Furthermore, within the MIKCC group, six subclasses can be further distinguished.

The functional divergence of homologous GPAT9 genes contributes to the erucic acid content of Brassica napus seeds

BackgroundThe early allopolyploid Brassica napus was a hybrid of two Brassica species, that had undergone a whole genome duplication event followed by genome restructuring, including deletions and small scale duplications. A large number of homologous genes appeared functional divergence during species domestication. Due to the high conservation of de novo glycerolipid biosynthesis, multiple homologues of glycerol-3-phosphate acyltransferases (GPATs) have been found in B. napus. Moreover, the functional variances among these homologous GPAT-encoding genes are unclear.ResultsIn this study, four B. napus homologous genes encoding glycerol-3-phosphate acyltransferase 9 (BnaGPAT9) were characterized. Although a bioinformatics analysis indicated high protein sequence similarity, the homologues demonstrated tissue-specific expression patterns and functional divergence. Yeast genetic complementation assays revealed that BnaGPAT9-A1/C1 homologues but not BnaGPAT9-A10/C9 homologues encoded functional GPAT enzymes. Furthermore, a single nucleotide polymorphism of BnaGPAT9-C1 that occurred during the domestication process was associated with enzyme activity and contributed to the fatty acid composition. The seed-specific expression of BnGPAT9-C11124A increased the erucic acid content in the transformant seeds.ConclusionsThis study revealed that BnaGPAT9 gene homologues evolved into functionally divergent forms with important roles in erucic acid biosynthesis.

Hiding in plain sight: The discovery of complete genomes of 11 hypothetical spindle-shaped viruses that putatively infect mesophilic ammonia-oxidizing archaea.

The genome of a putative Nitrosopumilaceae virus with a hypothetical spindle-shaped particle morphology was identified in the Yangshan Harbour metavirome from the East China Sea through protein similarity comparison and structure analysis. This discovery was accompanied by a set of 10 geographically dispersed close relatives found in the environmental virus datasets from typical locations of ammonia-oxidizing archaeon distribution. Its host prediction was supported by iPHoP prediction and protein sequence similarity. The structure of the predicted major capsid protein, together with the overall N-glycosylation site, the transmembrane helices prediction, the hydrophilicity profile, and the docking simulation of the major capsid proteins, indicate that these viruses resemble spindle-shaped viruses. It suggests a similarly assembled structure and, consequently, a possibly spindle-shaped morphology of these newly discovered archaeal viruses.

Molecular evolution analysis of MYB5 in Brassicaceae with specific focus on seed coat color of Brassica napus

BackgroundMYB transcription factors are splay a vital role in plant biology, with previous research highlighting the significant impact of the R2R3-MYB-like transcription factor MYB5 on seed mucilage biosynthesis, trichome branching, and seed coat development. However, there is a dearth of studies investigating its role in the regulation of proanthocyanidin (PA) biosynthesis.ResultsIn this study, a total of 51 MYB5 homologous genes were identified across 31 species belonging to the Brassicaceae family, with particular emphasis on Brassica napus for subsequent investigation. Through phylogenetic analysis, these genes were categorized into four distinct subclasses. Protein sequence similarity and identity analysis demonstrated a high degree of conservation of MYB5 among species within the Brassicaceae family. Additionally, the examination of selection pressure revealed that MYB5 predominantly underwent purifying selection during its evolutionary history, as indicated by the Ka/Ks values of all MYB5 homologous gene pairs being less than one. Notably, we observed a higher rate of non-synonymous mutations in orthologous genes compared to paralogous genes, and the Ka/Ks value displayed a stronger correlation with Ka. In B. napus, an examination of expression patterns in five tissues revealed that MYB5 exhibited particularly high expression in the black seed coat. The findings from the WGCNA demonstrated a robust correlation between MYB5 and BAN(ANR) associated with PA biosynthesis in the black seed coat, providing further evidence of their close association and co-expression. Furthermore, the results obtained from of the analysis of protein interaction networks offer supplementary support for the proposition that MYB5 possesses the capability to interact with transcriptional regulatory proteins, specifically TT8 and TT2, alongside catalytic enzymes implicated in the synthesis of PAs, thereby making a contribution to the biosynthesis of PAs. These findings imply a plausible and significant correlation between the nuique expression pattern of MYB5 and the pigmentation of rapeseed coats. Nevertheless, additional research endeavors are imperative to authenticate and substantiate these findings.ConclusionsThis study offers valuable insights into the genetic evolution of Brassicaceae plants, thereby serving as a significant reference for the genetic enhancement of Brassicaceae seed coat color.

GSLAlign: community detection and local PPI network alignment

High throughput protein–protein interaction (PPI) profiling and computational techniques have resulted in generating a large amount of PPI network data. The study of PPI networks helps in understanding the biological processes of the proteins. The comparative study of the PPI networks helps in identifying the conserved interactions across the species. This article presents a novel local PPI network aligner ‘GSLAlign’ that consists of two stages. It first detects the communities from the PPI networks by applying the GraphSAGE algorithm using gene expression data. In the second stage, the detected communities are aligned using a community aligner that is based on protein sequence similarity. The community detection algorithm produces more separable and biologically accurate communities as compared to previous community detection algorithms. Moreover, the proposed community alignment algorithm achieves 3–8% better results in terms of semantic similarity as compared to previous local aligners. The average connectivity and coverage of the proposed algorithm are also better than the existing aligners. Communicated by Ramaswamy H. Sarma

SARS-CoV-2 helicase might interfere with cellular nonsense-mediated RNA decay: insights from a bioinformatics study

BackgroundViruses employ diverse strategies to interfere with host defense mechanisms, including the production of proteins that mimic or resemble host proteins. This study aimed to analyze the similarities between SARS-CoV-2 and human proteins, investigate their impact on virus-host interactions, and elucidate underlying mechanisms.ResultsComparing the proteins of SARS-CoV-2 with human and mammalian proteins revealed sequence and structural similarities between viral helicase with human UPF1. The latter is a protein that is involved in nonsense-mediated RNA decay (NMD), an mRNA surveillance pathway which also acts as a cellular defense mechanism against viruses. Protein sequence similarities were also observed between viral nsp3 and human Poly ADP-ribose polymerase (PARP) family of proteins. Gene set enrichment analysis on transcriptomic data derived from SARS-CoV-2 positive samples illustrated the enrichment of genes belonging to the NMD pathway compared with control samples. Moreover, comparing transcriptomic data from SARS-CoV-2-infected samples with transcriptomic data derived from UPF1 knockdown cells demonstrated a significant overlap between datasets.ConclusionsThese findings suggest that helicase/UPF1 sequence and structural similarity might have the ability to interfere with the NMD pathway with pathogenic and immunological implications.

Humanization reveals pervasive incompatibility of yeast and human kinetochore components.

Kinetochores assemble on centromeres to drive chromosome segregation in eukaryotic cells. Humans and budding yeast share most of the structural subunits of the kinetochore, whereas protein sequences have diverged considerably. The conserved centromeric histone H3 variant, CenH3 (CENP-A in humans and Cse4 in budding yeast), marks the site for kinetochore assembly in most species. A previous effort to complement Cse4 in yeast with human CENP-A was unsuccessful; however, co-complementation with the human core nucleosome was not attempted. Previously, our lab successfully humanized the core nucleosome in yeast; however, this severely affected cellular growth. We hypothesized that yeast Cse4 is incompatible with humanized nucleosomes and that the kinetochore represented a limiting factor for efficient histone humanization. Thus, we argued that including the human CENP-A or a Cse4-CENP-A chimera might improve histone humanization and facilitate kinetochore function in humanized yeast. The opposite was true: CENP-A expression reduced histone humanization efficiency, was toxic to yeast, and disrupted cell cycle progression and kinetochore function in wild-type (WT) cells. Suppressors of CENP-A toxicity included gene deletions of subunits of 3 conserved chromatin remodeling complexes, highlighting their role in CenH3 chromatin positioning. Finally, we attempted to complement the subunits of the NDC80 kinetochore complex, individually and in combination, without success, in contrast to a previous study indicating complementation by the human NDC80/HEC1 gene. Our results suggest that limited protein sequence similarity between yeast and human components in this very complex structure leads to failure of complementation.

Developing similarity matrices for antibody-protein binding interactions.

The inventions of AlphaFold and RoseTTAFold are revolutionizing computational protein science due to their abilities to reliably predict protein structures. Their unprecedented successes are due to the parallel consideration of several types of information, one of which is protein sequence similarity information. Sequence homology has been studied for many decades and depends on similarity matrices to define how similar or different protein sequences are to one another. A natural extension of predicting protein structures is predicting the interactions between proteins, but similarity matrices for protein-protein interactions do not exist. This study conducted a mutational analysis of 384 non-redundant antibody-protein antigen complexes to calculate antibody-protein interaction similarity matrices. Every important residue in each antibody and each antigen was mutated to each of the other 19 commonly occurring amino acids and the percentage changes in interaction energies were calculated using three force fields: CHARMM, Amber, and Rosetta. The data were used to construct six interaction similarity matrices, one for antibodies and another for antigens using each force field. The matrices exhibited both commonalities, such as mutations of aromatic and charged residues being the most detrimental, and differences, such as Rosetta predicting mutations of serines to be better tolerated than either Amber or CHARMM. A comparison to nine previously published similarity matrices for protein sequences revealed that the new interaction matrices are more similar to one another than they are to any of the previous matrices. The created similarity matrices can be used in force field specific applications to help guide decisions regarding mutations in protein-protein binding interfaces.

PC_ali: A tool for improved multiple alignments and evolutionary inference based on a hybrid protein sequence and structure similarity score.

Evolutionary inference depends crucially on the quality of multiple sequence alignments (MSA), which is problematic for distantly related proteins. Since protein structure is more conserved than sequence, it seems natural to use structure alignments for distant homologs. However, structure alignments may not be suitable for inferring evolutionary relationships. Here we examined four protein similarity measures that depend on sequence and structure (fraction of aligned residues, sequence identity, fraction of superimposed residues and contact overlap), finding that they are intimately correlated but none of them provides a complete and unbiased picture of conservation in proteins. Therefore, we propose the new hybrid protein sequence and structure similarity score PC_sim based on their main Principal Component. The corresponding divergence measure PC_div shows the strongest correlation with divergences obtained from individual similarities, suggesting that it infers accurate evolutionary divergences. We developed the program PC_ali that constructs protein MSAs either de novo or modifying an input MSA, using a similarity matrix based on PC_sim. tructs a starting MSA based on the maximal cliques of the graph of these PAs and it refines it through progressive alignments along the tree reconstructed with PC_div. Compared with eight state-of-the-art multiple structure or sequence alignment tools, PC_ali achieve higher or equal aligned fraction and structural scores, sequence identity higher than structure aligners although lower than sequence aligners, highest score PC_sim and highest similarity with the MSAs produced by other tools and with the reference MSA Balibase. https://github.com/ugobas/PC_ali. Supplementary data are available at Bioinformatics online.