Vesicle formation induced by thermal fluctuations.

The process of fission and vesicle formation depends on the geometry of the membrane that will split. For instance, a flat surface finds it difficult to form vesicles because of the lack of curved regions where to start the process. Here we show that vesicle formation can be promoted by temperature, by using a membrane phase field model with Gaussian curvature. We find a phase transition between fluctuating and vesiculation phases that depends on temperature, spontaneous curvature, and the ratio between bending and Gaussian moduli. We analysed the energy dynamical behaviour of these processes and found that the main driving ingredient is the Gaussian energy term, although the curvature energy term usually helps with the process as well. We also found that the chemical potential can be used to investigate the temperature of the system. Finally we address how temperature changes the condition for spontaneous vesiculation for all geometries, making it happen in a wider range of values of the Gaussian modulus.

Treatment with low-dose prednisone in refractory obstetric antiphospholipid syndrome: A retrospective cohort study and meta-analysis.

BACKGROUND: Glucocorticoids have been suggested as a potential therapy in refractory obstetric antiphospholipid syndrome (oAPS). Our aims were to describe a cohort of patients with oAPS treated with low-dose glucocorticoids and to perform a systematic review and meta-analysis evaluating the effects of additional glucocorticoids on the pregnancy outcomes in oAPS patients. METHODS: Retrospective study that included 11 women diagnosed with primary antiphospholipid syndrome. The meta-analysis was conducted by fitting random effects models and was checked for heterogeneity. RESULTS: All women had suffered from early pregnancy losses and two also had a history of fetal deaths. We studied 47 pregnancies that resulted in 32 abortions (68.1%) and 3 fetal deaths (6.4%). Twenty-six pregnancies were under treatment, mainly LDA and LMWH. Low-dose glucocorticoids were indicated in 13 pregnancies (always in association with LDA and LMWH). There was a decrease in pregnancy loss in those patients treated with LDA and LMWH. Treatment with glucocorticoids significantly increased the rate of successful pregnancy (38.5% abortions in treated vs 85.3% abortions in non-treated pregnancies; p=0.003). After multivariate GEE analysis, only glucocorticoids remained inversely associated with pregnancy loss (OR=0.157, (CI 0.025-0.968, p=0.046)). The meta-analysis showed that glucocorticoids tended to improve the frequency of successful pregnancy (OR= 0.509 (0.252-1.028), p=0.06). Three cases of gestational diabetes and one of preeclampsia were observed in our cohort. The meta-analysis, which mostly included studies using high-dose steroids, showed that glucocorticoids increased not only the frequency of preeclampsia and gestational diabetes, but also the rate of pre-term birth. CONCLUSIONS: The efficacy of low-dose glucocorticoids in addition to the standard therapy in patients with refractory oAPS should be confirmed in well-designed clinical trials. However, high doses of steroids significantly increase the frequency of maternal and fetal morbidities, making their use strongly inadvisable.

Coherent directed movement toward food modeled in Trichoplax, a ciliated animal lacking a nervous system.

Trichoplax adhaerens is a small, ciliated marine animal that glides on surfaces grazing upon algae, which it digests externally. It has no muscles or nervous system and only six cell types, all but two of which are embedded in its epithelium. The epithelial cells are joined by apical adherens junctions; neither tight junctions nor gap junctions are present. Monociliated epithelial cells on the lower surface propel gliding. The cilia beat regularly, but asynchronously, and transiently contact the substrate with each stroke. The animal moves in random directions in the absence of food. We show here that it exhibits chemotaxis, moving preferentially toward algae embedded in a disk of agar. We present a mathematical model to explain how coherent, directional movements could arise from the collective actions of a set of ciliated epithelial cells, each independently sensing and responding to a chemoattractant gradient. The model incorporates realistic values for viscoelastic properties of cells and produces coordinated movements and changes in body shape that resemble the actual movements of the animal. The model demonstrates that an animal can move coherently in search of food without any need for chemical signaling between cells and introduces a different approach to modeling behavior in primitive multicellular organisms.

Dipole-dipole interactions control the interfacial rheological response of cyclodextrin/surfactant solutions.

