The role of transcription factories-mediated interchromosomal contacts in the organization of nuclear architecture.

Using numerical simulations, we investigate the underlying physical effects responsible for the overall organization of chromosomal territories in interphase nuclei. In particular, we address the following three questions: (i) why are chromosomal territories with relatively high transcriptional activity on average, closer to the centre of cell's nucleus than those with the lower activity? (ii) Why are actively transcribed genes usually located at the periphery of their chromosomal territories? (iii) Why are pair-wise contacts between active and inactive genes less frequent than those involving only active or only inactive genes? We show that transcription factories-mediated contacts between active genes belonging to different chromosomal territories are instrumental for all these features of nuclear organization to emerge spontaneously due to entropic effects arising when chromatin fibres are highly crowded.

Effects of supercoiling on enhancer-promoter contacts.

Using Brownian dynamics simulations, we investigate here one of possible roles of supercoiling within topological domains constituting interphase chromosomes of higher eukaryotes. We analysed how supercoiling affects the interaction between enhancers and promoters that are located in the same or in neighbouring topological domains. We show here that enhancer-promoter affinity and supercoiling act synergistically in increasing the fraction of time during which enhancer and promoter stay in contact. This stabilizing effect of supercoiling only acts on enhancers and promoters located in the same topological domain. We propose that the primary role of recently observed supercoiling of topological domains in interphase chromosomes of higher eukaryotes is to assure that enhancers contact almost exclusively their cognate promoters located in the same topological domain and avoid contacts with very similar promoters but located in neighbouring topological domains.

Modeling complex wells with the multi-scale finite-volume method

In this paper, an extension of the multi-scale finite-volume (MSFV) method is devised, which allows to Simulate flow and transport in reservoirs with complex well configurations. The new framework fits nicely into the data Structure of the original MSFV method,and has the important property that large patches covering the whole well are not required. For each well. an additional degree of freedom is introduced. While the treatment of pressure-constraint wells is trivial (the well-bore reference pressure is explicitly specified), additional equations have to be solved to obtain the unknown well-bore pressure of rate-constraint wells. Numerical Simulations of test cases with multiple complex wells demonstrate the ability of the new algorithm to capture the interference between the various wells and the reservoir accurately. (c) 2008 Elsevier Inc. All rights reserved.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Fas (CD95/APO-1) limits the expansion of T lymphocytes in an environment of limited T-cell antigen receptor/MHC contacts.

Fas-deficient mice (Fas(lpr/lpr)) and humans have profoundly dysregulated T lymphocyte homeostasis, which manifests as an accumulation of CD4(+) and CD8(+) T cells as well as an unusual population of CD4(-)CD8(-)TCRαβ(+) T cells. To date, no unifying model has explained both the increased T-cell numbers and the origin of the CD4(-)CD8(-)TCRαβ(+) T cells. As Fas(lpr/lpr) mice raised in a germ-free environment still manifest lymphadenopathy, we considered that this process is primarily driven by recurrent low-avidity TCR signaling in response to self-peptide/MHC as occurs during homeostatic proliferation. In these studies, we developed two independent systems to decrease the number of self-peptide/MHC contacts. First, expression of MHC class I was reduced in OT-I TCR transgenic mice. Although OT-I Fas(lpr/lpr) mice did not develop lymphadenopathy characteristic of Fas(lpr/lpr) mice, in the absence of MHC class I, OT-I Fas(lpr/lpr) T cells accumulated as both CD8(+) and CD4(-)CD8(-) T cells. In the second system, re-expression of β(2)m limited to thymic cortical epithelial cells of Fas(lpr/lpr) β(2)m-deficient mice yielded a model in which polyclonal CD8(+) thymocytes entered a peripheral environment devoid of MHC class I. These mice accumulated significantly greater numbers of CD4(-)CD8(-)TCRαβ(+) T cells than conventional Fas(lpr/lpr) mice. Thus, Fas shapes the peripheral T-cell repertoire by regulating the survival of a subset of T cells proliferating in response to limited self-peptide/MHC contacts.

A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

Importance of polyethylene thickness in total shoulder arthroplasty: a finite element analysis.

BACKGROUND: Articular surfaces reconstruction is essential in total shoulder arthroplasty. Because of the limited glenoid bone support, thin glenoid component could improve anatomical reconstruction, but adverse mechanical effects might appear. METHODS: With a numerical musculoskeletal shoulder model, we analysed and compared three values of thickness of a typical all-polyethylene glenoid component: 2, 4 (reference) and 6mm. A loaded movement of abduction in the scapular plane was simulated. We evaluated the humeral head translation, the muscle moment arms, the joint force, the articular contact pattern, and the polyethylene and cement stress. Findings Decreasing polyethylene thickness from 6 to 2mm slightly increased humeral head translation and muscle moment arms. This induced a small decreased of the joint reaction force, but important increase of stress within the polyethylene and the cement mantel. Interpretation The reference thickness of 4mm seems a good compromise to avoid stress concentration and joint stuffing.

An operator formulation of the multiscale finite-volume method with correction function

The multiscale finite-volume (MSFV) method has been derived to efficiently solve large problems with spatially varying coefficients. The fine-scale problem is subdivided into local problems that can be solved separately and are coupled by a global problem. This algorithm, in consequence, shares some characteristics with two-level domain decomposition (DD) methods. However, the MSFV algorithm is different in that it incorporates a flux reconstruction step, which delivers a fine-scale mass conservative flux field without the need for iterating. This is achieved by the use of two overlapping coarse grids. The recently introduced correction function allows for a consistent handling of source terms, which makes the MSFV method a flexible algorithm that is applicable to a wide spectrum of problems. It is demonstrated that the MSFV operator, used to compute an approximate pressure solution, can be equivalently constructed by writing the Schur complement with a tangential approximation of a single-cell overlapping grid and incorporation of appropriate coarse-scale mass-balance equations.