Search

Understanding the mechanisms of the cellular aging processes is crucial for attempting to extend organismal lifespan and for studying age-related degenerative diseases. Yeast cells divide through budding, providing a classical biological model for studying cellular aging. With their powerful genetics, relatively short cell cycle, and well-established signaling pathways also found in animals, yeast cells offer valuable insights into the aging process. Recent experiments suggested the existence of two aging modes in yeast characterized by nucleolar and mitochondrial declines, respectively. By analyzing experimental data, this study shows that cells evolving into those two aging modes behave differently when they are young. While buds grow linearly in both modes, cells that consistently generate spherical buds throughout their lifespan demonstrate greater efficacy in controlling bud size and growth rate at young ages. A three-dimensional multiscale chemical-mechanical model was developed and used to suggest and test hypothesized impacts of aging on bud morphogenesis. Experimentally calibrated model simulations showed that during the early stage of budding, tubular bud shape in one aging mode could be generated by locally inserting new materials at the bud tip, a process guided by the polarized Cdc42 signal. Furthermore, the aspect ratio of the tubular bud could be stabilized during the late stage as observed in experiments in this work. The model simulation results suggest that the localization of new cell surface material insertion, regulated by chemical signal polarization, could be weakened due to cellular aging in yeast and other cell types, leading to the change and stabilization of the bud aspect ratio.

In plants, the robust maintenance of tissue structure is crucial to supporting its functionality. The multi-layered shoot apical meristem (SAM) of Arabidopsis, containing stem cells, is an approximately radially symmetric tissue whose shape and structure is maintained throughout the life of the plant. In this paper, a new biologically calibrated pseudo-three-dimensional (P3D) computational model of a longitudinal section of the SAM is developed. It includes anisotropic expansion and division of cells out of the cross-section plane, as well as representation of tension experienced by the SAM epidermis. Results from the experimentally calibrated P3D model provide new insights into maintenance of the structure of the SAM epidermal cell monolayer under tension and quantify dependence of epidermal and subepidermal cell anisotropy on the amount of tension. Moreover, the model simulations revealed that out-of-plane cell growth is important in offsetting cell crowding and regulating mechanical stresses experienced by tunica cells. Predictive model simulations show that tension-determined cell division plane orientation in the apical corpus may be regulating cell and tissue shape distributions needed for maintaining structure of the wild-type SAM. This suggests that cells' responses to local mechanical cues may serve as a mechanism to regulate cell- and tissue-scale patterning.

While blood clot formation has been relatively well studied, little is known about the mechanisms underlying the subsequent structural and mechanical clot remodeling called contraction or retraction. Impairment of the clot contraction process is associated with both life-threatening bleeding and thrombotic conditions, such as ischemic stroke, venous thromboembolism, and others. Recently, blood clot contraction was observed to be hindered in patients with COVID-19. A three-dimensional multiscale computational model is developed and used to quantify biomechanical mechanisms of the kinetics of clot contraction driven by platelet-fibrin pulling interactions. These results provide important biological insights into contraction of platelet filopodia, the mechanically active thin protrusions of the plasma membrane, described previously as performing mostly a sensory function. The biomechanical mechanisms and modeling approach described can potentially apply to studying other systems in which cells are embedded in a filamentous network and exert forces on the extracellular matrix modulated by the substrate stiffness.

The exact mechanism controlling cell growth remains a grand challenge in developmental biology and regenerative medicine. The Drosophila wing disc tissue serves as an ideal biological model to study mechanisms involved in growth regulation. Most existing computational models for studying tissue growth focus specifically on either chemical signals or mechanical forces. Here we developed a multiscale chemical-mechanical model to investigate the growth regulation mechanism based on the dynamics of a morphogen gradient. By comparing the spatial distribution of dividing cells and the overall tissue shape obtained in model simulations with experimental data of the wing disc, it is shown that the size of the domain of the Dpp morphogen is critical in determining tissue size and shape. A larger tissue size with a faster growth rate and more symmetric shape can be achieved if the Dpp gradient spreads in a larger domain. Together with Dpp absorbance at the peripheral zone, the feedback regulation that downregulates Dpp receptors on the cell membrane allows for further spreading of the morphogen away from its source region, resulting in prolonged tissue growth at a more spatially homogeneous growth rate.

