Contributions of Glycosaminoglycans to Collagen Fiber Recruitment in Constitutive Modeling of Arterial Mechanics

The contribution of glycosaminoglycans (GAGs) to the biological and mechanical functions of biological tissue has emerged as an important area of research. GAGs provide structural basis for the organization and assembly of extracellular matrix (ECM). The mechanics of tissue with low GAG content can be indirectly affected by the interaction of GAGs with collagen fibers, which have long been known to be one of the primary contributors to soft tissue mechanics. Our earlier study showed that enzymatic GAG depletion results in straighter collagen fibers that are recruited at lower levels of stretch, and a corresponding shift in earlier arterial stiffening (Mattson et al., 2016). In this study, the effect of GAGs on collagen fiber recruitment was studied through a structure-based constitutive model. The model incorporates structural information, such as fiber orientation distribution, content, and recruitment of medial elastin, medial collagen, and adventitial collagen fibers. The model was first used to study planar biaxial tensile stress-stretch behavior of porcine descending thoracic aorta. Changes in elastin and collagen fiber orientation distribution, and collagen fiber recruitment were then incorporated into the model in order to predict the stress-stretch behavior of GAG depleted tissue. Our study shows that incorporating early collagen fiber recruitment into the model predicts the stress-stretch response of GAG depleted tissue reasonably well (rms = 0.141); considering further changes of fiber orientation distribution does not improve the predicting capability (rms = 0.149). Our study suggests an important role of GAGs in arterial mechanics that should be considered in developing constitutive models.     

Mathematical modelling of bone adaptation of the metacarpal subchondral bone in racehorses

In Thoroughbred racehorses, fractures of the distal limb are commonly catastrophic. Most of these fractures occur due to the accumulation of fatigue damage from repetitive loading, as evidenced by microdamage at the predilection sites for fracture. Adaptation of the bone in response to training loads is important for fatigue resistance. In order to better understand the mechanism of subchondral bone adaptation to its loading environment, we utilised a square root function defining the relationship between bone volume fraction \((f_{BM} )\) and specific surface \((S_v )\) of the subchondral bone of the lateral condyles of the third metacarpal bone (MCIII) of the racehorse, and using this equation, developed a mathematical model of subchondral bone that adapts to loading conditions observed in vivo. The model is expressed as an ordinary differential equation incorporating a formation rate that is dependent on strain energy density. The loading conditions applied to a selected subchondral region, i.e. volume of interest, were estimated based on joint contact forces sustained by racehorses in training. For each of the initial conditions of \(f_{BM} \) we found no difference between subsequent homoeostatic \(f_{BM} \) at any given loading condition, but the time to reach equilibrium differed by initial \(f_{BM} \) and loading condition. We found that the observed values for \(f_{BM} \) from the mathematical model output were a good approximation to the existing data for racehorses in training or at rest. This model provides the basis for understanding the effect of changes to training strategies that may reduce the risk of racehorse injury.     

Arterial Extracellular Matrix: A Mechanobiological Study of the Contributions and Interactions of Elastin and Collagen

The complex network structure of elastin and collagen extracellular matrix (ECM) forms the primary load bearing components in the arterial wall. The structural and mechanobiological interactions between elastin and collagen are important for properly functioning arteries. Here, we examined the elastin and collagen organization, realignment, and recruitment by coupling mechanical loading and multiphoton imaging. Two-photon excitation fluorescence and second harmonic generation methods were performed with a multiphoton video-rate microscope to capture real time changes to the elastin and collagen structure during biaxial deformation. Enzymatic removal of elastin was performed to assess the structural changes of the remaining collagen structure. Quantitative analysis of the structural changes to elastin and collagen was made using a combination of two-dimensional fast Fourier transform and fractal analysis, which allows for a more complete understanding of structural changes. Our study provides new quantitative evidence, to our knowledge on the sequential engagement of different arterial ECM components in response to mechanical loading. The adventitial collagen exists as large wavy bundles of fibers that exhibit fiber engagement after 20% strain. The medial collagen is engaged throughout the stretching process, and prominent elastic fiber engagement is observed up to 20% strain after which the engagement plateaus. The fiber orientation distribution functions show remarkably different changes in the ECM structure in response to mechanical loading. The medial collagen shows an evident preferred circumferential distribution, however the fiber families of adventitial collagen are obscured by their waviness at no or low mechanical strains. Collagen fibers in both layers exhibit significant realignment in response to unequal biaxial loading. The elastic fibers are much more uniformly distributed and remained relatively unchanged due to loading. Removal of elastin produces similar structural changes in collagen as mechanical loading. Our study suggests that the elastic fibers are under tension and impart an intrinsic compressive stress on the collagen.     

