Acute Repair of Meniscus Root Tear Partially Restores Joint Displacements as Measured With Magnetic Resonance Images and Loading in a Cadaveric Porcine Knee

The meniscus serves important load-bearing functions and protects the underlying articular cartilage. Unfortunately, meniscus tears are common and impair the ability of the meniscus to distribute loads, increasing the risk of developing osteoarthritis. Therefore, surgical repair of the meniscus is a frequently performed procedure; however, repair does not always prevent osteoarthritis. This is hypothesized to be due to altered joint loading post-injury and repair, where the functional deficit of the meniscus prevents it from performing its role of distributing forces. The objective of this study was to quantify joint kinematics in an intact joint, after a meniscus root tear, and after suture repair in cadaveric porcine knees, a frequently used in vivo model. We utilized an magnetic resonance images-compatible loading device and novel use of a T1 vibe sequence to measure meniscus and femur displacements under physiological axial loads. We found that anterior root tear led to large meniscus displacements under physiological axial loading and that suture anchor repair reduced these displacements but did not fully restore intact joint kinematics. After tear and repair, the anterior region of the meniscus moved posteriorly and medially as it was forced out of the joint space under loading, while the posterior region had small displacements as the posterior attachment acted as a hinge about which the meniscus pivoted in the axial plane. Methods from this study can be applied to assess altered joint kinematics following human knee injuries and evaluate repair strategies aimed to restore joint kinematics.

Overload in a Rat In Vivo Model of Synergist Ablation Induces Tendon Multiscale Structural and Functional Degeneration

Tendon degeneration is typically described as an overuse injury with little distinction made between magnitude of load (overload) and number of cycles (overuse). Further, in vivo, animal models of tendon degeneration are mostly overuse models, where tendon damage is caused by a high number of load cycles. As a result, there is a lack of knowledge of how isolated overload leads to degeneration in tendons. A surgical model of synergist ablation (SynAb) overloads the target tendon, plantaris, by ablating its synergist tendon, Achilles. The objective of this study was to evaluate the structural and functional changes that occur following overload of plantaris tendon in a rat SynAb model. Tendon cross-sectional area (CSA) and shape changes were evaluated by longitudinal MR imaging up to 8 weeks postsurgery. Tissue-scale structural changes were evaluated by semiquantified histology and second harmonic generation microscopy. Fibril level changes were evaluated with serial block face scanning electron microscopy (SBF-SEM). Functional changes were evaluated using tension tests at the tissue and microscale using a custom testing system allowing both video and microscopy imaging. At 8 weeks, overloaded plantaris tendons exhibited degenerative changes including increases in CSA, cell density, collagen damage area fraction (DAF), and fibril diameter, and decreases in collagen alignment, modulus, and yield stress. To interpret the differences between overload and overuse in tendon, we introduce a new framework for tendon remodeling and degeneration that differentiates between the inputs of overload and overuse. In summary, isolated overload induces multiscale degenerative structural and functional changes in plantaris tendon.

Volume Loss and Recovery in Bovine Knee Meniscus Loaded in Circumferential Tension

Load-induced volume change is an important aspect of knee meniscus function because volume loss creates fluid pressure, which minimizes friction and helps support compressive loads. The knee meniscus is unusual amongst cartilaginous tissues in that it is loaded not only in axial compression, but also in circumferential tension between its tibial attachments. Despite the physiologic importance of the knee meniscus' tensile properties, its volumetric strain in tension has never been directly measured, and predictions of volume strain in the scientific literature are inconsistent. In this study, we apply uniaxial tension to bovine knee meniscus and use biplanar imaging to directly observe the resulting three-dimensional volume change and unloaded recovery, revealing that tension causes volumetric contraction. Compression is already known to also cause contraction; therefore, all major physiologic loads compress and pressurize the meniscus, inducing fluid outflow. Although passive unloaded recovery is often described as slow relative to loaded loss, here we show that at physiologic strains the volume recovery rate in the meniscus upon unloading is faster than the rate of volume loss. These measurements of volumetric strain are an important step toward a complete theory of knee meniscus fluid flow and load support.

Acute repair of meniscus root tear partially restores joint displacements as measured with MRI and loading in a porcine knee

The meniscus serves important load-bearing functions and protects the underlying articular cartilage. Unfortunately, meniscus tears are common and impair the ability of the meniscus to distribute loads, greatly increasing the risk for developing osteoarthritis. Therefore, surgical repair of the meniscus is a frequently performed procedure; however, this repair does not always prevent osteoarthritis. This is hypothesized to be due to altered joint loading post injury and repair, where the functional deficit of the meniscus prevents it from performing its role of distributing forces. However, many studies of meniscus function required opening the joint, which alters kinematics. The objective of this study was to use novel MRI methods to image the intact joint under axial load and measure the acute meniscus and femur displacements in an intact joint, after a meniscus root tear, and after suture repair in the porcine knee, a frequently used in vivo model. We found that anterior root tear led to large meniscus and femur displacements under physiological axial loading, and that suture anchor repair reduced these displacements, but did not fully restore intact joint kinematics. After tear and repair, the anterior region of the meniscus moved posteriorly and medially as it was forced out of the joint space under loading, while the posterior region had small displacements as the posterior attachment acted as a hinge about which the meniscus rotated in the axial plane. This technique can be applied to evaluate the effect of knee injuries and to develop improved repair strategies to restore joint kinematics.

MRI ‐based measurement of in vivo disc mechanics in a young population due to flexion, extension, and diurnal loading

Background

Intervertebral disc degeneration is often implicated in low back pain; however, discs with structural degeneration often do not cause pain. It may be that disc mechanics can provide better diagnosis and identification of the pain source. In cadaveric testing, the degenerated disc has altered mechanics, but in vivo, disc mechanics remain unknown. To measure in vivo disc mechanics, noninvasive methods must be developed to apply and measure physiological deformations.

Aim

Thus, this study aimed to develop methods to measure disc mechanical function via noninvasive MRI during flexion and extension and after diurnal loading in a young population. This data will serve as baseline disc mechanics to later compare across ages and in patients.

Materials & Methods