Gene expression specificity of homeobox transcription factors has remained paradoxical. WUSCHEL activates and represses CLAVATA3 transcription at lower and higher concentrations, respectively. We use computational modeling and experimental analysis to investigate the properties of the cis-regulatory module. We find that intrinsically each cis-element can only activate CLAVATA3 at a higher WUSCHEL concentration. However, together, they repress CLAVATA3 at higher WUSCHEL and activate only at lower WUSCHEL, showing that the concentration-dependent interactions among cis-elements regulate both activation and repression. Biochemical experiments show that two adjacent functional cis-elements bind WUSCHEL with higher affinity and dimerize at relatively lower levels. Moreover, increasing the distance between cis-elements prolongs WUSCHEL monomer binding window, resulting in higher CLAVATA3 activation. Our work showing a constellation of optimally spaced cis-elements of defined affinities determining activation and repression thresholds in regulating CLAVATA3 transcription provides a previously unknown mechanism of cofactor-independent regulation of transcription factor binding in mediating gene expression specificity.

Fibrin deformation and interaction of fibrin with other blood components play critical roles in hemostasis and thrombosis. In this review, computational and mathematical biomechanical models of fibrin network deformation and contraction at different spatio-temporal scales as well as challenges in developing and calibrating multiscale models are discussed. There are long standing challenges. For instance, applicability of models to identify and test potential mechanisms of the biomechanical processes mediating interactions between platelets and fiber networks in blood clot stretching and contraction needs to be examined carefully. How the structural and mechanical properties of major blood clot components influences biomechanical responses of the entire clot subjected to external forces, such as blood flow or vessel wall deformations needs to be investigated thoroughly.

Microtubules (MTs) are cytoskeletal fibers that undergo dynamic instability (DI), a remarkable process involving phases of growth and shortening separated by stochastic transitions called catastrophe and rescue. Dissecting DI mechanism(s) requires first characterizing and quantifying these dynamics, a subjective process that often ignores complexity in MT behavior. We present a Statistical Tool for Automated Dynamic Instability Analysis (STADIA) that identifies and quantifies not only growth and shortening, but also a category of intermediate behaviors that we term "stutters." During stutters, the rate of MT length change tends to be smaller in magnitude than during typical growth or shortening phases. Quantifying stutters and other behaviors with STADIA demonstrates that stutters precede most catastrophes in our in vitro experiments and dimer-scale MT simulations, suggesting that stutters are mechanistically involved in catastrophes. Related to this idea, we show that the anticatastrophe factor CLASP2γ works by promoting the return of stuttering MTs to growth. STADIA enables more comprehensive and data-driven analysis of MT dynamics compared with previous methods. The treatment of stutters as distinct and quantifiable DI behaviors provides new opportunities for analyzing mechanisms of MT dynamics and their regulation by binding proteins.

Stem cell maintenance in multilayered shoot apical meristems (SAMs) of plants requires strict regulation of cell growth and division. Exactly how the complex milieu of chemical and mechanical signals interact in the central region of the SAM to regulate cell division plane orientation is not well understood. In this paper, simulations using a newly developed multiscale computational model are combined with experimental studies to suggest and test three hypothesized mechanisms for the regulation of cell division plane orientation and the direction of anisotropic cell expansion in the corpus. Simulations predict that in the Apical corpus, WUSCHEL and cytokinin regulate the direction of anisotropic cell expansion, and cells divide according to tensile stress on the cell wall. In the Basal corpus, model simulations suggest dual roles for WUSCHEL and cytokinin in regulating both the direction of anisotropic cell expansion and cell division plane orientation. Simulation results are followed by a detailed analysis of changes in cell characteristics upon manipulation of WUSCHEL and cytokinin in experiments that support model predictions. Moreover, simulations predict that this layer-specific mechanism maintains both the experimentally observed shape and structure of the SAM as well as the distribution of WUSCHEL in the tissue. This provides an additional link between the roles of WUSCHEL, cytokinin, and mechanical stress in regulating SAM growth and proper stem cell maintenance in the SAM.

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell. [ABSTRACT FROM AUTHOR]

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case-mechanical or biochemical-could help define improved strategies for thrombosis prevention tailored for each patient.

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case-mechanical or biochemical-could help define improved strategies for thrombosis prevention tailored for each patient.

