Biomedical subjects

C C van Donkelaar

Publications and source records attributed to C C van Donkelaar.

At least 19 recordsLinked to original sources

PubMed2006

Prediction of collagen orientation in articular cartilage by a collagen remodeling algorithm.

OBJECTIVE: Tissue engineering is a promising method to treat damaged cartilage. So far it has not been possible to create tissue-engineered cartilage with an appropriate structural organization. It is envisaged that cartilage tissue engineering will significantly benefit from knowledge of how the collagen fiber orientation is directed by mechanical conditions. The goal of the present study is to evaluate whether a collagen remodeling algorithm based on mechanical loading can be corroborated by the collagen orientation in healthy cartilage. METHODS: According to the remodeling algorithm, collagen fibrils align with a preferred fibril direction, situated between the positive principal strain directions. The remodeling algorithm was implemented in an axisymmetric finite element model of the knee joint. Loading as a result of typical daily activities was represented in three different phases: rest, standing and gait. RESULTS: In the center of the tibial plateau the collagen fibrils run perpendicular to the subchondral bone. Just below the articular surface they bend over to merge with the articular surface. Halfway between the center and the periphery, the collagen fibrils bend over earlier, resulting in a thicker superficial and transitional zones. Near the periphery fibrils in the deep zone run perpendicular to the articular surface and slowly bend over to angles of -45 degrees and +45 degrees with the articular surface. CONCLUSION: The collagen structure as predicted with the collagen remodeling algorithm corresponds very well with the collagen structure in healthy knee joints. This remodeling algorithm is therefore considered to be a valuable tool for developing loading protocols for tissue engineering of articular cartilage.

Algorithms↗

PubMed2006

Depth-dependent compressive equilibrium properties of articular cartilage explained by its composition.

For this study, we hypothesized that the depth-dependent compressive equilibrium properties of articular cartilage are the inherent consequence of its depth-dependent composition, and not the result of depth-dependent material properties. To test this hypothesis, our recently developed fibril-reinforced poroviscoelastic swelling model was expanded to include the influence of intra- and extra-fibrillar water content, and the influence of the solid fraction on the compressive properties of the tissue. With this model, the depth-dependent compressive equilibrium properties of articular cartilage were determined, and compared with experimental data from the literature. The typical depth-dependent behavior of articular cartilage was predicted by this model. The effective aggregate modulus was highly strain-dependent. It decreased with increasing strain for low strains, and increases with increasing strain for high strains. This effect was more pronounced with increasing distance from the articular surface. The main insight from this study is that the depth-dependent material behavior of articular cartilage can be obtained from its depth-dependent composition only. This eliminates the need for the assumption that the material properties of the different constituents themselves vary with depth. Such insights are important for understanding cartilage mechanical behavior, cartilage damage mechanisms and tissue engineering studies.

Cartilage, Articular↗

PubMed2006

The PTHrP-Ihh feedback loop in the embryonic growth plate allows PTHrP to control hypertrophy and Ihh to regulate proliferation.

Growth plate and long bone development is governed by biochemical signaling pathways of which the PTHrP-Ihh system is the best known. Other factors, such as BMPs, FGFs and mechanical loading, may interact with this system. This study aims at elucidating the relative importance of PTHrP and Ihh for controlling proliferation, and hypertrophy in fetal growth plate cartilage. We assessed the question why reduced Ihh expression leads to more pronounced effects on the number of non-hypertrophic cells and total bone formation, compared to PTHrP down-regulation. Using few basic equations, constituted from literature data, this paper shows how the PTHrP-Ihh feedback system can control different aspects of tissue differentiation at distinct locations. In particular, it is shown that (mechanical or biochemical) perturbations will affect proliferation via Ihh-related parameters, whereas changes in PTHrP-related parameters selectively interact with hypertrophy. This is contra-intuitive, since PTHrP acts to keep cells proliferating. In this context, the critical PTHrP level for keeping cells proliferating has been reconsidered. In addition, an explanation is provided for the aforementioned difference in effect between reduced Ihh and PTHrP expression.

Animals↗

PubMed2006

A composition-based cartilage model for the assessment of compositional changes during cartilage damage and adaptation.