Arterial Extracellular Matrix: A Mechanobiological Study of the Contributions and Interactions of Elastin and Collagen

The complex network structure of elastin and collagen extracellular matrix (ECM) forms the primary load bearing components in the arterial wall. The structural and mechanobiological interactions between elastin and collagen are important for properly functioning arteries. Here, we examined the elastin and collagen organization, realignment, and recruitment by coupling mechanical loading and multiphoton imaging. Two-photon excitation fluorescence and second harmonic generation methods were performed with a multiphoton video-rate microscope to capture real time changes to the elastin and collagen structure during biaxial deformation. Enzymatic removal of elastin was performed to assess the structural changes of the remaining collagen structure. Quantitative analysis of the structural changes to elastin and collagen was made using a combination of two-dimensional fast Fourier transform and fractal analysis, which allows for a more complete understanding of structural changes. Our study provides new quantitative evidence, to our knowledge on the sequential engagement of different arterial ECM components in response to mechanical loading. The adventitial collagen exists as large wavy bundles of fibers that exhibit fiber engagement after 20% strain. The medial collagen is engaged throughout the stretching process, and prominent elastic fiber engagement is observed up to 20% strain after which the engagement plateaus. The fiber orientation distribution functions show remarkably different changes in the ECM structure in response to mechanical loading. The medial collagen shows an evident preferred circumferential distribution, however the fiber families of adventitial collagen are obscured by their waviness at no or low mechanical strains. Collagen fibers in both layers exhibit significant realignment in response to unequal biaxial loading. The elastic fibers are much more uniformly distributed and remained relatively unchanged due to loading. Removal of elastin produces similar structural changes in collagen as mechanical loading. Our study suggests that the elastic fibers are under tension and impart an intrinsic compressive stress on the collagen.     

Toward translating near-infrared spectroscopy oxygen saturation data for the non-invasive prediction of spatial and temporal hemodynamics during exercise

Image-based computational fluid dynamics (CFD) studies conducted at rest have shown that atherosclerotic plaque in the thoracic aorta (TA) correlates with adverse wall shear stress (WSS), but there is a paucity of such data under elevated flow conditions. We developed a pedaling exercise protocol to obtain phase contrast magnetic resonance imaging (PC-MRI) blood flow measurements in the TA and brachiocephalic arteries during three-tiered supine pedaling at 130, 150, and 170 % of resting heart rate (HR), and relate these measurements to non-invasive tissue oxygen saturation \((\hbox {StO}_{2})\) acquired by near-infrared spectroscopy (NIRS) while conducting the same protocol. Local quantification of WSS indices by CFD revealed low time-averaged WSS on the outer curvature of the ascending aorta and the inner curvature of the descending aorta (dAo) that progressively increased with exercise, but that remained low on the anterior surface of brachiocephalic arteries. High oscillatory WSS observed on the inner curvature of the aorta persisted during exercise as well. Results suggest locally continuous exposure to potentially deleterious indices of WSS despite benefits of exercise. Linear relationships between flow distributions and tissue oxygen extraction calculated from \(\hbox {StO}_{2}\) were found between the left common carotid versus cerebral tissue \((r^{2}=0.96)\) and the dAo versus leg tissue \((r^{2}=0.87)\). A resulting six-step procedure is presented to use NIRS data as a surrogate for exercise PC-MRI when setting boundary conditions for future CFD studies of the TA under simulated exercise conditions. Relationships and ensemble-averaged PC-MRI inflow waveforms are provided in an online repository for this purpose.     