Regulation of the homeodomain transcription factor WUSCHEL concentration is critical for stem cell homeostasis in Arabidopsis shoot apical meristems. WUSCHEL regulates the transcription of CLAVATA3 through a concentration-dependent activation-repression switch. CLAVATA3, a secreted peptide, activates receptor kinase signaling to repress WUSCHEL transcription. Considering the revised regulation, CLAVATA3 mediated repression of WUSCHEL transcription alone will lead to an unstable system. Here we show that CLAVATA3 signaling regulates nuclear-cytoplasmic partitioning of WUSCHEL to control nuclear levels and its diffusion into adjacent cells. Our work also reveals that WUSCHEL directly interacts with EXPORTINS via EAR-like domain which is also required for destabilizing WUSCHEL in the cytoplasm. We develop a combined experimental and computational modeling approach that integrates CLAVATA3-mediated transcriptional repression of WUSCHEL and post-translational control of nuclear levels with the WUSCHEL concentration-dependent regulation of CLAVATA3. We show that the dual control by the same signal forms a seamless connection between de novo WUSCHEL synthesis and sub-cellular partitioning in providing robustness to the WUSCHEL gradient.

Budding yeast, Saccharomyces cerevisiae, serves as a prime biological model to study mechanisms underlying asymmetric growth. Previous studies have shown that prior to bud emergence, polarization of a conserved small GTPase Cdc42 must be established on the cell membrane of a budding yeast. Additionally, such polarization contributes to the delivery of cell wall remodeling enzymes and hydrolase from cytosol through the membrane, to change the mechanical properties of the cell wall. This leads to the hypothesis that Cdc42 and its associated proteins at least indirectly regulate cell surface mechanical properties. However, how the surface mechanical properties in the emerging bud are changed and whether such change is important are not well understood. To test several hypothesised mechanisms, a novel three-dimensional coarse-grained particle-based model has been developed which describes inhomogeneous mechanical properties of the cell surface. Model simulations predict alternation of the levels of stretching and bending stiffness of the cell surface in the bud region by the polarized Cdc42 signals is essential for initiating bud formation. Model simulations also suggest that bud shape depends strongly on the distribution of the polarized signaling molecules while the neck width of the emerging bud is strongly impacted by the mechanical properties of the chitin and septin rings. Moreover, the temporal change of the bud mechanical properties is shown to affect the symmetry of the bud shape. The 3D model of asymmetric cell growth can also be used for studying viral budding and other vegetative reproduction processes performed via budding, as well as detailed studies of cell growth.

Experimental measurements or computational model predictions of the post-translational regulation of enzymes needed in a metabolic pathway is a difficult problem. Consequently, regulation is mostly known only for well-studied reactions of central metabolism in various model organisms. In this study, we use two approaches to predict enzyme regulation policies and investigate the hypothesis that regulation is driven by the need to maintain the solvent capacity in the cell. The first predictive method uses a statistical thermodynamics and metabolic control theory framework while the second method is performed using a hybrid optimization-reinforcement learning approach. Efficient regulation schemes were learned from experimental data that either agree with theoretical calculations or result in a higher cell fitness using maximum useful work as a metric. As previously hypothesized, regulation is herein shown to control the concentrations of both immediate and downstream product concentrations at physiological levels. Model predictions provide the following two novel general principles: (1) the regulation itself causes the reactions to be much further from equilibrium instead of the common assumption that highly non-equilibrium reactions are the targets for regulation; and (2) the minimal regulation needed to maintain metabolite levels at physiological concentrations maximizes the free energy dissipation rate instead of preserving a specific energy charge. The resulting energy dissipation rate is an emergent property of regulation which may be represented by a high value of the adenylate energy charge. In addition, the predictions demonstrate that the amount of regulation needed can be minimized if it is applied at the beginning or branch point of a pathway, in agreement with common notions. The approach is demonstrated for three pathways in the central metabolism of E. coli (gluconeogenesis, glycolysis-tricarboxylic acid (TCA) and pentose phosphate-TCA) that each require different regulation schemes. It is shown quantitatively that hexokinase, glucose 6-phosphate dehydrogenase and glyceraldehyde phosphate dehydrogenase, all branch points of pathways, play the largest roles in regulating central metabolism.

Epithelial sheets define organ architecture during development. Here, we employed an iterative multiscale computational modeling and quantitative experimental approach to decouple direct and indirect effects of actomyosin-generated forces, nuclear positioning, extracellular matrix, and cell-cell adhesion in shaping Drosophila wing imaginal discs. Basally generated actomyosin forces generate epithelial bending of the wing disc pouch. Surprisingly, acute pharmacological inhibition of ROCK-driven actomyosin contractility does not impact the maintenance of tissue height or curved shape. Computational simulations show that ECM tautness provides only a minor contribution to modulating tissue shape. Instead, passive ECM pre-strain serves to maintain the shape independent from actomyosin contractility. These results provide general insight into how the subcellular forces are generated and maintained within individual cells to induce tissue curvature. Thus, the results suggest an important design principle of separable contributions from ECM prestrain and actomyosin tension during epithelial organogenesis and homeostasis.

