A new approach of fitting biomass dynamics models to data

A non-traditional approach of fitting dynamic resource biomass models to data is developed in this paper. A variational adjoint technique is used for dynamic parameter estimation. In the variational formulation, a cost function measuring the distance between the model solution and the observations is minimized. The data assimilation method provides a novel and computationally efficient procedure for combining all available information, i.e., the data and the model in the analysis of a resource system. This technique will be used to analyze data for the North-east Arctic cod stock. Two alternative population growth models: the logistic and the Gompertz model are used for estimating parameters of simple bioeconomic models by the method of constrained least squares. Estimates of the parameters of the models dynamics are reasonable and can be accepted. The main inference from the work is that the average fishing mortality is found to be significantly above the maximum sustainable yield value.     

Parameter identification of the calvin photosynthesis cycle

Summary This article is concerned with the determination of kinetic parameters of the Calvin photosynthesis cycle which is described by seventeen nonlinear ordinary differential equations. It is shown that the task requires dynamic data for several sets of initial conditions. The numerical technique is based upon an algorithm for non-linear optimization and Gear's numerical integration scheme for stiff systems of differential equations. The sensitivity of the parameters to noise in the data is tested with a method adapted from Rosenbrook and Storey. A preliminary set of parameters has been obtained from a preliminary set of experimental data. The numerical methods are then tested with synthetic data derived from these parameters. The mathematical model and the results obtained in the simulation are used as an aid in designing new experiments.     

Program Code Generator for Cardiac Electrophysiology Simulation with Automatic PDE Boundary Condition Handling

Clinical and experimental studies involving human hearts can have certain limitations. Methods such as computer simulations can be an important alternative or supplemental tool. Physiological simulation at the tissue or organ level typically involves the handling of partial differential equations (PDEs). Boundary conditions and distributed parameters, such as those used in pharmacokinetics simulation, add to the complexity of the PDE solution. These factors can tailor PDE solutions and their corresponding program code to specific problems. Boundary condition and parameter changes in the customized code are usually prone to errors and time-consuming. We propose a general approach for handling PDEs and boundary conditions in computational models using a replacement scheme for discretization. This study is an extension of a program generator that we introduced in a previous publication. The program generator can generate code for multi-cell simulations of cardiac electrophysiology. Improvements to the system allow it to handle simultaneous equations in the biological function model as well as implicit PDE numerical schemes. The replacement scheme involves substituting all partial differential terms with numerical solution equations. Once the model and boundary equations are discretized with the numerical solution scheme, instances of the equations are generated to undergo dependency analysis. The result of the dependency analysis is then used to generate the program code. The resulting program code are in Java or C programming language. To validate the automatic handling of boundary conditions in the program code generator, we generated simulation code using the FHN, Luo-Rudy 1, and Hund-Rudy cell models and run cell-to-cell coupling and action potential propagation simulations. One of the simulations is based on a published experiment and simulation results are compared with the experimental data. We conclude that the proposed program code generator can be used to generate code for physiological simulations and provides a tool for studying cardiac electrophysiology.     

A versatile digital computer program for non-linear regression analysis

1.

1. A general purpose digital computer program is described for application to biological experiments that require a non-linear regression analysis. The mathematical function, or model, to be fitted to a given set of experimental data is written as a section within the program. Given initial estimates for the parameters of the function, the program uses an iterative procedure to adjust the parameters until the sum of squares of residuals has converged to a minimum.     

Model parameterization to represent processes at unresolved scales and changing properties of evolving systems

