There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning.

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

LFA-1 nanoclusters integrate TCR stimulation strength to tune T-cell cytotoxic activity.

T-cell cytotoxic function relies on the cooperation between the highly specific but poorly adhesive T-cell receptor (TCR) and the integrin LFA-1. How LFA-1-mediated adhesion may scale with TCR stimulation strength is ill-defined. Here, we show that LFA-1 conformation activation scales with TCR stimulation to calibrate human T-cell cytotoxicity. Super-resolution microscopy analysis reveals that >1000 LFA-1 nanoclusters provide a discretized platform at the immunological synapse to translate TCR engagement and density of the LFA-1 ligand ICAM-1 into graded adhesion. Indeed, the number of high-affinity conformation LFA-1 nanoclusters increases as a function of TCR triggering strength. Blockade of LFA-1 conformational activation impairs adhesion to target cells and killing. However, it occurs at a lower TCR stimulation threshold than lytic granule exocytosis implying that it licenses, rather than directly controls, the killing decision. We conclude that the organization of LFA-1 into nanoclusters provides a calibrated system to adjust T-cell killing to the antigen stimulation strength.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

Protein overexpression can induce the elongation of cell membrane nanodomains.

In cell membranes, proteins and lipids are organized into submicrometric nanodomains of varying sizes, shapes, and compositions, performing specific functions. Despite their biological importance, the detailed morphology of these nanodomains remains unknown. Not only can they hardly be observed by conventional microscopy due to their small size, but there is no full consensus on the theoretical models to describe their structuring and their shapes. Here, we use a combination of analytical calculations and Monte Carlo simulations based upon a model coupling membrane composition and shape to show that increasing protein concentration leads to an elongation of membrane nanodomains. The results are corroborated by single-particle tracking measurements on HIV receptors, whose level of expression in the membrane of specifically designed living cells can be tuned. These findings highlight that protein abundance can modulate nanodomain shape and potentially their biological function. Beyond biomembranes, this mesopatterning mechanism is of relevance in several soft-matter systems because it relies on generic physical arguments.

Role of the membrane anchor in the regulation of Lck activity.

Theoretical work suggests that collective spatiotemporal behavior of integral membrane proteins should be modulated by boundary lipids sheathing their membrane anchors. Here, we show evidence for this prediction while investigating the mechanism for maintaining a steady amount of the active form of integral membrane protein Lck kinase (LckA) by Lck trans-autophosphorylation regulated by the phosphatase CD45. We used super-resolution microscopy, flow cytometry, and pharmacological and genetic perturbation to gain insight into the spatiotemporal context of this process. We found that LckA is generated exclusively at the plasma membrane, where CD45 maintains it in a ceaseless dynamic equilibrium with its unphosphorylated precursor. Steady LckA shows linear dependence, after an initial threshold, over a considerable range of Lck expression levels. This behavior fits a phenomenological model of trans-autophosphorylation that becomes more efficient with increasing LckA. We then challenged steady LckA formation by genetically swapping the Lck membrane anchor with structurally divergent ones, such as that of Src or the transmembrane domains of LAT, CD4, palmitoylation-defective CD4 and CD45 that were expected to drastically modify Lck boundary lipids. We observed small but significant changes in LckA generation, except for the CD45 transmembrane domain that drastically reduced LckA due to its excessive lateral proximity to CD45. Comprehensively, LckA formation and maintenance can be best explained by lipid bilayer critical density fluctuations rather than liquid-ordered phase-separated nanodomains, as previously thought, with "like/unlike" boundary lipids driving dynamical proximity and remoteness of Lck with itself and with CD45.

Boundary fluctuation dynamics of a phase-separated domain in planar geometry.