The concept of critical concentration (CC) is central to understanding the behavior of microtubules (MTs) and other cytoskeletal polymers. Traditionally, these polymers are understood to have one CC, measured in multiple ways and assumed to be the subunit concentration necessary for polymer assembly. However, this framework does not incorporate dynamic instability (DI), and there is work indicating that MTs have two CCs. We use our previously established simulations to confirm that MTs have (at least) two experimentally relevant CCs and to clarify the behavior of individuals and populations relative to the CCs. At free subunit concentrations above the lower CC (CCElongation), growth phases of individual filaments can occur transiently; above the higher CC (CCNetAssembly), the population's polymer mass will increase persistently. Our results demonstrate that most experimental CC measurements correspond to CCNetAssembly, meaning that "typical" DI occurs below the concentration traditionally considered necessary for polymer assembly. We report that [free tubulin] at steady state does not equal CCNetAssembly, but instead approaches CCNetAssembly asymptotically as [total tubulin] increases, and depends on the number of stable MT nucleation sites. We show that the degree of separation between CCElongation and CCNetAssembly depends on the rate of nucleotide hydrolysis. This clarified framework helps explain and unify many experimental observations.

Platelets play a key role in the formation of hemostatic clots and obstructive thrombi as well as in other biological processes. In response to physiological stimulants, including thrombin, platelets change shape, express adhesive molecules, aggregate, and secrete bioactive substances, but their subsequent fate is largely unknown. Here we examined late-stage structural, metabolic, and functional consequences of thrombin-induced platelet activation. Using a combination of confocal microscopy, scanning and transmission electron microscopy, flow cytometry, biochemical and biomechanical measurements, we showed that thrombin-induced activation is followed by time-dependent platelet dysfunction and disintegration. After ~30 minutes of incubation with thrombin, unlike with collagen or ADP, human platelets disintegrated into cellular fragments containing organelles, such as mitochondria, glycogen granules, and vacuoles. This platelet fragmentation was preceded by Ca2+ influx, integrin αIIbβ3 activation and phosphatidylserine exposure (activation phase), followed by mitochondrial depolarization, generation of reactive oxygen species, metabolic ATP depletion and impairment of platelet contractility along with dramatic cytoskeletal rearrangements, concomitant with platelet disintegration (death phase). Coincidentally with the platelet fragmentation, thrombin caused calpain activation but not activation of caspases 3 and 7. Our findings indicate that the late functional and structural damage of thrombin-activated platelets comprise a calpain-dependent platelet death pathway that shares some similarities with the programmed death of nucleated cells, but is unique to platelets, therefore representing a special form of cellular destruction. Fragmentation of activated platelets suggests that there is an underappreciated pathway of enhanced elimination of platelets from the circulation in (pro)thrombotic conditions once these cells have performed their functions.

Fibrin is a viscoelastic proteinaceous polymer that determines the deformability and integrity of blood clots and fibrin-based biomaterials in response to biomechanical forces. Here, a previously unnoticed structural mechanism of fibrin clots' mechanical response to external tensile loads is tested using high-resolution confocal microscopy and recently developed three-dimensional computational model. This mechanism, underlying local strain-stiffening of individual fibers as well as global stiffening of the entire network, is based on previously neglected nascent cohesive pairwise interactions between individual fibers (crisscrossing) in fibrin networks formed under tensile load. Existence of fiber-fiber crisscrossings of reoriented fibers was confirmed using 3D imaging of experimentally obtained stretched fibrin clots. The computational model enabled us to study structural details and quantify mechanical effects of the fiber-fiber cohesive crisscrossing during stretching of fibrin gels at various spatial scales. The contribution of the fiber-fiber cohesive contacts to the elasticity of stretched fibrin networks was characterized by changes in individual fiber stiffness, the length, width, and alignment of fibers, as well as connectivity and density of the entire network. The results show that the nascent cohesive crisscrossing of fibers in stretched fibrin networks comprise an underappreciated important structural mechanism underlying the mechanical response of fibrin to (patho)physiological stresses that determine the course and outcomes of thrombotic and hemostatic disorders, such as heart attack and ischemic stroke. STATEMENT OF SIGNIFICANCE: Fibrin is a viscoelastic proteinaceous polymer that determines the deformability and integrity of blood clots and fibrin-based biomaterials in response to biomechanical forces. In this paper, a novel structural mechanism of fibrin clots' mechanical response to external tensile loads is tested using high-resolution confocal microscopy and newly developed computational model. This mechanism, underlying local strain-stiffening of individual fibers as well as global stiffening of the entire network, is based on previously neglected nascent cohesive pairwise interactions between individual fibers (crisscrossing) in fibrin networks formed under tensile load. Cohesive crisscrossing is an important structural mechanism that influences the mechanical response of blood clots and which can determine the outcomes of blood coagulation disorders, such as heart attacks and strokes.

