Model simulation of heat and water transport dynamics in an airway

Heat and water transport processes in the respiratory tract depend on environmental conditions, breathing patterns, and the physiological state of the respiratory system. To study these processes, we have developed a mathematical model of the dynamics of temperature and water vapor in the radial and axial directions of an idealized trachea. The model is expressed as two implicit finite-difference equations and solved using an alternating-direction algorithm. Using these equations, we simulated the effects of inspired gas temperature and humidity, velocity profile, and flow rate on heat and water transport between the gas and airway wall. Under inspired gas conditions of low temperature or high relative humidity, supersaturation occurs. Increasing either the velocity gradient at the wall or the flow rate increases the heat and water transport rates. However, these rates change by only 10 percent when the velocity gradient is doubled, and by about 35 percent when flow rate undergoes a two-fold change. The model can be used with in-vivo data from the trachea to test hypotheses concerning normal and abnormal heat and water transport.     

Extracellular noise-induced stochastic synchronization in heterogeneous quorum sensing network

Quorum sensing is a bacterial mechanism used to synchronize the coordinated response of a microbial population. Because quorum sensing in Gram-negative bacteria depends on release and detection of a diffusible signaling molecule (autoinducer) among a multicellular group, it is considered a simple form of cell-cell communication for the purposes of mathematical analysis. Stochastic equation systems have provided a common approach to model biochemical or biophysical processes. Recently, the effect of noise to synchronize a specific homogeneous quorum sensing network was successfully modeled using a stochastic equation system with fixed parameters. The question remains of how to model quorum sensing networks in a general setting. To address this question, we first set a stochastic equation system as a general model for a heterogeneous quorum sensing network. Then, using two relevant biophysical characteristics of Gram-negative bacteria (the permeability of the cell membrane to the autoinducer and the symmetry of autoinducer diffusion) we construct the solution of the stochastic equation system at an abstract level. The solution indicates that stable synchronization of a quorum sensing network is robustly induced by an environment with a heterogenous distribution of extracellular and intracellular noise. The synchronization is independent of the initial state of the system and is solely the result of the connectivity of the cell network established through the effects of extracellular noise.     

Discriminative responses of squid (Loligo pealeii) photoreceptors to polarized light

Cephalopods behaviorally respond to polarized light. Electrophysiology experiments with the squid, Loligo pealeii, demonstrated that spike responses from individual photoreceptors are a cosine2 function of the e-vector orientation of a polarized stimulus. The discrimination limit to this polarization sensitivity depended upon the difference between the orientation of a polarized stimulus with a preferred e-vector. The limit ranged from 2 degrees to 9.2 degrees with a direct stimulus in the dark or 4.8 degrees -22.1 degrees with non-directed background illumination and the cells were least discriminative at the preferred orientations. This limit can be explained partly by the variability in anatomical alignment of microvilli in the photoreceptors around a dominant axis. A few light-sensitive retinal fibers showed no polarization sensitivity. The coding of polarization information suggests that light intensity is transformed into an average spike rate. This average results from silent periods interspersed between bursts of spikes, each burst possessing a consistent interspike interval. The variations in the length and frequency of silent periods depend upon the difference between the polarization e-vector and a preferred e-vector orientation. The minimal discriminated orientation of a squid photoreceptor agrees well with the minimum behavioral discrimination of polarized light by another cephalopod, the octopus.     

Characterization of the Aerial Escape Response of the African Butterfly Fish, Pantodon buchholzi Peters

We studied the startle response of the African butterfly fish, Pantodon buchholzi (Osteoglossomorpha, Osteoglossoidea). It is an upward movement, mediated by abduction of the pectoral fins, and is elicited by mechanical and visual stimuli. Because this fish inhabits the first few centimeters beneath the water surface, its startle response results in an aerial excursion that may be described as ballistic-like, following a motion as defined by linear acceleration. We show that the aerial excursion is well-modeled by a parabola. On average, a fish jumps no more than twice its height and travels horizontally about five times its standard length. The fish may exhibit variable in-flight trunk and fin movements, but neither increases the travel distance in air following the initial in-water propulsive event. Similar vertical jumps also occur entirely within the water column suggesting that this motor behavior of Pantodon is a general escape behavior analogous to a Mauthner neuron-induced escape response. The variability in its posture in air and its direction of motion after reentering the water enhances this act of vertical flight as a step in this fish's escape behavior. The aerial aspect of its escape behavior is only a consequence of its position in the water column.     

