N-ethyl-N-nitrosourea (ENU) induced mutations within the Klotho gene lead to ectopic calcification and reduced lifespan in mouse models

Ectopic calcification (EC), which is the pathological deposition of calcium and phosphate in extra-skeletal tissues, may be associated with hypercalcaemic and hyperphosphataemic disorders, or it may occur in the absence of metabolic abnormalities. In addition, EC may be inherited as part of several monogenic disorders and studies of these have provided valuable insights into the metabolic pathways regulating mineral metabolism. For example, studies of tumoural calcinosis, a disorder characterised by hyperphosphataemia and progressive EC, have revealed mutations of fibroblast growth factor 23 (FGF23), polypeptide N-acetyl galactosaminyltransferase 3 (GALNT3) and klotho (KL), which are all part of a phosphate-regulating pathway. However, such studies in humans are limited by the lack of available large families with EC, and to facilitate such studies we assessed the progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) for EC. This identified two mutants with autosomal recessive forms of EC, and reduced lifespan, designated Ecalc1 and Ecalc2. Genetic mapping localized the Ecalc1 and Ecalc2 loci to a 11.0 Mb region on chromosome 5 that contained the klotho gene (Kl), and DNA sequence analysis identified nonsense (Gln203Stop) and missense (Ile604Asn) Kl mutations in Ecalc1 and Ecalc2 mice, respectively. The Gln203Stop mutation, located in KL1 domain, was severely hypomorphic and led to a 17-fold reduction of renal Kl expression. The Ile604Asn mutation, located in KL2 domain, was predicted to impair klotho protein stability and in vitro expression studies in COS-7 cells revealed endoplasmic reticulum retention of the Ile604Asn mutant. Further phenotype studies undertaken in Ecalc1 (kl203X/203X) mice demonstrated elevations in plasma concentrations of phosphate, FGF23 and 1,25-dihydroxyvitamin D. Thus, two allelic variants of Kl that develop EC and represent mouse models for tumoural calcinosis have been established. © 2015 Esapa et al.

Populations of models, experimental designs and coverage of parameter space by Latin Hypercube and Orthogonal Sampling

In this paper we have used simulations to make a conjecture about the coverage of a t-dimensional subspace of a d-dimensional parameter space of size n when performing k trials of Latin Hypercube sampling. This takes the form P(k,n,d,t) = 1 - e^(-k/n^(t-1)). We suggest that this coverage formula is independent of d and this allows us to make connections between building Populations of Models and Experimental Designs. We also show that Orthogonal sampling is superior to Latin Hypercube sampling in terms of allowing a more uniform coverage of the t-dimensional subspace at the sub-block size level. These ideas have particular relevance when attempting to perform uncertainty quantification and sensitivity analyses.

Revisiting the Task-Capability interface model for incorporating human factors into car-following models

Human factors such as distraction, fatigue, alcohol and drug use are generally ignored in car-following (CF) models. Such ignorance overestimates driver capability and leads to most CF models’ inability in realistically explaining human driving behaviors. This paper proposes a novel car-following modeling framework by introducing the difficulty of driving task measured as the dynamic interaction between driving task demand and driver capability. Task difficulty is formulated based on the famous Task Capability Interface (TCI) model, which explains the motivations behind driver’s decision making. The proposed method is applied to enhance two popular CF models: Gipps’ model and IDM, and named as TDGipps and TDIDM respectively. The behavioral soundness of TDGipps and TDIDM are discussed and their stabilities are analyzed. Moreover, the enhanced models are calibrated with the vehicle trajectory data, and validated to explain both regular and human factor influenced CF behavior (which is distraction caused by hand-held mobile phone conversation in this paper). Both the models show better performance than their predecessors, especially in presence of human factors.

Predicting surgical outcomes for paediatric spinal deformity patients using subject specific finite element models

INTRODUCTION Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. Our novel patient specific modelling software creates individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We have recently applied the model to investigate the influence of increasing magnitudes of surgically applied corrective force on predicted deformity correction...

Visible light-driven cross-coupling reactions at lower temperatures using a photocatalyst of palladium and gold alloy nanoparticles

