Impact of matrix metalloproteinases on inhibition of mineralization by fetuin.

BACKGROUND AND OBJECTIVE: Human subjects affected by inflammatory diseases, such as periodontitis, may be at increased risk for the development of cardiovascular diseases and calcification of atheromas; however, the potential mechanisms have not been defined. Alpha-2-Heremans Schmid glycoprotein (fetuin A) is an abundant serum glycoprotein of ~49 kDa that inhibits ectopic arterial calcification. We examined whether matrix metalloproteinases (MMPs), which are increased in inflammatory diseases, including periodontitis, bind and degrade fetuin and alter its ability to inhibit calcification in vitro. MATERIAL AND METHODS: Binding and cleavage of fetuin by MMPs were assessed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, in-silico analyses and mass spectrometry. The effects of intact and MMP-degraded human fetuin on mineralization were measured in a cell-free assay. RESULTS: From in-silico analyses and literature review, we found that only MMP-3 (stromelysin) and MMP-7 (matrilysin) were predicted to cleave human fetuin at levels that were physiologically relevant. In-vitro assays showed that MMP-7, and, to a lesser extent, MMP-3, degraded human fetuin in a time- and dose-dependent manner. Fetuin peptides generated by MMP-7 cleavage were identified and sequenced by mass spectrometry; novel cleavage sites were found. Hydroxyapatite mineralization in vitro was strongly inhibited by fetuin (> 1 µm), as was MMP-3-cleaved fetuin, while MMP-7-cleaved fetuin was threefold less effective in blocking mineralization. CONCLUSION: MMP-7 and, to a lesser extent, MMP-3, affect the ability of fetuin to inhibit the formation of hydroxyapatite in vitro. These data suggest that the MMPs increased in inflammatory diseases, such as periodontitis, could affect regulation of mineralization and potentially enhance the risk of calcified atheroma formation.

Delivery of gelfoam-enabled cells and vectors into the pericardial space using a percutaneous approach in a porcine model.

Intrapericardial drug delivery is a promising procedure, with the ability to localize therapeutics with the heart. Gelfoam particles are nontoxic, inexpensive, nonimmunogenic and biodegradable compounds that can be used to deliver therapeutic agents. We developed a new percutaneous approach method for intrapericardial injection, puncturing the pericardial sac safely under fluoroscopy and intravascular ultrasound (IVUS) guidance. In a porcine model of myocardial infarction (MI), we deployed gelfoam particles carrying either (a) autologous mesenchymal stem cells (MSCs) or (b) an adenovirus encoding enhanced green fluorescent protein (eGFP) 48 h post-MI. The presence of MSCs and viral infection at the infarct zone was confirmed by immunoflourescence and PCR. Puncture was performed successfully in 16 animals. Using IVUS, we successfully determined the size of the pericardial space before the puncture, and safely accessed that space in setting of pericardial effusion and also adhesions induced by the MI. Intrapericardial injection of gelfoam was safe and reliable. Presence of the MSCs and eGFP expression from adenovirus in the myocardium were confirmed after delivery. Our novel percutaneous approach to deliver (stem-) cells or adenovirus was safe and efficient in this pre-clinical model. IVUS-guided delivery is a minimally invasive procedure that seems to be a promising new strategy to deliver therapeutic agents locally to the heart.

Mechanical properties of bovine articular cartilage under microscale indentation loading from atomic force microscopy.

Atomic force microscopy (AFM) techniques have been increasingly used for investigating the mechanical properties of articular cartilage. According to the previous studies reporting the microscale Young's modulus under AFM indentation tests, the Hertz contact model has been employed with a sharp conical tip indenter. However, the non-linear microscale behaviour of articular cartilage could not be resolved by the standardized Hertz analysis using small and sharp atomic force microscope tips. Therefore, the objective of this study was to evaluate the microscale Young's modulus of articular cartilage more accurately through a non-Hertzian approach with a spherical tip of 5 microm diameter, and to characterize its microscale mechanical behaviour. This methodology adopted in the present study was proved by the consistent values between the microscale (2 per cent, about 9.3 kPa; 3 per cent, about 17.5kPa) and macroscale (2 per cent, about 8.3kPa; 3 per cent, about 18.3kPa) Young's moduli for 2 per cent and 3 per cent agarose gel (n = 100). Therefore, the microscale Young's modulus evaluated in this study is representative of more accurate measurements of cartilage stiffness at the 600 nm deformation level and corresponds to approximately 30.9 kPa (n = 100). Furthermore, on this level of the microscale deformation, articular cartilage showed depth-dependent and frequency-independent behaviour under AFM indentation loading. These findings reveal the microscale mechanical behaviour of articular cartilage more accurately and can be employed further to design microscale structures of chondrocyte-seeded scaffolds and tissue-engineered cartilage by evaluating their microscale properties.