Modeling has become an indispensable tool for scientific research. However, models generate great uncertainty when they are used to predict or forecast ecosystem responses to global change. This uncertainty is partly due to parameterization, which is an essential procedure for model specification via defining parameter values for a model. The classic doctrine of parameterization is that a parameter is constant. However, it is commonly known from modeling practice that a model that is well calibrated for its parameters at one site may not simulate well at another site unless its parameters are tuned again. This common practice implies that parameter values have to vary with sites. Indeed, parameter values that are estimated using a statistically rigorous approach, that is, data assimilation, vary with time, space, and treatments in global change experiments. This paper illustrates that varying parameters is to account for both processes at unresolved scales and changing properties of evolving systems. A model, no matter how complex it is, could not represent all the processes of one system at resolved scales. Interactions of processes at unresolved scales with those at resolved scales should be reflected in model parameters. Meanwhile, it is pervasively observed that properties of ecosystems change over time, space, and environmental conditions. Parameters, which represent properties of a system under study, should change as well. Tuning has been practiced for many decades to change parameter values. Yet this activity, unfortunately, did not contribute to our knowledge on model parameterization at all. Data assimilation makes it possible to rigorously estimate parameter values and, consequently, offers an approach to understand which, how, how much, and why parameters vary. To fully understand those issues, extensive research is required. Nonetheless, it is clear that changes in parameter values lead to different model predictions even if the model structure is the same.     

Computer simulation of the P31 NMR equations governing the creatine kinase reaction

A method for solving the Bloch equations that govern magnetization transfer NMR experiments is presented. It requires the numerical evaluation of a matrix exponential and lends itself to computer simulation. It is a simple, versatile method for evaluating P31-NMR magnetization transfer experiments designed to measure biochemical exchange rates. We apply the method to the saturation and inversion transfer experiments and find that the "initial slope" method of determining flux is subject to at least two interpretations. This verifies the accepted concept that biochemical models representing compartmentation and competing reaction hypotheses cannot be reliably distinguished by simply selecting values for NMR and biochemical parameters that give a "best fit" to experimental data. However, our results do indicate that a controlled manipulation of biochemical exchange rate may distinguish between these two models.     

The Routine Fitting of Kinetic Data to Models: A Mathematical Formalism for Digital Computers

A mathematical formalism is presented for use with digital computers to permit the routine fitting of data to physical and mathematical models. Given a set of data, the mathematical equations describing a model, initial conditions for an experiment, and initial estimates for the values of model parameters, the computer program automatically proceeds to obtain a least squares fit of the data by an iterative adjustment of the values of the parameters. When the experimental measures are linear combinations of functions, the linear coefficients for a least squares fit may also be calculated. The values of both the parameters of the model and the coefficients for the sum of functions may be unknown independent variables, unknown dependent variables, or known constants. In the case of dependence, only linear dependencies are provided for in routine use. The computer program includes a number of subroutines, each one of which performs a special task. This permits flexibility in choosing various types of solutions and procedures. One subroutine, for example, handles linear differential equations, another, special non-linear functions, etc. The use of analytic or numerical solutions of equations is possible.     

The effect of time-varying mortality and carbon assimilation on models of carbon allocation in annual plants

Summary Models of optimal carbon allocation schedules have influenced the way plant ecologists think about life history evolution, particularly for annual plants. The present study asks (1) how, within the framework of these models, are their predictions affected by within-season variation in mortality and carbon assimilation rates?; and (2) what are the consequences of these prediction changes for empirical tests of the models? A companion paper examines the basic assumptions of the models themselves. I conducted a series of numerical experiments with a simple carbon allocation model. Results suggest that both qualitative and quantitative predictions can sometimes be sensitive to parameter values for net assimilation rate and mortality: for some parameter values, both the time and size at onset of reproduction, as well as the number of reproductive intervals, vary considerably as a result of small variations in these parameters. For other parameter values, small variations in the parameters result in only small changes in predicted phenotype, but these have very large fitness consequences. Satisfactory empirical tests are thus likely to require much accuracy in parameter estimates. The effort required for parameter estimation imposes a practical constraint on empirical tests, making large multipopulation comparisons impractical. It may be most practical to compare the predicted and observed fitness consequences of variation in the timing of onset of reproduction.     

