Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Mechanisms of species divergence through visual adaptation and sexual selection: Perspectives from a cichlid model system

The theory of ecological speciation suggests that assortative mating evolves most easily when mating preferences are;directly linked to ecological traits that are subject to divergent selection. Sensory adaptation can play a major role in this process,;because selective mating is often mediated by sexual signals: bright colours, complex song, pheromone blends and so on. When;divergent sensory adaptation affects the perception of such signals, mating patterns may change as an immediate consequence.;Alternatively, mating preferences can diverge as a result of indirect effects: assortative mating may be promoted by selection;against intermediate phenotypes that are maladapted to their (sensory) environment. For Lake Victoria cichlids, the visual environment;constitutes an important selective force that is heterogeneous across geographical and water depth gradients. We investigate;the direct and indirect effects of this heterogeneity on the evolution of female preferences for alternative male nuptial colours;(red and blue) in the genus Pundamilia. Here, we review the current evidence for divergent sensory drive in this system, extract;general principles, and discuss future perspectives

Use of the Supermolecule Approach To Model the Syn and Anti Conformations of Solvated Cyclic 3‘,5‘-Adenosine Monophosphate

The supermolecule approach has been used to model the hydration of cyclic 3‘,5‘-adenosine monophosphate, cAMP. Model building combined with PM3 optimizations predict that the anti conformer of cAMP is capable of hydrogen bonding to an additional solvent water molecule compared to the syn conformer. The addition of one water to the syn superstructure with concurrent rotation of the base about the glycosyl bond to form the anti superstructure leads to an additional enthalpy of stabilization of approximately −6 kcal/mol at the PM3 level. This specific solute−solvent interaction is an example of a large solvent effect, as the method predicts that cAMP has a conformational preference for the anti isomer in solution. This conformational preference results from a change in the number of specific solute−solvent interactions in this system. This prediction could be tested by NMR techniques. The number of waters predicted to be in the first hydration sphere around cAMP is in agreement with the results of hydration studies of nucleotides in DNA. In addition, the detailed picture of solvation about this cyclic nucleotide is in agreement with infrared experimental results.

Hydraulic Conductivity of Model Soil-Bentonite Backfills Subjected to Wet-Dry Cycling

The potential for changes in hydraulic conductivity, k, of two model soil-bentonite (SB) backfills subjected to wet-dry cycling was investigated. The backfills were prepared with the same base soil (clean, fine sand) but different bentonite contents (2.7 and 5.6 dry wt %). Saturation (S), volume change, and k of consolidated backfill specimens (effective stress = 24 kPa) were evaluated over three to seven cycles in which the matric suction, Ym, in the drying stage ranged from 50 to 700 kPa. Both backfills exhibited susceptibility to degradation in k caused by wet-dry cycling. Mean values of k for specimens dried at Ym = 50 kPa (S = 30-60 % after drying) remained low after two cycles, but increased by 5- to 300-fold after three or more cycles. Specimens dried at Ym ≥ 150 kPa (S < 30 % after drying) were less resilient and exhibited 500- to 10 000-fold increases in k after three or more cycles. The greater increases in k for these specimens correlated with greater vertical shrinkage upon drying. The findings suggest that increases in hydraulic conductivity due to wet-dry cycling may be a concern for SB vertical barriers located within the zone of a fluctuating groundwater table.

Tracing of 15N in a mountain spruce forest: comparison of experimental data with th TRACE model

