Search

Extrapolating in vivo hepatic clearance from in vitro uptake data is a known challenge, especially for organic anion-transporting polypeptide transporter (OATP) substrates, and the well-stirred model (WSM) commonly yields systematic underpredictions for those anionic drugs, hypothetically due to "albumin-mediated hepatic drug uptake". In the present study, we demonstrate that the WSM incorporating the dynamic free fraction ( f D ), a measure of drug protein binding affinity, performs reasonably well in predicting hepatic clearance of OATP substrates. For a selection of anionic drugs, including atorvastatin, fluvastatin, pravastatin, rosuvastatin, pitavastatin, cerivastatin, and repaglinide, this dynamic well-stirred model (dWSM) correctly predicts hepatic plasma clearance within 2-fold error for six out of seven OATP substrates examined. The geometric mean of clearance ratios between the predicted and the observed values falls in the range of 1.21-1.38. As expected, the WSM with unbound fraction ( f u ) systematically underpredicts hepatic clearance with greater than 2-fold error for five out of seven drugs, and the geometric mean of clearance ratios between the predicted and the observed values is in the range of 0.20-0.29. The results suggest that, despite its simplicity, the dWSM operates well for transporter-mediated uptake clearance, and that clearance under-prediction of OATP substrates may not necessarily be associated with the chemical class of the anionic drugs, nor is it a result of albumin-mediated hepatic drug uptake as currently hypothesized. Instead, the superior prediction power of the dWSM confirms the utility of the dynamic free fraction in clearance prediction and the importance of drug plasma binding kinetics in hepatic uptake clearance. SIGNIFICANCE STATEMENT: The traditional well-stirred model (WSM) consistently underpredicts organin anion-transporting polypeptide transporter (OATP)-mediated hepatic uptake clearance, hypothetically due to the albumin-mediated hepatic drug uptake. In this manuscript, we apply the dynamic WSM to extrapolate hepatic clearance of the OATP substrates, and our results show significant improvements in clearance prediction without assuming albumin-mediated hepatic drug uptake., (Copyright © 2024 by The American Society for Pharmacology and Experimental Therapeutics.)

Most nutrient uptake problems are modeled by the convection-diffusion equation (CDE) abiding by Fick's law. Because nutrients needed by plants exist in the soil solution as a form of ions and the soil is a typical fractal structure of heterogeneity, it makes the solute transport appear anomalous diffusion in soil. Taking anomalous diffusion as a transport process, we propose time and space fractional nutrient uptake models based on the classic Nye-Tinker-Barber model. There does not appear apparent sub-diffusion of nitrate in the time fractional model until four months and the time fractional models are appropriate for describing long-term dynamics and slow sorption reaction; the space fractional model can capture super-diffusion in short term and it is suitable for describing nonlocal phenomena and daily variations driven by transpiration and metabolism; the anomalous diffusion more apparently appears near the root surface in the modeling simulation., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Homeostasis in living cells refers to the steady state of internal, physical, and chemical conditions. It is sustained by self-regulation of the dynamic cellular system. To gain insight into the homeostatic mechanisms that maintain cytosolic nutrient concentrations in plant cells within a homeostatic range, we performed computational cell biology experiments. We mathematically modeled membrane transporter systems and simulated their dynamics. Detailed analyses of 'what-if' scenarios demonstrated that a single transporter type for a nutrient, irrespective of whether it is a channel or a cotransporter, is not sufficient to calibrate a desired cytosolic concentration. A cell cannot flexibly react to different external conditions. Rather, at least two different transporter types for the same nutrient, which are energized differently, are required. The gain of flexibility in adjusting a cytosolic concentration was accompanied by the establishment of energy-consuming cycles at the membrane, suggesting that these putatively "futile" cycles are not as futile as they appear. Accounting for the complex interplay of transporter networks at the cellular level may help design strategies for increasing nutrient use efficiency of crop plants., (© American Society of Plant Biologists 2021. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

The acidification of plant vacuoles is of great importance for various physiological processes, as a multitude of secondary active transporters utilize the proton gradient established across the vacuolar membrane. Vacuolar-type H + -translocating ATPases and a pyrophosphatase are thought to enable vacuoles to accumulate protons against their electrochemical potential. However, recent studies pointed to the ATPase located at the trans-Golgi network/early endosome (TGN/EE) to contribute to vacuolar acidification in a manner not understood as of now. Here, we combined experimental data and computational modeling to test different hypotheses for vacuolar acidification mechanisms. For this, we analyzed different models with respect to their ability to describe existing experimental data. To better differentiate between alternative acidification mechanisms, new experimental data have been generated. By fitting the models to the experimental data, we were able to prioritize the hypothesis in which vesicular trafficking of Ca 2+ /H + -antiporters from the TGN/EE to the vacuolar membrane and the activity of ATP-dependent Ca 2+ -pumps at the tonoplast might explain the residual acidification observed in Arabidopsis mutants defective in vacuolar proton pump activity. The presented modeling approach provides an integrative perspective on vacuolar pH regulation in Arabidopsis and holds potential to guide further experimental work., (© 2021 Society for Experimental Biology and John Wiley & Sons Ltd.)

This review is a summary of factors affecting the drug pharmacokinetics (PK) of dogs versus humans. Identifying these interspecies differences can facilitate canine-human PK extrapolations while providing mechanistic insights into species-specific drug in vivo behavior. Such a cross-cutting perspective can be particularly useful when developing therapeutics targeting diseases shared between the two species such as cancer, diabetes, cognitive dysfunction, and inflammatory bowel disease. Furthermore, recognizing these differences also supports a reverse PK extrapolations from humans to dogs. To appreciate the canine-human differences that can affect drug absorption, distribution, metabolism, and elimination, this review provides a comparison of the physiology, drug transporter/enzyme location, abundance, activity, and specificity between dogs and humans. Supplemental material provides an in-depth discussion of certain topics, offering additional critical points to consider. Based upon an assessment of available state-of-the-art information, data gaps were identified. The hope is that this manuscript will encourage the research needed to support an understanding of similarities and differences in human versus canine drug PK.

We study the dynamics of a model of membrane vesicle transport into dendritic spines, which are bulbous intracellular compartments in neurons driven by molecular motors. We reduce the lubrication model proposed in Fai et al. (Phys Rev Fluids 2:113601, 2017) to a fast-slow system, yielding an analytically and numerically tractable equation equivalent to the original model in the overdamped limit. The model's key parameters include: (1) the ratio of motors that prefer to push toward the head of the dendritic spine to the motors that prefer to push in the opposite direction, and (2) the viscous drag exerted on the vesicle by the spine constriction. We perform a numerical bifurcation analysis in these parameters and find that steady-state vesicle velocities appear and disappear through several saddle-node bifurcations. This process allows us to identify the region of parameter space in which multiple stable velocities exist. We show by direct calculations that there can only be unidirectional motion for sufficiently close vesicle-to-spine diameter ratios. Our analysis predicts the critical vesicle-to-spine diameter ratio, at which there is a transition from unidirectional to bidirectional motion, consistent with experimental observations of vesicle trajectories in the literature.

