Streaming potential of intact wet bone

The conversion of mechanical loads to bioelectrical signals in bone have been suggested to control repair and remodeling. These signals in wet bone are attributed to the electrokinetic behavior where mechanical forces cause electrical signals due to motion of an ion carrying extracellular fluid in the bone matrix (streaming potentials). Streaming potential experiments were performed on control and chemically treated intact wet bone plugs in aphosphate and phosphate buffers to examine the contribution of bone constituents to the electrokinetic behavior of bone tissue. Data indicate that the organic constituents of bone dominate streaming potentials. Slopes of streaming potential vs pressure are related to the electrokinetic (zeta) potential. The slopes should be analyzed in the low pressure region where data is mainly linear. Comparisons of estimated zeta potentials from streaming potentials with existing data obtained by particle electrophoresis showed similar trends.     

Three-dimensional inhomogeneous triphasic finite-element analysis of physical signals and solute transport in human intervertebral disc under axial compression

A 3D inhomogeneous finite-element model for charged hydrated soft tissues containing charged/uncharged solutes was developed and applied to analyze the mechanical, chemical, and electrical signals within the human intervertebral disc during an axial unconfined compression. The effects of tissue properties and boundary conditions on the physical signals and the transport of fluid and solute were investigated. The numerical simulation showed that, during disc compression, the fluid pressurization and the effective (von Misses) solid stress were more pronounced in the annulus fibrosus (AF) region near the interface between AF and nucleus pulposus (NP). In NP, the distributions of the fluid pressure, effective stress, and electrical potential were more uniform than those in AF. The electrical signals were very sensitive to fixed charge density. Changes in material properties of NP (water content, fixed charge density, and modulus) affected fluid pressure, electrical potential, effective stress, and solute transport in the disc. This study is important for understanding disc biomechanics, disc nutrition, and disc mechanobiology.     

Mechanically induced electrical potentials of articular cartilage

While there is increasing evidence that chondrocytes are affected by mechanically induced stimuli, endogenous force-related electrical potentials within articular cartilage have been so far observed only in-vitro. Using a porcine ex-vivo model (German Land Race), 8 knee joints were explanted and exposed to mechanical force (up to 800 N) using a special device. Electrodes were inserted into the cartilage matrix. With an amplifier and an A/D transducer the changes of electrical voltage between the electrodes as well as those of the force were recorded online and simultaneously on a computer. Additionally, we located one pair of electrodes on the surface of the cartilage tissue to detect electrical fields outside the cartilage tissue. In relation to the applied force we observed that electrical potentials derived from inside and outside the articular cartilage showed a correspondence. When an alternating force with an amplitude of 360 N and a frequency of about 0.2 Hz was periodically applied, we measured peak amplitudes ranging from 2.1 to 5.5 mV within the cartilage tissue with electrical negativity within the weight bearing area of the cartilage tissue. The measured voltages depended on the applied force, the location of the electrodes, and on anatomical variations. We found an almost linear relation between the magnitude of the applied force and the recorded voltage. With the help of the electrodes located outside and within the cartilage tissue, we were able to show that force dependent fields are generated inside the cartilage. There are several theories explaining the origin of these electrical phenomena, many of them focusing on the negative charges of the proteoglycans in relation to the flow of interstitial fluid and ions under compression. However, the consequences of these phenomena are yet not clear.     

A dynamical study of the mechanical stimuli and tissue differentiation within a CaP scaffold based on micro-CT finite element models

The control of the mechanical stimuli transmitted to the cells is critical for the design of functional scaffolds for tissue engineering. The objective of this study was to investigate the dynamics of the mechanical stimuli transmitted to the cells during tissue differentiation in an irregular morphology scaffold under compressive load and perfusion flow. A calcium phosphate-based glass porous scaffold was used. The solid phase and the fluid flow within the pores were modeled as linear elastic solid material and Newtonian fluid, respectively. In the fluid model, different levels of viscosity were used to simulate tissue differentiation. Compressive strain of 0.5% and fluid flow with constant inlet velocity of 10 μm/s or constant inlet pressure of 3 Pa were applied. Octahedral shear strain and fluid shear stress were used as mechano-regulatory stimuli. For constant inlet velocity, stimuli equivalent to bone were predicted in 80% of pore volume for the case of low tissue viscosity. For the cases of high viscosity, fluctuations between stimuli equivalent to tissue formation and cell death were predicted due to the increase in the fluid shear stress when tissue started to fill pores. When constant pressure was applied, stimuli equivalent to bone were predicted in 62% of pore volume when low tissue viscosity was used and 42% when high tissue viscosity was used. This study predicted critical variations of fluid shear stress when cells differentiated. If these variations are not controlled in vitro, they can impede the formation of new matured tissue.     

