Specimen specific parameter identification of ovine lumbar intervertebral discs: On the influence of fibre-matrix and fibre-fibre shear interactions.

Numerical models of the intervertebral disc, which address mechanical questions commonly make use of the difference in water content between annulus and nucleus, and thus fluid and solid parts are separated. Despite this simplification, models remain complex due to the anisotropy and nonlinearity of the annulus and regional variations of the collagen fibre density. Additionally, it has been shown that cross-links make a large contribution to the stiffness of the annulus. Because of this complex composite structure, it is difficult to reproduce several sets of experimental data with one single set of material parameters. This study addresses the question to which extent the ultrastructure of the intervertebral disc should be modelled so that its moment-angle behaviour can be adequately described. Therefore, a hyperelastic constitutive law, based on continuum mechanical principles was derived, which does not only consider the anisotropy from the collagen fibres, but also interactions among the fibres and between the fibres and the ground substance. Eight ovine lumbar intervertebral discs were tested on a custom made spinal loading simulator in flexion/extension, lateral bending and axial rotation. Specimen-specific geometrical models were generated using CT images and T2 maps to distinguish between annulus fibrosus and nucleus pulposus. For the identification of the material parameters the annulus fibrosus was described with two scenarios: with and without fibre-matrix and fibre-fibre interactions. Both scenarios showed a similar behaviour on a load displacement level. Comparing model predictions to the experimental data, the mean RMS of all specimens and all load cases was 0.54±0.15° without the interaction and 0.54±0.19° when the fibre-matrix and fibre-fibre interactions were included. However, due to the increased stiffness when cross-links effects were included, this scenario showed more physiological stress-strain relations in uniaxial and biaxial stress states. Thus, the present study suggests that fibre-matrix and fibre-fibre interactions should be considered in the constitutive law when the model addresses questions concerning the stress field of the annulus fibrosus.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Object Kinetic Monte Carlo calculations of irradiated Fe-Cr dilute alloys: The effect of the interaction radius between substitutional Cr and self-interstitial Fe

ObjectKineticMonteCarlo models allow for the study of the evolution of the damage created by irradiation to time scales that are comparable to those achieved experimentally. Therefore, the essential ObjectKineticMonteCarlo parameters can be validated through comparison with experiments. However, this validation is not trivial since a large number of parameters is necessary, including migration energies of point defects and their clusters, binding energies of point defects in clusters, as well as the interactionradii. This is particularly cumbersome when describing an alloy, such as the Fe–Cr system, which is of interest for fusion energy applications. In this work we describe an ObjectKineticMonteCarlo model for Fe–Cr alloys in the dilute limit. The parameters used in the model come either from density functional theory calculations or from empirical interatomic potentials. This model is used to reproduce isochronal resistivity recovery experiments of electron irradiateddiluteFe–Cr alloys performed by Abe and Kuramoto. The comparison between the calculated results and the experiments reveal that an important parameter is the capture radius between substitutionalCr and self-interstitialFe atoms. A parametric study is presented on the effect of the capture radius on the simulated recovery curves.

Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients

Optimal sampling times are found for a study in which one of the primary purposes is to develop a model of the pharmacokinetics of itraconazole in patients with cystic fibrosis for both capsule and solution doses. The optimal design is expected to produce reliable estimates of population parameters for two different structural PK models. Data collected at these sampling times are also expected to provide the researchers with sufficient information to reasonably discriminate between the two competing structural models.

Interaction of metal salts with polyethers