A recent surface rheological study has shown that aqueous solutions of α-cyclodextrin (αCD) with anionic surfactants (S) display a remarkable viscoelasticity at the liquid/air interface, which has not been observed in similar systems. The dilatational modulus is various orders of magnitude larger than those for the binary mixtures αCD + water and S + water. The rheological response has been qualitatively related to the bulk distribution of species, the 2 : 1 inclusion complexes (αCD2 : S) playing a fundamental role. In this work, we have developed a model that considers dipole-dipole interactions between 2 : 1 inclusion complexes ordered on the liquid/air interface. When the model is applied to the specific experimental conditions, the dependencies on concentration and temperature of the dilatational modulus and the surface tension were found to be in excellent agreement with the data, indicating clearly that dipole-dipole interactions determine and control the rheological behavior of the interface.

A model for social spreading of Covid-19: Cases of Mexico, Finland and Iceland.

The shocking severity of the Covid-19 pandemic has woken up an unprecedented interest and accelerated effort of the scientific community to model and forecast epidemic spreading to find ways to control it regionally and between regions. Here we present a model that in addition to describing the dynamics of epidemic spreading with the traditional compartmental approach takes into account the social behaviour of the population distributed over a geographical region. The region to be modelled is defined as a two-dimensional grid of cells, in which each cell is weighted with the population density. In each cell a compartmental SEIRS system of delay difference equations is used to simulate the local dynamics of the disease. The infections between cells are modelled by a network of connections, which could be terrestrial, between neighbouring cells, or long range, between cities by air, road or train traffic. In addition, since people make trips without apparent reason, noise is considered to account for them to carry contagion between two randomly chosen distant cells. Hence, there is a clear separation of the parameters related to the biological characteristics of the disease from the ones that represent the spatial spread of infections due to social behaviour. We demonstrate that these parameters provide sufficient information to trace the evolution of the pandemic in different situations. In order to show the predictive power of this kind of approach we have chosen three, in a number of ways different countries, Mexico, Finland and Iceland, in which the pandemics have followed different dynamic paths. Furthermore we find that our model seems quite capable of reproducing the path of the pandemic for months with few initial data. Unlike similar models, our model shows the emergence of multiple waves in the case when the disease becomes endemic.

Dynamics of hierarchical weighted networks of van der Pol oscillators.

We investigate the dynamics of regular fractal-like networks of hierarchically coupled van der Pol oscillators. The hierarchy is imposed in terms of the coupling strengths or link weights. We study the low frequency modes, as well as frequency and phase synchronization, in the network by a process of repeated coarse-graining of oscillator units. At any given stage of this process, we sum over the signals from the oscillator units of a clique to obtain a new oscillating unit. The frequencies and the phases for the coarse-grained oscillators are found to progressively synchronize with the number of coarse-graining steps. Furthermore, the characteristic frequency is found to decrease and finally stabilize to a value that can be tuned via the parameters of the system. We compare our numerical results with those of an approximate analytic solution and find good qualitative agreement. Our study on this idealized model shows how oscillations with a precise frequency can be obtained in systems with heterogeneous couplings. It also demonstrates the effect of imposing a hierarchy in terms of link weights instead of one that is solely topological, where the connectivity between oscillators would be the determining factor, as is usually the case.

A physico-genetic module for the polarisation of auxin efflux carriers PIN-FORMED (PIN).

Intracellular polarisation of auxin efflux carriers is crucial for understanding how auxin gradients form in plants. The polarisation dynamics of auxin efflux carriers PIN-FORMED (PIN) depends on both biomechanical forces as well as chemical, molecular and genetic factors. Biomechanical forces have shown to affect the localisation of PIN transporters to the plasma membrane. We propose a physico-genetic module of PIN polarisation that integrates biomechanical, molecular, and cellular processes as well as their non-linear interactions. The module was implemented as a discrete Boolean model and then approximated to a continuous dynamic system, in order to explore the relative contribution of the factors mediating PIN polarisation at the scale of single cell. Our models recovered qualitative behaviours that have been experimentally observed and enable us to predict that, in the context of PIN polarisation, the effects of the mechanical forces can predominate over the activity of molecular factors such as the GTPase ROP6 and the ROP-INTERACTIVE CRIB MOTIF-CONTAINING PROTEIN RIC1.

