An investigation into the cytotoxic effects of microbubbles and their constituents on osteosarcoma and bone marrow stromal cells.

BACKGROUND: Ultrasound-responsive microbubbles offer a means of achieving minimally invasive, localised drug delivery in applications including regenerative medicine. To facilitate their use, however, it is important to determine any cytotoxic effects they or their constituents may have. The aim of this study was to test the hypothesis that phospholipid-shelled microbubbles are non-toxic to human bone-derived cells at biologically-relevant concentrations. METHODS: Microbubbles were fabricated using combinations of 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC), 1,2-dibehenoyl-sn-glycero-3-phosphocholine (DBPC), polyoxyethylene(40) stearate (PEG40S) and 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[methoxy(polyethylene-glycol)-2000] (DSPE-PEG2000). Microbubble size and concentration were measured as a function of time and temperature by optical microscopy. Effects on MG63 osteosarcoma and human bone marrow stromal cells (BMSCs) were measured for up to 72 h by assay for viability, metabolic activity and proliferation. RESULTS: DBPC:DSPE-PEG2000 microbubbles were significantly more stable than DSPC:PEG40S microbubbles under all conditions tested. Serum-containing medium had no detrimental effect on microbubble stability, but storage at 37 °C compared to at 4 °C reduced stability for both preparations, with almost complete dissolution of microbubbles at times ≥24 h. DSPC:PEG40S microbubbles had greater inhibitory effects on cell metabolism and growth than DBPC:DSPE-PEG2000 microbubbles, with PEG40S found to be the principle inhibitory component. These effects were only evident at high microbubble concentrations (≥20% (v/v)) or with prolonged culture (≥24 h). Increasing cell-microbubble contact by inversion culture in a custom-built device had no inhibitory effect on metabolism. CONCLUSIONS: These data indicate that, over a broad range of concentrations and incubation times, DBPC:DSPE-PEG2000 and DSPC:PEG40S microbubbles have little effect on osteoblastic cell viability and growth, and that PEG40S is the principle inhibitory component in the formulations investigated.

Predicting QRS and PR interval prolongations in humans using nonclinical data.

BACKGROUND AND PURPOSE: Risk of cardiac conduction slowing (QRS/PR prolongations) is assessed prior to clinical trials using in vitro and in vivo studies. Understanding the quantitative translation of these studies to the clinical situation enables improved risk assessment in the nonclinical phase. EXPERIMENTAL APPROACH: Four compounds that prolong QRS and/or PR (AZD1305, flecainide, quinidine and verapamil) were characterized using in vitro (sodium/calcium channels), in vivo (guinea pigs/dogs) and clinical data. Concentration-matched translational relationships were developed based on in vitro and in vivo modelling, and the in vitro to clinical translation of AZD1305 was quantified using an in vitro model. KEY RESULTS: Meaningful (10%) human QRS/PR effects correlated with low levels of in vitro Nav 1.5 block (3-7%) and Cav 1.2 binding (13-21%) for all compounds. The in vitro model developed using AZD1305 successfully predicted QRS/PR effects for the remaining drugs. Meaningful QRS/PR changes in humans correlated with small effects in guinea pigs and dogs (QRS 2.3-4.6% and PR 2.3-10%), suggesting that worst-case human effects can be predicted by assuming four times greater effects at the same concentration from dog/guinea pig data. CONCLUSION AND IMPLICATIONS: Small changes in vitro and in vivo consistently translated to meaningful PR/QRS changes in humans across compounds. Assuming broad applicability of these approaches to assess cardiovascular safety risk for non-arrhythmic drugs, this study provides a means of predicting human QRS/PR effects of new drugs from effects observed in nonclinical studies.

A mathematical model characterising Achilles tendon dynamics in flexion.

