Theory of commensurable magnetic structures in holmium

The tendency for the period of the helically ordered moments in holmium to lock into values which are commensurable with the lattice is studied theoretically as a function of temperature and magnetic field. The commensurable effects are derived in the mean-field approximation from numerical calculations of the free energy of various commensurable structures, and the results are compared with the extensive experimental evidence collected during the last ten years on the magnetic structures in holmium. In general the stability of the different commensurable structures is found to be in accord with the experiments, except for the tau=5/18 structure observed a few degrees below T_N in a b-axis field. The trigonal coupling recently detected in holmium is found to be the interaction required to explain the increased stability of the tau=1/5 structure around 42 K, and of the tau=1/4 structure around 96 K, when a field is applied along the c-axis.Comment: REVTEX, 31 pages, 7 postscript figure

Mathematical modelling of contact dermatitis from nickel and chromium

Dermal exposure to metal allergens can lead to irritant and allergic contact dermatitis (ACD). In this paper we present a mathematical model of the absorption of metal ions, hexavalent chromium and nickel, into the viable epidermis and compare the localised irritant and T-lymphocyte (T-cell) mediated immune responses. The model accounts for the spatial-temporal variation of skin health, extra and intracellular allergen concentrations, innate immune cells, T-cells, cytokine signalling and lymph node activity up to about 6 days after contact with these metals; repair processes associated with withdrawal of exposure to both metals is not considered in the current model, being assumed secondary during the initial phases of exposure. Simulations of the resulting system of PDEs are studied in one-dimension, i.e. across skin depth, and three-dimensional scenarios with the aim of comparing the responses to the two ions in the cases of first contact (no T-cells initially present) and second contact (T-cells initially present). The results show that on continuous contact, chromium ions elicit stronger skin inflammation, but for nickel, subsequent re-exposure stimulates stronger responses due to an accumulation of cytotoxic T-cell mediated responses which characterise ACD. Furthermore, the surface area of contact to these metals has little effect on the speed of response, whilst sensitivity is predicted to increase with the thickness of skin. The modelling approach is generic and should be applicable to describe contact dermatitis from a wide range of allergens

Survey of US public attitudes toward pharmacogenetic testing

To assess public attitudes and interest in pharmacogenetic (PGx) testing, we conducted a random-digit-dial telephone survey of U.S. adults, achieving a response rate of 42% (n=1139). Most respondents expressed interest in PGx testing to predict mild or serious side effects (73% ±3.29% and 85% ±2.91%, respectively), guide dosing (91%) and assist with drug selection (92%). Younger individuals (ages 18–34) were more likely to be interested in PGx testing to predict serious side effects (vs. ages 55+), as well as Whites, those with a college degree, and who had experienced side effects from medications. However, most respondents (78% ±3.14%) were not likely to have a PGx test if there was a risk that their DNA sample or test result could be shared without their permission. Given differences in interest among some groups, providers should clearly discuss the purpose of testing, alternative testing options (if available), and policies to protect patient privacy and confidentiality

Survey of genetic counselors and clinical geneticists' use and attitudes toward pharmacogenetic testing

Pharmacogenetic (PGx) testing aims to improve therapeutic outcomes through tailoring treatment based on a patient’s genetic risk for non-response and/or an adverse event. Given their expertise, geneticists could facilitate the use of PGx testing; however, , the preparedness and perceived role of the clinical genetics community is unclear. To assess the attitudes, preparedness, and perceived roles of geneticists in the delivery of PGx testing, we conducted a survey of 1500 randomly selected board-certified genetic counselors and clinical geneticists in the U.S (response rate: 37.8% (n=516)). Twelve percent of genetic counselors and 41% of clinical geneticists indicated that they had ordered or coordinated patient care for PGx testing, a seemingly high proportion at this early stage of adoption. Almost all respondents had some education on pharmacogenetics, though only 28% of counselors and 58% of clinical geneticists indicated they felt well-informed about PGx testing. About half of counselors (52%) and clinical geneticists (46%) felt they would play ‘some’ role in the delivery of PGx testing; 17% and 19%, respectively, felt that they would play ‘no’ or ‘a little’ role. At this early stage of PGx testing, the role of geneticists and genetic counselors is unclear. However, their experience may aid in readying PGx testing and informing delivery strategies into clinical practice

Adaptive response and enlargement of dynamic range

Many membrane channels and receptors exhibit adaptive, or desensitized, response to a strong sustained input stimulus, often supported by protein activity-dependent inactivation. Adaptive response is thought to be related to various cellular functions such as homeostasis and enlargement of dynamic range by background compensation. Here we study the quantitative relation between adaptive response and background compensation within a modeling framework. We show that any particular type of adaptive response is neither sufficient nor necessary for adaptive enlargement of dynamic range. In particular a precise adaptive response, where system activity is maintained at a constant level at steady state, does not ensure a large dynamic range neither in input signal nor in system output. A general mechanism for input dynamic range enlargement can come about from the activity-dependent modulation of protein responsiveness by multiple biochemical modification, regardless of the type of adaptive response it induces. Therefore hierarchical biochemical processes such as methylation and phosphorylation are natural candidates to induce this property in signaling systems.Comment: Corrected typos, minor text revision

Characterization of growth and metabolism of the haloalkaliphile Natronomonas pharaonis