Effects of deception in social networks.

Honesty plays a crucial role in any situation where organisms exchange information or resources. Dishonesty can thus be expected to have damaging effects on social coherence if agents cannot trust the information or goods they receive. However, a distinction is often drawn between prosocial lies ('white' lies) and antisocial lying (i.e. deception for personal gain), with the former being considered much less destructive than the latter. We use an agent-based model to show that antisocial lying causes social networks to become increasingly fragmented. Antisocial dishonesty thus places strong constraints on the size and cohesion of social communities, providing a major hurdle that organisms have to overcome (e.g. by evolving counter-deception strategies) in order to evolve large, socially cohesive communities. In contrast, white lies can prove to be beneficial in smoothing the flow of interactions and facilitating a larger, more integrated network. Our results demonstrate that these group-level effects can arise as emergent properties of interactions at the dyadic level. The balance between prosocial and antisocial lies may set constraints on the structure of social networks, and hence the shape of society as a whole.

Cell patterns emerge from coupled chemical and physical fields with cell proliferation dynamics: the Arabidopsis thaliana root as a study system.

A central issue in developmental biology is to uncover the mechanisms by which stem cells maintain their capacity to regenerate, yet at the same time produce daughter cells that differentiate and attain their ultimate fate as a functional part of a tissue or an organ. In this paper we propose that, during development, cells within growing organs obtain positional information from a macroscopic physical field that is produced in space while cells are proliferating. This dynamical interaction triggers and responds to chemical and genetic processes that are specific to each biological system. We chose the root apical meristem of Arabidopsis thaliana to develop our dynamical model because this system is well studied at the molecular, genetic and cellular levels and has the key traits of multicellular stem-cell niches. We built a dynamical model that couples fundamental molecular mechanisms of the cell cycle to a tension physical field and to auxin dynamics, both of which are known to play a role in root development. We perform extensive numerical calculations that allow for quantitative comparison with experimental measurements that consider the cellular patterns at the root tip. Our model recovers, as an emergent pattern, the transition from proliferative to transition and elongation domains, characteristic of stem-cell niches in multicellular organisms. In addition, we successfully predict altered cellular patterns that are expected under various applied auxin treatments or modified physical growth conditions. Our modeling platform may be extended to explicitly consider gene regulatory networks or to treat other developmental systems.

Modelling the creation of friends and foes groups in small real social networks.

Although friendship networks have been extensively studied, few models and studies are available to understand the reciprocity of friendship and foes. Here a model is presented to explain the directed friendship and foes network formation observed in experiments of Mexican and Hungarian schools. Within the presented model, each agent has a private opinion and a public one that shares to the group. There are two kinds of interactions between agents. The first kind represent interactions with the neighbors while the other represents the attitude of an agent to the overall public available information. Links between agents evolve as a combination of the public and private information available. Friendship is defined using a fitness function according to the strength of the agent's bonds, clustering coefficient, betweenness centrality and degree. Enmity is defined as very negative links. The model allows us to reproduce the distribution of mentions for friends and foes observed in the experiments, as well as the topology of the directed networks.

Socio-economic pandemic modelling: case of Spain.

A global disaster, such as the recent Covid-19 pandemic, affects every aspect of our lives and there is a need to investigate these highly complex phenomena if one aims to diminish their impact in the health of the population, as well as their socio-economic stability. In this paper we present an attempt to understand the role of the governmental authorities and the response of the rest of the population facing such emergencies. We present a mathematical model that takes into account the epidemiological features of the pandemic and also the actions of people responding to it, focusing only on three aspects of the system, namely, the fear of catching this serious disease, the impact on the economic activities and the compliance of the people to the mitigating measures adopted by the authorities. We apply the model to the specific case of Spain, since there are accurate data available about these three features. We focused on tourism as an example of the economic activity, since this sector of economy is one of the most likely to be affected by the restrictions imposed by the authorities, and because it represents an important part of Spanish economy. The results of numerical calculations agree with the empirical data in such a way that we can acquire a better insight of the different processes at play in such a complex situation, and also in other different circumstances.

