COMMIT: Convex Optimization Modeling for Microstructure Informed Tractography.

Tractography is a class of algorithms aiming at in vivo mapping the major neuronal pathways in the white matter from diffusion magnetic resonance imaging (MRI) data. These techniques offer a powerful tool to noninvasively investigate at the macroscopic scale the architecture of the neuronal connections of the brain. However, unfortunately, the reconstructions recovered with existing tractography algorithms are not really quantitative even though diffusion MRI is a quantitative modality by nature. As a matter of fact, several techniques have been proposed in recent years to estimate, at the voxel level, intrinsic microstructural features of the tissue, such as axonal density and diameter, by using multicompartment models. In this paper, we present a novel framework to reestablish the link between tractography and tissue microstructure. Starting from an input set of candidate fiber-tracts, which are estimated from the data using standard fiber-tracking techniques, we model the diffusion MRI signal in each voxel of the image as a linear combination of the restricted and hindered contributions generated in every location of the brain by these candidate tracts. Then, we seek for the global weight of each of them, i.e., the effective contribution or volume, such that they globally fit the measured signal at best. We demonstrate that these weights can be easily recovered by solving a global convex optimization problem and using efficient algorithms. The effectiveness of our approach has been evaluated both on a realistic phantom with known ground-truth and in vivo brain data. Results clearly demonstrate the benefits of the proposed formulation, opening new perspectives for a more quantitative and biologically plausible assessment of the structural connectivity of the brain.

An introductory review on gravitational-deformation induced structures, fabrics and modeling

Recent studies have pointed out a similarity between tectonics and slope tectonic-induced structures. Numerous studies have demonstrated that structures and fabrics previously interpreted as of purely geodynamical origin are instead the result of large slope deformation, and this led in the past to erroneous interpretations. Nevertheless, their limit seems not clearly defined, but it is somehow transitional. Some studies point out continuity between failures developing at surface with upper crust movements. In this contribution, the main studies which examine the link between rock structures and slope movements are reviewed. The aspects regarding model and scale of observation are discussed together with the role of pre-existing weaknesses in the rock mass. As slope failures can develop through progressive failure, structures and their changes in time and space can be recognized. Furthermore, recognition of the origin of these structures can help in avoiding misinterpretations of regional geology. This also suggests the importance of integrating different slope movement classifications based on distribution and pattern of deformation and the application of structural geology techniques. A structural geology approach in the landslide community is a tool that can greatly support the hazard quantification and related risks, because most of the physical parameters, which are used for landslide modeling, are derived from geotechnical tests or the emerging geophysical approaches.

Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

Historical review: Negative efficacy and the constitutive activity of G-protein-coupled receptors.

The idea that a receptor can produce signalling without agonist intervention and that several antagonists can be 'active' in repressing such spontaneous activity is contained in the concept of ligand-induced conformational changes. Yet, this idea was neglected by pharmacologists for many years. In this article, we review the events that brought inverse agonism and constitutive activity to general attention and made this phenomenon a topic of current research. We also suggest a classification of antagonists based on the cooperativity that links their primary site of interaction with other functional domains of the receptor.

GLUT2 surface expression and intracellular transport via the constitutive pathway in pancreatic beta cells and insulinoma: evidence for a block in trans-Golgi network exit by brefeldin A.

The biosynthesis, intracellular transport, and surface expression of the beta cell glucose transporter GLUT2 was investigated in isolated islets and insulinoma cells. Using a trypsin sensitivity assay to measure cell surface expression, we determined that: (a) greater than 95% of GLUT2 was expressed on the plasma membrane; (b) GLUT2 did not recycle in intracellular vesicles; and (c) after trypsin treatment, reexpression of the intact transporter occurred with a t1/2 of approximately 7 h. Kinetics of intracellular transport of GLUT2 was investigated in pulse-labeling experiments combined with glycosidase treatment and the trypsin sensitivity assay. We determined that transport from the endoplasmic reticulum to the trans-Golgi network (TGN) occurred with a t1/2 of 15 min and that transport from the TGN to the plasma membrane required a similar half-time. When added at the start of a pulse-labeling experiment, brefeldin A prevented exit of GLUT2 from the endoplasmic reticulum. When the transporter was first accumulated in the TGN during a 15-min period of chase, but not following a low temperature (22 degrees C) incubation, addition of brefeldin A (BFA) prevented subsequent surface expression of the transporter. This indicated that brefeldin A prevented GLUT2 exit from the TGN by acting at a site proximal to the 22 degrees C block. Together, these data demonstrate that GLUT2 surface expression in beta cells is via the constitutive pathway, that transport can be blocked by BFA at two distinct steps and that once on the surface, GLUT2 does not recycle in intracellular vesicles.

