Some novel techniques of parameter estimation for the dynamical models in biological systems

Inverse problems based on using experimental data to estimate unknown parameters of a system often arise in biological and chaotic systems. In this paper, we consider parameter estimation in systems biology involving linear and non-linear complex dynamical models, including the Michaelis–Menten enzyme kinetic system, a dynamical model of competence induction in Bacillus subtilis bacteria and a model of feedback bypass in B. subtilis bacteria. We propose some novel techniques for inverse problems. Firstly, we establish an approximation of a non-linear differential algebraic equation that corresponds to the given biological systems. Secondly, we use the Picard contraction mapping, collage methods and numerical integration techniques to convert the parameter estimation into a minimization problem of the parameters. We propose two optimization techniques: a grid approximation method and a modified hybrid Nelder–Mead simplex search and particle swarm optimization (MH-NMSS-PSO) for non-linear parameter estimation. The two techniques are used for parameter estimation in a model of competence induction in B. subtilis bacteria with noisy data. The MH-NMSS-PSO scheme is applied to a dynamical model of competence induction in B. subtilis bacteria based on experimental data and the model for feedback bypass. Numerical results demonstrate the effectiveness of our approach.

Exploring the parameter space of a rabbit ventricular action potential model to investigate the effect of variation on action potential and calcium transients

Computational models for cardiomyocyte action potentials (AP) often make use of a large parameter set. This parameter set can contain some elements that are fitted to experimental data independently of any other element, some elements that are derived concurrently with other elements to match experimental data, and some elements that are derived purely from phenomenological fitting to produce the desired AP output. Furthermore, models can make use of several different data sets, not always derived for the same conditions or even the same species. It is consequently uncertain whether the parameter set for a given model is physiologically accurate. Furthermore, it is only recently that the possibility of degeneracy in parameter values in producing a given simulation output has started to be addressed. In this study, we examine the effects of varying two parameters (the L-type calcium current (I(CaL)) and the delayed rectifier potassium current (I(Ks))) in a computational model of a rabbit ventricular cardiomyocyte AP on both the membrane potential (V(m)) and calcium (Ca(2+)) transient. It will subsequently be determined if there is degeneracy in this model to these parameter values, which will have important implications on the stability of these models to cell-to-cell parameter variation, and also whether the current methodology for generating parameter values is flawed. The accuracy of AP duration (APD) as an indicator of AP shape will also be assessed.

Parameter estimation for a phenomenological model of the cardiac action potential

The action potential (ap) of a cardiac cell is made up of a complex balance of ionic currents which flow across the cell membrane in response to electrical excitation of the cell. Biophysically detailed mathematical models of the ap have grown larger in terms of the variables and parameters required to model new findings in subcellular ionic mechanisms. The fitting of parameters to such models has seen a large degree of parameter and module re-use from earlier models. An alternative method for modelling electrically exciteable cardiac tissue is a phenomenological model, which reconstructs tissue level ap wave behaviour without subcellular details. A new parameter estimation technique to fit the morphology of the ap in a four variable phenomenological model is presented. An approximation of a nonlinear ordinary differential equation model is established that corresponds to the given phenomenological model of the cardiac ap. The parameter estimation problem is converted into a minimisation problem for the unknown parameters. A modified hybrid Nelder–Mead simplex search and particle swarm optimization is then used to solve the minimisation problem for the unknown parameters. The successful fitting of data generated from a well known biophysically detailed model is demonstrated. A successful fit to an experimental ap recording that contains both noise and experimental artefacts is also produced. The parameter estimation method’s ability to fit a complex morphology to a model with substantially more parameters than previously used is established.

Examining the black box between stimulating regions and individual entrepreneurship

This paper seeks to investigate the link between the objective regional opportunity structure (captured by regional data) and individuals’ engagement in different stages in the venture creation process (intention to start a business and engagement in nascent entrepreneurship). We further investigate pathways through which a favourable regional environment could affect entrepreneurial intentions and the propensity to be a nascent entrepreneur. We combine individual level GEM-data for Western Germany with regional level data from the statistical office and use multi-level analysis to test our hypotheses. We find support for our contention that a favourable regional opportunity structure affects entrepreneurial intentions and engagement. As pathways between the region and individual behaviour serve the individual perception of founding opportunities and the individual social capital.

Age-related loss of brain volume and T2 relaxation time in youth with type 1 diabetes

