Elisa Domínguez-Hüttinger

Resistance to Water Diffusion in the Stratum Corneum Is Depth-Dependent

The stratum corneum (SC) provides a permeability barrier that limits the inflow and outflow of water. The permeability barrier is continuously and dynamically formed, maintained, and degraded along the depth, from the bottom to the top, of the SC. Naturally, its functioning and structure also change dynamically in a depth-dependent manner. While transepidermal water loss is typically used to assess the function of the SC barrier, it fails to provide any information about the dynamic mechanisms that are responsible for the depth-dependent characteristics of the permeability barrier. This paper aims to quantitatively characterize the depth-dependency of the permeability barrier using in vivo non-invasive measurement data for understanding the underlying mechanisms for barrier formation, maintenance, and degradation. As a framework to combine existing experimental data, we propose a mathematical model of the SC, consisting of multiple compartments, to explicitly address and investigate the depth-dependency of the SC permeability barrier. Using this mathematical model, we derive a measure of the water permeability barrier, i.e. resistance to water diffusion in the SC, from the measurement data on transepidermal water loss and water concentration profiles measured non-invasively by Raman spectroscopy. The derived resistance profiles effectively characterize the depth-dependency of the permeability barrier, with three distinct regions corresponding to formation, maintenance, and degradation of the barrier. Quantitative characterization of the obtained resistance profiles allows us to compare and evaluate the permeability barrier of skin with different morphology and physiology (infants vs adults, different skin sites, before and after application of oils) and elucidates differences in underlying mechanisms of processing barriers. The resistance profiles were further used to predict the spatial-temporal effects of skin treatments by in silico experiments, in terms of spatial-temporal dynamics of percutaneous water penetration.

Risk factor-dependent dynamics of atopic dermatitis: modelling multi-scale regulation of epithelium homeostasis

Epithelial tissue provides the body with its first layer of protection against harmful environmental stimuli by enacting the regulatory interplay between a physical barrier preventing the influx of external stimuli and an inflammatory response to the infiltrating stimuli. Importantly, this interdependent regulation occurs on different time scales: the tissue-level barrier permeability is regulated over the course of hours, whereas the cellular-level enzymatic reactions leading to inflammation take place within minutes. This multi-scale regulation is key to the epitheliums function and its dysfunction leads to various diseases. This paper presents a mathematical model of regulatory mechanisms in the epidermal epithelium that includes processes on two different time scales at the cellular and tissue levels. We use this model to investigate the essential regulatory interactions between epidermal barrier integrity and skin inflammation and how their dysfunction leads to atopic dermatitis (AD). Our model exhibits a structure of dual (positive and negative) control at both cellular and tissue levels. We also determined how the variation induced by well-known risk factors for AD can break the balance of the dual control. Our model analysis based on time-scale separation suggests that each risk factor leads to qualitatively different dynamic behaviours of different severity for AD, and that the coincidence of multiple risk factors dramatically increases the fragility of the epitheliums function. The proposed mathematical framework should also be applicable to other inflammatory diseases that have similar time-scale separation and control architectures.

Mathematical modeling of atopic dermatitis reveals “double-switch” mechanisms underlying 4 common disease phenotypes

Background The skin barrier acts as the first line of defense against constant exposure to biological, microbial, physical, and chemical environmental stressors. Dynamic interplay between defects in the skin barrier, dysfunctional immune responses, and environmental stressors are major factors in the development of atopic dermatitis (AD). A systems biology modeling approach can yield significant insights into these complex and dynamic processes through integration of prior biological data. Objective We sought to develop a multiscale mathematical model of AD pathogenesis that describes the dynamic interplay between the skin barrier, environmental stress, and immune dysregulation and use it to achieve a coherent mechanistic understanding of the onset, progression, and prevention of AD. Methods We mathematically investigated synergistic effects of known genetic and environmental risk factors on the dynamic onset and progression of the AD phenotype, from a mostly asymptomatic mild phenotype to a severe treatment‐resistant form. Results Our model analysis identified a “double switch,” with 2 concatenated bistable switches, as a key network motif that dictates AD pathogenesis: the first switch is responsible for the reversible onset of inflammation, and the second switch is triggered by long‐lasting or frequent activation of the first switch, causing irreversible onset of systemic TH2 sensitization and worsening of AD symptoms. Conclusions Our mathematical analysis of the bistable switch predicts that genetic risk factors decrease the threshold of environmental stressors to trigger systemic TH2 sensitization. This analysis predicts and explains 4 common clinical AD phenotypes from a mild and reversible phenotype through to severe and recalcitrant disease and provides a mechanistic explanation for clinically demonstrated preventive effects of emollient treatments against development of AD. Graphical abstract Figure. No caption available.

Gene regulatory network underlying the immortalization of epithelial cells

