Cephalosporin C production by immobilized Cephalosporium acremonium cells in a repeated batch tower bioreactor

The industrial production of antibiotics with filamentous fungi is usually carried out in conventional aerated and agitated tank fermentors. Highly viscous non-Newtonian broths are produced and a compromise must be found between convenient shear stress and adequate oxygen transfer. In this work, cephalosporin C production by bioparticles of immobilized cells of Cephalosporium acremonium ATCC 48272 was studied in a repeated batch tower bioreactor as an alternative to the conventional process. Also, gas-liquid oxygen transfer volumetric coefficients, k(L)a, were determined at various air flow-rates and alumina contents in the bioparticle. The bioparticles were composed of calcium alginate (2.0% w/w), alumina (<44 micra), cells, and water. A model describing the cell growth, cephalosporin C production, oxygen, glucose, and sucrose consumption was proposed. To describe the radial variation of oxygen concentration within the pellet, the reaction-diffusion model forecasting a dead core bioparticle was adopted. The k(L)a measurements with gel beads prepared with 0.0, 1.0, 1.5, and 2.0% alumina showed that a higher k(L)a value is attained with 1.5 and 2.0%. An expression relating this coefficient to particle density, liquid density, and air velocity was obtained and further utilized in the simulation of the proposed model. Batch, followed by repeated batch experiments, were accomplished by draining the spent medium, washing with saline solution, and pouring fresh medium into the bioreactor. Results showed that glucose is consumed very quickly, within 24 h, followed by sucrose consumption and cephalosporin C production. Higher productivities were attained during the second batch, as cell concentration was already high, resulting in rapid glucose consumption and an early derepression of cephalosporin C synthesizing enzymes. The model incorporated this improvement predicting higher cephalosporin C productivity. (C) 2004 Wiley Periodicals, Inc.

Studies on the respiration rate of free and immobilized cells of Cephalosporium acremonium in cephalosporin C production

Bioprocesses using filamentous fungi immobilized in inert supports present many advantages when compared to conventional free cell processes. However, assessment of the real advantages of the unconventional process demands a rigorous study of the limitations to diffusional mass transfer of the reagents, especially concerning oxygen. In this work, a comparative study was carried out on the cephalosporin C production process in defined medium containing glucose and sucrose as main carbon and energy sources, by free and immobilized cells of Cephalosporium acremonium ATCC 48272 in calcium alginate gel beads containing alumina. The effective diffusivity of oxygen through the gel beads and the effectiveness factors related to the respiration rate of the microorganism were determined experimentally. By applying Monod kinetics, the respiration kinetics parameters were experimentally determined in independent experiments in a complete production medium. The effectiveness factor experimental values presented good agreement with the theoretical values of the approximated zero-order effectiveness factor, considering the dead core model. Furthermore, experimental results obtained with immobilized cells in a 1.7-L tower bioreactor were compared with those obtained in 5-L conventional fermenter with free cells. It could be concluded that it is possible to attain rather high production rates working with relatively large diameter gel beads (ca. 2.5 mm) and sucrose consumption-based productivity was remarkably higher with immobilized cells, i.e., 0.33 gCPC/kg sucrose/h against 0.24 gCPC/kg sucrose/h in the aerated stirred tank bioreactor process. (C) 1999 John Wiley & Sons, Inc.

Description of inclusive meson spectra in e+e- annihilation reactions by a quark cascade model

We present a model to describe inclusive meson production in e+e- reactions based on a quark cascade approach whose formulation is put in terms of diffusion equations for three quark flavors (u, d, s). These equations are solved by using a formalism previously developed for the problem of the electromagnetic cascade generated in the atmosphere by cosmicray interactions. The obtained solutions are given in terms of a combination of power-law functions whose profiles are adequate to describe the characteristics observed in the inclusive spectrum of mesons.

Allee effects and population spread in patchy landscapes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

