A mechanistic model of the actin cycle.

We have derived a broad, deterministic model of the steady-state actin cycle that includes its major regulatory mechanisms. Ours is the first model to solve the complete nucleotide profile within filaments, a feature that determines the dynamics and geometry of actin networks at the leading edges of motile cells, and one that has challenged investigators developing models to interpret steady-state experiments. We arrived at the nucleotide profile through analytic and numerical approaches that completely agree. Our model reproduces behaviors seen in numerous experiments with purified proteins, but allows a detailed inspection of the concentrations and fluxes that might exist in these experiments. These inspections provide new insight into the mechanisms that determine the rate of actin filament treadmilling. Specifically, we find that mechanisms for enhancing Pi release from the ADP.Pi intermediate on filaments, for increasing the off rate of ADP-bound subunits at pointed ends, and the multiple, simultaneous functions of profilin, make unique and essential contributions to increased treadmilling. In combination, these mechanisms have a theoretical capacity to increase treadmilling to levels limited only by the amount of available actin. This limitation arises because as the cycle becomes more dynamic, it tends toward the unpolymerized state.

Regulation of the actin cycle in vivo by actin filament severing.

Cycling of actin subunits between monomeric and filamentous phases is essential for cell crawling behavior. We investigated actin filament turnover rates, length, number, barbed end exposure, and binding of cofilin in bovine arterial endothelial cells moving at different speeds depending on their position in a confluent monolayer. Fast-translocating cells near the wound edge have short filament lifetimes compared with turnover values that proportionately increase in slower moving cells situated at increasing distances from the wound border. Contrasted with slow cells exhibiting slow actin filament turnover speeds, fast cells have less polymerized actin, shorter actin filaments, more free barbed ends, and less actin-associated cofilin. Cultured primary fibroblasts manifest identical relationships between speed and actin turnover as the endothelial cells, and fast fibroblasts expressing gelsolin have higher actin turnover rates than slow fibroblasts that lack this actin-severing protein. These results implicate actin filament severing as an important control mechanism for actin cycling in cells.

The human physiome as an information environment.

Modern biology is rapidly laying the foundation necessary for integrated modeling of physiological processes in living organisms. The human physiome project attempts to model interactions between biochemicals, cellular organelles, cells, tissues, and organs within whole organisms. One of the first challenges that this project faces is the development of a database environment flexible enough to accommodate the diversity in structure and content of physiological data. This paper reviews the current state of database technology, presents our understanding of the physiome database problem, and proposes a preliminary strategy for addressing it.

Three-dimensional reconstruction of the actin cytoskeleton from stereo images.

In eucaryotic cells, actin filaments are abundant components in the cytoskeleton where they form a complex three dimensional (3D) structural network that provides the cell with its shape and mechanical properties. However, understanding the structural and mechanical properties of actin filaments composing the cell cytoskeleton is often hampered by the inability to faithfully reconstruct the three-dimensional geometric relationships. This paper presents a vision-based reconstruction approach that automatically reconstitutes the three-dimensional structures of cytoskeletal polymers from stereo image pairs taken at the different tilt angles. The approach finds corresponding points between two images and recovers the depth information about the structures. The computational process consists of three major procedures: feature representation, stereo matching, and disparity refinement, implemented in a multi-resolution manner based on a coarse-to-fine strategy. The reconstruction depicts the three-dimensional structure of cytoskeletal polymers and their geometric relationships. New and useful information becomes available and allows quantitative analysis of the structure. Measurement of the cytoskeleton geometrical properties and the filament concentration in a defined volume are obtained by direct calculation.

Vascular endothelial cells respond to spatial gradients in fluid shear stress by enhanced activation of transcription factors.