Monitoring Growth of Microorganisms using the Methods of Mass-Energy Balance and Utilization of Computer on Line with Fermenter in Real Time Processing (in Russian)

The effective means of microbial culture monitoring is the measurement of low-inertial parameters (respiration rate, rates of supply of alkali for pH maintenance and the limiting substrate) and utilization of computer on line with fermenter for recalculation of these rates into the instant values of mass and energy cell yields, specific rates of cell growth and substrate and oxygen consumption, using the method of mass-energy balance. In this paper, the equations of mass-energy balance are presented both in general form and in the form of numerical algorythms for computer programming. The installation for automation of microbial cultivation experiment is described. Experimental data are presented which indicate the effectiveness of the method of indirect measurement of cell biomass yield and specific rates of physiological processes.     

A model of the predator-prey relationship

A general nonlinear model of the predator-prey relationship is proposed. The model is justified on biological grounds and its parameters are constrained to reasonable biological values. Analytical and graphical methods are used to classify and to analyze the equilibrium points of the model. Finally, numerical integration on a digital computer illustrates the nature of the model solutions.     

Parameter estimation using Simulated Annealing for S-system models of biochemical networks

MOTIVATION: High-throughput technologies now allow the acquisition of biological data, such as comprehensive biochemical time-courses at unprecedented rates. These temporal profiles carry topological and kinetic information regarding the biochemical network from which they were drawn. Retrieving this information will require systematic application of both experimental and computational methods. RESULTS: S-systems are non-linear mathematical approximative models based on the power-law formalism. They provide a general framework for the simulation of integrated biological systems exhibiting complex dynamics, such as genetic circuits, signal transduction and metabolic networks. We describe how the heuristic optimization technique simulated annealing (SA) can be effectively used for estimating the parameters of S-systems from time-course biochemical data. We demonstrate our methods using three artificial networks designed to simulate different network topologies and behavior. We then end with an application to a real biochemical network by creating a working model for the cadBA system in Escherichia coli. AVAILABILITY: The source code written in C++ is available at http://www.engg.upd.edu.ph/~naval/bioinformcode.html. All the necessary programs including the required compiler are described in a document archived with the source code. SUPPLEMENTARY INFORMATION: Supplementary material is available at Bioinformatics online.     

A model of cell killing by low-dose-rate radiation including repair of sublethal damage, G2 block, and cell division

A computer model that simulates the killing of exponentially growing cells by low-dose-rate radiation is described. The model incorporates cell killing by single-hit damage and double-hit (sublethal) damage, as well as repair of sublethal damage, delay of cell cycle progression, blockage and increased sensitivity of cells in the G2 phase of the cell cycle, and cell division. Seven cellular parameters determine the rate of cell killing. Initial estimates of most of these parameters can be made from independent experiments. Parameters were obtained that gave the best fit to the data for four cell lines, using constant or variable dose rates, and using as end points either the fraction of single cells forming colonies or the total number of clonogenic cells in a mass culture. Some of the parameters were determined to be insignificant or similar for the four cell lines. The main differences between the cell lines in patterns of cell killing over a range of dose rates appeared to be determined by differences in the values of four of the parameters.     

Estimating kinetic constants from single channel data.

The process underlying the opening and closing of ionic channels in biological or artificial lipid membranes can be modeled kinetically as a time-homogeneous Markov chain. The elements of the chain are kinetic states that can be either open or closed. A maximum likelihood procedure is described for estimating the transition rates between these states from single channel data. The method has been implemented for linear kinetic schemes of fewer than six states, and is suitable for nonstationary data in which one or more independent channels are functioning simultaneously. It also provides standard errors for all estimates of rate constants and permits testing of smoothly parameterized subhypotheses of a general model. We have illustrated our approach by analysis of single channel data simulated on a computer and have described a procedure for analysis of experimental data.     