Membrane transport is generally thought to occur via an alternating access mechanism in which the transporter adopts at least two states, accessible from two different sides of the membrane to exchange substrates from the extracellular environment and the cytoplasm or from the cytoplasm and the intracellular matrix of the organelles (only in eukaryotes). In recent years, a number of high resolution structures have supported this general framework for a wide class of transport molecules, although additional states along the transport pathway are emerging as critically important. Given that substrate binding is often weak in order to enhance overall transport rates, there exists the distinct possibility that transporters may transport the incorrect substrate. This is certainly the case for many pharmaceutical compounds that are absorbed in the gut or cross the blood brain barrier through endogenous transporters. Docking studies on the bacterial sugar transporter vSGLT reveal that many highly toxic compounds are compatible with binding to the orthosteric site, further motivating the selective pressure for additional modes of selectivity. Motivated by recent work in which we observed failed substrate delivery in a molecular dynamics simulation where the energized ion still goes down its concentration gradient, we hypothesize that some transporters evolved to harness this 'slip' mechanism to increase substrate selectivity and reduce the uptake of toxic molecules. Here, we test this idea by constructing and exploring a kinetic transport model that includes a slip pathway. While slip reduces the overall productive flux, when coupled with a second toxic molecule that is more prone to slippage, the overall substrate selectivity dramatically increases, suppressing the accumulation of the incorrect compound. We show that the mathematical framework for increased substrate selectivity in our model is analogous to the classic proofreading mechanism originally proposed for tRNA synthase; however, because the transport cycle is reversible we identified conditions in which the selectivity is essentially infinite and incorrect substrates are exported from the cell in a 'detoxification' mode. The cellular consequences of proofreading and membrane slippage are discussed as well as the impact on future drug development., Competing Interests: The authors have declared that no competing interests exist.

The cortical collecting duct (CCD) of the mammalian kidney plays a major role in the maintenance of total body electrolyte, acid/base, and fluid homeostasis by tubular reabsorption and excretion. The mammalian CCD is heterogeneous, composed of Na + -absorbing principal cells (PCs) and acid-base-transporting intercalated cells (ICs). Perturbations in luminal flow rate alter hydrodynamic forces to which these cells in the cylindrical tubules are exposed. However, most studies of tubular ion transport have been performed in cell monolayers grown on or epithelial sheets affixed to a flat support, since analysis of transepithelial transport in native tubules by in vitro microperfusion requires considerable expertise. Here, we report on the generation and characterization of an in vitro, perfusable three-dimensional kidney CCD model (3D CCD), in which immortalized mouse PC-like mpkCCD cells are seeded within a cylindrical channel embedded within an engineered extracellular matrix and subjected to luminal fluid flow. We find that a tight epithelial barrier composed of differentiated and polarized PCs forms within 1 wk. Immunofluorescence microscopy reveals the apical epithelial Na + channel ENaC and basolateral Na + /K + -ATPase. On cessation of luminal flow, benzamil-inhibitable cell doming is observed within these 3D CCDs consistent with the presence of ENaC-mediated Na + absorption. Our 3D CCD provides a geometrically and microphysiologically relevant platform for studying the development and physiology of renal tubule segments.

A combined experimental and computational model approach was developed to assess heat effects on drug delivery from transdermal delivery systems (TDSs) in vitro and nicotine was the model drug. A Franz diffusion cell system was modified to allow close control of skin temperature when heat was applied from an infrared lamp in vitro. The effects of different heat application regimens on nicotine fluxes from two commercial TDSs across human cadaver skin were determined. Results were interpreted in terms of transport parameters estimated using a computational heat and mass transport model. Steady-state skin surface temperature was obtained rapidly after heat application. Increasing skin surface temperature from 32 to 42°C resulted in an approximately 2-fold increase in average nicotine flux for both TDSs, with maximum flux observed during early heat application. ANOVA statistical analyses of the in vitro permeation data identified TDS differences, further evidenced by the need for a two-layer model to describe one of the TDSs. Activation energies associated with these data suggest similar temperature effects on nicotine transport across the skin despite TDS design differences. Model simulations based on data obtained from continuous heat application were able to predict system response to intermittent heat application, as shown by the agreement between the simulation results and experimental data of nicotine fluxes under four different heat application regimens. The combination of in vitro permeation testing and a computational model provided a parameter-based heat and mass transport approach to evaluate heat effects on nicotine TDS delivery.

Cell polarity is an important cellular process that cells use for various cellular functions such as asymmetric division, cell migration, and directionality determination. In asymmetric cell division, a mother cell creates multiple polarities of various proteins simultaneously within her membrane and cytosol to generate two different daughter cells. The formation of multiple polarities in asymmetric cell division has been found to be controlled via the regulatory system by upstream polarity of the membrane to downstream polarity of the cytosol, which is involved in not only polarity establishment but also polarity positioning. However, the mechanism for polarity positioning remains unclear. In this study, we found a general mechanism and mathematical structure for the multiple streams of polarities to determine their relative position via conceptional models based on the biological example of the asymmetric cell division process of C. elegans embryo. Using conceptional modeling and model reductions, we show that the positional relation of polarities is determined by a contrasting role of regulation by upstream polarity proteins on the transition process of diffusion dynamics of downstream proteins. We analytically prove that our findings hold under the general mathematical conditions, suggesting that the mechanism of relative position between upstream and downstream dynamics could be understood without depending on a specific type of bio-chemical reaction, and it could be the universal mechanism in multiple streams of polarity dynamics of the cell.

Creatinine is widely used as a biomarker of glomerular filtration, and, hence, renal function. However, transporter-mediated secretion also contributes to its renal clearance, albeit to a lesser degree. Inhibition of these transporters causes transient serum creatinine elevation, which can be mistaken as impaired renal function. The current study developed mechanistic models of creatinine kinetics within physiologically based framework accounting for multiple transporters involved in creatinine renal elimination, assuming either unidirectional or bidirectional-OCT2 transport (driven by electrochemical gradient). Robustness of creatinine models was assessed by predicting creatinine-drug interactions with 10 perpetrators; performance evaluation accounted for 5% intra-individual variability in serum creatinine. Models showed comparable predictive performances of the maximum steady-state effect regardless of OCT2 directionality assumptions. However, only the bidirectional-OCT2 model successfully predicted the minimal effect of ranitidine. The dynamic nature of models provides clear advantage to static approaches and most advanced framework for evaluating interplay between multiple processes in creatinine renal disposition., (© 2020 The Authors. CPT: Pharmacometrics & Systems Pharmacology published by Wiley Periodicals, Inc. on behalf of the American Society for Clinical Pharmacology and Therapeutics.)