To accomplish this, subjects were imaged in the morning in a reference supine position, in flexion, in extension, and at the end of the day in a supine position. Disc deformations and vertebral motions were used to quantify disc axial strain, changes in wedge angle, and anterior-posterior (A-P) shear displacement. T₂ weighted MRI was also used to evaluate disc degeneration via Pfirrmann grading and T₂ time. All measures were then tested for effect of sex and disc level.

Results

We found that flexion and extension caused level-dependent strains in the anterior and posterior of the disc, changes in wedge angle, and A-P shear displacements. Flexion had higher magnitude changes overall. Diurnal loading did not cause level-dependent strains but did cause small level-dependent changes in wedge angle and A-P shear displacements.

Discussion

Correlations between disc degeneration and mechanics were largest in flexion, likely due to the smaller contribution of the facet joints in this condition.

Conclusion

In summary, this study established methods to measure in vivo disc mechanical function via noninvasive MRI and established a baseline in a young population that may be compared to older subjects and clinical disorders in the future.

Multiaxial validation of a finite element model of the intervertebral disc with multigenerational fibers to establish residual strain

Finite element models of the intervertebral disc are used to address research questions that cannot be tested through typical experimentation. A disc model requires complex geometry and tissue properties to be accurately defined to mimic the physiological disc. The physiological disc possesses residual strain in the annulus fibrosus (AF) due to osmotic swelling and due to inherently pre-strained fibers. We developed a disc model with residual contributions due to swelling-only, and a multigeneration model with residual contributions due to both swelling and AF fiber pre-strain and validated it against organ-scale uniaxial, quasi-static and multiaxial, dynamic mechanical tests. In addition, we demonstrated the models' ability to mimic the opening angle observed following radial incision of bovine discs. Both models were validated against organ-scale experimental data. While the swelling only model responses were within the experimental 95% confidence interval, the multigeneration model offered outcomes closer to the experimental mean and had a bovine model opening angle within one SD of the experimental mean. The better outcomes for the multigeneration model, which allowed for the inclusion of inherently pre-strained fibers in AF, is likely due to its uniform fiber contribution throughout the AF. We conclude that the residual contribution of pre-strained fibers in the AF should be included to best simulate the physiological disc and its behaviors.

Structure, function, and defect tolerance with maturation of the radial tie fiber network in the knee meniscus

The knee menisci are comprised of two orthogonal collagenous networks-circumferential and radial-that combine to enable efficient load bearing by the tissue in adults. Here, we assessed how the structural and functional characteristics of these networks developed over the course of skeletal maturation and determined the role of these fiber networks in defect tolerance with tissue injury. Imaging of the radial tie fiber (RTF) collagen structure in medial bovine menisci from fetal, juvenile, and adult specimens showed increasing heterogeneity, anisotropy, thickness, and density with skeletal development. The mechanical analysis showed that the tensile modulus in the radial direction did not change with skeletal development, though the resilience (in the radial direction) increased and the tolerance to defects in the circumferential direction decreased, in adult compared to fetal tissues. This loss of defect tolerance correlated with increased order in the RTF network in adult tissue. These data provide new insights into the role of the radial fiber network in meniscus function, will lead to improved clinical decision-making in the presence of a tear and may improve engineering efforts to reproduce this critical load-bearing structure in the knee.

Impact of pulse sequence, analysis method, and signal to noise ratio on the accuracy of intervertebral disc T 2 measurement

Noninvasive assessments of intervertebral disc health and degeneration are critical for addressing disc degeneration and low back pain. Magnetic resonance imaging (MRI) is exceptionally sensitive to tissue with high water content, and measurement of the MR transverse relaxation time, T 2, has been applied as a quantitative, continuous, and objective measure of disc degeneration that is linked to the water and matrix composition of the disc. However, T 2 measurement is susceptible to inaccuracies due to Rician noise, T 1 contamination, and stimulated echo effects. These error generators can all be controlled for with proper data collection and fitting methods. The objective of this study was to identify sequence parameters to appropriately acquire MR data and to establish curve fitting methods to accurately calculate disc T 2 in the presence of noise by correcting for Rician noise. To do so, we compared T 2 calculated from the typical monoexponential (MONO) fits and noise corrected exponential (NCEXP) fits. We examined how the selected sequence parameters altered the calculated T 2 in silico and in vivo. Typical MONO fits were frequently poor due to Rician noise, and NCEXP fits were more likely to provide accurate T 2 calculations. NCEXP is particularly less biased and less uncertain at low SNR. This study showed that the NCEXP using sequences with data from 20 echoes out to echo times of ~300 ms is the best method for calculating T 2 of discs. By acquiring signal data out to longer echo times and accounting for Rician noise, the curve fitting is more robust in calculating T 2 despite the noise in the data. This is particularly important when considering degenerate discs or AF tissue because the SNR of these regions is lower.

Helical fibrillar microstructure of tendon using serial block-face scanning electron microscopy and a mechanical model for interfibrillar load transfer

Tendon's hierarchical structure allows for load transfer between its fibrillar elements at multiple length scales. Tendon microstructure is particularly important, because it includes the cells and their surrounding collagen fibrils, where mechanical interactions can have potentially important physiological and pathological contributions. However, the three-dimensional (3D) microstructure and the mechanisms of load transfer in that length scale are not known. It has been postulated that interfibrillar matrix shear or direct load transfer via the fusion/branching of small fibrils are responsible for load transfer, but the significance of these mechanisms is still unclear. Alternatively, the helical fibrils that occur at the microstructural scale in tendon may also mediate load transfer; however, these structures are not well studied due to the lack of a three-dimensional visualization of tendon microstructure. In this study, we used serial block-face scanning electron microscopy to investigate the 3D microstructure of fibrils in rat tail tendon. We found that tendon fibrils have a complex architecture with many helically wrapped fibrils. We studied the mechanical implications of these helical structures using finite-element modelling and found that frictional contact between helical fibrils can induce load transfer even in the absence of matrix bonding or fibril fusion/branching. This study is significant in that it provides a three-dimensional view of the tendon microstructure and suggests friction between helically wrapped fibrils as a mechanism for load transfer, which is an important aspect of tendon biomechanics.

Off-axis response due to mechanical coupling across all six degrees of freedom in the human disc