A comparison between mechano-electrochemical and biphasic swelling theories for soft hydrated tissues

Biological tissues like intervertebral discs and articular cartilage primarily consist of interstitial fluid, collagen fibrils and negatively charged proteoglycans. Due to the fixed charges of the proteoglycans, the total ion concentration inside the tissue is higher than in the surrounding synovial fluid (cation concentration is higher and the anion concentration is lower). This excess of ion particles leads to an osmotic pressure difference, which causes swelling of the tissue. In the last decade several mechano-electrochemical models, which include this mechanism, have been developed. As these models are complex and computationally expensive, it is only possible to analyze geometrically relatively small problems. Furthermore, there is still no commercial finite element tool that includes such a mechano-electrochemical theory. Lanir (Biorheology, 24, pp. 173-187, 1987) hypothesized that electrolyte flux in articular cartilage can be neglected in mechanical studies. Lanir's hypothesis implies that the swelling behavior of cartilage is only determined by deformation of the solid and by fluid flow. Hence, the response could be described by adding a deformation-dependent pressure term to the standard biphasic equations. Based on this theory we developed a biphasic swelling model. The goal of the study was to test Lanir's hypothesis for a range of material properties. We compared the deformation behavior predicted by the biphasic swelling model and a full mechano-electrochemical model for confined compression and 1D swelling. It was shown that, depending on the material properties, the biphasic swelling model behaves largely the same as the mechano-electrochemical model, with regard to stresses and strains in the tissue following either mechanical or chemical perturbations. Hence, the biphasic swelling model could be an alternative for the more complex mechano-electrochemical model, in those cases where the ion flux itself is not the subject of the study. We propose thumbrules to estimate the correlation between the two models for specific problems.     

The influence of the fixed negative charges on mechanical and electrical behaviors of articular cartilage under unconfined compression

Unconfined compression test has been frequently used to study the mechanical behaviors of articular cartilage, both theoretically and experimentally. It has also been used in explant and gel-cell-complex studies in tissue engineering. In biphasic and poroelastic theories, the effect of charges fixed on the proteoglycan macromolecules in articular cartilage is embodied in the apparent compressive Young's modulus and the apparent Poisson's ratio of the tissue, and the fluid pressure is considered to be the portion above the osmotic pressure. In order to understand how proteoglycan fixed charges might affect the mechanical behaviors of articular cartilage, and in order to predict the osmotic pressure and electric fields inside the tissue in this experimental configuration, it is necessary to use a model that explicitly takes into account the charged nature of the tissue and the flow of ions within its porous interstices. In this paper, we used a finite element model based on the triphasic theory to study how fixed charges in the porous-permeable soft tissue can modulate its mechanical and electrochemical responses under a step displacement in unconfined compression. The results from finite element calculations showed that: 1) A charged tissue always supports a larger load than an uncharged tissue of the same intrinsic elastic moduli. 2) The apparent Young's modulus (the ratio of the equilibrium axial stress to the axial strain) is always greater than the intrinsic Young's modulus of an uncharged tissue. 3) The apparent Poisson's ratio (the negative ratio of the lateral strain to the axial strain) is always larger than the intrinsic Poisson's ratio of an uncharged tissue. 4) Load support derives from three sources: intrinsic matrix stiffness, hydraulic pressure and osmotic pressure. Under the unconfined compression, the Donnan osmotic pressure can constitute between 13%-22% of the total load support at equilibrium. 5) During the stress-relaxation process following the initial instant of loading, the diffusion potential (due to the gradient of the fixed charge density and the associated gradient of ion concentrations) and the streaming potential (due to fluid convection) compete against each other. Within the physiological range of material parameters, the polarity of the electric potential depends on both the mechanical properties and the fixed charge density (FCD) of the tissue. For softer tissues, the diffusion effects dominate the electromechanical response, while for stiffer tissues, the streaming potential dominates this response. 6) Fixed charges do not affect the instantaneous strain field relative to the initial equilibrium state. However, there is a sudden increase in the fluid pressure above the initial equilibrium osmotic pressure. These new findings are relevant and necessary for the understanding of cartilage mechanics, cartilage biosynthesis, electromechanical signal transduction by chondrocytes, and tissue engineering.     