OBJECTIVE: Childhood-onset type 1 diabetes is associated with neurocognitive deficits, but there is limited evidence to date regarding associated neuroanatomical brain changes and their relationship to illness variables such as age at disease onset. This report examines age-related changes in volume and T2 relaxation time (a fundamental parameter of magnetic resonance imaging that reflects tissue health) across the whole brain. RESEARCH DESIGN AND METHODS: Type 1 diabetes, N = 79 (mean age 20.32 ± 4.24 years), and healthy control participants, N = 50 (mean age 20.53 ± 3.60 years). There were no substantial group differences on socioeconomic status, sex ratio, or intelligence quotient. RESULTS: Regression analyses revealed a negative correlation between age and brain changes, with decreasing gray matter volume and T2 relaxation time with age in multiple brain regions in the type 1 diabetes group. In comparison, the age-related decline in the control group was small. Examination of the interaction of group and age confirmed a group difference (type 1 diabetes vs. control) in the relationship between age and brain volume/T2 relaxation time. CONCLUSIONS: We demonstrated an interaction between age and group in predicting brain volumes and T2 relaxation time such that there was a decline in these outcomes in type 1 diabetic participants that was much less evident in control subjects. Findings suggest the neurodevelopmental pathways of youth with type 1 diabetes have diverged from those of their healthy peers by late adolescence and early adulthood but the explanation for this phenomenon remains to be clarified.

Variational Bayes and the reduced dependence approximation for the autologistic model on an irregular grid with applications

Discrete Markov random field models provide a natural framework for representing images or spatial datasets. They model the spatial association present while providing a convenient Markovian dependency structure and strong edge-preservation properties. However, parameter estimation for discrete Markov random field models is difficult due to the complex form of the associated normalizing constant for the likelihood function. For large lattices, the reduced dependence approximation to the normalizing constant is based on the concept of performing computationally efficient and feasible forward recursions on smaller sublattices which are then suitably combined to estimate the constant for the whole lattice. We present an efficient computational extension of the forward recursion approach for the autologistic model to lattices that have an irregularly shaped boundary and which may contain regions with no data; these lattices are typical in applications. Consequently, we also extend the reduced dependence approximation to these scenarios enabling us to implement a practical and efficient non-simulation based approach for spatial data analysis within the variational Bayesian framework. The methodology is illustrated through application to simulated data and example images. The supplemental materials include our C++ source code for computing the approximate normalizing constant and simulation studies.

Effect of soil parameter uncertainty on seismic response of buried segmented pipeline

Pipelines are important lifeline facilities spread over a large area and they generally encounter a range of seismic hazards and different soil conditions. The seismic response of a buried segmented pipe depends on various parameters such as the type of buried pipe material and joints, end restraint conditions, soil characteristics, burial depths, and earthquake ground motion, etc. This study highlights the effect of the variation of geotechnical properties of the surrounding soil on seismic response of a buried pipeline. The variations of the properties of the surrounding soil along the pipe are described by sampling them from predefined probability distribution. The soil-pipe interaction model is developed in OpenSEES. Nonlinear earthquake time-history analysis is performed to study the effect of soil parameters variability on the response of pipeline. Based on the results, it is found that uncertainty in soil parameters may result in significant response variability of the pipeline.

New exact solutions for Hele-Shaw flow in doubly connected regions

Radial Hele-Shaw flows are treated analytically using conformal mapping techniques. The geometry of interest has a doubly-connected annular region of viscous fluid surrounding an inviscid bubble that is either expanding or contracting due to a pressure difference caused by injection or suction of the inviscid fluid. The zero-surface-tension problem is ill-posed for both bubble expansion and contraction, as both scenarios involve viscous fluid displacing inviscid fluid. Exact solutions are derived by tracking the location of singularities and critical points in the analytic continuation of the mapping function. We show that by treating the critical points, it is easy to observe finite-time blow-up, and the evolution equations may be written in exact form using complex residues. We present solutions that start with cusps on one interface and end with cusps on the other, as well as solutions that have the bubble contracting to a point. For the latter solutions, the bubble approaches an ellipse in shape at extinction.

Early local functional changes in the human diabetic retina : a global flash multifocal electroretinogram study

Purpose: To investigate early functional changes of local retinal defects in type II diabetic patients using the global flash multifocal electroretinogram (MOFO mfERG). Methods: Thirty-eight diabetic patients and 14 age-matched controls were recruited. Nine of the diabetics were free from diabetic retinopathy (DR), while the remainder had mild to moderate non-proliferative diabetic retinopathy. The MOFO mfERG was performed at high (98%) and low (46%) contrast levels. MfERG responses were grouped into 35 regions for comparison with DR classification at those locations. Z-scores of the regional mfERG responses were compared across different types of DR defects. Results: The mfERG waveform consisted of the direct component (DC) and the induced component (IC). Local reduction in DC and IC amplitudes were found in diabetic patients with and without DR. With increasing severity of retinopathy, there was a further deterioration in amplitude of both components. Under MOFO mfERG paradigm, amplitude was a useful screening parameter. Conclusion: The MOFO mfERG can help in detecting early functional anomalies before the appearance of visible signs, and may assist in monitoring further functional deterioration in diabetic patients.

Territorial strategies and regions