OBJECTIVE: The composition of articular cartilage changes with progression of osteoarthritis. Since compositional changes are associated with changes in the mechanical properties of the tissue, they are relevant for understanding how mechanical loading induces progression. The objective of this study is to present a computational model of articular cartilage which enables to study the interaction between composition and mechanics. METHODS: Our previously developed fibril-reinforced poroviscoelastic swelling model for articular cartilage was combined with our tissue composition-based model. In the combined model both the depth- and strain-dependencies of the permeability are governed by tissue composition. All local mechanical properties in the combined model are directly related to the local composition of the tissue, i.e., to the local amounts of proteoglycans and collagens and to tissue anisotropy. RESULTS: Solely based on the composition of the cartilage, we were able to predict the equilibrium and transient response of articular cartilage during confined compression, unconfined compression, indentation and two different 1D-swelling tests, simultaneously. CONCLUSION: Since both the static and the time-dependent mechanical properties have now become fully dependent on tissue composition, the model allows assessing the mechanical consequences of compositional changes seen during osteoarthritis without further assumptions. This is a major step forward in quantitative evaluations of osteoarthritis progression.

Adaptation, Physiological↗

PubMed2005

Can the growth factors PTHrP, Ihh and VEGF, together regulate the development of a long bone?

Endochondral ossification is the process of differentiation of cartilaginous into osseous tissue. Parathyroid hormone related protein (PTHrP), Indian hedgehog (Ihh) and vascular endothelial growth factor (VEGF), which are synthesized in different zones of the growth plate, were found to have crucial roles in regulating endochondral ossification. The aim of this study was to evaluate whether the three growth factors PTHrP, Ihh and VEGF, together, could regulate longitudinal growth in a normal human, fetal femur. For this purpose, a one-dimensional finite element (FE) model, incorporating growth factor signaling, was developed of the human, distal, femoral growth plate. It included growth factor synthesis in the relevant zones, their transport and degradation and their effects. Simulations ran from initial hypertrophy in the center of the bone until secondary ossification starts at approximately 3.5 months postnatal. For clarity, we emphasize that no mechanical stresses were considered. The FE model showed a stable growth plate in which the bone growth rate was constant and the number of cells per zone oscillated around an equilibrium. Simulations incorporating increased and decreased PTHrP and Ihh synthesis rates resulted, respectively, in more and less cells per zone and in increased and decreased bone growth rates. The FE model correctly reflected the development of a growth plate and the rate of bone growth in the femur. Simulations incorporating increased and decreased PTHrP and Ihh synthesis rates reflected growth plate pathologies and growth plates in PTHrP-/- and Ihh-/- mice. The three growth factors, PTHrP, Ihh and VEGF, could potentially together regulate tissue differentiation.

Animals↗

PubMed2005

Determining diffusion coefficients in inhomogeneous tissues using fluorescence recovery after photobleaching.

Diffusion plays an important role in the transport of nutrients and signaling molecules in cartilaginous tissues. Diffusion coefficients can be measured by fluorescence recovery after photobleaching (FRAP). Available methods to analyze FRAP data, however, assume homogeneity in the environment of the bleached area and neglect geometrical restrictions to diffusion. Hence, diffusion coefficients in inhomogeneous materials, such as most biological tissues, cannot be assessed accurately. In this study, a new method for analyzing data from FRAP measurements has been developed, which is applicable to inhomogeneous tissues. It is based on a fitting procedure of the intensity recovery after photobleaching with a two-dimensional finite element analysis, which includes Fick's law for diffusion. The finite element analysis can account for distinctive diffusivity in predefined zones, which allows determining diffusion coefficients in inhomogeneous samples. The method is validated theoretically and experimentally in both homogeneous and inhomogeneous tissues and subsequently applied to the proliferation zone of the growth plate. Finally, the importance of accounting for inhomogeneities, for appropriate assessment of diffusivity in inhomogeneous tissues, is illustrated.

Algorithms↗

PubMed2005

A fibril-reinforced poroviscoelastic swelling model for articular cartilage.

