Limitations of field-theory simulation for exploring phase separation: The role of repulsion in a lattice protein model.

Field-theory simulation by the complex Langevin method offers an alternative to conventional sampling techniques for exploring the forces driving biomolecular liquid-liquid phase separation. Such simulations have recently been used to study several polyampholyte systems. Here, we formulate a field theory corresponding to the hydrophobic/polar (HP) lattice protein model, with finite same-site repulsion and nearest-neighbor attraction between HH bead pairs. By direct comparison with particle-based Monte Carlo simulations, we show that complex Langevin sampling of the field theory reproduces the thermodynamic properties of the HP model only if the same-site repulsion is not too strong. Unfortunately, the repulsion has to be taken weaker than what is needed to prevent condensed droplets from assuming an artificially compact shape. Analysis of a minimal and analytically solvable toy model hints that the sampling problems caused by repulsive interaction may stem from loss of ergodicity.

Cholesterol binding to the sterol-sensing region of Niemann Pick C1 protein confines dynamics of its N-terminal domain.

Lysosomal accumulation of cholesterol is a hallmark of Niemann Pick type C (NPC) disease caused by mutations primarily in the lysosomal membrane protein NPC1. NPC1 contains a transmembrane sterol-sensing domain (SSD), which is supposed to regulate protein activity upon cholesterol binding, but the mechanisms underlying this process are poorly understood. Using atomistic simulations, we show that in the absence of cholesterol in the SSD, the luminal domains of NPC1 are highly dynamic, resulting in the disengagement of the NTD from the rest of the protein. The disengaged NPC1 adopts a flexed conformation that approaches the lipid bilayer, and could represent a conformational state primed to receive a sterol molecule from the soluble lysosomal cholesterol carrier NPC2. The binding of cholesterol to the SSD of NPC1 allosterically suppresses the conformational dynamics of the luminal domains resulting in an upright NTD conformation. The presence of an additional 20% cholesterol in the membrane has negligible impact on this process. The additional presence of an NTD-bound cholesterol suppresses the flexing of the NTD. We propose that cholesterol acts as an allosteric effector, and the modulation of NTD dynamics by the SSD-bound cholesterol constitutes an allosteric feedback mechanism in NPC1 that controls cholesterol abundance in the lysosomal membrane.

A continuous growth model for plant tissue.

Morphogenesis in plants and animals involves large irreversible deformations. In plants, the response of the cell wall material to internal and external forces is determined by its mechanical properties. An appropriate model for plant tissue growth must include key features such as anisotropic and heterogeneous elasticity and cell dependent evaluation of mechanical variables such as turgor pressure, stress and strain. In addition, a growth model needs to cope with cell divisions as a necessary part of the growth process. Here we develop such a growth model, which is capable of employing not only mechanical signals but also morphogen signals for regulating growth. The model is based on a continuous equation for updating the resting configuration of the tissue. Simultaneously, material properties can be updated at a different time scale. We test the stability of our model by measuring convergence of growth results for a tissue under the same mechanical and material conditions but with different spatial discretization. The model is able to maintain a strain field in the tissue during re-meshing, which is of particular importance for modeling cell division. We confirm the accuracy of our estimations in two and three-dimensional simulations, and show that residual stresses are less prominent if strain or stress is included as input signal to growth. The approach results in a model implementation that can be used to compare different growth hypotheses, while keeping residual stresses and other mechanical variables updated and available for feeding back to the growth and material properties.

Stress and strain provide positional and directional cues in development.

The morphogenesis of organs necessarily involves mechanical interactions and changes in mechanical properties of a tissue. A long standing question is how such changes are directed on a cellular scale while being coordinated at a tissular scale. Growing evidence suggests that mechanical cues are participating in the control of growth and morphogenesis during development. We introduce a mechanical model that represents the deposition of cellulose fibers in primary plant walls. In the model both the degree of material anisotropy and the anisotropy direction are regulated by stress anisotropy. We show that the finite element shell model and the simpler triangular biquadratic springs approach provide equally adequate descriptions of cell mechanics in tissue pressure simulations of the epidermis. In a growing organ, where circumferentially organized fibers act as a main controller of longitudinal growth, we show that the fiber direction can be correlated with both the maximal stress direction and the direction orthogonal to the maximal strain direction. However, when dynamic updates of the fiber direction are introduced, the mechanical stress provides a robust directional cue for the circumferential organization of the fibers, whereas the orthogonal to maximal strain model leads to an unstable situation where the fibers reorient longitudinally. Our investigation of the more complex shape and growth patterns in the shoot apical meristem where new organs are initiated shows that a stress based feedback on fiber directions is capable of reproducing the main features of in vivo cellulose fiber directions, deformations and material properties in different regions of the shoot. In particular, we show that this purely mechanical model can create radially distinct regions such that cells expand slowly and isotropically in the central zone while cells at the periphery expand more quickly and in the radial direction, which is a well established growth pattern in the meristem.

The MOSS Physcomitrella patens reproductive organ development is highly organized, affected by the two SHI/STY genes and by the level of active auxin in the SHI/STY expression domain.

In order to establish a reference for analysis of the function of auxin and the auxin biosynthesis regulators SHORT INTERNODE/STYLISH (SHI/STY) during Physcomitrella patens reproductive development, we have described male (antheridial) and female(archegonial) development in detail, including temporal and positional information of organ initiation. This has allowed us to define discrete stages of organ morphogenesis and to show that reproductive organ development in P. patens is highly organized and that organ phyllotaxis differs between vegetative and reproductive development. Using the PpSHI1 and PpSHI2 reporter and knockout lines, the auxin reporters GmGH3(pro):GUS and PpPINA(pro):GFP-GUS, and the auxin-conjugating transgene PpSHI2(pro):IAAL, we could show that the PpSHI genes, and by inference also auxin, play important roles for reproductive organ development in moss. The PpSHI genes are required for the apical opening of the reproductive organs, the final differentiation of the egg cell, and the progression of canal cells into a cell death program. The apical cells of the archegonium, the canal cells, and the egg cell are also sites of auxin responsiveness and are affected by reduced levels of active auxin, suggesting that auxin mediates PpSHI function in the reproductive organs.