Evidence evaluation in fingerprint comparison and automated fingerprint identification systems: modeling between finger variability

In the context of the investigation of the use of automated fingerprint identification systems (AFIS) for the evaluation of fingerprint evidence, the current study presents investigations into the variability of scores from an AFIS system when fingermarks from a known donor are compared to fingerprints that are not from the same source. The ultimate goal is to propose a model, based on likelihood ratios, which allows the evaluation of mark-to-print comparisons. In particular, this model, through its use of AFIS technology, benefits from the possibility of using a large amount of data, as well as from an already built-in proximity measure, the AFIS score. More precisely, the numerator of the LR is obtained from scores issued from comparisons between impressions from the same source and showing the same minutia configuration. The denominator of the LR is obtained by extracting scores from comparisons of the questioned mark with a database of non-matching sources. This paper focuses solely on the assignment of the denominator of the LR. We refer to it by the generic term of between-finger variability. The issues addressed in this paper in relation to between-finger variability are the required sample size, the influence of the finger number and general pattern, as well as that of the number of minutiae included and their configuration on a given finger. Results show that reliable estimation of between-finger variability is feasible with 10,000 scores. These scores should come from the appropriate finger number/general pattern combination as defined by the mark. Furthermore, strategies of obtaining between-finger variability when these elements cannot be conclusively seen on the mark (and its position with respect to other marks for finger number) have been presented. These results immediately allow case-by-case estimation of the between-finger variability in an operational setting.

Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

The cover basement contact beneath the Rawil axial depression (Western Alps) - True amplitude seismic processing, petrophysical properties, and modeling

Since 1986, several near-vertical seismic reflection profiles have been recorded in Switzerland in order to map the deep geologic structure of the Alps. One objective of this endeavour has been to determine the geometries of the autochthonous basement and of the external crystalline massifs, important elements for understanding the geodynamics of the Alpine orogeny. The PNR-20 seismic line W1, located in the Rawil depression of the western Swiss Alps, provides important information on this subject. It extends northward from the `'Penninic front'' across the Helvetic nappes to the Prealps. The crystalline massifs do not outcrop along this profile. Thus, the interpretation of `'near-basement'' reflections has to be constrained by down-dip projections of surface geology, `'true amplitude'' processing, rock physical property studies and modelling. 3-D seismic modelling has been used to evaluate the seismic response of two alternative down-dip projection models. To constrain the interpretation in the southern part of the profile, `'true amplitude'' processing has provided information on the strength of the reflections. Density and velocity measurements on core samples collected up-dip from the region of the seismic line have been used to evaluate reflection coefficients of typical lithologic boundaries in the region. The cover-basement contact itself is not a source of strong reflections, but strong reflections arise from within the overlaying metasedimentary cover sequence, allowing the geometry of the top of the basement to be determined on the basis of `'near-basement'' reflections. The front of the external crystalline massifs is shown to extend beneath the Prealps, about 6 km north of the expected position. A 2-D model whose seismic response shows reflection patterns very similar to the observed is proposed.

The alpha 1a and alpha 1b-adrenergic receptor subtypes: molecular mechanisms of receptor activation and of drug action.

In this chapter we summarize some aspects of the structure-functional relationship of the alpha 1a and alpha 1b-adrenergic receptor subtypes related to the receptor activation process as well as the effect of different alpha-blockers on the constitutive activity of the receptor. Molecular modeling of the alpha 1a and alpha 1b-adrenergic receptor subtypes and computational simulation of receptor dynamics were useful to interpret the experimental findings derived from site directed mutagenesis studies.

Characterisation of "fou2", a constitutive wound-activated mutant and analysis of the proteome and leaf physiology in the region proximal to wounds in Arabidopsis

