A Systems Approach to Biomechanics, Mechanobiology, and Biotransport

The human body represents a collection of interacting systems that range in scale from nanometers to meters. Investigations from a systems perspective focus on how the parts work together to enact changes across spatial scales, and further our understanding of how systems function and fail. Here, we highlight systems approaches presented at the 2022 Summer Biomechanics, Bio-engineering, and Biotransport Conference in the areas of solid mechanics; fluid mechanics; tissue and cellular engineering; biotransport; and design, dynamics, and rehabilitation; and biomechanics education. Systems approaches are yielding new insights into human biology by leveraging state-of-the-art tools, which could ultimately lead to more informed design of therapies and medical devices for preventing and treating disease as well as rehabilitating patients using strategies that are uniquely optimized for each patient. Educational approaches can also be designed to foster a foundation of systems-level thinking.

Effects of Porosity on Piezoelectric Characteristics of Polyvinylidene Fluoride Films for Biomedical Applications

Objective: The objective of this work is to study the effects of porosity on mechanical and piezoelectric properties of polyvinylidene fluoride (PVDF) films for biomedical applications. Impact Statement: By investigating the piezoelectric properties of PVDF and the porosity effect on its electromechanical performance, there is potential for further development of PVDF as a hemodynamic sensor that can lead to further technological advancements in the biomedical field, benefiting patients and physicians alike. Introduction: PVDF thin films have shown potential in the application of hemodynamic flow sensing and monitoring the effects on blood flow caused by prosthetic valve implantation via the transcatheter aortic valve replacement operation. The piezoelectric performance of PVDF films can be influenced by the porosity of the material. Methods: In this study, strain tracking was performed on thin film PVDF specimens with various levels of porosity and pore sizes to determine the mechanical properties of the specimens. The mechanical properties were used to model the PVDF material in COMSOL multiphysics software, in which compression test simulations were performed to determine the piezoelectric coefficient d33 of the PVDF. Results: A decline in the elastic modulus was found to be highly inversely correlated with porosity of the specimens and the simulation results show that elastic modulus had a much greater effect on the piezoelectric properties than Poisson's ratio. Conclusion: A combination of experimental and computational techniques was able to characterize and correlate the mechanical properties of PVDF films of varying porosities to their piezoelectric properties.

Elucidating the signal for contact guidance contained in aligned fibrils with a microstructural-mechanical model

Despite its importance in physiological processes and tissue engineering, the mechanism underlying cell contact guidance in an aligned fibrillar network has defied elucidation due to multiple interdependent signals that such a network presents to cells, namely, anisotropy of adhesion, porosity and mechanical behaviour. A microstructural-mechanical model of fibril networks was used to assess the relative magnitudes of these competing signals in networks of varied alignment strength based on idealized cylindrical pseudopods projected into the aligned and orthogonal directions and computing the anisotropy of metrics chosen for adhesion, porosity and mechanical behaviour: cylinder-fibre contact area for adhesion, persistence length of pores for porosity and total force to displace fibres from the cylindrical volume as well as network stiffness experienced upon cylinder retraction for mechanical behaviour. The signals related to mechanical anisotropy are substantially higher than adhesion and porosity anisotropy, especially at stronger network alignments, although their signal to noise (S/N) values are substantially lower. The former finding is consistent with a recent report that fibroblasts can sense fibril alignment via anisotropy of network mechanical resistance, and the model reveals this can be due to either mechanical resistance to pseudopod protrusion or retraction given their signal and S/N values are similar.

Responses of ticks to immersion in hot bathing water: Effect of surface type, water temperature, and soap on tick motor control