Fueled by breakthrough technology developments, the biological, biomedical, and behavioral sciences are now collecting more data than ever before. There is a critical need for time- and cost-efficient strategies to analyze and interpret these data to advance human health. The recent rise of machine learning as a powerful technique to integrate multimodality, multifidelity data, and reveal correlations between intertwined phenomena presents a special opportunity in this regard. However, machine learning alone ignores the fundamental laws of physics and can result in ill-posed problems or non-physical solutions. Multiscale modeling is a successful strategy to integrate multiscale, multiphysics data and uncover mechanisms that explain the emergence of function. However, multiscale modeling alone often fails to efficiently combine large datasets from different sources and different levels of resolution. Here we demonstrate that machine learning and multiscale modeling can naturally complement each other to create robust predictive models that integrate the underlying physics to manage ill-posed problems and explore massive design spaces. We review the current literature, highlight applications and opportunities, address open questions, and discuss potential challenges and limitations in four overarching topical areas: ordinary differential equations, partial differential equations, data-driven approaches, and theory-driven approaches. Towards these goals, we leverage expertise in applied mathematics, computer science, computational biology, biophysics, biomechanics, engineering mechanics, experimentation, and medicine. Our multidisciplinary perspective suggests that integrating machine learning and multiscale modeling can provide new insights into disease mechanisms, help identify new targets and treatment strategies, and inform decision making for the benefit of human health.

Energetic Variational Approach is used to derive a novel thermodynamically consistent three-phase model of a mixture of Newtonian and visco-elastic fluids. The model which automatically satisfies the energy dissipation law and is Galilean invariant, consists of coupled Navier-Stokes and Cahn-Hilliard equations. Modified General Navier Boundary Condition with fluid elasticity taken into account is also introduced for using the model to study moving contact line problems. Energy stable numerical scheme is developed to solve system of model equations efficiently. Convergence of the numerical scheme is verified by simulating a droplet sliding on an inclined plane under gravity. The model can be applied for studying various biological or biophysical problems. Predictive abilities of the model are demonstrated by simulating deformation of venous blood clots with different visco-elastic properties and experimentally observed internal structures under different biologically relevant shear blood flow conditions.

Blood clot contraction plays an important role in prevention of bleeding and in thrombotic disorders. Here, we unveil and quantify the structural mechanisms of clot contraction at the level of single platelets. A key elementary step of contraction is sequential extension-retraction of platelet filopodia attached to fibrin fibers. In contrast to other cell-matrix systems in which cells migrate along fibers, the "hand-over-hand" longitudinal pulling causes shortening and bending of platelet-attached fibers, resulting in formation of fiber kinks. When attached to multiple fibers, platelets densify the fibrin network by pulling on fibers transversely to their longitudinal axes. Single platelets and aggregates use actomyosin contractile machinery and integrin-mediated adhesion to remodel the extracellular matrix, inducing compaction of fibrin into bundled agglomerates tightly associated with activated platelets. The revealed platelet-driven mechanisms of blood clot contraction demonstrate an important new biological application of cell motility principles.

Thromboembolism, one of the leading causes of morbidity and mortality worldwide, is characterized by formation of obstructive intravascular clots (thrombi) and their mechanical breakage (embolization). A novel two-dimensional multi-phase computational model is introduced that describes active interactions between the main components of the clot, including platelets and fibrin, to study the impact of various physiologically relevant blood shear flow conditions on deformation and embolization of a partially obstructive clot with variable permeability. Simulations provide new insights into mechanisms underlying clot stability and embolization that cannot be studied experimentally at this time. In particular, model simulations, calibrated using experimental intravital imaging of an established arteriolar clot, show that flow-induced changes in size, shape and internal structure of the clot are largely determined by two shear-dependent mechanisms: reversible attachment of platelets to the exterior of the clot and removal of large clot pieces. Model simulations predict that blood clots with higher permeability are more prone to embolization with enhanced disintegration under increasing shear rate. In contrast, less permeable clots are more resistant to rupture due to shear rate-dependent clot stiffening originating from enhanced platelet adhesion and aggregation. These results can be used in future to predict risk of thromboembolism based on the data about composition, permeability and deformability of a clot under specific local haemodynamic conditions.