The kinematics of the intervertebral disc are defined by six degrees of freedom (DOF): three translations (Tz: axial compression, Tx: lateral shear, and Ty: anterior-posterior shear) and three rotations (Rz: torsion, Rx: flexion-extension, and Ry: lateral bending). There is some evidence that the six DOFs are mechanically coupled, such that loading in one DOF affects the mechanics of the other five "off-axis" DOFs, however, most studies have not controlled and/or measured all six DOFs simultaneously. Additionally, the relationships between disc geometry and disc mechanics are important for evaluation of data from different sized donor and patient discs. The objectives of this study were to quantify the mechanical behavior of the intervertebral disc in all six degrees of freedom (DOFs), measure the coupling between the applied motion in each DOF with the resulting off-axis motions, and test the hypothesis that disc geometry influences these mechanical behaviors. All off-axis displacements and rotations were significantly correlated with the applied DOF and were of similar magnitude as physiologically relevant motion, confirming that off-axis coupling is an important mechanical response. Interestingly, there were pairs of DOFs that were especially strongly coupled: lateral shear (Tx) and lateral bending (Ry), anterior-posterior shear (Ty) and flexion-extension (Rx), and compression (Tz) and torsion (Rz). Large off-axis shears may contribute to injury risk in bending and flexion. In addition, the disc responded to shear (Tx, Ty) and rotational loading (Rx, Ry, and Rz) by increasing in disc height in order to maintain the applied compressive load. Quantifying these mechanical behaviors across all six DOF are critical for designing and testing disc therapies, such as implants and tissue engineered constructs, and also for validating finite element models.

Short cracks in knee meniscus tissue cause strain concentrations, but do not reduce ultimate stress, in single-cycle uniaxial tension

Tears are central to knee meniscus pathology and, from a mechanical perspective, are crack-like defects (cracks). In many materials, cracks create stress concentrations that cause progressive local rupture and reduce effective strength. It is currently unknown if cracks in meniscus have these consequences; if they do, this would have repercussions for management of meniscus pathology. The objective of this study was to determine if a short crack in meniscus tissue, which mimics a preclinical meniscus tear, (a) causes crack growth and reduces effective strength, (b) creates a near-tip strain concentration and (c) creates unloaded regions on either side of the crack. Specimens with and without cracks were tested in uniaxial tension and compared in terms of macroscopic stress-strain curves and digital image correlation strain fields. The strain fields were used as an indicator of stress concentrations and unloaded regions. Effective strength was found to be insensitive to the presence of a crack (potential effect < 0.86 s.d.; β = 0.2), but significant strain concentrations, which have the potential to lead to long-term accumulation of tissue or cell damage, were observed near the crack tip.

Strain Distribution of Intact Rat Rotator Cuff Tendon-to-Bone Attachments and Attachments With Defects

This study aimed to experimentally track the tissue-scale strains of the tendon-bone attachment with and without a localized defect. We hypothesized that attachments with a localized defect would develop strain concentrations and would be weaker than intact attachments. Uniaxial tensile tests and digital image correlation were performed on rat infraspinatus tendon-to-bone attachments with defects (defect group) and without defects (intact group). Biomechanical properties were calculated, and tissue-scale strain distributions were quantified for superior and inferior fibrous and calcified regions. At the macroscale, the defect group exhibited reduced stiffness (31.3±3.7 N/mm), reduced ultimate load (24.7±3.8 N), and reduced area under the curve at ultimate stress (3.7±1.5 J/m2) compared to intact attachments (42.4±4.3 N/mm, 39.3±3.7 N, and 5.6±1.4 J/m2, respectively). Transverse strain increased with increasing axial load in the fibrous region of the defect group but did not change for the intact group. Shear strain of the superior fibrous region was significantly higher in the defect group compared to intact group near yield load. This work experimentally identified that attachments may resist failure by distributing strain across the interface and that strain concentrations develop near attachment defects. By establishing the tissue-scale deformation patterns of the attachment, we gained insight into the micromechanical behavior of this interfacial tissue and bolstered our understanding of the deformation mechanisms associated with its ability to resist failure.

Advances in Quantification of Meniscus Tensile Mechanics Including Nonlinearity, Yield, and Failure

The meniscus provides crucial knee function and damage to it leads to osteoarthritis of the articular cartilage. Accurate measurement of its mechanical properties is therefore important, but there is uncertainty about how the test procedure affects the results, and some key mechanical properties are reported using ad hoc criteria (modulus) or not reported at all (yield). This study quantifies the meniscus' stress-strain curve in circumferential and radial uniaxial tension. A fiber recruitment model was used to represent the toe region of the stress-strain curve, and new reproducible and objective procedures were implemented for identifying the yield point and measuring the elastic modulus. Patterns of strain heterogeneity were identified using strain field measurements. To resolve uncertainty regarding whether rupture location (i.e., midsubstance rupture versus at-grip rupture) influences the measured mechanical properties, types of rupture were classified in detail and compared. Dogbone (DB)-shaped specimens are often used to promote midsubstance rupture; to determine if this is effective, we compared DB and rectangle (R) specimens in both the radial and circumferential directions. In circumferential testing, we also compared expanded tab (ET) specimens under the hypothesis that this shape would more effectively secure the meniscus' curved fibers and thus produce a stiffer response. The fiber recruitment model produced excellent fits to the data. Full fiber recruitment occurred approximately at the yield point, strongly supporting the model's physical interpretation. The strain fields, especially shear and transverse strain, were extremely heterogeneous. The shear strain field was arranged in pronounced bands of alternating positive and negative strain in a pattern similar to the fascicle structure. The site and extent of failure showed great variation, but did not affect the measured mechanical properties. In circumferential tension, ET specimens underwent earlier and more rapid fiber recruitment, had less stretch at yield, and had greater elastic modulus and peak stress. No significant differences were observed between R and DB specimens in either circumferential or radial tension. Based on these results, ET specimens are recommended for circumferential tests and R specimens for radial tests. In addition to the data obtained, the procedural and modeling advances made in this study are a significant step forward for meniscus research and are applicable to other fibrous soft tissues.

MRI quantification of human spine cartilage endplate geometry: Comparison with age, degeneration, level, and disc geometry