Preventing bites from undetected ticks through bathing practices would benefit public health, but the effects of these practices have been researched minimally. We immersed nymphal and adult hard ticks of species common in the eastern United States in tap water, using temperatures and durations that are realistic for human hot bathing. The effect of (a) different skin-equivalent surfaces (silicone and pig skin), and (b) water temperature was tested on Amblyomma americanum, Dermacentor variabilis and Ixodes scapularis nymphs. Overall, the type of surface had a much larger effect on the nymphs’ tendency to stay in contact with the surface than water temperature did. Most nymphs that separated from the surface did so within the first 10 s of immersion, with the majority losing contact due to the formation of an air bubble between their ventral side and the test surface. In addition, adult Ixodes scapularis were tested for the effect of immersion time, temperature, and soap on tick responsiveness. Some individual adults moved abnormally or stopped moving as a result of longer or hotter immersion, but soap had little effect on responsiveness. Taken together, our results suggest that the surface plays a role in ticks’ tendency to stay in contact; the use of different bath additives warrants further research. While water temperature did not have a significant short-term effect on tick separation, ticks that have not attached by their mouth parts may be rendered unresponsive and eventually lose contact with a person’s skin in a hot bath. It should be noted that our research did not consider potential temperature effects on the pathogens themselves, as previous research suggests that some tickborne pathogens may become less hazardous even if the tick harboring them survives hot-water exposures and later bites the bather after remaining undetected.

Pandemic-Driven Online Teaching-The Natural Setting for a Flipped Classroom?

As the COVID-19 pandemic forced a sudden shift to online teaching and learning in April 2020, one of the more significant challenges faced by instructors is encouraging and maintaining student engagement in their online classes. This paper describes my experience of flipping an online classroom for a core Chemical Engineering Fluid Mechanics class to promote student engagement and collaboration in an online setting. Comparing exam scores with prior semesters involving in-person, traditional lecture-style classes suggests that students need a certain degree of adjustment to adapt to this new learning mode. A decrease in student rating of teaching (SRT) scores indicates that students largely prefer in-person, traditional lectures over an online flipped class, even though written comments in the SRT contained several responses favorable to flipping the class in an online setting. Overall, SRT scores on a department level also showed a similar decrease, which suggests students were less satisfied with the quality of teaching overall throughout the department, with this flipped method of instruction neither improving nor worsening student sentiment toward online learning. In addition, whereas most students liked the prerecorded lecture videos, they were less enthusiastic about using breakout rooms to encourage student collaboration and discussion. Further thought and discussion on best practices to facilitate online student interaction and collaboration are recommended, as online learning will likely continue to grow in popularity even when in-person instruction resumes after the pandemic.

Swelling of Collagen-Hyaluronic Acid Co-Gels: An In Vitro Residual Stress Model

Tissue-equivalents (TEs), simple model tissues with tunable properties, have been used to explore many features of biological soft tissues. Absent in most formulations however, is the residual stress that arises due to interactions among components with different unloaded levels of stress, which has an important functional role in many biological tissues. To create a pre-stressed model system, co-gels were fabricated from a combination of hyaluronic acid (HA) and reconstituted Type-I collagen (Col). When placed in solutions of varying osmolarity, HA-Col co-gels swell as the HA imbibes water, which in turn stretches (and stresses) the collagen network. In this way, co-gels with residual stress (i.e., collagen fibers in tension and HA in compression) were fabricated. When the three gel types tested here were immersed in hypotonic solutions, pure HA gels swelled the most, followed by HA-Col co-gels; no swelling was observed in pure collagen gels. The greatest swelling rates and swelling ratios occurred in the lowest salt concentration solutions. Tension on the collagen component of HA-Col co-gels was calculated from a stress balance and increased nonlinearly as swelling increased. The swelling experiment results were in good agreement with the stress predicted by a fibril network + non-fibrillar interstitial matrix computational model.

Multiscale Mechanical Model of the Pacinian Corpuscle Shows Depth and Anisotropy Contribute to the Receptor's Characteristic Response to Indentation