The dielectric properties of pure low to medium molecular weight poly(ethylene glycol) and poly(propylene glycol) and a variety of their salt complexes have been studied through the measurement of the dielectric permittivity and dielectric loss over a range of frequency and temperature. The major proportion of this study has been concerned with the examination of the nature of the interaction between mercuric chloride and poly(propylene glycol) (PPG). Other salt-poly-ether combinations have also been considered such as cobalt chloride-PPG cadmium chloride-PPG zinc chloride-PPG and ferric chloride-PEG (polyethylene glycol). Some of this work was also supported by chemical shift and spin-lattice Nuclear Magnetic Resonance (N.M.R.) spectroscopy. The dielectric permittivity data were analysed using the Onsager relation to calculate the mean dipole moment per dipolar unit. This approach was employed in the discussion of various models proposed for the structure of salt-polyether complexes. The effect of mercuric chloride on the statistical conformations of poly(propylene-glycol) was studied in a quantitative manner using the relationships of marchal-Benoit. The dielectric relaxation activation energy and mean energy difference between gauche and trans conformations of poly(propylene glycol) in the presence of mercuric chloride, both showed a distinct minimum when the concentration of mercuric chloride was close to 5 mole %. Opposite behaviour was observed for the Cole-Cole parameter. It was concluded that the majority of the dielectric data could be rationalised in terms of a 5-membered cyclic complex formed between mercuric chloride and PPG in which the complexed segment of the polyether-(OMeCH2CH2O)- adopted either gauche or cis conformations.

Spatial and temporal interaction between sediment and microphytobenthos in a temperate estuarine macro-intertidal bay

Intertidal flats of the estuarine macro-intertidal Baie des Veys (France) were investigated to identify spatial features of sediment and microphytobenthos (MPB) in April 2003. Gradients occurred within the domain, and patches were identified close to vegetated areas or within the oyster-farming areas where calm physical conditions and biodeposition altered the sediment and MPB landscapes. Spatial patterns of chl a content were explained primarily by the influence of sediment features, while bed elevation and compaction brought only minor insights into MPB distribution regulation. The smaller size of MPB patches compared to silt patches revealed the interplay between physical structure defining the sediment landscape, the biotic patches that they contain, and that median grain-size is the most important parameter in explaining the spatial pattern of MPB. Small-scale temporal dynamics of sediment chl a content and grain-size distribution were surveyed in parallel during 2 periods of 14 d to detect tidal and seasonal variations. Our results showed a weak relationship between mud fraction and MPB biomass in March, and this relationship fully disappeared in July. Tidal exposure was the most important parameter in explaining the summer temporal dynamics of MPB. This study reveals the general importance of bed elevation and tidal exposure in muddy habitats and that silt content was a prime governing physical factor in winter. Biostabilisation processes seemed to behave only as secondary factors that could only amplify the initial silt accumulation in summer rather than primary factors explaining spatial or long-term trends of sediment changes.

Thermodynamic Study Of The Micellization Of Zwitterionic Surfactants And Their Interaction With Polymers In Water By Isothermal Titration Calorimetry.

The micellization of a homologous series of zwitterionic surfactants, a group of sulfobetaines, was studied using isothermal titration calorimetry (ITC) in the temperature range from 15 to 65 °C. The increase in both temperature and the alkyl chain length leads to more negative values of ΔGmic(0) , favoring the micellization. The entropic term (ΔSmic(0)) is predominant at lower temperatures, and above ca. 55-65 °C, the enthalpic term (ΔHmic(0)) becomes prevalent, figuring a jointly driven process as the temperature increases. The interaction of these sulfobetaines with different polymers was also studied by ITC. Among the polymers studied, only two induced the formation of micellar aggregates at lower surfactant concentration: poly(acrylic acid), PAA, probably due to the formation of hydrogen bonds between the carboxylic group of the polymer and the sulfonate group of the surfactant, and poly(sodium 4-styrenesulfonate), PSS, probably due to the incorporation of the hydrophobic styrene group into the micelles. The prevalence of the hydrophobic and not the electrostatic contributions to the interaction between sulfobetaine and PSS was confirmed by an increased interaction enthalpy in the presence of electrolytes (NaCl) and by the observation of a significant temperature dependence, the latter consistent with the proposed removal of hydrophobic groups from water.

