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.     

Trend analysis of variations in carbon stock using stock big data

Changes in land use affect the terrestrial carbon stock through changes in the land cover. Research on land use and analysis of variations in carbon stock have practical applications in the optimization of land use and the mitigation of climate change effects. This study was conducted in Baixiang and Julu counties in the Taihang Piedmont by employing the trend analysis method to characterize the variation in county land use and carbon stock. The findings show that in both counties, agricultural and unused land areas decreased while built-up land area increased, and the reduction in cropland was the main reason behind the agricultural land reduction. An inflection point appeared on the cropland curves of Julu, because the cropland area decreased by 1576.97 hm\(^{2}\) from 2004 to 2006. Cropland area in Baixiang decreased from 1996 to 1998 by a total of 129.89 hm\(^{2}\) and then remained relatively stable after 1998. The total carbon storage and variation in land use in the two counties displayed similar trends. Total carbon reserves in Julu increased by 2.76 \(\times \) 10\(^{4}\) tC (carbon equivalent), while those in Baixiang decreased by 0.63 \(\times \) 10\(^{4}\) tC. Carbon stock of built-up land in Julu and Baixiang increased by 2.44 \(\times \) 10\(^{4}\) and 1.22 \(\times \) 10\(^{4}\) tC, respectively.     

Simulated bone remodeling around tilted dental implants in the anterior maxilla

Dental implants have to be placed with the long axis in different angulations due to the change in bone morphology. The objective of this study was to investigate the different bone remodeling response induced by the tilted dental implants and to assess whether it could lead to bone loss and implant failure. In this study, bone remodeling due to palato-labially inclined dental implants placed in the anterior maxillary incisor region was simulated. CT-based finite element models of a maxillary bone with dental implants were created herein. Five dental implants were placed at \(+10^{\circ }\), \(+5^{\circ }\), \(0^{\circ }\), \(-5^{\circ }\) and \(-10^{\circ }\), respectively. The remodeling progression was recorded and compared. Model \(-10^{\circ }\) (palatal side) shows the highest bone density values, but the inclined implant at \(+10^{\circ }\) (labial side) leads to significant bone loss. From a biomechanical perspective, it is speculated that a palatally inclined implant is more likely to enhance the bone density in the maxillary anterior region, but labial inclination of implant could jeopardize its stability.     

Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur

A micro-finite element-based method to estimate the bone loading history based on bone architecture was recently presented in the literature. However, a thorough investigation of the parameter sensitivity and plausibility of this method to predict joint loads is still missing. The goals of this study were (1) to analyse the parameter sensitivity of the joint load predictions at one proximal femur and (2) to assess the plausibility of the results by comparing load predictions of ten proximal femora to in vivo hip joint forces measured with instrumented prostheses (available from www.orthoload.com). Joint loads were predicted by optimally scaling the magnitude of four unit loads (inclined \(-20^{\circ }\) to \(100^{\circ }\) with respect to the vertical axis) applied to micro-finite element models created from high-resolution computed tomography scans (\(30.3~\upmu \)m voxel size). Parameter sensitivity analysis was performed by varying a total of nine parameters and showed that predictions of the peak load directions (range 10\(^{\circ }\)–\(30^{\circ }\)) are more robust than the predicted peak load magnitudes (range 2344.8–4689.5 N). Comparing the results of all ten femora with the in vivo loading data of ten subjects showed that peak loads are plausible both in terms of the load direction (in vivo: \(18.2\pm 2.0^{\circ }\), predicted: \(20.0^{\circ }\)) and magnitude (in vivo: \(2707.6\pm 443.3~\hbox {N}\), predicted: \(3372.2\pm 597.9~\hbox {N}\)). Overall, this study suggests that micro-finite element-based joint load predictions are both plausible and robust in terms of the predicted peak load direction, but predicted load magnitudes should be interpreted with caution.     

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of \(^{15}\)N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of \(^{15}\)N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2–4\(\times\) due to relatively low “compressibility” of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of \(t_1\)/\(t_{relax}\) is sampled in a way similar to NUS of \(t_1\)/\(t_2\) domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional \(^{15}\)N HSQC spectra.     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