The purpose of this study is to acquire mechanistic knowledge of the gastrocnemius muscle-Achilles tendon complex behaviour during specific movements in humans through mathematical modelling. Analysis of this muscle-tendon complex was performed to see if already existing muscle-tendon models of other parts of the body could be applied to the leg muscles, especially the gastrocnemius muscle-Achilles tendon complex, and whether they could adequately characterise its behaviour. Five healthy volunteers were asked to take part in experiments where dorsiflexion and plantar flexion of the foot were studied. A model of the Achilles tendon-gastrocnemius muscle was developed, incorporating assumptions regarding the mechanical properties of the muscle fibres and the tendinous tissue in series. Ultrasound images of the volunteers, direct measurements and additional mathematical calculations were used to parameterise the model. Ground reaction forces, forces on specific joints and moments and angles for the ankle were obtained from a Vicon 3D motion capture system. Model validation was performed from the experimental data captured for each volunteer and from reconstruction of the movements of specific trajectories of the joints, muscles and tendons involved in those movements.

Solubility constraints affecting the migration of selenium through the cementitious backfill of a geological disposal facility.

This work presents the study of the solubility of selenium under cementitious conditions and its diffusion, as SeO3(2-), through monolithic cement samples. The solubility studies were carried out under alkaline conditions similar to those anticipated in the near-field of a cement-based repository for low- and intermediate-level radioactive waste. Experiments were conducted in NaOH solution, 95%-saturated Ca(OH)2, water equilibrated with a potential backfill material (Nirex reference vault backfill, NRVB) and in solutions containing cellulose degradation products, with and without reducing agents. The highest selenium concentrations were found in NaOH solution. In the calcium-containing solutions, analysis of the precipitates suggests that the solubility controlling phase is Ca2SeO3(OH)2·2H2O, which appears as euhedral rhombic crystals. The presence of cellulose degradation products caused an increase in selenium concentration, possibly due to competitive complexation, thereby, limiting the amount of calcium available for precipitation. Iron coupons had a minor effect on selenium solubility in contrast to Na2S2O4, suggesting that effective reduction of Se(IV) occurs only at Eh values below -300mV. Radial through-diffusion experiments on NRVB and in a fly ash cement showed no evidence of selenium breakthrough after one year. However, autoradiography of the exposed surfaces indicated that some migration had occurred and that selenium was more mobile in the higher porosity backfill than in the fly ash cement.

Effect of anthropogenic organic complexants on the solubility of Ni, Th, U(IV) and U(VI).

The influence of anthropogenic organic complexants (citrate, EDTA and DTPA from 0.005 to 0.1M) on the solubility of nickel(II), thorium(IV) and uranium (U(IV) and U(VI)) has been studied. Experiments were carried out in 95%-saturated Ca(OH)2 solutions, representing the high pH conditions anticipated in the near field of a cementitious intermediate level radioactive waste repository. Results showed that Ni(II) solubility increased by 2-4 orders of magnitude in the presence of EDTA and DTPA and from 3 to 4 orders of magnitude in the case of citrate. Citrate had the greatest effect on the solubility of Th(IV) and U(IV)/(VI). XRD and SEM analyses indicate that the precipitates are largely amorphous; only in the case of Ni(II), is there some evidence of incipient crystallinity, in the form of Ni(OH)2 (theophrastite). A study of the effect of calcium suggests that U(VI) and Ni(II) may form metal-citrate-OH complexes stabilised by Ca(2+). Thermodynamic modelling underestimates the concentrations in solution in the presence of the ligands for all the elements considered here. Further investigation of the behaviour of organic ligands under hyperalkaline conditions is important because of the use of the thermodynamic constants in preparing the safety case for the geological disposal of radioactive wastes.

Kinetic modelling of haemodialysis removal of myoglobin in rhabdomyolysis patients.

An extended two compartment model is proposed to describe the dynamics of myoglobin in rhabdomyolysis patients undergoing dialysis. Before using clinical data to estimate the model's unknown parameters, structural identifiability analysis was performed to determine the parameters uniqueness given certain clinical observations. A Taylor series expansion method was implemented which found that the model was structurally globally/uniquely identifiable for both on- and off-dialysis phases. The fitted model was then used in a predictive capacity showing that the use of Theralite high cut-off (HCO) or HCO 1100 dialyser gave a significant reduction in myoglobin renal exposure compared to standard haemodialysis (HD).

Describing the effectiveness of immunosuppression drugs and apheresis in the treatment of transplant patients.

