Mesoscopic lateral diffusion in lipid bilayers

The lateral diffusion in bilayers is modeled with a multiscale mesoscopic simulation. The methodology consists of two simulations, where the first employs atomistic models to obtain exact results for the mesoscopic model. The second simulation takes the results obtained from the first to parameterize an effective force field that is employed in a new coarse-grained model. The multiscale aspect of this scheme occurs at the point where the microscopic time-averaged results of the first simulation are employed to parameterize the second simulation that operates in a higher spatial and temporal domain. The results of both simulation schemes give quantitative information on the details associated with lipid lateral diffusion.     

A quantitative kinetic model for ATP-induced intracellular Ca2+ oscillations

A quantitative kinetic model is proposed to simulate the ATP-induced intracellular Ca(2+) oscillations. The quantitative effect of ATP concentration upon the oscillations was successfully simulated. Our simulation results support previous experimental explanations that the Ca(2+) oscillations are mainly due to interaction of Ca(2+) release from the endoplasmic reticulum (ER) and the ATP-dependent Ca(2+) pump back into the ER, and the oscillations are prolonged by extracellular Ca(2+) entry that maintains the constant Ca(2+) supplies to its intracellular stores. The model is also able to simulate the sudden disappearance phenomenon of the Ca(2+) oscillations observed in some cell types by taking into account of the biphasic characteristic of the Ca(2+) release from the endoplasmic reticulum (ER). Moreover, the model simulation results for the Ca(2+) oscillations characteristics such as duration, peak [Ca(2+)](cyt), and average interval, etc., lead to prediction of some possible factors responsible for the variations of Ca(2+) oscillations in different types of cells.     

Nucleosome dynamics between tension-induced states

We studied the dynamical behavior of a mononucleosome under tension using a theoretical model that takes into account the nucleosomal geometry, DNA elasticity, nonspecific DNA-protein binding, and effective repulsion between the two DNA turns. Using a dynamical Monte-Carlo simulation algorithm, we demonstrate that this model shows a behavior that for an appropriate set of parameters is in quantitative agreement with data from micromanipulation experiments on individual nucleosomes. All of the parameters of the model follow from the data obtained from two types of pulling experiments, namely, constant force and constant loading rate ensembles.     

Effects of organelle shape on fluorescence recovery after photobleaching

The determination of diffusion coefficients from fluorescence recovery data is often complicated by geometric constraints imposed by the complex shapes of intracellular compartments. To address this issue, diffusion of proteins in the lumen of the endoplasmic reticulum (ER) is studied using cell biological and computational methods. Fluorescence recovery after photobleaching (FRAP) experiments are performed in tissue culture cells expressing GFP-KDEL, a soluble, fluorescent protein, in the ER lumen. The three-dimensional (3D) shape of the ER is determined by confocal microscopy and computationally reconstructed. Within these ER geometries diffusion of solutes is simulated using the method of particle strength exchange. The simulations are compared to experimental FRAP curves of GFP-KDEL in the same ER region. Comparisons of simulations in the 3D ER shapes to simulations in open 3D space show that the constraints imposed by the spatial confinement result in two- to fourfold underestimation of the molecular diffusion constant in the ER if the geometry is not taken into account. Using the same molecular diffusion constant in different simulations, the observed speed of fluorescence recovery varies by a factor of 2.5, depending on the particular ER geometry and the location of the bleached area. Organelle shape considerably influences diffusive transport and must be taken into account when relating experimental photobleaching data to molecular diffusion coefficients. This novel methodology combines experimental FRAP curves with high accuracy computer simulations of diffusion in the same ER geometry to determine the molecular diffusion constant of the solute in the particular ER lumen.     

Rapid, endoplasmic reticulum-independent diffusion of the mitotic Golgi haze