The vascular endothelium is exposed to a spectrum of fluid mechanical forces generated by blood flow; some of these, such as fluid shear stress, can directly modulate endothelial gene expression. Previous work by others and in our laboratory, using an in vitro uniform laminar shear stress model, has identified various shear stress response elements (SSREs) within the promoters of certain endothelial genes that regulate their expression by interacting with various transcription factors, including nuclear factor-kappaB (NF-kappaB), early growth response-1 (Egr-1), and activator protein-1 (AP-1, composed of c-Jun/c-Jun and c-Jun/c-Fos protein dimers). In the current study, we have examined the topographical patterns of NF-kappaB, Egr-1, c-Jun, and c-Fos activation in a specially designed in vitro disturbed laminar shear stress model, which incorporates regions of significant spatial shear stress gradients similar to those found in atherosclerosis-prone arterial geometries in vivo (eg, arterial bifurcations, curvatures, ostial openings). Using newly developed quantitative image analysis techniques, we demonstrate that endothelial cells subjected to disturbed laminar shear stress exhibit increased levels of nuclear localized NF-kappaB, Egr-1, c-Jun, and c-Fos, compared with cells exposed to uniform laminar shear stress or maintained under static conditions. In addition, individual cells display a heterogeneity in responsiveness to disturbed flow, as measured by the amount of NF-kappaB, Egr-1, c-Jun, and c-Fos in their nuclei. This differential regulation of transcription factor expression by disturbed versus uniform laminar shear stress indicates that regional differences in blood flow patterns in vivo-in particular, the occurrence of spatial shear stress gradients-may represent important local modulators of endothelial gene expression at anatomic sites predisposed for atherosclerotic development.

Measuring actin dynamics in endothelial cells.

Cytoplasmic actin distributes between monomeric and filamentous phases in cells. As cells crawl, actin polymerizes near the plasma membrane of expanding peripheral cytoplasm and depolymerizes elsewhere. Thus, the finite actin filament lifetime, the diffusivity of actin monomer, and the distribution of actin between the polymer and monomer phases are key parameters in cell motility. The dynamics of cellular actin can be determined by following the evolution of fluorescence in the techniques of photoactivated fluorescence (PAF) or fluorescence recovery after photobleaching (FRAP) of microinjected actin derivatives. A mathematical model is discussed that measures monomer diffusion coefficients, filament turnover rates, and the fraction of actin polymerized from measurements of the evolution of fluorescence from a photoactivated band [Tardy et al. (1995) Biophys. J., 69:1674-1682; McGrath et al. (1998) Biophys. J., in press]. Applying this model to subconfluent endothelial cells shows that approximately 40% of the actin is polymer and that these filaments turn over on average every 6 minutes. This report discusses how PAF and FRAP can be combined with more traditional biochemistry to probe actin cytoskeleton remodeling in endothelial cells.

Simultaneous measurements of actin filament turnover, filament fraction, and monomer diffusion in endothelial cells.

The analogous techniques of photoactivation of fluorescence (PAF) and fluorescence recovery after photobleaching (FRAP) have been applied previously to the study of actin dynamics in living cells. Traditionally, separate experiments estimate the mobility of actin monomer or the lifetime of actin filaments. A mathematical description of the dynamics of the actin cytoskeleton, however, predicts that the evolution of fluorescence in PAF and FRAP experiments depends simultaneously on the diffusion coefficient of actin monomer, D, the fraction of actin in filaments, FF, and the lifetime of actin filaments, tau (, Biophys. J. 69:1674-1682). Here we report the application of this mathematical model to the interpretation of PAF and FRAP experiments in subconfluent bovine aortic endothelial cells (BAECs). The following parameters apply for actin in the bulk cytoskeleton of subconfluent BAECs. PAF: D = 3.1 +/- 0.4 x 10(-8) cm2/s, FF = 0.36 +/- 0.04, tau = 7.5 +/- 2.0 min. FRAP: D = 5.8 +/- 1.2 x 10(-8) cm2/s, FF = 0.5 +/- 0.04, tau = 4.8 +/- 0.97 min. Differences in the parameters are attributed to differences in the actin derivatives employed in the two studies and not to inherent differences in the PAF and FRAP techniques. Control experiments confirm the modeling assumption that the evolution of fluorescence is dominated by the diffusion of actin monomer, and the cyclic turnover of actin filaments, but not by filament diffusion. The work establishes the dynamic state of actin in subconfluent endothelial cells and provides an improved framework for future applications of PAF and FRAP.