Despite numerous studies about nitrogen-cycling in forest ecosystems, many uncertainties remain, especially regarding the longer-term nitrogen accumulation. To contribute to filling this gap, the dynamic process-based model TRACE, with the ability to simulate 15N tracer redistribution in forest ecosystems was used to study N cycling processes in a mountain spruce forest of the northern edge of the Alps in Switzerland (Alptal, SZ). Most modeling analyses of N-cycling and C-N interactions have very limited ability to determine whether the process interactions are captured correctly. Because the interactions in such a system are complex, it is possible to get the whole-system C and N cycling right in a model without really knowing if the way the model combines fine-scale interactions to derive whole-system cycling is correct. With the possibility to simulate 15N tracer redistribution in ecosystem compartments, TRACE features a very powerful tool for the validation of fine-scale processes captured by the model. We first adapted the model to the new site (Alptal, Switzerland; long-term low-dose N-amendment experiment) by including a new algorithm for preferential water flow and by parameterizing of differences in drivers such as climate, N deposition and initial site conditions. After the calibration of key rates such as NPP and SOM turnover, we simulated patterns of 15N redistribution to compare against 15N field observations from a large-scale labeling experiment. The comparison of 15N field data with the modeled redistribution of the tracer in the soil horizons and vegetation compartments shows that the majority of fine-scale processes are captured satisfactorily. Particularly, the model is able to reproduce the fact that the largest part of the N deposition is immobilized in the soil. The discrepancies of 15N recovery in the LF and M soil horizon can be explained by the application method of the tracer and by the retention of the applied tracer by the well developed moss layer, which is not considered in the model. Discrepancies in the dynamics of foliage and litterfall 15N recovery were also observed and are related to the longevity of the needles in our mountain forest. As a next step, we will use the final Alptal version of the model to calculate the effects of climate change (temperature, CO2) and N deposition on ecosystem C sequestration in this regionally representative Norway spruce (Picea abies) stand.

Neuronal expression of the ubiquitin ligase Nedd4-2 in rat dorsal root ganglia: modulation in the spared nerve injury model of neuropathic pain

Neuronal hyperexcitability following peripheral nerve lesions may stem from altered activity of voltage-gated sodium channels (VGSCs), which gives rise to allodynia or hyperalgesia. In vitro, the ubiquitin ligase Nedd4-2 is a negative regulator of VGSC α-subunits (Na(v)), in particular Na(v)1.7, a key actor in nociceptor excitability. We therefore studied Nedd4-2 in rat nociceptors, its co-expression with Na(v)1.7 and Na(v)1.8, and its regulation in pathology. Adult rats were submitted to the spared nerve injury (SNI) model of neuropathic pain or injected with complete Freund's adjuvant (CFA), a model of inflammatory pain. L4 dorsal root ganglia (DRG) were analyzed in sham-operated animals, seven days after SNI and 48 h after CFA with immunofluorescence and Western blot. We observed Nedd4-2 expression in almost 50% of DRG neurons, mostly small and medium-sized. A preponderant localization is found in the non-peptidergic sub-population. Additionally, 55.7 ± 2.7% and 55.0 ± 3.6% of Nedd4-2-positive cells are co-labeled with Na(v)1.7 and Na(v)1.8 respectively. SNI significantly decreases the proportion of Nedd4-2-positive neurons from 45.9 ± 1.9% to 33.5 ± 0.7% (p<0.01) and the total Nedd4-2 protein to 44% ± 0.13% of its basal level (p<0.01, n=4 animals in each group, mean ± SEM). In contrast, no change in Nedd4-2 was found after peripheral inflammation induced by CFA. These results indicate that Nedd4-2 is present in nociceptive neurons, is downregulated after peripheral nerve injury, and might therefore contribute to the dysregulation of Na(v)s involved in the hyperexcitability associated with peripheral nerve injuries.

Maternal mental health in the first 3-week postpartum: the impact of caregiver support and the subjective experience of childbirth - a longitudinal path model

Acute stress reactions (ASR) and postpartum depressive symptoms (PDS) are frequent after childbirth. The present study addresses the change and overlap of ASR and PDS from the 1- to 3-week postpartum and examines the interplay of caregiver support and subjective birth experience with regard to the development of ASR/PDS within a longitudinal path model.

Rapid climate change during the Weichselian Lateglacial in Ireland: Chironomid-inferred summer temperatures from Fiddaun, Co. Galway

A lacustrine sediment core from Fiddaun, western Ireland was studied to reconstruct summer temperature changes during the Weichselian Lateglacial. This site is located close to the Atlantic Ocean; and so is potentially sensitive to climatic changes associated with changes in ocean circulation. The record, comprising the end of the Weichselian Pleniglacial to the early Holocene, was analysed for fossil chironomids, lithology, and oxygen and carbon isotopes in the sedimentary carbonates. These proxies clearly show rapid warming at the onset of the Lateglacial Interstadial, relatively high summer temperatures during the Interstadial, pronounced cooling during the Younger Dryas, and subsequent warming at the transition to the Holocene. Chironomid-inferred mean July air temperatures for the Interstadial are ~12.5–14.5 °C, ~7.5 °C for the Younger Dryas, and ~15.0 °C for the early Holocene. Furthermore, this research provides evidence for at least two cold events during the Interstadial. These more moderate temperature oscillations can be correlated to Greenland Interstadial events 1b and 1d, on the basis of the age-depth model for the Fiddaun sequence. Based on multiple proxies, the first cold oscillation (GI-1d) was the more severe of the two in Ireland.