Wear particle diffusion and tissue differentiation in TKA implant fibrous interfaces

In the context of mechanical loosening, we studied the hypothesis that wear-particle migration in the fibrous membrane under tibial plateaus after total knee arthroplasty can be explained by the pumping effects of the interstitial fluid in the tissue. Further, as a secondary objective we investigated the possibility that interface-tissue differentiation is influenced by interstitial fluid flow and strain, as mechanical effects of interface motions. For comparative reasons, we analyzed a previously published simplified two-dimensional finite-element model, this time assuming biphasic tissue properties. We wanted to determine hydrostatic pressure and flow velocities in the fluid phase, in addition to stresses and strains, for time-dependent loading of the plateau. We found that fluid flow in the interface was extremely slow, except in the periphery. Hence, loosening due to particle-induced bone resorption appears improbable. The results, however, do support the idea that particles migrate with fluid flow, when such flow occurs. Where fibrous tissue tends to be prominent in reality, the fluid is repeatedly extruded and reabsorbed in the model. Where these values are low, fibrocartilage is commonly found. When material properties were varied to subsequently represent fibrocartilage and two stages of mineralization, the strains and fluid velocities is reduced. Fluid pressure, however, did not change. Our results refute the hypothesis that wear particles are pumped through the interface tissue below a TKA but support the hypothesis that interface tissue type and loosening processes are influenced by mechanical tissue variables such as tissue strain and interstitial fluid velocity.     

A finite element study of mechanical stimuli in scaffolds for bone tissue engineering

Mechanical stimuli are one of the factors that affect cell proliferation and differentiation in the process of bone tissue regeneration. Knowledge on the specific deformation sensed by cells at a microscopic level when mechanical loads are applied is still missing in the development of biomaterials for bone tissue engineering. The objective of this study was to analyze the behavior of the mechanical stimuli within some calcium phosphate-based scaffolds in terms of stress and strain distributions in the solid material phase and fluid velocity, fluid pressure and fluid shear stress distributions in the pores filled of fluid, by means of micro computed tomographed (CT)-based finite element (FE) models. Two samples of porous materials, one of calcium phosphate-based cement and another of biodegradable glass, were used. Compressive loads equivalent to 0.5% of compression applied to the solid material phase and interstitial fluid flows with inlet velocities of 1, 10 and 100 microm/s applied to the interconnected pores were simulated, changing also the inlet side and the viscosity of the medium. Similar strain distributions for both materials were found, with compressive and tensile strain maximal values of 1.6% and 0.6%, respectively. Mean values were consistent with the applied deformation. When 10 microm/s of inlet fluid velocity and 1.45 Pas viscosity, maximal values of fluid velocity were 12.76 mm/s for CaP cement and 14.87 mm/s for glass. Mean values were consistent with the inlet ones applied, and mean values of shear stress were around 5 x 10(-5)Pa. Variations on inlet fluid velocity and fluid viscosity produce proportional and independent changes in fluid velocity, fluid shear stress and fluid pressure. This study has shown how mechanical loads and fluid flow applied on the scaffolds cause different levels of mechanical stimuli within the samples according to the morphology of the materials.     

理解生物组织细胞力学信号传导机制的必要分析：以关节软骨为例

关节软骨是覆盖于骨关节面的一薄层低摩擦、耐磨损、负重的水化组织。这种功能通过散在镶嵌于软骨组织深层致密胞外组织内的软骨细胞的代谢和生物合成作用来维持。这种代谢和生物合成进程很大一部分是由物理因素,如应力、张力、电流及电势、液压及渗透压等来调控的。这两者都存在于这一带电的、渗透性的基质当中。本文主要讲述关节软骨及软骨细胞的机械一电化学行为及理论模型的最近进展,同时着眼于软骨细胞生物合成活动的物理调控及其对于组织维持、功能组织工程软骨修复及再生医学的意义。     

