Editorial: Machine learning, epistasis, and protein engineering: From sequence-structure-function relationships to regulation of metabolic pathways
Frederic Cadet
Emma Saavedra
Per-Olof Syren
Brigitte Gontero- 1Laboratory of Excellence LABEX GR, DSIMB, Inserm UMR S1134, University of Paris City and University of Reunion, Paris, France
- 2PEACCEL, Artificial Intelligence Department, Paris, France
- 3Department of Biochemistry, Instituto Nacional de Cardiología Ignacio Chávez, Mexico City, Mexico
- 4Science for Life Laboratory, School of Engineering Sciences in Chemistry, Biotechnology, and Health, KTH Royal Institute of Technology, Stockholm, Sweden
- 5Department of Fibre and Polymer Technology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, Stockholm, Sweden
- 6Aix Marseille University, CNRS, UMR7281 Bioénergétique et Ingénierie des Protéines, Marseille, France
Editorial on the Research Topic
Machine learning, epistasis, and protein engineering: From sequence-structure-function relationships to regulation of metabolic pathways
Epistasis is a term originating from genomics and describes the non-additivity of effects of gene interactions on functional parameters (Phillips, 2008). Both within and between genes, epistasis plays a fundamental role in the ability of protein sequences to evolve (Breen et al., 2012).
In protein sciences, epistasis reflects non-linearity effects, ie., non-additive impacts, resulting from interactions between mutations within a protein sequence (Reetz, 2013). The evolution of proteins cannot be understood without the knowledge of the epistasis phenomena (non-additive mutational effects) that take place within them. Indeed, epistasis can reverse the effect of a mutation from beneficial to deleterious. Likewise, thanks to epistasis phenomena, the conservation of a neutral mutation during evolution can lead to beneficial effects of greater amplitude. Interactions between mutated amino acids, as well as intramolecular interaction networks that can be set up following mutations, condition the function (Acevedo-Rocha et al., 2021). While it is clear that understanding sequence-structure-function relationships, and in particular intramolecular epistasis, is of paramount importance in protein engineering, the means to predict these phenomena are currently limited. Non-linear interactions are still poorly understood, making the design of networks of interacting amino acid residues serving to introduce desired functionality in an enzyme a bottleneck.
The prediction of protein structure has been an on-going challenge for computational methods, including artificial intelligence. However, the recent development of structure prediction deep learning (DL) tools such as alphafold2 (Jumper et al., 2021), ESMfold (Lin et al., 2022) or ProteinMPNN (Dauparas et al., 2022), has the potential to revolutionize this area (de Brevern, 2022; Goulet and Cambillau, 2022). Nevertheless, these DL tools are not suitable for predicting how individual amino acid changes alter protein structure and function (Eisenstein, 2021): they can’t predict epistatic effects.
After protein folding powered by DeepMind, Meta and/or Baker’s team, the next challenge is to accurately predict epistasis ie., the impact of non-linear interactions of mutations within the protein sequence. A recent review of machine learning and deep learning strategies examines how epistatic effects influence the success rate of protein engineering projects by comparing fifteen state-of-the-art approaches and provides a general workflow for non-experts when using such learning strategies (Cadet et al., 2022).
The aim of this Research Topic is to give an overview of recent advances and to discuss the current understanding of epistatic phenomena. Particular attention is paid to better understanding and modeling of epistatic phenomena that can impair the prediction of a property of interest.
The articles in this Research Topic illustrate the role of non-linear interactions between players through four particularly interesting examples: negative epistasis for an RNA enzyme; higher-order interactions in complex phenotypes revealed by the BowSaw tool; synergistic effect between Lysine acetylation levels, mSWI/SNF activity and BRD9 inhibition in many cellular contexts, and positive epistasis during the acquisition of resistance of pathogenic bacteria to antibiotics.
Each of these articles is summarized below
Machine learning approaches based on large experimental datasets allow the behavior of some complex biological systems to be predicted. As we know, the rational design of RNA enzyme (ribozyme) activity is challenging, and many ribozyme-based systems are only engineered or improved by random mutagenesis and selection (in vitro evolution). Moreover, extensive pairwise and higher-order epistasis prevent straightforward prediction of the effect of multiple mutations that is needed for rational design. Beck et al. used high-throughput experimental data from variants of a self-cleaving ribozyme to train a predictive model through machine learning approaches. Using only sequence and activity data, they showed that a machine learning approach can be used for RNA design even for RNA molecules with unknown structures. This work is very important because self-cleaving ribozymes can be used to engineer control of gene expression thanks to their capacity to alter RNA processing and stability. Beck et al. unveiled negative epistasis in the RNA data suggesting that additional information, eg. thermodynamic stability of helices, might be necessary for increasing accuracy when predicting effects of long-range distance mutations. This will enable a more comprehensive understanding of RNA fitness landscapes for studying evolution and for guiding RNA-based engineering efforts.
Machine learning also enables the inference and classification of cellular, organismal and ecological phenotypes based on large datasets, e.g., from genomic, transcriptomic and metagenomic analyses. DiMucci et al. have developed a suite of algorithms, named BowSaw, that takes into account different combinations of variables also called “rules”. They applied this approach to study the role of the gut microbiome in Crohn’s disease. They found a previously unreported combination of microbial taxa that is broadly and precisely associated with Crohn’s disease samples. This is an important work whose current implementation shows that BowSaw can be applied to other datasets and used to uncover patterns associated with other diseases. This new suit of algorithms utilizes variable interactions in a trained Random Forest (RF) model in order to extract multiple candidate explanatory rules. Applying BowSaw to a study on the role of the gut microbiome in Crohn’s disease, DiMucci et al. shows that it can find a previously unreported combination of microbial taxa that is broadly and precisely associated with Crohn’s disease samples in the data set: they unveil higher-order interactions in complex phenotypes.
Protein-protein interaction networks are often affected by post-translational modifications. Since proteins rarely act alone, these interactions and their regulation play key roles in a cell. For example, acetylated histones not only affect chromatin condensation but also act as anchor points for bromodomain (BRD)-containing adapter proteins. These domains, and the proteins containing them, can act as chromatin scaffolds that organize large interaction networks which regulate transcription. Loehr et al. have created a cellular model allowing the study of lysine acetylation-dependent protein interaction networks. Using cell lines in which histone acetylation was dependent on acetate supplementation, they showed that the loss of lysine acetylation remodels the composition of the chromatin and the protein-protein interactions of BRDs. Their model may help to understand the functional characterization of BRD-containing proteins across distinct cellular contexts. Synergy between Kac levels, mSWI/SNF activity and BRD9 inhibition appears to be present in many cellular contexts. Loehr et al. created an inexpensive and flexible cellular model allowing the study of Kac-dependent protein interaction networks.
The increasing resistance of pathogenic bacteria to antibiotics is a major challenge to the health of humans, livestock and wildlife. One approach to overcome resistance is to use drug combinations. However, Mehta et al. showed that in the case of the pathogen Francisella tularensis, exposure to two different drugs, sequentially or in combination, resulted in a generalist mutation followed by further mutations that alternated between adaptation to one drug or the other, hence overcoming resistance. Mehta et al. showed that clonal interference, weak pleiotropy and positive epistasis also contributed to combinatorial evolution. The finding suggests that the use of this non-interacting drug pair against F. tularensis may render both drugs ineffective because of mutational switch-backs that accelerate evolution of dual resistance. A better understanding of how drug combinations affect adaptation to multi-drug resistance is required to help overcome antibiotic resistance.
Our understanding of how to account for the effects of genetic variation and in particular the impact of non-linear interactions both within and between genes, is still in its infancy (Cadet et al., 2022; Wittmund et al., 2022). As shown in the articles cited above, the fundamental question of epistasis must be approached in a transdisciplinary manner mixing theoretical approaches, data driven modeling and wet laboratory experiments.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Conflict of interest
Author FC is a member of the scientific board of Peaccel.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Acevedo-Rocha, C. G., Li, A., D’Amore, L., Hoebenreich, S., Sanchis, J., Lubrano, P., et al. (2021). Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics. Nat. Commun. 12, 1621. doi:10.1038/s41467-021-21833-w
Breen, M. S., Kemena, C., Vlasov, P. K., Notredame, C., and Kondrashov, F. A. (2012). Epistasis as the primary factor in molecular evolution. Nature 490, 535–538. doi:10.1038/nature11510
Cadet, X. F., Gelly, J. C., van Noord, A., Cadet, F., and Acevedo-Rocha, C. G. (2022). “Learning strategies in protein directed evolution,” in Directed evolution, methods in molecular biology. Editors A. Currin, and N. Swainston (New York, NY: Springer US), 225–275. doi:10.1007/978-1-0716-2152-3_15
Dauparas, J., Anishchenko, I., Bennett, N., Bai, H., Ragotte, R. J., Milles, L. F., et al. (2022). Robust deep learning–based protein sequence design using ProteinMPNN. Science 378, 49–56. doi:10.1126/science.add2187
de Brevern, A. G. (2022). An agnostic analysis of the human AlphaFold2 proteome using local protein conformations. Biochimie. doi:10.1016/j.biochi.2022.11.009
Eisenstein, M. (2021). Artificial intelligence powers protein-folding predictions. Nature 599, 706–708. doi:10.1038/d41586-021-03499-y
Goulet, A., and Cambillau, C. (2022). Present impact of AlphaFold2 revolution on structural biology, and an illustration with the structure prediction of the bacteriophage J-1 host adhesion device. Front. Mol. Biosci. 9, 907452. doi:10.3389/fmolb.2022.907452
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589. doi:10.1038/s41586-021-03819-2
Lin, Z., Akin, H., Rao, R., Hie, B., Zhu, Z., Lu, W., et al. (2022). Evolutionary-scale prediction of atomic level protein structure with a language model (preprint). Synth. Biol. doi:10.1101/2022.07.20.500902
Phillips, P. C. (2008). Epistasis — The essential role of gene interactions in the structure and evolution of genetic systems. Nat. Rev. Genet. 9, 855–867. doi:10.1038/nrg2452
Reetz, M. T. (2013). The importance of additive and non-additive mutational effects in protein engineering. Angew. Chem. Int. Ed. Engl. 52, 2658–2666. doi:10.1002/anie.201207842
Keywords: epistasis, non-linear interactions, machine learning, artificial intelligence, RNA enzyme, Crohn ‘s disease, antibiotic resistance, Protein-protein interaction networks
Citation: Cadet F, Saavedra E, Syren P-O and Gontero B (2022) Editorial: Machine learning, epistasis, and protein engineering: From sequence-structure-function relationships to regulation of metabolic pathways. Front. Mol. Biosci. 9:1098289. doi: 10.3389/fmolb.2022.1098289
Received: 14 November 2022; Accepted: 25 November 2022;
Published: 02 December 2022.
Edited and reviewed by:
Andrea Mozzarelli, University of Parma, ItalyCopyright © 2022 Cadet, Saavedra, Syren and Gontero. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Frederic Cadet, frederic.cadet.run@gmail.com