The structure of tubular transport networks is thought to underlie much of biological regularity, from individuals to ecosystems. A core assumption of transport network models is either area-preserving or area-increasing branching, such that the summed cross-sectional area of all child branches is equal to or greater than the cross-sectional area of their respective parent branch. For insects, the most diverse group of animals, the assumption of area-preserving branching of tracheae is, however, based on measurements of a single individual and an assumption of gas exchange by diffusion. Here we show that ants exhibit neither area-preserving nor area-increasing branching in their abdominal tracheal systems. We find for 20 species of ants that the sum of child tracheal cross-sectional areas is typically less than that of the parent branch (area-decreasing). The radius, rather than the area, of the parent branch is conserved across the sum of child branches. Interpretation of the tracheal system as one optimized for the release of carbon dioxide, while readily catering to oxygen demand, explains the branching pattern. Our results, together with widespread demonstration that gas exchange in insects includes, and is often dominated by, convection, indicate that for generality, network transport models must include consideration of systems with different architectures., Competing Interests: The authors have declared that no competing interests exist.

The human central nervous system (CNS) vasculature expresses a distinctive barrier phenotype, the blood-brain barrier (BBB). As the BBB contributes to low efficiency in CNS pharmacotherapy by restricting drug transport, the development of an in vitro human BBB model has been in demand. Here, we present a microfluidic model of CNS angiogenesis having three-dimensional (3D) lumenized vasculature in concert with perivascular cells. We confirmed the necessity of the angiogenic tri-culture system (brain endothelium in direct interaction with pericytes and astrocytes) to attain essential phenotypes of BBB vasculature, such as minimized vessel diameter and maximized junction expression. In addition, lower vascular permeability is achieved in the tri-culture condition compared to the monoculture condition. Notably, we focussed on reconstituting the functional efflux transporter system, including p-glycoprotein (p-gp), which is highly responsible for restrictive drug transport. By conducting the calcein-AM efflux assay on our 3D perfusable vasculature after treatment of efflux transporter inhibitors, we confirmed the higher efflux property and prominent effect of inhibitors in the tri-culture model. Taken together, we designed a 3D human BBB model with functional barrier properties based on a developmentally inspired CNS angiogenesis protocol. We expect the model to contribute to a deeper understanding of pathological CNS angiogenesis and the development of effective CNS medications., (© 2019 Wiley Periodicals, Inc.)

