Cartilage stress-relaxation proceeds slower at higher compressive strains

Articular cartilage is the connective tissue which covers bone surfaces and deforms during in vivo activity. Previous research has investigated flow-dependent cartilage viscoelasticity, but relatively few studies have investigated flow-independent mechanisms. This study investigated polymer dynamics as an explanation for the molecular basis of flow-independent cartilage viscoelasticity. Polymer dynamics predicts that stress-relaxation will proceed more slowly at higher volumetric concentrations of polymer. Stress-relaxation tests were performed on cartilage samples after precompression to different strain levels. Precompression increases the volumetric concentration of cartilage biopolymers, and polymer dynamics predicts an increase in relaxation time constant. Stress-relaxation was slower for greater precompression. There was a significant correlation between the stress-relaxation time constant and cartilage volumetric concentration. Estimates of the flow-dependent timescale suggest that flow-dependent relaxation occurs on a longer timescale than presently observed. These results are consistent with polymer dynamics as a mechanism of cartilage viscoelasticity.     

Singular perturbation analysis of the nonlinear, flow-dependent compressive stress relaxation behavior of articular cartilage

The dominant mechanism giving rise to the viscoelastic response of articular cartilage during compression is the nonlinear diffusive interaction of the fluid and solid phases of the tissue as they flow relative to one another. The present study is concerned with the role of this interaction under uniaxial stress relaxation in compression. The model is a biphasic mixture of fluid and solid which incorporates the strain-dependent permeability found earlier from permeation experiments. When a ramp-displacement is imposed on the articular surface, simple, but accurate, asymptotic approximations are derived for the deformation and stress fields in the tissue for slow and moderately fast rates of compression. They are shown to agree very well with experiment and they provide a simple means for determining the material parameters. Moreover, they lead to important insights into the role of the flow-dependent viscoelastic nature of articular cartilage and other hydrated biological tissues.     

Tribology of polymer brush: microscale modelling and simulation

In order to examine the tribological behaviour of the polymer brush, we constructed a mesoscale particle model of a polymer brush system, based on a Brownian dynamics scheme. The polymer model consists of coarse-grained beads connected with harmonic springs. The Lennard-Jones type interaction is assumed between beads. The flow velocity is affected by the local packing fraction of beads. With this model, we executed a series of molecular dynamics simulation to investigate the mechanical responses of brush to shear flow.     

A comparison of cartilage stress-relaxation models in unconfined compression: QLV and stretched exponential in combination with fluid flow

Cartilage exhibits nonlinear viscoelastic behaviour. Various models have been proposed to explain cartilage stress relaxation, but it is unclear whether explicit modelling of fluid flow in unconfined compression is needed. This study compared Fung's quasi-linear viscoelastic (QLV) model with a stretched-exponential model of cartilage stress relaxation and examined each of these models both alone and in combination with a fluid-flow model in unconfined compression. Cartilage explants were harvested from bovine calf patellofemoral joints and equilibrated in tissue culture for 5 days before stress-relaxation testing in unconfined compression at 5% nominal strain. The stretched exponential models fit as well as the QLV models. Furthermore, the average stretched exponential relaxation time determined by this model lies within the range of experimentally measured relaxation times for extracted proteoglycan aggregates, consistent with the hypothesis that the stretched exponential model represents polymeric mechanisms of cartilage viscoelasticity.     

Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone

The molecular interaction between common polymer chains and the cell membrane is unknown. Molecular dynamics simulations offer an emerging tool to characterise the nature of the interaction between common degradable polymer chains used in biomedical applications, such as polycaprolactone, and model cell membranes. Herein we characterise with all-atomistic and coarse-grained molecular dynamics simulations the interaction between single polycaprolactone chains of varying chain lengths with a phospholipid membrane. We find that the length of the polymer chain greatly affects the nature of interaction with the membrane, as well as the membrane properties. Furthermore, we next utilise advanced sampling techniques in molecular dynamics to characterise the two-dimensional free energy surface for the interaction of varying polymer chain lengths (short, intermediate, and long) with model cell membranes. We find that the free energy minimum shifts from the membrane-water interface to the hydrophobic core of the phospholipid membrane as a function of chain length. Finally, we perform coarse-grained molecular dynamics simulations of slightly larger membranes with polymers of the same length and characterise the results as compared with all-atomistic molecular dynamics simulations. These results can be used to design polymer chain lengths and chemistries to optimise their interaction with cell membranes at the molecular level.     

Stress relaxation behaviors of articular cartilages in porcine temporomandibular joint