Forecasting virus outbreaks with social media data via neural ordinary differential equations.

During the Covid-19 pandemic, real-time social media data could in principle be used as an early predictor of a new epidemic wave. This possibility is examined here by employing a neural ordinary differential equation (neural ODE) trained to forecast viral outbreaks in a specific geographic region. It learns from multivariate time series of signals derived from a novel set of large online polls regarding COVID-19 symptoms. Once trained, the neural ODE can capture the dynamics of interconnected local signals and effectively estimate the number of new infections up to two months in advance. In addition, it may predict the future consequences of changes in the number of infected at a certain period, which might be related with the flow of individuals entering or exiting a region. This study provides persuasive evidence for the predictive ability of widely disseminated social media surveys for public health applications.

Modelling of chemotactic sprouting endothelial cells through an extracellular matrix.

Sprouting angiogenesis is a core biological process critical to vascular development. Its accurate simulation, relevant to multiple facets of human health, is of broad, interdisciplinary appeal. This study presents an in-silico model replicating a microfluidic assay where endothelial cells sprout into a biomimetic extracellular matrix, specifically, a large-pore, low-concentration fibrin-based porous hydrogel, influenced by chemotactic factors. We introduce a novel approach by incorporating the extracellular matrix and chemotactic factor effects into a unified term using a single parameter, primarily focusing on modelling sprouting dynamics and morphology. This continuous model naturally describes chemotactic-induced sprouting with no need for additional rules. In addition, we extended our base model to account for matrix sensing and degradation, crucial aspects of angiogenesis. We validate our model via a hybrid in-silico experimental method, comparing the model predictions with experimental results derived from the microfluidic setup. Our results underscore the intricate relationship between the extracellular matrix structure and angiogenic sprouting, proposing a promising method for predicting the influence of the extracellular matrix on angiogenesis.

From ABC genes to regulatory networks, epigenetic landscapes and flower morphogenesis: making biological sense of theoretical approaches.

The ABC model postulates that expression combinations of three classes of genes (A, B and C) specify the four floral organs at early stages of flower development. This classic model provides a solid framework to study flower development and has been the foundation for multiple studies in different plant species, as well as for new evolutionary hypotheses. Nevertheless, it has been shown that in spite of being necessary, these three gene classes are not sufficient for flower organ specification. Rather, flower organ specification depends on complex interactions of several genes, and probably other non-genetic factors. Being useful to study systems of complex interactions, mathematical and computational models have enlightened the origin of the A, B and C stereotyped and robust expression patterns and the process of early flower morphogenesis. Here, we present a brief introduction to basic modeling concepts and techniques and review the results that these models have rendered for the particular case of the Arabidopsis thaliana flower organ specification. One of the main results is the uncovering of a robust functional module that is sufficient to recover the gene configurations characterizing flower organ primordia. Another key result is that the temporal sequence with which such gene configurations are attained may be recovered only by modeling the aforementioned functional module as a noisy or stochastic system. Finally, modeling approaches enable testable predictions regarding the role of non-genetic factors (noise, mechano-elastic forces, etc.) in development. These predictions, along with some perspectives for future work, are also reviewed and discussed.

Dynamical shapes of droplets of cyclodextrin-surfactant solutions.

We present a series of experiments with droplets of aqueous cyclodextrin-surfactant solutions, in which the volume is reduced after the equilibrium spherical shape is reached. The final shape of the drop after this perturbation is found to be dependent on the concentration of inclusion complexes in the bulk of the solution. These inclusion complexes are formed by two cyclodextrin molecules and one surfactat molecule. We propose a model to describe these dynamical processes. Dipole-dipole interactions on the surface of the drop trigger a competition between water surface tension and dipole-dipole interaction energies. The results of the model reproduce the spherical and rod-like shapes found in the experiments.

Membrane rigidity regulates E. coli proliferation rates.

