A partial velocity approach to subcycling structural dynamics

Subcycling, or the use of different timesteps at different nodes, can be an effective way of improving the computational efficiency of explicit transient dynamic structural solutions. The method that has been most widely adopted uses a nodal partition. extending the central difference method, in which small timestep updates are performed interpolating on the displacement at neighbouring large timestep nodes. This approach leads to narrow bands of unstable timesteps or statistical stability. It also can be in error due to lack of momentum conservation on the timestep interface. The author has previously proposed energy conserving algorithms that avoid the first problem of statistical stability. However, these sacrifice accuracy to achieve stability. An approach to conserve momentum on an element interface by adding partial velocities is considered here. Applied to extend the central difference method. this approach is simple. and has accuracy advantages. The method can be programmed by summing impulses of internal forces, evaluated using local element timesteps, in order to predict a velocity change at a node. However, it is still only statistically stable, so an adaptive timestep size is needed to monitor accuracy and to be adjusted if necessary. By replacing the central difference method with the explicit generalized alpha method. it is possible to gain stability by dissipating the high frequency response that leads to stability problems. However. coding the algorithm is less elegant, as the response depends on previous partial accelerations. Extension to implicit integration, is shown to be impractical due to the neglect of remote effects of internal forces acting across a timestep interface. (C) 2002 Elsevier Science B.V. All rights reserved.

The tree cut and merge algorithm for estimation of network reliability

This article presents Monte Carlo techniques for estimating network reliability. For highly reliable networks, techniques based on graph evolution models provide very good performance. However, they are known to have significant simulation cost. An existing hybrid scheme (based on partitioning the time space) is available to speed up the simulations; however, there are difficulties with optimizing the important parameter associated with this scheme. To overcome these difficulties, a new hybrid scheme (based on partitioning the edge set) is proposed in this article. The proposed scheme shows orders of magnitude improvement of performance over the existing techniques in certain classes of network. It also provides reliability bounds with little overhead.

GLUT4 recycles via a trans-Golgi network (TGN) subdomain enriched in Syntaxins 6 and 16 but not TGN38: Involvement of an acidic targeting motif

Insulin stimulates glucose transport in fat and muscle cells by triggering exocytosis of the glucose transporter GLUT4. To define the intracellular trafficking of GLUT4, we have studied the internalization of an epitope-tagged version of GLUT4 from the cell surface. GLUT4 rapidly traversed the endosomal system en route to a perinuclear location. This perinuclear GLUT4 compartment did not colocalize with endosomal markers (endosomal antigen I protein, transferrin) or TGN38, but showed significant overlap with the TGN target (t)-soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) Syntaxins 6 and 16. These results were confirmed by vesicle immunoisolation. Consistent with a role for Syntaxins 6 and 16 in GLUT4 trafficking we found that their expression was up-regulated significantly during adipocyte differentiation and insulin stimulated their movement to the cell surface. GLUT4 trafficking between endosomes and trans-Golgi network was regulated via an acidic targeting motif in the carboxy terminus of GLUT4, because a mutant lacking this motif was retained in endosomes. We conclude that GLUT4 is rapidly transported from the cell surface to a subdomain of the trans-Golgi network that is enriched in the t-SNAREs Syntaxins 6 and 16 and that an acidic targeting motif in the C-terminal tail of GLUT4 plays an important role in this process.

Dynamics of Gray Matter Loss in Alzheimer's Disease

We detected and mapped a dynamically spreading wave of gray matter loss in the brains of patients with Alzheimer's disease (AD). The loss pattern was visualized in four dimensions as it spread over time from temporal and limbic cortices into frontal and occipital brain regions, sparing sensorimotor cortices. The shifting deficits were asymmetric (left hemisphere >right hemisphere) and correlated with progressively declining cognitive status ( p 15% loss). The maps distinguished different phases of AD and differentiated AD from normal aging. Local gray matter loss rates (5.3 +/- 2.3% per year in AD v 0.9 +/- 0.9% per year in controls) were faster in the left hemisphere ( p < 0.029) than the right. Transient barriers to disease progression appeared at limbic/frontal boundaries. This degenerative sequence, observed in vivo as it developed, provides the first quantitative, dynamic visualization of cortical atrophic rates in normal elderly populations and in those with dementia.

The trans-membrane protein p25 forms highly specialized domains that regulate membrane composition and dynamics