Robust dynamic experiments for the precise estimation of respiration and fermentation parameters of fruit and vegetables

Dynamic models based on non-linear differential equations are increasingly being used in many biological applications. Highly informative dynamic experiments are valuable for the identification of these dynamic models. The storage of fresh fruit and vegetables is one such application where dynamic experimentation is gaining momentum. In this paper, we construct optimal O₂ and CO₂ gas input profiles to estimate the respiration and fermentation kinetics of pear fruit. The optimal input profiles, however, depend on the true values of the respiration and fermentation parameters. Locally optimal design of input profiles, which uses a single initial guess for the parameters, is the traditional method to deal with this issue. This method, however, is very sensitive to the initial values selected for the model parameters. Therefore, we present a robust experimental design approach that can handle uncertainty on the model parameters.     

Analysis of single gene multitrait effects in livestock by the use of Gibbs sampling

The paper presents a method of multivariate data analysis described by a model which involves fixed effects, additive polygenic individual effects and the effects of a major gene. To find the estimates of model parameters, the maximization of likelihood function method is applied. The maximum of likelihood function is computed by the use of the Gibbs sampling approach. In this approach, following the conditional posterior distributions, values of all unknown parameters are generated. On the basis of the obtained samples the marginal posterior densities as well as the estimates of fixed effects, gene frequency, genotypic values, major gene, polygenic and error (co)variances are calculated. A numerical example, supplemented to theoretical considerations, deals with data simulated according to the considered model.     

Data-based identifiability analysis of non-linear dynamical models

MOTIVATION: Mathematical modelling of biological systems is becoming a standard approach to investigate complex dynamic, non-linear interaction mechanisms in cellular processes. However, models may comprise non-identifiable parameters which cannot be unambiguously determined. Non-identifiability manifests itself in functionally related parameters, which are difficult to detect. RESULTS: We present the method of mean optimal transformations, a non-parametric bootstrap-based algorithm for identifiability testing, capable of identifying linear and non-linear relations of arbitrarily many parameters, regardless of model size or complexity. This is performed with use of optimal transformations, estimated using the alternating conditional expectation algorithm (ACE). An initial guess or prior knowledge concerning the underlying relation of the parameters is not required. Independent, and hence identifiable parameters are determined as well. The quality of data at disposal is included in our approach, i.e. the non-linear model is fitted to data and estimated parameter values are investigated with respect to functional relations. We exemplify our approach on a realistic dynamical model and demonstrate that the variability of estimated parameter values decreases from 81 to 1% after detection and fixation of structural non-identifiabilities.     

A new framework for the estimation of control parameters in metabolic pathways using lin-log kinetics.

The control properties of biochemical pathways can be described by control coefficients and elasticities, as defined in the framework of metabolic control analysis. The determination of these parameters using the traditional metabolic control analysis relationships is, however, limited by experimental difficulties (e.g. realizing and measuring small changes in biological systems) and lack of appropriate mathematical procedures (e.g. when the more practical large changes are made). In this paper, the recently developed lin-log approach is proposed to avoid the above-mentioned problems and is applied to estimate control parameters from measurements obtained in steady state experiments. The lin-log approach employs approximative linear-logarithmic kinetics parameterized by elasticities and provides analytical solutions for fluxes and metabolite concentrations when large changes are made. Published flux and metabolite concentration data are used, obtained from a reconstructed section of glycolysis converting 3-phosphoglycerate to pyruvate [Giersch, C. (1995) Eur. J. Biochem. 227, 194-201]. With the lin-log approach, all data from different experiments can be combined to give realistic elasticity and flux control coefficient estimates by linear regression. Despite the large changes, a good agreement of fluxes and metabolite concentrations is obtained between the measured and calculated values according to the lin-log model. Furthermore, it is shown that the lin-log approach allows a rigorous statistical evaluation to identify the optimal reference state and the optimal model structure assumption. In conclusion, the lin-log approach addresses practical problems encountered in the traditional metabolic control analysis-based methods by introducing suitable nonlinear kinetics, thus providing a novel framework with improved procedures for the estimation of elasticities and control parameters from large perturbation experiments.     