A theoretical analysis of water transport through chondrocytes.

Because of the avascular nature of adult cartilage, nutrients and waste products are transported to and from the chondrocytes by diffusion and convection through the extracellular matrix. The convective interstitial fluid flow within and around chondrocytes is poorly understood. This theoretical study demonstrates that the incorporation of a semi-permeable membrane when modeling the chondrocyte leads to the following findings: under mechanical loading of an isolated chondrocyte the intracellular fluid pressure is on the order of tens of Pascals and the transmembrane fluid outflow, on the order of picometers per second, takes several days to subside; consequently, the chondrocyte behaves practically as an incompressible solid whenever the loading duration is on the order of minutes or hours. When embedded in its extracellular matrix (ECM), the chondrocyte response is substantially different. Mechanical loading of the tissue leads to a fluid pressure difference between intracellular and extracellular compartments on the order of tens of kilopascals and the transmembrane outflow, on the order of a nanometer per second, subsides in about 1 h. The volume of the chondrocyte decreases concomitantly with that of the ECM. The interstitial fluid flow in the extracellular matrix is directed around the cell, with peak values on the order of tens of nanometers per second. The viscous fluid shear stress acting on the cell surface is several orders of magnitude smaller than the solid matrix shear stresses resulting from the ECM deformation. These results provide new insight toward our understanding of water transport in chondrocytes.

Removal of the superficial zone of bovine articular cartilage does not increase its frictional coefficient.

OBJECTIVE: To investigate the role of the superficial zone in regulating the frictional response of articular cartilage. This zone contains the superficial protein (SZP), a proteoglycan synthesized exclusively by superficial zone chondrocytes and implicated in reducing the friction coefficient of cartilage. DESIGN: Unconfined compression creep tests with sliding of cartilage against glass in saline were carried out on fresh bovine cylindrical plugs (slashed circle Ø6 mm, n=35) obtained from 16 bovine shoulder joints (ages 1-3 months). In the first two experiments, friction tests were carried out before and after removal of the superficial zone ( approximately 100 microm), in a control and treatment group, using two different applied load magnitudes (4.4 N and 22.2 N). In the third experiment, friction tests were conducted on intact surfaces and the corresponding microtomed deep zone of the same specimen. RESULTS: In all tests the friction coefficient exhibited a transient response, increasing from a minimum value (mu(min)) to a near-equilibrium final value (micro(eq)). No statistical change (P>0.5) was found in micro(min) before and after removal of the superficial zone in both experiments 1 and 2. However, micro(eq) was observed to decrease significantly (P<0.001) after removal of the surface zone. Results from the third experiment confirm that micro(eq) is even lower at the deep zone. Surface roughness measurements with atomic force microscopy (AFM) revealed an increase in surface roughness after microtoming. Immunohistochemical staining confirmed the presence of SZP in intact specimens and its removal in microtomed specimens. CONCLUSIONS: The topmost ( approximately 100 microm) superficial zone of articular cartilage does not have special properties which enhances its frictional response.

Relating myocardial laminar architecture to shear strain and muscle fiber orientation.

Cardiac myofibers are organized into laminar sheets about four cells thick. Recently, it has been suggested that these layers coincide with the plane of maximum shear during systole. In general, there are two such planes, which are oriented at +/-45 degrees to the main principal strain axes. These planes do not necessarily contain the fiber axis. In the present study, we explicitly added the constraint that the sheet planes should also contain the muscle fiber axis. In a mathematical analysis of previously measured three-dimensional transmural systolic strain distributions in six dogs, we computed the planes of maximum shear, adding the latter constraint by using the also-measured muscle fiber axis. Generally, for such planes two solutions were found, suggesting that two populations of sheet orientation may exist. The angles at which the predicted sheets intersected transmural tissue slices, cut along left ventricular short- or long-axis planes, were strikingly similar to experimentally measured values. In conclusion, sheets coincide with planes of maximum systolic shear subject to the constraint that the muscle fiber axis is contained in this plane. Sheet orientation is not a unique function of the transmural location but occurs in two distinct populations.

Mechanism underlying mechanical dysfunction in the border zone of left ventricular aneurysm: a finite element model study.

