Simulation methods with extended stability for stiff biochemical Kinetics

Background  

With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows.     

Climatic sensitivity of species’ vegetative and reproductive phenology in a Hawaiian montane wet forest

Understanding how tropical tree phenology (i.e., the timing and amount of seed and leaf production) responds to climate is vital for predicting how climate change may alter ecological functioning of tropical forests. We examined the effects of temperature, rainfall, and photosynthetically active radiation (PAR) on seed phenology of four dominant species and community-level leaf phenology in a montane wet forest on the island of Hawaiʻi using monthly data collected over ~ 6 years. We expected that species phenologies would be better explained by variation in temperature and PAR than rainfall because rainfall at this site is not limiting. The best-fit model for all four species included temperature, rainfall, and PAR. For three species, including two foundational species of Hawaiian forests (Acacia koa and Metrosideros polymorpha), seed production declined with increasing maximum temperatures and increased with rainfall. Relationships with PAR were the most variable across all four species. Community-level leaf litterfall decreased with minimum temperatures, increased with rainfall, and showed a peak at PAR of ~ 400 μmol/m²s⁻¹. There was considerable variation in monthly seed and leaf production not explained by climatic factors, and there was some evidence for a mediating effect of daylength. Thus, the impact of future climate change on this forest will depend on how climate change interacts with other factors such as daylength, biotic, and/or evolutionary constraints. Our results nonetheless provide insight into how climate change may affect different species in unique ways with potential consequences for shifts in species distributions and community composition.     

Aberrant Mitochondria in a Bethlem Myopathy Patient with a Homozygous Amino Acid Substitution That Destabilizes the Collagen VI α2(VI) Chain

Bethlem myopathy and Ullrich congenital muscular dystrophy (UCMD) sit at opposite ends of a clinical spectrum caused by mutations in the extracellular matrix protein collagen VI. Bethlem myopathy is relatively mild, and patients remain ambulant in adulthood while many UCMD patients lose ambulation by their teenage years and require respiratory interventions. Dominant and recessive mutations are found across the entire clinical spectrum; however, recessive Bethlem myopathy is rare, and our understanding of the molecular pathology is limited. We studied a patient with Bethlem myopathy. Electron microscopy of his muscle biopsy revealed abnormal mitochondria. We identified a homozygous COL6A2 p.D871N amino acid substitution in the C-terminal C2 A-domain. Mutant α2(VI) chains are unable to associate with α1(VI) and α3(VI) and are degraded by the proteasomal pathway. Some collagen VI is assembled, albeit more slowly than normal, and is secreted. These molecules contain the minor α2(VI) C2a splice form that has an alternative C terminus that does include the mutation. Collagen VI tetramers containing the α2(VI) C2a chain do not assemble efficiently into microfibrils and there is a severe collagen VI deficiency in the extracellular matrix. We expressed wild-type and mutant α2(VI) C2 domains in mammalian cells and showed that while wild-type C2 domains are efficiently secreted, the mutant p.D871N domain is retained in the cell. These studies shed new light on the protein domains important for intracellular and extracellular collagen VI assembly and emphasize the importance of molecular investigations for families with collagen VI disorders to ensure accurate diagnosis and genetic counseling.     

Chemokine Receptor Ccr6 Deficiency Alters Hepatic Inflammatory Cell Recruitment and Promotes Liver Inflammation and Fibrosis

Chronic liver diseases are characterized by a sustained inflammatory response in which chemokines and chemokine-receptors orchestrate inflammatory cell recruitment. In this study we investigated the role of the chemokine receptor CCR6 in acute and chronic liver injury. In the absence of liver injury Ccr6 ^-/- mice presented a higher number of hepatic macrophages and increased expression of pro-inflammatory cytokines and M1 markers Tnf-α, Il6 and Mcp1. Inflammation and cell recruitment were increased after carbon tetrachloride-induced acute liver injury in Ccr6 ^-/- mice. Moreover, chronic liver injury by carbon tetrachloride in Ccr6 ^-/- mice was associated with enhanced inflammation and fibrosis, altered macrophage recruitment, enhanced CD4⁺ cells and a reduction in Th17 (CD4⁺IL17⁺) and mature dendritic (MHCII⁺CD11c⁺) cells recruitment. Clodronate depletion of macrophages in Ccr6 ^-/- mice resulted in a reduction of hepatic pro-inflammatory and pro-fibrogenic markers in the absence and after liver injury. Finally, increased CCR6 hepatic expression in patients with alcoholic hepatitis was found to correlate with liver expression of CCL20 and severity of liver disease. In conclusion, CCR6 deficiency affects hepatic inflammatory cell recruitment resulting in the promotion of hepatic inflammation and fibrosis.     

Pattern Selection by Dynamical Biochemical Signals

The development of multicellular organisms involves cells to decide their fate upon the action of biochemical signals. This decision is often spatiotemporally coordinated such that a spatial pattern arises. The dynamics that drive pattern formation usually involve genetic nonlinear interactions and positive feedback loops. These complex dynamics may enable multiple stable patterns for the same conditions. Under these circumstances, pattern formation in a developing tissue involves a selection process: why is a certain pattern formed and not another stable one? Herein we computationally address this issue in the context of the Notch signaling pathway. We characterize a dynamical mechanism for developmental selection of a specific pattern through spatiotemporal changes of the control parameters of the dynamics, in contrast to commonly studied situations in which initial conditions and noise determine which pattern is selected among multiple stable ones. This mechanism can be understood as a path along the parameter space driven by a sequence of biochemical signals. We characterize the selection process for three different scenarios of this dynamical mechanism that can take place during development: the signal either 1) acts in all the cells at the same time, 2) acts only within a cluster of cells, or 3) propagates along the tissue. We found that key elements for pattern selection are the destabilization of the initial pattern, the subsequent exploration of other patterns determined by the spatiotemporal symmetry of the parameter changes, and the speeds of the path compared to the timescales of the pattern formation process itself. Each scenario enables the selection of different types of patterns and creates these elements in distinct ways, resulting in different features. Our approach extends the concept of selection involved in cellular decision-making, usually applied to cell-autonomous decisions, to systems that collectively make decisions through cell-to-cell interactions.     

Recombinant Candida rugosa LIP2 expression in Pichia pastoris under the control of the AOX1 promoter

The LIP2 isoenzyme gene from Candida rugosa has been completely synthesised and functionally expressed under the AOX1 promoter control in Pichia pastoris. The on-line monitoring and control of methanol, the key inducer carbon source in fed-batch cultures, has enhanced the yield product/biomass 7.8-fold and the productivity 12.8-fold compared to the best batch cultivation with the Pichia system and, 10-fold compared to the fed-batch cultivation process using the native C. rugosa strain.Nevertheless, the high ionic strength of culture broth favoured aggregation of Lip2, leading to total loss of lipolytic activity. After cultivation, a diaultrafiltration process was implemented to diminish ionic strength, allowing for the recovery of lipolytic activity in the diaultrafiltrate. The developed bioprocess resulted into a reproducible product in terms of quality and productivity.