Examining the reaction between antioxidant compounds and 2,2-diphenyl-1-picrylhydrazyl (DPPH) through a computational investigation

In this work, we present a computational investigation on the reactions between two well-known antioxidants (quercetin and morin) and 2,2-diphenyl-1-picrylhydrazyl (DPPH). A density functional theory (DFT) approach with the B3LYP functional and the 6-31G(d,p) basis set was used for the simulations. The structural and energetic parameters (Gibbs free-energy, ΔG, and Gibbs free-energy of activation, ΔG⁺⁺) were determined to provide information on the antioxidant activity as well as to evaluate the contributions of each hydroxyl group to the referred property. According to the results obtained, quercetin presented three hydroxyls as being thermodynamically spontaneous in the reaction with DPPH (4\(^{\prime }\)-ArOH, 3\(^{\prime }\)-ArOH, and 3-ArOH, with ΔG = -4.93 kcal/mol, -2.89 kcal/mol, and -1.87 kcal/mol, respectively) against only two in the case of morin (2\(^{\prime }\)-ArOH and 3-ArOH, with ΔG = -7.56 kcal/mol and -4.57 kcal/mol, respectively). Hence, quercetin was found to be a more efficient antioxidant, which is in agreement with different experimental and computational investigations of bond dissociation enthalpies (BDEs). However, the order of contribution of the OH groups of each compound to the antioxidant potential present some differences when compared to what was seen in the previous investigations, especially for morin. These findings are in contrast to what was observed in studies based on the determinations of BDEs. Therefore, experimental investigations on the hydrogen-atom transfer mechanism (HAT) for both compounds are encouraged in order to clarify these observations.     

Linking chloroplast relocation to different responses of photosynthesis to blue and red radiation in low and high light-acclimated leaves of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> (L.)

Low light (LL) and high light (HL)-acclimated plants of A. thaliana were exposed to blue (BB) or red (RR) light or to a mixture of blue and red light (BR) of incrementally increasing intensities. The light response of photosystem II was measured by pulse amplitude-modulated chlorophyll fluorescence and that of photosystem I by near infrared difference spectroscopy. The LL but not HL leaves exhibited blue light-specific responses which were assigned to relocation of chloroplasts from the dark to the light-avoidance arrangement. Blue light (BB and BR) decreased the minimum fluorescence (\(F_{0}^{\prime }\)) more than RR light. This extra reduction of the \(F_{0}^{\prime }\) was stronger than theoretically predicted for \(F_{0}^{\prime }\) quenching by energy dissipation but actual measurement and theory agreed in RR treatments. The extra \(F_{0}^{\prime }\) reduction was assigned to decreased light absorption of chloroplasts in the avoidance position. A maximum reduction of 30% was calculated. Increasing intensities of blue light affected the fluorescence parameters NPQ and q_P to a lesser degree than red light. After correcting for the optical effects of chloroplast relocation, the NPQ responded similarly to blue and red light. The same correction method diminished the color-specific variations in q_P but did not abolish it; thus strongly indicating the presence of another blue light effect which also moderates excitation pressure in PSII but cannot be ascribed to absorption variations. Only after RR exposure, a post-illumination overshoot of \(F_{0}^{\prime }\) and fast oxidation of PSI electron acceptors occurred, thus, suggesting an electron flow from stromal reductants to the plastoquinone pool.     

Towards the prediction of flow-induced shear stress distributions experienced by breast cancer cells in the lymphatics

