Source author record

Anne-Florence Bitbol

Anne-Florence Bitbol appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Correlations from structure and phylogeny combine constructively in the inference of protein partners from sequences

Inferring protein-protein interactions from sequences is an important task in computational biology. Recent methods based on Direct Coupling Analysis (DCA) or Mutual Information (MI) allow to find interaction partners among paralogs of two protein families. Does successful inference mainly rely on correlations from structural contacts or from phylogeny, or both? Do these two types of signal combine constructively or hinder each other? To address these questions, we generate and analyze synthetic data produced using a minimal model that allows us to control the amounts of structural constraints and phylogeny. We show that correlations from these two sources combine constructively to increase the performance of partner inference by DCA or MI. Furthermore, signal from phylogeny can rescue partner inference when signal from contacts becomes less informative, including in the realistic case where inter-protein contacts are restricted to a small subset of sites. We also demonstrate that DCA-inferred couplings between non-contact pairs of sites improve partner inference in the presence of strong phylogeny, while deteriorating it otherwise. Moreover, restricting to non-contact pairs of sites preserves inference performance in the presence of strong phylogeny. In a natural data set, as well as in realistic synthetic data based on it, we find that non-contact pairs of sites contribute positively to partner inference performance, and that restricting to them preserves performance, evidencing an important role of phylogeny.

preprint2021arXiv

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

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 to 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.

preprint2020arXiv

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

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 cristae, 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.

preprint2020arXiv

Resist or perish: fate of a microbial population subjected to a periodic presence of antimicrobial

The evolution of antimicrobial resistance can be strongly affected by variations of antimicrobial concentration. Here, we study the impact of periodic alternations of absence and presence of antimicrobial on resistance evolution in a microbial population, using a stochastic model that includes variations of both population composition and size, and fully incorporates stochastic population extinctions. We show that fast alternations of presence and absence of antimicrobial are inefficient to eradicate the microbial population and strongly favor the establishment of resistance, unless the antimicrobial increases enough the death rate. We further demonstrate that if the period of alternations is longer than a threshold value, the microbial population goes extinct upon the first addition of antimicrobial, if it is not rescued by resistance. We express the probability that the population is eradicated upon the first addition of antimicrobial, assuming rare mutations. Rescue by resistance can happen either if resistant mutants preexist, or if they appear after antimicrobial is added to the environment. Importantly, the latter case is fully prevented by perfect biostatic antimicrobials that completely stop division of sensitive microorganisms. By contrast, we show that the parameter regime where treatment is efficient is larger for biocidal drugs than for biostatic drugs. This sheds light on the respective merits of different antimicrobial modes of action.

preprint2020arXiv

Statistical physics of interacting proteins: impact of dataset size and quality assessed in synthetic sequences

Identifying protein-protein interactions is crucial for a systems-level understanding of the cell. Recently, algorithms based on inverse statistical physics, e.g. Direct Coupling Analysis (DCA), have allowed to use evolutionarily related sequences to address two conceptually related inference tasks: finding pairs of interacting proteins, and identifying pairs of residues which form contacts between interacting proteins. Here we address two underlying questions: How are the performances of both inference tasks related? How does performance depend on dataset size and the quality? To this end, we formalize both tasks using Ising models defined over stochastic block models, with individual blocks representing single proteins, and inter-block couplings protein-protein interactions; controlled synthetic sequence data are generated by Monte-Carlo simulations. We show that DCA is able to address both inference tasks accurately when sufficiently large training sets are available, and that an iterative pairing algorithm (IPA) allows to make predictions even without a training set. Noise in the training data deteriorates performance. In both tasks we find a quadratic scaling relating dataset quality and size that is consistent with noise adding in square-root fashion and signal adding linearly when increasing the dataset. This implies that it is generally good to incorporate more data even if its quality is imperfect, thereby shedding light on the empirically observed performance of DCA applied to natural protein sequences.

preprint2019arXiv

Antibody-mediated cross-linking of gut bacteria hinders the spread of antibiotic resistance

The body is home to a diverse microbiota, mainly in the gut. Resistant bacteria are selected for by antibiotic treatments, and once resistance becomes widespread in a population of hosts, antibiotics become useless. Here, we develop a multiscale model of the interaction between antibiotic use and resistance spread in a host population, focusing on an important aspect of within-host immunity. Antibodies secreted in the gut enchain bacteria upon division, yielding clonal clusters of bacteria. We demonstrate that immunity-driven bacteria clustering can hinder the spread of a novel resistant bacterial strain in a host population. We quantify this effect both in the case where resistance pre-exists and in the case where acquiring a new resistance mutation is necessary for the bacteria to spread. We further show that the reduction of spread by clustering can be countered when immune hosts are silent carriers, and are less likely to get treated, and/or have more contacts. We demonstrate the robustness of our findings to including stochastic within-host bacterial growth, a fitness cost of resistance, and its compensation. Our results highlight the importance of interactions between immunity and the spread of antibiotic resistance, and argue in the favor of vaccine-based strategies to combat antibiotic resistance.