Shear stress gradients remodel endothelial monolayers in vitro via a cell proliferation-migration-loss cycle.

Wall shear stress has been implicated in the genesis of atherosclerosis because a strong correlation exists between the location of developing arterial lesions and regions where particular gradients in stress occur. Studying the behavior of endothelial cells in such regions may contribute to our understanding of the disease etiology. We report the detailed migratory history of endothelial cells subjected to large shear stress gradients caused by a surface protuberance in an in vitro model system. The history of cell migration, cell division, and cell loss from the surface was continuously monitored in confluent human umbilical vein endothelial cell monolayers for 48 hours after the onset of flow. Individual cells were tracked using time-lapse video microscopy. In contrast to a uniform laminar flow field in which cells were observed to continually rearrange their relative position with no net migration, in a disturbed flow field there was a net migration directed away from the region of high shear gradient. This organized migration pattern under disturbed flow conditions was accompanied by more than a twofold increase in cell motility. In addition, cell division increased in the vicinity of the flow separation (maximum shear stress gradient of 34 dyne/cm2 per mm) whereas cell loss was increased upstream and downstream in the regions where the shear gradient diminishes. These data suggest a steady cell proliferation-migration-loss cycle and indicate that local shear stress gradient may play a key role in the morphological remodeling of the vascular endothelium in vivo.

Mechanical remodeling of the endothelial surface and actin cytoskeleton induced by fluid flow.

OBJECTIVE: The mechanism by which cultured endothelial cells respond to shear stress is controversial. The cell surface and cytoskeleton are involved, but their roles are undefined. In this study, previously unknown changes in the surface detail and actin cytoskeleton of bovine aortic endothelial cells were identified. METHODS: Actin filament content and filament number in resting and flow-oriented cells were determined by biochemical assays. The three-dimensional organization of the actin cytoskeleton in cells was defined in the confocal microscope and in the electron microscope after rapid-freezing, freeze-drying, and metal coating of detergent-permeabilized cells. RESULTS: Endothelial cells have smooth apical membranes in situ. However, cultured cells exhibit surface microvilli which increase the apical surface area, exposing the ruffled surface to forces from fluid flow and potentially enhancing cell interactions with blood-borne white cells. Stereoscopic micrographs show that stress fibers are integrated into a complex distributed cytoplasmic structural actin network (DCSA). This lattice is formed by actin filaments that frequently cross and connect to each other, stress fibers, and microfilaments and microtubules. The cytoskeletons of cells cultured in static media lack apparent order when compared to cytoskeletons from cells which have been exposed to 24 hours of laminar flow. CONCLUSIONS: The DCSA physically connects the apical and basal cell membranes and fills the volume between nucleus and membrane, providing a mechanism for transmitting mechanical forces across cells and a signaling pathway from membrane to nucleus. Stress fibers increase the mechanical modulus of the DCSA, although this increase is probably unnecessary to withstand the increase in shear stress caused by blood flow in vivo. This implies that actin rearrangements are not required for mechanical integrity, but serve an alternate function.

Prospects for telediagnosis using ultrasound.

Ultrasound imaging is currently used as a primary diagnostic tool in cardiology, abdominal disorders, pulmonary medicine, trauma, and obstetrics. Because of its relatively low capital and operating costs as well as its growth potential, it represents one of the major diagnostic modalities of future health care. However, the use of ultrasonography as a mobile and powerful modality is controlled by the availability of a highly skilled technician to acquire the images and an experienced physician to interpret them. This paper discusses the technology required to increase the availability of a diagnosing physician by employing telerobotics. With this technology, the physician can guide the motion of the transducer by the technician from a remote location. Thus, the physician controls the examination and renders the diagnosis. It is shown that communication lines at 1.5 Mbits/s (T-1 speed) can, with appropriate compression, support both real-time viewing of the ultrasound images and telerobotic manipulation of the transducer. The incremental costs of telediagnosis for an examination are estimated to be a small fraction of the base charges and significantly less than the expense of bringing a physician to a remote location or transporting a patient to a regional medical center. Telediagnosis can, in addition, provide benefits from immediate interpretation and consultation that cannot be duplicated using store-and-forward scenarios.

