Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, Biological”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

A universal steady state I-V relationship for membrane current

A purely electrical mechanism for the gating of membrane ionic channel gives rise to a simple I-V relationship for membrane current. Our approach is based on the known presence of gating charge, which is an established property of the membrane channel gating. The gating charge is systematically treated as a polarization of the channel protein which varies with the external electric field and modifies the effective potential through which the ions migrate in the channel. Two polarization effects have been considered: 1) the up or down shift of the whole potential function, and 2) the change in the effective electric field inside the channel which is due to familiar effect of the effective reduction of the electric field inside a dielectric body because of the presence of surface charges on its surface. Both effects are linear in the channel polarization. The ionic current is described by a steady state solution of the Nernst-Planck equation with the potential directly controlled by the gating charge system. The solution describes reasonably well the steady state and peak-current I-V relationships for different channels, and when applied adiabatically, explains the time lag between the gating charge current and the rise of the ionic current. The approach developed can be useful as an effective way to model the ionic currents in axons, cardiac cells and other excitable tissues.

NASA Discipline Regulatory Physiology↗

A case for bone canaliculi as the anatomical site of strain generated potentials

We address the question of determining the anatomical site that is the source of the experimentally observed strain generated potentials (SGPs) in bone tissue. There are two candidates for the anatomical site that is the SGP source, the collagen-hydroxyapatite porosity and the larger size lacunar-canalicular porosity. In the past it has been argued, on the basis of experimental data and a reasonable model, that the site of the SGPs in bone is the collagen-hydroxyapatite porosity. The theoretically predicted pore radius necessary for the SGPs to reside in this porosity is 16 nm, which is somewhat larger than the pore radii estimated from gas adsorption data where the preponderance of the pores were estimated to be in the range 5-12.5 nm. However, this pore size is significantly larger than the 2 nm size of the small tracer, microperoxidase, which appears to be excluded from the mineralized matrix. In this work a similar model, but one in which the effects of fluid dynamic drag of the cell surface matrix in the bone canaliculi are included, is used to show that it is possible for the generation of SGPs to be associated with the larger size lacunar-canalicular porosity when the hydraulic drag and electrokinetic contribution of the bone fluid passage through the cell coat (glycocalyx) is considered. The consistency of the SGP data with this model is demonstrated. A general boundary condition is introduced to allow for current leakage at the bone surface. The results suggest that the current leakage is small for the in vitro studies in which the strain generated potentials have been measured.

Non-NASA Center↗

An equivalent body surface charge model representing three-dimensional bioelectrical activity

A new surface-source model has been developed to account for the bioelectrical potential on the body surface. A single-layer surface-charge model on the body surface has been developed to equivalently represent bioelectrical sources inside the body. The boundary conditions on the body surface are discussed in relation to the surface-charge in a half-space conductive medium. The equivalent body surface-charge is shown to be proportional to the normal component of the electric field on the body surface just outside the body. The spatial resolution of the equivalent surface-charge distribution appears intermediate between those of the body surface potential distribution and the body surface Laplacian distribution. An analytic relationship between the equivalent surface-charge and the surface Laplacian of the potential was found for a half-space conductive medium. The effects of finite spatial sampling and noise on the reconstruction of the equivalent surface-charge were evaluated by computer simulations. It was found through computer simulations that the reconstruction of the equivalent body surface-charge from the body surface Laplacian distribution is very stable against noise and finite spatial sampling. The present results suggest that the equivalent body surface-charge model may provide an additional insight to our understanding of bioelectric phenomena.

NASA Discipline Regulatory Physiology↗

A fiber matrix model for fluid flow and streaming potentials in the canaliculi of an osteon

A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a magnitude that lies within the middle of the range, 6-30 dynes/cm2, where fluid shear stresses in tissue culture studies with osteoblast monolayers have elicited an intracellular Ca++ response (31). The implications of these results for intracellular electrical communication are discussed.

Non-NASA Center↗

Physiological time-series analysis: what does regularity quantify?