BackgroundTumorigenic transformation of human epithelial cells in vitro has been described experimentally as the potential result of spontaneous immortalization. This process is characterized by a series of cell–state transitions, in which normal epithelial cells acquire first a senescent state which is later surpassed to attain a mesenchymal stem–like phenotype with a potentially tumorigenic behavior. In this paper we aim to provide a system–level mechanistic explanation to the emergence of these cell types, and to the time–ordered transition patterns that are common to neoplasias of epithelial origin. To this end, we first integrate published functional and well–curated molecular data of the components and interactions that have been found to be involved in such cell states and transitions into a network of 41 molecular components. We then reduce this initial network by removing simple mediators (i.e., linear pathways), and formalize the resulting regulatory core into logical rules that govern the dynamics of each of the network components as a function of the states of its regulators.ResultsComputational dynamic analysis shows that our proposed Gene Regulatory Network model recovers exactly three attractors, each of them defined by a specific gene expression profile that corresponds to the epithelial, senescent, and mesenchymal stem–like cellular phenotypes, respectively. We show that although a mesenchymal stem–like state can be attained even under unperturbed physiological conditions, the likelihood of converging to this state is increased when pro–inflammatory conditions are simulated, providing a systems–level mechanistic explanation for the carcinogenic role of chronic inflammatory conditions observed in the clinic. We also found that the regulatory core yields an epigenetic landscape that restricts temporal patterns of progression between the steady states, such that recovered patterns resemble the time–ordered transitions observed during the spontaneous immortalization of epithelial cells, both in vivo and in vitro.ConclusionOur study strongly suggests that the in vitro tumorigenic transformation of epithelial cells, which strongly correlates with the patterns observed during the pathological progression of epithelial carcinogenesis in vivo, emerges from underlying regulatory networks involved in epithelial trans–differentiation during development.

Mathematical Modeling of Streptococcus pneumoniae Colonization, Invasive Infection and Treatment

Streptococcus pneumoniae (Sp) is a commensal bacterium that normally resides on the upper airway epithelium without causing infection. However, factors such as co-infection with influenza virus can impair the complex Sp-host interactions and the subsequent development of many life-threatening infectious and inflammatory diseases, including pneumonia, meningitis or even sepsis. With the increased threat of Sp infection due to the emergence of new antibiotic resistant Sp strains, there is an urgent need for better treatment strategies that effectively prevent progression of disease triggered by Sp infection, minimizing the use of antibiotics. The complexity of the host-pathogen interactions has left the full understanding of underlying mechanisms of Sp-triggered pathogenesis as a challenge, despite its critical importance in the identification of effective treatments. To achieve a systems-level and quantitative understanding of the complex and dynamically-changing host-Sp interactions, here we developed a mechanistic mathematical model describing dynamic interplays between Sp, immune cells, and epithelial tissues, where the host-pathogen interactions initiate. The model serves as a mathematical framework that coherently explains various in vitro and in vitro studies, to which the model parameters were fitted. Our model simulations reproduced the robust homeostatic Sp-host interaction, as well as three qualitatively different pathogenic behaviors: immunological scarring, invasive infection and their combination. Parameter sensitivity and bifurcation analyses of the model identified the processes that are responsible for qualitative transitions from healthy to such pathological behaviors. Our model also predicted that the onset of invasive infection occurs within less than 2 days from transient Sp challenges. This prediction provides arguments in favor of the use of vaccinations, since adaptive immune responses cannot be developed de novo in such a short time. We further designed optimal treatment strategies, with minimal strengths and minimal durations of antibiotics, for each of the three pathogenic behaviors distinguished by our model. The proposed mathematical framework will help to design better disease management strategies and new diagnostic markers that can be used to inform the most appropriate patient-specific treatment options.

Computational design of treatment strategies for proactive therapy on atopic dermatitis using optimal control theory

Atopic dermatitis (AD) is a common chronic skin disease characterized by recurrent skin inflammation and a weak skin barrier, and is known to be a precursor to other allergic diseases such as asthma. AD affects up to 25% of children worldwide and the incidence continues to rise. There is still uncertainty about the optimal treatment strategy in terms of choice of treatment, potency, duration and frequency. This study aims to develop a computational method to design optimal treatment strategies for the clinically recommended ‘proactive therapy’ for AD. Proactive therapy aims to prevent recurrent flares once the disease has been brought under initial control. Typically, this is done by using an anti-inflammatory treatment such as a potent topical corticosteroid intensively for a few weeks to ‘get control’, followed by intermittent weekly treatment to suppress subclinical inflammation to ‘keep control’. Using a hybrid mathematical model of AD pathogenesis that we recently proposed, we computationally derived the optimal treatment strategies for individual virtual patient cohorts, by recursively solving optimal control problems using a differential evolution algorithm. Our simulation results suggest that such an approach can inform the design of optimal individualized treatment schedules that include application of topical corticosteroids and emollients, based on the disease status of patients observed on their weekly hospital visits. We demonstrate the potential and the gaps of our approach to be applied to clinical settings. This article is part of the themed issue ‘Mathematical methods in medicine: neuroscience, cardiology and pathology’.

Medical Systems Biology

The aim of this volume is to encourage the use of systems-level methodologies to contribute to the improvement of human-health . We intend to motivate biomedical researchers to complement their current theoretical and empirical practice with up-to-date systems biology conceptual approaches. Our perspective is based on the deep understanding of the key biomolecular regulatory mechanisms that underlie health, as well as the emergence and progression of human-disease . We strongly believe that the contemporary systems biology perspective opens the door to the effective development of novel methodologies to the improvement of prevention . This requires a deeper and integrative understanding of the involved underlying systems-level mechanisms. In order to explain our proposal in a simple way, in this chapter we privilege the conceptual exposition of our chosen framework over formal considerations. The formal exposition of our proposal will be expanded and discussed later in the next chapters.