Tumour metastasis in the lymphatics is a crucial step in the progression of breast cancer. The dynamics by which breast cancer cells (BCCs) travel in the lymphatics remains poorly understood. The goal of this work is to develop a model capable of predicting the shear stresses metastasising BCCs experience using numerical and experimental techniques. This paper models the fluidic transport of large particles (\(\eta =d_{\mathrm{p}}/W=0.1-0.4\) where \(d_{\mathrm{p}}\) is the particle diameter and W is the channel width) subjected to lymphatic flow conditions (\({ Re}=0.04\)), in a \(100\times 100\,\upmu \hbox {m}\) microchannel. The feasibility of using the dynamic fluid body interaction (DFBI) method to predict particle motion was assessed, and particle tracking experiments were performed. The experiments found that particle translational velocity decreased from the undisturbed fluid velocity with increasing particle size (5–14% velocity lag for \(\eta =0.1-0.3\)). DFBI simulations were found to better predict particle behaviour than theoretical predictions; however, mesh restrictions in the near-wall region (\(0.2\,\mathrm{W}>y>0.8\,\mathrm{W}\)) result in computationally expensive models. The simulations were in good agreement with the experiments (\(<12\%\) difference) across the channel (\(0.2\,\mathrm{W}\le y\le 0.8\,\mathrm{W}\)), with differences up to 25% in the near-wall region. Particles experience a range of shear stresses (0.002–0.12 Pa) and spatial shear gradients (\(0.004-0.137\,\hbox {Pa}/\upmu \hbox {m}\)) depending on their size and radial position. The predicted shear gradients are far in excess of values associated with BCC apoptosis (\(0.004-0.023\,\hbox {Pa}/\upmu \hbox {m}\)). Increasing our understanding of the shear stress magnitudes and gradients experienced by BCCs could be leveraged to elucidate whether a particular BCC size or location exists that encourages metastasis within the lymphatics.     

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.     

Optimal Therapy Scheduling Based on a Pair of Collaterally Sensitive Drugs

Despite major strides in the treatment of cancer, the development of drug resistance remains a major hurdle. One strategy which has been proposed to address this is the sequential application of drug therapies where resistance to one drug induces sensitivity to another drug, a concept called collateral sensitivity. The optimal timing of drug switching in these situations, however, remains unknown. To study this, we developed a dynamical model of sequential therapy on heterogeneous tumors comprised of resistant and sensitive cells. A pair of drugs (DrugA, DrugB) are utilized and are periodically switched during therapy. Assuming resistant cells to one drug are collaterally sensitive to the opposing drug, we classified cancer cells into two groups, \(A_\mathrm{R}\) and \(B_\mathrm{R}\), each of which is a subpopulation of cells resistant to the indicated drug and concurrently sensitive to the other, and we subsequently explored the resulting population dynamics. Specifically, based on a system of ordinary differential equations for \(A_\mathrm{R}\) and \(B_\mathrm{R}\), we determined that the optimal treatment strategy consists of two stages: an initial stage in which a chosen effective drug is utilized until a specific time point, T, and a second stage in which drugs are switched repeatedly, during which each drug is used for a relative duration (i.e., \(f \Delta t\)-long for DrugA and \((1-f) \Delta t\)-long for DrugB with \(0 \le f \le 1\) and \(\Delta t \ge 0\)). We prove that the optimal duration of the initial stage, in which the first drug is administered, T, is shorter than the period in which it remains effective in decreasing the total population, contrary to current clinical intuition. We further analyzed the relationship between population makeup, \(\mathcal {A/B} = A_\mathrm{R}/B_\mathrm{R}\), and the effect of each drug. We determine a critical ratio, which we term \(\mathcal {(A/B)}^{*}\), at which the two drugs are equally effective. As the first stage of the optimal strategy is applied, \(\mathcal {A/B}\) changes monotonically to \(\mathcal {(A/B)}^{*}\) and then, during the second stage, remains at \(\mathcal {(A/B)}^{*}\) thereafter. Beyond our analytic results, we explored an individual-based stochastic model and presented the distribution of extinction times for the classes of solutions found. Taken together, our results suggest opportunities to improve therapy scheduling in clinical oncology.     

Failure properties of vena cava tissue due to deep penetration during filter insertion