Approximate entropy (ApEn) is a recently developed statistic quantifying regularity and complexity that appears to have potential application to a wide variety of physiological and clinical time-series data. The focus here is to provide a better understanding of ApEn to facilitate its proper utilization, application, and interpretation. After giving the formal mathematical description of ApEn, we provide a multistep description of the algorithm as applied to two contrasting clinical heart rate data sets. We discuss algorithm implementation and interpretation and introduce a general mathematical hypothesis of the dynamics of a wide class of diseases, indicating the utility of ApEn to test this hypothesis. We indicate the relationship of ApEn to variability measures, the Fourier spectrum, and algorithms motivated by study of chaotic dynamics. We discuss further mathematical properties of ApEn, including the choice of input parameters, statistical issues, and modeling considerations, and we conclude with a section on caveats to ensure correct ApEn utilization.

NASA Discipline Cardiopulmonary↗

Overview of spaceflight immunology studies

The effects of spaceflight and analogues of spaceflight are discussed here and in nine accompanying articles. In this summary we present spaceflight studies with human subjects, animal subjects, and cell cultures and we review ground-based systems used to model the observed effects of spaceflight on the immune system. Human paradigms include bed rest, academic or psychological stress, physical stress, hypobaric or high altitude stress, and confinement. Animal models include antiorthostatic and orthostatic suspension, hypobarism, and confinement. The ten manuscripts in this collection were selected to provide a summary that should give the reader an overview of the various activities of spaceflight immunology researchers throughout the history of space travel. This manuscript identifies the major contributors to the study of spaceflight immunology, explains what types of studies have been conducted, and how they have changed over the years. Also presented is a discussion of the unusual limitations associated with spaceflight research and the efforts to develop appropriate ground-based surrogate model systems. Specific details, data, and mechanistic speculations will be held to a minimum, because they will be discussed in depth in the other articles in the collection.

NASA Discipline Number 18-10↗

Influence of simulated microgravity on the VO2 max of nontrained and trained rats

Head-down suspension (HDS) of rats has evolved as a useful model for the simulation of a microgravity environment. Previous HDS experiments with rats have shown an impaired capacity to perform aerobic exercise as demonstrated by reductions in maximum oxygen consumption (VO2 max), treadmill run time (RT), and mechanical efficiency (ME) of treadmill running at submaximal conditions. To determine whether endurance training (TR) before HDS would modify exercise performance, male Sprague-Dawley rats were assigned to nontrained (NT) or TR groups for 6 wk and exposed to HDS or cage control (CC) conditions for 29 days. The rats were tested for VO2 max, RT, and ME before treatment and on days 7, 14, 21, and 28. In addition, water and electrolyte excretion was measured on days 1 and 21 of the experimental period. Before HDS, the TR rats had significantly higher measures of VO2 max (15%) and RT (22%) than the NT rats. On day 28, HDS was associated with significant reductions in absolute VO2 max (ml/min) in TR (-30%) and NT (-14%) rats. Relative VO2 max (ml.min-1.kg-1) was significantly reduced in TR (-15%) but not NT rats. Similar reductions in RT occurred in TR (-37%) and NT (-35%) rats by day 28. ME was reduced 22% in both TR and NT rats after 28 days of suspension. HDS elicited diuresis, natriuresis, and kaliuresis in TR rats after 21 days but not after 24 h. In contrast, HDS-NT rats exhibited no diuretic, natriuretic, or kaliuretic responses.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

Decomposition of superimposed ground reaction forces into left and right force profiles

