Diversity in times of adversity: probabilistic strategies in microbial survival games

Population diversification strategies are ubiquitous among microbes, encompassing random phase-variation (RPV) of pathogenic bacteria, viral latency as observed in some bacteriophage and HIV, and the non-genetic diversity of bacterial stress responses. Precise conditions under which these diversification strategies confer an advantage have not been well defined. We develop a model of population growth conditioned on dynamical environmental and cellular states. Transitions among cellular states, in turn, may be biased by possibly noisy readings of the environment from cellular sensors. For various types of environmental dynamics and cellular sensor capability, we apply game-theoretic analysis to derive the evolutionarily stable strategy (ESS) for an organism and determine when that strategy is diversification. We find that: (1) RPV, effecting a sort of Parrondo paradox wherein random alternations between losing strategies produce a winning strategy, is selected when transitions between different selective environments cannot be sensed, (2) optimal RPV cell switching rates are a function of environmental lifecycle asymmetries and environmental autocorrelation, (3) probabilistic diversification upon entering a new environment is selected when sensors can detect environmental transitions but have poor precision in identifying new environments, and (4) in the presence of excess additive noise, low-pass filtering is required for evolutionary stability. We show that even when RPV is not the ESS, it may minimize growth rate variance and the risk of extinction due to 'unlucky' environmental dynamics.     

Saccharomyces cerevisiae in the stationary phase as a model organism--characterization at cellular and proteome level

The yeast Saccharomyces cerevisiae has been used as a model organism to investigate responses to different environmental stressors. The importance of their conclusions has been expanded to human cells. The experiments were done with exponentially growing cells, which do not resemble human cells. Human and other eukaryotic cells spend the greater part of their lives in a quiescent state, known as G0 corresponding to the yeast stationary phase. Providing energy, which comes from mitochondrial respiration, is also common. Thus, in the present study S. cerevisiae was used in the stationary phase for characterization at the cellular and proteome levels. At the cellular level, optical density, cell viability, glycogen content, intracellular oxidation and cell energy metabolic activity were measured, while at the proteome level, protein profiles were analyzed using two-dimensional electrophoresis. The data obtained at both levels provide better insight into quiescence program state, which still remains poorly understood. At their base, optimal time period reflecting a stable metabolic and oxidative state of the yeast was determined. Consequently, this period is the appropriate to study changes in cell oxidant status and energy metabolic activity in response to different environmental stressors.     

Autophagy, nutrition and immunology

Turnover of cellular components in lysosomes or autophagy is an essential mechanism for cellular quality control. Added to this cleaning role, autophagy has recently been shown to participate in the dynamic interaction of cells with the surrounding environment by acting as a point of integration of extracellular cues. In this review, we focus on the relationship between autophagy and two types of environmental factors: nutrients and pathogens. We describe their direct effect on autophagy and discuss how the autophagic reaction to these stimuli allows cells to accommodate the requirements of the cellular response to stress, including those specific to the immune responses.     

Energetic adaptations: Metabolic control of endocytic membrane traffic

Endocytic membrane traffic controls the access of myriad cell surface proteins to the extracellular milieu, and thus gates nutrient uptake, ion homeostasis, signaling, adhesion and migration. Coordination of the regulation of endocytic membrane traffic with a cell's metabolic needs represents an important facet of maintenance of homeostasis under variable conditions of nutrient availability and metabolic demand. Many studies have revealed intimate regulation of endocytic membrane traffic by metabolic cues, from the specific control of certain receptors or transporters, to broader adaptation or remodeling of the endocytic membrane network. We examine how metabolic sensors such as AMP‐activated protein kinase, mechanistic target of rapamycin complex 1 and hypoxia inducible factor 1 determine sufficiency of various metabolites, and in turn modulate cellular functions that includes control of endocytic membrane traffic. We also examine how certain metabolites can directly control endocytic traffic proteins, such as the regulation of specific protein glycosylation by limiting levels of uridine diphosphate N‐acetylglucosamine (UDP‐GlcNAc) produced by the hexosamine biosynthetic pathway. From these ideas emerge a growing appreciation that endocytic membrane traffic is orchestrated by many intrinsic signals derived from cell metabolism, allowing alignment of the functions of cell surface proteins with cellular metabolic requirements. Endocytic membrane traffic determines how cells interact with their environment, thus defining many aspects of nutrient uptake and energy consumption. We examine how intrinsic signals that reflect metabolic status of a cell regulate endocytic traffic of specific proteins, and, in some cases, exert broad control of endocytic membrane traffic phenomena. Hence, endocytic traffic is versatile and adaptable and can be modulated to meet the changing metabolic requirements of a cell.     