preprint2019arXiv

Phylogenetic correlations can suffice to infer protein partners from sequences

Determining which proteins interact together is crucial to a systems-level understanding of the cell. Recently, algorithms based on Direct Coupling Analysis (DCA) pairwise maximum-entropy models have allowed to identify interaction partners among paralogous proteins from sequence data. This success of DCA at predicting protein-protein interactions could be mainly based on its known ability to identify pairs of residues that are in contact in the three-dimensional structure of protein complexes and that coevolve to remain physicochemically complementary. However, interacting proteins possess similar evolutionary histories. What is the role of purely phylogenetic correlations in the performance of DCA-based methods to infer interaction partners? To address this question, we employ controlled synthetic data that only involve phylogeny and no interactions or contacts. We find that DCA accurately identifies the pairs of synthetic sequences that share evolutionary history. While phylogenetic correlations confound the identification of contacting residues by DCA, they are thus useful to predict interacting partners among paralogs. We find that DCA performs as well as phylogenetic methods to this end, and slightly better than them with large and accurate training sets. Employing DCA or phylogenetic methods within an Iterative Pairing Algorithm (IPA) allows to predict pairs of evolutionary partners without a training set. We demonstrate the ability of these various methods to correctly predict pairings among real paralogous proteins with genome proximity but no known physical interaction, illustrating the importance of phylogenetic correlations in natural data. However, for physically interacting and strongly coevolving proteins, DCA and mutual information outperform phylogenetic methods. We discuss how to distinguish physically interacting proteins from those only sharing evolutionary history.

preprint2016arXiv

Inferring interaction partners from protein sequences

Specific protein-protein interactions are crucial in the cell, both to ensure the formation and stability of multi-protein complexes, and to enable signal transduction in various pathways. Functional interactions between proteins result in coevolution between the interaction partners, causing their sequences to be correlated. Here we exploit these correlations to accurately identify which proteins are specific interaction partners from sequence data alone. Our general approach, which employs a pairwise maximum entropy model to infer couplings between residues, has been successfully used to predict the three-dimensional structures of proteins from sequences. Thus inspired, we introduce an iterative algorithm to predict specific interaction partners from two protein families whose members are known to interact. We first assess the algorithm's performance on histidine kinases and response regulators from bacterial two-component signaling systems. We obtain a striking 0.93 true positive fraction on our complete dataset without any a priori knowledge of interaction partners, and we uncover the origin of this success. We then apply the algorithm to proteins from ATP-binding cassette (ABC) transporter complexes, and obtain accurate predictions in these systems as well. Finally, we present two metrics that accurately distinguish interacting protein families from non-interacting ones, using only sequence data.

preprint2015arXiv

Fundamental constraints on the abundances of chemotaxis proteins

Flagellated bacteria, such as Escherichia coli, perform directed motion in gradients of concentration of attractants and repellents in a process called chemotaxis. The E. coli chemotaxis signaling pathway is a model for signal transduction, but it has unique features. We demonstrate that the need for fast signaling necessitates high abundances of the proteins involved in this pathway. We show that further constraints on the abundances of chemotaxis proteins arise from the requirements of self-assembly, both of flagellar motors and of chemoreceptor arrays. All these constraints are specific to chemotaxis, and published data confirm that chemotaxis proteins tend to be more highly expressed than their homologs in other pathways. Employing a chemotaxis pathway model, we show that the gain of the pathway at the level of the response regulator CheY increases with overall chemotaxis protein abundances. This may explain why, at least in one E. coli strain, the abundance of all chemotaxis proteins is higher in media with lower nutrient content. We also demonstrate that the E. coli chemotaxis pathway is particularly robust to abundance variations of the motor protein FliM.

preprint2014arXiv

Interplay of packing and flip-flop in local bilayer deformation. How phosphatidylglycerol could rescue mitochondrial function in a cardiolipin-deficient yeast mutant