Geometry is an important indicator of disc mechanical function and degeneration. While the geometry and associated degenerative changes in the nucleus pulposus and the annulus fibrosus are well-defined, the geometry of the cartilage endplate (CEP) and its relationship to disc degeneration are unknown. The objectives of this study were to quantify CEP geometry in three dimensions using an MRI FLASH imaging sequence and evaluate relationships between CEP geometry and age, degeneration, spinal level, and overall disc geometry. To do so, we assessed the MRI-based measurements for accuracy and repeatability. Next, we measured CEP geometry across a larger sample set and correlated CEP geometric parameters to age, disc degeneration, level, and disc geometry. The MRI-based measures resulted in thicknesses (0.3-1 mm) that are comparable to prior measurements of CEP thickness. CEP thickness was greatest at the anterior/posterior (A/P) margins and smallest in the center. The CEP A/P thickness, axial area, and lateral width decreased with age but were not related to disc degeneration. Age-related, but not degeneration-related, changes in geometry suggest that the CEP may not follow the progression of disc degeneration. Ultimately, if the CEP undergoes significant geometric changes with aging and if these can be related to low back pain, a clinically feasible translation of the FLASH MRI-based measurement of CEP geometry presented in this study may prove a useful diagnostic tool. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1410-1417, 2016.

Novel human intervertebral disc strain template to quantify regional three-dimensional strains in a population and compare to internal strains predicted by a finite element model

Tissue strain is an important indicator of mechanical function, but is difficult to noninvasively measure in the intervertebral disc. The objective of this study was to generate a disc strain template, a 3D average of disc strain, of a group of human L4-L5 discs loaded in axial compression. To do so, magnetic resonance images of uncompressed discs were used to create an average disc shape. Next, the strain tensors were calculated pixel-wise by using a previously developed registration algorithm. Individual disc strain tensor components were then transformed to the template space and averaged to create the disc strain template. The strain template reduced individual variability while highlighting group trends. For example, higher axial and circumferential strains were present in the lateral and posterolateral regions of the disc, which may lead to annular tears. This quantification of group-level trends in local 3D strain is a significant step forward in the study of disc biomechanics. These trends were compared to a finite element model that had been previously validated against the disc-level mechanical response. Depending on the strain component, 81-99% of the regions within the finite element model had calculated strains within one standard deviation of the template strain results. The template creation technique provides a new measurement technique useful for a wide range of studies, including more complex loading conditions, the effect of disc pathologies and degeneration, damage mechanisms, and design and evaluation of treatments. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1264-1273, 2016.

A comparison of stress in cracked fibrous tissue specimens with varied crack location, loading, and orientation using finite element analysis

Cracks in fibrous soft tissue, such as intervertebral disc annulus fibrosus and knee meniscus, cause pain and compromise joint mechanics. A crack concentrates stress at its tip, making further failure and crack extension (fracture) more likely. Ex vivo mechanical testing is an important tool for studying the loading conditions required for crack extension, but prior work has shown that it is difficult to reproduce crack extension. Most prior work used edge crack specimens in uniaxial tension, with the crack 90° to the edge of the specimen. This configuration does not necessarily represent the loading conditions that cause in vivo crack extension. To find a potentially better choice for experiments aiming to reproduce crack extension, we used finite element analysis to compare, in factorial combination, (1) center crack vs. edge crack location, (2) biaxial vs. uniaxial loading, and (3) crack-fiber angles ranging from 0° to 90°. The simulated material was annulus fibrosus fibrocartilage with a single fiber family. We hypothesized that one of the simulated test cases would produce a stronger stress concentration than the commonly used uniaxially loaded 90° crack-fiber angle edge crack case. Stress concentrations were compared between cases in terms of fiber-parallel stress (representing risk of fiber rupture), fiber-perpendicular stress (representing risk of matrix rupture), and fiber shear stress (representing risk of fiber sliding). Fiber-perpendicular stress and fiber shear stress concentrations were greatest in edge crack specimens (of any crack-fiber angle) and center crack specimens with a 90° crack-fiber angle. However, unless the crack is parallel to the fiber direction, these stress components alone are insufficient to cause crack opening and extension. Fiber-parallel stress concentrations were greatest in center crack specimens with a 45° crack-fiber angle, either biaxially or uniaxially loaded. We therefore recommend that the 45° center crack case be tried in future experiments intended to study crack extension by fiber rupture.

Internal three-dimensional strains in human intervertebral discs under axial compression quantified noninvasively by magnetic resonance imaging and image registration

Study objectives were to develop, validate, and apply a method to measure three-dimensional (3D) internal strains in intact human discs under axial compression. A custom-built loading device applied compression and permitted load-relaxation outside of the magnet while also maintaining compression and hydration during imaging. Strain was measured through registration of 300 μm isotropic resolution images. Excellent registration accuracy was achieved, with 94% and 65% overlap of disc volume and lamellae compared to manual segmentation, and an average Hausdorff, a measure of distance error, of 0.03 and 0.12 mm for disc volume and lamellae boundaries, respectively. Strain maps enabled qualitative visualization and quantitative regional annulus fibrosus (AF) strain analysis. Axial and circumferential strains were highest in the lateral AF and lowest in the anterior and posterior AF. Radial strains were lowest in the lateral AF, but highly variable. Overall, this study provided new methods that will be valuable in the design and evaluation surgical procedures and therapeutic interventions.

Human L3L4 intervertebral disc mean 3D shape, modes of variation, and their relationship to degeneration

Intervertebral disc mechanics are affected by both disc shape and disc degeneration, which in turn each affect the other; disc mechanics additionally have a role in the etiology of disc degeneration. Finite element analysis (FEA) is a favored tool to investigate these relationships, but limited data for intervertebral disc 3D shape has forced the use of simplified or single-subject geometries, with the effect of inter-individual shape variation investigated only in specialized studies. Similarly, most data on disc shape variation with degeneration is based on 2D mid-sagittal images, which incompletely define 3D shape changes. Therefore, the objective of this study was to quantify inter-individual disc shape variation in 3D, classify this variation into independently-occurring modes using a statistical shape model, and identify correlations between disc shape and degeneration. Three-dimensional disc shapes were obtained from MRI of 13 human male cadaver L3L4 discs. An average disc shape and four major modes of shape variation (representing 90% of the variance) were identified. The first mode represented disc axial area and was significantly correlated to degeneration (R(2)=0.44), indicating larger axial area in degenerate discs. Disc height variation occurred in three distinct modes, each also involving non-height variation. The statistical shape model provides an average L3L4 disc shape for FEA that is fully defined in 3D, and makes it convenient to generate a set of shapes with which to represent aggregate inter-individual variation. Degeneration grade-specific shapes can also be generated. To facilitate application, the model is included in this paper׳s supplemental content.