In this study, we tested the compressive stress relaxation behaviors of the mandibular condylar and temporal cartilages in the porcine temporomandibular joint (TMJ). The aim was to determine the quantitative and qualitative similarities and differences of compressive stress relaxation behaviors between the two cartilages. Ten porcine TMJs were used; the articular surface was divided into 5 regions: anterior, central, posterior, lateral and medial. Compressive relaxation test was carried out at a strain level of 5% in each region of the two cartilages. The stress relaxation was monitored over a period of 5 min. In all the regions of the two cartilages, the time-dependent stress relaxation curves showed a marked drop in stress within the initial 10 s, which can be fitted by a standard linear viscoelastic model. The instantaneous moduli in the temporal cartilage were dominantly larger than those in the condylar cartilage, while the condylar cartilage had slightly larger relaxation moduli than the temporal cartilage except for the medial region. The both cartilages showed the regional differences in the compressive stress relaxation behavior, and in the temporal cartilage the lateral and medial regions revealed the largest values for the instantaneous and relaxation moduli. The present results demonstrate that the viscoelastic properties of compressive stress relaxation in both cartilages are region-specific, which might have an important implication for stress distribution and transmission along with the TMJ disc.     

The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation

The relative importance of fluid-dependent and fluid-independent transient mechanical behavior in articular cartilage was examined for tensile and unconfined compression testing using a fibril reinforced model. The collagen matrix of articular cartilage was modeled as viscoelastic using a quasi-linear viscoelastic formulation with strain-dependent elastic modulus, while the proteoglycan matrix was considered as linearly elastic. The collagen viscoelastic properties were obtained by fitting experimental data from a tensile test. These properties were used to investigate unconfined compression testing, and the sensitivity of the properties was also explored. It was predicted that the stress relaxation observed in tensile tests was not caused by fluid pressurization at the macroscopic level. A multi-step tensile stress relaxation test could be approximated using a hereditary integral in which the elastic fibrillar modulus was taken to be a linear function of the fibrillar strain. Applying the same formulation to the radial fibers in unconfined compression, stress relaxation could not be simulated if fluid pressurization were absent. Collagen viscoelasticity was found to slightly weaken fluid pressurization in unconfined compression, and this effect was relatively more significant at moderate strain rates. Therefore, collagen viscoelasticity appears to play an import role in articular cartilage in tensile testing, while fluid pressurization dominates the transient mechanical behavior in compression. Collagen viscoelasticity plays a minor role in the mechanical response of cartilage in unconfined compression if significant fluid flow is present.     

A continuum theory and an experiment for the ion-induced swelling behavior of articular cartilage

Swelling of normal bovine articular cartilage equilibrated in NaCl solutions was dimensionally measured in thin strips of tissue. The ion-induced strains show that free swelling of articular cartilage is anisotropic and inhomogeneous. For the molar concentrations used, contraction increased linearly with concentration, defining a "coefficient of chemical contraction" (alpha c). Isometrically constrained specimens registered a rise in tensile force followed by stress relaxation. An extension of the biphasic theory incorporating this ion-induced strain is proposed. This theory can describe the equilibrium anisotropic swelling behavior of cartilage and explain the transient force history observed in the isometric experiment.     

Nuclear magnetic relaxation dispersion study of the dynamics in solid homopolypeptides

The (1)H nuclear magnetic relaxation dispersion profiles were measured from 10 kHz to 30 MHz as a function of temperature for polyglycine, polyalanine, polyvaline, and polyphenylalanine to examine the contributions of different side chain motions to the polypeptide proton relaxation rate constants. The spin-fracton theory for (1)H relaxation is modified to account for high frequency motions of side chains that are dynamically connected to the linear polymer backbone. The (1)H relaxation is dominated by propagation of rare disturbances along the backbone of the polymer. The side-chain dynamics cause an off-set in the field dependence of the (1)H spin-lattice relaxation rate constants which obey a power law in the Larmor frequency in the limit of low and high magnetic field strength.     

Determinants of Fluidlike Behavior and Effective Viscosity in Cross-Linked Actin Networks

The actin cortex has a well-documented ability to rapidly remodel and flow while maintaining long-range connectivity, but how this is achieved remains poorly understood. Here, we use computer simulations to explore how stress relaxation in cross-linked actin networks subjected to extensional stress depends on the interplay between network architecture and turnover. We characterize a regime in which a network response is nonaffine and stress relaxation is governed by the continuous dissipation of elastic energy via cyclic formation, elongation, and turnover of tension-bearing elements. Within this regime, for a wide range of network parameters, we observe a constant deformation (creep) rate that is linearly proportional to the rate of filament turnover, leading to a constant effective viscosity that is inversely proportional to turnover rate. Significantly, we observe a biphasic dependence of the creep rate on applied stress: below a critical stress threshold, the creep rate increases linearly with applied stress; above that threshold, the creep rate becomes independent of applied stress. We show that this biphasic stress dependence can be understood in terms of the nonlinear force-extension behavior of individual force-transmitting network elements. These results have important implications for understanding the origins and control of viscous flows both in the cortex of living cells and in other polymer networks.     