The formation of platelet thrombi is determined by the integrin αIIbβ3-mediated interactions of platelets with fibrinogen and fibrin. Blood clotting in vivo is catalyzed by thrombin, which simultaneously induces fibrinogen binding to αIIbβ3 and converts fibrinogen to fibrin. Thus, after a short time, thrombus formation is governed by αIIbβ3 binding to fibrin fibers. Surprisingly, there is little understanding of αIIbβ3 interaction with fibrin polymers. Here we used an optical trap-based system to measure the binding of single αIIbβ3 molecules to polymeric fibrin and compare it to αIIbβ3 binding to monomeric fibrin and fibrinogen. Like αIIbβ3 binding to fibrinogen and monomeric fibrin, we found that αIIbβ3 binding to polymeric fibrin can be segregated into two binding regimes, one with weaker rupture forces of 30-60 pN and a second with stronger rupture forces >60 pN that peaked at 70-80 pN. However, we found that the mechanical stability of the bimolecular αIIbβ3-ligand complexes had the following order: fibrin polymer > fibrin monomer > fibrinogen. These quantitative differences reflect the distinct specificity and underlying molecular mechanisms of αIIbβ3-mediated reactions, implying that targeting platelet interactions with fibrin could increase the therapeutic indices of antithrombotic agents by focusing on the destabilization of thrombi rather than the prevention of platelet aggregation.

A technological revolution in both light and electron microscopy imaging now allows unprecedented views of clotting, especially in animal models of hemostasis and thrombosis. However, our understanding of three-dimensional high-resolution clot structure remains incomplete since most of our recent knowledge has come from studies of relatively small clots or thrombi, due to the optical impenetrability of clots beyond a few cell layers in depth. Here, we developed an optimized optical clearing method termed cCLOT that renders large whole blood clots transparent and allows confocal imaging as deep as one millimeter inside the clot. We have tested this method by investigating the 3D structure of clots made from reconstituted pre-labeled blood components yielding new information about the effects of clot contraction on erythrocytes. Although it has been shown recently that erythrocytes are compressed to form polyhedrocytes during clot contraction, observations of this phenomenon have been impeded by the inability to easily image inside clots. As an efficient and non-destructive method, cCLOT represents a powerful research tool in studying blood clot structure and mechanisms controlling clot morphology. Additionally, cCLOT optical clearing has the potential to facilitate imaging of ex vivo clots and thrombi derived from healthy or pathological conditions.

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.

Swarming groups of bacteria coordinate their behavior by self-organizing as a population to move over surfaces in search of nutrients and optimal niches for colonization. Many open questions remain about the cues used by swarming bacteria to achieve this self-organization. While chemical cue signaling known as quorum sensing is well-described, swarming bacteria often act and coordinate on time scales that could not be achieved via these extracellular quorum sensing cues. Here, cell-cell contact-dependent protein exchange is explored as a mechanism of intercellular signaling for the bacterium Myxococcus xanthus. A detailed biologically calibrated computational model is used to study how M. xanthus optimizes the connection rate between cells and maximizes the spread of an extracellular protein within the population. The maximum rate of protein spreading is observed for cells that reverse direction optimally for swarming. Cells that reverse too slowly or too fast fail to spread extracellular protein efficiently. In particular, a specific range of cell reversal frequencies was observed to maximize the cell-cell connection rate and minimize the time of protein spreading. Furthermore, our findings suggest that predesigned motion reversal can be employed to enhance the collective behavior of biological synthetic active systems.

A novel multi-component model is introduced for studying interaction between blood flow and deforming aortic wall with intramural hematoma (IMH). The aortic wall is simulated by a composite structure submodel representing material properties of the three main wall layers. The IMH is described by a poroelasticity submodel which takes into account both the pressure inside hematoma and its deformation. The submodel of the hematoma is fully coupled with the aortic submodel as well as with the submodel of the pulsatile blood flow. Model simulations are used to investigate the relation between the peak wall stress, hematoma thickness and permeability in patients of different age. The results indicate that an increase in hematoma thickness leads to larger wall stress, which is in agreement with clinical data. Further simulations demonstrate that a hematoma with smaller permeability results in larger wall stress, suggesting that blood coagulation in hematoma might increase its mechanical stability. This is in agreement with previous experimental observations of coagulation having a beneficial effect on the condition of a patient with the IMH.