Cutaneous mechanoreceptors transduce different tactile stimuli into neural signals that produce distinct sensations of touch. The Pacinian corpuscle (PC), a cutaneous mechanoreceptor located deep within the dermis of the skin, detects high frequency vibrations that occur within its large receptive field. The PC is comprised of lamellae that surround the nerve fiber at its core. We hypothesized that a layered, anisotropic structure, embedded deep within the skin, would produce the nonlinear strain transmission and low spatial sensitivity characteristic of the PC. A multiscale finite-element model was used to model the equilibrium response of the PC to indentation. The first simulation considered an isolated PC with fiber networks aligned with the PC's surface. The PC was subjected to a 10 μm indentation by a 250 μm diameter indenter. The multiscale model captured the nonlinear strain transmission through the PC, predicting decreased compressive strain with proximity to the receptor's core, as seen experimentally by others. The second set of simulations considered a single PC embedded epidermally (shallow) or dermally (deep) to model the PC's location within the skin. The embedded models were subjected to 10 μm indentations at a series of locations on the surface of the skin. Strain along the long axis of the PC was calculated after indentation to simulate stretch along the nerve fiber at the center of the PC. Receptive fields for the epidermis and dermis models were constructed by mapping the long-axis strain after indentation at each point on the surface of the skin mesh. The dermis model resulted in a larger receptive field, as the calculated strain showed less indenter location dependence than in the epidermis model.

A coupled fiber-matrix model demonstrates highly inhomogeneous microstructural interactions in soft tissues under tensile load

A soft tissue's macroscopic behavior is largely determined by its microstructural components (often a collagen fiber network surrounded by a nonfibrillar matrix (NFM)). In the present study, a coupled fiber-matrix model was developed to fully quantify the internal stress field within such a tissue and to explore interactions between the collagen fiber network and nonfibrillar matrix (NFM). Voronoi tessellations (representing collagen networks) were embedded in a continuous three-dimensional NFM. Fibers were represented as one-dimensional nonlinear springs and the NFM, meshed via tetrahedra, was modeled as a compressible neo-Hookean solid. Multidimensional finite element modeling was employed in order to couple the two tissue components and uniaxial tension was applied to the composite representative volume element (RVE). In terms of the overall RVE response (average stress, fiber orientation, and Poisson's ratio), the coupled fiber-matrix model yielded results consistent with those obtained using a previously developed parallel model based upon superposition. The detailed stress field in the composite RVE demonstrated the high degree of inhomogeneity in NFM mechanics, which cannot be addressed by a parallel model. Distributions of maximum/minimum principal stresses in the NFM showed a transition from fiber-dominated to matrix-dominated behavior as the matrix shear modulus increased. The matrix-dominated behavior also included a shift in the fiber kinematics toward the affine limit. We conclude that if only gross averaged parameters are of interest, parallel-type models are suitable. If, however, one is concerned with phenomena, such as individual cell-fiber interactions or tissue failure that could be altered by local variations in the stress field, then the detailed model is necessary in spite of its higher computational cost.

A multiscale approach to modeling the passive mechanical contribution of cells in tissues

In addition to their obvious biological roles in tissue function, cells often play a significant mechanical role through a combination of passive and active behaviors. This study focused on the passive mechanical contribution of cells in tissues by improving our multiscale model via the addition of cells, which were treated as dilute spherical inclusions. The first set of simulations considered a rigid cell, with the surrounding ECM modeled as (1) linear elastic, (2) Neo-Hookean, and (3) a fiber network. Comparison with the classical composite theory for rigid inclusions showed close agreement at low cell volume fraction. The fiber network case exhibited nonlinear stress-strain behavior and Poisson's ratios larger than the elastic limit of 0.5, characteristics similar to those of biological tissues. The second set of simulations used a fiber network for both the cell (simulating cytoskeletal filaments) and matrix, and investigated the effect of varying relative stiffness between the cell and matrix, as well as the effect of a cytoplasmic pressure to enforce incompressibility of the cell. Results showed that the ECM network exerted negligible compression on the cell, even when the stiffness of fibers in the network was increased relative to the cell. Introduction of a cytoplasmic pressure significantly increased the stresses in the cell filament network, and altered how the cell changed its shape under tension. Findings from this study have implications on understanding how cells interact with their surrounding ECM, as well as in the context of mechanosensation.

Mechanical behavior of collagen-fibrin co-gels reflects transition from series to parallel interactions with increasing collagen content