When any foreign object is found in the human body antibodies are generated that mark it for removal by the immune system. In most cases these are natural and healthy responses; however, when considering organ transplants the immune response to the implanted organ must be kept to a minimum to avoid host rejection. To reduce the host's immune response to the implant, clinicians are able to manipulate the antibody dynamics through drug therapy, to minimise the antibody synthesis (immunosuppression), and by the removal of antibodies directly from the patients' blood, a process known as apheresis. In this paper models are presented that describe the in vivo kinetics of three immune complexes which are routinely measured pre- and post-operatively in implant patients, namely IgA, IgG and IgM. These models are then used to analyse the effective clearance rates of different apheresis methods (plasmapheresis, plasma absorption or plasma exchange) and to quantify the impact immune-suppression drugs have on the underlying antibody synthesis. It is hoped that the simplicity of the mathematical models, and associated implementation, will allow the translation of knowledge gained of the process dynamics to positively impact future patient diagnosis and treatment.

Dose Adjustment for Normal Eating (DAFNE) in routine clinical practice: who benefits?

AIMS: To explore the effectiveness of Dose Adjustment for Normal Eating in routine clinical practice in the UK. METHODS: Participants were 124 adults with Type 1 diabetes who had completed a Dose Adjustment for Normal Eating course. Data were collected before the course and again 1 year later on a variety of biological, psychological and social measures. RESULTS: There were a range of significant benefits consistent with Dose Adjustment for Normal Eating aims, including: better control among those with baseline HbA(1c) ≥ 81 mmol/mol (9.6%) (z = -2.8, P = 0.004); reduced number of participants reporting severe hypoglycaemia (χ² = 4.27, P = 0.039); total eradication of diabetic ketoacidosis (χ² = 4.17, P = 0.041) and lower diabetes-related distress (z = -4.5, P < 0.001). The most deprived of the clinic population were significantly under-represented (χ² = 17.8, P = 0.001) and the levels of clinical depression were unusually low. CONCLUSIONS: These results indicate that Dose Adjustment for Normal Eating delivered in routine clinical practice is associated with a range of benefits and that certain clinical and psychosocial characteristics are associated with better outcomes.

Indistinguishability and identifiability of kinetic models for the MurC reaction in peptidoglycan biosynthesis.

An important question in Systems Biology is the design of experiments that enable discrimination between two (or more) competing chemical pathway models or biological mechanisms. In this paper analysis is performed between two different models describing the kinetic mechanism of a three-substrate three-product reaction, namely the MurC reaction in the cytoplasmic phase of peptidoglycan biosynthesis. One model involves ordered substrate binding and ordered release of the three products; the competing model also assumes ordered substrate binding, but with fast release of the three products. The two versions are shown to be distinguishable; however, if standard quasi-steady-state assumptions are made distinguishability cannot be determined. Once model structure uniqueness is ensured the experimenter must determine if it is possible to successfully recover rate constant values given the experiment observations, a process known as structural identifiability. Structural identifiability analysis is carried out for both models to determine which of the unknown reaction parameters can be determined uniquely, or otherwise, from the ideal system outputs. This structural analysis forms an integrated step towards the modelling of the full pathway of the cytoplasmic phase of peptidoglycan biosynthesis.

The use of a formal sensitivity analysis on epidemic models with immune protection from maternally acquired antibodies.

This paper considers the outcome of a formal sensitivity analysis on a series of epidemic model structures developed to study the population level effects of maternal antibodies. The analysis is used to compare the potential influence of maternally acquired immunity on various age and time domain observations of infection and serology, with and without seasonality. The results of the analysis indicate that time series observations are largely insensitive to variations in the average duration of this protection, and that age related empirical data are likely to be most appropriate for estimating these characteristics.

Kinetic modelling of the role of the aldehyde dehydrogenase enzyme and the breast cancer resistance protein in drug resistance and transport.

A compartmental model for the in vitro uptake kinetics of the anti-cancer agent topotecan (TPT) has been extended from a previously published model. The extended model describes the drug activity and delivery of the pharmacologically active form to the DNA target as well as the catalysis of the aldehyde dehydrogenase (ALDH) enzyme and the elimination of drug from the cytoplasm via the efflux pump. Verification of the proposed model is achieved using scanning-laser microscopy data from live human breast cancer cells. Before estimating the unknown model parameters from the experimental in vitro data it is essential to determine parameter uniqueness (or otherwise) from this imposed output structure. This is formally performed as a structural identifiability analysis, which demonstrates that all of the unknown model parameters are uniquely determined by the output structure corresponding to the experiment.