Cell segmentation and motion tracking in time-lapse images are fundamental problems in computer vision, and are also crucial for various biomedical studies. Myxococcus xanthus is a type of rod-like cells with highly coordinated motion. The segmentation and tracking of M. xanthus are challenging, because cells may touch tightly and form dense swarms that are difficult to identify individually in an accurate manner. The known cell tracking approaches mainly fall into two frameworks, detection association and model evolution, each having its own advantages and disadvantages. In this paper, we propose a new hybrid framework combining these two frameworks into one and leveraging their complementary advantages. Also, we propose an active contour model based on the Ribbon Snake, which is seamlessly integrated with our hybrid framework. Evaluated by 10 different datasets, our approach achieves considerable improvement over the state-of-the-art cell tracking algorithms on identifying complete cell trajectories, and higher segmentation accuracy than performing segmentation in individual 2D images.

The rheological properties of fibrin networks have been of long-standing interest. As such there is a wealth of studies of their shear and tensile responses, but their compressive behavior remains unexplored. Here, by characterization of the network structure with synchronous measurement of the fibrin storage and loss moduli at increasing degrees of compression, we show that the compressive behavior of fibrin networks is similar to that of cellular solids. A nonlinear stress-strain response of fibrin consists of three regimes: (1) an initial linear regime, in which most fibers are straight, (2) a plateau regime, in which more and more fibers buckle and collapse, and (3) a markedly nonlinear regime, in which network densification occurs by bending of buckled fibers and inter-fiber contacts. Importantly, the spatially non-uniform network deformation included formation of a moving "compression front" along the axis of strain, which segregated the fibrin network into compartments with different fiber densities and structure. The Young's modulus of the linear phase depends quadratically on the fibrin volume fraction while that in the densified phase depends cubically on it. The viscoelastic plateau regime corresponds to a mixture of these two phases in which the fractions of the two phases change during compression. We model this regime using a continuum theory of phase transitions and analytically predict the storage and loss moduli which are in good agreement with the experimental data. Our work shows that fibrin networks are a member of a broad class of natural cellular materials which includes cancellous bone, wood and cork.

A novel biofilm model is described which systemically couples bacteria, extracellular polymeric substances (EPS) and solvent phases in biofilm. This enables the study of contributions of rheology of individual phases to deformation of biofilm in response to fluid flow as well as interactions between different phases. The model, which is based on first and second laws of thermodynamics, is derived using an energetic variational approach and phase-field method. Phase-field coupling is used to model structural changes of a biofilm. A newly developed unconditionally energy-stable numerical splitting scheme is implemented for computing the numerical solution of the model efficiently. Model simulations predict biofilm cohesive failure for the flow velocity between [Formula: see text] and [Formula: see text] m s(-1) which is consistent with experiments. Simulations predict biofilm deformation resulting in the formation of streamers for EPS exhibiting a viscous-dominated mechanical response and the viscosity of EPS being less than [Formula: see text]. Higher EPS viscosity provides biofilm with greater resistance to deformation and to removal by the flow. Moreover, simulations show that higher EPS elasticity yields the formation of streamers with complex geometries that are more prone to detachment. These model predictions are shown to be in qualitative agreement with experimental observations.

Bacterial surface motility, such as swarming, is commonly examined in the laboratory using plate assays that necessitate specific concentrations of agar and sometimes inclusion of specific nutrients in the growth medium. The preparation of such explicit media and surface growth conditions serves to provide the favorable conditions that allow not just bacterial growth but coordinated motility of bacteria over these surfaces within thin liquid films. Reproducibility of swarm plate and other surface motility plate assays can be a major challenge. Especially for more "temperate swarmers" that exhibit motility only within agar ranges of 0.4%-0.8% (wt/vol), minor changes in protocol or laboratory environment can greatly influence swarm assay results. "Wettability", or water content at the liquid-solid-air interface of these plate assays, is often a key variable to be controlled. An additional challenge in assessing swarming is how to quantify observed differences between any two (or more) experiments. Here we detail a versatile two-phase protocol to prepare and image swarm assays. We include guidelines to circumvent the challenges commonly associated with swarm assay media preparation and quantification of data from these assays. We specifically demonstrate our method using bacteria that express fluorescent or bioluminescent genetic reporters like green fluorescent protein (GFP), luciferase (lux operon), or cellular stains to enable time-lapse optical imaging. We further demonstrate the ability of our method to track competing swarming species in the same experiment.

