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.

Sensitivity of the NMR density matrix to pulse sequence parameters : a simplified analytic approach

We present a formalism for the analysis of sensitivity of nuclear magnetic resonance pulse sequences to variations of pulse sequence parameters, such as radiofrequency pulses, gradient pulses or evolution delays. The formalism enables the calculation of compact, analytic expressions for the derivatives of the density matrix and the observed signal with respect to the parameters varied. The analysis is based on two constructs computed in the course of modified density-matrix simulations: the error interrogation operators and error commutators. The approach presented is consequently named the Error Commutator Formalism (ECF). It is used to evaluate the sensitivity of the density matrix to parameter variation based on the simulations carried out for the ideal parameters, obviating the need for finite-difference calculations of signal errors. The ECF analysis therefore carries a computational cost comparable to a single density-matrix or product-operator simulation. Its application is illustrated using a number of examples from basic NMR spectroscopy. We show that the strength of the ECF is its ability to provide analytic insights into the propagation of errors through pulse sequences and the behaviour of signal errors under phase cycling. Furthermore, the approach is algorithmic and easily amenable to implementation in the form of a programming code. It is envisaged that it could be incorporated into standard NMR product-operator simulation packages.

Optimization of human papillomavirus type 16 (HPV-16) L1 expression in plants: Comparison of the suitability of different HPV-16 L1 gene variants and different cell-compartment localization

Virus-like particle-based vaccines for high-risk human papillomaviruses (HPVs) appear to have great promise; however, cell culture-derived vaccines will probably be very expensive. The optimization of expression of different codon-optimized versions of the HPV-16 L1 capsid protein gene in plants has been explored by means of transient expression from a novel suite of Agrobacterium tumefaciens binary expression vectors, which allow targeting of recombinant protein to the cytoplasm, endoplasmic reticulum (ER) or chloroplasts. A gene resynthesized to reflect human codon usage expresses better than the native gene, which expresses better than a plant-optimized gene. Moreover, chloroplast localization allows significantly higher levels of accumulation of L1 protein than does cytoplasmic localization, whilst ER retention was least successful. High levels of L1 (>17% total soluble protein) could be produced via transient expression: the protein assembled into higher-order structures visible by electron microscopy, and a concentrated extract was highly immunogenic in mice after subcutaneous injection and elicited high-titre neutralizing antibodies. Transgenic tobacco plants expressing a human codon-optimized gene linked to a chloroplast-targeting signal expressed L1 at levels up to 11% of the total soluble protein. These are the highest levels of HPV L1 expression reported for plants: these results, and the excellent immunogenicity of the product, significantly improve the prospects of making a conventional HPV vaccine by this means. © 2007 SGM.

A multivariate approach to the identification of surrogate parameters for heavy metals in stormwater

Stormwater is a potential and readily available alternative source for potable water in urban areas. However, its direct use is severely constrained by the presence of toxic pollutants, such as heavy metals (HMs). The presence of HMs in stormwater is of concern because of their chronic toxicity and persistent nature. In addition to human health impacts, metals can contribute to adverse ecosystem health impact on receiving waters. Therefore, the ability to predict the levels of HMs in stormwater is crucial for monitoring stormwater quality and for the design of effective treatment systems. Unfortunately, the current laboratory methods for determining HM concentrations are resource intensive and time consuming. In this paper, applications of multivariate data analysis techniques are presented to identify potential surrogate parameters which can be used to determine HM concentrations in stormwater. Accordingly, partial least squares was applied to identify a suite of physicochemical parameters which can serve as indicators of HMs. Datasets having varied characteristics, such as land use and particle size distribution of solids, were analyzed to validate the efficacy of the influencing parameters. Iron, manganese, total organic carbon, and inorganic carbon were identified as the predominant parameters that correlate with the HM concentrations. The practical extension of the study outcomes to urban stormwater management is also discussed.

Resources for improving journalists' understanding of the economic and financial parameters of Australia's agricultural and commodity producers

The research addresses how an understanding of the fundamentals of economics will better inform general journalists who report on issues or events affecting rural and regional Australia. The research draws on practice-based experience of the author, formal economics studies, interviews with news editors from Australian television news organisations, and interviews from leading economists. A guidebook has also been written to help journalists apply economic theories to their reporting. The guidebook enables reporters to think strategically and consider the 'big picture' when they inform society about policies, commodity trade, the environment, or any issues involving rural and regional Australia.

The essential parameters of a resource-based carrying capacity assessment model : an Australian case study

Carrying capacity assessments model a population’s potential self-sufficiency. A crucial first step in the development of such modelling is to examine the basic resource-based parameters defining the population’s production and consumption habits. These parameters include basic human needs such as food, water, shelter and energy together with climatic, environmental and behavioural characteristics. Each of these parameters imparts land-usage requirements in different ways and varied degrees so their incorporation into carrying capacity modelling also differs. Given that the availability and values of production parameters may differ between locations, no two carrying capacity models are likely to be exactly alike. However, the essential parameters themselves can remain consistent so one example, the Carrying Capacity Dashboard, is offered as a case study to highlight one way in which these parameters are utilised. While examples exist of findings made from carrying capacity assessment modelling, to date, guidelines for replication of such studies in other regions and scales have largely been overlooked. This paper addresses such shortcomings by describing a process for the inclusion and calibration of the most important resource-based parameters in a way that could be repeated elsewhere.