A Biomechanical Triphasic Approach to the Transport of Nondilute Solutions in Articular Cartilage

Biomechanical models for biological tissues such as articular cartilage generally contain an ideal, dilute solution assumption. In this article, a biomechanical triphasic model of cartilage is described that includes nondilute treatment of concentrated solutions such as those applied in vitrification of biological tissues. The chemical potential equations of the triphasic model are modified and the transport equations are adjusted for the volume fraction and frictional coefficients of the solutes that are not negligible in such solutions. Four transport parameters, i.e., water permeability, solute permeability, diffusion coefficient of solute in solvent within the cartilage, and the cartilage stiffness modulus, are defined as four degrees of freedom for the model. Water and solute transport in cartilage were simulated using the model and predictions of average concentration increase and cartilage weight were fit to experimental data to obtain the values of the four transport parameters. As far as we know, this is the first study to formulate the solvent and solute transport equations of nondilute solutions in the cartilage matrix. It is shown that the values obtained for the transport parameters are within the ranges reported in the available literature, which confirms the proposed model approach.     

A study of absorption of energy of the extremely high frequency electromagnetic radiation in the rat skin by various dosimetric methods and approaches

Using experimental and theoretical methods of dosimetry, the energy absorption of extremely high-frequency electromagnetic radiation (EHF EMR) in the skin of laboratory rats was analyzed. Specific absorption rate (SAR) in the skin was determined on the basis of both microthermometric measurements of initial rates of temperature rise in rat skin induced by the exposure and microcalorimetric measurements of specific heat of the skin. Theoretical calculations of SAR in the skin were performed with consideration for dielectric parameters of rat skin obtained from the measurements of the standing wave ratio upon reflection of electromagnetic waves from the skin surface and for the effective area of stationary overheating measured by infrared thermography. A numerical method was developed to determine electromagnetic wave energy reflected, absorbed, and transmitted in the model of flat layers. The algorithm of the method was realized in a computer program and used to calculate SAR in the skin on the basis of the complex dielectric constant of rat skin. The SAR values obtained from experimental measurements, theoretical calculations and numerical analysis are in good mutual correspondence and make about 220-280 W/kg at a frequency of 42.25 GHz and a power of 20 mW at the radiator output. The results obtained can be used for dosimetric supply of biomedical experiments on studying the physicochemical mechanisms of the biological effects of EHF EMR.     

Flux module decomposition for parameter estimation in a multiple-feedback loop model of biochemical networks

Computer simulation is an important technique to capture the dynamics of biochemical networks. Since few quantitative values are measured in vivo, the values for unmeasured parameters should be estimated so that the simulation agrees with the experimental data. Considering the sparsity and error rates of experimentally measured data, the first thing is not to find a numerically exact and global solution but to explore a variety of the plausible parameter solutions. To find many plausible parameter solutions without any biases, we developed the two-phase search (TPS) method. However, calculation complexity makes it hard for TPS to optimize a large-scale dynamic model. In this study divide-and-conquer methods are used to solve this problem. The flux module decomposition (FMD) is first proposed that separates a complex, large-scale dynamic model into multiple flux modules without deteriorating its basic control architectures. FMD is combined with TPS, named FMD-TPS, to find many plausible parameter solutions for a dynamic model. To demonstrate the feasibility of FMD-TPS, it is applied to the E. coli ammonia assimilation system that consists of multiple-feedback loops. The variability of the solutions is verified by measuring the space distribution of the parameter solution vectors and by defining the binary vectors checking the consistency with biological behaviors. Compared with non-decomposition methods, FMD-TPS efficiently explored a variety of plausible parameter solutions that reproduce the dynamic behaviors in vivo.