Theoretical estimates of mechanical properties of the endothelial cell cytoskeleton.

Current modeling of endothelial cell mechanics does not account for the network of F-actin that permeates the cytoplasm. This network, the distributed cytoplasmic structural actin (DCSA), extends from apical to basal membranes, with frequent attachments. Stress fibers are intercalated within the network, with similar frequent attachments. The microscopic structure of the DCSA resembles a foam, so that the mechanical properties can be estimated with analogy to these well-studied systems. The moduli of shear and elastic deformations are estimated to be on the order of 10(5) dynes/cm2. This prediction agrees with experimental measurements of the properties of cytoplasm and endothelial cells reported elsewhere. Stress fibers can potentially increase the modulus by a factor of 2-10, depending on whether they act in series or parallel to the network in transmitting surface forces. The deformations produced by physiological flow fields are of insufficient magnitude to disrupt cell-to-cell or DCSA cross-linkages. The questions raised by this paradox, and the ramifications of implicating the previously unreported DCSA as the primary force transmission element are discussed.

Interpreting photoactivated fluorescence microscopy measurements of steady-state actin dynamics.

A continuum model describing the steady-state actin dynamics of the cytoskeleton of living cells has been developed to aid in the interpretation of photoactivated fluorescence experiments. In a simplified cell geometry, the model assumes uniform concentrations of cytosolic and cytoskeletal actin throughout the cell and no net growth of either pool. The spatiotemporal evolution of the fluorescent actin population is described by a system of two coupled linear partial-differential equations. An analytical solution is found using a Fourier-Laplace transform and important limiting cases relevant to the design of experiments are discussed. The results demonstrate that, despite being a complex function of the parameters, the fluorescence decay in photoactivated fluorescence experiments has a biphasic behavior featuring a short-term decay controlled by monomer diffusion and a long-term decay governed by the monomer exchange rate between the polymerized and unpolymerized actin pools. This biphasic behavior suggests a convenient mechanism for extracting the parameters governing the fluorescence decay from data records. These parameters include the actin monomer diffusion coefficient, filament turnover rate, and ratio of polymerized to unpolymerized actin.

Shear stress selectively upregulates intercellular adhesion molecule-1 expression in cultured human vascular endothelial cells.

Hemodynamic forces induce various functional changes in vascular endothelium, many of which reflect alterations in gene expression. We have recently identified a cis-acting transcriptional regulatory element, the shear stress response element (SSRE), present in the promoters of several genes, that may represent a common pathway by which biomechanical forces influence gene expression. In this study, we have examined the effect of shear stress on endothelial expression of three adhesion molecules: intercellular adhesion molecule-1 (ICAM-1), which contains the SSRE in its promoter, and E-selectin (ELAM-1) and vascular cell adhesion molecule-1 (VCAM-1), both of which lack the SSRE. Cultured human umbilical vein endothelial cells, subjected to a physiologically relevant range of laminar shear stresses (2.5-46 dyn/cm2) in a cone and plate apparatus for up to 48 h, showed time-dependent but force-independent increases in surface immunoreactive ICAM-1. Upregulated ICAM-1 expression was correlated with increased adhesion of the JY lymphocytic cell line. Northern blot analysis revealed increased ICAM-1 transcript as early as 2 h after the onset of shear stress. In contrast, E-selectin and vascular cell adhesion molecule-1 transcript and cell-surface protein were not upregulated at any time point examined. This selective regulation of adhesion molecule expression in vascular endothelium suggests that biomechanical forces, in addition to humoral stimuli, may contribute to differential endothelial gene expression and thus represent pathophysiologically relevant stimuli in inflammation and atherosclerosis.