Contribution of left ventricular residual stress by myocytes and collagen: existence of inter-constituent mechanical interaction

We quantify the contribution of myocytes, collagen fibers and their interactions to the residual stress field found in the left ventricle (LV) using both experimental and theoretical methods. Ring tissue samples extracted from normal rat, male and female, LV were treated with collagenase and decellularization to isolate myocytes and collagen fibers, respectively. Opening angle tests were then performed on these samples as well as intact tissue samples containing both constituents that served as control. Our results show that the collagen fibers are the main contributor to the residual stress fields found in the LV. Specifically, opening angle measured in collagen-only samples (106.45\(^\circ \) ± 23.02\(^\circ \)) and myocytes-only samples (21.00\(^\circ \) ± 4.37\(^\circ \)) was significantly higher and lower than that of the control (57.88\(^\circ \) ± 12.29\(^\circ \)), respectively. A constrained mixture (CM) modeling framework was then used to infer these experimental results. We show that the framework cannot reproduce the opening angle found in the intact tissue with measurements made on the collagen-only and myocytes-only samples. Given that the CM framework assumes that each constituent contributes to the overall mechanics simply by their mere presence, this result suggests the existence of some myocyte–collagen mechanical interaction that cannot be ignored in the LV. We then propose an extended CM formulation that takes into account of the inter-constituent mechanical interaction in which constituents are deformed additionally when they are physically combined into a mixture. We show that the intact tissue opening angle can be recovered in this framework.     

Local distribution of collagen fibers determines crack initiation site and its propagation direction during aortic rupture

Although elucidation of the mechanism of aortic aneurysm rupture is important, the characteristics of crack initiation and propagation sites remain unknown. To determine the microscopic properties of these sites, the characteristics of local strains and constituents at crack initiation and propagation sites were investigated during biaxial stretching of porcine thoracic aortas (PTAs). PTAs were sliced into approximately 50-\(\upmu \hbox {m}\)-thick sections, and the center of the sections was made especially thin using our previously developed technique. Alpha-elastin and cell nuclei were fluorescently labeled as indices of local elastin density and as a strain marker, respectively. Birefringence and second harmonic generation (SHG) light images were used to determine local collagen distributions. The specimens were then stretched biaxially with a laboratory-made tensile tester under a fluorescent microscope equipped with a birefringence imaging system. Local strains were calculated from the local displacement of the cell nuclei. The degree of alignment and density of local collagen fibers were measured from retardance and SHG images. The strain distributions, specifically the first and second principal, and maximum shear strains, fluorescent intensity of \(\upalpha \)-elastin, and degree of alignment of collagen fibers, showed insignificant differences between the crack initiation sites and other sites. The retardance and intensity of SHG light at the crack initiation sites were significantly lower than those at other sites for all (\(n = 6\)) specimens. Cracks tended to propagate along the local direction of the collagen fibers. These results indicate that the local density and direction of collagen fibers play an important role in aorta rupture.     

A finite-element model of mechanosensation by a Pacinian corpuscle cluster in human skin

The Pacinian corpuscle (PC) is the cutaneous mechanoreceptor responsible for sensation of high-frequency (20–1000 Hz) vibrations. PCs lie deep within the skin, often in multicorpuscle clusters with overlapping receptive fields. We developed a finite-element mechanical model of one or two PCs embedded within human skin, coupled to a multiphysics PC model to simulate action potentials elicited by each PC. A vibration was applied to the skin surface, and the resulting mechanical signal was analyzed using two metrics: the deformation amplitude ratio (\({\rho }_{\mathrm{1S}} \), \({\rho }_{\mathrm{2S}} )\) and the phase shift of the vibration (\({\delta }_{\mathrm{S}1}^{\mathrm{mech}} \), \({\delta }_{\mathrm{S}2}^{\mathrm{mech}} )\) between the stimulus and the PC. Our results showed that the amplitude attenuation and phase shift at a PC increased with distance from the stimulus to the PC. Differences in amplitude (\(\rho _{12} )\) and phase shift (\({\delta }_{12}^{\mathrm{mech}} )\) between the two PCs in simulated clusters directly affected the interspike interval between the action potentials elicited by each PC (\({\delta }_{12}^{\mathrm{spike}} )\). While \({\delta }_{12}^{\mathrm{mech}} \) had a linear relationship with \({\delta }_{12}^{\mathrm{spike}} \), \(\rho _{12} \)’s effect on \({\delta }_{12}^{\mathrm{spike}} \) was greater for lower values of \(\rho _{12} \). In our simulations, the separation between PCs and the distance of each PC from the stimulus location resulted in differences in amplitude and phase shift at each PC that caused \({\delta }_{12}^{\mathrm{spike}} \) to vary with PC location. Our results suggest that PCs within a cluster receive different mechanical stimuli which may enhance source localization of vibrotactile stimuli, drawing parallels to sound localization in binaural hearing.     