Résumé : Les jasmonates (JA), une famille d'hor1none végétale, jouent un rôle central dans la réponse à la blessure, et aux attaques d'insectes et de pathogènes. Les JA sont principalement dérivés d'un acide gras, l'acide linolénique. L'addition par une lipoxygénase d'une molécule d'oxygène à l'acide linolénique initie la synthèse de JA. Cependant les mécanismes régulant l'activation de la biosynthèse de JA ne sont pas encore connus. C'est pour cette raison que dans ce travail, nous avons caractérisé chez Arabidopsis thaliana (l'Arabette des Dames) un mutant fou2 dont l'activité lipoxygénase est plus élevée que celle d'une plante sauvage. Les niveaux de JA sont constitutivement plus élevés et l'activation de la synthèse de JA après blessure est fortement plus induite chez fou2 que chez le type sauvage. En outre, fou2 est plus résistant au pathogène Botrytis cinerea et à la chenille Spodoptera littoralis. Afin de comprendre quel mécanisme chez fou2 génére ce phénotype, nous avons cloné le gène responsable du phénotype de fou2. Le mutant fou2 porte une mutation dans le gène d'un canal à deux pores transportant probablement du potassium, du lumen de la vacuole végétale vers le compartiment cytosolique. L'analyse du protéome de fou2 a permis d'identifier une expression plus élevée de sept protéines régulées par les JA ou le stress. La découverte de l'implication d'un canal dans le phénotype de fou2 renforce l'hypothèse que les flux de cations pourraient être impliqués dans les étapes précoces de la synthèse des JA. Nous avons également étudié le protéome et la physiologie d'une feuille blessée, Pour évaluer les changements d'expression protéique en réponse à la blessure et contrôlés par les JA, nous avons quantifié l'expression de 5937 protéines chez une plante d'Arabidopsis sauvage et chez un mutant incapable de synthétiser des JA. Parmi ces 5937 protéines, nous avons identifié 99 protéines régulées par la blessure chez le type sauvage. Nous avons observé pour 65% des protéines dont l'expression protéique changeait après blessure une bonne corrélation entre la quantité de transcrits et de protéines. Plusieurs enzymes de la voie des chorismates impliquées dans la biosynthèse des acides aminés phénoliques étaient induites par les JA après blessure. Une quantification des acides aminés a montré que les niveaux d'acides aminés phénoliques augmentaient significativement après blessure. La blessure induisait aussi des changements dans l'expression de protéines impliquées dans la réponse au stress et particulièrement au stress oxydatif. Nous avons quantifié l'état réduit et oxydé du glutathion, un tripeptide qui, sous sa forme réduite, est l'antioxydant majeur des cellules. Nous avons trouvé une quantité significativement plus élevée de glutathion oxydé chez le type sauvage blessé que chez la plante aus blessée. Ce résultat suggère que la génération d'un stress oxydatif et la proportion relative de glutathions réduits et oxydés sont contrôlés par les JA après blessure. Abstract : Plants possess a family of potent fatty acid-derived wound-response and developmental regulators: the jasmonates. These compounds are derived from the tri?unsaturated fatty acid a-linolenic-acid (18:3). Addition of an oxygen molecule to 18:3 by 13-lipoxygenases (13-LOX) initiates JA biosynthesis. Actually components regulating the activation of JA biosynthesis are poorly defined. Therefore we characterized in Arabidopsis thaliana the fatty acid Qxygenation upregulated 2 (fou2) mutant, which was previously isolated in a screen for mutants with an enhanced 13-LOX activity. As a consequence of this increased 13-LOX activity, JA levels in fou2 are higher than in wild type (WT) and wounding strongly increased JA biosynthesis compared to WT. fou2 was more resistant to the fungus Botrytis cinerea and the generalist caterpillar Spodaptera littomlis, The fou2 mutant carries a missense mutation in the Two Pore Channel 1 gene (TPCJ), which encodes a vacuolar cation channel transporting probably K* into the cytosol. Patchclamp analysis of fou2 vacuolar membranes showed faster time-dependent conductivity and activation of the mutated channel at lower membrane potentials than wild-type. Proteomic analysis of fou2 leaves identified increased levels of seven biotic stress- and JA- inducible proteins. The discovery of the implication of a channel in the fou2 phenotype strenghtens the hypothesis that cation fluxes might be implicated in early steps of JA synthesis. We further concentrated on the proteome and leaf physiology in the region proximal to wounds in Arabidopsis using the WT and the aos JA-biosynthesis deficient mutant in order to find JA- induced proteins changes. We used two successive proteomic methods to assess protein changes in response to wounding Arabidopsis leaves, two dimensional electrophoresis (2DE) and linear trap quadrupole ion-trap mass spectrometry. In total 5937 proteins were quantified. We identified 99 wound-regulated proteins in the WT. Most these proteins were also wound-regulated at the transcript level showing a good correlation between transcript and protein abundance. We identified several wound-regulated enzymes involved in amino acid biosynthesis and confirmed this result by amino acid quantification. Proteins involved in stress reponses were upregulated, particularly in redox species regulation. We found a significantly higher quantity of oxidized glutathione in wounded WT relative to wounded aos leaves. This result suggests that levels of reduced glutathione are controlled by JA after wounding.

Lifelong dynamics of human CD4+CD25+ regulatory T cells: insights from in vivo data and mathematical modeling.

Despite their limited proliferation capacity, regulatory T cells (T(regs)) constitute a population maintained over the entire lifetime of a human organism. The means by which T(regs) sustain a stable pool in vivo are controversial. Using a mathematical model, we address this issue by evaluating several biological scenarios of the origins and the proliferation capacity of two subsets of T(regs): precursor CD4(+)CD25(+)CD45RO(-) and mature CD4(+)CD25(+)CD45RO(+) cells. The lifelong dynamics of T(regs) are described by a set of ordinary differential equations, driven by a stochastic process representing the major immune reactions involving these cells. The model dynamics are validated using data from human donors of different ages. Analysis of the data led to the identification of two properties of the dynamics: (1) the equilibrium in the CD4(+)CD25(+)FoxP3(+)T(regs) population is maintained over both precursor and mature T(regs) pools together, and (2) the ratio between precursor and mature T(regs) is inverted in the early years of adulthood. Then, using the model, we identified three biologically relevant scenarios that have the above properties: (1) the unique source of mature T(regs) is the antigen-driven differentiation of precursors that acquire the mature profile in the periphery and the proliferation of T(regs) is essential for the development and the maintenance of the pool; there exist other sources of mature T(regs), such as (2) a homeostatic density-dependent regulation or (3) thymus- or effector-derived T(regs), and in both cases, antigen-induced proliferation is not necessary for the development of a stable pool of T(regs). This is the first time that a mathematical model built to describe the in vivo dynamics of regulatory T cells is validated using human data. The application of this model provides an invaluable tool in estimating the amount of regulatory T cells as a function of time in the blood of patients that received a solid organ transplant or are suffering from an autoimmune disease.

Protocol for the modeling the epidemiologic transition study: a longitudinal observational study of energy balance and change in body weight, diabetes and cardiovascular disease risk.

ABSTRACT: BACKGROUND: The prevalence of obesity has increased in societies of all socio-cultural backgrounds. To date, guidelines set forward to prevent obesity have universally emphasized optimal levels of physical activity. However there are few empirical data to support the assertion that low levels of energy expenditure in activity is a causal factor in the current obesity epidemic are very limited. METHODS: The Modeling the Epidemiologic Transition Study (METS) is a cohort study designed to assess the association between physical activity levels and relative weight, weight gain and diabetes and cardiovascular disease risk in five population-based samples at different stages of economic development. Twenty-five hundred young adults, ages 25-45, were enrolled in the study; 500 from sites in Ghana, South Africa, Seychelles, Jamaica and the United States. At baseline, physical activity levels were assessed using accelerometry and a questionnaire in all participants and by doubly labeled water in a subsample of 75 per site. We assessed dietary intake using two separate 24-h recalls, body composition using bioelectrical impedance analysis, and health history, social and economic indicators by questionnaire. Blood pressure was measured and blood samples collected for measurement of lipids, glucose, insulin and adipokines. Full examination including physical activity using accelerometry, anthropometric data and fasting glucose will take place at 12 and 24 months. The distribution of the main variables and the associations between physical activity, independent of energy intake, glucose metabolism and anthropometric measures will be assessed using cross-section and longitudinal analysis within and between sites. DISCUSSION: METS will provide insight on the relative contribution of physical activity and diet to excess weight, age-related weight gain and incident glucose impairment in five populations' samples of young adults at different stages of economic development. These data should be useful for the development of empirically-based public health policy aimed at the prevention of obesity and associated chronic diseases.

Influence of structural heterogeneities and of large scale topography on imbricate gravitational rock slope failures: New insights from 3-D physical modeling and geomorphological analysis

Initial topography and inherited structural discontinuities are known to play a dominant role in rock slope stability. Previous 2-D physical modeling results demonstrated that even if few preexisting fractures are activated/propagated during gravitational failure all of those heterogeneities had a great influence on mobilized volume and its kinematics. The question we address in the present study is to determine if such a result is also observed in 3-D. As in 2-D previous models we examine geologically stable model configuration, based upon the well documented landslide at Randa, Switzerland. The 3-D models consisted of a homogeneous material in which several fracture zones were introduced in order to study simplified but realistic configurations of discontinuities (e.g. based on natural example rather than a parametric study). Results showed that the type of gravitational failure (deep-seated landslide or sequential failure) and resulting slope morphology evolution are the result of the interplay of initial topography and inherited preexisting fractures (orientation and density). The three main results are i) the initial topography exerts a strong control on gravitational slope failure. Indeed in each tested configuration (even in the isotropic one without fractures) the model is affected by a rock slide, ii) the number of simulated fracture sets greatly influences the volume mobilized and its kinematics, and iii) the failure zone involved in the 1991 event is smaller than the results produced by the analog modeling. This failure may indicate that the zone mobilized in 1991 is potentially only a part of a larger deep-seated landslide and/or wider deep seated gravitational slope deformation.