In vivo, eukaryotic cells are embedded in a matrix environment, where they grow and develop. Generally, this extracellular matrix (ECM) is an anisotropic fibrous structure, through which macromolecules and biochemical signaling molecules at the nanometer scale diffuse. The ECM is continuously remodeled by cells, via mechanical interactions, which lead to a potential link between biomechanical and biochemical cell-cell interactions. Here, it is studied how cell-induced forces applied on the ECM impact the biochemical transport of molecules between distant cells. It is experimentally observed that cells remodel the ECM by increasing fiber alignment and density of the matrix between them over time. Using random walk simulations on a 3D lattice, elongated fixed obstacles are implemented that mimic the fibrous ECM structure. Both diffusion of a tracer molecule and the mean first-passage time a molecule secreted from one cell takes to reach another cell are measured. The model predicts that cell-induced remodeling can lead to a dramatic speedup in the transport of molecules between cells. Fiber alignment and densification cause reduction of the transport dimensionality from a 3D to a much more rapid 1D process. Thus, a novel mechanism of mechano-biochemical feedback in the regulation of long-range cell-cell communication is suggested., (© 2019 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Regulation of molecular transport via intercellular channels called plasmodesmata (PDs) is important for both coordinating developmental and environmental responses among neighbouring cells, and isolating (groups of) cells to execute distinct programs. Cell-to-cell mobility of fluorescent molecules and PD dimensions (measured from electron micrographs) are both used as methods to predict PD transport capacity (i.e., effective symplasmic permeability), but often yield very different values. Here, we build a theoretical bridge between both experimental approaches by calculating the effective symplasmic permeability from a geometrical description of individual PDs and considering the flow towards them. We find that a dilated central region has the strongest impact in thick cell walls and that clustering of PDs into pit fields strongly reduces predicted permeabilities. Moreover, our open source multi-level model allows to predict PD dimensions matching measured permeabilities and add a functional interpretation to structural differences observed between PDs in different cell walls., Competing Interests: ED, BM, YB No competing interests declared, (© 2019, Deinum et al.)

For many years, double-layer phospholipid membrane vesicles, released by most cells, were not considered to be of biological significance. This stance has dramatically changed with the recognition of extracellular vesicles (EVs) as carriers of biologically active molecules that can traffic to local or distant targets and execute defined biological functions. The dimensionality of the field has expanded with the appreciation of diverse types of EVs and the complexity of vesicle biogenesis, cargo loading, release pathways, targeting mechanisms, and vesicle processing. With the expanded interest in the field and the accelerated rate of publications on EV structure and function in diverse biomedical fields, it has become difficult to distinguish between well-established biological features of EV and the untested hypotheses or speculative assumptions that await experimental proof. With the growing interest despite the limited evidence, we sought in this essay to formulate a set of unsolved mysteries in the field, sort out established data from fascinating hypotheses, and formulate several challenging questions that must be answered for the field to advance., Competing Interests: Dr. Margolis has declared that no competing interests exist. Dr. Sadovsky is a member of a Clinical Expert Advisory Panel to Illumina, Inc.

It is often assumed in biophysical studies that when multiple identical molecular motors interact with two parallel microtubules, the microtubules will be crosslinked and locked together. The aim of this study is to examine this assumption mathematically. We model the forces and movements generated by motors with a time-continuous Markov process and find that, counter-intuitively, a tug-of-war results from opposing actions of identical motors bound to different microtubules. The model shows that many motors bound to the same microtubule generate a great force applied to a smaller number of motors bound to another microtubule, which increases detachment rate for the motors in minority, stabilizing the directional sliding. However, stochastic effects cause occasional changes of the sliding direction, which has a profound effect on the character of the long-term microtubule motility, making it effectively diffusion-like. Here, we estimate the time between the rare events of switching direction and use them to estimate the effective diffusion coefficient for the microtubule pair. Our main result is that parallel microtubules interacting with multiple identical motors are not locked together, but rather slide bidirectionally. We find explicit formulae for the time between directional switching for various motor numbers.

Cancers exhibit spatially heterogeneous, unique vascular architectures across individual samples, cell-lines and patients. This inherently disorganised collection of leaky blood vessels contribute significantly to suboptimal treatment efficacy. Preclinical tools are urgently required which incorporate the inherent variability and heterogeneity of tumours to optimise and engineer anti-cancer therapies. In this study, we present a novel computational framework which incorporates whole, realistic tumours extracted ex vivo to efficiently simulate vascular blood flow and interstitial fluid transport in silico for validation against in vivo biomedical imaging. Our model couples Poiseuille and Darcy descriptions of vascular and interstitial flow, respectively, and incorporates spatially heterogeneous blood vessel lumen and interstitial permeabilities to generate accurate predictions of tumour fluid dynamics. Our platform enables highly-controlled experiments to be performed which provide insight into how tumour vascular heterogeneity contributes to tumour fluid transport. We detail the application of our framework to an orthotopic murine glioma (GL261) and a human colorectal carcinoma (LS147T), and perform sensitivity analysis to gain an understanding of the key biological mechanisms which determine tumour fluid transport. Finally we mimic vascular normalization by modifying parameters, such as vascular and interstitial permeabilities, and show that incorporating realistic vasculatures is key to modelling the contrasting fluid dynamic response between tumour samples. Contrary to literature, we show that reducing tumour interstitial fluid pressure is not essential to increase interstitial perfusion and that therapies should seek to develop an interstitial fluid pressure gradient. We also hypothesise that stabilising vessel diameters and permeabilities are not key responses following vascular normalization and that therapy may alter interstitial hydraulic conductivity. Consequently, we suggest that normalizing the interstitial microenvironment may provide a more effective means to increase interstitial perfusion within tumours., Competing Interests: The authors have declared that no competing interests exist.

Polar auxin transport lies at the core of many self-organizing phenomena sustaining continuous plant organogenesis. In angiosperms, the shoot apical meristem is a potentially unique system in which the two main modes of auxin-driven patterning-convergence and canalization-co-occur in a coordinated manner and in a fully three-dimensional geometry. In the epidermal layer, convergence points form, from which auxin is canalized towards inner tissue. Each of these two patterning processes has been extensively investigated separately, but the integration of both in the shoot apical meristem remains poorly understood. We present here a first attempt of a three-dimensional model of auxin-driven patterning during phyllotaxis. We base our simulations on a biochemically plausible mechanism of auxin transport proposed by Cieslak et al. (2015) which generates both convergence and canalization patterns. We are able to reproduce most of the dynamics of PIN1 polarization in the meristem, and we explore how the epidermal and inner cell layers act in concert during phyllotaxis. In addition, we discuss the mechanism by which initiating veins connect to the already existing vascular system., Competing Interests: The authors have declared that no competing interests exist.

Motivation: Single cells often show stochastic behaviour and variations in the physiological state of individual cells affect the behaviour observed in cell populations. This may be partially explained by variations in the concentration and spatial location of molecules within and in the vicinity of each cell., Methods: This paper introduces an agent-based model that represents single-molecule transport through the cellular envelope of Escherichia coli at the micrometre scale. This model enables broader observation of molecular transport throughout the different membrane layers and the study of the effect of molecular concentration in cellular noise. Simulations considered various low molecular weight molecules, i.e. ampicillin, bosentan, coumarin, saquinavir, and terbutaline, and a gradient of molecular concentrations. The model ensured stochasticity in the location of the agents, using diffusing spherical particles with physical dimensions., Results: Simulation results were validated against theoretical and experimental data. For example, theoretically, ampicillin molecules take 0.6 s to cross the entire cell envelope, and computational simulations took 0.68 s, 0.68 s, 0.70 s, and 0.69 s, for concentrations of 1.44 μM, 13.21 μM, 26.4 μM and 105.61 μM, respectively. Replicate standard deviation decreased with growing initial concentrations of the molecules. In turn, no clear relationship could be observed between molecular size and variability., Conclusions: This work presented a novel agent-based model to study the effect of the initial concentration of low molecular weight molecules on cellular noise. Cellular noise during molecule diffusion was found to be concentration-dependent and size-independent. The new model holds considerable potential for future, more complex analyses, when emerging experimental data may enable modelling of membrane transport mechanisms., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Many microorganisms secrete molecules that interact with resources outside of the cell. This includes, for example, enzymes that degrade polymers like chitin, and chelators that bind trace metals like iron. In contrast to direct uptake via the cell surface, such release strategies entail the risk of losing the secreted molecules to environmental sinks, including 'cheating' genotypes. Nevertheless, such secretion strategies are widespread, even in the well-mixed marine environment. Here, we investigate the benefits of a release strategy whose efficiency has frequently been questioned: iron uptake in the ocean by secretion of iron chelators called siderophores. We asked the question whether the release itself is essential for the function of siderophores, which could explain why this risky release strategy is widespread. We developed a reaction-diffusion model to determine the impact of siderophore release on iron uptake from the predominant iron sources in marine environments, colloidal or particulate iron, formed due to poor iron solubility. We found that release of siderophores is essential to accelerate iron uptake, as secreted siderophores transform slowly diffusing large iron particles to small, quickly diffusing iron-siderophore complexes. In addition, we found that cells can synergistically share their siderophores, depending on their distance and the size of the iron sources. Our study helps understand why release of siderophores is so widespread: even though a large fraction of siderophores is lost, the solubilization of iron through secreted siderophores can efficiently increase iron uptake, especially if siderophores are produced cooperatively by several cells. Overall, resource uptake mediated via release of molecules transforming their substrate could be essential to overcome diffusion limitation specifically in the cases of large, aggregated resources. In addition, we find that including the reaction of the released molecule with the substrate can impact the result of cooperative and competitive interactions, making our model also relevant for release-based uptake of other substrates.

Cells are neither flat nor smooth, which has serious implications for prevailing plasma membrane models and cellular processes like cell signalling, adhesion and molecular clustering. Using probability distributions from diffusion simulations, we demonstrate that 2D and 3D Euclidean distance measurements substantially underestimate diffusion on non-flat surfaces. Intuitively, the shortest within surface distance (SWSD), the geodesic distance, should reduce this problem. The SWSD is accurate for foldable surfaces but, although it outperforms 2D and 3D Euclidean measurements, it still underestimates movement on deformed surfaces. We demonstrate that the reason behind the underestimation is that topographical features themselves can produce both super- and subdiffusion, i.e. the appearance of anomalous diffusion. Differentiating between topography-induced and genuine anomalous diffusion requires characterising the surface by simulating Brownian motion on high-resolution cell surface images and a comparison with the experimental data., Competing Interests: The authors declare no competing interests.

The retinal pigment epithelium (RPE) is the outermost cell layer of the retina. It has several important physiological functions, among which is removal of excess fluid from the sub-retinal space by pumping it isotonically towards the choroid. Failure of this pumping leads to fluid accumulation, which is closely associated with several pathological conditions, such as age-related macular degeneration, macular oedema and retinal detachment. In the present work we study mechanisms responsible for fluid transport across the RPE with the aim of understanding how fluid accumulation can be prevented. We focus on two possible mechanisms, osmosis and electroosmosis, and develop a spatially resolved mathematical model that couples fluid and ion transport across the epithelium, accounting for the presence of Na + ,K + and Cl - ions. Our model predicts spatial variability of ion concentrations and the electrical potential along the cleft gap between two adjacent cells, which osmotically drives the flow across the lateral membranes. This flow is directed from the sub-retinal space to the choroid and has a magnitude close to measured values. Electroosmosis is subdominant by three orders of magnitude to osmosis and has an opposite direction, suggesting that local osmosis is the main driving mechanism for water transport across the RPE., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Current theory and supporting research suggests that radial transport is the most limiting factor to root water uptake, raising the question whether only absorbing root length and radial conductivity matter to water uptake. Here, we extended the porous pipe analytical model of root water uptake to entire root networks in 3D and analysed the relative importance of axial and radial characteristics to total uptake over parameter ranges reported in the literature. We found that network conductance can be more sensitive to axial than radial conductance of absorbing roots. When axial transport limits uptake, more dichotomous topology, especially towards the base of the network, increases water uptake efficiency, while the effect of root length is reduced. Whole root system conductance was sensitive to radial transport and length in model lupin (Lupinus angustifolius L.), but to axial transport and topology in wheat (Triticum aestivum L.), suggesting the root habit niche space of monocots may be constrained by their loss of secondary growth. A deep tap root calibrated to oak (Quercus fusiformis J. Buchholz) hydraulic parameters required 15 times more xylem volume to transport comparable amounts of water once recalibrated to parameters from juniper (Juniperus ashei Small 1901), showing that anatomical constraints on axial conductance can lead to significant trade-offs in woody roots as well. Root system water uptake responds to axial transport and can be limited by it in a biologically meaningful way., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Vessel length is an important functional trait for plant hydraulics, because it determines the ratio of flow resistances posed by lumen and pit membranes and hence controls xylem hydraulic efficiency. The most commonly applied methods to estimate vessel lengths are based on the injection of silicon or paint into cut-off stem segments. The number of stained vessels in a series of cross-sections in increasing distance from the injection point is then counted. The resulting infusion profiles are used to estimate the vessel length distribution using one of several statistical algorithms. However, the basis of these algorithms has not been systematically analysed using probability theory. We derive a general mathematical expression for the expected shape of the infusion profile for a given vessel length distribution, provide analytic solutions for five candidate distributions (exponential, Erlang(2), gamma, Weibull, and log-normal), and present maximum likelihood estimators for the parameters of these distributions including implementations in R based on two potential sampling schemes (counting all injected vessels or counting the injected and empty vessels in a random subset of each cross-section). We then explore the performance of these estimators relative to other methods with Monte Carlo experiments. Our analysis demonstrates that most published methods estimate the conditional length distribution of vessels that cross an injection point, which is a size-biased version of the overall length distribution in the stem. We show the mathematical relationship between these distributions and provide methods to estimate either of them. According to our simulation experiments, vessel length distribution was best described by the more flexible models, especially the Weibull distribution. In simulations, the estimators were able to recover the parameters of the vessel length distribution if its functional form was known, achieving an overlap of 90% or more between the true and predicted length distribution when counting no more than 500 injected vessels in 10 cross-sections. This sample size nowadays can easily be reached with the help of automated image analysis., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

The primary objective of this study is to introduce a simple and flexible mathematical approach which models transport processes in skin using compartments. The main feature of the presented approach is that the rate constants for exchange between compartments are derived from physiologically relevant diffusional transport parameters. This allows for better physical interpretation of the rate constants, and limits the number of parameters for the compartmental model. The resulting compartmental solution is in good agreement with previously published solutions for the diffusion model of skin when ten or more compartments are used. It was found that the new compartmental model with three compartments provided a better fit of the previously publish water penetration data than the diffusion model. Two special cases for which it is difficult to implement the diffusion model were considered using our compartmental approach. In both cases the compartmental model predictions agreed well with the diffusion model., (Copyright © 2018 Elsevier B.V. All rights reserved.)

One of the basic and vital processes in living organisms is mass exchange, which occurs on several levels: it goes from blood vessels to cells and organelles within cells. On that path, molecules, as oxygen, metabolic products, drugs, etc. Traverse different macro and micro environments - blood, extracellular/intracellular space, and interior of organelles; and also biological barriers such as walls of blood vessels and membranes of cells and organelles. Many aspects of this mass transport remain unknown, particularly the biophysical mechanisms governing drug delivery. The main research approach relies on laboratory and clinical investigations. In parallel, considerable efforts have been directed to develop computational tools for additional insight into the intricate process of mass exchange and transport. Along these lines, we have recently formulated a composite smeared finite element (CSFE) which is composed of the smeared continuum pressure and concentration fields of the capillary and lymphatic system, and of these fields within tissue. The element offers an elegant and simple procedure which opens up new lines of inquiry and can be applied to large systems such as organs and tumors models. Here, we extend this concept to a multiscale scheme which concurrently couples domains that span from large blood vessels, capillaries and lymph, to cell cytosol and further to organelles of nanometer size. These spatial physical domains are coupled by the appropriate connectivity elements representing biological barriers. The composite finite element has "degrees of freedom" which include pressures and concentrations of all compartments of the vessels-tissue assemblage. The overall model uses the standard, measurable material properties of the continuum biological environments and biological barriers. It can be considered as a framework into which we can incorporate various additional effects (such as electrical or biochemical) for transport through membranes or within cells. This concept and the developed FE software within our package PAK offers a computational tool that can be applied to whole-organ systems, while also including specific domains such as tumors. The solved examples demonstrate the accuracy of this model and its applicability to large biological systems., (Copyright © 2018. Published by Elsevier Ltd.)

Recent advances in molecular biology and fluorescence microscopy imaging have made possible the inference of the dynamics of molecules in living cells. Such inference allows us to understand and determine the organization and function of the cell. The trajectories of particles (e.g., biomolecules) in living cells, computed with the help of object tracking methods, can be modeled with diffusion processes. Three types of diffusion are considered: (i) free diffusion, (ii) subdiffusion, and (iii) superdiffusion. The mean-square displacement (MSD) is generally used to discriminate the three types of particle dynamics. We propose here a nonparametric three-decision test as an alternative to the MSD method. The rejection of the null hypothesis, i.e., free diffusion, is accompanied by claims of the direction of the alternative (subdiffusion or superdiffusion). We study the asymptotic behavior of the test statistic under the null hypothesis and under parametric alternatives which are currently considered in the biophysics literature. In addition, we adapt the multiple-testing procedure of Benjamini and Hochberg to fit with the three-decision-test setting, in order to apply the test procedure to a collection of independent trajectories. The performance of our procedure is much better than the MSD method as confirmed by Monte Carlo experiments. The method is demonstrated on real data sets corresponding to protein dynamics observed in fluorescence microscopy.

A three-dimensional diffusion model with either hexagonal or cylindrical symmetry has been constructed to simulate desorption profiles of hydrophilic chemicals from the topmost layer of human skin (the stratum corneum) as measured in ex vivo studies. The tissue is pierced by skin appendages - sweat glands and hair follicles - which in this particular scenario are considered to be perfect sinks. Desorption profiles of nine test permeants covering a wide range of lipophilicity were analyzed. By optimizing transverse and lateral diffusion coefficients to match these profiles, it was found that the lateral diffusivity values exceeded the transverse values by average factors ranging from 45 (hexagon model) to 71 (cylinder model). However, transverse clearance exceeded lateral clearance by factors ranging from 8 to 27 (cylinder model); these values were strongly influenced by the thickness of the individual tissue samples, as expected. The results confirm the validity of earlier estimates of transverse diffusivity of hydrophilic compounds in human stratum corneum based on purely one-dimensional models. They furthermore confirm that transcellular transport is an important component of the stratum corneum's polar pathway, in addition to the already-recognized appendageal transport mechanism., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Although various drugs, environmental pollutants and nanoparticles (NP) can cross the human placental barrier and may harm the developing fetus, knowledge on predictive placental transfer rates and the underlying transport pathways is mostly lacking. Current available in vitro placental transfer models are often inappropriate for translocation studies of macromolecules or NPs and do not consider barrier function of placental endothelial cells (EC). Therefore, we developed a human placental in vitro co-culture transfer model with tight layers of trophoblasts (BeWo b30) and placental microvascular ECs (HPEC-A2) on a low-absorbing, 3 µm porous membrane. Translocation studies with four model substances and two polystyrene (PS) NPs across the individual and co-culture layers revealed that for most of these compounds, the trophoblast and the EC layer both demonstrate similar, but not additive, retention capacity. Only the paracellular marker Na-F was substantially more retained by the BeWo layer. Furthermore, simple shaking, which is often applied to mimic placental perfusion, did not alter translocation kinetics compared to static exposure. In conclusion, we developed a novel placental co-culture model, which provides predictive values for translocation of a broad variety of molecules and NPs and enables valuable mechanistic investigations on cell type-specific placental barrier function.

To model non-Markovian fluctuations arising in biomolecular transport, we introduce a stochastic process with memory where Brownian motion is modulated sinusoidally. The probability density function and moments of this non-Markovian process are evaluated analytically as Hida stochastic functional integrals. Comparison of graphs of computed variance vis-á-vis empirical data for protein diffusion coefficients closely match with both exhibiting emergent superdiffusive then subdiffusive behavior for longer proteins., (Copyright © 2018 Elsevier Inc. All rights reserved.)

In this work, we present structured capillaries that were inspired by the microstructures of the external scent efferent system as found in different European true bug species (Pentatomidae and Cydnidae). These make use of small, orientated structures in order to facilitate fluid movement towards desired areas where defensive substances are evaporated. Gland channels and microstructures were investigated by means of scanning electron microscopy and abstracted into three-dimensional models. We used these models to create scent channel replicas from different technical substrates (steel and polymers) by means of laser ablation, laser structuring and casting. Video analysis of conducted fluid-flow experiments showed that bug-inspired, artificial scent fluid channels can indeed transport different fluids (water solutions and oils/lubricants) passively in one direction (velocities of up to 1 mm s -1 ), while halting the fluid movement in the opposite direction. At the end of this contribution, we present a physical theory that explains the observed fluid transport and sets the rules for performance optimization in future work., (© 2018 The Authors.)

Cells typically vary in their response to extracellular ligands. Receptor transport processes modulate ligand-receptor induced signal transduction and impact the variability in cellular responses. Here, we quantitatively characterized cellular variability in erythropoietin receptor (EpoR) trafficking at the single-cell level based on live-cell imaging and mathematical modeling. Using ensembles of single-cell mathematical models reduced parameter uncertainties and showed that rapid EpoR turnover, transport of internalized EpoR back to the plasma membrane, and degradation of Epo-EpoR complexes were essential for receptor trafficking. EpoR trafficking dynamics in adherent H838 lung cancer cells closely resembled the dynamics previously characterized by mathematical modeling in suspension cells, indicating that dynamic properties of the EpoR system are widely conserved. Receptor transport processes differed by one order of magnitude between individual cells. However, the concentration of activated Epo-EpoR complexes was less variable due to the correlated kinetics of opposing transport processes acting as a buffering system.

Phloem transport is the process by which carbohydrates produced by photosynthesis in the leaves get distributed in a plant. According to Münch, the osmotically generated hydrostatic phloem pressure is the force driving the long-distance transport of photoassimilates. Following Thompson and Holbrook[35]'s approach, we develop a mathematical model of coupled water-carbohydrate transport. It is first proven that the model presented here preserves the positivity. The model is then applied to simulate the flow of phloem sap for an organic tree shape, on a 3D surface and in a channel with orthotropic hydraulic properties. Those features represent an significant advance in modelling the pathway for carbohydrate transport in trees.

We propose a mechanism for the homeostatic control of synapses along the ventral cord of Caenorhabditis elegans during development, based on a form of Turing pattern formation on a growing domain. C. elegans is an important animal model for understanding cellular mechanisms underlying learning and memory. Our mathematical model consists of two interacting chemical species, where one is passively diffusing and the other is actively trafficked by molecular motors, which switch between forward and backward moving states (bidirectional transport). This differs significantly from the standard mechanism for Turing pattern formation based on the interaction between fast and slow diffusing species. We derive evolution equations for the chemical concentrations on a slowly growing one-dimensional domain, and use numerical simulations to demonstrate the insertion of new concentration peaks as the length increases. Taking the passive component to be the protein kinase CaMKII and the active component to be the glutamate receptor GLR-1, we interpret the concentration peaks as sites of new synapses along the length of C. elegans, and thus show how the density of synaptic sites can be maintained.

Stem cells are identical in many scales, they share the same molecular composition, DNA, genes, and genetic networks, yet they should acquire different properties to form a functional tissue. Therefore, they must interact and get some external information from their environment, either spatial (dynamical fields) or temporal (lineage). In this paper we test to what extent coupled chemical and physical fields can underlie the cell's positional information during development. We choose the root apical meristem of Arabidopsis thaliana to model the emergence of cellular patterns. We built a model to study the dynamics and interactions between the cell divisions, the local auxin concentration, and physical elastic fields. Our model recovers important aspects of the self-organized and resilient behavior of the observed cellular patterns in the Arabidopsis root, in particular, the reverse fountain pattern observed in the auxin transport, the PIN-FORMED (protein family of auxin transporters) polarization pattern and the accumulation of auxin near the region of maximum curvature in a bent root. Our model may be extended to predict altered cellular patterns that are expected under various applied auxin treatments or modified physical growth conditions.

The sites of evaporation within leaves are unknown, but they have drawn attention for decades due to their perceived implications for many factors, including patterns of leaf isotopic enrichment, the maintenance of mesophyll water status, stomatal regulation, and the interpretation of measured stomatal and leaf hydraulic conductances. We used a spatially explicit model of coupled water and heat transport outside the xylem, MOFLO 2.0, to map the distribution of net evaporation across leaf tissues in relation to anatomy and environmental parameters. Our results corroborate earlier predictions that most evaporation occurs from the epidermis at low light and moderate humidity but that the mesophyll contributes substantially when the leaf center is warmed by light absorption, and more so under high humidity. We also found that the bundle sheath provides a significant minority of evaporation (15% in darkness and 18% in high light), that the vertical center of amphistomatous leaves supports net condensation, and that vertical temperature gradients caused by light absorption vary over 10-fold across species, reaching 0.3°C. We show that several hypotheses that depend on the evaporating sites require revision in light of our findings, including that experimental measurements of stomatal and hydraulic conductances should be affected directly by changes in the location of the evaporating sites. We propose a new conceptual model that accounts for mixed-phase water transport outside the xylem. These conclusions have far-reaching implications for inferences in leaf hydraulics, gas exchange, water use, and isotope physiology., (© 2017 American Society of Plant Biologists. All Rights Reserved.)

Optimization of fluid transport in the slime mold Physarum polycephalum has been the subject of several modeling efforts in recent literature. Existing models assume that the tube adaptation mechanism in P. polycephalum's tubular network is controlled by the sheer amount of fluid flow through the tubes. We put forward the hypothesis that the controlling variable may instead be the flow's pressure gradient along the tube. We carry out the stability analysis of such a revised mathematical model for a parallel-edge network, proving that the revised model supports the global flow-optimizing behavior of the slime mold for a substantially wider class of response functions compared to previous models. Simulations also suggest that the same conclusion may be valid for arbitrary network topologies.

The lymphatic system is responsible for transporting interstitial fluid back to the bloodstream, but unlike the cardiovascular system, lacks a centralized pump-the heart-to drive flow. Instead, each collecting lymphatic vessel can individually contract and dilate producing unidirectional flow enforced by intraluminal check valves. Due to the large number and spatial distribution of such pumps, high-level coordination would be unwieldy. This leads to the question of how each segment of lymphatic vessel responds to local signals that can contribute to the coordination of pumping on a network basis. Beginning with elementary fluid mechanics and known cellular behaviors, we show that two complementary oscillators emerge from i) mechanical stretch with calcium ion transport and ii) fluid shear stress induced nitric oxide production (NO). Using numerical simulation and linear stability analysis we show that the newly identified shear-NO oscillator shares similarities with the well-known Van der Pol oscillator, but has unique characteristics. Depending on the operating conditions, the shear-NO process may i) be inherently stable, ii) oscillate spontaneously in response to random disturbances or iii) synchronize with weak periodic stimuli. When the complementary shear-driven and stretch-driven oscillators interact, either may dominate, producing a rich family of behaviors similar to those observed in vivo., Competing Interests: The authors have declared that no competing interests exist.