The process of collecting ground reaction force data by mounting a forceplate beneath a treadmill belt has the advantage that numerous walking trials can be analyzed without the problem of subjects 'targeting' their footsteps. However, a potential problem is that the measured forces represent a summation of bilateral force profiles during the double support phase of walking. To address this issue, an algorithm is described for decomposing superimposed ground reaction force data into individual left and right profiles. It is based on an examination of the side-to-side oscillations of the measured center of pressure (CoP). Whenever the measured CoP exceeds a certain threshold, it is assumed that the person is being supported by a single limb, and the measured GRF data reflect the forces under that limb. Conversely, when the measured CoP indicates that both feet are on the treadmill, it is assumed that the location of the individual CoP under each foot is given by wL2 and wR2. These quantities reflect the greatest excursion of the measured CoP towards the left and right sides of the forceplate, respectively. With this assumption, individual GRF profiles can be calculated by means of solving two simultaneous equations--one describing the equilibrium of forces in the vertical direction, and one describing the equilibrium of moments about an antero-posterior axis of the forceplate. The algorithm describing this procedure is simple enough to be implemented on a spreadsheet and yields estimates for average force, impulse, peak force and stance time that are typically within 3% of the true values.

Non-NASA Center↗

Alveolar ventilation to perfusion heterogeneity and diffusion impairment in a mathematical model of gas exchange

This study describes a two-compartment model of pulmonary gas exchange in which alveolar ventilation to perfusion (VA/Q) heterogeneity and impairment of pulmonary diffusing capacity (D) are simultaneously taken into account. The mathematical model uses as input data measurements usually obtained in the lung function laboratory. It consists of two compartments and an anatomical shunt. Each compartment receives fractions of alveolar ventilation and blood flow. Mass balance equations and integration of Fick's law of diffusion are used to compute alveolar and blood O2 and CO2 values compatible with input O2 uptake and CO2 elimination. Two applications are presented. The first is a method to partition O2 and CO2 alveolar-arterial gradients into VA/Q and D components. The technique is evaluated in data of patients with chronic obstructive pulmonary disease (COPD). The second is a theoretical analysis of the effects of blood flow variation in alveolar and blood O2 partial pressures. The results show the importance of simultaneous consideration of D to estimate VA/Q heterogeneity in patients with diffusion impairment. This factor plays an increasing role in gas alveolar-arterial gradients as severity of COPD increases. Association of VA/Q heterogeneity and D may produce an increase of O2 arterial pressure with decreasing QT which would not be observed if only D were considered. We conclude that the presented computer model is a useful tool for description and interpretation of data from COPD patients and for performing theoretical analysis of variables involved in the gas exchange process.

NASA Discipline Environmental Health↗

Mass balance approaches for estimating the intestinal absorption and metabolism of peptides and analogues: theoretical development and applications

A theoretical analysis for estimating the extent of intestinal peptide and peptide analogue absorption was developed on the basis of a mass balance approach that incorporates convection, permeability, and reaction. The macroscopic mass balance analysis (MMBA) was extended to include chemical and enzymatic degradation. A microscopic mass balance analysis, a numerical approach, was also developed and the results compared to the MMBA. The mass balance equations for the fraction of a drug absorbed and reacted in the tube were derived from the general steady state mass balance in a tube: [formula: see text] where M is mass, z is the length of the tube, R is the tube radius, Pw is the intestinal wall permeability, kr is the reaction rate constant, C is the concentration of drug in the volume element over which the mass balance is taken, VL is the volume of the tube, and vz is the axial velocity of drug. The theory was first applied to the oral absorption of two tripeptide analogues, cefaclor (CCL) and cefatrizine (CZN), which degrade and dimerize in the intestine. Simulations using the mass balance equations, the experimental absorption parameters, and the literature stability rate constants yielded a mean estimated extent of CCL (250-mg dose) and CZN (1000-mg dose) absorption of 89 and 51%, respectively, which was similar to the mean extent of absorption reported in humans (90 and 50%). It was proposed previously that 15% of the CCL dose spontaneously degraded systematically; however, our simulations suggest that significant CCL degradation occurs (8 to 17%) presystemically in the intestinal lumen.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

A multidimensional model of the effect of gravity on the spatial orientation of the monkey