From a mechanical point of view, the most relevant components of articular cartilage are the tight and highly organized collagen network together with the charged proteoglycans. Due to the fixed charges of the proteoglycans, the cation concentration inside the tissue is higher than in the surrounding synovial fluid. This excess of ion particles leads to an osmotic pressure difference, which causes swelling of the tissue. The fibrillar collagen network resists straining and swelling pressures. This combination makes cartilage a unique, highly hydrated and pressurized tissue, enforced with a strained collagen network. Many theories to explain articular cartilage behavior under loading, expressed in computational models that either include the swelling behavior or the properties of the anisotropic collagen structure, can be found in the literature. The most common tests used to determine the mechanical quality of articular cartilage are those of confined compression, unconfined compression, indentation and swelling. All theories currently available in the literature can explain the cartilage response occurring in some of the above tests, but none of them can explain these for all of the tests. We hypothesized that a model including simultaneous mathematical descriptions of (1) the swelling properties due to the fixed-change densities of the proteoglycans and (2) the anisotropic viscoelastic collagen structure, can explain all these test simultaneously. To study this hypothesis we extended our fibril-reinforced poroviscoelastic finite element model with our biphasic swelling model. We have shown that the newly developed fibril-reinforced poroviscoelastic swelling (FPVES) model for articular cartilage can simultaneously account for the reaction force during swelling, confined compression, indentation and unconfined compression as well as the lateral deformation during unconfined compression. Using this theory it is possible to analyze the link between the collagen network and the swelling properties of articular cartilage.

Animals↗

PubMed2005

The role of computational models in the search for the mechanical behavior and damage mechanisms of articular cartilage.

Articular cartilage plays a vital role in the function of diarthrodial joints. Due to osteoarthritis degeneration of articular cartilage occurs. The initial event that triggers the pathological process of cartilage degeneration is still unknown. Cartilage damage due to osteoarthritis is believed to be mechanically induced. Hence, to investigate the initiation of osteoarthritis the stresses and strains in the cartilage must be determined. So far the most common method to accomplish that is finite element analysis. This paper provides an overview of computational descriptions developed for this purpose, and what they can be used for. Articular cartilage composition and structure are discussed in relation with degenerative changes, and how these affect mechanical properties.

Animals↗

PubMed2005

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.

Animals↗

PubMed2004

Stresses in the local collagen network of articular cartilage: a poroviscoelastic fibril-reinforced finite element study.

Osteoarthritis (OA) is a multifactorial disease, resulting in diarthrodial joint wear and eventually destruction. Swelling of cartilage, which is proportional to the amount of collagen damage, is an initial event of cartilage degeneration, so damage to the collagen fibril network is likely to be one of the earliest signs of OA cartilage degeneration. We propose that the local stresses and strains in the collagen fibrils, which cause the damage, cannot be determined dependably without taking the local arcade-like collagen-fibril structure into account. We investigate this using a poroviscoelastic fibril-reinforced FEA model. The constitutive fibril properties were determined by fitting numerical data to experimental results of unconfined compression and indentation tests on samples of bovine patellar articular cartilage. It was demonstrated that with this model the stresses and strains in the collagen fibrils can be calculated. It was also exhibited that fibrils with different orientations at the same location can be loaded differently, depending on the local architecture of the collagen network. To the best of our knowledge, the present model is the first that can account for these features. We conclude that the local stresses and strains in the articular cartilage are highly influenced by the local morphology of the collagen-fibril network.

Animals↗

PubMed2003

An ionised/non-ionised dual porosity model of intervertebral disc tissue.

The volume of the intrafibrillar water space--i.e. the water contained inside the collagen fibres--is a key parameter that is relevant to concepts of connective tissue structure and function. Confined compression and swelling experiments on annulus fibrosus samples are interpreted in terms of a dual porosity model that distinguishes between a non-ionised intrafibrillar porosity and an ionised extrafibrillar porosity. Both porosities intercommunicate and are saturated with a monovalent ionic solution, i.c. NaCl. The extrafibrillar fixed charge density of the samples is assessed using radiotracer techniques and the collagen content is evaluated by measurement of hydroxyproline concentration. The interpretation of the experimental data yields values for the intrafibrillar water content, the average activity coefficient of the ions, the Donnan osmotic coefficient, the fraction of intrafibrillar water, the stress-free deformation state, and an effective stress-strain relationship as a function of the radial position in the disc. A linear fit between the second Piola-Kirchhoff effective stress and Green-Lagrange strain yielded an effective stiffness: H(e)=1.087 +/- 0.657 MPa. The average fraction of intrafibrillar water was 1.16 g/g collagen. The results were sensitive to changes in the activity and osmotic coefficients and the fraction of intrafibrillar water. The fixed charge density increased with distance from the outer edge of the annulus, whereas the hydroxyproline decreased.