In this work, we use an in-vitro mechanical test to explore the resistance of biaxially stretched vena cava tissue against deep perforation and a methodology which integrates experimental and numerical modeling to identify constitutive fracture properties of the vena cava. Six sheep vena cava were harvested just after killing, and cyclic uniaxial tension tests in longitudinal and circumferential directions and biaxial deep penetration tests were performed. After that, we use a nonlinear finite element model to simulate in vitro penetration of the cava tissue in order to fit the fracture properties under penetration of the vena cava by defining a cohesive fracture zone. An iterative process was developed in order to fit the fracture properties of the vena cava using the previously obtained experimental results. The proposed solutions were obtained with fracture energy of 0.22 or 0.33 N/mm. In comparison with the experimental data, the simulation using \(\delta _{0}=0.01\,\hbox {mm}\), \(\delta _{r}=0.35\,\hbox {mm}\), and \(K=220\, \hbox {N}/\hbox {mm}^{3}\) parameters (\(F_{\hbox {max}}=0.92\)) is in good agreement with results from penetration experiments of cava tissue. It is noticeable that the parameter estimation process of the fracture behavior is more accurate than the estimation process of the elastic behavior for the toe region of the curve.     

Cellular mechanosensitivity to substrate stiffness decreases with increasing dissimilarity to cell stiffness

Computational modelling has received increasing attention to investigate multi-scale coupled problems in micro-heterogeneous biological structures such as cells. In the current study, we investigated for a single cell the effects of (1) different cell-substrate attachment (2) and different substrate modulus \(\textit{E}_\mathrm{s}\) on intracellular deformations. A fibroblast was geometrically reconstructed from confocal micrographs. Finite element models of the cell on a planar substrate were developed. Intracellular deformations due to substrate stretch of \(\lambda =1.1\), were assessed for: (1) cell-substrate attachment implemented as full basal contact (FC) and 124 focal adhesions (FA), respectively, and \(\textit{E}_\mathrm{s}\,=\,\)140 KPa and (2) \(\textit{E}_\mathrm{s}\,=\,10\), 140, 1000, and 10,000 KPa, respectively, and FA attachment. The largest strains in cytosol, nucleus and cell membrane were higher for FC (1.35\(\text {e}^{-2}\), 0.235\(\text {e}^{-2}\) and 0.6\(\text {e}^{-2}\)) than for FA attachment (0.0952\(\text {e}^{-2}\), 0.0472\(\text {e}^{-2}\) and 0.05\(\text {e}^{-2}\)). For increasing \(\textit{E}_\mathrm{s}\), the largest maximum principal strain was 4.4\(\text {e}^{-4}\), 5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\) and 5.3\(\text {e}^{-4}\) in the membrane, 9.5\(\text {e}^{-4}\), 1.1\(\text {e}^{-4}\), 1.2\(\text {e}^{-3}\) and 1.2\(\text {e}^{-3}\) in the cytosol, and 4.5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\), 5.7\(\text {e}^{-4}\) and 5.7\(\text {e}^{-4}\) in the nucleus. The results show (1) the importance of representing FA in cell models and (2) higher cellular mechanical sensitivity for substrate stiffness changes in the range of cell stiffness. The latter indicates that matching substrate stiffness to cell stiffness, and moderate variation of the former is very effective for controlled variation of cell deformation. The developed methodology is useful for parametric studies on cellular mechanics to obtain quantitative data of subcellular strains and stresses that cannot easily be measured experimentally.     

Constitutive modeling of an electrospun tubular scaffold used for vascular tissue engineering

In this study, we sought to model the mechanical behavior of an electrospun tubular scaffold previously reported for vascular tissue engineering with hyperelastic constitutive equations. Specifically, the scaffolds were made by wrapping electrospun polycaprolactone membranes that contain aligned fibers around a mandrel in such a way that they have microstructure similar to the native arterial media. The biaxial stress-stretch data of the scaffolds made of moderately or highly aligned fibers with three different off-axis fiber angles \(\alpha \) (30\(^\circ \), 45\(^\circ \), and 60\(^\circ \)) were fit by a phenomenological Fung model and a series of structurally motivated models considering fiber directions and fiber angle distributions. In particular, two forms of fiber strain energy in the structurally motivated model for a linear and a nonlinear fiber stress–strain relation, respectively, were tested. An isotropic neo-Hookean strain energy function was also added to the structurally motivated models to examine its contribution. The two forms of fiber strain energy did not result in significantly different goodness of fit for most groups of the scaffolds. The absence of the neo-Hookean term in the structurally motivated model led to obvious nonlinear stress-stretch fits at a greater axial stretch, especially when fitting data from the scaffolds with a small \(\alpha \). Of the models considered, the Fung model had the overall best fitting results; its applications are limited because of its phenomenological nature. Although a structurally motivated model using the nonlinear fiber stress–strain relation with the neo-Hookean term provided fits comparably as good as the Fung model, the values of its model parameters exhibited large within-group variations. Prescribing the dispersion of fiber orientation in the structurally motivated model, however, reduced the variations without compromising the fits and was thus considered to be the best structurally motivated model for the scaffolds. It appeared that the structurally motivated models could be further improved for fitting the mechanical behavior of the electrospun scaffold; fiber interactions are suggested to be considered in future models.     