BACKGROUND: The global left ventricular dysfunction characteristic of left ventricular aneurysm is associated with muscle fiber stretching in the adjacent noninfarcted (border zone) region during isovolumic systole. The mechanism of this regional dysfunction is poorly understood. METHODS: An anteroapical transmural myocardial infarct was created by coronary arterial ligation in an adult Dorset sheep and was allowed to mature into left ventricular aneurysm for 10 weeks. The animal was imaged subsequently using magnetic resonance imaging with simultaneous recording of intraventricular pressures. A realistic mathematical model of the three-dimensional ovine left ventricle with an anteroapical aneurysm was constructed from multiple short-axis and long-axis magnetic resonance imaging slices at the beginning of diastolic filling. RESULTS: Three model simulations are presented: (1) normal border zone contractility and normal aneurysmal material properties; (2) greatly reduced border zone contractility (by 50%) and normal aneurysmal material properties; and (3) greatly reduced border zone contractility (by 50%) and stiffened aneurysmal material properties (by 1000%). Only the latter two simulations were able to reproduce experimentally observed stretching of border zone fibers during isovolumic systole. CONCLUSIONS: The mechanism underlying mechanical dysfunction in the border zone region of left ventricular aneurysm is primarily the result of myocardial contractile dysfunction rather than increased wall stress in this region.

Analysis of indentation: implications for measuring mechanical properties with atomic force microscopy.

Indentation using the atomic force microscope (AFM) has potential to measure detailed micromechanical properties of soft biological samples. However, interpretation of the results is complicated by the tapered shape of the AFM probe tip, and its small size relative to the depth of indentation. Finite element models (FEMs) were used to examine effects of indentation depth, tip geometry, and material nonlinearity and heterogeneity on the finite indentation response. Widely applied infinitesimal strain models agreed with FEM results for linear elastic materials, but yielded substantial errors in the estimated properties for nonlinear elastic materials. By accounting for the indenter geometry to compute an apparent elastic modulus as a function of indentation depth, nonlinearity and heterogeneity of material properties may be identified. Furthermore, combined finite indentation and biaxial stretch may reveal the specific functional form of the constitutive law--a requirement for quantitative estimates of material constants to be extracted from AFM indentation data.

Laminar fiber architecture and three-dimensional systolic mechanics in canine ventricular myocardium.

Previous studies suggest that the laminar architecture of left ventricular myocardium may be critical for normal ventricular mechanics. However, systolic three-dimensional deformation of the laminae has never been measured. Therefore, end-systolic finite strains relative to end diastole, from biplane radiography of transmural markers near the apex and base of the anesthetized open-chest canine anterior left ventricular free wall (n = 6), were referred to three-dimensional laminar microstructural axes reconstructed from histology. Whereas fiber shortening was uniform [-0.07 +/- 0.04 (SD)], radial wall thickening increased from base (0. 10 +/- 0.09) to apex (0.14 +/- 0.13). Extension of the laminae transverse to the muscle fibers also increased from base (0.08 +/- 0. 07) to apex (0.11 +/- 0.08), and interlaminar shear changed sign [0. 05 +/- 0.07 (base) and -0.07 +/- 0.09 (apex)], reflecting variations in laminar architecture. Nevertheless, the apex and base were similar in that at each site laminar extension and shear contributed approximately 60 and 40%, respectively, of mean transmural thickening. Kinematic considerations suggest that these dual wall-thickening mechanisms may have distinct ultrastructural origins.

Three-dimensional residual strain in midanterior canine left ventricle.

All previous studies of residual strain in the ventricular wall have been based on one- or two-dimensional measurements. Transmural distributions of three-dimensional (3-D) residual strains were measured by biplane radiography of columns of lead beads implanted in the midanterior free wall of the canine left ventricle (LV). 3-D bead coordinates were reconstructed with the isolated arrested LV in the zero-pressure state and again after local residual stress had been relieved by excising a transmural block of tissue. Nonhomogeneous 3-D residual strains were computed by finite element analysis. Mean +/- SD (n = 8) circumferential residual strain indicated that the intact unloaded myocardium was prestretched at the epicardium (0.07 +/- 0.06) and compressed in the subendocardium (-0.04 +/- 0.05). Small but significant longitudinal shortening and torsional shear residual strains were also measured. Residual fiber strain was tensile at the epicardium (0.05 +/- 0.06) and compressive in the subendocardium (-0.01 +/- 0.04), with residual extension and shortening, respectively, along structural axes parallel and perpendicular to the laminar myocardial sheets. Relatively small residual shear strains with respect to the myofiber sheets suggest that prestretching in the plane of the myocardial laminae may be a primary mechanism of residual stress in the LV.

A three-dimensional finite element method for large elastic deformations of ventricular myocardium: I--Cylindrical and spherical polar coordinates.