Fluid flow modulates vascular endothelial cytosolic calcium responses to adenine nucleotides.

OBJECTIVE: To determine whether fluid flow influences the action of soluble vasoactive agonists on vascular endothelium. METHODS: Confluent monolayers of bovine aortic endothelial cells (BAEC) were cultured on glass coverslips, prelabeled with the Ca(2+)-sensitive dye fura-2, and placed in a parallel-plate flow chamber designed to generate defined laminar fluid flow. Cytosolic free Ca2+ concentration ([Ca2+]i) in individual BAEC was monitored during perfusion with medium containing adenine nucleotide under defined flow conditions. RESULTS: Continuous perfusion with ATP (0.3-3.0 microM) or ADP (0.1-1.0 microM) evoked repetitive oscillations in [Ca2+]i in individual BAEC. The frequency of the [Ca2+]i oscillations was dependent on both nucleotide concentration and levels of applied shear stress; at constant bulk concentration of nucleotide, the frequency increased with shear stress. Stopping flow in the continuous presence of agonists immediately extinguished the oscillatory response. Elimination of extracellular Ca2+ did not inhibit the [Ca2+]i oscillations. In the presence of nonhydrolyzable nucleotide analog, ATP gamma S or ADP beta S, application of flow resulted in similar shear-dependent [Ca2+]i oscillations, suggesting that flow modulation of the [Ca2+]i response was not simply due to depletion of ATP or ADP in the vicinity of BAEC monolayers as a result of hydrolysis of nucleotides by ectonucleotidases. CONCLUSIONS: These findings suggest that local hemodynamic conditions may modulate the action of vasoactive agents on the vascular endothelium in vivo.

Platelet-derived growth factor B chain promoter contains a cis-acting fluid shear-stress-responsive element.

The endothelial lining of blood vessels is constantly exposed to fluid mechanical forces generated by flowing blood. In vitro application of fluid shear stresses to cultured endothelial cells influences the expression of multiple genes, as reflected by changes in their steady-state mRNA levels. We have utilized the B chain of platelet-derived growth factor (PDGF-B) as a model to investigate the mechanisms of shear-stress-induced gene regulation in cultured bovine aortic endothelial cells (BAECs). Northern blot analysis revealed elevated endogenous PDGF-B transcript levels in BAECs, after exposure to a physiological level of laminar shear stress (10 dynes/cm2; 1 dyne = 100 mN) for 4 h. A transfected reporter gene, consisting of a 1.3-kb fragment of the human PDGF-B promoter coupled to chloramphenicol acetyltransferase (CAT), indicated a direct effect on transcriptional activity. Transfection of a series of PDGF-B-CAT deletion mutants led to the characterization of a cis-acting component within the PDGF-B promoter that was necessary for shear-stress responsiveness. In gel-shift assays, overlapping oligonucleotide probes of this region formed several protein-DNA complexes with nuclear extracts prepared from both static and shear-stressed BAECs. A 12-bp component (CTCTCAGAGACC) was identified that formed a distinct pattern of complexes with nuclear proteins extracted from shear-stressed BAECs. This shear-stress-responsive element does not encode binding sites for any known transcription factor but does contain a core binding sequence (GAGACC), as defined by deletion mutation in gel-shift assays. Interestingly, this putative transcription factor binding site is also present in the promoters of certain other endothelial genes, including tissue plasminogen activator, intercellular adhesion molecule 1, and transforming growth factor beta 1, that also are induced by shear stress. Thus, the expression of PDGF-B and other pathophysiologically relevant genes in vascular endothelium appears to be regulated, in part, by shear-stress-induced transcription factors interacting with a common promoter element.

Regulation of adenine nucleotide concentration at endothelium-fluid interface by viscous shear flow.