Les flux d’eau dans le tendon sous-contrainte

It is necessary before studying tendon mechanobiology to understand its liquid phase behaviour under stress. As observed experimentally (particularly with NMR) the tendon viscoelastic characteristics depend on water content and tissue geometric modifications cause physicochemical changes, which induce fluid diffusion and determine its properties (pressure, flow directions, diffusion pathways). All these experimental data have been collected with our current knowledge of tendinous microstructure in a model that describes the tendon fluid flow under stress.     

Mechanical regulation of signaling pathways in bone

A wide range of cell types depend on mechanically induced signals to enable appropriate physiological responses. The skeleton is particularly dependent on mechanical information to guide the resident cell population towards adaptation, maintenance and repair. Research at the organ, tissue, cell and molecular levels has improved our understanding of how the skeleton can recognize the functional environment, and how these challenges are translated into cellular information that can site-specifically alter phenotype. This review first considers those cells within the skeleton that are responsive to mechanical signals, including osteoblasts, osteoclasts, osteocytes and osteoprogenitors. This is discussed in light of a range of experimental approaches that can vary parameters such as strain, fluid shear stress, and pressure. The identity of mechanoreceptor candidates is approached, with consideration of integrins, pericellular tethers, focal adhesions, ion channels, cadherins, connexins, and the plasma membrane including caveolar and non-caveolar lipid rafts and their influence on integral signaling protein interactions. Several mechanically regulated intracellular signaling cascades are detailed including activation of kinases (Akt, MAPK, FAK), β-catenin, GTPases, and calcium signaling events. While the interaction of bone cells with their mechanical environment is complex, an understanding of mechanical regulation of bone signaling is crucial to understanding bone physiology, the etiology of diseases such as osteoporosis, and to the development of interventions to improve bone strength.     

The electrical potential difference through the foot epithelium of the snail Achatina achatina, Lameere during mechanical and chemical stimulation

An important electrophysiological variable--the transepithelial potential difference reflects the electrogenic transepithelial ion currents, which are produced and modified by ion transport processes in polarized cells of epithelium. These processes result from coordinated function of transporters in apical and basolateral cell membranes and have been observed in all epithelial tissues studied so far. The experiments were performed on isolated specimens of snail foot. In the experiments, the baseline transepithelial electrical potential difference--PD, changes of transepithelial difference during mechanical stimulation--dPD and the transepithelial resistance were measured with an Ussing apparatus. A total of 60 samples of foot ventral surface of 28 snails were studied. The transepithelial electrical potential difference of isolated foot ranged from -6.0 to 10.0 mV under different experimental conditions. Mechanical stimulation of foot ventral surface caused changes of electrogenic ion transport, observed as transient hyperpolarization (electrical potential difference became more positive). When the transepithelial electrical potential difference decreased during stimulation, the reaction was described as depolarization. When amiloride and bumetanide were added to the stimulating fluid so that the sodium and chloride ion transport pathways were inhibited, prolonged depolarization occurred. Under the influence of different stimuli: mechanical (gentle rinsing), chemical (changes of ion concentrations) and pharmacological (application of ion inhibitors), transient changes of potential difference (dPD) were evoked, ranging from about -0.7 to almost 2.0 mV. Changes in transepithelial potential difference of the pedal surface of the snail's foot related to these physiological stimuli are probably involved in the locomotion of the animal and are under control of the part of the nervous system in which tachykinin related peptides (TRP) act as transmitters.     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Time-Resolved Analysis of Signals Involved in Systemic Induction of Pin2 Gene Expression