Early in mitosis, the mammalian Golgi apparatus disassembles, and fluorescence microscopy reveals Golgi clusters and an extensive, nonresolvable haze that either represents scattered vesicles or a merged endoplasmic reticulum (ER)-Golgi compartment. To help decide between these alternatives, we have carried out a combined microscopic and pharmacological analysis, by using a BS-C-1 cell line stably coexpressing ER and Golgi markers. Video fluorescence microscopy showed that these two organelles were morphologically distinguishable at all stages of mitosis, and photobleaching experiments showed that diffusion of the Golgi marker was unaffected by the presence of the ER. Fragmentation of the ER by using filipin III completely blocked diffusion of the ER marker but had no effect on the Golgi marker, unless it was first relocated to the ER by using brefeldin A. The Golgi haze was also studied using BODIPY ceramide. Its diffusion was slower in mitotic Golgi than in mitotic ER, but similar to that of a Golgi enzyme marker in the mitotic Golgi haze or in Golgi vesicles generated by ilimaquinone. Together, these results support the idea that the Golgi and the ER remain separate during mitosis and strongly suggest that Golgi markers move by vesicle diffusion, as opposed to lateral diffusion in continuous membranes.     

The making and breaking of the endoplasmic reticulum

The endoplasmic reticulum (ER) is a dynamic organelle central to many essential cellular functions. It is an important calcium store, which functions in cellular signal transduction cascades. It is also the site of entry for secreted proteins into the secretory pathway. Lumenal enzymes will fold and glycosylate these proteins, and if a protein is destined to be secreted, it will be packaged into membrane vesicles that bud off from the ER. The ER is also the site where most cellular lipids are synthesized. It is contiguous with the nuclear envelope, which serves as a diffusion barrier to control entry into and out of the nucleus. In the life cycle of a cell, the ER is in a constant flux of membrane traffic. What maintains the ER in the shape of an intact reticulum among this constant flux of material? We discuss the mechanisms that contribute to the biogenesis of the ER, the maintenance of the organelle, as well as processes that give the ER its characteristic shape and pattern of inheritance.     

Simulations of (an)isotropic diffusion on curved biological surfaces

We present a computational particle method for the simulation of isotropic and anisotropic diffusion on curved biological surfaces that have been reconstructed from image data. The method is capable of handling surfaces of high curvature and complex shape, which are often encountered in biology. The method is validated on simple benchmark problems and is shown to be second-order accurate in space and time and of high parallel efficiency. It is applied to simulations of diffusion on the membrane of endoplasmic reticula (ER) in live cells. Diffusion simulations are conducted on geometries reconstructed from real ER samples and are compared to fluorescence recovery after photobleaching experiments in the same ER samples using the transmembrane protein tsO45-VSV-G, C-terminally tagged with green fluorescent protein. Such comparisons allow derivation of geometry-corrected molecular diffusion constants for membrane components from fluorescence recovery after photobleaching data. The results of the simulations indicate that the diffusion behavior of molecules in the ER membrane differs significantly from the volumetric diffusion of soluble molecules in the lumen of the same ER. The apparent speed of recovery differs by a factor of approximately 4, even when the molecular diffusion constants of the two molecules are identical. In addition, the specific shape of the membrane affects the recovery half-time, which is found to vary by a factor of approximately 2 in different ER samples.     

Why fluorescent probes for endoplasmic reticulum are selective: an experimental and QSAR-modelling study

The basis of the selectivity of fluorochromes routinely used to visualize the endoplasmic reticulum (ER) in live cells remains obscure. To clarify this, interactions of living cells with fluorochromes of varied physicochemical properties were analyzed experimentally and numerically using a quantitative structure activity relationship analysis (QSAR). Routine selective ER probes were found to be amphipathic, lipophilic cations with moderate-sized conjugated systems. The moderately lipophilic character permits probe uptake by passive diffusion without nonspecific accumulation in biomembranes. The moderately amphipathic character favors uptake into the ER, perhaps owing to its high concentration of zwitterionic lipid head-groups. The QSAR model rationalizes the impractical character of some ER probes mentioned in the literature, and could permit design of novel ER probes with different emission colors. The possibility of using the QSAR model as a tool to predict the accumulation of xenobiotics in the ER of living cells is illustrated by the localization of certain antipsychotic drugs in cultured cells.     

A diffusion model for mesoderm induction in amphibian embryos