Using phase-ordering kinetics and of renormalization group theories, we derive analytically the relaxation times of the long wavelength fluctuations of a phase-separated domain boundary in the vicinity of (and below) the critical temperature in the planar Ising universality class. For a conserved order parameter, the relaxation time grows like Λ^{3} at wavelength Λ and can be expressed in terms of parameters relevant at the microscopic scale: lattice spacing, bulk diffusion coefficient of the minority phase, and temperature. These results are supported by numerical simulations of 2D Ising models, enabling in addition calculating the nonuniversal numerical prefactor. We discuss the applications of these findings to the determination of the real timescale associated with elementary Monte Carlo moves from the measurement of long wavelength relaxation times on experimental systems or molecular dynamics simulations.

Physical constraints in polymer modeling of chromatin associations with the nuclear periphery at kilobase scale.

Interactions of chromatin with the nuclear lamina imposes a radial genome distribution important for nuclear functions. How physical properties of chromatin affect these interactions is unclear. We used polymer simulations to model how physical parameters of chromatin affect its interaction with the lamina. Impact of polymer stiffness is greater than stretching on its configurations at the lamina; these are manifested as trains describing extended interactions, and loops describing desorbed regions . Conferring an attraction potential leads to persistent interaction and adsorption-desorption regimes manifested by fluctuations between trains and loops. These are modulated by polymer stiffness and stretching, with a dominant impact of stiffness on resulting structural configurations. We infer that flexible euchromatin is more prone to stochastic interactions with lamins than rigid heterochromatin characterizing constitutive LADs. Our models provide insights on the physical properties of chromatin as a polymer which affect the dynamics and patterns of interactions with the nuclear lamina.

Domain formation in bicomponent vesicles induced by composition-curvature coupling.

Lipid vesicles composed of a mixture of two types of lipids are studied by intensive Monte Carlo numerical simulations. The coupling between the local composition and the membrane shape is induced by two different spontaneous curvatures of the components. We explore the various morphologies of these biphasic vesicles coupled to the observed patterns such as nano-domains or labyrinthine mesophases. The effect of the difference in curvatures, the surface tension, and the interaction parameter between components is thoroughly explored. Our numerical results quantitatively agree with the previous analytical results obtained by Gueguen et al. [Eur. Phys. J. E 37, 76 (2014)] in the disordered (high temperature) phase. Numerical simulations allow us to explore the full parameter space, especially close to and below the critical temperature, where analytical results are not accessible. Phase diagrams are constructed and domain morphologies are quantitatively studied by computing the structure factor and the domain size distribution. This mechanism likely explains the existence of nano-domains in cell membranes as observed by super-resolution fluorescence microscopy.

Dynamical control of denaturation bubble nucleation in supercoiled DNA minicircles.

We examine the behavior of supercoiled DNA minicircles containing between 200 and 400 base-pairs, also named microDNA, in which supercoiling favors thermally assisted DNA denaturation bubbles of nanometer size and controls their lifetime. Mesoscopic modeling and accelerated dynamics simulations allow us to overcome the limitations of atomistic simulations encountered in such systems, and offer detailed insight into the thermodynamic and dynamical properties associated with the nucleation and closure mechanisms of long-lived thermally assisted denaturation bubbles which do not stem from bending- or torque-driven stress. Suitable tuning of the degree of supercoiling and size of specifically designed microDNA is observed to lead to the control of opening characteristic times in the millisecond range, and closure characteristic times ranging over well distinct timescales, from microseconds to several minutes. We discuss how our results can be seen as a dynamical bandwidth which might enhance selectivity for specific DNA binding proteins.

The conical shape of DIM lipids promotes Mycobacterium tuberculosis infection of macrophages.