The notion of territorial strategy emerged in the 1990s and has become more and more popular since. It refers to that combination of factors purposely assembled by governments, private and public companies, universities, and industrial associations to exploit a specific geographic competitive advantage in order to boost economic growth through the development of entrepreneurial activity and innovation. Three factors are generally considered to be the building blocks of a territorial strategy: natural resources, human capital, and industrial capabilities. Natural resources derive from environ­ mental conditions and represent raw materials or land available in a region. The presence of natural resources characterizes the typology of an industry (related to tourism, oil, wood, fish, and so forth) that exists or could exist in a certain area. Human capital refers to the stock of competences available in a certain region resulting from education and work experience. Industrial capabilities relate to complex constructs of specialized expertise, the confidence to apply knowledge and skills in various contexts and under changing conditions, and an ability repeatedly to improve methods and processes in a specific industry.

From slum to suburbia : examining a spectrum of planned and unplanned urban morphology in humid sub-tropical regions of the world

The current rapid urban growth throughout the world manifests in various ways and historically cities have grown, similarly, alternately or simultaneously between planned extensions and organic informal settlements (Mumford, 1989). Within cities different urban morphological regions can reveal different contexts of economic growth and/or periods of dramatic social/technological change (Whitehand, 2001, 105). Morpho-typological study of alternate contexts can present alternative models and contribute to the present discourse which questions traditional paradigms of urban planning and design (Todes et al, 2010). In this study a series of cities are examined as a preliminary exploration into the urban morphology of cities in ‘humid subtropical’ climates. From an initial set of twenty, six cities were selected: Sao Paulo, Brazil; Jacksonville, USA; Maputo, Mozambique; Kanpur, India; Hong Kong, China; and Brisbane, Australia. The urban form was analysed from satellite imagery at a constant scale. Urban morphological regions (types) were identified as those demonstrating particular consistant characteristics of form (density, typology and pattern) different to their surroundings when examined at a constant scale. This analysis was correlated against existing data and literature discussing the proliferation of two types of urban development, ‘informal settlement’ (defined here as self-organised communities identifiable but not always synonymous with ‘slums’) and ‘suburbia’ (defined here as master planned communities of generally detached houses prevalent in western society) - the extreme ends of a hypothetical spectrum from ‘planned’ to ‘spontaneous’ urban development. Preliminary results show some cities contain a wide variety of urban form ranging from the highly organic ‘self-organised’ type to the highly planned ‘master planned community’ (in the case of Sao Paulo) while others tend to fall at one end of the planning spectrum or the other (more planned in the cases of Brisbane and Jacksonville; and both highly planned and highly organic in the case of Maputo). Further research will examine the social, economical and political drivers and controls which lead to this diversity or homogeneity of urban form and speculates on the role of self-organisation as a process for the adaptation of urban form.

Inherent errors in pollutant build-up estimation in considering urban land use as a lumped parameter

Stormwater quality modelling results is subject to uncertainty. The variability of input parameters is an important source of overall model error. An in-depth understanding of the variability associated with input parameters can provide knowledge on the uncertainty associated with these parameters and consequently assist in uncertainty analysis of stormwater quality models and the decision making based on modelling outcomes. This paper discusses the outcomes of a research study undertaken to analyse the variability related to pollutant build-up parameters in stormwater quality modelling. The study was based on the analysis of pollutant build-up samples collected from 12 road surfaces in residential, commercial and industrial land uses. It was found that build-up characteristics vary appreciably even within the same land use. Therefore, using land use as a lumped parameter would contribute significant uncertainties in stormwater quality modelling. Additionally, it was also found that the variability in pollutant build-up can also be significant depending on the pollutant type. This underlines the importance of taking into account specific land use characteristics and targeted pollutant species when undertaking uncertainty analysis of stormwater quality models or in interpreting the modelling outcomes.

An assessment tool for selection of appropriate daylighting solutions for buildings in tropical and subtropical regions : Validation using radiance simulation

A new decision-making tool that will assist designers in the selection of appropriate daylighting solutions for buildings in tropical locations has been previously proposed by the authors. Through an evaluation matrix that prioritizes the parameters that best respond to the needs of tropical climates (e.g. reducing solar gain and protection from glare) the tool determines the most appropriate devices for specific climate and building inputs. The tool is effective in demonstrating the broad benefits and limitations of the different daylight strategies for buildings in the tropics. However for thorough analysis and calibration of the tool, validation is necessary. This paper presents a first step in the validation process. RADIANCE simulations were conducted to compare simulation performance with the performance predicted by the tool. To this end, an office building case study in subtropical Brisbane, Australia, and five different daylighting devices including openings, light guiding systems and light transport systems were simulated. Illuminance, light uniformity, daylight penetration and glare analysis were assessed for each device. The results indicate the tool can appropriately rank and recommend daylighting strategies based on specific building inputs for tropical and subtropical regions, making it a useful resource for designers.