In this paper we try to answer the question whether diffusion is a possible mechanism to explain mesoderm induction in Amphibians. First the embryological data are discussed and a hypothesis for mesoderm formation is set forth. The blastula being essentially a hollow sphere, we assume that the induction mechanism in an embryo at the blastula stage can be simulated by diffusion-reaction processes on spherical surfaces. A model is constructed for the simple case when the source is held constant with respect to time, the decay proportional to the concentration and the diffusion coefficient a constant, From simulation we find a (best) value for the decay constant to be 6 × 10^–5/sec and for the diffusion constant to be 0.24 × 10^{– 6} cm²/sec. The relation between the parameters is derived from an analytic solution for the diffusion process on a spherical surface with a continuously producing point source and the concentration proportional to the decay. The form and regulative properties of the steady concentration gradient are discussed.     

Why fluorescent probes for endoplasmic reticulum are selective: an experimental and QSAR-modelling study

The basis of the selectivity of fluorochromes routinely used to visualize the endoplasmic reticulum (ER) in live cells remains obscure. To clarify this, interactions of living cells with fluorochromes of varied physicochemical properties were analyzed experimentally and numerically using a quantitative structure activity relationship analysis (QSAR). Routine selective ER probes were found to be amphipathic, lipophilic cations with moderate-sized conjugated systems. The moderately lipophilic character permits probe uptake by passive diffusion without nonspecific accumulation in biomembranes. The moderately amphipathic character favors uptake into the ER, perhaps owing to its high concentration of zwitterionic lipid head-groups. The QSAR model rationalizes the impractical character of some ER probes mentioned in the literature, and could permit design of novel ER probes with different emission colors. The possibility of using the QSAR model as a tool to predict the accumulation of xenobiotics in the ER of living cells is illustrated by the localization of certain antipsychotic drugs in cultured cells.     

A quantitative method for detection of spliced X-box binding protein-1 (XBP1) mRNA as a measure of endoplasmic reticulum (ER) stress

Endoplasmic reticulum (ER) stress is increasingly recognized as an important mechanism in a wide range of diseases including cystic fibrosis, alpha-1 antitrypsin deficiency, Parkinson's and Alzheimer's disease. Therefore, there is an increased need for reliable and quantitative markers for detection of ER stress in human tissues and cells. Accumulation of unfolded or misfolded proteins in the endoplasmic reticulum can cause ER stress, which leads to the activation of the unfolded protein response (UPR). UPR signaling involves splicing of X-box binding protein-1 (XBP1) mRNA, which is frequently used as a marker for ER stress. In most studies, the splicing of the XBP1 mRNA is visualized by gel electrophoresis which is laborious and difficult to quantify. In the present study, we have developed and validated a quantitative real-time RT-PCR method to detect the spliced form of XBP1 mRNA.     

Photoactivation of GFP reveals protein dynamics within the endoplasmic reticulum membrane

Components of the plant cell secretory pathway, including the endoplasmic reticulum and Golgi apparatus, are in constant motion. The photoactivation of GFP has been used to determine that proteins within the membrane of the ER flow as the ER is remodelled. Measurement of the rate at which activated GFP moves away from the activation spot shows that this motion is much faster than would be expected if membrane components moved simply by diffusion. Treatment with latrunculin to depolymerize the actin cytoskeleton stops ER remodelling and reduces the rate of GFP movement to that expected from diffusion alone. This suggests that myosin binds directly or indirectly to ER membrane proteins and actively moves them around over the actin scaffold. Tracking of Golgi body movement was used to demonstrate that they move at the same rate and in the same direction as do photoactivated ER surface proteins. Golgi bodies, therefore, move with, and not over, the surface of the ER. These observations support the current theory of continuity between Golgi bodies and discrete ER exit sites in the ER membrane.     

Estrogen receptor binding to a DNA response element in vitro is not dependent upon estradiol

Gel shift assays were employed to distinguish between the contribution of 17 beta-estradiol (E2) and a short heating step to the ability of the rat uterine cytosolic estrogen receptor (ER) to bind to the estrogen response element (ERE) from the vitellogenin A2 gene (vitERE). Despite the popularity of models in which the ER is a ligand-activated DNA-binding protein, these studies find that estrogen does not significantly contribute to receptor-DNA complex formation. An avidin-biotin complex with DNA (ABCD) assay was utilized to obtain quantitative measurement of the affinities of the ER for the vitERE and a mutant sequence. Scatchard analysis gave a dissociation constant of 390 +/- 40 pM for the E2-occupied, heated ER to the vitERE. The data fit a one-site model and evidence for cooperatively was not observed. A dissociation constant of 450 +/- 170 pM was obtained for the unoccupied, heated ER, leading to the conclusion that estrogen was not necessary for specific binding to DNA. The percentage of ER capable of binding vitERE varied with each cytosol preparation, ranging from 60 to 100% and estrogen did not appear to affect this variation. Competition against the vitERE with a 2-bp mutant sequence showed a 250-fold lower relative binding affinity of the receptor for the mutant over the vitERE sequence. This ability of the ER to discriminate between target and nonspecific DNA sequences was also not dependent on the presence of estrogen.     