Determinants of Fluidlike Behavior and Effective Viscosity in Cross-Linked Actin Networks

The actin cortex has a well-documented ability to rapidly remodel and flow while maintaining long-range connectivity, but how this is achieved remains poorly understood. Here, we use computer simulations to explore how stress relaxation in cross-linked actin networks subjected to extensional stress depends on the interplay between network architecture and turnover. We characterize a regime in which a network response is nonaffine and stress relaxation is governed by the continuous dissipation of elastic energy via cyclic formation, elongation, and turnover of tension-bearing elements. Within this regime, for a wide range of network parameters, we observe a constant deformation (creep) rate that is linearly proportional to the rate of filament turnover, leading to a constant effective viscosity that is inversely proportional to turnover rate. Significantly, we observe a biphasic dependence of the creep rate on applied stress: below a critical stress threshold, the creep rate increases linearly with applied stress; above that threshold, the creep rate becomes independent of applied stress. We show that this biphasic stress dependence can be understood in terms of the nonlinear force-extension behavior of individual force-transmitting network elements. These results have important implications for understanding the origins and control of viscous flows both in the cortex of living cells and in other polymer networks.     

Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method

To relate the subcellular molecular events to organ level physiology in heart, we have developed a three-dimensional finite-element-based simulation program incorporating the cellular mechanisms of excitation-contraction coupling and its propagation, and simulated the fluid-structure interaction involved in the contraction and relaxation of the human left ventricle. The FitzHugh-Nagumo model and four-state model representing the cross-bridge kinetics were adopted for cellular model. Both ventricular wall and blood in the cavity were modeled by finite element mesh. An arbitrary Lagrangian Eulerian finite element method with automatic mesh updating has been formulated for large domain changes, and a strong coupling strategy has been taken. Using electrical analog of pulmonary circulation and left atrium as a preload and the windkessel model as an afterload, dynamics of ventricular filling as well as ejection was simulated. We successfully reproduced the biphasic filling flow consisting of early rapid filling and atrial contraction similar to that reported in clinical observation. Furthermore, fluid-structure analysis enabled us to analyze the wave propagation velocity of filling flow. This simulator can be a powerful tool for establishing a link between molecular abnormality and the clinical disorder at the macroscopic level.     

Enzymatic digestion of articular cartilage results in viscoelasticity changes that are consistent with polymer dynamics mechanisms

Background  

Cartilage degeneration via osteoarthritis affects millions of elderly people worldwide, yet the specific contributions of matrix biopolymers toward cartilage viscoelastic properties remain unknown despite 30 years of research. Polymer dynamics theory may enable such an understanding, and predicts that cartilage stress-relaxation will proceed faster when the average polymer length is shortened.     

Calculation and Visualization of Atomistic Mechanical Stresses in Nanomaterials and Biomolecules

Many biomolecules have machine-like functions, and accordingly are discussed in terms of mechanical properties like force and motion. However, the concept of stress, a mechanical property that is of fundamental importance in the study of macroscopic mechanics, is not commonly applied in the biomolecular context. We anticipate that microscopical stress analyses of biomolecules and nanomaterials will provide useful mechanistic insights and help guide molecular design. To enable such applications, we have developed Calculator of Atomistic Mechanical Stress (CAMS), an open-source software package for computing atomic resolution stresses from molecular dynamics (MD) simulations. The software also enables decomposition of stress into contributions from bonded, nonbonded and Generalized Born potential terms. CAMS reads GROMACS topology and trajectory files, which are easily generated from AMBER files as well; and time-varying stresses may be animated and visualized in the VMD viewer. Here, we review relevant theory and present illustrative applications.     

Viscoelastic properties of proteoglycan subunits and aggregates in varying solution concentrations

Using a cone-on-plate mechanical spectrometer, we have measured the linear and non-linear rheological properties of cartilage proteoglycan solutions at concentrations similar to those found in situ. Solutions of bovine nasal cartilage proteoglycan subunits (22S) and aggregates (79S) were studied at concentrations ranging from 10 to 50 mg ml-1. We determined: (1) the complex viscoelastic shear modulus G (omega) under small amplitude (0.02 radians) oscillatory excitation at frequencies (omega) ranging from 1.0 to 20.0 Hz, (2) the non-linear shear rate (gamma) dependent apparent viscosity napp (gamma) in continuous shear, and (3) the non-linear shear rate dependent primary normal stress difference sigma 1 (gamma) in continuous shear. Both the apparent viscosity and normal stress difference were measured over four decades of shear rates ranging from 0.25 to 250 s-1. Analysis of the experimental results were performed using a variety of materially objective non-linear viscoelastic constitutive laws. We found that the non-linear, four-coefficient Oldroyd rate-type model was most effective for describing the measured flow characteristics of proteoglycan subunit and aggregate solutions. Values of the relaxation time lambda 1, retardation time lambda 2, zero shear viscosity no, and nonlinear viscosity parameter muo were computed for the aggregate and subunit solutions at all of the solute concentrations used. The four independent material coefficients showed marked dependence on the two different molecular conformations, i.e. aggregate or subunit, of proteoglycans in solution.     