A coupled drug kinetics-cell cycle model to analyse the response of human cells to intervention by topotecan.

A model describing the response of the growth of single human cells in the absence and presence of the anti-cancer agent topotecan (TPT) is presented. The model includes a novel coupling of both the kinetics of TPT and cell cycle responses to the agent. By linking the models in this way, rather than using separate (disjoint) approaches, it is possible to illustrate how the drug perturbs the cell cycle. The model is compared to experimental in vitro cell cycle response data (comprising single cell descriptors for molecular and behavioural events), showing good qualitative agreement for a range of TPT dose levels.

Competitive effect of iron(III) on metal complexation by humic substances: characterisation of ageing processes.

Aiming at an assessment of counteractive effects on colloid-borne migration of actinides in the event of release from an underground repository, competition by Fe(III) in respect of metal complexation by dissolved organic matter was investigated for the example of Eu(III) as an analogue of trivalent actinides. Complexation with different humic materials was examined in cation exchange experiments, using (59)Fe and (152)Eu as radioactive tracers for measurements in dilute systems as encountered in nature. Competitive effects proved to be significant when Fe is present at micromolar concentrations. Flocculation as a limiting process was attributed to charge compensation of humic colloids. Fe fractions bound to humic acids (HA) were higher than 90%, exceeding the capacity of binding sites at high Fe concentrations. It is thus concluded that the polynuclear structure of hydrolysed Fe(III) is maintained when bound to HA, which is also inferred from UV-Vis spectrometry. The competitive effect was found to be enhanced if Fe and HA were in contact before Eu was added. Depending on the time of Fe/HA pre-equilibration, Eu complexation decreased asymptotically over a time period of several weeks, the amount of bound Fe being unchanged. Time-dependent observations of UV-Vis spectra and pH values revealed that the ageing effect was due to a decline in Fe hydrolysis rather than structural changes within HA molecules. Fe polycations are slowly degraded in contact with humic colloids, and more binding sites are occupied as a consequence of dispersion. The extent of degradation as derived from pH shifts depended on the Fe/HA ratio.

Glucose-dependent changes in NAD(P)H-related fluorescence lifetime of adipocytes and fibroblasts in vitro: potential for non-invasive glucose sensing in diabetes mellitus.

The aim of this study was to test the hypothesis that glucose can be monitored non-invasively by measuring NAD(P)H-related fluorescence lifetime of cells in an in vitro cell culture model. Autofluorescence decay functions were measured in 3T3-L1 adipocytes by time-correlated single-photon counting (excitation 370nm, emission 420-480nm). Free NADH had a two-exponential decay but cell autofluorescence fitted best to a three-exponential decay. Addition of 30mM glucose caused a 29% increase in autofluorescence intensity, a significantly shortened mean lifetime (from 7.23 to 6.73ns), and an increase in the relative amplitude and fractional intensity of the short-lifetime component at the expense of the two longer-lifetime components. Similar effects were seen with rotenone, an agent that maximizes mitochondrial NADH. 3T3-L1 fibroblasts stained with the fluorescent mitochondrial marker, rhodamine 123 showed a 16% quenching of fluorescence intensity when exposed to 30mM glucose, and an increase in the relative amplitude and fractional intensity of the short lifetime at the expense of the longer lifetime component. We conclude that, though the effect size is relatively small, glucose can be measured non-invasively in cells by monitoring changes in the lifetimes of cell autofluorescence or of a dye marker of mitochondrial metabolism. Further investigation and development of fluorescence intensity and lifetime sensing is therefore indicated for possible non-invasive metabolic monitoring in human diabetes.

The structural identifiability and parameter estimation of a multispecies model for the transmission of mastitis in dairy cows with postmilking teat disinfection.