Metabolic rate does not scale with body mass in cultured mammalian cells

The allometric scaling of metabolic rate with organism body mass can be partially accounted for by differences in cellular metabolic rates. For example, hepatocytes isolated from horses consume almost 10-fold less oxygen per unit time as mouse hepatocytes [Porter and Brand, Am J Physiol Regul Integr Comp Physiol 269: R226-R228, 1995]. This could reflect a genetically programmed, species-specific, intrinsic metabolic rate set point, or simply the adaptation of individual cells to their particular in situ environment (i.e., within the organism). We studied cultured cell lines derived from 10 mammalian species with donor body masses ranging from 5 to 600,000 g to determine whether cells propagated in an identical environment (media) exhibited metabolic rate scaling. Neither metabolic rate nor the maximal activities of key enzymes of oxidative or anaerobic metabolism scaled significantly with donor body mass in cultured cells, indicating the absence of intrinsic, species-specific, cellular metabolic rate set points. Furthermore, we suggest that changes in the metabolic rates of isolated cells probably occur within 24 h and involve a reduction of cellular metabolism toward values observed in lower metabolic rate organisms. The rate of oxygen delivery has been proposed to limit cellular metabolic rates in larger organisms. To examine the effect of oxygen on steady-state cellular respiration rates, we grew cells under a variety of physiologically relevant oxygen regimens. Long-term exposure to higher medium oxygen levels increased respiration rates of all cells, consistent with the hypothesis that higher rates of oxygen delivery in smaller mammals might increase cellular metabolic rates.     

A computational model to describe the collagen orientation in statically cultured engineered tissues

Collagen provides cardiovascular tissues with the ability to withstand haemodynamic loads. A similar network is essential to obtain in tissue-engineered (TE) samples of the same nature. Yet, the mechanism of collagen orientation is not fully understood. Typically collagen remodelling is linked to mechanical loading. However, TE constructs also show an oriented collagen network when developed under static culture. Experiments under these conditions also indicate that the tissue gradually compacts due to contractile stresses developed in the α-actin fibres of the cells. Therefore, it is hypothesised that cellular contractile stresses are responsible for collagen orientation. A model describing the cellular α-actin turnover and the stresses developed by them is integrated in a structural constitutive model describing the mechanical behaviour of collagen fibres. Results show that the model can successfully capture the sample compaction, tissue stress generation and its heterogeneous collagen arrangement.     

Local actin dynamics couple speed and persistence in a cellular Potts model of cell migration

Cell migration is astoundingly diverse. Molecular signatures, cell-cell interactions, and environmental structures each play their part in shaping cell motion, yielding numerous morphologies and migration modes. Nevertheless, in recent years, a simple unifying law was found to describe cell migration across many different cell types and contexts: faster cells turn less frequently. This universal coupling between speed and persistence (UCSP) was explained by retrograde actin flow from front to back, but it remains unclear how this mechanism generalizes to cells with complex shapes and cells migrating in structured environments, which may not have a well-defined front-to-back orientation. Here, we present an in-depth characterization of an existing cellular Potts model, in which cells polarize dynamically from a combination of local actin dynamics (stimulating protrusions) and global membrane tension along the perimeter (inhibiting protrusions). We first show that the UCSP emerges spontaneously in this model through a cross talk of intracellular mechanisms, cell shape, and environmental constraints, resembling the dynamic nature of cell migration in vivo. Importantly, we find that local protrusion dynamics suffice to reproduce the UCSP—even in cases in which no clear global, front-to-back polarity exists. We then harness the spatial nature of the cellular Potts model to show how cell shape dynamics limit both the speed and persistence a cell can reach and how a rigid environment such as the skin can restrict cell motility even further. Our results broaden the range of potential mechanisms underlying the speed-persistence coupling that has emerged as a fundamental property of migrating cells.     