A three-dimensional Galerkin finite element method was developed for large deformations of ventricular myocardium and other incompressible, nonlinear elastic, anisotropic materials. Cylindrical and spherical elements were used to solve axisymmetric problems with r.m.s. errors typically less than 2 percent. Isochoric interpolation and pressure boundary constraint equations enhanced low-order curvilinear elements under special circumstances (69 percent savings in degrees of freedom, 78 percent savings in solution time for inflation of a thick-walled cylinder). Generalized tensor products of linear Lagrange and cubic Hermite polynomials permitted custom elements with improved performance, including 52 percent savings in degrees of freedom and 66 percent savings in solution time for compression of a circular disk. Such computational efficiencies become significant for large scale problems such as modeling the heart.

A three-dimensional finite element method for large elastic deformations of ventricular myocardium: II--Prolate spheroidal coordinates.

A three-dimensional finite element method for nonlinear finite elasticity is presented using prolate spheroidal coordinates. For a thick-walled ellipsoidal model of passive anisotropic left ventricle, a high-order (cubic Hermite) mesh with 3 elements gave accurate continuous stresses and strains, with a 69 percent savings in degrees of freedom (dof) versus a 70-element standard low-order model. A custom mixed-order model offered 55 percent savings in dof and 39 percent savings in solution time compared with the low-order model. A nonsymmetric 3D model of the passive canine LV was solved using 16 high-order elements. Continuous nonhomogeneous stresses and strains were obtained within 1 hour on a laboratory workstation, with an estimated solution time of less than 4 hours to model end-systole. This method represents the first practical opportunity to solve large-scale anatomically detailed models for cardiac stress analysis.

Finite element stress analysis of left ventricular mechanics in the beating dog heart.

A three-dimensional finite element model was used to explore whether or not transmural distributions of end-diastolic and end-systolic fiber stress are uniform from the apex to the base of the canine left ventricular wall. An elastance model for active fiber stress was incorporated in an axisymmetric model that accurately represented the geometry and fiber angle distribution of the anterior free wall. The nonlinear constitutive equation for the resting myocardium was transversely isotropic with respect to the local fiber axis. Transmural distributions of end-diastolic fiber stress became increasingly nonuniform from midventricle toward the apex or the base. At a typical diastolic left ventricular pressure (1 kPa), the differences between largest and smallest fiber stresses were only 0.5 kPa near midventricle, compared with 4.6 kPa at the apex, and 3.3 kPa at the base. Transmural fiber stress differences at end-systole (14 kPa) were relatively small in regions from the base to the midventricle (13-22 kPa), but were larger between midventricle and the apex (30-43 kPa). All six three-dimensional end-diastolic strain components were within or very close to one standard deviation of published measurements through the midanterior left ventricular free wall of the passive canine heart [Omens et al., Am. J. Physiol. 261, H918-H928 (1991)]. End-systolic in-plane normal and shear strains also agreed closely with published experimental measurements in the beating dog heart [Waldman et al., Circ. Res. 63, 550-562 (1988)]. The results indicate that, unlike in the midventricle region that has been studied most fully, there may be significant regional nonhomogeneity of fiber stress in the normal left ventricle associated with regional variations in shape and fiber angle.

Physiological interpretations based on lumped element models fit to respiratory impedance data: use of forward-inverse modeling.

Respiratory impedance (Zrs) data at lower (less than 4 Hz) and higher (greater than 32 Hz) frequencies require more complicated inverse models than the standard series combination of a respiratory resistance, inertance, and compliance. In this paper, a forward-inverse modeling approach was used to provide insight on how the parameters in these more complicated inverse models reflect the true physiological system. Forward models are set up to incorporate explicit physiological and anatomical detail. Simulated forward data are then fit with identifiable inverse models and the parameter estimates related to the known detail in the forward model. It is shown that inverse fitting of low frequency data alone will not allow a distinction between frequency dependence due to airway inhomogeneities and frequency dependence due to tissue viscoelasticity. With higher frequency data, a forward model based on an asymmetric branching airways network was used to simulate Zrs from 0.1-128 Hz with increasing amounts of nonuniform peripheral airway obstruction. Here, inverse modeling is more amenable to sensibly separating estimates of airway and tissue properties. A key result, however, is that changes in the tissue parameters of an inverse model (which provides an excellent fit to Zrs data) will appropriately occur in response to inhomogeneous alterations in airway diameters only. The apparent altered tissue properties reflect the decreased communication of some tissue segments with the airway opening and not an explicit change at the tissue level. These phenomena present a substantial problem for the inverse modeler. Finally, inverse model fitting of low and high frequency Zrs data simultaneously with a single model is not helpful for extracting additional physiological detail. Instead, separate models should be applied to each frequency range.