Phthiocerol dimycocerosate (DIM) is a major virulence factor of the pathogen Mycobacterium tuberculosis (Mtb). While this lipid promotes the entry of Mtb into macrophages, which occurs via phagocytosis, its molecular mechanism of action is unknown. Here, we combined biophysical, cell biology, and modeling approaches to reveal the molecular mechanism of DIM action on macrophage membranes leading to the first step of Mtb infection. Matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry showed that DIM molecules are transferred from the Mtb envelope to macrophage membranes during infection. Multiscale molecular modeling and 31P-NMR experiments revealed that DIM adopts a conical shape in membranes and aggregates in the stalks formed between 2 opposing lipid bilayers. Infection of macrophages pretreated with lipids of various shapes uncovered a general role for conical lipids in promoting phagocytosis. Taken together, these results reveal how the molecular shape of a mycobacterial lipid can modulate the biological response of macrophages.

Statistical physics and mesoscopic modeling to interpret tethered particle motion experiments.

Tethered particle motion experiments are versatile single-molecule techniques enabling one to address in vitro the molecular properties of DNA and its interactions with various partners involved in genetic regulations. These techniques provide raw data such as the tracked particle amplitude of movement, from which relevant information about DNA conformations or states must be recovered. Solving this inverse problem appeals to specific theoretical tools that have been designed in the two last decades, together with the data pre-processing procedures that ought to be implemented to avoid biases inherent to these experimental techniques. These statistical tools and models are reviewed in this paper.

Dependence of DNA Persistence Length on Ionic Strength and Ion Type.

Even though the persistence length L_{P} of double-stranded DNA plays a pivotal role in cell biology and nanotechnologies, its dependence on ionic strength I lacks a consensual description. Using a high-throughput single-molecule technique and statistical physics modeling, we measure L_{P} in the presence of monovalent (Li^{+}, Na^{+}, K^{+}) and divalent (Mg^{2+}, Ca^{2+}) metallic and alkyl ammonium ions, over a large range 0.5 mM≤I≤5 M. We show that linear Debye-Hückel-type theories do not describe even part of these data. By contrast, the Netz-Orland and Trizac-Shen formulas, two approximate theories including nonlinear electrostatic effects and the finite DNA radius, fit our data with divalent and monovalent ions, respectively, over the whole I range. Furthermore, the metallic ion type does not influence L_{P}(I), in contrast to alkyl ammonium monovalent ions at high I.

A Rationale for Mesoscopic Domain Formation in Biomembranes.

Cell plasma membranes display a dramatically rich structural complexity characterized by functional sub-wavelength domains with specific lipid and protein composition. Under favorable experimental conditions, patterned morphologies can also be observed in vitro on model systems such as supported membranes or lipid vesicles. Lipid mixtures separating in liquid-ordered and liquid-disordered phases below a demixing temperature play a pivotal role in this context. Protein-protein and protein-lipid interactions also contribute to membrane shaping by promoting small domains or clusters. Such phase separations displaying characteristic length-scales falling in-between the nanoscopic, molecular scale on the one hand and the macroscopic scale on the other hand, are named mesophases in soft condensed matter physics. In this review, we propose a classification of the diverse mechanisms leading to mesophase separation in biomembranes. We distinguish between mechanisms relying upon equilibrium thermodynamics and those involving out-of-equilibrium mechanisms, notably active membrane recycling. In equilibrium, we especially focus on the many mechanisms that dwell on an up-down symmetry breaking between the upper and lower bilayer leaflets. Symmetry breaking is an ubiquitous mechanism in condensed matter physics at the heart of several important phenomena. In the present case, it can be either spontaneous (domain buckling) or explicit, i.e., due to an external cause (global or local vesicle bending properties). Whenever possible, theoretical predictions and simulation results are confronted to experiments on model systems or living cells, which enables us to identify the most realistic mechanisms from a biological perspective.

Energy required to pinch a DNA plectoneme.