How population loss through habitat boundaries determines the dynamics of a predator-prey system

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

An introductory informatics theoretical framework, based on an integrable lattice model using quantized function algebras, for studying DNA-ionized gas interaction system

Usually we observe that Bio-physical systems or Bio-chemical systems are many a time based on nanoscale phenomenon in different host environments, which involve many particles can often not be solved explicitly. Instead a physicist, biologist or a chemist has to rely either on approximate or numerical methods. For a certain type of systems, called integrable in nature, there exist particular mathematical structures and symmetries which facilitate the exact and explicit description. Most integrable systems, we come across are low-dimensional, for instance, a one-dimensional chain of coupled atoms in DNA molecular system with a particular direction or exist as a vector in the environment. This theoretical research paper aims at bringing one of the pioneering ‘Reaction-Diffusion’ aspects of the DNA-plasma material system based on an integrable lattice model approach utilizing quantized functional algebras, to disseminate the new developments, initiate novel computational and design paradigms.

Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Convergence results for stochastic particle systems with social interaction

We consider stochastic individual-based models for social behaviour of groups of animals. In these models the trajectory of each animal is given by a stochastic differential equation with interaction. The social interaction is contained in the drift term of the SDE. We consider a global aggregation force and a short-range repulsion force. The repulsion range and strength gets rescaled with the number of animals N. We show that for N tending to infinity stochastic fluctuations disappear and a smoothed version of the empirical process converges uniformly towards the solution of a nonlinear, nonlocal partial differential equation of advection-reaction-diffusion type. The rescaling of the repulsion in the individual-based model implies that the corresponding term in the limit equation is local while the aggregation term is non-local. Moreover, we discuss the effect of a predator on the system and derive an analogous convergence result. The predator acts as an repulsive force. Different laws of motion for the predator are considered.

Efficiency limits for silicon/perovskite tandem solar cells: a theoretical model

The primary goal of this work is related to the extension of an analytic electro-optical model. It will be used to describe single-junction crystalline silicon solar cells and a silicon/perovskite tandem solar cell in the presence of light-trapping in order to calculate efficiency limits for such a device. In particular, our tandem system is composed by crystalline silicon and a perovskite structure material: metilammoniumleadtriiodide (MALI). Perovskite are among the most convenient materials for photovoltaics thanks to their reduced cost and increasing efficiencies. Solar cell efficiencies of devices using these materials increased from 3.8% in 2009 to a certified 20.1% in 2014 making this the fastest-advancing solar technology to date. Moreover, texturization increases the amount of light which can be absorbed through an active layer. Using Green’s formalism it is possible to calculate the photogeneration rate of a single-layer structure with Lambertian light trapping analytically. In this work we go further: we study the optical coupling between the two cells in our tandem system in order to calculate the photogeneration rate of the whole structure. We also model the electronic part of such a device by considering the perovskite top cell as an ideal diode and solving the drift-diffusion equation with appropriate boundary conditions for the silicon bottom cell. We have a four terminal structure, so our tandem system is totally unconstrained. Then we calculate the efficiency limits of our tandem including several recombination mechanisms such as Auger, SRH and surface recombination. We focus also on the dependence of the results on the band gap of the perovskite and we calculare an optimal band gap to optimize the tandem efficiency. The whole work has been continuously supported by a numerical validation of out analytic model against Silvaco ATLAS which solves drift-diffusion equations using a finite elements method. Our goal is to develop a simpler and cheaper, but accurate model to study such devices.

Invariant barriers to reactive front propagation in fluid flows

We present theory and experiments on the dynamics of reaction fronts in two-dimensional, vortex-dominated flows, for both time-independent and periodically driven cases. We find that the front propagation process is controlled by one-sided barriers that are either fixed in the laboratory frame (time-independent flows) or oscillate periodically (periodically driven flows). We call these barriers burning invariant manifolds (BIMs), since their role in front propagation is analogous to that of invariant manifolds in the transport and mixing of passive impurities under advection. Theoretically, the BIMs emerge from a dynamical systems approach when the advection-reaction-diffusion dynamics is recast as an ODE for front element dynamics. Experimentally, we measure the location of BIMs for several laboratory flows and confirm their role as barriers to front propagation.