The systemic induction of proteinase inhibitor genes in tomato plants is mediated either by electrical signals, hydraulic signals or chemical messengers. In the present study we analyzed the effects of mechanical wounding, heat treatment and electrical current application on wild-type tomato plants (Lycopersicon esculentum Mill, cv Moneymaker) and ABA-deficient mutants of tomato (sitiens). Kinetic studies revealed that systemic Pin2 gene expression could be slightly induced by the fast transient membrane potential change which left the damaged leaf within 30–60s after wounding. Moreover, a signal leaving the damaged tissue between 2 and 4 minutes after wounding was responsible for a significant amplification of Pin2 gene expression. This signal could either be a decrease in turgor pressure, which occurred 3–4min after treatment, or a slow electrical transient. In addition, mechanical wounding and electrical current seem to involve ABA to induce changes in membrane potential and to promote Pin2 gene expression. In contrast, heat triggers fast and slow electrical transients leading to an induction of Pin2 gene expression within the plant independently of ABA. Turgor pressure, in turn, is presumably adjusted in relation to ionic movements across the membrane, elucidated by membrane potential recordings. In conclusion, wound-induced changes in membrane potential seem to be dependent on the endogenous level of ABA. These shifts in membrane potentials, in turn, are involved in regulation of turgor pressure within the plant.     

Extracellular matrix, mechanotransduction and structural hierarchies in heart tissue engineering

The spatial and temporal scales of cardiac organogenesis and pathogenesis make engineering of artificial heart tissue a daunting challenge. The temporal scales range from nanosecond conformational changes responsible for ion channel opening to fibrillation which occurs over seconds and can lead to death. Spatial scales range from nanometre pore sizes in membrane channels and gap junctions to the metre length scale of the whole cardiovascular system in a living patient. Synchrony over these scales requires a hierarchy of control mechanisms that are governed by a single common principle: integration of structure and function. To ensure that the function of ion channels and contraction of muscle cells lead to changes in heart chamber volume, an elegant choreography of metabolic, electrical and mechanical events are executed by protein networks composed of extracellular matrix, transmembrane integrin receptors and cytoskeleton which are functionally connected across all size scales. These structural control networks are mechanoresponsive, and they process mechanical and chemical signals in a massively parallel fashion, while also serving as a bidirectional circuit for information flow. This review explores how these hierarchical structural networks regulate the form and function of living cells and tissues, as well as how microfabrication techniques can be used to probe this structural control mechanism that maintains metabolic supply, electrical activation and mechanical pumping of heart muscle. Through this process, we delineate various design principles that may be useful for engineering artificial heart tissue in the future.     

Modeling viscous dissipation during vocal fold contact: the influence of tissue viscosity and thickness with implications for hydration

The mechanics of vocal fold contact during phonation is known to play a crucial role in both normal and pathological speech production, though the underlying physics is not well understood. Herein, a viscoelastic model of the stresses during vocal fold contact is developed. This model assumes the cover to be a poroelastic structure wherein interstitial fluid translocates in response to mechanical squeezing. The maximum interstitial fluid pressure is found to generally increase with decreasing viscous dissipation and/or decreasing tissue elasticity. A global minimum in the total contact stress, comprising interstitial fluid pressure and elastic stress in the tissue, is observed over the studied dimensionless parameter range. Interestingly, physiologically reasonable estimates for the governing parameters fall within this global minimum region. The model is validated against prior experimental and computational work, wherein the predicted contact stress magnitude and impact duration agree well with published results. Lastly, observations of the potential relationship between vocal fold hydration and increased risk of tissue damage are discussed based upon model predictions of stress as functions of cover layer thickness and viscosity.     

Mineral and matrix contributions to rigidity in fracture healing

The purpose of this study was to investigate the relationships among selected properties of fracture callus: bending rigidity, tissue density, mineral density, matrix density and mineral-to-matrix ratio. The experimental model was an osteotomized canine radius in which the development of the fracture callus was modified by electrical stimulation with various levels of direct current. This resulted in a range of values for the selected properties of the callus, determined post mortem at 7 weeks after osteotomy. We found that the rigidity (R) of the bone-callus combination obeyed relationships of the form R = axb, where x is the tissue density, mineral density, matrix density or the mineral-to-matrix ratio of the repair tissue. These are analogous to power-law relationships found in studies of compact and cancellous bone. The results suggest that fracture callus at 7 weeks after osteotomy in canine radius behaves more like immature compact bone than cancellous bone in its mechanical and physicochemical properties. The present study demonstrates the feasibility of developing non-invasive in vivo densitometric methods to monitor fracture healing, since models may be developed that can predict mechanical properties from densitometric data. Further studies are needed to develop a refined model based on experimental data on the mechanical and physicochemical properties and microstructure of fracture callus at different stages of healing.     