DNA supercoiling plays an important role from a biological point of view. One of its consequences at the supramolecular level is the formation of DNA superhelices named plectonemes. Normally separated by a distance on the order of 10 nm, the two opposite double strands of a DNA plectoneme must be brought closer if a protein or protein complex implicated in genetic regulation is to be bound simultaneously to both strands, as if the plectoneme was locally pinched. We propose an analytic calculation of the energetic barrier, of elastic nature, required to bring closer the two loci situated on the opposed double strands. We examine how this energy barrier scales with the DNA supercoiling. For physically relevant values of elastic parameters and of supercoiling density, we show that the energy barrier is in the k_{B}T range under physiological conditions, thus demonstrating that the limiting step to loci encounter is more likely the preceding plectoneme slithering bringing the two loci side by side.

The Wiskott-Aldrich Syndrome Protein Contributes to the Assembly of the LFA-1 Nanocluster Belt at the Lytic Synapse.

T lymphocyte cytotoxicity relies on a synaptic ring of lymphocyte function-associated antigen 1 (LFA-1), which permits polarized delivery of lytic granules. How LFA-1 organization is controlled by underlying actin cytoskeleton dynamics is poorly understood. Here, we explored the contribution of the actin cytoskeleton regulator WASP to the topography of LFA-1 using a combination of microscopy modalities. We uncover that the reduced cytotoxicity of Wiskott-Aldrich syndrome patient-derived CD8+ T lymphocytes lacking WASP is associated with reduced LFA-1 activation, unstable synapse, and delayed lethal hit. At the nanometric scale, WASP constrains high-affinity LFA-1 into dense nanoclusters located in actin meshwork interstices. At the cellular scale, WASP is required for the assembly of a radial belt composed of hundreds of LFA-1 nanoclusters and for lytic granule docking within this belt. Our study unravels the nanoscale topography of LFA-1 at the lytic synapse and identifies WASP as a molecule controlling individual LFA-1 cluster density and LFA-1 nanocluster belt integrity.

How does temperature impact the conformation of single DNA molecules below melting temperature?

The double stranded DNA molecule undergoes drastic structural changes during biological processes such as transcription during which it opens locally under the action of RNA polymerases. Local spontaneous denaturation could contribute to this mechanism by promoting it. Supporting this idea, different biophysical studies have found an unexpected increase in the flexibility of DNA molecules with various sequences as a function of the temperature, which would be consistent with the formation of a growing number of locally denatured sequences. Here, we take advantage of our capacity to detect subtle changes occurring on DNA by using high throughput tethered particle motion to question the existence of bubbles in double stranded DNA under physiological salt conditions through their conformational impact on DNA molecules ranging from several hundreds to thousands of base pairs. Our results strikingly differ from previously published ones, as we do not detect any unexpected change in DNA flexibility below melting temperature. Instead, we measure a bending modulus that remains stable with temperature as expected for intact double stranded DNA.

Fluctuation tension and shape transition of vesicles: renormalisation calculations and Monte Carlo simulations.

It has been known for long that the fluctuation surface tension of membranes r, computed from the height fluctuation spectrum, is not equal to the bare surface tension σ, which is introduced in the theory either as a Lagrange multiplier to conserve the total membrane area or as an external constraint. In this work we relate these two surface tensions both analytically and numerically. They are also compared to the Laplace tension γ, and the mechanical frame tension τ. Using the Helfrich model and one-loop renormalisation calculations, we obtain, in addition to the effective bending modulus κeff, a new expression for the effective surface tension σeff = σ - εkBT/(2ap) where kBT is the thermal energy, ap the projected cut-off area, and ε = 3 or 1 according to the allowed configurations that keep either the projected area or the total area constant. Moreover we show that the crumpling transition for an infinite planar membrane occurs for σeff = 0, and also that it coincides with vanishing Laplace and frame tensions. Using extensive Monte Carlo (MC) simulations, triangulated membranes of vesicles made of N = 100-2500 vertices are simulated within the Helfrich theory. As compared to alternative numerical models, no local constraint is applied and the shape is only controlled by the constant volume, the spontaneous curvature and σ. It is shown that the numerical fluctuation surface tension r is equal to σeff both with radial MC moves (ε = 3) and with corrected MC moves locally normal to the fluctuating membrane (ε = 1). For finite vesicles of typical size R, two different regimes are defined: a tension regime for [small sigma, Greek, circumflex]eff = σeffR2/κeff > 0 and a bending one for -1 < [small sigma, Greek, circumflex]eff < 0. A shape transition from a quasi-spherical shape imposed by the large surface energy, to more deformed shapes only controlled by the bending energy, is observed numerically at [small sigma, Greek, circumflex]eff ≃ 0. We propose that the buckling transition, observed for planar supported membranes in the literature, occurs for [small sigma, Greek, circumflex]eff ≃ -1, the associated negative frame tension playing the role of a compressive force. Hence, a precise control of the value of σeff in simulations cannot but enhance our understanding of shape transitions of vesicles and cells.