Uncertainty quantification and sensitivity analysis of an arterial wall mechanics model for evaluation of vascular drug therapies

Quantification of the uncertainty in constitutive model predictions describing arterial wall mechanics is vital towards non-invasive assessment of vascular drug therapies. Therefore, we perform uncertainty quantification to determine uncertainty in mechanical characteristics describing the vessel wall response upon loading. Furthermore, a global variance-based sensitivity analysis is performed to pinpoint measurements that are most rewarding to be measured more precisely. We used previously published carotid diameter–pressure and intima–media thickness (IMT) data (measured in triplicate), and Holzapfel–Gasser–Ogden models. A virtual data set containing 5000 diastolic and systolic diameter–pressure points, and IMT values was generated by adding measurement error to the average of the measured data. The model was fitted to single-exponential curves calculated from the data, obtaining distributions of constitutive parameters and constituent load bearing parameters. Additionally, we (1) simulated vascular drug treatment to assess the relevance of model uncertainty and (2) evaluated how increasing the number of measurement repetitions influences model uncertainty. We found substantial uncertainty in constitutive parameters. Simulating vascular drug treatment predicted a 6% point reduction in collagen load bearing (\(L_\mathrm {coll}\)), approximately 50% of its uncertainty. Sensitivity analysis indicated that the uncertainty in \(L_{\mathrm {coll}}\) was primarily caused by noise in distension and IMT measurements. Spread in \(L_{\mathrm {coll}}\) could be decreased by 50% when increasing the number of measurement repetitions from 3 to 10. Model uncertainty, notably that in \(L_{\mathrm {coll}}\), could conceal effects of vascular drug therapy. However, this uncertainty could be reduced by increasing the number of measurement repetitions of distension and wall thickness measurements used for model parameterisation.     

A Mathematical Model for the Macrophage Response to Respiratory Viral Infection in Normal and Asthmatic Conditions

Respiratory viral infections are common in the general population and one of the most important causes of asthma aggravation and exacerbation. Despite many studies, it is not well understood how viral infections cause more severe symptoms and exacerbations in asthmatics. We develop a mathematical model of two types of macrophages that play complementary roles in fighting viral infection: classically \((\hbox {CA}\)-\(\hbox {M}\Phi )\) and alternatively activated macrophages \((\hbox {AA}\)-\(\hbox {M}\Phi )\). \(\hbox {CA}\)-\(\hbox {M}\Phi \) destroy infected cells and tissues to remove viruses, while \(\hbox {AA}\)-\(\hbox {M}\Phi \) repair damaged tissues. We show that a higher viral load or longer duration of infection provokes a stronger immune response from the macrophage system. By adjusting the parameters, we model the differences in response to respiratory viral infection in normal and asthmatic subjects and show how this skews the system toward a response that generates more severe symptoms in asthmatic patients.     

Multiphoton microscopy observations of 3D elastin and collagen fiber microstructure changes during pressurization in aortic media

