Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Prediction of Leishmania major Key Proteins Via Topological Analysis of Protein-Protein Interaction Network

Background:

Although leishmaniasis is regarded as a public health problem, no effective vaccine or decisive treatment has been introduced for this disease. Therefore, representing novel therapeutic proteins is essential. Protein-protein Interaction network analysis is a suitable tool to discover novel drug targets for leishmania major. To this aim, gene and protein expression data is used for instructing protein network and the key proteins are highlighted.

Materials and Methods:

In this computational and bioinformatics study, the protein/gene expression data related to leishmania major were studied, and 252 candidate proteins were extracted. Then, the protein networks of these proteins were explored and visualized by using String database and Cytoscape software. Finally, clustering and gene ontology were performed by MCODE and PANTHER databases, respectively.

Results:

Based on gene ontology analysis, most of the leishmania major proteins were located in cell compartments and membrane. Catalytic activity and binding were regarded as the relevant molecular functions and metabolic and cellular processes were the significant biological process. In this network analysis, UB-EP52, EF-2, chaperonin, Hsp70.4, Hsp60, tubulin alpha and beta chain, and ENOL and LACK were introduced as hub-bottleneck proteins. Based on clustering analysis, Lmjf.32.3270, ENOL and Lmjf.13.0290 were determined as seed proteins in each cluster.

Conclusion:

The results indicated that hub proteins play a significant role in pathogenesis and life cycle of leishmania major. Further studies of hubs will provide a better understanding of leishmaniasis mechanisms. Finally, these key hub proteins could be a suitable and helpful potential for drug targets and treating leishmaniasis by considering their validation.