The goal of this study was to investigate water and solute transport, with a focus on sodium transport (T Na ) and metabolism along individual nephron segments under differing physiological and pathophysiological conditions. To accomplish this goal, we developed a computational model of solute transport and oxygen consumption (Q O 2 ) along different nephron populations of a rat kidney. The model represents detailed epithelial and paracellular transport processes along both the superficial and juxtamedullary nephrons, with the loop of Henle of each model nephron extending to differing depths of the inner medulla. We used the model to assess how changes in T Na may alter Q O 2 in different nephron segments and how shifting the T Na sites alters overall kidney Q O 2 Under baseline conditions, the model predicted a whole kidney T Na /Q O 2 , which denotes the number of moles of Na + reabsorbed per moles of O 2 consumed, of ∼15, with T Na efficiency predicted to be significantly greater in cortical nephron segments than in medullary segments. The T Na /Q O 2 ratio was generally similar among the superficial and juxtamedullary nephron segments, except for the proximal tubule, where T Na /Q O 2 was ∼20% higher in superficial nephrons, due to the larger luminal flow along the juxtamedullary proximal tubules and the resulting higher, flow-induced transcellular transport. Moreover, the model predicted that an increase in single-nephron glomerular filtration rate does not significantly affect T Na /Q O 2 in the proximal tubules but generally increases T Na /Q O 2 along downstream segments. The latter result can be attributed to the generally higher luminal [Na + ], which raises paracellular T Na Consequently, vulnerable medullary segments, such as the S3 segment and medullary thick ascending limb, may be relatively protected from flow-induced increases in Q O 2 under pathophysiological conditions., (Copyright © 2016 the American Physiological Society.)