Evaluation of strengthening mechanisms in calcite single crystals from mollusk shells

Biogenic single-crystal calcite is often reported to be harder and tougher than geologic calcite in the form of Iceland spar. However, the mechanistic origins of the superior mechanical properties of the biogenic materials are still debated. We investigate the hardness and modulus of biogenic calcite from the prismatic layer of the mollusk Atrina rigida compared with a pure geologic calcite, Iceland spar. On the {001} face, biogenic calcite is found to be 50-70% harder than geologic calcite. This range is due to the fact that changes in azimuthal angle of the indenter tip lead to a hardness variation of ∼20% in A. rigida but only ∼7% in Iceland spar. The higher hardness and increased anisotropy of biogenic calcite could be accounted for by hardening mechanisms based on hindered dislocation motion rather than crack deflection.

Biaxial tensile testing and constitutive modeling of human supraspinatus tendon

The heterogeneous composition and mechanical properties of the supraspinatus tendon offer an opportunity for studying the structure-function relationships of fibrous musculoskeletal connective tissues. Previous uniaxial testing has demonstrated a correlation between the collagen fiber angle distribution and tendon mechanics in response to tensile loading both parallel and transverse to the tendon longitudinal axis. However, the planar mechanics of the supraspinatus tendon may be more appropriately characterized through biaxial tensile testing, which avoids the limitation of nonphysiologic traction-free boundary conditions present during uniaxial testing. Combined with a structural constitutive model, biaxial testing can help identify the specific structural mechanisms underlying the tendon's two-dimensional mechanical behavior. Therefore, the objective of this study was to evaluate the contribution of collagen fiber organization to the planar tensile mechanics of the human supraspinatus tendon by fitting biaxial tensile data with a structural constitutive model that incorporates a sample-specific angular distribution of nonlinear fibers. Regional samples were tested under several biaxial boundary conditions while simultaneously measuring the collagen fiber orientations via polarized light imaging. The histograms of fiber angles were fit with a von Mises probability distribution and input into a hyperelastic constitutive model incorporating the contributions of the uncrimped fibers. Samples with a wide fiber angle distribution produced greater transverse stresses than more highly aligned samples. The structural model fit the longitudinal stresses well (median R(2) ≥ 0.96) and was validated by successfully predicting the stress response to a mechanical protocol not used for parameter estimation. The transverse stresses were fit less well with greater errors observed for less aligned samples. Sensitivity analyses and relatively affine fiber kinematics suggest that these errors are not due to inaccuracies in measuring the collagen fiber organization. More likely, additional strain energy terms representing fiber-fiber interactions are necessary to provide a closer approximation of the transverse stresses. Nevertheless, this approach demonstrated that the longitudinal tensile mechanics of the supraspinatus tendon are primarily dependent on the moduli, crimp, and angular distribution of its collagen fibers. These results add to the existing knowledge of structure-function relationships in fibrous musculoskeletal tissue, which is valuable for understanding the etiology of degenerative disease, developing effective tissue engineering design strategies, and predicting outcomes of tissue repair.

Age-related intimal stiffening enhances endothelial permeability and leukocyte transmigration

Age is the most significant risk factor for atherosclerosis; however, the link between age and atherosclerosis is poorly understood. During both aging and atherosclerosis progression, the blood vessel wall stiffens owing to alterations in the extracellular matrix. Using in vitro and ex vivo models of vessel wall stiffness and aging, we show that stiffening of extracellular matrix within the intima promotes endothelial cell permeability--a hallmark of atherogenesis. When cultured on hydrogels fabricated to match the elasticity of young and aging intima, endothelial monolayers exhibit increased permeability and disrupted cell-cell junctions on stiffer matrices. In parallel experiments, we showed a corresponding increase in cell-cell junction width with age in ex vivo aortas from young (10 weeks) and old (21 to 25 months) healthy mice. To investigate the mechanism by which matrix stiffening alters monolayer integrity, we found that cell contractility increases with increased matrix stiffness, mechanically destabilizing cell-cell junctions. This increase in endothelial permeability results in increased leukocyte extravasation, which is a critical step in atherosclerotic plaque formation. Mild inhibition of Rho-dependent cell contractility using Y-27632, an inhibitor of Rho-associated kinase, or small interfering RNA restored monolayer integrity in vitro and in vivo. Our results suggest that extracellular matrix stiffening alone, which occurs during aging, can lead to endothelial monolayer disruption and atherosclerosis pathogenesis. Because previous therapeutics designed to decrease vascular stiffness have been met with limited success, our findings could be the basis for the design of therapeutics that target the Rho-dependent cellular contractile response to matrix stiffening, rather than stiffness itself, to more effectively prevent atherosclerosis progression.

Pulsed and parallel-polarization EPR characterization of the photosystem II oxygen-evolving complex

Photosystem II uses visible light to drive the oxidation of water, resulting in bioactivated electrons and protons, with the production of molecular oxygen as a byproduct. This water-splitting reaction is carried out by a manganese cluster/tyrosine radial ensemble, the oxygen -evolving complex. Although conventional continuous-wave, perpendicular -polarization electron paramagnetic resonance (EPR) spectroscopy has significantly advanced our knowledge of the structure and function of the oxygen-evolving complex, significant additional information can be obtained with the application of additional EPR methodologies. Specifically, parallel-polarization EPR spectroscopy can be use to obtain highly resolved EPR spectra of integer spin Mn species, and pulsed EPR spectroscopy with electron spin echo-based sequences, such as electron spin echo envelope modulation and electron spin echo-electron nuclear double resonance, can be used to measure weak interactions obscured in continuous-wave spectroscopy by inhomogeneous broadening.

EPR/ENDOR characterization of the physical and electronic structure of the OEC Mn cluster