In a previous work, we have shown that a spatially localized transmembrane pH gradient, produced by acid micro-injection near the external side of cardiolipin-containing giant unilamellar vesicles, leads to the formation of tubules that retract after the dissipation of this gradient. These tubules have morphologies similar to mitochondrial cristae. The tubulation effect is due to direct phospholipid packing modification in the outer leaflet that is promoted by protonation of cardiolipin headgroups. Here we compare the case of cardiolipin-containing giant unilamellar vesicles with that of phosphatidylglycerol-containing giant unilamellar vesicles. Local acidification also promotes formation of tubules in the latter. However, compared to cardiolipin-containing giant unilamellar vesicles the tubules are longer, exhibit a visible pearling and have a much longer lifetime after acid micro-injection is stopped. We attribute these differences to an additional mechanism that increases monolayer surface imbalance, namely inward PG flip-flop promoted by the local transmembrane pH-gradient. Simulations using a fully non-linear membrane model as well as geometrical calculations are in agreement with this hypothesis. Interestingly, among yeast mutants deficient in cardiolipin biosynthesis, only the crd1-null mutant, which accumulates phosphatidylglycerol, displays significant mitochondrial activity. Our work provides a possible explanation of such a property and further emphasizes the salient role of specific lipids in mitochondrial function.

preprint2014arXiv

Quantifying the role of population subdivision in evolution on rugged fitness landscapes

Natural selection drives populations towards higher fitness, but crossing fitness valleys or plateaus may facilitate progress up a rugged fitness landscape involving epistasis. We investigate quantitatively the effect of subdividing an asexual population on the time it takes to cross a fitness valley or plateau. We focus on a generic and minimal model that includes only population subdivision into equivalent demes connected by global migration, and does not require significant size changes of the demes, environmental heterogeneity or specific geographic structure. We determine the optimal speedup of valley or plateau crossing that can be gained by subdivision, if the process is driven by the deme that crosses fastest. We show that isolated demes have to be in the sequential fixation regime for subdivision to significantly accelerate crossing. Using Markov chain theory, we obtain analytical expressions for the conditions under which optimal speedup is achieved: valley or plateau crossing by the subdivided population is then as fast as that of its fastest deme. We verify our analytical predictions through stochastic simulations. We demonstrate that subdivision can substantially accelerate the crossing of fitness valleys and plateaus in a wide range of parameters extending beyond the optimal window. We study the effect of varying the degree of subdivision of a population, and investigate the trade-off between the magnitude of the optimal speedup and the width of the parameter range over which it occurs. Our results also hold for weakly beneficial intermediate mutations. We extend our work to the case of a population connected by migration to one or several smaller islands. Our results demonstrate that subdivision with migration alone can significantly accelerate the crossing of fitness valleys and plateaus, and shed light onto the quantitative conditions necessary for this to occur.

preprint2012arXiv

Bilayer elasticity at the nanoscale: the need for new terms

Continuum elastic models that account for membrane thickness variations are especially useful in the description of nanoscale deformations due to the presence of membrane proteins with hydrophobic mismatch. We show that terms involving the gradient and the Laplacian of the area per lipid are significant and must be retained in the effective Hamiltonian of the membrane. We reanalyze recent numerical data, as well as experimental data on gramicidin channels, in light of our model. This analysis yields consistent results for the term stemming from the gradient of the area per molecule. The order of magnitude we find for the associated amplitude, namely 13-60 mN/m, is in good agreement with the 25 mN/m contribution of the interfacial tension between water and the hydrophobic part of the membrane. The presence of this term explains a systematic variation in previously published numerical data.

preprint2012arXiv

Membrane properties revealed by spatiotemporal response to a local inhomogeneity

We study theoretically the spatiotemporal response of a lipid membrane submitted to a local chemical change of its environment, taking into account the time-dependent profile of the reagent concentration due to diffusion in the solution above the membrane. We show that the effect of the evolution of the reagent concentration profile becomes negligible after some time. It then becomes possible to extract interesting properties of the membrane response to the chemical modification. We find that a local density asymmetry between the two monolayers relaxes by spreading diffusively in the whole membrane. This behavior is driven by intermonolayer friction. Moreover, we show how the ratio of the spontaneous curvature change to the equilibrium density change induced by the chemical modification can be extracted from the dynamics of the local membrane deformation. Such information cannot be obtained by analyzing the equilibrium vesicle shapes that exist in different membrane environments in light of the area-difference elasticity model.

preprint2012arXiv

Pairwise summation approximation for Casimir potentials and its limitations

