Metabolic engineering with multi-objective optimization of kinetic models

Kinetic models have a great potential for metabolic engineering applications. They can be used for testing which genetic and regulatory modifications can increase the production of metabolites of interest, while simultaneously monitoring other key functions of the host organism. This work presents a methodology for increasing productivity in biotechnological processes exploiting dynamic models. It uses multi-objective dynamic optimization to identify the combination of targets (enzymatic modifications) and the degree of up- or down-regulation that must be performed in order to optimize a set of pre-defined performance metrics subject to process constraints. The capabilities of the approach are demonstrated on a realistic and computationally challenging application: a large-scale metabolic model of Chinese Hamster Ovary cells (CHO), which are used for antibody production in a fed-batch process. The proposed methodology manages to provide a sustained and robust growth in CHO cells, increasing productivity while simultaneously increasing biomass production, product titer, and keeping the concentrations of lactate and ammonia at low values. The approach presented here can be used for optimizing metabolic models by finding the best combination of targets and their optimal level of up/down-regulation. Furthermore, it can accommodate additional trade-offs and constraints with great flexibility.

Reproductive characteristics of characid fish species (Teleostei, Characiformes) and their relationship with body size and phylogeny

In this study, I investigated the reproductive biology of fish species from the family Characidae of the order Characiformes. I also investigated the relationship between reproductive biology and body weight and interpreted this relationship in a phylogenetic context. The results of the present study contribute to the understanding of the evolution of the reproductive strategies present in the species of this family. Most larger characid species and other characiforms exhibit a reproductive pattern that is generally characterized by a short seasonal reproductive period that lasts one to three months, between September and April. This is accompanied by total spawning, an extremely high fecundity, and, in many species, a reproductive migration. Many species with lower fecundity exhibit some form of parental care. Although reduction in body size may represent an adaptive advantage, it may also require evolutionary responses to new biological problems that arise. In terms of reproduction, smaller species have a tendency to reduce the number of oocytes that they produce. Many small characids have a reproductive pattern similar to that of larger characiforms. On the other hand they may also exhibit a range of modifications that possibly relate to the decrease in body size and the consequent reduction in fecundity. Examples of changes in the general reproductive pattern include the following: reduction in the size of mature oocytes; increase in fecundity; production of several batches of oocytes; an extended reproductive period or even continuous reproduction that allows individuals to reproduce more than once a year; high growth rates; rapid recruitment of juveniles; presence of more than one reproductive cohort that increases the sexually active population; and multiple independent development of insemination as a reproductive strategy. These changes are possibly associated with adaptive pressures that are related to the reduction in body size. In addition, such reproductive characteristics or novelties may reflect the phylogenetic history of a given species.

Melanin-based coloration covaries with ovary size in an age-specific manner in the barn owl.

While the adaptive function of black eumelanin-based coloration is relatively well known, the function of reddish-brown pheomelanin-based coloration is still unclear. Only a few studies have shown or suggested that the degree of reddish-brownness is associated with predator-prey relationships, reproductive parameters, growth rate and immunity. To gain insight into the physiological correlates of melanin-based coloration, I collected barn owl (Tyto alba) cadavers and examined the covariation between this colour trait and ovary size, an organ that increases in size before reproduction. A relationship is expected because melanin-based coloration often co-varies with sexual activity. The results showed that reddish-brown juveniles had larger ovaries than whiter juveniles particularly in individuals in poor condition and outside the breeding season, while in birds older than 2 years lightly coloured females had larger ovaries than reddish-brown conspecifics. As barn owls become less reddish-brown between the first and second year of age, the present study suggests that reddish-brown pheomelanic and whitish colorations are associated with juvenile- and adult-specific adaptations, respectively.

Dpp signaling activity requires Pentagone to scale with tissue size in the growing Drosophila wing imaginal disc.

The wing of the fruit fly, Drosophila melanogaster, with its simple, two-dimensional structure, is a model organ well suited for a systems biology approach. The wing arises from an epithelial sac referred to as the wing imaginal disc, which undergoes a phase of massive growth and concomitant patterning during larval stages. The Decapentaplegic (Dpp) morphogen plays a central role in wing formation with its ability to co-coordinately regulate patterning and growth. Here, we asked whether the Dpp signaling activity scales, i.e. expands proportionally, with the growing wing imaginal disc. Using new methods for spatial and temporal quantification of Dpp activity and its scaling properties, we found that the Dpp response scales with the size of the growing tissue. Notably, scaling is not perfect at all positions in the field and the scaling of target gene domains is ensured specifically where they define vein positions. We also found that the target gene domains are not defined at constant concentration thresholds of the downstream Dpp activity gradients P-Mad and Brinker. Most interestingly, Pentagone, an important secreted feedback regulator of the pathway, plays a central role in scaling and acts as an expander of the Dpp gradient during disc growth.

Perturbing Hele-Shaw flow with a small gap gradient

A controlled perturbation is introduced into the Saffman-Taylor flow problem by adding a gradient to the gap of a Hele-Shaw cell. The stability of the single-finger steady state was found to be strongly affected by such a perturbation. Compared with patterns in a standard Hele-Shaw cell, the single Saffman-Taylor finger was stabilized or destabilized according to the sign of the gap gradient. While a linear stability analysis shows that this perturbation should have a negligible effect on the early-stage pattern formation, the experimental data indicate that the characteristic length for the initial breakup of a flat interface has been changed by the perturbation.