Palladium (Pd)-catalyzed cross-coupling reactions are among the most important methods in organic synthesis. We report the discovery of highly efficient and green photocatalytic processes by which cross-coupling reactions, including Sonogashira, Stille, Hiyama, Ullmann, and Buchwald–Hartwig reactions, can be driven with visible light at temperatures slightly above room temperature using alloy nanoparticles of gold and Pd on zirconium oxide, thus achieving high yields. The alloy nanoparticles absorb visible light, and their conduction electrons gain energy, which is available at the surface Pd sites. Results of the density functional theory calculations indicate that transfer of the light excited electrons from the nanoparticle surface to the reactant molecules adsorbed on the nanoparticle surface activates the reactants. When the light intensity was increased, a higher reaction rate was observed, because of the increased population of photoexcited electrons. The irradiation wavelength also has an important impact on the reaction rates. Ultraviolet irradiation can drive some reactions with the chlorobenzene substrate, while visible light irradiation failed to, and substantially improve the yields of the reactions with the bromobenzene substrate. The discovery reveals the possibility of using low-energy and -density sources such as sunlight to drive chemical transformations.

Adjustable stiffness, external fixator for the rat femur osteotomy and segmental bone defect models

The mechanical environment around the healing of broken bone is very important as it determines the way the fracture will heal. Over the past decade there has been great clinical interest in improving bone healing by altering the mechanical environment through the fixation stability around the lesion. One constraint of preclinical animal research in this area is the lack of experimental control over the local mechanical environment within a large segmental defect as well as osteotomies as they heal. In this paper we report on the design and use of an external fixator to study the healing of large segmental bone defects or osteotomies. This device not only allows for controlled axial stiffness on the bone lesion as it heals, but it also enables the change of stiffness during the healing process in vivo. The conducted experiments have shown that the fixators were able to maintain a 5 mm femoral defect gap in rats in vivo during unrestricted cage activity for at least 8 weeks. Likewise, we observed no distortion or infections, including pin infections during the entire healing period. These results demonstrate that our newly developed external fixator was able to achieve reproducible and standardized stabilization, and the alteration of the mechanical environment of in vivo rat large bone defects and various size osteotomies. This confirms that the external fixation device is well suited for preclinical research investigations using a rat model in the field of bone regeneration and repair.

Copper nanoparticles on graphene support: An efficient photocatalyst for coupling of nitroaromatics in visible light

Copper is a low-cost plasmonic metal. Efficient photocatalysts of copper nanoparticles on graphene support are successfully developed for controllably catalyzing the coupling reactions of aromatic nitro compounds to the corresponding azoxy or azo compounds under visible-light irradiation. The coupling of nitrobenzene produces azoxybenzene with a yield of 90 % at 60 °C, but azobenzene with a yield of 96 % at 90 °C. When irradiated with natural sunlight (mean light intensity of 0.044 W cm−2) at about 35 °C, 70 % of the nitrobenzene is converted and 57 % of the product is azobenzene. The electrons of the copper nanoparticles gain the energy of the incident light through a localized surface plasmon resonance effect and photoexcitation of the bound electrons. The excited energetic electrons at the surface of the copper nanoparticles facilitate the cleavage of the NO bonds in the aromatic nitro compounds. Hence, the catalyzed coupling reaction can proceed under light irradiation and moderate conditions. This study provides a green photocatalytic route for the production of azo compounds and highlights a potential application for graphene.

Fractional diffusion models of cardiac electrical propagation: Role of structural heterogeneity in dispersion of repolarization

Impulse propagation in biological tissues is known to be modulated by structural heterogeneity. In cardiac muscle, improved understanding on how this heterogeneity influences electrical spread is key to advancing our interpretation of dispersion of repolarization. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a means of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, analysed against in vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies relevant characteristics of cardiac electrical propagation at tissue level. These include conduction effects on action potential (AP) morphology, the shortening of AP duration along the activation pathway and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media.

Mimicking breast cancer-induced bone metastasis in vivo: Current transplantation models and advanced humanized strategies

Bone metastasis is a complication that occurs in 80 % of women with advanced breast cancer. Despite the prevalence of bone metastatic disease, the avenues for its clinical management are still restricted to palliative treatment options. In fact, the underlying mechanisms of breast cancer osteotropism have not yet been fully elucidated due to a lack of suitable in vivo models that are able to recapitulate the human disease. In this work, we review the current transplantation-based models to investigate breast cancer-induced bone metastasis and delineate the strengths and limitations of the use of different grafting techniques, tissue sources, and hosts. We further show that humanized xenograft models incorporating human cells or tissue grafts at the primary tumor site or the metastatic site mimic more closely the human disease. Tissue-engineered constructs are emerging as a reproducible alternative to recapitulate functional humanized tissues in these murine models. The development of advanced humanized animal models may provide better platforms to investigate the mutual interactions between human cancer cells and their microenvironment and ultimately improve the translation of preclinical drug trials to the clinic.

Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