Elastin and collagen fibers play important roles in the mechanical properties of aortic media. Because knowledge of local fiber structures is required for detailed analysis of blood vessel wall mechanics, we investigated 3D microstructures of elastin and collagen fibers in thoracic aortas and monitored changes during pressurization. Using multiphoton microscopy, autofluorescence images from elastin and second harmonic generation signals from collagen were acquired in media from rabbit thoracic aortas that were stretched biaxially to restore physiological dimensions. Both elastin and collagen fibers were observed in all longitudinal–circumferential plane images, whereas alternate bright and dark layers were observed along the radial direction and were recognized as elastic laminas (ELs) and smooth muscle-rich layers (SMLs), respectively. Elastin and collagen fibers are mainly oriented in the circumferential direction, and waviness of collagen fibers was significantly higher than that of elastin fibers. Collagen fibers were more undulated in longitudinal than in radial direction, whereas undulation of elastin fibers was equibiaxial. Changes in waviness of collagen fibers during pressurization were then evaluated using 2-dimensional fast Fourier transform in mouse aortas, and indices of waviness of collagen fibers decreased with increases in intraluminal pressure. These indices also showed that collagen fibers in SMLs became straight at lower intraluminal pressures than those in EL, indicating that SMLs stretched more than ELs. These results indicate that deformation of the aorta due to pressurization is complicated because of the heterogeneity of tissue layers and differences in elastic properties of ELs, SMLs, and surrounding collagen and elastin.     

Elasticity–density and viscoelasticity–density relationships at the tibia mid-diaphysis assessed from resonant ultrasound spectroscopy measurements

Cortical bone tissue is an anisotropic material characterized by typically five independent elastic coefficients (for transverse isotropy) governing shear and longitudinal deformations in the different anatomical directions. It is well established that the Young’s modulus in the direction of the bone axis of long bones has a strong relationship with mass density. It is not clear, however, whether relationships of similar strength exist for the other elastic coefficients, for they have seldom been investigated, and the results available in the literature are contradictory. The objectives of the present work were to document the anisotropic elastic properties of cortical bone at the tibia mid-diaphysis and to elucidate their relationships with mass density. Resonant ultrasound spectroscopy (RUS) was used to measure the transverse isotropic stiffness tensor of 55 specimens from 19 donors. Except for Poisson’s ratios and the non-diagonal stiffness coefficient, strong linear correlations between the different elastic coefficients \((0.7 < {r^{2}} < 0.99)\) and between these coefficients and density \((0.79 < {r^{2}} < 0.89)\) were found. Comparison with previously published data from femur specimens suggested that the strong correlations evidenced in this study may not only be valid for the mid-tibia. RUS also measures the viscous part of the stiffness tensor. An anisotropy ratio close to two was found for damping coefficients. Damping increased as the mass density decreased. The data suggest that a relatively accurate estimation of all the mid-tibia elastic coefficients can be derived from mass density. This is of particular interest (1) to design organ-scale bone models in which elastic coefficients are mapped according to Hounsfield values from computed tomography scans as a surrogate for mass density and (2) to model ultrasound propagation at the mid-tibia, which is an important site for the in vivo assessment of bone status with axial transmission techniques.     

Tissue growth constrained by extracellular matrix drives invagination during optic cup morphogenesis

In the early embryo, the eyes form initially as relatively spherical optic vesicles (OVs) that protrude from both sides of the brain tube. Each OV grows until it contacts and adheres to the overlying surface ectoderm (SE) via an extracellular matrix (ECM) that is secreted by the SE and OV. The OV and SE then thicken and bend inward (invaginate) to create the optic cup (OC) and lens vesicle, respectively. While constriction of cell apices likely plays a role in SE invagination, the mechanisms that drive OV invagination are poorly understood. Here, we used experiments and computational modeling to explore the hypothesis that the ECM locally constrains the growing OV, forcing it to invaginate. In chick embryos, we examined the need for the ECM by (1) removing SE at different developmental stages and (2) exposing the embryo to collagenase. At relatively early stages of invagination (Hamburger–Hamilton stage HH14\(-\)), removing the SE caused the curvature of the OV to reverse as it ‘popped out’ and became convex, but the OV remained concave at later stages (HH15) and invaginated further during subsequent culture. Disrupting the ECM had a similar effect, with the OV popping out at early to mid-stages of invagination (HH14\(-\) to HH14\(+\)). These results suggest that the ECM is required for the early stages but not the late stages of OV invagination. Microindentation tests indicate that the matrix is considerably stiffer than the cellular OV, and a finite-element model consisting of a growing spherical OV attached to a relatively stiff layer of ECM reproduced the observed behavior, as well as measured temporal changes in OV curvature, wall thickness, and invagination depth reasonably well. Results from our study also suggest that the OV grows relatively uniformly, while the ECM is stiffer toward the center of the optic vesicle. These results are consistent with our matrix-constraint hypothesis, providing new insight into the mechanics of OC (early retina) morphogenesis.     