Quantitative study of polymer conformation and dynamics by single-particle tracking

We present a new method for analyzing the dynamics of conformational fluctuations of individual flexible polymer molecules. In single-particle tracking (SPT), one end of the polymer molecule is tethered to an immobile substratum. A microsphere attached to the other end serves as an optical marker. The conformational fluctuations of the polymer molecule can be measured by optical microscopy via the motion of the microsphere. The bead-and-spring theory for polymer dynamics is further developed to account for the microsphere, and together the measurement and the theory yield quantitative information about molecular conformations and dynamics under nonperturbing conditions. Applying the method to measurements carried out on DNA molecules provides information complementary to recent studies of single DNA molecules under extensional force. Combining high precision measurements with the theoretical analysis presented here creates a powerful tool for studying conformational dynamics of biological and synthetic macromolecules at the single-molecule level.     

Dynamic response of immature bovine articular cartilage in tension and compression, and nonlinear viscoelastic modeling of the tensile response

Very limited information is currently available on the constitutive modeling of the tensile response of articular cartilage and its dynamic modulus at various loading frequencies. The objectives of this study were to (1) formulate and experimentally validate a constitutive model for the intrinsic viscoelasticity of cartilage in tension, (2) confirm the hypothesis that energy dissipation in tension is less than in compression at various loading frequencies, and (3) test the hypothesis that the dynamic modulus of cartilage in unconfined compression is dependent upon the dynamic tensile modulus. Experiment 1: Immature bovine articular cartilage samples were tested in tensile stress relaxation and cyclical loading. A proposed reduced relaxation function was fitted to the stress-relaxation response and the resulting material coefficients were used to predict the response to cyclical loading. Adjoining tissue samples were tested in unconfined compression stress relaxation and cyclical loading. Experiment 2: Tensile stress relaxation experiments were performed at varying strains to explore the strain-dependence of the viscoelastic response. The proposed relaxation function successfully fit the experimental tensile stress-relaxation response, with R2 = 0.970+/-0.019 at 1% strain and R2 = 0.992+/-0.007 at 2% strain. The predicted cyclical response agreed well with experimental measurements, particularly for the dynamic modulus at various frequencies. The relaxation function, measured from 2% to 10% strain, was found to be strain dependent, indicating that cartilage is nonlinearly viscoelastic in tension. Under dynamic loading, the tensile modulus at 10 Hz was approximately 2.3 times the value of the equilibrium modulus. In contrast, the dynamic stiffening ratio in unconfined compression was approximately 24. The energy dissipation in tension was found to be significantly smaller than in compression (dynamic phase angle of 16.7+/-7.4 deg versus 53.5+/-12.8 deg at 10(-3) Hz). A very strong linear correlation was observed between the dynamic tensile and dynamic compressive moduli at various frequencies (R2 = 0.908+/-0.100). The tensile response of cartilage is nonlinearly viscoelastic, with the relaxation response varying with strain. A proposed constitutive relation for the tensile response was successfully validated. The frequency response of the tensile modulus of cartilage was reported for the first time. Results emphasize that fluid-flow dependent viscoelasticity dominates the compressive response of cartilage, whereas intrinsic solid matrix viscoelasticity dominates the tensile response. Yet the dynamic compressive modulus of cartilage is critically dependent upon elevated values of the dynamic tensile modulus.     

Quasielastic light-scattering study of polyadenylic acid solutions. II

Quasielastic light scattering is used to study saline solutions of polyadenylic acid with varying polymer concentrations and molecular masses. These experiments clearly show the existence of two relaxation times. For dilute solutions, when the chains are mutually independent, the fast mode is due to the free diffusion of the polymer chains. For concentrations above the overlap concentration C*, the fast mode is due to the propagation of collective excitations of the pseudolattice of polymer chains. The slow modes are observed when the polymer concentration is in the vicinity of the overlap concentration C*. A series of experiments shows that both their relaxation time and amplitude depend only on the polymer concentration and not on the polymer molecular mass. This result rules out any previous explanation based on individual chain motion. Furthermore, since the amplitudes depend on the time elapsed from the preparation of the solution, the slow modes are due to the diffusion of concentration inhomogeneities in the pseudolattice.