Nanodomains in Biomembranes with Recycling.

Cell membranes are out of thermodynamic equilibrium notably because of membrane recycling, i.e., active exchange of material with the cytosol. We propose an analytically tractable model of biomembrane predicting the effects of recycling on the size of protein nanodomains also called protein clusters. The model includes a short-range attraction between proteins and a weaker long-range repulsion which ensures the existence of so-called cluster phases in equilibrium, where monomeric proteins coexist with finite-size domains. Our main finding is that, when taking recycling into account, the typical cluster size at steady state increases logarithmically with the recycling rate at fixed protein concentration. Using physically realistic model parameters, the predicted 2-fold increase due to recycling in living cells is most likely experimentally measurable with the help of super-resolution microscopy.

Where Biology Meets Physics--A Converging View on Membrane Microdomain Dynamics.

For several decades, the phenomenon of membrane component segregation into microdomains has been a well-known and highly debated subject, and varying concepts including the raft hypothesis, the fence-and-picket model, hydrophobic-mismatch, and specific protein-protein interactions have been offered as explanations. Here, we review the level of insight into the molecular architecture of membrane domains one is capable of obtaining through biological experimentation. Using SNARE proteins as a paradigm, comprehensive data suggest that several dozens of molecules crowd together into almost circular spots smaller than 100 nm. Such clusters are highly dynamical as they constantly capture and lose molecules. The organization has a strong influence on the functional availability of proteins and likely provides a molecular scaffold for more complex protein networks. Despite this high level of insight, fundamental open questions remain, applying not only to SNARE protein domains but more generally to all types of membrane domains. In this context, we explain the view of physical models and how they are beneficial in advancing our concept of micropatterning. While biological models generally remain qualitative and descriptive, physics aims towards making them quantitative and providing reproducible numbers, in order to discriminate between different mechanisms which have been proposed to account for experimental observations. Despite the fundamental differences in biological and physical approaches as far as cell membrane microdomains are concerned, we are able to show that convergence on common points of views is in reach.

Probing a label-free local bend in DNA by single molecule tethered particle motion.

Being capable of characterizing DNA local bending is essential to understand thoroughly many biological processes because they involve a local bending of the double helix axis, either intrinsic to the sequence or induced by the binding of proteins. Developing a method to measure DNA bend angles that does not perturb the conformation of the DNA itself or the DNA-protein complex is a challenging task. Here, we propose a joint theory-experiment high-throughput approach to rigorously measure such bend angles using the Tethered Particle Motion (TPM) technique. By carefully modeling the TPM geometry, we propose a simple formula based on a kinked Worm-Like Chain model to extract the bend angle from TPM measurements. Using constructs made of 575 base-pair DNAs with in-phase assemblies of one to seven 6A-tracts, we find that the sequence CA6CGG induces a bend angle of 19° ± 4°. Our method is successfully compared to more theoretically complex or experimentally invasive ones such as cyclization, NMR, FRET or AFM. We further apply our procedure to TPM measurements from the literature and demonstrate that the angles of bends induced by proteins, such as Integration Host Factor (IHF) can be reliably evaluated as well.