Calcium waves driven by "sensitization" wave-fronts

OBJECTIVE: Cellular Ca(2+) waves are understood as reaction-diffusion systems sustained by Ca(2+)-induced Ca(2+) release (CICR) from Ca(2+) stores. Given the recently discovered sensitization of Ca(2+) release channels (ryanodine receptors; RyRs) of the sarcoplasmic reticulum (SR) by luminal SR Ca(2+), waves could also be driven by RyR sensitization, mediated by SR overloading via Ca(2+) pump (SERCA), acting in tandem with CICR. METHODS: Confocal imaging of the Ca(2+) indicator fluo-3 was combined with UV-flash photolysis of caged compounds and the whole-cell configuration of the patch clamp technique to carry out these experiments in isolated guinea pig ventricular cardiomyocytes. RESULTS: Upon sudden slowing of the SERCA in cardiomyocytes with a photoreleased inhibitor, waves indeed decelerated immediately. No secondary changes of Ca(2+) signaling or SR Ca(2+) content due to SERCA inhibition were observed in the short time-frame of these experiments. CONCLUSIONS: Our findings are consistent with Ca(2+) loading resulting in a zone of RyR 'sensitization' traveling within the SR, but inconsistent with CICR as the predominant mechanism driving the Ca(2+) waves. This alternative mode of RyR activation is essential to fully conceptualize cardiac arrhythmias triggered by spontaneous Ca(2+) release.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Design of a Semi-analytical Method to Describe Plasma-Wall Interactions

Electric probes are objects immersed in the plasma with sharp boundaries which collect of emit charged particles. Consequently, the nearby plasma evolves under abrupt imposed and/or naturally emerging conditions. There could be localized currents, different time scales for plasma species evolution, charge separation and absorbing-emitting walls. The traditional numerical schemes based on differences often transform these disparate boundary conditions into computational singularities. This is the case of models using advection-diffusion differential equations with source-sink terms (also called Fokker-Planck equations). These equations are used in both, fluid and kinetic descriptions, to obtain the distribution functions or the density for each plasma species close to the boundaries. We present a resolution method grounded on an integral advancing scheme by using approximate Green's functions, also called short-time propagators. All the integrals, as a path integration process, are numerically calculated, what states a robust grid-free computational integral method, which is unconditionally stable for any time step. Hence, the sharp boundary conditions, as the current emission from a wall, can be treated during the short-time regime providing solutions that works as if they were known for each time step analytically. The form of the propagator (typically a multivariate Gaussian) is not unique and it can be adjusted during the advancing scheme to preserve the conserved quantities of the problem. The effects of the electric or magnetic fields can be incorporated into the iterative algorithm. The method allows smooth transitions of the evolving solutions even when abrupt discontinuities are present. In this work it is proposed a procedure to incorporate, for the very first time, the boundary conditions in the numerical integral scheme. This numerical scheme is applied to model the plasma bulk interaction with a charge-emitting electrode, dealing with fluid diffusion equations combined with Poisson equation self-consistently. It has been checked the stability of this computational method under any number of iterations, even for advancing in time electrons and ions having different time scales. This work establishes the basis to deal in future work with problems related to plasma thrusters or emissive probes in electromagnetic fields.

Twist grain boundaries in three-dimensional lamellar Turing structures

Steady spatial self-organization of three-dimensional chemical reaction-diffusion systems is discussed with the emphasis put on the possible defects that may alter the Turing patterns. It is shown that one of the stable defects of a three-dimensional lamellar Turing structure is a twist grain boundary embedding a Scherk minimal surface.