Fibrin and collagen, biopolymers occurring naturally in the body, are biomaterials commonly-used as scaffolds for tissue engineering. How collagen and fibrin interact to confer macroscopic mechanical properties in collagen-fibrin composite systems remains poorly understood. In this study, we formulated collagen-fibrin co-gels at different collagen-to-fibrin ratios to observe changes in the overall mechanical behavior and microstructure. A modeling framework of a two-network system was developed by modifying our micro-scale model, considering two forms of interaction between the networks: (a) two interpenetrating but noninteracting networks ("parallel"), and (b) a single network consisting of randomly alternating collagen and fibrin fibrils ("series"). Mechanical testing of our gels show that collagen-fibrin co-gels exhibit intermediate properties (UTS, strain at failure, tangent modulus) compared to those of pure collagen and fibrin. The comparison with model predictions show that the parallel and series model cases provide upper and lower bounds, respectively, for the experimental data, suggesting that a combination of such interactions exists between the collagen and fibrin in co-gels. A transition from the series model to the parallel model occurs with increasing collagen content, with the series model best describing predominantly fibrin co-gels, and the parallel model best describing predominantly collagen co-gels.

Mechanics of a fiber network within a non-fibrillar matrix: model and comparison with collagen-agarose co-gels

While collagen is recognized as the predominant mechanical component of soft connective tissues, the role of the non-fibrillar matrix (NFM) is less well understood. Even model systems, such as the collagen-agarose co-gel, can exhibit complex behavior, making it difficult to identify relative contributions of specific tissue constituents. In the present study, we developed a two-component microscale model of collagen-agarose tissue analogs and used it to elucidate the interaction between collagen and NFM in uniaxial tension. Collagen fibers were represented with Voronoi networks, and the NFM was modeled as a neo-Hookean solid. Model predictions of total normal stress and Poisson's ratio matched experimental observations well (including high Poisson's values of ~3), and the addition of NFM led to composition-dependent decreases in volume change and increases in fiber stretch. Because the NFM was more resistant to volume change than the fiber network, extension of the composite led to pressurization of the NFM. Within a specific range of parameter values (low shear modulus and moderate Poisson's ratio), the magnitude of the reaction force decreased relative to this pressurization component resulting in a negative (compressive) NFM stress in the loading direction, even though the composite tissue was in tension.

Identification of ubiquitin ligases required for skeletal muscle atrophy

Skeletal muscle adapts to decreases in activity and load by undergoing atrophy. To identify candidate molecular mediators of muscle atrophy, we performed transcript profiling. Although many genes were up-regulated in a single rat model of atrophy, only a small subset was universal in all atrophy models. Two of these genes encode ubiquitin ligases: Muscle RING Finger 1 (MuRF1), and a gene we designate Muscle Atrophy F-box (MAFbx), the latter being a member of the SCF family of E3 ubiquitin ligases. Overexpression of MAFbx in myotubes produced atrophy, whereas mice deficient in either MAFbx or MuRF1 were found to be resistant to atrophy. These proteins are potential drug targets for the treatment of muscle atrophy.

Accelerated mammary tumor development in mutant polyomavirus middle T transgenic mice expressing elevated levels of either the Shc or Grb2 adapter protein

The Grb2 and Shc adapter proteins play critical roles in coupling activated growth factor receptors to several cellular signaling pathways. To assess the role of these molecules in mammary epithelial development and tumorigenesis, we have generated transgenic mice which individually express the Grb2 and Shc proteins in the mammary epithelium. Although mammary epithelial cell-specific expression of Grb2 or Shc accelerated ductal morphogenesis, mammary tumors were rarely observed in these strains. To explore the potential role of these adapter proteins in mammary tumorigenesis, mice coexpressing either Shc or Grb2 and a mutant form of polyomavirus middle T (PyV mT) antigen in the mammary epithelium were generated. Coexpression of either Shc or Grb2 with the mutant PyV mT antigen resulted in a dramatic acceleration of mammary tumorigenesis compared to parental mutant PyV mT strain. The increased rate of tumor formation observed in these mice was correlated with activation of the epidermal growth factor receptor family and mitogen-activated protein kinase pathway. These observations suggest that elevated levels of the Grb2 or Shc adapter protein can accelerate mammary tumor progression by sensitizing the mammary epithelial cell to growth factor receptor signaling.

Identification of residues that control specific binding of the Shc phosphotyrosine-binding domain to phosphotyrosine sites

The Shc adaptor protein contains two phosphotyrosine [Tyr(P)]binding modules--an N-terminal Tyr(P) binding (PTB) domain and a C-terminal Src homology 2 (SH2) domain. We have compared the ability of the Shc PTB domain to bind the receptors for nerve growth factor and insulin, both of which contain juxtamembrane Asn-Pro-Xaa-Tyr(P) motifs implicated in PTB binding. The Shc PTB domain binds with high affinity to a phosphopeptide corresponding to the nerve growth factor receptor Tyr-490 autophosphorylation site. Analysis of individual residues within this motif indicates that the Asn at position -3 [with respect to Tyr(P)], in addition to Tyr(P), is critical for PTB binding, while the Pro at position -2 plays a less significant role. A hydrophobic amino acid 5 residues N-terminal to the Tyr(P) is also essential for high-affinity binding. In contrast, the Shc PTB domain does not bind stably to the Asn-Pro-Xaa-Tyr(P) site at Tyr-960 in the activated insulin receptor, which has a polar residue (Ser) at position -5. Substitution of this Ser at position -5 with Ile markedly increased binding of the insulin receptor Tyr-960 phosphopeptide to the PTB domain. These results suggest that while the Shc PTB domain recognizes a core sequence of Asn-Pro-Xaa-Tyr(P), its binding affinity is modulated by more N-terminal residues in the ligand, which therefore contribute to the specificity of PTB-receptor interactions. An analysis of residues in the Shc PTB domain required for binding to Tyr(P) sites identified a specific and evolutionarily conserved Arg (Arg-175) that is uniquely important for ligand binding and is potentially involved in Tyr(P) recognition.

A conserved amino-terminal Shc domain binds to phosphotyrosine motifs in activated receptors and phosphopeptides

BACKGROUND

Signal transduction by growth factor receptor protein-tyrosine kinases is generally initiated by autophosphorylation on tyrosine residues following ligand binding. Phosphotyrosines within activated receptors form binding sites for the Src homology 2 (SH2) domains of cytoplasmic signalling proteins. One such protein, Shc, is tyrosine phosphorylated in response to a large number of growth factors and cytokines. Phosphorylation of Shc on tyrosine residue Y317 allows binding to the SH2 domain of Grb2, and hence stimulation of the Ras pathway. Shc is therefore implicated as an adaptor protein able to couple normal and oncogenic protein-tyrosine kinases to Ras activation. Shc itself contains an SH2 domain at its carboxyl terminus, but the function of the amino-terminal half of the protein is unknown.

RESULTS

We have found that the Shc amino-terminal region binds to a number of tyrosine-phosphorylated proteins in v-src-transformed cells. This domain also bound directly to the activated epidermal growth factor (EGF) receptor. A phosphotyrosine (pY)-containing peptide modeled after the Shc-binding site in polyoma middle T antigen (LLSNPTpYSVMRSK) was able to compete efficiently with the activated EGF receptor for binding to the Shc amino terminus. This competition was dependent on phosphorylation of the tyrosine residue within the peptide, and was abrogated by deletion of the leucine residue at position -5. The Shc amino-terminal domain also bound to the autophosphorylated nerve growth factor receptor (Trk), but bound significantly less well to a mutant receptor in which tyrosine Y490 in the receptor's Shc-binding site had been substituted by phenylalanine.

CONCLUSION

These data implicate the amino-terminal region of Shc in binding to activated receptors and other tyrosine-phosphorylated proteins. Binding appears to be specific for phosphorylated tyrosine residues within the sequence NPXpY, which is conserved in many Shc-binding sites. The Shc amino-terminal region bears only very limited sequence identify to known SH2 domains, suggesting that it represents a new class of phosphotyrosine-binding modules. Consistent with this view, the amino-terminal Shc domain is highly conserved in a Drosophila Shc homologue. Binding of Shc to activated receptors through its amino terminus could leave the carboxy-terminal SH2 domain free for other interactions. In this way, Shc may function as an adaptor protein to bring two tyrosine-phosphorylated proteins together.