A fiber-based constitutive model predicts changes in amount and organization of matrix proteins with development and disease in the mouse aorta

Decreased elastin in mice (Eln+/?) yields a functioning vascular system with elevated blood pressure and increased arterial stiffness that is morphologically distinct from wild-type mice (WT). Yet, function is retained enough that there is no appreciable effect on life span and some mechanical properties are maintained constant. It is not understood how the mouse modifies the normal developmental process to produce a functioning vascular system despite a deficiency in elastin. To quantify changes in mechanical properties, we have applied a fiber-based constitutive model to mechanical data from the ascending aorta during postnatal development of WT and Eln+/? mice. Results indicate that the fiber-based constitutive model is capable of distinguishing elastin amounts and identifying trends during development. We observe an increase in predicted circumferential stress contribution from elastin with age, which correlates with increased elastin amounts from protein quantification data. The model also predicts changes in the unloaded collagen fiber orientation with age, which must be verified in future work. In Eln+/? mice, elastin amounts are decreased at each age, along with the predicted circumferential stress contribution of elastin. Collagen amounts in Eln+/? aorta are comparable to WT, but the predicted circumferential stress contribution of collagen is increased. This may be due to altered organization or structure of the collagen fibers. Relating quantifiable changes in arterial mechanics with changes in extracellular matrix (ECM) protein amounts will help in understanding developmental remodeling and in producing treatments for human diseases affecting ECM proteins.     

Fluorescence microscopic distinction between elastin and collagen

Summary Distinction between elastin and collagen in arteriosclerotic lesions is difficult because immature and incompletely cross-linked collagen bind so-called elastica stains; furthermore, abnormal collagen can lack cross-striation and thus resemble elastin in electron microscopy. However, collagen and elastin differ significantly in their content of basic amino acids and hence in their affinity for heteropolyacids. This chemical difference was utilized for the development of a fluorescence microscopic method for distinction between collagen and elastin.Paraffin sections of human autopsy material were treated with a 1% aqueous solution of phosphomolybdic acid (PMA) for five minutes, rinsed in distilled water, dehydrated and mounted. Other series were treated with the PMA-molybdenum blue reaction and with various special stains.Only elastic membranes of aorta, the elastica interna and externa of sizable arteries, and true elastic fibers remained strongly fluorescent; the autofluorescence of collagen, reticulum fibers, basement membranes, pseudo-elastic fibers, and elastic membranes in small arteries was quenched. In other series PMA abolished the fluorescence of basic fluorochromes.Correlation of fluorescence and direct light microscopic observations with chemical and electron microscopic data showed that the PMA-fluorescence method permits distinction between elastin and various types of collagen.     

An improved Four-Russians method and sparsified Four-Russians algorithm for RNA folding

Background

The basic RNA secondary structure prediction problem or single sequence folding problem (SSF) was solved 35 years ago by a now well-known \(O(n^3)\)-time dynamic programming method. Recently three methodologies—Valiant, Four-Russians, and Sparsification—have been applied to speedup RNA secondary structure prediction. The sparsification method exploits two properties of the input: the number of subsequence Z with the endpoints belonging to the optimal folding set and the maximum number base-pairs L. These sparsity properties satisfy \(0 \le L \le n / 2\) and \(n \le Z \le n^2 / 2\), and the method reduces the algorithmic running time to O(LZ). While the Four-Russians method utilizes tabling partial results.

Results