Trans-membrane proteins of the p24 family are abundant, oligomeric proteins predominantly found in cis-Golgi membranes. They are not easily studied in vivo and their functions are controversial. We found that p25 can be targeted to the plasma membrane after inactivation of its canonical KKXX motif (KK to SS, p25SS), and that p25SS causes the co-transport of other p24 proteins beyond the Golgi complex, indicating that wild-type p25 plays a crucial role in retaining p24 proteins in cis-Golgi membranes. We then made use of these observations to study the intrinsic properties of these proteins, when present in a different membrane context. At the cell surface, the p25SS mutant segregates away from both the transferrin receptor and markers of lipid rafts, which are enriched in cholesterol and glycosphingolipids. This suggests that p25SS localizes to, or contributes to form, specialized membrane domains, presumably corresponding to oligomers of p25SS and other p24 proteins. Once at the cell surface, p25SS is endocytosed, together with other p24 proteins, and eventually accumulates in late endosomes, where it remains confined to well-defined membrane regions visible by electron microscopy. We find that this p25SS accumulation causes a concomitant accumulation of cholesterol in late endosomes, and an inhibition of their motility - two processes that are functionally linked. Yet, the p25SS-rich regions themselves seem to-exclude not only Lamp1 but also accumulated cholesterol. One may envision that p25SS accumulation, by excluding cholesterol from oligomers, eventually overloads neighboring late endosomal membranes with cholesterol beyond their capacity (see Discussion). In any case, our data show that p25 and presumably other p24 proteins are endowed with the intrinsic capacity to form highly specialized domains that control membrane composition and dynamics. We propose that p25 and other p24 proteins control the fidelity of membrane transport by maintaining cholesterol-poor membranes in the Golgi complex.

GRIP domain-mediated targeting of two new coiled-coil proteins, GCC88 and GCC185, to subcompartments of the trans-Golgi network

The GRIP domain is a targeting sequence found in a family of coiled-coil peripheral Golgi proteins. Previously we demonstrated that the GRIP domain of p230/golgin245 is specifically recruited to tubulovesicular structures of the traps-Golgi network (TGN). Here we have characterized two novel Golgi proteins with functional GRIP domains, designated GCC88 and GCC185. GCC88 cDNA encodes a protein of 88 kDa, and GCC185 cDNA encodes a protein of 185 kDa. Both molecules are brefeldin A-sensitive peripheral membrane proteins and are predicted to have extensive coiled-coil regions with the GRIP domain at the C terminus. By immunofluorescence and immunoelectron microscopy GCC88 and GCC185, and the GRIP protein golgin97, are all localized to the TGN of Hela cells. Overexpression of full-length GCC88 leads to the formation of large electron dense structures that extend from the traps-Golgi. These de novo structures contain GCC88 and co-stain for the TGN markers syntaxin 6 and TGN38 but not for alpha2,6-sialyltransferase, beta-COP, or cis-Golgi GM130. The formation of these abnormal structures requires the N-terminal domain of GCC88. TGN38, which recycles between the TGN and plasma membrane, was transported into and out of the GCC88 decorated structures. These data introduce two new GRIP domain proteins and implicate a role for GCC88 in the organization of a specific TGN subcompartment involved with membrane transport.

GAIP participates in budding of membrane carriers at the trans-Golgi network

Galpha interacting protein (GAIP) is a regulator of G protein signaling protein that associates dynamically with vesicles and has been implicated in membrane trafficking, although its specific role is not yet known. Using an in vitro budding assay, we show that GAIP is recruited to a specific population of trans-Golgi network-derived vesicles and that these are distinct from coatomer or clathrin-coated vesicles. A truncation mutant (NT-GAIP) encoding only the N-terminal half of GAIP is recruited to trans -Golgi network membranes during the formation of vesicle carriers. Overexpression of NT-GAIP induces the formation of long, coated tubules, which are stabilized by microtubules. Results from the budding assay and from imaging in live cells show that these tubules remain attached to the Golgi stack rather than being released as carrier vesicles. NT-GAIP expression blocks membrane budding and results in the accumulation of tubular carrier intermediates. NT-GAIP-decorated tubules are competent to load vesicular stomatitis virus protein G-green fluorescent protein as post-Golgi, exocytic cargo and in cells expressing NT-GAIP there is reduced surface delivery of vesicular stomatitis virus protein G-green fluorescent protein. We conclude that GAIP functions as an essential part of the membrane budding machinery for a subset of post-Golgi exocytic carriers derived from the trans-Golgi network.

Influence of hysteresis on tidal capillary fringe dynamics in a well-sorted sand