Thermodynamic Investigation of the Mechanism of Heat Production During Membrane Depolarization.

When an action potential passes through a neuron, heat is first produced and then reabsorbed by the neuronal membrane, resulting in a small measurable temperature spike. Here, we describe the thermodynamics and molecular features of the heat production using a coarse-grained molecular dynamics approach. We study a simple unicomponent lipid bilayer membrane surrounded by physiological salt solution with and without an external electric field, which represents an imbalanced charge across the membrane. We show that the temperature increases significantly upon removal of the electric field under constant pressure conditions. The potential energy converted to heat is initially stored mainly in the imbalanced ion distribution across the membrane and the elastic energy of the membrane has only a minor role to play. We demonstrate that the mechanism of heat production involves interaction between ions as well as lipid headgroup dipoles while the interactions between polar water molecules and lipid headgroup dipoles absorbs a considerable portion of such produced heat upon removal of the electric field. Our data provide novel thermodynamic insights into the molecular processes governing membrane reorganization upon discharging of lipid membranes and insight into energy metabolism in nerves.

Anisotropic growth is achieved through the additive mechanical effect of material anisotropy and elastic asymmetry.

Fast directional growth is a necessity for the young seedling; after germination, it needs to quickly penetrate the soil to begin its autotrophic life. In most dicot plants, this rapid escape is due to the anisotropic elongation of the hypocotyl, the columnar organ between the root and the shoot meristems. Anisotropic growth is common in plant organs and is canonically attributed to cell wall anisotropy produced by oriented cellulose fibers. Recently, a mechanism based on asymmetric pectin-based cell wall elasticity has been proposed. Here we present a harmonizing model for anisotropic growth control in the dark-grown Arabidopsis thaliana hypocotyl: basic anisotropic information is provided by cellulose orientation) and additive anisotropic information is provided by pectin-based elastic asymmetry in the epidermis. We quantitatively show that hypocotyl elongation is anisotropic starting at germination. We present experimental evidence for pectin biochemical differences and wall mechanics providing important growth regulation in the hypocotyl. Lastly, our in silico modelling experiments indicate an additive collaboration between pectin biochemistry and cellulose orientation in promoting anisotropic growth.

A Model Analysis of Mechanisms for Radial Microtubular Patterns at Root Hair Initiation Sites.

Plant cells have two main modes of growth generating anisotropic structures. Diffuse growth where whole cell walls extend in specific directions, guided by anisotropically positioned cellulose fibers, and tip growth, with inhomogeneous addition of new cell wall material at the tip of the structure. Cells are known to regulate these processes via molecular signals and the cytoskeleton. Mechanical stress has been proposed to provide an input to the positioning of the cellulose fibers via cortical microtubules in diffuse growth. In particular, a stress feedback model predicts a circumferential pattern of fibers surrounding apical tissues and growing primordia, guided by the anisotropic curvature in such tissues. In contrast, during the initiation of tip growing root hairs, a star-like radial pattern has recently been observed. Here, we use detailed finite element models to analyze how a change in mechanical properties at the root hair initiation site can lead to star-like stress patterns in order to understand whether a stress-based feedback model can also explain the microtubule patterns seen during root hair initiation. We show that two independent mechanisms, individually or combined, can be sufficient to generate radial patterns. In the first, new material is added locally at the position of the root hair. In the second, increased tension in the initiation area provides a mechanism. Finally, we describe how a molecular model of Rho-of-plant (ROP) GTPases activation driven by auxin can position a patch of activated ROP protein basally along a 2D root epidermal cell plasma membrane, paving the way for models where mechanical and molecular mechanisms cooperate in the initial placement and outgrowth of root hairs.

Auxin Acts through MONOPTEROS to Regulate Plant Cell Polarity and Pattern Phyllotaxis.

The periodic formation of plant organs such as leaves and flowers gives rise to intricate patterns that have fascinated biologists and mathematicians alike for hundreds of years [1]. The plant hormone auxin plays a central role in establishing these patterns by promoting organ formation at sites where it accumulates due to its polar, cell-to-cell transport [2-6]. Although experimental evidence as well as modeling suggest that feedback from auxin to its transport direction may help specify phyllotactic patterns [7-12], the nature of this feedback remains unclear [13]. Here we reveal that polarization of the auxin efflux carrier PIN-FORMED 1 (PIN1) is regulated by the auxin response transcription factor MONOPTEROS (MP) [14]. We find that in the shoot, cell polarity patterns follow MP expression, which in turn follows auxin distribution patterns. By perturbing MP activity both globally and locally, we show that localized MP activity is necessary for the generation of polarity convergence patterns and that localized MP expression is sufficient to instruct PIN1 polarity directions non-cell autonomously, toward MP-expressing cells. By expressing MP in the epidermis of mp mutants, we further show that although MP activity in a single-cell layer is sufficient to promote polarity convergence patterns, MP in sub-epidermal tissues helps anchor these polarity patterns to the underlying cells. Overall, our findings reveal a patterning module in plants that determines organ position by orienting transport of the hormone auxin toward cells with high levels of MP-mediated auxin signaling. We propose that this feedback process acts broadly to generate periodic plant architectures.