The action of adenine nucleotides on vascular endothelial cells is apparently mediated by the local flow conditions. Because nucleotides are sequentially degraded from ATP-->ADP-->AMP-->adenosine by ecto-enzymes at the endothelial surface, it has been hypothesized that the observed flow effect is caused by the flow-dependent change of nucleotide concentration at the cell surface. In this study, we have calculated the concentration profiles of adenine nucleotides at the cell surface under flow conditions encountered in an in vitro parallel-plate flow system, as has been used in several related experimental studies. When medium containing uniformly distributed ATP is perfused over endothelial monolayers, our results show that ATP concentration in the cell vicinity gradually decreases in the streamwise direction as a result of enzymatic degradation. This hydrolysis of ATP results in the generation of ADP, and ADP concentration in turn gradually increases at the cell surface. The concentration profiles of nucleotides are dependent on the levels of applied wall shear rate. As the corresponding shear stress increases from 0.1 to 30 dynes/cm2, ATP concentration at the cell surface at the center of coverslip increases from 0.66 to 0.93. Under no-flow conditions, our model predicts a steady decline of ATP concentration and a transient increase of ATP-derived ADP, comparable to the published results of previous experiments. These numerical results, combined with our recent experimental data, provide insights into the cellular mechanisms by which hemodynamic flow modulates the effects of vasoactive agents on endothelium.

Vascular endothelium responds to fluid shear stress gradients.

In vitro investigations of the responses of vascular endothelium to fluid shear stress have typically been conducted under conditions where the time-mean shear stress is uniform. In contrast, the in vitro experiments reported here have re-created the large gradients in surface fluid shear stress found near arterial branches in vivo; specifically, we have produced a disturbed-flow region that includes both flow separation and reattachment. Near reattachment regions, shear stress is small but its gradient is large. Cells migrate away from this region, predominantly in the downstream direction. Those that remain divide at a rate that is high compared with that of cells subjected to uniform shear. We speculate that large shear stress gradients can induce morphological and functional changes in the endothelium in regions of disturbed flow in vivo and thus may contribute to the formation of atherosclerotic lesions.

The distribution of fluid forces on model arterial endothelium using computational fluid dynamics.

Numerical calculations are used in conjunction with linear perturbation theory to analyze the problem of laminar flow of an incompressible fluid over a wavy surface which approximates a monolayer of vascular endothelial cells. These calculations model flow conditions in an artery very near the vessel wall at any instant in time, providing a description of the velocity field with detail that would be difficult to identify experimentally. The surface pressure and shear stress distributions are qualitatively similar for linear theory and numerical computations. However, the results diverge as the amplitude of surface undulation is increased. The shear stress gradient along the cell model surface is reduced for geometries which correspond to aligned endothelial cells (versus nonaligned geometries).

Fluid shear stress modulates cytosolic free calcium in vascular endothelial cells.

Cytosolic free Ca2+ concentration ([Ca2+]i) was monitored in single and groups of fura-2-loaded bovine aortic endothelial cells (BAEC) during exposure to laminar fluid shear stress. Application of a step increase in shear stress from 0.08 to 8 dyn/cm2 to confluent BAEC monolayers resulted in a transient increase in [Ca2+]i, which attained a peak value in 15-40 s, followed by a decline to baseline within 40-80 s. The magnitude of the [Ca2+]i responses increased with applied shear stress over the range of 0.2-4 dyn/cm2 and reached a maximum at greater than 4 dyn/cm2. Transient oscillations in [Ca2+]i with gradually diminishing amplitude were observed in individual cells subjected to continuous high shear stress. Elimination of extracellular Ca2+ with ethylene glycol-bis(beta-aminoethyl ether)-N,N,N',N'-tetraacetic acid, blockade of Ca2+ entry with lanthanum, depolarization of the cell membrane with high K+, and preconditioning of BAEC in steady laminar flow had little effect on the [Ca2+]i response. In the presence of ATP or ADP, application of shear stress caused repetitive oscillations in [Ca2+]i in single BAEC, whose frequency was dependent on both agonist concentration and the magnitude of applied shear stress. However, apyrase, an ATPase and ADPase, did not inhibit the shear-induced [Ca2+]i responses in standard medium (no added ATP or ADP), suggesting that the shear-induced [Ca2+]i response is not due to ATP released by endothelial cells.