Nielsen and Perrochet [Adv. Water Resour. 23 (2000) 503] presented experimental data for cyclic water movement in the vadose zone above an oscillating watertable. The response of the watertable to cyclic forcing was characterised by the ratios of the forcing head to watertable amplitudes and their associated phase lag. They found that their non-hysteretic Richards' equation model failed to represent the observed behaviour of these parameters. This paper explores the effect on the simulated capillary fringe dynamics (in terms of these parameters) of including varying degrees of hysteresis in the moisture retention curve used in a numerical model of their experiment. It is clear that hysteresis can indeed account for observed discrepancies between simulation and experiment and that the effect of hysteresis varies with the frequency of oscillation. The use of a single-valued mean retention curve, as advocated by some authors, fails to provide a match between the simulated and observed behaviour of the Nielsen and Perrochet parameters, but is shown to be adequate for predicting time-averaged soil moisture profiles. (C) 2003 Elsevier Ltd. All rights reserved.

A reduced load approximation accounting for link interactions in a loss network

This paper is concerned with evaluating the performance of loss networks. Accurate determination of loss network performance can assist in the design and dimen- sioning of telecommunications networks. However, exact determination can be difficult and generally cannot be done in reasonable time. For these reasons there is much interest in developing fast and accurate approximations. We develop a reduced load approximation that improves on the famous Erlang fixed point approximation (EFPA) in a variety of circumstances. We illustrate our results with reference to a range of networks for which the EFPA may be expected to perform badly.

Genotype and environment effects on dynamics of harvest index during grain filling in sorghum

An approach based on a linear rate of increase in harvest index (141) with time after anthesis has been used as a simple means-to predict grain growth and yield in many crop simulation models. When applied to diverse situations, however, this approach has been found to introduce significant error in grain yield predictions. Accordingly, this study was undertaken to examine the stability of the HI approach for yield prediction in sorghum [Sorghum bicolor (L.) Moench]. Four field experiments were conducted under nonlimiting water. and N conditions. The experiments were sown at times that ensured a broad range in temperature and radiation conditions. Treatments consisted of two population densities and three genotypes varying in maturity. Frequent sequential harvests were used to monitor crop growth, yield, and the dynamics of 111. Experiments varied greatly in yield and final HI. There was also a tendency for lower HI with later maturity. Harvest index dynamics also varied among experiments and, to a lesser extent, among treatments within experiments. The variation was associated mostly with the linear rate of increase in HI and timing of cessation of that increase. The average rate of HI increase was 0.0198 d(-1), but this was reduced considerably (0.0147) in one experiment that matured in cool conditions. The variations found in IN dynamics could be largely explained by differences in assimilation during grain filling and remobilization of preanthesis assimilate. We concluded that this level of variation in HI dynamics limited the general applicability of the HI approach in yield prediction and suggested a potential alternative for testing.

SPAnDX: a process-based population dynamics model to explore management and climate change impacts on an invasive alien plant, Acacia nilotica

This paper describes a process-based metapopulation dynamics and phenology model of prickly acacia, Acacia nilotica, an invasive alien species in Australia. The model, SPAnDX, describes the interactions between riparian and upland sub-populations of A. nilotica within livestock paddocks, including the effects of extrinsic factors such as temperature, soil moisture availability and atmospheric concentrations of carbon dioxide. The model includes the effects of management events such as changing the livestock species or stocking rate, applying fire, and herbicide application. The predicted population behaviour of A. nilotica was sensitive to climate. Using 35 years daily weather datasets for five representative sites spanning the range of conditions that A. nilotica is found in Australia, the model predicted biomass levels that closely accord with expected values at each site. SPAnDX can be used as a decision-support tool in integrated weed management, and to explore the sensitivity of cultural management practices to climate change throughout the range of A. nilotica. The cohort-based DYMEX modelling package used to build and run SPAnDX provided several advantages over more traditional population modelling approaches (e.g. an appropriate specific formalism (discrete time, cohort-based, process-oriented), user-friendly graphical environment, extensible library of reusable components, and useful and flexible input/output support framework). (C) 2003 Published by Elsevier Science B.V.

Effect of pore-network connectivity on multicomponent adsorption of large molecules

The effect of pore-network connectivity on binary liquid-phase adsorption equilibria using the ideal adsorbed solution theory (LAST) was studied. The liquid-phase binary adsorption experiments used ethyl propionate, ethyl butyrate, and ethyl isovalerate as the adsorbates and commercial activated carbons Filtrasorb-400 and Norit ROW 0.8 as adsorbents. As the single-component isotherm, a modified Dubinin-Radushkevich equation was used. A comparison with experimental data shows that incorporating the connectivity of the pore network and considering percolation processes associated with different molecular sizes of the adsorptives in the mixture, as well as their different corresponding accessibility, can improve the prediction of binary adsorption equilibria using the LAST Selectivity of adsorption for the larger molecule in binary systems increases with an increase in the pore-network coordination number, as well with an increase in the mean pore width and in the spread of the pore-size distribution.