Active translocon complexes labeled with GFP-Dad1 diffuse slowly as large polysome arrays in the endoplasmic reticulum

In the ER, the translocon complex (TC) functions in the translocation and cotranslational modification of proteins made on membrane-bound ribosomes. The oligosaccharyltransferase (OST) complex is associated with the TC, and performs the cotranslational N-glycosylation of nascent polypeptide chains. Here we use a GFP-tagged subunit of the OST complex (GFP-Dad1) that rescues the temperature-sensitive (ts) phenotype of tsBN7 cells, where Dad1 is degraded and N-glycosylation is inhibited, to study the lateral mobility of the TC by FRAP. GFP-Dad1 that is functionally incorporated into TCs diffuses extremely slow, exhibiting an effective diffusion constant (Deff) about seven times lower than that of GFP-tagged ER membrane proteins unhindered in their lateral mobility. Termination of protein synthesis significantly increases the lateral mobility of GFP-Dad1 in the ER membranes, but to a level that is still lower than that of free GFP-Dad1. This suggests that GFP-Dad1 as part of the OST remains associated with inactive TCs. Our findings that TCs assembled into membrane-bound polysomes diffuse slowly within the ER have mechanistic implications for the segregation of the ER into smooth and rough domains.     

Dynamics of the Ins(1,4,5)P3 receptor during polarization of MDCK cells

BACKGROUND INFORMATION: The uneven distribution of the Ins(1,4,5)P3R [Ins(1,4,5)P3 receptor] within the ER (endoplasmic reticulum) membrane generates spatially complex Ca2+ signals. The ER is a dynamic network, which allows the rapid diffusion of membrane proteins from one part of the cell to another. However, little is known about the localization and the dynamics of the Ins(1,4,5)P3R in the ER of living cells. We have used a MDCK (Madin-Darby canine kidney) clone stably expressing the Ins(1,4,5)P3R1-GFP (where GFP stands for green fluorescent protein) to investigate the effect of cell polarity on the lateral mobility of the Ins(1,4,5)P3R. RESULTS: In non-confluent MDCK cells, the chimaera is homogeneously distributed throughout the ER and the nuclear envelope. FRAP (fluorescence recovery after photobleaching) experiments showed that the receptor can move freely in the ER with a diffusion constant (D=0.01 microm2/s) approx. ten times lower than other ER membrane proteins. In confluent polarized cells, two populations of receptor can be defined: one population is distributed in the cytoplasm and is mobile but with a slower diffusion constant (D=0.004 microm2/s) compared with non-confluent cells, whereas the other population is concentrated at the periphery of the cells and is apparently immobile. CONCLUSIONS: The observed differences in the mobility of the Ins(1,4,5)P3R are most probably due to its interactions with stable protein complexes that form at the periphery of the polarized cells.     

Anisotropic nutrient transport in three-dimensional single species bacterial biofilms

The ability for a biofilm to grow and function is critically dependent on the nutrient availability, and this in turn is dependent on the structure of the biofilm. This relationship is therefore an important factor influencing biofilm maturation. Nutrient transport in bacterial biofilms is complex; however, mathematical models that describe the transport of particles within biofilms have made three simplifying assumptions: the effective diffusion coefficient (EDC) is constant, the EDC is that of water, and/or the EDC is isotropic. Using a Monte Carlo simulation, we determined the EDC, both parallel to and perpendicular to the substratum, within 131 real, single species, three-dimensional biofilms that were constructed from confocal laser scanning microscopy images. Our study showed that diffusion within bacterial biofilms was anisotropic and depth dependent. The heterogeneous distribution of bacteria varied between and within species, reducing the rate of diffusion of particles via steric hindrance. In biofilms with low porosity, the EDCs for nutrient transport perpendicular to the substratum were significantly lower than the EDCs for nutrient transport parallel to the substratum. Here, we propose a reaction-diffusion model to describe the nutrient concentration within a bacterial biofilm that accounts for the depth dependence of the EDC.     