A "sensory conflict" model of spatial orientation was developed. This mathematical model was based on concepts derived from observer theory, optimal observer theory, and the mathematical properties of coordinate rotations. The primary hypothesis is that the central nervous system of the squirrel monkey incorporates information about body dynamics and sensory dynamics to develop an internal model. The output of this central model (expected sensory afference) is compared to the actual sensory afference, with the difference defined as "sensory conflict." The sensory conflict information is, in turn, used to drive central estimates of angular velocity ("velocity storage"), gravity ("gravity storage"), and linear acceleration ("acceleration storage") toward more accurate values. The model successfully predicts "velocity storage" during rotation about an earth-vertical axis. The model also successfully predicts that the time constant of the horizontal vestibulo-ocular reflex is reduced and that the axis of eye rotation shifts toward alignment with gravity following postrotatory tilt. Finally, the model predicts the bias, modulation, and decay components that have been observed during off-vertical axis rotations (OVAR).

NASA Discipline Neuroscience↗

Experimental tests of a superposition hypothesis to explain the relationship between the vestibuloocular reflex and smooth pursuit during horizontal combined eye-head tracking in humans

1. We used a modeling approach to test the hypothesis that, in humans, the smooth pursuit (SP) system provides the primary signal for cancelling the vestibuloocular reflex (VOR) during combined eye-head tracking (CEHT) of a target moving smoothly in the horizontal plane. Separate models for SP and the VOR were developed. The optimal values of parameters of the two models were calculated using measured responses of four subjects to trials of SP and the visually enhanced VOR. After optimal parameter values were specified, each model generated waveforms that accurately reflected the subjects' responses to SP and vestibular stimuli. The models were then combined into a CEHT model wherein the final eye movement command signal was generated as the linear summation of the signals from the SP and VOR pathways. 2. The SP-VOR superposition hypothesis was tested using two types of CEHT stimuli, both of which involved passive rotation of subjects in a vestibular chair. The first stimulus consisted of a "chair brake" or sudden stop of the subject's head during CEHT; the visual target continued to move. The second stimulus consisted of a sudden change from the visually enhanced VOR to CEHT ("delayed target onset" paradigm); as the vestibular chair rotated past the angular position of the stationary visual stimulus, the latter started to move in synchrony with the chair. Data collected during experiments that employed these stimuli were compared quantitatively with predictions made by the CEHT model. 3. During CEHT, when the chair was suddenly and unexpectedly stopped, the eye promptly began to move in the orbit to track the moving target. Initially, gaze velocity did not completely match target velocity, however; this finally occurred approximately 100 ms after the brake onset. The model did predict the prompt onset of eye-in-orbit motion after the brake, but it did not predict that gaze velocity would initially be only approximately 70% of target velocity. One possible explanation for this discrepancy is that VOR gain can be dynamically modulated and, during sustained CEHT, it may assume a lower value. Consequently, during CEHT, a smaller-amplitude SP signal would be needed to cancel the lower-gain VOR. This reduction of the SP signal could account for the attenuated tracking response observed immediately after the brake. We found evidence for the dynamic modulation of VOR gain by noting differences in responses to the onset and offset of head rotation in trials of the visually enhanced VOR.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

A new animal model for modulating myosin isoform expression by altered mechanical activity

The purpose of this study was to develop a new rodent model that is capable of delineating the importance of mechanical loading on myosin heavy chain (MHC) isoform expression of the plantar and dorsi flexor muscles of the ankle. The essential components of this system include 1) stimulating electrodes that are chronically implanted into a muscle, allowing for the control of the activation pattern of the target muscle(s); 2) a training apparatus that translates the moment of the ankle into a linear force; and 3) a computer-controlled Cambridge 310 ergometer. The isovelocity profile of the ergometer ensured that the medial gastrocnemius (MG) produced forces that were > 90% of maximal isometric force (Po), and the eccentric contractions of the tibialis anterior (TA) were typically 120% of Po. Both the concentric and eccentric training programs produced statistically significant increases in the muscle mass of the MG (approximately 15%) and TA (approximately 7%) as well as a decrease in myofibrillar adenosinetriphosphatase activity. Both the white and red regions of the MG and TA exhibited significant increases in the relative content of the type IIa MHC and concomitant decreases in type IIb MHC expression. Although the red regions of the MG and red TA contained approximately 10% type I MHC, the training programs did not affect this isoform. It appears that when a fast-twitch muscle is stimulated at a high frequency (100 Hz) and required to contract either concentrically or eccentrically under high loading conditions, the expression of the type IIa MHC isoform will be upregulated, whereas that of the type IIb MHC will be concomitantly downregulated.