Combining single cell experiments, population dynamics and theoretical methods of membrane mechanics, we put forward that the rate of cell proliferation in E. coli colonies can be regulated by modifiers of the mechanical properties of the bacterial membrane. Bacterial proliferation was modelled as mediated by cell division through a membrane constriction divisome based on FtsZ, a mechanically competent protein at elastic interaction against membrane rigidity. Using membrane fluctuation spectroscopy in the single cells, we revealed either membrane stiffening when considering hydrophobic long chain fatty substances, or membrane softening if short-chained hydrophilic molecules are used. Membrane stiffeners caused hindered growth under normal division in the microbial cultures, as expected for membrane rigidification. Membrane softeners, however, altered regular cell division causing persistent microbes that abnormally grow as long filamentous cells proliferating apparently faster. We invoke the concept of effective growth rate under the assumption of a heterogeneous population structure composed by distinguishable individuals with different FtsZ-content leading the possible forms of cell proliferation, from regular division in two normal daughters to continuous growing filamentation and budding. The results settle altogether into a master plot that captures a universal scaling between membrane rigidity and the divisional instability mediated by FtsZ at the onset of membrane constriction.

On Gaussian curvature and membrane fission.

We propose a three-dimensional mathematical model to describe dynamical processes of membrane fission. The model is based on a phase field equation that includes the Gaussian curvature contribution to the bending energy. With the addition of the Gaussian curvature energy term numerical simulations agree with the predictions that tubular shapes can break down into multiple vesicles. A dispersion relation obtained with linear analysis predicts the wavelength of the instability and the number of formed vesicles. Finally, a membrane shape diagram is obtained for the different Gaussian and bending modulus, showing different shape regimes.

Behavioral role of PACAP signaling reflects its selective distribution in glutamatergic and GABAergic neuronal subpopulations.

The neuropeptide PACAP, acting as a co-transmitter, increases neuronal excitability, which may enhance anxiety and arousal associated with threat conveyed by multiple sensory modalities. The distribution of neurons expressing PACAP and its receptor, PAC1, throughout the mouse nervous system was determined, in register with expression of glutamatergic and GABAergic neuronal markers, to develop a coherent chemoanatomical picture of PACAP role in brain motor responses to sensory input. A circuit role for PACAP was tested by observing Fos activation of brain neurons after olfactory threat cue in wild-type and PACAP knockout mice. Neuronal activation and behavioral response, were blunted in PACAP knock-out mice, accompanied by sharply downregulated vesicular transporter expression in both GABAergic and glutamatergic neurons expressing PACAP and its receptor. This report signals a new perspective on the role of neuropeptide signaling in supporting excitatory and inhibitory neurotransmission in the nervous system within functionally coherent polysynaptic circuits.

Universal model for the skin colouration patterns of neotropical catfishes of the genus Pseudoplatystoma.

Fish skin colouration has been widely studied because it involves a variety of processes that are important to the broad field of the developmental biology. Mathematical modelling of fish skin patterning first predicted the existence of morphogens and helped to elucidate the mechanisms of pattern formation. The catfishes of the genus Pseudoplatystoma offer a good biological study model, since its species exhibit the most spectacular and amazing variations of colour patterns on the skin. They present labyrinths, closed loops (or cells), alternate spots and stripes, only spots and combinations of these. We have extended a well known mathematical model to study the skin of Pseudoplatystoma. The basic model is a two component, non-linear reaction diffusion system that presents a richness of bifurcations. The extended model assumes that there are two interacting cell/tissue layers in which morphogens diffuse and interact giving rise to the skin colouration pattern. We have found that by varying only two parameters we are able to accurately reproduce the distinct patterns found in all species of Pseudoplatystoma. The histological analysis of skin samples of two species of this genus, with different patterns, revealed differences on the disposition of the colouration cells that are consistent with our theoretical predictions.

Modeling the skin pattern of fishes.

Complicated patterns showing various spatial scales have been obtained in the past by coupling Turing systems in such a way that the scales of the independent systems resonate. This produces superimposed patterns with different length scales. Here we propose a model consisting of two identical reaction-diffusion systems coupled together in such a way that one of them produces a simple Turing pattern of spots or stripes, and the other traveling wave fronts that eventually become stationary. The basic idea is to assume that one of the systems becomes fixed after some time and serves as a source of morphogens for the other system. This mechanism produces patterns very similar to the pigmentation patterns observed in different species of stingrays and other fishes. The biological mechanisms that support the realization of this model are discussed.