Ratcheting in post-translational protein translocation: a mathematical model

We have developed a non-steady-state mathematical model describing post-translational protein translocation across the endoplasmic reticulum membrane. Movement of the polypeptide chain through the channel in the endoplasmic reticulum membrane is considered to be a stochastic process which is biased at the lumenal side of the channel by the binding of BiP (Kar2p), a member of the Hsp70 family of ATPases (ratcheting model). Assuming that movement of the chain through the channel is caused by passive diffusion (Brownian ratchet), the model describes all available experimental data. The optimum set of model parameters indicates that the ratcheting mechanism functions at near-maximum rate, being relatively insensitive to variations of the association or dissociation rate constants of BiP or its concentration. The estimated rate constant for diffusion of a polypeptide inside the channel indicates that the chain makes contact with the walls of the channel. Since fitting of the model to the data required that the backward rate constant be larger than the forward constant during early diffusion steps, translocation must occur against a force. The latter may arise, for example, from the unfolding of the polypeptide chain in the cytosol. Our results indicate that the ratchet can transport polypeptides against a free energy of about 25 kJ/mol without significant retardation of translocation. The modeling also suggests that the BiP ratchet is optimized, allowing fast translocation to be coupled with minimum consumption of ATP and rapid dissociation of BiP in the lumen of the ER. Finally, we have estimated the maximum hydrophobicity of a polypeptide segment up to which lateral partitioning from the channel into the lipid phase does not result in significant retardation of translocation.     

Maximum likelihood techniques for the mapping and analysis of quantitative trait loci with the aid of genetic markers

A method is presented to estimate the biometric parameters of a quantitative trait locus linked to a genetic marker when both loci are segregating in the F-2 generation of a cross between two inbred lines. The method, which assumes underlying normal distributions, is a combination of maximum likelihood and moments methods and uses the statistics of the genetic marker genotype samples for the quantitative trait to estimate the recombination frequency between the two loci and the means and variances of the genotypes of the quantitative trait locus. With this method, the genetic parameters of a locus affecting plant height linked to an electrophoretic marker for esterase were accurately estimated from a sample of 1596 F-2 progeny of a cross between two species of Lycopersicon (tomato). Linkage distance between the two loci was 38 map units and the effect of the quantitative trait locus was 1.6 phenotypic standard deviation units. Accurate estimates of the genetic parameters and linkage distance for populations of 2000 individuals simulated with a segregating codominant locus with an effect of 1.63 standard deviations linked to a genetic marker with .2 recombination were also derived by this method. The method is not effective in distinguishing between complete and partial linkage in samples of only 500 individuals or for quantitative loci with effects less than a phenotypic standard deviation. The method is more effective for codominant than for dominant loci.     

Calcium regulates the association between mitochondria and a smooth subdomain of the endoplasmic reticulum

Association between the ER and mitochondria has long been observed, and the formation of close contacts between ER and mitochondria is necessary for the ER-mediated sequestration of cytosolic calcium by mitochondria. Autocrine motility factor receptor (AMF-R) is a marker for a smooth subdomain of the ER, shown here by confocal microscopy to be distinct from, yet closely associated with the calnexin- or calreticulin-labeled ER. By EM, smooth ER AMF-R tubules exhibit direct interactions with mitochondria, identifying them as a mitochondria-associated smooth ER subdomain. In digitonin-permeabilized MDCK cells, the addition of rat liver cytosol stimulates the dissociation of smooth ER and mitochondria under conditions of low calcium. Using BAPTA chelators of various affinities and CaEGTA buffers of defined free Ca(2+) concentrations and quantitative confocal microscopy, we show that free calcium concentrations <100 nM favor dissociation, whereas those >1 microM favor close association between these two organelles. Therefore, we describe a cellular mechanism that facilitates the close association of this smooth ER subdomain and mitochondria when cytosolic free calcium rises above physiological levels.