Transport of solutes helps to regulate normal physiology and proper function of cartilage in diarthrodial joints. Multiple studies have shown the effects of characteristic parameters such as concentration of proteoglycans and collagens and the orientation of collagen fibrils on the diffusion process. However, not much quantitative information and accurate models are available to help understand how the characteristics of the fluid surrounding articular cartilage influence the diffusion process. In this study, we used a combination of micro-computed tomography experiments and biphasic-solute finite element models to study the effects of three parameters of the overlying bath on the diffusion of neutral solutes across cartilage zones. Those parameters include bath size, degree of stirring of the bath, and the size and concentration of the stagnant layer that forms at the interface of cartilage and bath. Parametric studies determined the minimum of the finite bath size for which the diffusion behavior reduces to that of an infinite bath. Stirring of the bath proved to remarkably influence neutral solute transport across cartilage zones. The well-stirred condition was achieved only when the ratio of the diffusivity of bath to that of cartilage was greater than ≈1000. While the thickness of the stagnant layer at the cartilage-bath interface did not significantly influence the diffusion behavior, increase in its concentration substantially elevated solute concentration in cartilage. Sufficient stirring attenuated the effects of the stagnant layer. Our findings could be used for efficient design of experimental protocols aimed at understanding the transport of molecules across articular cartilage., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The development of indefinitely propagating human 'mini-guts' has led to a rapid advance in gastrointestinal research related to transport physiology, developmental biology, pharmacology, and pathophysiology. These mini-guts, also called enteroids or colonoids, are derived from LGR5 + intestinal stem cells isolated from the small intestine or colon. Addition of WNT3A and other growth factors promotes stemness and results in viable, physiologically functional human intestinal or colonic cultures that develop a crypt-villus axis and can be differentiated into all intestinal epithelial cell types. The success of research using human enteroids has highlighted the limitations of using animals or in vitro, cancer-derived cell lines to model transport physiology and pathophysiology. For example, curative or preventive therapies for acute enteric infections have been limited, mostly due to the lack of a physiological human intestinal model. However, the human enteroid model enables specific functional studies of secretion and absorption in each intestinal segment as well as observations of the earliest molecular events that occur during enteric infections. This Review describes studies characterizing these human mini-guts as a physiological model to investigate intestinal transport and host-pathogen interactions., Competing Interests: statement The authors declare no competing interests.