Natronomonas pharaonis is an archaeon adapted to two extreme conditions: high salt concentration and alkaline pH. It has become one of the model organisms for the study of extremophilic life. Here, we present a genome-scale, manually curated metabolic reconstruction for the microorganism. The reconstruction itself represents a knowledge base of the haloalkaliphile's metabolism and, as such, would greatly assist further investigations on archaeal pathways. In addition, we experimentally determined several parameters relevant to growth, including a characterization of the biomass composition and a quantification of carbon and oxygen consumption. Using the metabolic reconstruction and the experimental data, we formulated a constraints-based model which we used to analyze the behavior of the archaeon when grown on a single carbon source. Results of the analysis include the finding that Natronomonas pharaonis, when grown aerobically on acetate, uses a carbon to oxygen consumption ratio that is theoretically near-optimal with respect to growth and energy production. This supports the hypothesis that, under simple conditions, the microorganism optimizes its metabolism with respect to the two objectives. We also found that the archaeon has a very low carbon efficiency of only about 35%. This inefficiency is probably due to a very low P/O ratio as well as to the other difficulties posed by its extreme environment

Overview of mathematical approaches used to model bacterial chemotaxis II: bacterial populations

We review the application of mathematical modeling to understanding the behavior of populations of chemotactic bacteria. The application of continuum mathematical models, in particular generalized Keller–Segel models, is discussed along with attempts to incorporate the microscale (individual) behavior on the macroscale, modeling the interaction between different species of bacteria, the interaction of bacteria with their environment, and methods used to obtain experimentally verified parameter values. We allude briefly to the role of modeling pattern formation in understanding collective behavior within bacterial populations. Various aspects of each model are discussed and areas for possible future research are postulated

Overview of mathematical approaches used to model bacterial chemotaxis I: the single cell

Mathematical modeling of bacterial chemotaxis systems has been influential and insightful in helping to understand experimental observations. We provide here a comprehensive overview of the range of mathematical approaches used for modeling, within a single bacterium, chemotactic processes caused by changes to external gradients in its environment. Specific areas of the bacterial system which have been studied and modeled are discussed in detail, including the modeling of adaptation in response to attractant gradients, the intracellular phosphorylation cascade, membrane receptor clustering, and spatial modeling of intracellular protein signal transduction. The importance of producing robust models that address adaptation, gain, and sensitivity are also discussed. This review highlights that while mathematical modeling has aided in understanding bacterial chemotaxis on the individual cell scale and guiding experimental design, no single model succeeds in robustly describing all of the basic elements of the cell. We conclude by discussing the importance of this and the future of modeling in this area

Mathematical description of bacterial traveling pulses

The Keller-Segel system has been widely proposed as a model for bacterial waves driven by chemotactic processes. Current experiments on {\em E. coli} have shown precise structure of traveling pulses. We present here an alternative mathematical description of traveling pulses at a macroscopic scale. This modeling task is complemented with numerical simulations in accordance with the experimental observations. Our model is derived from an accurate kinetic description of the mesoscopic run-and-tumble process performed by bacteria. This model can account for recent experimental observations with {\em E. coli}. Qualitative agreements include the asymmetry of the pulse and transition in the collective behaviour (clustered motion versus dispersion). In addition we can capture quantitatively the main characteristics of the pulse such as the speed and the relative size of tails. This work opens several experimental and theoretical perspectives. Coefficients at the macroscopic level are derived from considerations at the cellular scale. For instance the stiffness of the signal integration process turns out to have a strong effect on collective motion. Furthermore the bottom-up scaling allows to perform preliminary mathematical analysis and write efficient numerical schemes. This model is intended as a predictive tool for the investigation of bacterial collective motion

A combined model reduction algorithm for controlled biochemical systems

Background: Systems Biology continues to produce increasingly large models of complex biochemical reaction networks. In applications requiring, for example, parameter estimation, the use of agent-based modelling approaches, or real-time simulation, this growing model complexity can present a significant hurdle. Often, however, not all portions of a model are of equal interest in a given setting. In such situations methods of model reduction offer one possible approach for addressing the issue of complexity by seeking to eliminate those portions of a pathway that can be shown to have the least effect upon the properties of interest. Methods: In this paper a model reduction algorithm bringing together the complementary aspects of proper lumping and empirical balanced truncation is presented. Additional contributions include the development of a criterion for the selection of state-variable elimination via conservation analysis and use of an ‘averaged’ lumping inverse. This combined algorithm is highly automatable and of particular applicability in the context of ‘controlled’ biochemical networks. Results: The algorithm is demonstrated here via application to two examples; an 11 dimensional model of bacterial chemotaxis in Escherichia coli and a 99 dimensional model of extracellular regulatory kinase activation (ERK) mediated via the epidermal growth factor (EGF) and nerve growth factor (NGF) receptor pathways. In the case of the chemotaxis model the algorithm was able to reduce the model to 2 state-variables producing a maximal relative error between the dynamics of the original and reduced models of only 2.8% whilst yielding a 26 fold speed up in simulation time. For the ERK activation model the algorithm was able to reduce the system to 7 state-variables, incurring a maximal relative error of 4.8%, and producing an approximately 10 fold speed up in the rate of simulation. Indices of controllability and observability are additionally developed and demonstrated throughout the paper. These provide insight into the relative importance of individual reactants in mediating a biochemical system’s input-output response even for highly complex networks. Conclusions: Through application, this paper demonstrates that combined model reduction methods can produce a significant simplification of complex Systems Biology models whilst retaining a high degree of predictive accuracy. In particular, it is shown that by combining the methods of proper lumping and empirical balanced truncation it is often possible to produce more accurate reductions than can be obtained by the use of either method in isolation