Pseudomonas aeruginosa is a ubiquitous bacterium that survives in many environments, including as an acute and chronic pathogen in humans. Substantial evidence shows that P. aeruginosa behavior is affected by its motility, and appendages known as flagella and type IV pili (TFP) are known to confer such motility. The role these appendages play when not facilitating motility or attachment, however, is unclear. Here we discern a passive intercellular role of TFP during flagellar-mediated swarming of P. aeruginosa that does not require TFP extension or retraction. We studied swarming at the cellular level using a combination of laboratory experiments and computational simulations to explain the resultant patterns of cells imaged from in vitro swarms. Namely, we used a computational model to simulate swarming and to probe for individual cell behavior that cannot currently be otherwise measured. Our simulations showed that TFP of swarming P. aeruginosa should be distributed all over the cell and that TFP-TFP interactions between cells should be a dominant mechanism that promotes cell-cell interaction, limits lone cell movement, and slows swarm expansion. This predicted physical mechanism involving TFP was confirmed in vitro using pairwise mixtures of strains with and without TFP where cells without TFP separate from cells with TFP. While TFP slow swarm expansion, we show in vitro that TFP help alter collective motion to avoid toxic compounds such as the antibiotic carbenicillin. Thus, TFP physically affect P. aeruginosa swarming by actively promoting cell-cell association and directional collective motion within motile groups to aid their survival.

Links between cell division and other cellular processes are poorly understood. It is difficult to simultaneously examine division and function in most cell types. Most of the research probing aspects of cell division has experimented with stationary or immobilized cells or distinctly asymmetrical cells. Here we took an alternative approach by examining cell division events within motile groups of cells growing on solid medium by time-lapse microscopy. A total of 558 cell divisions were identified among approximately 12,000 cells. We found an interconnection of division, motility, and polarity in the bacterium Myxococcus xanthus. For every division event, motile cells stop moving to divide. Progeny cells of binary fission subsequently move in opposing directions. This behavior involves M. xanthus Frz proteins that regulate M. xanthus motility reversals but is independent of type IV pilus "S motility." The inheritance of opposing polarity is correlated with the distribution of the G protein RomR within these dividing cells. The constriction at the point of division limits the intracellular distribution of RomR. Thus, the asymmetric distribution of RomR at the parent cell poles becomes mirrored at new poles initiated at the site of division.

When a blood vessel ruptures or gets inflamed, the human body responds by rapidly forming a clot to restrict the loss of blood. Platelets aggregation at the injury site of the blood vessel occurring via platelet-platelet adhesion, tethering and rolling on the injured endothelium is a critical initial step in blood clot formation. A novel three-dimensional multi-scale model is introduced and used in this paper to simulate receptor-mediated adhesion of deformable platelets at the site of vascular injury under different shear rates of blood flow. The novelty of the model is based on a new approach of coupling submodels at three biological scales crucial for the early clot formation: novel hybrid cell membrane submodel to represent physiological elastic properties of a platelet, stochastic receptor-ligand binding submodel to describe cell adhesion kinetics and lattice Boltzmann submodel for simulating blood flow. The model implementation on the GPU cluster significantly improved simulation performance. Predictive model simulations revealed that platelet deformation, interactions between platelets in the vicinity of the vessel wall as well as the number of functional GPIbα platelet receptors played significant roles in platelet adhesion to the injury site. Variation of the number of functional GPIbα platelet receptors as well as changes of platelet stiffness can represent effects of specific drugs reducing or enhancing platelet activity. Therefore, predictive simulations can improve the search for new drug targets and help to make treatment of thrombosis patient-specific.

Fibrin is a protein polymer that forms a 3D filamentous network, a major structural component of protective physiological blood clots as well as life threatening pathological thrombi. It plays an important role in wound healing, tissue regeneration and is widely employed in surgery as a sealant and in tissue engineering as a scaffold. The goal of this study was to establish correlations between structural changes and mechanical responses of fibrin networks exposed to compressive loads. Rheological measurements revealed nonlinear changes of fibrin network viscoelastic properties under dynamic compression, resulting in network softening followed by its dramatic hardening. Repeated compression/decompression enhanced fibrin clot stiffening. Combining fibrin network rheology with simultaneous confocal microscopy provided direct evidence of structural modulations underlying nonlinear viscoelasticity of compressed fibrin networks. Fibrin clot softening in response to compression strongly correlated with fiber buckling and bending, while hardening was associated with fibrin network densification. Our results suggest a complex interplay of entropic and enthalpic mechanisms accompanying structural changes and accounting for the nonlinear mechanical response in fibrin networks undergoing compressive deformations. These findings provide new insight into the fibrin clot structural mechanics and can be useful for designing fibrin-based biomaterials with modulated viscoelastic properties.