The root locus method: application to linear stability analysis and design of cooperative car-following models

A global framework for linear stability analyses of traffic models, based on the dispersion relation root locus method, is presented and is applied taking the example of a broad class of car-following (CF) models. This approach is able to analyse all aspects of the dynamics: long waves and short wave behaviours, phase velocities and stability features. The methodology is applied to investigate the potential benefits of connected vehicles, i.e. V2V communication enabling a vehicle to send and receive information to and from surrounding vehicles. We choose to focus on the design of the coefficients of cooperation which weights the information from downstream vehicles. The coefficients tuning is performed and different ways of implementing an efficient cooperative strategy are discussed. Hence, this paper brings design methods in order to obtain robust stability of traffic models, with application on cooperative CF models

Predictive models for water sources with high susceptibility for bromine-containing disinfection by-product formation: Implications for water treatment

This study examines a matrix of synthetic water samples designed to include conditions that favour brominated disinfection by-product (Br-DBP) formation, in order to provide predictive models suitable for high Br-DBP forming waters such as salinity-impacted waters. Br-DBPs are known to be more toxic than their chlorinated analogues, in general, and their formation may be favoured by routine water treatment practices such as coagulation/flocculation under specific conditions; therefore, circumstances surrounding their formation must be understood. The chosen factors were bromide concentration, mineral alkalinity, bromide to dissolved organic carbon (Br/DOC) ratio and Suwannee River natural organic matter concentration. The relationships between these parameters and DBP formation were evaluated by response surface modelling of data generated using a face-centred central composite experimental design. Predictive models for ten brominated and/or chlorinated DBPs are presented, as well as models for total trihalomethanes (tTHMs) and total dihaloacetonitriles (tDHANs), and bromide substitution factors for the THMs and DHANs classes. The relationships described revealed that increasing alkalinity and increasing Br/DOC ratio were associated with increasing bromination of THMs and DHANs, suggesting that DOC lowering treatment methods that do not also remove bromide such as enhanced coagulation may create optimal conditions for Br-DBP formation in waters in which bromide is present.

Stochastic dynamics of interacting haematopoietic stem cell niche lineages

Since we still know very little about stem cells in their natural environment, it is useful to explore their dynamics through modelling and simulation, as well as experimentally. Most models of stem cell systems are based on deterministic differential equations that ignore the natural heterogeneity of stem cell populations. This is not appropriate at the level of individual cells and niches, when randomness is more likely to affect dynamics. In this paper, we introduce a fast stochastic method for simulating a metapopulation of stem cell niche lineages, that is, many sub-populations that together form a heterogeneous metapopulation, over time. By selecting the common limiting timestep, our method ensures that the entire metapopulation is simulated synchronously. This is important, as it allows us to introduce interactions between separate niche lineages, which would otherwise be impossible. We expand our method to enable the coupling of many lineages into niche groups, where differentiated cells are pooled within each niche group. Using this method, we explore the dynamics of the haematopoietic system from a demand control system perspective. We find that coupling together niche lineages allows the organism to regulate blood cell numbers as closely as possible to the homeostatic optimum. Furthermore, coupled lineages respond better than uncoupled ones to random perturbations, here the loss of some myeloid cells. This could imply that it is advantageous for an organism to connect together its niche lineages into groups. Our results suggest that a potential fruitful empirical direction will be to understand how stem cell descendants communicate with the niche and how cancer may arise as a result of a failure of such communication.

Estimates on the coverage of parameter space using populations of models

In this paper we provide estimates for the coverage of parameter space when using Latin Hypercube Sampling, which forms the basis of building so-called populations of models. The estimates are obtained using combinatorial counting arguments to determine how many trials, k, are needed in order to obtain specified parameter space coverage for a given value of the discretisation size n. In the case of two dimensions, we show that if the ratio (Ø) of trials to discretisation size is greater than 1, then as n becomes moderately large the fractional coverage behaves as 1-exp-ø. We compare these estimates with simulation results obtained from an implementation of Latin Hypercube Sampling using MATLAB.