Cell Volume Regulation in the Proximal Tubule of Rat Kidney

We developed a dynamic model of a rat proximal convoluted tubule cell in order to investigate cell volume regulation mechanisms in this nephron segment. We examined whether regulatory volume decrease (RVD), which follows exposure to a hyposmotic peritubular solution, can be achieved solely via stimulation of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. We also determined whether regulatory volume increase (RVI), which follows exposure to a hyperosmotic peritubular solution under certain conditions, may be accomplished by activating basolateral \(\hbox {Na}^+\)/H\(^+\) exchangers. Model predictions were in good agreement with experimental observations in mouse proximal tubule cells assuming that a 10% increase in cell volume induces a fourfold increase in the expression of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. Our results also suggest that in response to a hyposmotic challenge and subsequent cell swelling, \(\hbox {Na}^+\)–\(\hbox {HCO}^-_3\) cotransporters are more efficient than basolateral K\(^+\) and \(\hbox {Cl}^-\) channels at lowering intracellular osmolality and reducing cell volume. Moreover, both RVD and RVI are predicted to stabilize net transcellular \(\hbox {Na}^+\) reabsorption, that is, to limit the net \(\hbox {Na}^+\) flux decrease during a hyposmotic challenge or the net \(\hbox {Na}^+\) flux increase during a hyperosmotic challenge.     

Investigating the effects of membrane deformability on artificial capsule adhesion to the functionalized surface

Understanding, manipulating and controlling cellular adhesion processes can be critical in developing biomedical technologies. Adhesive mechanisms can be used to the target, pattern and separate cells such as leukocytes from whole blood for biomedical applications. The deformability response of the cell directly affects the rolling and adhesion behavior under viscous linear shear flow conditions. To that end, the primary objective of the present study was to investigate numerically the influence of capsule membrane’s nonlinear material behavior (i.e. elastic-plastic to strain hardening) on the rolling and adhesion behavior of representative artificial capsules. Specifically, spherical capsules with radius of \(3.75\, \upmu \hbox {m}\) were represented using an elastic membrane governed by a Mooney–Rivlin strain energy functions. The surfaces of the capsules were coated with P-selectin glycoprotein-ligand-1 to initiate binding interaction with P-selectin-coated planar surface with density of \(150\,\upmu \hbox {m}^{-2}\) under linear shear flow varying from 100 to \(400\,\hbox {s}^{-1}\). The numerical model is based on the Immersed Boundary Method for rolling of deformable capsule in shear flow coupled with Monte Carlo simulation for receptor/ligand interaction modeled using Bell model. The results reveal that the mechanical properties of the capsule play an important role in the rolling behavior and the binding kinetics between the capsule contact surface and the substrate. The rolling behavior of the strain hardening capsules is relatively smoother and slower compared to the elastic-plastic capsules. The strain hardening capsules exhibits higher contact area at any given shear rate compared to elastic-plastic capsules. The increase in contact area leads to decrease in rolling velocity. The capsule contact surface is not in complete contact with the substrate because of thin lubrication film that is trapped between the capsule and substrate. This creates a concave shape on the bottom surface of the capsule that is referred to as a dimple. In addition, the present study demonstrates that the average total bond force from the capsules lifetime increases by 37 % for the strain hardening capsules compared to elastic-plastic capsules at shear rate of \(400\,\hbox {s}^{-1}\). Finally, the model demonstrates the effect of finite membrane deformation on the coupling between hydrodynamic and receptor/ligand interaction.