High-resolution Transcript Profiling Of The Atypical Biotrophic Interaction Between Theobroma Cacao And The Fungal Pathogen Moniliophthora Perniciosa.

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

Clustering Of Water Bodies In Unpolluted And Polluted Environments Based On Escherichia Coli Phylogroup Abundance Using A Simple Interaction Database.

Different types of water bodies, including lakes, streams, and coastal marine waters, are often susceptible to fecal contamination from a range of point and nonpoint sources, and have been evaluated using fecal indicator microorganisms. The most commonly used fecal indicator is Escherichia coli, but traditional cultivation methods do not allow discrimination of the source of pollution. The use of triplex PCR offers an approach that is fast and inexpensive, and here enabled the identification of phylogroups. The phylogenetic distribution of E. coli subgroups isolated from water samples revealed higher frequencies of subgroups A1 and B23 in rivers impacted by human pollution sources, while subgroups D1 and D2 were associated with pristine sites, and subgroup B1 with domesticated animal sources, suggesting their use as a first screening for pollution source identification. A simple classification is also proposed based on phylogenetic subgroup distribution using the w-clique metric, enabling differentiation of polluted and unpolluted sites.

Human Mitochondrial Hsp70 (mortalin): Shedding Light On Atpase Activity, Interaction With Adenosine Nucleotides, Solution Structure And Domain Organization.

The human mitochondrial Hsp70, also called mortalin, is of considerable importance for mitochondria biogenesis and the correct functioning of the cell machinery. In the mitochondrial matrix, mortalin acts in the importing and folding process of nucleus-encoded proteins. The in vivo deregulation of mortalin expression and/or function has been correlated with age-related diseases and certain cancers due to its interaction with the p53 protein. In spite of its critical biological roles, structural and functional studies on mortalin are limited by its insoluble recombinant production. This study provides the first report of the production of folded and soluble recombinant mortalin when co-expressed with the human Hsp70-escort protein 1, but it is still likely prone to self-association. The monomeric fraction of mortalin presented a slightly elongated shape and basal ATPase activity that is higher than that of its cytoplasmic counterpart Hsp70-1A, suggesting that it was obtained in the functional state. Through small angle X-ray scattering, we assessed the low-resolution structural model of monomeric mortalin that is characterized by an elongated shape. This model adequately accommodated high resolution structures of Hsp70 domains indicating its quality. We also observed that mortalin interacts with adenosine nucleotides with high affinity. Thermally induced unfolding experiments indicated that mortalin is formed by at least two domains and that the transition is sensitive to the presence of adenosine nucleotides and that this process is dependent on the presence of Mg2+ ions. Interestingly, the thermal-induced unfolding assays of mortalin suggested the presence of an aggregation/association event, which was not observed for human Hsp70-1A, and this finding may explain its natural tendency for in vivo aggregation. Our study may contribute to the structural understanding of mortalin as well as to contribute for its recombinant production for antitumor compound screenings.

Reduced Lrp6 Expression And Increase In The Interaction Of Gsk3β With P53 Contribute To Podocyte Apoptosis In Diabetes Mellitus And Are Prevented By Green Tea.