Electron paramagnetic resonance (EPR) spectroscopy has often played a crucial role in characterizing the various cofactors and processes of photosynthesis, and photosystem II and its oxygen evolving chemistry is no exception. Until recently, the application of EPR spectroscopy to the characterization of the oxygen evolving complex (OEC) has been limited to the S2-state of the Kok cycle. However, in the past few years, continuous wave-EPR signals have been obtained for both the S0- and S1-state as well as for the S2 (radical)(Z)-state of a number of inhibited systems. Furthermore, the pulsed EPR technique of electron spin echo electron nuclear double resonance spectroscopy has been used to directly probe the 55Mn nuclei of the manganese cluster. In this review, we discuss how the EPR data obtained from each of these states of the OEC Kok cycle are being used to provide insight into the physical and electronic structure of the manganese cluster and its interaction with the key tyrosine, Y(Z).

Computational and experimental analyses reveal previously undetected coding exons of the KRIT1 (CCM1) gene

A notable difficulty in annotating genomic sequence is identifying the correct start codon in a gene. An important such case has been found with KRIT1, the cerebral cavernous malformation type 1 (CCM1) gene. Analysis of human and mouse genomic sequence encompassing the region containing KRIT1/Krit1 using exon/gene-prediction and comparative alignment programs revealed putative exons upstream of the previously described first exon. These additional candidate exons show significant matches to mouse and human ESTs that are contiguous with and extend upstream from the previously designated 5' end of the KRIT1 cDNA sequence. RT-PCR and 5'RACE experiments confirm the presence of four additional upstream coding exons that encode an additional 207 amino acids. Importantly, a novel frameshift mutation in one of these newly identified KRIT1 exons has been found in a CCM1 family. These data establish the authentic KRIT1 amino acid sequence and suggest that the additional KRIT1 exons may harbor mutations in other CCM1 families. In addition, these results provide another example of the utility of rigorous computational and comparative sequence analysis for refining gene structure.

ESEEM studies of succinate:ubiquinone reductase from Paracoccus denitrificans

Electron spin-echo envelope modulation (ESEEM) spectroscopy has been performed in order to obtain structural information about the environment of the reduced [2Fe-2S] cluster (S-1 center), the oxidized [3Fe-4S] cluster (S-3 center), and the flavin semiquinone radical in purified succinate:ubiquinone reductase from Paracoccus denitrificans. Spectral simulations of the ESEEM data from the reduced [2Fe-2S] yielded nuclear quadrupole interaction parameters that are indicative of peptide nitrogens. We also observed a weak interaction between the oxidized [3Fe-4S] cluster and a peptide 14N. There was no evidence for coordination of any of the Fe atoms to 14N atoms of imidazole rings. The ESEEM data from the flavin semiquinone radical were more complicated. Here, evidence was obtained for interactions between the unpaired electron and only the two nitrogen atoms in the flavin ring.

Histidine 332 of the D1 polypeptide modulates the magnetic and redox properties of the manganese cluster and tyrosine Y(Z) in photosystem II

An electron spin-echo envelope modulation study [Tang, X.-S., Diner, B. A., Larsen, B. S., Gilchrist, M. L., Jr., Lorigan, G. A., and Britt, R. D. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 704-708] and a recent Fourier transform infrared study [Noguchi, T., Inoue, Y., and Tang, X.-S. (1999) Biochemistry 38, 10187-10195], both conducted with [(15)N]histidine-labeled photosystem II particles, show that at least one histidine residue coordinates the O(2)-evolving Mn cluster in photosystem II. Evidence obtained from site-directed mutagenesis studies suggests that one of these residues may be His332 of the D1 polypeptide. The mutation D1-H332E is of particular interest because cells of the cyanobacterium Synechocystis sp. PCC 6803 that contain this mutation evolve no O(2) but appear to assemble Mn clusters in nearly all photosystem II reaction centers [Chu, H.-A., Nguyen, A. P. , and Debus, R. J. (1995) Biochemistry 34, 5859-5882]. Photosystem II particles isolated from the Synechocystis D1-H332E mutant are characterized in this study. Intact D1-H332E photosystem II particles exhibit an altered S(2) state multiline EPR signal that has more hyperfine lines and narrower splittings than the S(2) state multiline EPR signal observed in wild-type PSII particles. However, the quantum yield for oxidizing the S(1) state Mn cluster is very low, corresponding to an 8000-fold slowing of the rate of Mn oxidation by Y(Z)(*), and the temperature threshold for forming the S(2) state is approximately 100 K higher than in wild-type PSII preparations. Furthermore, the D1-H332E PSII particles are unable to advance beyond the Y(Z)(*)S(2) state, as shown by the accumulation of a narrow "split" EPR signal under multiple turnover conditions. In Mn-depleted photosystem II particles, charge recombination between Q(A)(*)(-) and Y(Z)(*) in D1-H332E is accelerated in comparison to wild-type, showing that the mutation alters the redox properties of Y(Z) in addition to those of the Mn cluster. These results are consistent with D1-His332 being located near the Mn-Y(Z) complex and perhaps ligating Mn.

An electron spin-echo envelope modulation (ESEEM) study of the QA binding pocket of PS II reaction centers from spinach and Synechocystis

We have used electron spin-echo envelope modulation spectroscopy (ESEEM) to characterize the protein-cofactor interactions present in the QA- binding pocket of PS II centers isolated from spinach and Synechocystis. We conclude that the ESEEM spectrum of QA- is the result of interactions of the S = 1/2 electron spin of QA- with the I = 1 nuclear spins of the peptide nitrogens of two different amino acids. One peptide nitrogen has ESEEM peaks near 0.7, 2.0, 2.85, and 5.0 MHz with isotropic and dipolar hyperfine couplings of Aiso = 2.0 MHz and Adip = 0.25 MHz, respectively. On the basis of these hyperfine couplings we predict the existence of a strong hydrogen bond between QA- and the peptide nitrogen with a hydrogen bond distance of about 2 A. We have not identified the amino acid origin of this peptide nitrogen. By using amino acid specific isotopic labeling in conjunction with site-directed mutagenesis, we demonstrate that the second peptide nitrogen is that of D2-Ala260, with ESEEM peaks near 0.6 and 1.5 MHz and an isotropic hyperfine coupling, Aiso, less than 0.2 MHz. This small isotropic coupling suggests that the D2-Ala260 peptide nitrogen at best forms a weak hydrogen bond with QA-.

The 23 and 17 kDa extrinsic proteins of photosystem II modulate the magnetic properties of the S1-state manganese cluster