A physics-based correction model for homogenizing sub-daily temperature series

A new physics-based technique for correcting inhomogeneities present in sub-daily temperature records is proposed. The approach accounts for changes in the sensor-shield characteristics that affect the energy balance dependent on ambient weather conditions (radiation, wind). An empirical model is formulated that reflects the main atmospheric processes and can be used in the correction step of a homogenization procedure. The model accounts for short- and long-wave radiation fluxes (including a snow cover component for albedo calculation) of a measurement system, such as a radiation shield. One part of the flux is further modulated by ventilation. The model requires only cloud cover and wind speed for each day, but detailed site-specific information is necessary. The final model has three free parameters, one of which is a constant offset. The three parameters can be determined, e.g., using the mean offsets for three observation times. The model is developed using the example of the change from the Wild screen to the Stevenson screen in the temperature record of Basel, Switzerland, in 1966. It is evaluated based on parallel measurements of both systems during a sub-period at this location, which were discovered during the writing of this paper. The model can be used in the correction step of homogenization to distribute a known mean step-size to every single measurement, thus providing a reasonable alternative correction procedure for high-resolution historical climate series. It also constitutes an error model, which may be applied, e.g., in data assimilation approaches.

A Cox Model for Biostatistics of the Future

Professor Sir David R. Cox (DRC) is widely acknowledged as among the most important scientists of the second half of the twentieth century. He inherited the mantle of statistical science from Pearson and Fisher, advanced their ideas, and translated statistical theory into practice so as to forever change the application of statistics in many fields, but especially biology and medicine. The logistic and proportional hazards models he substantially developed, are arguably among the most influential biostatistical methods in current practice. This paper looks forward over the period from DRC's 80th to 90th birthdays, to speculate about the future of biostatistics, drawing lessons from DRC's contributions along the way. We consider "Cox's model" of biostatistics, an approach to statistical science that: formulates scientific questions or quantities in terms of parameters gamma in probability models f(y; gamma) that represent in a parsimonious fashion, the underlying scientific mechanisms (Cox, 1997); partition the parameters gamma = theta, eta into a subset of interest theta and other "nuisance parameters" eta necessary to complete the probability distribution (Cox and Hinkley, 1974); develops methods of inference about the scientific quantities that depend as little as possible upon the nuisance parameters (Barndorff-Nielsen and Cox, 1989); and thinks critically about the appropriate conditional distribution on which to base infrences. We briefly review exciting biomedical and public health challenges that are capable of driving statistical developments in the next decade. We discuss the statistical models and model-based inferences central to the CM approach, contrasting them with computationally-intensive strategies for prediction and inference advocated by Breiman and others (e.g. Breiman, 2001) and to more traditional design-based methods of inference (Fisher, 1935). We discuss the hierarchical (multi-level) model as an example of the future challanges and opportunities for model-based inference. We then consider the role of conditional inference, a second key element of the CM. Recent examples from genetics are used to illustrate these ideas. Finally, the paper examines causal inference and statistical computing, two other topics we believe will be central to biostatistics research and practice in the coming decade. Throughout the paper, we attempt to indicate how DRC's work and the "Cox Model" have set a standard of excellence to which all can aspire in the future.

Protein buffering in model systems and in whole human saliva

The aim of this study was to quantify the buffer attributes (value, power, range and optimum) of two model systems for whole human resting saliva, the purified proteins from whole human resting saliva and single proteins. Two model systems, the first containing amyloglucosidase and lysozyme, and the second containing amyloglucosidase and alpha-amylase, were shown to provide, in combination with hydrogencarbonate and di-hydrogenphosphate, almost identical buffer attributes as whole human resting saliva. It was further demonstrated that changes in the protein concentration as small as 0.1% may change the buffer value of a buffer solution up to 15 times. Additionally, it was shown that there was a protein concentration change in the same range (0.16%) between saliva samples collected at the time periods of 13:00 and others collected at 9:00 am and 17:00. The mode of the protein expression changed between these samples corresponded to the change in basic buffer power and the change of the buffer value at pH 6.7. Finally, SDS Page and Ruthenium II tris (bathophenantroline disulfonate) staining unveiled a constant protein expression in all samples except for one 50 kDa protein band. As the change in the expression pattern of that 50 kDa protein band corresponded to the change in basic buffer power and the buffer value at pH 6.7, it was reasonable to conclude that this 50 kDa protein band may contain the protein(s) belonging to the protein buffer system of human saliva.