In diabetes mellitus (DM), podocyte apoptosis leads to albuminuria and nephropathy progression. Low-density lipoprotein receptor-related protein 6 (LRP6) is WNT pathway receptor that is involved in podocyte death, adhesion and motility. Glycogen synthase kinase 3 (GSK3) interaction with p53 (GSK3-p53) promotes apoptosis in carcinoma cells. It is unknown if GSK3-p53 contributes to podocyte apoptosis in DM. In experimental DM, green tea (GT) reduces albuminuria by an unknown mechanism. In the present study, we assessed the role of the GSK3β-p53 in podocyte apoptosis and the effects of GT on these abnormalities. In diabetic spontaneously hypertensive rats (SHRs), GT prevents podocyte's p-LRP6 expression reduction, increased GSK3β-p53 and high p53 levels. In diabetic SHR rats, GT reduces podocyte apoptosis, foot process effacement and albuminuria. In immortalized mouse podocytes (iMPs), high glucose (HG), silencing RNA (siRNA) or blocking LRP6 (DKK-1) reduced p-LRP6 expression, leading to high GSK3β-p53, p53 expression, apoptosis and increased albumin influx. GSK3β blockade by BIO reduced GSK3β-p53 and podocyte apoptosis. In iMPs under HG, GT reduced apoptosis and the albumin influx by blocking GSK3β-p53 following the rise in p-LRP6 expression. These effects of GT were prevented by LRP6 siRNA or DKK-1. In conclusion, in DM, WNT inhibition, via LRP6, increases GSK3β-p53 and podocyte apoptosis. Maneuvers that inactivate GSK3β-p53, such as GT, may be renoprotective in DM.

Characterization Of A Protein-protein Interaction Network Of The Cbl-interacting Protein Kinase 8 From Sugarcane.

Plants are sessile organisms and have evolved to tolerate a constantly changing environment. After the onset of different stress conditions, calcineurin B-like (CBL) proteins can sense calcium signals and activate CBL-interacting protein kinase (CIPK) proteins, which can phosphorylate downstream proteins to reestablish plant homeostasis. Previous studies in the bioenergy crop sugarcane showed that the ScCIPK8 gene is induced by drought stress and is also related to sucrose content. Here, we have characterized the protein-protein interactions of ScCIPK8 with six CBL proteins (ScCBL1, ScCBL2, ScCBL3, ScCBL6, ScCBL9, and ScCBL10). Yeast two-hybrid assays showed that ScCIPK8 interacts with ScCBL1, ScCBL3, and ScCBL6. Bimolecular fluorescence complementation assays confirmed in planta the interactions that were observed in yeast cells. These findings give insights on the regulatory networks related to sugar accumulation and drought stress responses in sugarcane.

New Interaction Partners For Nek4.1 And Nek4.2 Isoforms: From The Dna Damage Response To Rna Splicing.

Neks are serine-threonine kinases that are similar to NIMA, a protein found in Aspergillus nidulans which is essential for cell division. In humans there are eleven Neks which are involved in different biological functions besides the cell cycle control. Nek4 is one of the largest members of the Nek family and has been related to the primary cilia formation and in DNA damage response. However, its substrates and interaction partners are still unknown. In an attempt to better understand the role of Nek4, we performed an interactomics study to find new biological processes in which Nek4 is involved. We also described a novel Nek4 isoform which lacks a region of 46 amino acids derived from an insertion of an Alu sequence and showed the interactomics profile of these two Nek4 proteins. Isoform 1 and isoform 2 of Nek4 were expressed in human cells and after an immunoprecipitation followed by mass spectrometry, 474 interacting proteins were identified for isoform 1 and 149 for isoform 2 of Nek4. About 68% of isoform 2 potential interactors (102 proteins) are common between the two Nek4 isoforms. Our results reinforce Nek4 involvement in the DNA damage response, cilia maintenance and microtubule stabilization, and raise the possibility of new functional contexts, including apoptosis signaling, stress response, translation, protein quality control and, most intriguingly, RNA splicing. We show for the first time an unexpected difference between both Nek4 isoforms in RNA splicing control. Among the interacting partners, we found important proteins such as ANT3, Whirlin, PCNA, 14-3-3ε, SRSF1, SRSF2, SRPK1 and hNRNPs proteins. This study provides new insights into Nek4 functions, identifying new interaction partners and further suggests an interesting difference between isoform 1 and isoform 2 of this kinase. Nek4 isoform 1 may have similar roles compared to other Neks and these roles are not all preserved in isoform 2. Besides, in some processes, both isoforms showed opposite effects, indicating a possible fine controlled regulation.