An S1-state parallel polarization "multiline" EPR signal arising from the oxygen-evolving complex has been detected in spinach (PSII) membrane and core preparations depleted of the 23 and 17 kDa extrinsic polypeptides, but retaining the 33 kDa extrinsic protein. This S1-state multiline signal, with an effective g value of 12 and at least 18 hyperfine lines, has previously been detected only in PSII preparations from the cyanobacterium sp. Synechocystis sp. PCC6803 [Campbell, K. A., Peloquin, J. M., Pham, D. P., Debus, R. J., and Britt, R. D. (1998) J. Am. Chem. Soc. 120, 447-448]. It is absent in PSII spinach membrane and core preparations that either fully retain or completely lack the 33, 23, and 17 kDa extrinsic proteins. The S1-state multiline signal detected in spinach PSII cores and membranes has the same effective g value and hyperfine spacing as the signal detected in Synechocystis PSII particles. This signal provides direct evidence for the influence of the extrinsic PSII proteins on the magnetic properties of the Mn cluster.

Time-dependent thermodynamics during early electron transfer in reaction centers from Rhodobacter sphaeroides

The temperature dependence of fluorescence on the picosecond to nanosecond time scale from the reaction centers of Rhodobacter sphaeroides strain R-26 and two mutants with elevated P/P+ midpoint potentials has been measured with picosecond time resolution. In all three samples, the kinetics of the fluorescence decay is complex and can only be well described with four or more exponential decay terms spanning the picosecond to nanosecond time range. Multiexponential fits are needed at all temperatures between 295 and 20 K. The complex decay kinetics are explained in terms of a dynamic solvation model in which the charge-separated state is stabilized after formation by protein conformational changes. Many of these motions have not had time to occur on the time scale of initial electron transfer and/or are frozen out at low temperature. This results in a time- and temperature-dependent enthalpy change between the excited singlet state and the charge-separated state that is the dominant term in the free energy difference between these states. Long-lived fluorescence is still observed even at 20 K, particularly for the high-potential mutants. This implies that the driving force for electron transfer on the nanosecond time scale at low temperature is less than 200 cm-1 (25 meV) in R-26 reaction centers and even smaller on the picosecond time scale or in the high-potential mutants.(ABSTRACT TRUNCATED AT 250 WORDS)

Relationship between thermodynamics and mechanism during photoinduced charge separation in reaction centers from Rhodobacter sphaeroides

Detailed fast transient absorption measurements have been performed at low temperature on reaction centers from Rhodobacter sphaeroides strain R-26 and on a double mutant, [LH(L131) + LH-(M160)], in which the P/P+ oxidation potential is roughly 140 mV (1100 cm-1) above that of wild-type reaction centers. In both samples, the decay of the excited singlet state of the initial electron donor is not well described by a single-exponential decay term. This is particularly true for reaction centers from the double mutant where at least three exponential kinetic components are required to describe the decay, with time constants ranging from a few picoseconds to hundreds of picoseconds. However, singular value decomposition analysis of the time-dependent absorption change spectra indicates the presence of only two spectrally distinct states in reaction centers from both R-26 and the double mutant. Thus, the complex decay of P* at low temperature does not appear to be due to formation of either the state P+BA- as a distinct intermediate in electron transfer or P+BB- as an equilibrated side product of electron transfer. Instead, the decay kinetics are modeled by assuming dynamic solvation of the charge-separated state, as was done for the long-lived fluorescence decay in the accompanying paper [Peloquin, J. M., Williams, J. C., Lin, X., Alden, R. G., Taguchi, A. K. W., Allen, J.P., & Woodbury, N. W. (1994) Biochemistry 33, 8089-8100]. The results of assuming a static distribution of electron-transfer rates at early times followed by dynamic solvation of the charge-separated states on longer time scales are also presented. Regardless of which model is used to describe the early time kinetics of excited-state decay, the time-dependent excited-state population on the 100-ps or longer time scale is best described in terms of thermal repopulation of P* from the charge-separated state, even at 20 K. This results in a time- and temperature-dependent driving force estimated for initial electron transfer of less than 200 cm-1 on all time scales from picoseconds to nanoseconds. Assuming a nonzero internal reorganization energy associated with charge separation, the small driving force does not appear to be consistent with the lack of temperature dependence of electron transfer and the fact that a mutant with a P/P+ oxidation potential 140 mV (1100 cm-1) higher than wild type is still able to undergo electron transfer, even at low temperature.(ABSTRACT TRUNCATED AT 400 WORDS)

Mutations designed to modify the environment of the primary electron donor of the reaction center from Rhodobacter sphaeroides: phenylalanine to leucine at L167 and histidine to phenylalanine at L168

Two mutations, L168 His to Phe and L167 Phe to Leu, were made in residues near the primary electron donor, a bacteriochlorophyll dimer, of the reaction center from Rhodobacter sphaeroides. Blue shifts of 10-15 nm in the 865-nm band of the donor were observed in the optical absorption spectra of both of the mutant reaction centers. The rate of initial electron transfer was determined by measurement of the kinetics of the decay of the excited state of the donor, and the rate of charge recombination was determined by measurement of the recovery of the bleaching of the donor. The initial electron transfer time constant and the charge recombination time constant were determined to be 3.6 ps and 220 ms, respectively, in the L168 His to Phe mutant and 5.0 ps and 85 ms in the L167 Phe to Leu mutant, compared to 3.8 ps and 100 ms measured for the wild type. The oxidation potential of the donor measured by oxidation-reduction titrations was found to decrease by 80 mV in the L168 His to Phe mutant and increase by 25 mV in the L167 Phe to Leu mutant. Time-resolved fluorescence decay measurements indicated that the change in the oxidation potential of the donor in the L168 His to Phe mutant resulted in a change in the energies of the charge-separated states. The results show that an increase in the driving force does not increase the rate of the initial electron transfer reaction.(ABSTRACT TRUNCATED AT 250 WORDS)

Effects of mutations near the bacteriochlorophylls in reaction centers from Rhodobacter sphaeroides