Tissue-engineered bone regeneration

Bone lesions above a critical size become scarred rather than regenerated, leading to nonunion. We have attempted to obtain a greater degree of regeneration by using a resorbable scaffold with regeneration-competent cells to recreate an embryonic environment in injured adult tissues, and thus improve clinical outcome. We have used a combination of a coral scaffold with in vitro-expanded marrow stromal cells (MSC) to increase osteogenesis more than that obtained with the scaffold alone or the scaffold plus fresh bone marrow. The efficiency of the various combinations was assessed in a large segmental defect model in sheep. The tissue-engineered artificial bone underwent morphogenesis leading to complete recorticalization and the formation of a medullary canal with mature lamellar cortical bone in the most favorable cases. Clinical union never occurred when the defects were left empty or filled with the scaffold alone. In contrast, clinical union was obtained in three out of seven operated limbs when the defects were filled with the tissue-engineered bone.     

Mammalian target of rapamycin (mTOR): a central regulator of male fertility?

Mammalian target of rapamycin (mTOR) is a central regulator of cellular metabolic phenotype and is involved in virtually all aspects of cellular function. It integrates not only nutrient and energy-sensing pathways but also actin cytoskeleton organization, in response to environmental cues including growth factors and cellular energy levels. These events are pivotal for spermatogenesis and determine the reproductive potential of males. Yet, the molecular mechanisms by which mTOR signaling acts in male reproductive system remain a matter of debate. Here, we review the current knowledge on physiological and molecular events mediated by mTOR in testis and testicular cells. In recent years, mTOR inhibition has been explored as a prime strategy to develop novel therapeutic approaches to treat cancer, cardiovascular disease, autoimmunity, and metabolic disorders. However, the physiological consequences of mTOR dysregulation and inhibition to male reproductive potential are still not fully understood. Compelling evidence suggests that mTOR is an arising regulator of male fertility and better understanding of this atypical protein kinase coordinated action in testis will provide insightful information concerning its biological significance in other tissues/organs. We also discuss why a new generation of mTOR inhibitors aiming to be used in clinical practice may also need to include an integrative view on the effects in male reproductive system.     

Integrating gene expression and metabolic profiles

Recent advances in high throughput technologies have generated an abundance of biological information, such as gene expression, protein-protein interaction, and metabolic data. These various types of data capture different aspects of the cellular response to environmental factors. Integrating data from different measurements enhances the ability of modeling frameworks to predict cellular function more accurately and can lead to a more coherent reconstruction of the underlying regulatory network structure. Different techniques, newly developed and borrowed, have been applied for the purpose of extracting this information from experimental data. In this study, we developed a framework to integrate metabolic and gene expression profiles for a hepatocellular system. Specifically, we applied genetic algorithm and partial least square analysis to identify important genes relevant to a specific cellular function. We identified genes 1) whose expression levels quantitatively predict a metabolic function and 2) that play a part in regulating a hepatocellular function and reconstructed their role in the metabolic network. The framework 1) preprocesses the gene expression data using statistical techniques, 2) selects genes using a genetic algorithm and couples them to a partial least squares analysis to predict cellular function, and 3) reconstructs, with the assistance of a literature search, the pathways that regulate cellular function, namely intracellular triglyceride and urea synthesis. This provides a framework for identifying cellular pathways that are active as a function of the environment and in turn helps to uncover the interplay between gene and metabolic networks.     

Habitat selection in a variable environment

A Monte Carlo simulation scheme was utilized to determine optimal strategies of habitat utilization in a variable environment. The model allows for differences in quality among habitats at any one time and for varying levels of environmental variance and autocorrelation. When habitats are on the average equal in quality, tracking of temporal fluctuations in environment through variable habitat selection is universally advantageous with the gain in fitness limited by environmental variance, autocorrelation, and number of available habitats. Average differences in quality among habitats will restrict the advantage of variable habitat utilization (over complete usage of the average better habitat) to cases of high environmental autocorrelation or high ratios of enviromental variance to mean habitat separation. Extending an earlier prediction of Levins (1965), the average heterozygosity per individual in a natural population should increase with increasing environmental variance.     