Animals↗

PubMed2003

Pathways of load-induced cartilage damage causing cartilage degeneration in the knee after meniscectomy.

Results of both clinical and animal studies show that meniscectomy often leads to osteoarthritic degenerative changes in articular cartilage. It is generally assumed that this process of cartilage degeneration is due to changes in mechanical loading after meniscectomy. It is, however, not known why and where this cartilage degeneration starts. Load induced cartilage damage is characterized as either type (1)--damage without disruption of the underlying bone or calcified cartilage layer--or type (2), subchondral fracture with or without damage to the overlying cartilage. We asked the question whether cartilage degeneration after meniscectomy is likely to be initiated by type (1) and/or type (2) cartilage damage. To investigate that we applied an axisymmetric biphasic finite element analysis model of the knee joint. In this model the articular cartilage layers of the tibial and the femoral condyles, the meniscus and the bone underlying the articular cartilage of the tibia plateau were included. The model was validated with data from clinical studies, in which the effects of meniscectomy on contact areas and pressures were measured. It was found that both the maximal values and the distributions of the shear stress in the articular cartilage changed after meniscectomy, and that these changes could lead to both type (1) and type (2) cartilage damage. Hence it likely that the cartilage degeneration seen after meniscectomy is initiated by both type (1) and type (2) cartilage damage.

Anisotropy↗

PubMed2001

Spatial interaction between tissue pressure and skeletal muscle perfusion during contraction.

The vascular waterfall theory attributes decreased muscle perfusion during contraction to increased intramuscular pressure (P(IM)) and concomitant increase in venous resistance. Although P(IM) is distributed during contractions, this theory does not account for heterogeneity. This study hypothesises that pressure heterogeneity could affect the interaction between P(IM) rise and perfusion. Regional tissue perfusion during submaximum (100kPa) tetanic contraction is studied, using a finite element model of perfused contracting skeletal muscle. Capillary flow in muscles with one proximal artery and vein (SIM(1)) and with an additional distal artery and vein (SIM(2)) is compared. Blood flow and pressures at rest and P(IM) during contraction ( approximately 25kPa maximally) are similar between simulations, but capillary flow and venous pressure differ. In SIM(2), venous pressure and capillary flow correspond to P(IM) distribution, whereas capillary flow in SIM(1) is less than 10% of flow in SIM(2), in the muscle half without draining vein. This difference is caused by a high central P(IM), followed by central venous pressure rise, in agreement with the waterfall theory. The high central pressure (SIM(1)), obstructs outflow from the distal veins. Distal venous pressure rises until central blood pressure is reached, although local P(IM) is low. Adding a distal vein (SIM(2)) restores the perfusion. It is concluded that regional effects contribute to the interaction between P(IM) and perfusion during contraction. Unlike stated by the vascular waterfall theory, venous pressure may locally exceed P(IM). Although this can be explained by the principles of this theory, the theory does not include this phenomenon as such.

Animals↗

PubMed2001

Requirements for an artificial intervertebral disc.

Intervertebral disc degeneration is an important social and economic problem. Presently available artificial intervertebral discs (AIDs) are insufficient and the main surgical intervention is still spinal fusion. The objective of the present study is to present a list of requirements for the development of an AID which could replace the human lumbar intervertebral disc and restore its function. The list addresses geometry, stiffness, range of motion, strength, facet joint function, center of rotation, fixation, failsafety and implantation technique. Date are obtained from the literature, quantified where possible and checked for consistency. Existing AIDs are evaluated according to the presented list of requirements. Endplate size is a weak point in existing AIDs. These should be large and fit vertebral bodies to prevent migration. Disc height and wedge angle should be restored, unless this would overstretch ligaments. Finally, stiffness and range of motion in all directions should equal those of the healthy disc, except for the axial rotation to relieve the facet joints.

Biomechanical Phenomena↗

PubMed1999

Skeletal muscle transverse strain during isometric contraction at different lengths.