Differences in glycosylation and sperm-egg binding inhibition of pregnancy-related glycodelin

Glycodelin is a glycoprotein produced in many glands, particularly those of reproductive tissues. It appears as different glycoforms in amniotic fluid (glycodelin-A) and seminal plasma (glycodelin-S), but only glycodelin-A inhibits gamete adhesion. In the present study, glycodelin from secretory-phase endometrium, first-trimester pregnancy decidua, and midtrimester amniotic fluid was studied with respect to physicochemical properties, including glycosylation patterns and inhibitory activity of sperm-egg binding. Purified glycodelins from all these sources were similar in isoelectric focusing and in lectin immunoassays using lectins from Wisteria floribunda and Sambucus nigra. Likewise, the glycodelins inhibited sperm-egg binding in a dose-dependent manner, as measured by hemizona-binding assay. However, subtle quantitative physicochemical and biological differences were found between glycodelins from different sources as well as within the same tissue/fluid between different individuals. Differences were most pronounced between endometrial glycodelins from nonpregnancy and first-trimester pregnancy. The glycan structures studied by fast-atom bombardment mass spectrometry of individual amniotic fluid glycodelin-A samples also showed interindividual quantitative differences. In conclusion, glycodelins from different female reproductive tract tissues and amniotic fluid share substantial similarity, allowing all of them to be called glycodelin-A. However, these glycodelins exhibit quantitative physicochemical and functional differences between different sources and individuals.     

Physical signals and solute transport in human intervertebral disc during compressive stress relaxation: 3D finite element analysis

A 3D finite element model for charged hydrated soft tissues containing charged/uncharged solutes was developed based on the multi-phasic mechano-electrochemical mixture theory (Lai et al., J. Biomech. Eng. 113 (1991), 245-258; Gu et al., J. Biomech. Eng. 120 (1998), 169-180). This model was applied to analyze the mechanical, chemical and electrical signals within the human intervertebral disc during an unconfined compressive stress relaxation test. The effects of tissue composition [e.g., water content and fixed charge density (FCD)] on the physical signals and the transport rate of fluid, ions and nutrients were investigated. The numerical simulation showed that, during disc compression, the fluid pressurization was more pronounced at the center (nucleus) region of the disc while the effective (von Mises) stress was higher at the outer (annulus) region. Parametric analyses revealed that the decrease in initial tissue water content (0.7-0.8) increased the peak stress and relaxation time due to the reduction of permeability, causing greater fluid pressurization effect. The electrical signals within the disc were more sensitive to FCD than tissue porosity, and mechanical loading affected the large solute (e.g., growth factor) transport significantly, but not for small solute (e.g., glucose). Moreover, this study confirmed that the interstitial fluid pressurization plays an important role in the load support mechanism of IVD by sharing more than 40% of the total load during disc compression. This study is important for understanding disc biomechanics, disc nutrition and disc mechanobiology.     

A computational study of systemic hydration in vocal fold collision

Mechanical stresses develop within vocal fold (VF) soft tissues due to phonation-associated vibration and collision. These stresses in turn affect the hydration of VF tissue and thus influence voice health. In this paper, high-fidelity numerical computations are described, taking into account fully 3D geometry, realistic tissue and air properties, and high-amplitude vibration and collision. A segregated solver approach is employed, using sophisticated commercial solvers for both the VF tissue and glottal airflow domains. The tissue viscoelastic properties were derived from a biphasic formulation. Two cases were considered, whereby the tissue viscoelastic properties corresponded to two different volume fractions of the fluid phase of the VF tissue. For each case, hydrostatic stresses occurring as a result of vibration and collision were investigated. Assuming the VF tissue to be poroelastic, interstitial fluid movement within VF tissue was estimated from the hydrostatic stress gradient. Computed measures of overall VF dynamics (peak airflow velocity, magnitude of VF deformation, frequency of vibration and contact pressure) were well within the range of experimentally observed values. The VF motion leading to mechanical stresses within the VFs and their effect on the interstitial fluid flux is detailed. It is found that average deformation and vibration of VFs tend to increase the state of hydration of the VF tissue, whereas VF collision works to reduce hydration.