NASA Discipline Musculoskeletal↗

Modeling the effects of hypoxia on ATP turnover in exercising muscle

Most models of metabolic control concentrate on the regulation of ATP production and largely ignore the regulation of ATP demand. We describe a model, based on the results of Hogan et al. (J. Appl. Physiol. 73: 728-736, 1992), that incorporates the effects of ATP demand. The model is developed from the premise that a unique set of intracellular conditions can be measured at each level of ATP turnover and that this relationship is best described by energetic state. Current concepts suggest that cells are capable of maintaining oxygen consumption in the face of declines in the concentration of oxygen through compensatory changes in cellular metabolites. We show that these compensatory changes can cause significant declines in ATP demand and result in a decline in oxygen consumption and ATP turnover. Furthermore we find that hypoxia does not directly affect the rate of anaerobic ATP synthesis and associated lactate production. Rather, lactate production appears to be related to energetic state, whatever the PO2. The model is used to describe the interaction between ATP demand and ATP supply in determining final ATP turnover.

Non-NASA Center↗

Muscle fatigue in frog semitendinosus: alterations in contractile function

The purpose of this study was to characterize the contractile properties of the frog semitendinosus (ST) muscle before and during recovery from fatigue, to relate the observed functional changes to alterations in specific steps in the crossbridge model of muscle contraction, and to determine how fatigue affects the force-frequency relationship. The frog ST (22 degrees C) was fatigued by direct electrical stimulation with 100-ms 150-Hz trains at 1/s for 5 min. The fatigue protocol reduced peak twitch (Pt) and tetanic (Po) force to 32 and 8.5% of initial force, respectively. The decline in Pt was less than Po, in part due to a prolongation in the isometric contraction time (CT), which increased to 300% of the initial value. The isometric twitch duration was greatly prolonged as reflected by the lengthened CT and the 800% increase in the one-half relaxation time (1/2RT). Both Pt and Po showed a biphasic recovery, a rapid initial phase (2 min) followed by a slower (40 min) return to the prefatigue force. CT and 1/2RT also recovered in two phases, returning to 160 and 265% of control in the first 5 min. CT returned to the prefatigue value between 35 and 40 min, whereas even at 60 min 1/2RT was 133% of control. The maximal velocity of shortening, determined by the slack test, was significantly reduced [from 6.7 +/- 0.5 to 2.5 +/- 0.4 optimal muscle length/s] at fatigue. The force-frequency relationship was shifted to the left, so that optimal frequency for generating Po was reduced.(ABSTRACT TRUNCATED AT 250 WORDS).

NASA Discipline Musculoskeletal↗

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

Mathematical modelling of flow distribution in the human cardiovascular system

The paper presents a detailed model of the entire human cardiovascular system which aims to study the changes in flow distribution caused by external stimuli, changes in internal parameters, or other factors. The arterial-venous network is represented by 325 interconnected elastic segments. The mathematical description of each segment is based on equations of hydrodynamics and those of stress/strain relationships in elastic materials. Appropriate input functions provide for the pumping of blood by the heart through the system. The analysis employs the finite-element technique which can accommodate any prescribed boundary conditions. Values of model parameters are from available data on physical and rheological properties of blood and blood vessels. As a representative example, simulation results on changes in flow distribution with changes in the elastic properties of blood vessels are discussed. They indicate that the errors in the calculated overall flow rates are not significant even in the extreme case of arteries and veins behaving as rigid tubes.

NASA Center JSC↗