An important assumption in 2D numerical models of skeletal muscle contraction involves deformation in the third dimension of the included muscle section. The present paper studies the often used plane strain description. Therefore, 3D muscle surface deformation is measured from marker displacements during isometric contractions at various muscle lengths. Longitudinal strains at superficial muscle fibers ( - 14 +/- 2.6% at L0, n = 57) and aponeurosis (0.8 +/- 0.9% at L0) decrease with increasing muscle length. The same holds for transverse muscle surface strains in superficial muscle fibers and aponeurosis, which are comparable at intermediate muscle length, but differ at long and short muscle length. Because transverse strains during isometric contraction change with initial muscle length, it is concluded that the effect of muscle length on muscle deformation cannot be studied in plane strain models. These results do not counteract the use of these models to study deformation in contractions with approximately - 9 % longitudinal muscle fiber strain, as transverse strain in superficial muscle fibers and in aponeurosis tissue is minimal in that case. Aponeurosis surface area change decreases with increasing initial muscle length, but muscle fiber surface area change is - 11%, independent of muscle length. Assuming incompressible muscle material, this means that strain perpendicular to the muscle surface equals 11%. Taking the relationship between transverse and longitudinal muscle fiber strain into account, it is hypothesized that superficial muscle fibers flatten during isometric contractions.

Animals↗

PubMed1998

Mechanical blood-tissue interaction in contracting muscles: a model study.

A finite element (FE) model of blood perfused biological tissue has been developed. Blood perfusion is described by fluid flow through a series of 5 intercommunicating vascular compartments that are embedded in the tissue. Each compartment is characterized by a blood flow permeability tensor, blood volume fraction and vessel compliance. Local non-linear relationships between intra-extra vascular pressure difference and blood volume fraction, and between blood volume fraction and the permeability tensor, are included in the FE model. To test the implementation of these non-linear relations, FE results of blood perfusion in a piece of tissue that is subject to increased intramuscular pressure, are compared to results that are calculated with a lumped parameter (LP) model of blood perfusion. FE simulation of blood flow through a contracting rat calf muscle is performed. The FE model used in this simulation contains a transversely isotropic, non-linearly elastic description of deforming muscle tissue, in which local contraction stress is prescribed as a function of time. FE results of muscle tension, total arterial inflow and total venous outflow of the muscle during contraction, correspond to experimental results of an isometrically and tetanically contracting rat calf muscle.

Animals↗

PubMed1997

Finite-element simulation of blood perfusion in muscle tissue during compression and sustained contraction.

Mechanical interaction between tissue stress and blood perfusion in skeletal muscles plays an important role in blood flow impediment during sustained contraction. The exact mechanism of this interaction is not clear, and experimental investigation of this mechanism is difficult. We developed a finite-element model of the mechanical behavior of blood-perfused muscle tissue, which accounts for mechanical blood-tissue interaction in maximally vasodilated vasculature. Verification of the model was performed by comparing finite-element results of blood pressure and flow with experimental measurements in a muscle that is subject to well-controlled mechanical loading conditions. In addition, we performed simulations of blood perfusion during tetanic, isometric contraction and maximal vasodilation in a simplified, two-dimensional finite-element model of a rat calf muscle. A vascular waterfall in the venous compartment was identified as the main cause for blood flow impediment both in the experiment and in the finite-element simulations. The validated finite-element model offers possibilities for detailed analysis of blood perfusion in three-dimensional muscle models under complicated loading conditions.

Animals↗

PubMed1996

Three-dimensional reconstruction of the rat triceps surae muscle and finite element mesh generation of the gastrocnemius medialis muscle.

In order to simulate blood flow in skeletal muscle, our group has developed a finite element description of perfused skeletal muscle. This model requires input parameters concerning the vascular system, muscle contraction and the geometry of a muscle, including its aponeuroses. The objective of the present paper is to create a geometrical reconstruction of the rat gastrocnemius medialis muscle that can be incorporated in the finite element model. Since this muscle is connected to the plantaris and the gastrocnemius lateralis muscle, a detailed computer graphical reconstruction of the triceps surae muscle, based on histological cross-sections, has been accomplished first. Using this reconstruction, relevant sections were selected to create the finite element mesh of the gastrocnemius medialis muscle. Special attention was payed to the location of the aponeuroses. The mesh can be used in finite element simulations of perfused skeletal muscle.

Animals↗