A two dimensional numerical model is used to study the muco-ciliary transport process in human respiratory tract. Here, hybrid finite difference-lattice Boltzmann method is used to model the flow physics of the transport of mucus and periciliary liquid (PCL) layer in the airway surface liquid. The immersed boundary method is also used to implement the propulsive effect of the cilia and also the effects of the interface between the mucus and PCL layers. The main contribution of this study is on elucidating the role of the viscoelastic behavior of mucus on the muco-ciliary transport and for this purpose an Oldroyd-B model is used as the constitutive equation of mucus for the first time. Results show that the viscosity and viscosity ratio of mucus have an enormous effect on the muco-ciliary transport process. It is also seen that the mucus velocity is affected by mucus relaxation time when its value is less than 0.002 s. Results also indicate that the variation of these properties on the mucus velocity at lower values of viscosity ratio is more significant., (Copyright © 2015. Published by Elsevier Inc.)

One-dimensional models coupling a Langevin equation for the cargo position to stochastic stepping dynamics for the motors constitute a relevant framework for analyzing multiple-motor microtubule transport. In this work we explore the consistence of these models focusing on the effects of the thermal noise. We study how to define consistent stepping and detachment rates for the motors as functions of the local forces acting on them in such a way that the cargo velocity and run-time match previously specified functions of the external load, which are set on the base of experimental results. We show that due to the influence of the thermal fluctuations this is not a trivial problem, even for the single-motor case. As a solution, we propose a motor stepping dynamics which considers memory on the motor force. This model leads to better results for single-motor transport than the approaches previously considered in the literature. Moreover, it gives a much better prediction for the stall force of the two-motor case, highly compatible with the experimental findings. We also analyze the fast fluctuations of the cargo position and the influence of the viscosity, comparing the proposed model to the standard one, and we show how the differences on the single-motor dynamics propagate to the multiple motor situations. Finally, we find that the one-dimensional character of the models impede an appropriate description of the fast fluctuations of the cargo position at small loads. We show how this problem can be solved by considering two-dimensional models.