On CD-systems of stateless deterministic R-automata with window size one

We study cooperating distributed systems (CD-systems) of restarting automata that are very restricted: they are deterministic, they cannot rewrite, but only delete symbols, they restart immediately after performing a delete operation, they are stateless, and they have a read/write window of size 1 only, that is, these are stateless deterministic R(1)-automata. We study the expressive power of these systems by relating the class of languages that they accept by mode =1 computations to other well-studied language classes, showing in particular that this class only contains semi-linear languages, and that it includes all rational trace languages. In addition, we investigate the closure and non-closure properties of this class of languages and some of its algorithmic properties.

On Globally Deterministic CD-Systems of Stateless R-Automata with Window Size One

It is known that cooperating distributed systems (CD-systems) of stateless deterministic restarting automata with window size 1 accept a class of semi-linear languages that properly includes all rational trace languages. Although the component automata of such a CD-system are all deterministic, in general the CD-system itself is not, as in each of its computations, the initial component and the successor components are still chosen nondeterministically. Here we study CD-systems of stateless deterministic restarting automata with window size 1 that are themselves completely deterministic. In fact, we consider two such types of CD-systems, the strictly deterministic systems and the globally deterministic systems.

Stochastic resonance in a nonlinear model of a rotating, stratified shear flow, with a simple stochastic inertia-gravity wave parameterization

We report on a numerical study of the impact of short, fast inertia-gravity waves on the large-scale, slowly-evolving flow with which they co-exist. A nonlinear quasi-geostrophic numerical model of a stratified shear flow is used to simulate, at reasonably high resolution, the evolution of a large-scale mode which grows due to baroclinic instability and equilibrates at finite amplitude. Ageostrophic inertia-gravity modes are filtered out of the model by construction, but their effects on the balanced flow are incorporated using a simple stochastic parameterization of the potential vorticity anomalies which they induce. The model simulates a rotating, two-layer annulus laboratory experiment, in which we recently observed systematic inertia-gravity wave generation by an evolving, large-scale flow. We find that the impact of the small-amplitude stochastic contribution to the potential vorticity tendency, on the model balanced flow, is generally small, as expected. In certain circumstances, however, the parameterized fast waves can exert a dominant influence. In a flow which is baroclinically-unstable to a range of zonal wavenumbers, and in which there is a close match between the growth rates of the multiple modes, the stochastic waves can strongly affect wavenumber selection. This is illustrated by a flow in which the parameterized fast modes dramatically re-partition the probability-density function for equilibrated large-scale zonal wavenumber. In a second case study, the stochastic perturbations are shown to force spontaneous wavenumber transitions in the large-scale flow, which do not occur in their absence. These phenomena are due to a stochastic resonance effect. They add to the evidence that deterministic parameterizations in general circulation models, of subgrid-scale processes such as gravity wave drag, cannot always adequately capture the full details of the nonlinear interaction.

Shallow water flow with cylindrical symmetry

A numerical scheme is presented tor the solution of the shallow water equations in a single radial coordinate. This can prove useful when testing codes for the two-dimensional shallow water equations. The scheme is applied with success to problems involving converging and diverging bores.

Carbon fractions associated with silt-size particles in surface and subsurface soil horizons

Information on the distribution and behavior of C fractions in soil particle sizes is crucial for understanding C dynamics in soil. At present little is known about the behavior of the C associated with silt-size particles. We quantified the concentrations, distribution, and enrichment of total C (TC), readily oxidizable C (ROC), hotwater- extractable C (HWC), and cold-water-extractable C (CWC) fractions in coarse (63–20-mm), medium (20–6.3-mm), and fine (6.3–2-mm) silt-size subfractions and in coarse (2000–250 mm) and fine (250–63 mm) sand and clay (<2-mm) soil fractions isolated from bulk soil (<2 mm), and 2- to 4-mm aggregate-size fraction of surface (0–25 cm) and subsurface (25–55 cm) soils under different land uses. All measured C fractions varied significantly across all soil particle-size fractions. The highest C concentrations were associated with the <20-mm soil fractions and peaked in the medium (20–6.3-mm) and fine (6.3–2-mm) silt subfractions in most treatments. Carbon enrichment ratios (ERC) revealed the dual behavior of the C fractions associated with the medium silt-size fraction, demonstrating the simultaneous enrichment of TC and ROC, and the depletion of HWC and CWC fractions. The medium silt (20–6.3-mm) subfraction was identified in this study as a zone where the associated C fractions exhibit transitory qualities. Our results show that investigating subfractions within the silt-size particle fraction provides better understanding of the behavior of C fractions in this soil fraction.

Padé approximants for the acoustical properties of rigid frame porous media with pore size distributions

Expressions for the viscosity correction function, and hence bulk complex impedance, density, compressibility, and propagation constant, are obtained for a rigid frame porous medium whose pores are prismatic with fixed cross-sectional shape, but of variable pore size distribution. The lowand high-frequency behavior of the viscosity correction function is derived for the particular case of a log-normal pore size distribution, in terms of coefficients which can, in general, be computed numerically, and are given here explicitly for the particular cases of pores of equilateral triangular, circular, and slitlike cross-section. Simple approximate formulae, based on two-point Pade´ approximants for the viscosity correction function are obtained, which avoid a requirement for numerical integration or evaluation of special functions, and their accuracy is illustrated and investigated for the three pore shapes already mentioned

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.