Mutations were made in four residues near the bacteriochlorophyll cofactors of the photosynthetic reaction center from Rhodobacter sphaeroides. These mutations, L131 Leu to His and M160 Leu to His, near the dimer bacteriochlorophylls, and M203 Gly to Asp and L177 Ile to Asp, near the monomer bacteriochlorophylls, were designed to result in the placement of a hydrogen bond donor group near the ring V keto carbonyl of each bacteriochlorophyll. Perturbations of the electronic structures of the bacteriochlorophylls in the mutants are indicated by additional resolved transitions in the bacteriochlorophyll absorption bands in steady-state low-temperature and time-resolved room temperature spectra in three of the resulting mutant reaction centers. The major effect of the two mutations near the dimer was an increase up to 80 mV in the donor oxidation-reduction midpoint potential. Correspondingly, the calculated free energy difference between the excited state of the primary donor and the initial charge separated state decreased by up to 55 mV, the initial forward electron-transfer rate was up to 4 times slower, and the rate of charge recombination between the primary quinone and the donor was approximately 30% faster in these two mutants compared to the wild type. The two mutations near the monomer bacteriochlorophylls had minor changes of 25 mV or less in the donor oxidation-reduction potential, but the mutation close to the monomer bacteriochlorophyll on the active branch resulted in a roughly 3-fold decrease in the rate of the initial electron transfer.

Biochemical characterization and electron-transfer reactions of sym1, a Rhodobacter capsulatus reaction center symmetry mutant which affects the initial electron donor

A 51 bp section of the Rhodobacter capsulatus photosynthetic reaction center M subunit gene (nucleotides M562-M612 of the pufM structural sequence) encoding amino acids M187-M203 was replaced by the homologous region of the L subunit gene. This resulted in the symmetrization of much of the amino acid environment of the reaction center initial electron donor, P. This is the first in a series of large-scale symmetry mutations and is referred to as sym1. The sym1 mutant was able to grow photosynthetically, indicating that reaction center function was largely intact. Isolated reaction centers showed an approximately 10-nm blue shift in the QY band of P. The standard free energy change between P* and P+BphA- determined from analysis of the long-lived fluorescence from quinone-reduced reaction centers decreased from about -120 meV in the wild-type to about -75 meV in the sym1 mutant. A 65-70% quantum yield of electron transfer from P* to P+QA- was observed, most of the yield loss occurring between P* and P+BphA-. The decay of the stimulated emission from P* was about 3-fold slower in this mutant than in the wild-type. Time-resolved spectral analysis of the charge-separated intermediates formed in sym1 reaction centers indicated that the major product was P+BphA-. A model-dependent analysis of the observed rates and electron-transfer yields gave the following microscopic rate constants for sym1 reaction centers (wild-type values under the same conditions are given in parentheses): [formula: see text] Analysis of the sym1 mutant, mutants near P made by other groups, and interspecies variation of amino acids in the vicinity of P suggests that the protein asymmetry in the environment of the initial electron donor is important for optimizing the rate and yield of electron transfer, but is not strictly required for overall reaction center function.

Effects of pigment-protein interactions on the conformation of the primary electron acceptor in Rhodobacter capsulatus reaction centers

Resonance Raman spectra are reported for RCs from Rb. capsulatus in which the L104 glutamic acid is replaced by glutamine. The skeletal modes of the primary electron acceptor, BPhL, in these RCs undergo temperature-dependent frequency shifts that are identical to those observed for BPhL in RCs from wild-type. This observation suggests that the strength of the hydrogen bond between the L104 residue and the C9 keto group of BPhL is not a determinant of the temperature-dependent conformation of this pigment.

Resonance Raman studies of genetically modified reaction centers from Rhodobacter capsulatus

Resonance Raman (RR) spectra are reported for the photosynthetic reaction center (RC) proteins from Rhodobacter capsulatus wild type and the genetically modified systems GluL104----Leu and HisM200----Leu. The spectra were obtained with a variety of excitation wavelengths, spanning the UV, violet, and yellow-green regions of the absorption spectrum, and at temperatures of 30 and 200 K. The RR data indicate that the structures of the bacteriochlorin pigments in RCs from Rb. capsulatus wild type are similar to those in RCs from Rhodobacter sphaeroides wild type. The data also show that the amino acid modifications near the primary electron acceptor (GluL104----Leu) and special pair (HisM200----Leu) perturb only those bacteriochlorin pigments near the site of the mutation and do not influence the structures of the other pigments in the RC. In the case of the GluL104----Leu mutant, elimination of the hydrogen bond to the C9 keto group of BPhL results in frequency shifts of RR bands of certain skeletal modes of the macrocycle. This allows the assignment of bands to the individual BPhL and BPhM pigments. In the case of the HisM200----Leu mutant, in which the special pair is comprised of a bacteriochlorophyll (BChl)-bacteriopheophytin (BPh) heterodimer rather than the BChl2 unit bound in the wild type, certain skeletal vibrations due to the additional BPh unit are identified. The frequencies of these modes are similar to those of the analogous vibrations BPhL and BPhM, which indicates that the structure of the BPh in the heterodimer is not unusual in any discernible way.(ABSTRACT TRUNCATED AT 250 WORDS)

Temperature-dependent conformational changes in the bacteriopheophytins of Rhodobacter sphaeroides reaction centers

Resonance Raman (RR) spectra are reported for the photosynthetic reaction center (RC) protein from Rhodobacter sphaeroides 2.4.1. The spectra were obtained with a variety of excitation wavelengths, spanning the UV, violet, and yellow-green regions of the absorption spectrum, and at a number of temperatures ranging from 30 to 270 K. The RR data indicate that the frequencies of certain vibrational modes of the bacteriochlorin pigments in the RC shift with temperature. These shifts are reversible and do not depend on external factors such as solvent or detergent. The acetyl carbonyl bands exhibit the largest shifts with temperature. These shifts are attributed to thermal effects involving the torsional vibrations of the acetyl groups of several (or all) of the bacteriochlorins rather than to specific pigment-protein interactions. The frequency of the structure-sensitive skeletal mode near 1610 cm-1 of one of the two bacteriopheophytins (BPhs) in the RC is also sensitive to temperature. In contrast, no temperature sensitivity is observed for the analogous modes of the bacteriochlorophylls or other BPhs. Over the range 160-100 K, the skeletal mode of the BPh upshifts by approximately 4 cm-1. This upshift is attributed to a flattening of the macrocycle at low temperatures. It is suggested that the BPh active in the electron-transfer process is the pigment whose structure is temperature dependent. It is further suggested that such structural changes could be responsible in part for the temperature dependence of the electron-transfer rates in photosynthetic RCs.