We present and explore a kinetic model of ion transport across and between the membranes of an isolated plant cell with an emphasis on the cell's response to salt (Na(+)) stress. The vacuole, cytoplasm and apoplast are treated as concentric regions separated by tonoplast and plasma membranes. The model includes the transport of Na(+), K(+), Cl(-) and H(+) across both membranes via primary active proton pumps, secondary active antiporters and symporters, as well as passive ion channels. In addition, water transport is included, allowing us to investigate both the osmotic and ionic components of salt stress. The model's predictions of steady state and transient cytosolic pH and Na(+) concentrations were found to be quantitatively comparable to measured experimental values. Through an extensive simulation study we have identified and characterized scenarios in which individual transport processes (H(+) pumps, Na(+)/H(+) antiporters and channels involved in the transport of Na(+)) and their combinations have major effects on the level of Na(+) in each of the cell compartments. This systematic study emulates the effects of overexpressing and inhibiting transporter genes by genetic modification and hence we have compared our simulations with observations from experiments conducted on transgenic plants. The simulations suggest that overexpressing tonoplast Na(+)/H(+) antiporter genes and tonoplast H(+) pump genes lead to an increase in the storage of Na(+) in the vacuole (helping the cell to maintain water uptake under salt stress), with only a transient influence on the cytoplasmic Na(+) concentration. The model predicts effects of varying the expression of transporter genes (individually or in combination) which have yet to be investigated in experiments. For example, our findings indicate that simultaneously overexpressing plasma membrane and tonoplast Na(+)/H(+) antiporter genes would lead to improvements in both ionic and osmotic stress tolerance. The results demonstrate the importance of simultaneously modelling the transport of Na(+) across both the tonoplast and plasma membrane, a task not undertaken previously., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

The past decade has seen increased numbers of studies publishing ligand-based computational models for drug transporters. Although they generally use small experimental data sets, these models can provide insights into structure-activity relationships for the transporter. In addition, such models have helped to identify new compounds as substrates or inhibitors of transporters of interest. We recently proposed that many transporters are promiscuous and may require profiling of new chemical entities against multiple substrates for a specific transporter. Furthermore, it should be noted that virtually all of the published ligand-based transporter models are only accessible to those involved in creating them and, consequently, are rarely shared effectively. One way to surmount this is to make models shareable or more accessible. The development of mobile apps that can access such models is highlighted here. These apps can be used to predict ligand interactions with transporters using Bayesian algorithms. We used recently published transporter data sets (MATE1, MATE2K, OCT2, OCTN2, ASBT, and NTCP) to build preliminary models in a commercial tool and in open software that can deliver the model in a mobile app. In addition, several transporter data sets extracted from the ChEMBL database were used to illustrate how such public data and models can be shared. Predicting drug-drug interactions for various transporters using computational models is potentially within reach of anyone with an iPhone or iPad. Such tools could help prioritize which substrates should be used for in vivo drug-drug interaction testing and enable open sharing of models., (Copyright © 2015 by The American Society for Pharmacology and Experimental Therapeutics.)

Viruses and other cellular cargo that lack locomotion must rely on diffusion and cellular transport systems to navigate through a biological cell. Indeed, advances in single particle tracking have revealed that viral motion alternates between (a) diffusion in the cytoplasm and (b) active transport along microtubules. This intermittency makes quantitative analysis of trajectories difficult. Therefore, the purpose of this paper is to construct mathematical methods to approximate intermittent dynamics by effective stochastic differential equations. The coarse-graining method that we develop is more accurate than existing techniques and applicable to a wide range of intermittent transport models. In particular, we apply our method to two- and three-dimensional cell geometries (disk, sphere, and cylinder) and demonstrate its accuracy. In addition to these specific applications, we also explain our method in full generality for use on future intermittent models.

Man-made transport networks and their design are closely related to the shortest path problem and considered amongst the most debated problems of computational intelligence. Apart from using conventional or bio-inspired computer algorithms, many researchers tried to solve this kind of problem using biological computing substrates, gas-discharge solvers, prototypes of a mobile droplet, and hot ice computers. In this aspect, another example of biological computer is the plasmodium of acellular slime mould Physarum polycephalum (P. polycephalum), which is a large single cell visible by an unaided eye and has been proven as a reliable living substrate for implementing biological computing devices for computational geometry, graph-theoretical problems, and optimization and imitation of transport networks. Although P. polycephalum is easy to experiment with, computing devices built with the living slime mould are extremely slow; it takes slime mould days to execute a computation. Consequently, mapping key computing mechanisms of the slime mould onto silicon would allow us to produce efficient bio-inspired computing devices to tackle with hard to solve computational intelligence problems like the aforementioned. Toward this direction, a cellular automaton (CA)-based, Physarum-inspired, network designing model is proposed. This novel CA-based model is inspired by the propagating strategy, the formation of tubular networks, and the computing abilities of the plasmodium of P. polycephalum. The results delivered by the CA model demonstrate a good match with several previously published results of experimental laboratory studies on imitation of man-made transport networks with P. polycephalum. Consequently, the proposed CA model can be used as a virtual, easy-to-access, and biomimicking laboratory emulator that will economize large time periods needed for biological experiments while producing networks almost identical to the tubular networks of the real-slime mould.

Microtubule-based distribution of organelles/vesicles is crucial for the function of many types of eukaryotic cells and the molecular motor cytoplasmic dynein is required for transporting a variety of cellular cargos toward the microtubule minus ends. Early endosomes represent a major cargo of dynein in filamentous fungi, and dynein regulators such as LIS1 and the dynactin complex are both required for early endosome movement. In fungal hyphae, kinesin-3 and dynein drive bi-directional movements of early endosomes. Dynein accumulates at microtubule plus ends; this accumulation depends on kinesin-1 and dynactin, and it is important for early endosome movements towards the microtubule minus ends. The physical interaction between dynein and early endosome requires the dynactin complex, and in particular, its p25 component. The FTS-Hook-FHIP (FHF) complex links dynein-dynactin to early endosomes, and within the FHF complex, Hook interacts with dynein-dynactin, and Hook-early endosome interaction depends on FHIP and FTS.