An analysis of countercurrent exchange with emphasis on renal function

An analysis of countercurrent exchange in a U-tube is presented for a single-solute, constant-volume flow rate system with spatially varying source fluxes and permeabilities. Analytical solutions are given for the steady-state equations and numerical solutions for the unsteady-state equations. The solutions indicate that an external source of solute delivered to the stream flowing away from the U-tube bend can be distributed by the exchanger so that the concentration in both limbs increases toward the bend. In particular, there exist source fluxes whose magnitude decreases monotonically toward the bend for which the maximum solute concentration occurs at the bend. The point at which a concentration maximum occurs is governed principally by the solute permeability of the barrier separating the two limbs and by the volume flow rate through the exchanger. The system dynamics depend strongly on the relative cross-sectional areas of the two limbs or, equivalently, on the flow velocities within them. The model is used as a basis for discussion of various functional aspects of the renal vasa recta system.     

Countercurrent exchange in the inner renal medulla: Vasa recta-descending limb system

A model of countercurrent exchange has been developed to simulate transport of salt, urea and water among vasa recta and descending limbs of the loop of Henle in the inner medulla. These vessels are abstracted as three concentric cylinders: the inne one represents descending vasa recta, the middle one represents ascending vasa recta and the outer one represents descending limbs. The capillary plexus, which connects the ascending and descending vasa recta, is modeled as a series of well-mixe compartments. Multicomponent transport equations for the sytem are derived from steady state mass balances and simple passive flux relations. The resulting set of nonlinear equations are solved numerically by an iterative Gauss-Seidel algorithm with under-relaxation. Simulations yield the salt and urea concentrations as well as volume flow rates in all tubes and compartments. The simulations indicate that solute concentrations can increase monotonically toward the papillae even if all transport processes within the exchanger are passive and source fluxes decrease monotonically toward the papillae.     

System for dynamic measurements of membrane capacitance in intact epithelial monolayers.

Dynamic measurements of exocytosis have been difficult to perform in intact epithelial monolayers. We have designed a system that estimates with +/-1% accuracy (99% confidence) the total membrane capacitance of monolayers represented by a lumped model. This impedance measurement and analysis system operates through a conventional transepithelial electrophysiology clamp, performing all signal measurements as frequently as every 5 s. Total membrane capacitance (the series combination of apical and basolateral membranes) is the inverse of one of three unique coefficients that describe the monolayer impedance. These coefficients are estimated using a weighted, nonlinear, least-squares algorithm. Using the estimated coefficients, solution ranges for individual membrane parameters are calculated, frequently providing results within +/-20% of true values without additional electrophysiological measurements. We determined the measurement system specifications and statistical significance of estimated parameters using 1) analytical testing with circuit simulation software and equation-generated data; 2) a system noise analysis combined with Monte Carlo simulations; and 3) analog model circuits for calibration of the electronic system and to check equation-generated results. Finally, the time course of capacitance changes associated with purinergically stimulated mucin exocytosis are quantified in monolayers of the colonic goblet cell-like cell line HT29-CI.16E.     

Novel anti-prostate cancer scaffold identified by the combination of in silico and cell-based assays targeting the PI3K-AKT-mTOR pathway

Phenotypic assays were performed in prostate cancer cell lines to describe the biological activity of PI3K-AKT-mTOR pathway inhibitors retrieved from the virtual screening initiative. These novel chemicals share in common the aminopyridine scaffold, hitting PC-3 cells in macromolar range, with selectivity index over fibroblast cell lines. Moreover, a preliminary study of the mode of action by flow cytometry assay pointed out that these compounds had a rapamycin-like response for the PI3K-AKT-mTOR pathway modulation.