Anne-Florence Raphaëlle Bitbol

Related publications (42)

Combining phylogeny and coevolution improves the inference of interaction partners among paralogous proteins

Author summaryWhen two protein families interact, their sequences feature statistical dependencies. First, interacting proteins tend to share a common evolutionary history. Second, maintaining structure and interactions through the course of evolution yields coevolution, detectable via correlations in the amino-acid usage at contacting sites. Both signals can be used to computationally predict which proteins are specific interaction partners among the paralogs of two interacting protein families, starting just from their sequences. We show that combining them improves the performance of interaction partner inference, especially when the average number of potential partners is large and when the total data set size is modest. The resulting paired multiple-sequence alignments might be used as input to machine-learning algorithms to improve protein-complex structure prediction, as well as to understand interaction specificity in signaling pathways. Predicting protein-protein interactions from sequences is an important goal of computational biology. Various sources of information can be used to this end. Starting from the sequences of two interacting protein families, one can use phylogeny or residue coevolution to infer which paralogs are specific interaction partners within each species. We show that these two signals can be combined to improve the performance of the inference of interaction partners among paralogs. For this, we first align the sequence-similarity graphs of the two families through simulated annealing, yielding a robust partial pairing. We next use this partial pairing to seed a coevolution-based iterative pairing algorithm. This combined method improves performance over either separate method. The improvement obtained is striking in the difficult cases where the average number of paralogs per species is large or where the total number of sequences is modest.

PUBLIC LIBRARY SCIENCE2023

Generative power of a protein language model trained on multiple sequence alignments

Anne-Florence Raphaëlle Bitbol, Umberto Lupo, Damiano Sgarbossa

Computational models starting from large ensembles of evolutionarily related protein sequences capture a representation of protein families and learn constraints associated to protein structure and function. They thus open the possibility for generating novel sequences belonging to protein families. Protein language models trained on multiple sequence alignments, such as MSA Transformer, are highly attractive candidates to this end. We propose and test an iterative method that directly employs the masked language modeling objective to generate sequences using MSA Transformer. We demonstrate that the resulting sequences score as well as natural sequences, for homology, coevolution, and structure-based measures. For large protein families, our synthetic sequences have similar or better properties compared to sequences generated by Potts models, including experimentally validated ones. Moreover, for small protein families, our generation method based on MSA Transformer outperforms Potts models. Our method also more accurately reproduces the higher-order statistics and the distribution of sequences in sequence space of natural data than Potts models. MSA Transformer is thus a strong candidate for protein sequence generation and protein design.

eLIFE SCIENCES PUBL LTD2023

Impact of phylogeny on structural contact inference from protein sequence data

Anne-Florence Raphaëlle Bitbol, Nicola Dietler, Umberto Lupo

Local and global inference methods have been developed to infer structural contacts from multiple sequence alignments of homologous proteins. They rely on correlations in amino acid usage at contacting sites. Because homologous proteins share a common ancestry, their sequences also feature phylogenetic correlations, which can impair contact inference. We investigate this effect by generating controlled synthetic data from a minimal model where the importance of contacts and of phylogeny can be tuned. We demonstrate that global inference methods, specifically Potts models, are more resilient to phylogenetic correlations than local methods, based on covariance or mutual information. This holds whether or not phylogenetic corrections are used, and may explain the success of global methods. We analyse the roles of selection strength and of phylogenetic relatedness. We show that sites that mutate early in the phylogeny yield false positive contacts. We consider natural data and realistic synthetic data, and our findings generalize to these cases. Our results highlight the impact of phylogeny on contact prediction from protein sequences and illustrate the interplay between the rich structure of biological data and inference.

ROYAL SOC2023

Impact of population size on early adaptation in rugged fitness landscapes

Anne-Florence Raphaëlle Bitbol, Richard Marie Servajean

Owing to stochastic fluctuations arising from finite population size, known as genetic drift, the ability of a population to explore a rugged fitness landscape depends on its size. In the weak mutation regime, while the mean steady-state fitness increases with population size, we find that the height of the first fitness peak encountered when starting from a random genotype displays various behaviours versus population size, even among small and simple rugged landscapes. We show that the accessibility of the different fitness peaks is key to determining whether this height overall increases or decreases with population size. Furthermore, there is often a finite population size that maximizes the height of the first fitness peak encountered when starting from a random genotype. This holds across various classes of model rugged landscapes with sparse peaks, and in some experimental and experimentally inspired ones. Thus, early adaptation in rugged fitness landscapes can be more efficient and predictable for relatively small population sizes than in the large-size limit.This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

ROYAL SOC2023

Extracting phylogenetic dimensions of coevolution reveals hidden functional signals

Anne-Florence Raphaëlle Bitbol

Despite the structural and functional information contained in the statistical coupling between pairs of residues in a protein, coevolution associated with function is often obscured by artifactual signals such as genetic drift, which shapes a protein's phylogenetic history and gives rise to concurrent variation between protein sequences that is not driven by selection for function. Here, we introduce a background model for phylogenetic contributions of statistical coupling that separates the coevolution signal due to inter-clade and intra-clade sequence comparisons and demonstrate that coevolution can be measured on multiple phylogenetic timescales within a single protein. Our method, nested coevolution (NC), can be applied as an extension to any coevolution metric. We use NC to demonstrate that poorly conserved residues can nonetheless have important roles in protein function. Moreover, NC improved the structural-contact predictions of several coevolution-based methods, particularly in subsampled alignments with fewer sequences. NC also lowered the noise in detecting functional sectors of collectively coevolving residues. Sectors of coevolving residues identified after application of NC were more spatially compact and phylogenetically distinct from the rest of the protein, and strongly enriched for mutations that disrupt protein activity. Thus, our conceptualization of the phylogenetic separation of coevolution provides the potential to further elucidate relationships among protein evolution, function, and genetic diseases.