Distribution of Tight Junction Proteins in Adult Human Salivary Glands

Tight junctions (TJs) are an essential structure of fluid-secreting cells, such as those in salivary glands. Three major families of integral membrane proteins have been identified as components of the TJ: claudins, occludin, and junctional adhesion molecules (JAMs), plus the cytosolic protein zonula occludens (ZO). We have been working to develop an orally implantable artificial salivary gland that would be suitable for treating patients lacking salivary parenchymal tissue. To date, little is known about the distribution of TJ proteins in adult human salivary cells and thus what key molecular components might be desirable for the cellular component of an artificial salivary gland device. Therefore, the aim of this study was to determine the distribution of TJ proteins in human salivary glands. Salivary gland samples were obtained from 10 patients. Frozen and formalin-fixed paraffin-embedded sections were stained using IHC methods. Claudin-1 was expressed in ductal, endothelial, and ∼25% of serous cells. Claudins-2, -3, and -4 and JAM-A were expressed in both ductal and acinar cells, whereas claudin-5 was expressed only in endothelial cells. Occludin and ZO-1 were expressed in acinar, ductal, and endothelial cells. These results provide new information on TJ proteins in two major human salivary glands and should serve as a reference for future studies to assess the presence of appropriate TJ proteins in a tissue-engineered human salivary gland. (J Histochem Cytochem 56:1093–1098, 2008)     

How Biochemical Constraints of Cellular Growth Shape Evolutionary Adaptations in Metabolism

Evolutionary adaptations in metabolic networks are fundamental to evolution of microbial growth. Studies on unneeded-protein synthesis indicate reductions in fitness upon nonfunctional protein synthesis, showing that cell growth is limited by constraints acting on cellular protein content. Here, we present a theory for optimal metabolic enzyme activity when cells are selected for maximal growth rate given such growth-limiting biochemical constraints. We show how optimal enzyme levels can be understood to result from an enzyme benefit minus cost optimization. The constraints we consider originate from different biochemical aspects of microbial growth, such as competition for limiting amounts of ribosomes or RNA polymerases, or limitations in available energy. Enzyme benefit is related to its kinetics and its importance for fitness, while enzyme cost expresses to what extent resource consumption reduces fitness through constraint-induced reductions of other enzyme levels. A metabolic fitness landscape is introduced to define the fitness potential of an enzyme. This concept is related to the selection coefficient of the enzyme and can be expressed in terms of its fitness benefit and cost.     

A Multiscale Agent-Based Model for the Investigation of <Emphasis Type="Italic">E. coli</Emphasis> K12 Metabolic Response During Biofilm Formation

Bacterial biofilm formation is an organized collective response to biochemical cues that enables bacterial colonies to persist and withstand environmental insults. We developed a multiscale agent-based model that characterizes the intracellular, extracellular, and cellular scale interactions that modulate Escherichia coli MG1655 biofilm formation. Each bacterium’s intracellular response and cellular state were represented as an outcome of interactions with the environment and neighboring bacteria. In the intracellular model, environment-driven gene expression and metabolism were captured using statistical regression and Michaelis–Menten kinetics, respectively. In the cellular model, growth, death, and type IV pili- and flagella-dependent movement were based on the bacteria’s intracellular state. We implemented the extracellular model as a three-dimensional diffusion model used to describe glucose, oxygen, and autoinducer 2 gradients within the biofilm and bulk fluid. We validated the model by comparing simulation results to empirical quantitative biofilm profiles, gene expression, and metabolic concentrations. Using the model, we characterized and compared the temporal metabolic and gene expression profiles of sessile versus planktonic bacterial populations during biofilm formation and investigated correlations between gene expression and biofilm-associated metabolites and cellular scale phenotypes. Based on our in silico studies, planktonic bacteria had higher metabolite concentrations in the glycolysis and citric acid cycle pathways, with higher gene expression levels in flagella and lipopolysaccharide-associated genes. Conversely, sessile bacteria had higher metabolite concentrations in the autoinducer 2 pathway, with type IV pili, autoinducer 2 export, and cellular respiration genes upregulated in comparison with planktonic bacteria. Having demonstrated results consistent with in vitro static culture biofilm systems, our model enables examination of molecular phenomena within biofilms that are experimentally inaccessible and provides a framework for future exploration of how hypothesized molecular mechanisms impact bulk community behavior.     