A mathematical model for the transmission of two interacting classes of mastitis causing bacterial pathogens in a herd of dairy cows is presented and applied to a specific data set. The data were derived from a field trial of a specific measure used in the control of these pathogens, where half the individuals were subjected to the control and in the others the treatment was discontinued. The resultant mathematical model (eight non-linear simultaneous ordinary differential equations) therefore incorporates heterogeneity in the host as well as the infectious agent and consequently the effects of control are intrinsic in the model structure. A structural identifiability analysis of the model is presented demonstrating that the scope of the novel method used allows application to high order non-linear systems. The results of a simultaneous estimation of six unknown system parameters are presented. Previous work has only estimated a subset of these either simultaneously or individually. Therefore not only are new estimates provided for the parameters relating to the transmission and control of the classes of pathogens under study, but also information about the relationships between them. We exploit the close link between mathematical modelling, structural identifiability analysis, and parameter estimation to obtain biological insights into the system modelled.

The structural identifiability and parameter estimation of a multispecies model for the transmission of mastitis in dairy cows.

A structural identifiability analysis is performed on a mathematical model for the coupled transmission of two classes of pathogen. The pathogens, classified as major and minor, are aetiological agents of mastitis in dairy cows that interact directly and via the immunological reaction in their hosts. Parameter estimates are available from experimental data for all but four of the parameters in the model. Data from a longitudinal study of infection are used to estimate these unknown parameters. A novel approach and application of structural identifiability analysis is combined in this paper with the estimation of cross-protection parameters using epidemiological data.

A control theoretic approach to containing the spread of rabies.

Many problems in medicine and biology involve some kind of spatial spread, and quite often the need to control it. A large proportion of medical and biological systems distinguish themselves from those found in engineering by the way the control acts. We illustrate this by considering the specific example of the spread of rabies among foxes. We first give a brief description of a model proposed by Murray et al. (Murray, J. D., Stanley, E. A. & Brown, D. L., Proc. R. Soc. Lond., B 229, 111-150 (1986)), which we extend to include the control mechanism. The problem is to prevent the spread of rabies by vaccinating foxes via the distribution of bait in a region around an observed outbreak. The extended model can be formulated as a nonlinear time-varying control system described by partial differential equations. In contrast to most engineering type control problems, the control does not continuously affect the system but only acts through the initial distributions. We briefly outline a general theory developed for dealing with such nonlinear systems by the use of a fixed point theorem. The problem and the theory are illustrated by some numerical simulations.

An identifiability analysis of a parent-metabolite pharmacokinetic model for ivabradine.

The paper considers the structural identifiability of a parent-metabolite pharmacokinetic model for ivabradine and one of its metabolites. The model, which is linear, is considered initially for intravenous administration of ivabradine, and then for a combined intravenous and oral administration. In both cases, the model is shown to be unidentifiable. Simplification of the model (for both forms of administration) to that proposed by Duffull et al. (1) results in a globally structurally identifiable model. The analysis could be applied to the modeling of any drug undergoing first-pass metabolism, with plasma concentrations available from drug and metabolite.

Fabrication of specimens of metamorphic magnetite crystals for field ion microscopy and atom probe microanalysis.

Field ion specimens have been successfully fabricated from samples of metamorphic magnetite crystals (Fe3O4) extracted from a polymetamorphosed, granulite-facies marble with the use of a focused ion beam. These magnetite crystals contain nanometer-scale, disk-shaped inclusions making this magnetite particularly attractive for investigating the capabilities of atom probe field ion microscopy (APFIM) for geological materials. Field ion microscope images of these magnetite crystals were obtained in which the observed size and morphology of the precipitates agree with previous results. Samples were analyzed in the energy compensated optical position-sensitive atom probe. Mass spectra were obtained in which peaks for singly ionized 16O, 56Fe and 56FeO and doubly ionized 54Fe, 56Fe and 57Fe peaks were fully resolved. Manganese and aluminum were observed in a limited analysis of a precipitate in an energy compensated position sensitive atom probe.

Extensions to a procedure for generating locally identifiable reparameterisations of unidentifiable systems.

In this paper extensions to an existing procedure for generating locally identifiable reparameterisations of unidentifiable systems are presented. These extensions further formalise the constructive nature of the methodology and lend themselves to application within symbolic manipulation packages. The extended reparameterisation procedure is described in detail and is illustrated with application to two known non-trivial examples of unidentifiable systems of practical relevance.