Mechanical stimulation of the pro-angiogenic capacity of human fracture haematoma: Involvement of VEGF mechano-regulation

Compromised angiogenesis appears to be a major limitation in various suboptimal bone healing situations. Appropriate mechanical stimuli support blood vessel formation in vivo and improve healing outcomes. However, the mechanisms responsible for this association are unclear. To address this question, the paracrine angiogenic potential of early human fracture haematoma and its responsiveness to mechanical loading, as well as angiogenic growth factors involved, were investigated in vitro. Human haematomas were collected from healthy patients undergoing surgery within 72. h after bone fracture. The haematomas were embedded in a fibrin matrix, and cultured in a bioreactor resembling the in vivo conditions of the early phase of bone healing (20 compression, 1. Hz) over 3. days. Conditioned medium (CM) from the bioreactor was then analyzed. The matrices were also incubated in fresh medium for a further 24. h to evaluate the persistence of the effects. Growth factor (GF) concentrations were measured in the CM by ELISAs. In vitro tube formation assays were conducted on Matrigel with the HMEC-1 cell line, with or without inhibition of vascular endothelial growth factor receptor 2 (VEGFR2). Cell numbers were quantified using an MTS test. In vitro endothelial tube formation was enhanced by CM from haematomas, compared to fibrin controls. The angiogenesis regulators, vascular endothelial growth factor (VEGF) and transforming growth factor β1 (TGF-β1), were released into the haematoma CM, but not angiopoietins 1 or 2 (Ang1, 2), basic fibroblast growth factor (bFGF) or platelet-derived growth factor (PDGF). Mechanical stimulation of haematomas, but not fibrin controls, further increased the induction of tube formation by their CM. The mechanically stimulated haematoma matrices retained their elevated pro-angiogenic capacity for 24. h. The pro-angiogenic effect was cancelled by inhibition of VEGFR2 signalling. VEGF concentrations in CM tended to be elevated by mechanical stimulation; this was significant in haematomas from younger, but not from older patients. Other GFs were not mechanically regulated. In conclusion, the paracrine pro-angiogenic capacity of early human haematomas is enhanced by mechanical stimulation. This effect lasts even after removing the mechanical stimulus and appears to be VEGFR2-dependent.

Predicting the fatigue life of internal fracture fixation plates

Failures of fracture fixation plates, often related to fatigue fractures of the implants, have been reported (Banovetz et al, 1996). While metallurgical defects can usually be excluded, many of these fractures can be explained with the biomechanical situation. This study investigated the biomechanics of two clinical cases, both of which used a 14-hole locking compression plate. In the first case, a titanium plate was used in a rigid configuration with 12 screws resulting in plate breakage after 7 weeks (Sommer et al, 2003). In the second case, a stainless steel plate, which endured the entire healing process, was used in a flexible application with only 6 screws.

Numerical modelling and analysis of the blast performance of laminated glass panels and the influence of material parameters

This paper presents a comprehensive numerical procedure to treat the blast response of laminated glass (LG) panels and studies the influence of important material parameters. Post-crack behaviour of the LG panel and the contribution of the interlayer towards blast resistance are treated. Modelling techniques are validated by comparing with existing experimental results. Findings indicate that the tensile strength of glass considerably influences the blast response of LG panels while the interlayer material properties have a major impact on the response under higher blast loads. Initially, glass panes absorb most of the blast energy, but after the glass breaks, interlayer deforms further and absorbs most of the blast energy. LG panels should be designed to fail by tearing of the interlayer rather than failure at the supports to achieve a desired level of protection. From this aspect, material properties of glass, interlayer and sealant joints play important roles, but unfortunately they are not accounted for in the current design standards. The new information generated in this paper will enhance the capabilities of engineers to better design LG panels under blast loads and use better materials to improve the blast response of LG panels.

Improvement of the mode II interface fracture toughness of glass fiber reinforced plastics/aluminum laminates through vapor grown carbon fiber interleaves

The effects of acid treatment, vapor grown carbon fiber (VGCF) interlayer and the angle, i.e., 0° and 90°, between the rolling stripes of an aluminum (Al) plate and the fiber direction of glass fiber reinforced plastics (GFRP) on the mode II interlaminar mechanical properties of GFRP/Al laminates were investigated. The experimental results of an end notched flexure test demonstrate that the acid treatment and the proper addition of VGCF can effectively improve the critical load and mode II fracture toughness of GFRP/Al laminates. The specimens with acid treatment and 10 g m−2 VGCF addition possess the highest mode II fracture toughness, i.e., 269% and 385% increases in the 0° and 90° specimens, respectively compared to those corresponding pristine ones. Due to the induced anisotropy by the rolling stripes on the aluminum plate, the 90° specimens possess 15.3%–73.6% higher mode II fracture toughness compared to the 0° specimens. The improvement mechanisms were explored by the observation of crack propagation path and fracture surface with optical, laser scanning and scanning electron microscopies. Moreover, finite element analyses were carried out based on the cohesive zone model to verify the experimental fracture toughness and to predict the interface shear strength between the aluminum plates and GFRP laminates.