We investigate the error made by the pairwise summation (PWS) approximation in three geometries where the exact formula for the Casimir interaction is known: atom-slab, slab-slab and sphere-slab configurations. For each case the interactions are calculated analytically by summing the van der Waals interactions between the two objects. We show that the PWS result is incorrect even for an infinitely thin slab in the atom-slab configuration, because of local field effects, unless the material is infinitely dilute. In the experimentally relevant case of dielectric materials, in all considered geometries the error made by the PWS approximation is much higher than the well-known value obtained for perfect reflectors in the long-range regime. This error is maximized for permittivities close to the one of Silicon.

preprint2011arXiv

Dynamical membrane curvature instability controlled by intermonolayer friction

We study a dynamical curvature instability caused by a local chemical modification of a phospholipid membrane. In our experiments, a basic solution is microinjected close to a giant unilamellar vesicle, which induces a local chemical modification of some lipids in the external monolayer of the membrane. This modification causes a local deformation of the vesicle, which then relaxes. We present a theoretical description of this instability, taking into account both the change of the equilibrium lipid density and the change of the spontaneous membrane curvature induced by the chemical modification. We show that these two types of changes of the membrane properties yield different dynamics. In contrast, it is impossible to distinguish them when studying the equilibrium shape of a vesicle subjected to a global modification. In our model, the longest relaxation timescale is related to the intermonolayer friction, which plays an important part when there is a change in the equilibrium density in one monolayer. We compare our experimental results to the predictions of our model by fitting the measured time evolution of the deformation height to the solution of our dynamical equations. We obtain good agreement between theory and experiments. Our fits enable us to estimate the intermonolayer friction coefficient, yielding values that are consistent with previous measurements.

preprint2011arXiv

Forces exerted by a correlated fluid on embedded inclusions

We investigate the forces exerted on embedded inclusions by a fluid medium with long-range correlations, described by an effective scalar field theory. Such forces are the basis for the medium-mediated Casimir-like force. To study these forces beyond thermal average, it is necessary to define them in each microstate of the medium. Two different definitions of these forces are currently used in the literature. We study the assumptions underlying them. We show that only the definition that uses the stress tensor of the medium gives the sought-after force exerted by the medium on an embedded inclusion. If a second inclusion is embedded in the medium, the thermal average of this force gives the usual Casimir-like force between the two inclusions. The other definition can be used in the different physical case of an object that interacts with the medium without being embedded in it. We show in a simple example that the two definitions yield different results for the variance of the Casimir-like force.

preprint2011arXiv

Membrane stress tensor in the presence of lipid density and composition inhomogeneities

We derive the expression of the stress tensor for one and two-component lipid membranes with density and composition inhomogeneities. We first express the membrane stress tensor as a function of the free-energy density by means of the principle of virtual work. We then apply this general result to a monolayer model which is shown to be a local version of the area-difference elasticity (ADE) model. The resulting stress tensor expression generalizes the one associated with the Helfrich model, and can be specialized to obtain the one associated with the ADE model. Our stress tensor directly gives the force exchanged through a boundary in a monolayer with density and composition inhomogeneities. Besides, it yields the force density, which is also directly obtained in covariant formalism. We apply our results to study the forces induced in a membrane by a local perturbation.

preprint2011arXiv

Universal amplitudes of the Casimir-like interactions between four types of rods in fluid membranes

The fluctuation-induced, Casimir-like interaction between two parallel rods of length L adsorbed on a fluid membrane is calculated analytically at short separations d<<L. The rods are modeled as constraints imposed on the membrane curvature along a straight line. This allows to define four types of rods, according to whether the membrane can twist along the rod and/or curve across it. For stiff constraints, all the interaction potentials between the different types of rods are attractive and proportional to L/d. Two of the four types of rods are then equivalent, which yields six universal Casimir amplitudes. Repulsion can occur between different rods for soft constraints. Numerical results obtained for all ranges of d/L show that the attraction potential reaches kT for d/L\simeq0.2. At separations smaller than d_c \approx L(L/l_p)^(1/3), where l_p is the rod persistence length, two rods with fixed ends will bend toward each other and finally come into contact because of the Casimir interaction.

preprint2010arXiv

Fluctuations of the Casimir-like force between two membrane inclusions

Although Casimir forces are inseparable from their fluctuations, little is known about these fluctuations in soft matter systems. We use the membrane stress tensor to study the fluctuations of the membrane-mediated Casimir-like force. This method enables us to recover the Casimir force between two inclusions and to calculate its variance. We show that the Casimir force is dominated by its fluctuations. Furthermore, when the distance d between the inclusions is decreased from infinity, the variance of the Casimir force decreases as -1/d^2. This distance dependence shares a common physical origin with the Casimir force itself.