A coronary heart disease risk model for predicting the effect of potent antiretroviral therapy in HIV-1 infected men

BACKGROUND: Many HIV-infected patients on highly active antiretroviral therapy (HAART) experience metabolic complications including dyslipidaemia and insulin resistance, which may increase their coronary heart disease (CHD) risk. We developed a prognostic model for CHD tailored to the changes in risk factors observed in patients starting HAART. METHODS: Data from five cohort studies (British Regional Heart Study, Caerphilly and Speedwell Studies, Framingham Offspring Study, Whitehall II) on 13,100 men aged 40-70 and 114,443 years of follow up were used. CHD was defined as myocardial infarction or death from CHD. Model fit was assessed using the Akaike Information Criterion; generalizability across cohorts was examined using internal-external cross-validation. RESULTS: A parametric model based on the Gompertz distribution generalized best. Variables included in the model were systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, triglyceride, glucose, diabetes mellitus, body mass index and smoking status. Compared with patients not on HAART, the estimated CHD hazard ratio (HR) for patients on HAART was 1.46 (95% CI 1.15-1.86) for moderate and 2.48 (95% CI 1.76-3.51) for severe metabolic complications. CONCLUSIONS: The change in the risk of CHD in HIV-infected men starting HAART can be estimated based on typical changes in risk factors, assuming that HRs estimated using data from non-infected men are applicable to HIV-infected men. Based on this model the risk of CHD is likely to increase, but increases may often be modest, and could be offset by lifestyle changes.

Mechanisms of intrinsic beating variability in cardiac cell cultures and model pacemaker networks

Heart rate variability (HRV) exhibits fluctuations characterized by a power law behavior of its power spectrum. The interpretation of this nonlinear HRV behavior, resulting from interactions between extracardiac regulatory mechanisms, could be clinically useful. However, the involvement of intrinsic variations of pacemaker rate in HRV has scarcely been investigated. We examined beating variability in spontaneously active incubating cultures of neonatal rat ventricular myocytes using microelectrode arrays. In networks of mathematical model pacemaker cells, we evaluated the variability induced by the stochastic gating of transmembrane currents and of calcium release channels and by the dynamic turnover of ion channels. In the cultures, spontaneous activity originated from a mobile focus. Both the beat-to-beat movement of the focus and beat rate variability exhibited a power law behavior. In the model networks, stochastic fluctuations in transmembrane currents and stochastic gating of calcium release channels did not reproduce the spatiotemporal patterns observed in vitro. In contrast, long-term correlations produced by the turnover of ion channels induced variability patterns with a power law behavior similar to those observed experimentally. Therefore, phenomena leading to long-term correlated variations in pacemaker cellular function may, in conjunction with extracardiac regulatory mechanisms, contribute to the nonlinear characteristics of HRV.

The spatiotemporal pattern of auditory cortical responses during verbal hallucinations

Functional magnetic resonance imaging (fMRI) studies can provide insight into the neural correlates of hallucinations. Commonly, such studies require self-reports about the timing of the hallucination events. While many studies have found activity in higher-order sensory cortical areas, only a few have demonstrated activity of the primary auditory cortex during auditory verbal hallucinations. In this case, using self-reports as a model of brain activity may not be sensitive enough to capture all neurophysiological signals related to hallucinations. We used spatial independent component analysis (sICA) to extract the activity patterns associated with auditory verbal hallucinations in six schizophrenia patients. SICA decomposes the functional data set into a set of spatial maps without the use of any input function. The resulting activity patterns from auditory and sensorimotor components were further analyzed in a single-subject fashion using a visualization tool that allows for easy inspection of the variability of regional brain responses. We found bilateral auditory cortex activity, including Heschl's gyrus, during hallucinations of one patient, and unilateral auditory cortex activity in two more patients. The associated time courses showed a large variability in the shape, amplitude, and time of onset relative to the self-reports. However, the average of the time courses during hallucinations showed a clear association with this clinical phenomenon. We suggest that detection of this activity may be facilitated by examining hallucination epochs of sufficient length, in combination with a data-driven approach.