Absence of metabolic rate allometry in an ex vivo model of mammalian skeletal muscle

Within mammalian species, standard metabolic rate (SMR) increases disproportionately with body mass (Mb), such that the mass-specific SMR correlates negatively with Mb. This phenomenon can be explained in part by reduced cellular metabolic rates in larger species. To better understand the cause(s) of this cellular metabolic rate allometry we have used an ex vivo approach to isolate and identify potential contributors. Skeletal myoblasts from mammalian species ranging inMb from 30 g to over 300,000 g were isolated and differentiated into myotubes in vitro. Oxygen consumption rates, citrate synthase (CS) activity, and lactate dehydrogenase (LDH) activity were measured in myotubes under standardized conditions. No correlation of any of these parameters was observedwith speciesMb, suggesting that there is no genetic contribution to between-species differences in cellular metabolic rates. Myotubes were incubated in serum from species ranging from 30 g to 400,000 g to determine whether between-species differences in the levels of metabolically important hormones might produce allometric trends in the cultured cells. However, there was no observed effect of serum donor Mb on any of the metabolic characteristicsmeasured. Thus, there is no evidence for a relationship between skeletal muscle oxidative metabolism and Mb in an ex vivo model.     

Phenotype Switching and Mutations in Random Environments

Cell populations can benefit from changing phenotype when the environment changes. One mechanism for generating these changes is stochastic phenotype switching, whereby cells switch stochastically from one phenotype to another according to genetically determined rates, irrespective of the current environment, with the matching of phenotype to environment then determined by selective pressure. This mechanism has been observed in numerous contexts, but identifying the precise connection between switching rates and environmental changes remains an open problem. Here, we introduce a simple model to study the evolution of phenotype switching in a finite population subject to random environmental shocks. We compare the successes of competing genotypes with different switching rates, and analyze how the optimal switching rates depend on the frequency of environmental changes. If environmental changes are as rare as mutations, then the optimal switching rates mimic the rates of environmental changes. If the environment changes more frequently, then the optimal genotype either maximally favors fitness in the more common environment or has the maximal switching rate to each phenotype. Our results also explain why the optimum is relatively insensitive to fitness in each environment.     

A model of physiological adaptation in unicellular algae

A simple growth model for unicellular algae is used to show that environmentally induced changes in cellular composition can be explained in terms of controlled adjustments acting to maximize the specific growth rate. The model is based on a division of cellular carbon into four distinct compartments: carbon associated with the photosynthetic apparatus, carbon associated with those components engaged in macromolecular synthesis, carbon associated with structural components and stored carbon. Flows of material between compartments, and between cells and their environment, are defined in terms of the environmental conditions and the distribution of carbon amongst compartments. Given that growth is balanced under a specific set of environmental conditions, there exists a unique, optimal allocation of carbon for which the rate of growth is maximal. Changes in this optimal allocation of material induced by changes in light intensity, nutrient availability or temperature are qualitatively similar to compositional changes observed in a wide variety of algal species. Empirical estimates for each of the model parameters are derived and used to show that reasonable quantitative agreement between observed and predicted behaviour is attainable. The model and parameter set are also used to illustrate the influence of cell size on growth rate. Under a given set of environmental conditions, the function relating cell size to growth rate has a single maximum. The size at which growth rate is maximal varies inversely with light intensity and directly with nutrient availability and temperature. Such behaviour is consistent with some empirical observations on the influence of environmental factors on the size distribution of natural phytoplankton communities.