NATURE PORTFOLIO2022

Hydrodynamic flow and concentration gradients in the gut enhance neutral bacterial diversity

Anne-Florence Raphaëlle Bitbol

The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.

National Academy of Sciences2022

Protein language models trained on multiple sequence alignments learn phylogenetic relationships

Anne-Florence Raphaëlle Bitbol, Umberto Lupo, Damiano Sgarbossa

Self-supervised neural language models with attention have recently been applied to biological sequence data, advancing structure, function and mutational effect prediction. Some protein language models, including MSA Transformer and AlphaFold's EvoFormer, take multiple sequence alignments (MSAs) of evolutionarily related proteins as inputs. Simple combinations of MSA Transformer's row attentions have led to state-of-the-art unsupervised structural contact prediction. We demonstrate that similarly simple, and universal, combinations of MSA Transformer's column attentions strongly correlate with Hamming distances between sequences in MSAs. Therefore, MSA-based language models encode detailed phylogenetic relationships. We further show that these models can separate coevolutionary signals encoding functional and structural constraints from phylogenetic correlations reflecting historical contingency. To assess this, we generate synthetic MSAs, either without or with phylogeny, from Potts models trained on natural MSAs. We find that unsupervised contact prediction is substantially more resilient to phylogenetic noise when using MSA Transformer versus inferred Potts models. Protein language models taking multiple sequence alignments as inputs capture protein structure and mutational effects. Here, the authors show that these models also encode phylogenetic relationships, and can disentangle correlations due to structural constraints from those due to phylogeny.

NATURE PORTFOLIO2022

Evolutionary history expands the range of signaling interactions in hybrid multikinase networks

Mohamed Barakat, Anne-Florence Raphaëlle Bitbol

Two-component systems (TCSs) are ubiquitous signaling pathways, typically comprising a sensory histidine kinase (HK) and a response regulator, which communicate via intermolecular kinase-to-receiver domain phosphotransfer. Hybrid HKs constitute non-canonical TCS signaling pathways, with transmitter and receiver domains within a single protein communicating via intramolecular phosphotransfer. Here, we report how evolutionary relationships between hybrid HKs can be used as predictors of potential intermolecular and intramolecular interactions ('phylogenetic promiscuity'). We used domain-swap genes chimeras to investigate the specificity of phosphotransfer within hybrid HKs of the GacS-GacA multikinase network of Pseudomonas brassicacearum. The receiver domain of GacS was replaced with those from nine donor hybrid HKs. Three chimeras with receivers from other hybrid HKs demonstrated correct functioning through complementation of a gacS mutant, which was dependent on strains having a functional gacA. Formation of functional chimeras was predictable on the basis of evolutionary heritage, and raises the possibility that HKs sharing a common ancestor with GacS might remain components of the contemporary GacS network. The results also demonstrate that understanding the evolutionary heritage of signaling domains in sophisticated networks allows their rational rewiring by simple domain transplantation, with implications for the creation of designer networks and inference of functional interactions.

NATURE RESEARCH2021

Adapt or Perish: Evolutionary Rescue in a Gradually Deteriorating Environment

Anne-Florence Raphaëlle Bitbol, Loïc Marrec

We investigate the evolutionary rescue of a microbial population in a gradually deteriorating environment, through a combination of analytical calculations and stochastic simulations. We consider a population destined for extinction in the absence of mutants, which can only survive if mutants sufficiently adapted to the new environment arise and fix. We show that mutants that appear later during the environment deterioration have a higher probability to fix. The rescue probability of the population increases with a sigmoidal shape when the product of the carrying capacity and of the mutation probability increases. Furthermore, we find that rescue becomes more likely for smaller population sizes and/or mutation probabilities if the environment degradation is slower, which illustrates the key impact of the rapidity of environment degradation on the fate of a population. We also show that our main conclusions are robust across various types of adaptive mutants, including specialist and generalist ones, as well as mutants modeling antimicrobial resistance evolution. We further express the average time of appearance of the mutants that do rescue the population and the average extinction time of those that do not. Our methods can be applied to other situations with continuously variable fitnesses and population sizes, and our analytical predictions are valid in the weak-to-moderate mutation regime.

2020

Mitochondrial cristae modeled as an out-of-equilibrium membrane driven by a proton field

Anne-Florence Raphaëlle Bitbol

As the places where most of the fuel of the cell, namely, ATP, is synthesized, mitochondria are crucial organelles in eukaryotic cells. The shape of the invaginations of the mitochondria inner membrane, known as a crista, has been identified as a signature of the energetic state of the organelle. However, the interplay between the rate of ATP synthesis and the crista shape remains unclear. In this work, we investigate the crista membrane deformations using a pH-dependent Helfrich model, maintained out of equilibrium by a diffusive flux of protons. This model gives rise to shape changes of a cylindrical invagination, in particular to the formation of necks between wider zones under variable, and especially oscillating, proton flux.

2020