In this paper, we explore three different algorithmic speedups. We first expand the reformulate the single sequence folding Four-Russians \(\Theta \left(\frac{n^3}{\log ^2 n}\right)\)-time algorithm, to utilize an on-demand lookup table. Second, we create a framework that combines the fastest Sparsification and new fastest on-demand Four-Russians methods. This combined method has worst-case running time of \(O(\tilde{L}\tilde{Z})\), where \(\frac{{L}}{\log n} \le \tilde{L}\le min\left({L},\frac{n}{\log n}\right)\) and \(\frac{{Z}}{\log n}\le \tilde{Z} \le min\left({Z},\frac{n^2}{\log n}\right)\). Third we update the Four-Russians formulation to achieve an on-demand \(O( n^2/ \log ^2n )\)-time parallel algorithm. This then leads to an asymptotic speedup of \(O(\tilde{L}\tilde{Z_j})\) where \(\frac{{Z_j}}{\log n}\le \tilde{Z_j} \le min\left({Z_j},\frac{n}{\log n}\right)\) and \(Z_j\) the number of subsequence with the endpoint j belonging to the optimal folding set.

Conclusions

The on-demand formulation not only removes all extraneous computation and allows us to incorporate more realistic scoring schemes, but leads us to take advantage of the sparsity properties. Through asymptotic analysis and empirical testing on the base-pair maximization variant and a more biologically informative scoring scheme, we show that this Sparse Four-Russians framework is able to achieve a speedup on every problem instance, that is asymptotically never worse, and empirically better than achieved by the minimum of the two methods alone.

     

Extracellular matrix and the mechanics of large artery development

The large, elastic arteries, as their name suggests, provide elastic distention and recoil during the cardiac cycle in vertebrate animals. The arteries are distended from the pressure of ejecting blood during the active contraction of the left ventricle (LV) during systole and recoil to their original dimensions during relaxation of the LV during diastole. The cyclic distension occurs with minimal energy loss, due to the elastic properties of one of the major structural extracellular matrix (ECM) components, elastin. The maximum distension is limited to prevent damage to the artery by another major ECM component, collagen. The mix of ECM components in the wall largely determines the passive mechanical behavior of the arteries and the subsequent load on the heart during systole. While much research has focused on initial artery formation, there has been less attention on the continuing development of the artery to produce the mature composite wall complete with endothelial cells (ECs), smooth muscle cells (SMCs), and the necessary mix of ECM components for proper cardiovascular function. This review focuses on the physiology of large artery development, including SMC differentiation and ECM production. The effects of hemodynamic forces and ECM deposition on the evolving arterial structure and function are discussed. Human diseases and mouse models with genetic mutations in ECM proteins that affect large artery development are summarized. A review of constitutive models and growth and remodeling theories is presented, along with future directions to improve understanding of ECM and the mechanics of large artery development.     

Effect of collagen fibre orientation on intervertebral disc torsion mechanics

The intervertebral disc is a complex fibro-cartilaginous material, consisting of a pressurized nucleus pulposus surrounded by the annulus fibrosus, which has an angle-ply structure. Disc injury and degeneration are noted by significant changes in tissue structure and function, which significantly alters stress distribution and disc joint stiffness. Differences in fibre orientation are thought to contribute to changes in disc torsion mechanics. Therefore, the objective of this study was to evaluate the effect of collagen fibre orientation on internal disc mechanics under compression combined with axial rotation. We developed and validated a finite element model (FEM) to delineate changes in disc mechanics due to fibre orientation from differences in material properties. FEM simulations were performed with fibres oriented at \(\pm 30^{\circ }\) throughout the disc (uniform by region and fibre layer). The initial model was validated by published experimental results for two load conditions, including \(0.48\,\hbox {MPa}\) axial compression and \(10\,\hbox {Nm}\) axial rotation. Once validated, fibre orientation was rotated by \(4^{\circ }\) or \(8^{\circ }\) towards the horizontal plane, resulting in a decrease in disc joint torsional stiffness. Furthermore, we observed that axial rotation caused a sinusoidal change in disc height and radial bulge, which may be beneficial for nutrient transport. In conclusion, including anatomically relevant fibre angles in disc joint FEMs is important for understanding stress distribution throughout the disc and will be important for understanding potential causes for disc injury. Future models will include regional differences in fibre orientation to better represent the fibre architecture of the native disc.