Next Article in Journal
Biotransformation of d-Xylose-Rich Rice Husk Hydrolysate by a Rice Paddy Soil Bacterium, Priestia sp. Strain JY310, to Low Molecular Weight Poly(3-hydroxybutyrate)
Next Article in Special Issue
Transportan 10 Induces Perturbation and Pores Formation in Giant Plasma Membrane Vesicles Derived from Cancer Liver Cells
Previous Article in Journal
Mapping DNA Conformations Using Single-Molecule Conductance Measurements
Previous Article in Special Issue
Phosphatidic Acid Accumulates at Areas of Curvature in Tubulated Lipid Bilayers and Liposomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Perspective

Energy Dissipation in the Human Red Cell Membrane

by
Thomas M. Fischer
Laboratory for Red Cell Rheology, Krummer Weg 20, 52134 Herzogenrath, Germany
Biomolecules 2023, 13(1), 130; https://doi.org/10.3390/biom13010130
Submission received: 9 December 2022 / Revised: 4 January 2023 / Accepted: 6 January 2023 / Published: 9 January 2023
(This article belongs to the Special Issue Recent Developments in Biophysical Studies of Cell Membranes)

Abstract

:
The membrane of the human red cell consists of a lipid bilayer and a so-called membrane skeleton attached on the cytoplasmic side of the bilayer. Upon the deformation of red cells, energy is dissipated in their cytoplasm and their membrane. As to the membrane, three contributions can be distinguished: (i) A two-dimensional shear deformation with the membrane viscosity as the frictional parameter; (ii) A motion of the membrane skeleton relative to the bilayer; (iii) A relative motion of the two monolayers of the bilayer. The frictional parameter in contributions (ii) and (iii) is a frictional coefficient specific for the respective contribution. This perspective describes the history up to recent advances in the knowledge of these contributions. It reviews the mechanisms of energy dissipation on a molecular scale and suggests new ones, particularly for the first contribution. It proposes a parametric fitting expected to shed light on the discrepant values found for the membrane viscosity by different experimental approaches. It proposes strategies that could allow the determination of the frictional coefficients pertaining to the second and the third contribution. It highlights the consequences characteristic times have on the state of the red cell membrane in circulation as well as on the adaptation of computer models to the red cell history in an in vitro experiment.

1. Introduction

Compared to other cells in the human body, the structure of red cells is relatively simple. A closed membrane is filled with a concentrated solution of hemoglobin—from a mechanical point of view a Newtonian liquid. The shape is discoidal with average values of diameter and thickness of 8 µm and 2 µm, respectively. The membrane consists of a lipid bilayer and a so-called membrane skeleton attached to the cytoplasmic side of the bilayer. The lipid bilayer is referred to as a bilayer in the following. From a mechanical point of view, it can be considered a two-dimensional (2D) liquid. It provides bending stiffness to the membrane. Under the experimental situations considered here, the bilayer is essentially incompressible. So far, it also conserves the surface area of the membrane. Embedded in the bilayer are so-called intrinsic proteins spanning the bilayer. The membrane skeleton is referred to as the skeleton in the following. Its main constituents are the proteins spectrin and actin forming a 2D (loosely knitted) lattice. The skeleton provides shear elasticity to the membrane. It is connected to the intrinsic proteins Band 3 and glycophorin.
The coefficients describing the mechanical behavior of the membrane have been measured by comparison between experiments and adapted theoretical descriptions. Among these coefficients is the 2D viscosity of the red cell membrane. By comparison with the orders of magnitude lower value of the bilayer, the respective energy dissipation was considered to occur in the skeleton. Although the determination of the membrane viscosity became standard in the characterization of membrane mechanical properties, an attempt to explain the responsible energy dissipation on a molecular scale was never made.
With the advances in computational technology, membrane models of increasing sophistication were developed. (These models typically treat the membrane as a 2D sheet with no thickness. This has usually been justified by comparing red cell dimensions with the bilayer thickness of ~4 nm. By lowering a force-sensitive probe onto the membrane, Heinrich et al. [1] observed the first elastic recoil at a skeleton thickness of ~90 nm. Nevertheless, the 2D idealization is commonly adopted in modeling.) Among these, four recent studies [2,3,4,5] presented parametric fits of the membrane viscosity to existing experimental results. However, the results of these studies suggest that additional sources of energy dissipation have to be considered to describe the dynamic behavior of the red cell.
The paper is organized as follows. Section 2 reviews the collection of the knowledge with respect to the energy dissipation in the membrane by (i) Dynamic membrane shear with the membrane viscosity as the frictional parameter; (ii) Skeleton slip; i.e., the relative motion of bilayer and skeleton; (iii) Monolayer slip; i.e., the relative motion of the two monolayers of the bilayer. The frictional parameter in contributions (ii) and (iii) is a frictional coefficient specific to the respective contribution. Section 3 presents possible mechanisms on a molecular scale and illuminates some of them by reference to particular experiments. Section 4 proposes a parametric fitting expected to shed light on the discrepant values found for the membrane viscosity by different experimental approaches. Furthermore, strategies are proposed that could allow the determination of the friction coefficients pertaining to the second and the third contribution. Section 5 highlights the consequences the various characteristic times have on the state of the red cell membrane in the circulation as well as on the adaptation of computer models to the red cell history in an in vitro experiment. Section 6 ties the presented information together and ventures an outlook.

2. History of the Consideration of the Energy Dissipation in the Red Cell Membrane

Section 2.1 deals with the various methods used to determine the membrane viscosity. A matter of concern were and, to some extent, still are the different values obtained by different methods thus shedding doubt on the membrane viscosity being the sole contribution to the energy dissipation in the red cell membrane. The findings of Rezghi and Zhang [5] added to this suspicion. Section 2.2 and Section 2.3 describe the mechanisms of the two additional contributions to energy dissipation. No parametric fits were undertaken for these two cases. For the second contribution, theoretical estimates of the involved frictional coefficient differ by three orders of magnitude thus demonstrating the necessity of a direct determination by parametric fits to appropriate experiments.

2.1. Dynamic Membrane Shear—The First Contribution to the Energy Dissipation

The shear stiffness was evaluated by aspirating a membrane tongue into small bore glass pipettes and the shear modulus ( μ ) was determined as 7 µN/m [6]. In those days (1970s–1990s), the skeleton was considered to be incompressible—equivalent to local conservation of its surface area. The first account of the viscosity ( η mem ¯ , where the overline indicates that a single value was determined) used the time constant of relaxation of a membrane tongue after an aspirated cell was expelled from the pipette [7]. η mem ¯ = 10 6 Pa s m was given as an order of magnitude estimate. Since this value was three orders of magnitude above what was estimated to be the viscosity of the bilayer, η mem ¯ was considered to be the viscosity of the skeleton.
Chien et al. [8] observed the time courses during the aspiration of a membrane tongue into small bore pipettes and during its release. Several characteristic times were recorded, τ 1 and τ 2 during the entry of the tongue and τ 3 upon release of the aspiration pressure. The shear modulus was measured as μ = 4.2 µN/m and the membrane viscosity was calculated according to η mem ¯ = μ τ . From τ 1 , η mem ¯ was found in the range of (0.6–4) · 10 7 Pa s m varying inversely with the initial rate of entry which in turn depended on the final tongue length. Based on the τ 2 results, the membrane viscosity was 2 · 10 5 Pa s m. Using τ 3 , the membrane viscosity was estimated as 5.44 · 10 7 Pa s m after 20 s of deformation. Upon reduction of this time interval down to 2 s, τ 3 decreased to approach the low values of τ 1 . As an explanation for the complicated behavior, the authors suggested a shear thinning lasting over a certain period of time. Tözeren et al. [9] interpreted the two time scales upon entry as a creep phenomenon and presented a description in closed form. Fischer [10] ascribed τ 2 to skeleton slip introduced in Section 2.2.
Hochmuth et al. [11] observed the relaxation of a red cell after being released from a static extension at two diametrically opposed locations on the membrane. From the characteristic time of the relaxation ( τ = 0.1 0.13 s) and μ = 6 µN/m, η mem ¯ = 6–8 · 10 7 Pa s m was determined. Without going into detail, it is mentioned that this concept was later used in numerous studies by stretching the red cells using pipettes, optical tweezers, electrical fields or shear flows.
Tran-Son-Tay et al. [12] considered red cells suspended in a solution of viscosity ( η 0 ) of 35 mPa s and subjected to simple shear flow with shear rate ( γ ˙ ). In a steady state, the red cells display an elongated shape and a continuous flow of the membrane around that shape [13]. This particular motion of the membrane has been called the tank-tread (TT) motion [14]. The kinematics of the TT motion are sketched in Figure 1. Tran-Son-Tay et al. [12] collected the geometric data of the red cells and their TT frequency ( T T F ). Equating the power supplied to the red cell from the shear flow with the energy dissipated in the red cell, they deduced η mem ¯ . The value depended moderately on γ ˙ . At γ ˙ = 29 /s they found η mem ¯ = 1.03 · 10 7 Pa s m and at γ ˙ = 171 /s the value was 0.86 · 10 7 Pa s m, an order of magnitude below the value given by Hochmuth et al. [11].
Later this group refined the treatment and determined η mem depending on the location on the red cell membrane [16]. The range of these values even overlapped with that found for pure lipid bilayers [17].
With the appearance of discretized computer models, previous determinations of η mem ¯ were reevaluated by fitting these models to experimentally observed T T F s [18]. As in Tran-Son-Tay et al. [12], a decrease of η mem ¯ with the increase in γ ˙ was found. Tsubota [3] fitted data collected at η 0 = 54 mPa s and γ ˙ from 8/s to 260/s and found η mem ¯ to decrease from 4 to 1.5 · 10 7 Pa s m. Matteoli et al. [2] used data at η 0 = 13 mPa s from γ ˙ = 30 /s to 260/s and at η 0 = 54 mPa s from γ ˙ = 8 /s to 130/s and they found a common decrease of η mem ¯ from 24 to 2.8 · 10 8 Pa s m. It is surprising that the values of η mem ¯ at the same values of γ ˙ and η 0 are quite different in these two studies, although the axis ratios of the respective stress-free shapes were similar (Table 1).
The determination of η mem ¯ via the recovery time was revisited by Guglietta et al. [4]. Comparison of their parametric study with the relaxation data of Mills et al. [19] results in η mem ¯ 4.2 · 10 7 Pa s m. The difference to the respective values determined via TTing is smaller (line 4 and 5 in Table 2) than before (line 2 and 3 in Table 2) but is still existing.
Rezghi and Zhang [5] criticized the studies by Matteoli et al. [2] and Tsubota [3] in that the same η mem value was used for all membrane elements instead of determining η mem as a function of the local 2D membrane shear rate as one would expect for a non-Newtonian viscosity. Actually, this consideration had been adopted by Sutera et al. [16], where the membrane viscosity was assumed as a function of the second invariant of the surface strain rate. Rezghi and Zhang [5] fitted to T T F s collected at four values of η 0 (13, 29, 54, and 109 mPa s). In contrast to the authors in Table 1, they used the biconcave disc as the stress-free shape. A multitude of values for η mem ¯ was reported which for brevity are not repeated here. In agreement with the authors in Table 1, the trend in these data indicates a decrease of η mem ¯ with γ ˙ . On the other hand, comparison with experimental data on the angle of inclination during TTing [20] yielded a constant value for η mem ¯ . Furthermore, fitting the elongation of TTing red cells [21], even resulted in an increase of η mem ¯ with γ ˙ . Among other explanations, the authors conjectured that the fact that different blood samples had been used in the measurement of these three parameters might be responsible for the different behavior. However, such drastic discrepancies call for a widened view and to include additional sources of energy dissipation in the theoretical models.

2.2. Skeleton Slip—The Second Contribution to the Energy Dissipation

Let us go back to the early days when the skeleton was considered to be incompressible. Stokke et al. [22], considering the skeleton as an ionic gel, widened the view by claiming that the skeleton has finite compressibility. This means that besides the shear modulus a biaxial modulus describing the resistance of the skeleton to a change in the area has to be considered. If the red cell shape differs from the reference configuration of the skeleton, the minimum in both its shear and biaxial strain energy requires its rearrangement towards a non-uniform density on the membrane [23]. Provided the biaxial modulus is not much larger than the shear modulus, a non-uniform skeleton density becomes apparent in deformations with strong non-uniform local shear strains. An example is the aspiration of a membrane tongue into a small bore glass pipette. Fischer [10] predicted that the skeleton density would increase at the pipette tip and decrease at the tip of the aspirated membrane tongue. This redistribution was demonstrated experimentally by Discher et al. [24].
A dynamic redistribution gives rise to a second source of energy dissipation within the membrane. In mechanical terms, a redistribution is equivalent to a slip between the skeleton and the bilayer. Accordingly, its contribution to energy dissipation is called skeleton slip in the following. The mechanical effect of the skeleton slip can be represented by:
σ sb = c f · Δ v sb
where σ sb is shear stress opposing the slip between skeleton and bilayer, c f is the respective frictional coefficient, and Δ v sb is the difference in velocity between the skeleton and the bilayer.
Since the skeleton is connected to the bilayer via intrinsic proteins, skeleton slip involves a relative motion of these proteins and the bilayer. A first approach to quantify the lateral drag of these anchor proteins through the bilayer was made by Fischer [10]. The drag was estimated using a theoretical expression of the force necessary to move a circular cylinder embedded in a thin liquid sheet of finite size [25].
Peng et al. [26] estimated c f using the lateral diffusivity of the intrinsic anchor protein Band 3. They determined c f = ( ρ / ρ 0 ) 144 · 10 6 Ns/m³ where ρ is the actual and ρ 0 the undisturbed density of this protein. Accounting for both membrane shear and skeleton slip, Peng and Zhu [27] modeled TTing red cells. The viscosities of cytoplasm and suspending phase were both set at 13 mPa s, and the shear rate γ ˙ = 200/s and η mem ¯ = 5 · 10 8 Pa s m [12] were used. The results showed that it took about one minute or 392 TT cycles to reach the steady state in the distribution of skeleton density. As to the cell elongation, the initial fast rise was followed by a small decrease in parallel to the attainment of the steady state in skeleton distribution. The T T F remained basically constant.
Accounting for dynamic membrane shear alone, the steady state in red cell elongation is attained within a few TT cycles [26]. Accordingly, the initial fast rise in red cell elongation in [26] can be considered to be controlled by dynamic membrane shear. The slight decrease in red cell elongation can be attributed to a small decrease in energy dissipation after the redistribution of the skeleton has reached a steady state.
In their determination of c f , Peng et al. [26] used the larger of the two values determined for the diffusion coefficient by Kodippili et al. [28]. As discussed by the last authors, the attribution of the two values to the different fractions of Band 3 is presently unclear. Had Peng et al. [26] used the lower value, the redistribution would have been calculated to occur an order of magnitude more slowly. On the other hand, the use of the experimental diffusion coefficients has been criticized by Turlier et al. [29]. Based on an equation of Saffman and Delbrück [30], these authors calculated a value for the diffusion coefficient of free Band 3 two orders of magnitude larger (Their estimate was based on the diameter of the naked Band 3 molecule. After correction to the average diameter of the Band 3 complexes the ratio decreased somewhat but the two orders of magnitude remained.) than the value used by Peng et al. [26]. This discrepancy calls for an experimental determination of c f .

2.3. Monolayer Slip—The Third Contribution to the Energy Dissipation

In the early days of red cell mechanics studies, the bending stiffness of the red cell membrane was considered to result from the low compressibility of the two monolayers of the bilayer and their fixed distance [31]. Assuming that the biaxial modulus (K) of the bilayer is distributed equally between the two monolayers, the respective bending stiffness ( B bc ) reads [31]:
B bc = K h 2 4
where h is the distance between the neutral planes of the two monolayers. The motivation for the index “bc” will become apparent below.
Helfrich [32] on the other hand, by analogy to liquid crystals, claimed an intrinsic bending stiffness of each monolayer. The sum of the intrinsic bending stiffnesses of both monolayers would constitute the intrinsic bending stiffness of the bilayer. In the beginning, this stiffness was deduced from the observation of membrane flickering of phospholipid vesicles and tubes, e.g., [33,34]. Later Waugh and coworkers used pulling of so-called tethers of the bilayer to this end [35,36].
To distinguish these two mechanisms, the intrinsic bending stiffness is referred to as single-layer bending (slb) hereafter, and the mechanism described by Equation (2) is referred to as bilayer couple bending (bcb). In the literature, slb is often referred to as local bending because its elastic energy depends on the local curvature provided the lateral distribution of the different lipid species is uniform on the membrane. Bilayer couple bending is often referred to as global bending because in mechanical equilibrium its elastic energy depends on the average of all curvatures on the membrane. Here comes the third contribution to the energy dissipation into play.
If red cell deformations induce fast changes of local curvature, inner and outer monolayers are stressed in the opposite sense—in the above-introduced nomenclature: local bcb. At locations with positive curvatures, the outer layer is extended and the inner one compressed. As to the sign convention, a spherical membrane is defined to have a positive curvature. In a closed membrane, a positive change in curvature at some location requires a negative change at another location. This leads to a difference in lateral pressure in both monolayers which is relieved by a flow of their lipids in opposite directions—local bcb relaxes to global bcb. In the above-mentioned experiments on flickering [33,34], the changes in curvature are small. Obviously, bcb was essentially global during the experiments and the effects of slb were detected only.
The flow in opposite directions of the lipids in the two monolayers requires their flow relative to the intrinsic proteins which can be considered stationary if we conceptually assume that bending and relaxation occur consecutively. The respective frictional contribution was estimated [37] using calculations of the flow through regular arrays of parallel solid cylinders [38]. It was found that the influence of local bcb during tether extension [39] and during TTing γ rheo = 42 /s, in ref [40] cannot be neglected.
Besides the differential flow around the intrinsic proteins, the friction involved in the opposite motion of the interdigitating ends of the hydrocarbon chains of the inner and outer monolayers contributes as well. The respective frictional coefficient was estimated using two types of vesicle [41]. One consisted of a pure phospholipid. The other one was a binary mixture. A single order of magnitude was given for the frictional coefficient ( 10 8 Ns / m 3 ).
Using a coarse-grained molecular dynamics simulation of a bilayer with embedded proteins, Khoshnood et al. [42] found an order of magnitude value for the frictional coefficient of 10 6 Ns / m 3 . The boundary condition was shear stress acting on both sides of the simulated bilayer which is of course different from that driving monolayer slip. Furthermore, as mentioned by the authors, coarse-grained molecular dynamics simulations typically underestimate the actual values.

3. Mechanisms on a Molecular Scale

Historically, the origin of η mem ¯ was attributed to the skeleton proper and that solely by exclusion. Although its value was determined in numerous studies for almost half a century, a mechanism for dynamic membrane shear on a molecular scale has not been suggested so far. This gap is filled in Section 3.1. Most mechanisms for the second and third contributions have been described in Section 2. For the second contribution, an additional mechanism is put forward which, however, depends on the boundary conditions of the respective experiment (Section 3.2).

3.1. Dynamic Membrane Shear

Here are several mechanisms that could be responsible for the membrane shear viscosity:
1a 
Stretching single spectrin monomers resulted in a higher tension when the stretching rate was increased from 0.08 to 0.8 µm/s [43]. This indicates a viscous resistance against stretch. To set these rates in perspective to experiments, red cell TTing is considered. Using the images shown in [40], the stretch of a single spectrin dimer oriented at a right angle to the flow direction is estimated based on the difference of half the perimeter minus the width of the TTing red cell. Multiplication with 4 · T T F gives 0.2, 0.7, 1.5, and 4 µm/s at γ ˙ = 42, 125, 238, and 575/s, respectively.
1b 
It is conceivable that the flow of the cytoplasmic hemoglobin solution in the voids between the convoluted spectrins to accommodate the shear deformation may also contribute to membrane shear. The extent of this contribution could be estimated by comparing determinations of η mem ¯ in hemoglobin-free red cells (ghosts) and in red cells. The available experimental data do not provide a clear-cut picture. Waugh [44] found a 36% lower value for η mem ¯ for ghosts compared to red cells. Nash and Meiselman [45], on the other hand, found a decrease of only 6%.
By lowering a force-sensitive probe onto the ghost membrane, Heinrich et al. [1] probed the elastic thickness compressibility of the red cell membrane. The first elastic recoil was observed at a skeleton thickness of ~0.090 µm. An appreciable value was observed at ~0.050 µm. This value is in keeping with the results of an electron microscopic study [46]. The hardcore thickness was found by Heinrich et al. [1] at ~0.024 µm. Based on these numbers, the contribution of the envisaged mechanism could be estimated by applying the concept of flow through porous media.
1c 
For the sake of completeness, the analogous mechanism with respect to the glycocalyx on the outside of the membrane is listed here.
1d 
Consider a red cell after release from a static elongation—the standard experiment to determine η mem ¯ . The contracting skeleton drags the Band 3 complexes through the bilayer as illustrated schematically in Figure 2. Figure 2A shows the smallest unit of the idealized skeleton in the relaxed state. The circles indicate the cross sections of the intrinsic portions of Band 3 complexes and Band 3 dimers drawn to scale as adopted from Burton and Bruce [47]. Drawn are 4 dimers. The average number is 3.5. The edge length of the triangle is 75 nm [48].
Figure 2B shows two hexagons each composed of six triangles as in Figure 2A. For simplicity, the Band 3 complexes and free Band 3 dimers are not shown. The left hexagon is extended 1.5-fold. The right one is relaxed. The red arrows in the left hexagon indicate the direction in which the Band 3 complexes are dragged during the relaxation towards the positions shown in the right hexagon. Considering the motion of a single Band 3 complex, one is tempted to call the situation another type of skeleton slip. However as shown in the left hexagon, there is no net relative motion between the skeleton and bilayer. On average, the displacements of all Band 3 complexes cancel. This is a distinctive difference to skeleton slip which can be appreciated by comparing Figure 2B,C.
Figure 2C sketches the situation during the aspiration of a membrane tongue into a small bore pipette. Here, the hexagon moves as a whole relative to the bilayer whereas in Figure 2B the hexagon just changes its shape and does not move. Besides the relative motion, the surface area changes as sketched in Figure 2C. This change is linked to the skeleton slip but is not quoted in the context of the distinctive difference.

3.2. Skeleton Slip

2a 
The friction due to the drag of the intrinsic proteins through the bilayer has been described in Section 2.2.
2b 
For a contribution of cytoplasmic and ambient fluid flow through the layer of the skeletal proteins and the glycocalyx, the geometries of the respective experiments have to be considered. In TTing the shear stress effectively acts on the skeleton at the cytoplasmic side and on the glycocalyx on the outside. The 2D membrane flow would be slightly deformed due to the redistribution of the skeleton. However, the effect on energy dissipation is expected to be minute. During relaxation after release from an elongation, skeleton slip is directed from the attachment points towards the cell body. Again, the effect is expected to be small. The situation might be different during the aspiration of a membrane tongue. On the cytoplasmic side, the skeleton is squeezed like a sponge due to the increase in its local density. The same applies on the outside with the condensing glycocalyx. The last contribution could be estimated based on the observation of the fluid flow in the small gap between the membrane and the inner pipette wall [49].

3.3. Monolayer Slip

The two contributing mechanisms have been already described in Section 2.3 and are therefore just itemized here.
3a 
The flow of lipids around intrinsic proteins.
3b 
The friction due to the sliding of the hydrocarbon chains past each other.
An effect of these mechanisms might be observable during TTing where the membrane flows over the rim of an elongated red cell. The curvature at the rim is higher than in the rest of the cell thus producing gradients in the surface density of both monolayers which lead to a steady state flow of lipids relative to the skeleton.

4. Suggested Steps to Widen the Knowledge of Membrane Mechanical Parameters

The following items are the ideas of an experimentalist. A realization will depend on the ratio of necessary effort and possible results. Steps 1 to 3 are targeted to evolve η mem ¯ from a descriptive parameter depending on the respective experiment towards a universal physical coefficient. Step 5 attempts to pinpoint the molecular mechanism of dynamic membrane shear. Steps 4 and 6 suggest the first steps to put numbers to the respective frictional coefficients.
  • The discordant behavior shown in Table 2 could be revisited by using the data in [18]. For each value of the T T F , the axis dimensions L and W have been determined and are available upon request. (The respective raw data have been sent to a number of groups but have never been used in conjunction.) Using these data, both μ and η mem ¯ could be fitted. The value of the bending stiffness (in slb) should be sufficient to prevent buckling during TTing. Since the strain hardening behavior of the shear elasticity of red cell membrane is not known, μ should depend on the 2D shear strain ( ϵ ) of the membrane. The deformation and T T F in steady state should allow the determination of μ ( ϵ ) and η mem ¯ . The subsequent stop of the shear flow should allow to determine a second value for η mem ¯ from the recovery of the elongation. The double determination in the same model using the same data has not been done before.
    If the model uses the stress–strain law of Skalak et al. [50], the above calculations represent a test of a proposal of Dimitrakopoulos [51] who claimed that this law together with the small deformation value of μ measured by Lenormand et al. [52] describes various experiments sufficiently. Generally, the determined dependency μ ( ϵ ) indicates how close the constitutive equation of the respective model is to the real strain hardening behavior of the red cell membrane. The interpretation of the results requires some caution. At low red cell elongations, the membrane shear strains are small. As a consequence, the error in the fitted values of μ is larger the more the reference configuration of the model differs from reality. However, with increasing elongations, this error decreases.
  • Most likely, the trend of decreasing η mem ¯ with increasing γ ˙ as observed previously will persist in step 1. As a refinement, the calculations could be repeated with a membrane viscosity η mem depending on the 2D shear rate in the membrane.
  • The variation in the raw data of L and W is rather large. After excluding extreme values, the range normalized by the mean value is about 0.3. To exclude the error due to this variation, step 2 could be performed for every single cell.
  • In aspirating cell membrane into small bore pipettes, Chien et al. [8] found a second, much slower characteristic time ( τ 2 ) upon entry. The corresponding additional increase in tongue length was suggested to be due to skeleton slip [10]. If correct, it would allow a parametric fit of c f . Since no raw data from these experiments are available, the following procedure is suggested. Use a model with an incompressible skeleton. Determine the time course of tongue length in a typical pipette experiment of [8]; e.g., pressure difference 18 mm H 2 O , pipette radius 0.3 µm, μ = 4.2 µN/m, and η mem ¯ = 6 · 10 8 Pa s m. Add compressibility of the skeleton to the model and allow for skeleton slip. Repeat the calculation with c f = ( ρ / ρ 0 ) 144 · 10 6 Ns/m³. If the tongue length increases further, adjust c f to obtain a value of τ 2 = 5 s as found in [8]. This would constitute an experimental determination of c f .
    If the tongue length decreases, the origin of the observed increase is not related to skeleton slip as suggested [10]. However, a close look at the recorded trace of the tongue length in Figure 7 of [8] shows a decreasing trend 20 s after the aspiration. This trend might be due to the modeled decrease. In order to be able to perform a parametric fit, a complete recording of the time course of the tongue length would be required.
    Another useful information in this respect might be the observation that the redistribution of the skeleton was complete after < 1 min after aspiration of a membrane tongue into a pipette [49]. The time course was not given in this study probably because the required exposure time did not allow to record it. With 3 τ = 1 min, at least an upper bound for c f could be obtained.
  • As noted in Section 3.1, the friction due to the drag of a single Band 3 complex or Band 3 dimer through the bilayer is the same in skeleton slip and dynamic membrane shear. Therefore, it should be possible to express η mem ¯ as a function of c f found in step 4. If this value corresponds to the values determined in step 1, mechanism 1d (Section 3.1) can be considered dominant. Otherwise, additional mechanisms have to be taken into account.
  • As mentioned in Section 3.3, bcb may be local at the tips of strongly elongated and TTing red cells. To deduce a value of c f in monolayer slip, add slb with 2 · 10 19 J and bcb with 6 · 10 19 J (Considering each monolayer consisting of the same isotropic material, the ratio between the bending stiffnesses in bcb and slb is three. With a large variation, Waugh et al. [35] found a ratio of 3.4 in phospholipid vesicles.) as the respective stiffness parameters to the model in step 1. Choose the spontaneous curvature such that ratio of cell thickness in the dimple and the rim equals 0.55 [53,54]. Allow for monolayer slip during TTing and start with c f = 0 . If the curvatures at low elongations and low T T F s compare well with observations but deviate at high elongations and high T T F s, increase the frictional coefficient until you get agreement for all shapes. This would constitute an experimental determination of c f .

5. Consequences for the Modeling of Red Cell Experiments

A slow redistribution of skeleton density as found in [27] would have consequences for the reference configuration of the shear elasticity required in modeling. The circulation time through the vasculature is about 1 min. During a typical lap, the red cells are deformed in different shapes leaving not enough time to attain the steady state distribution of the skeleton in each of these shapes. In the long run, the distribution of the skeleton is therefore expected to be the mean of these (conceptual) steady state distributions weighted by their relative life time. This envisioned kind of average represents a modification of an equivalent suggestion that has been made earlier for the reference configuration in general [18,54,55].
The suggested average of skeleton density pertains to the situation in vivo. For parametric fits to experiments, the situation may be different. During handling of the blood after withdrawal as well as during the execution of an experiment, the envisioned average of skeleton density prevailing in the circulation may change. Based on the value of c f determined in Section 4 step 4, the time course of skeleton redistribution can be determined. Comparison of this time course to the history of the red cells in vitro could be used to choose/guess the appropriate reference configuration on which the modeling has to be based. The same argument applies to the relaxation of local bcb to global bcb based on the value of c f determined in Section 4 step 6.

6. Conclusions

Three contributions to the energy dissipation in the membrane under dynamic red cell deformations are identified. The first (dynamic membrane shear) corresponds to the well-known 2D membrane viscosity. Despite a lot of experimental and theoretical work, the underlying mechanisms on a molecular scale are not known and its frictional parameter ( η mem ¯ ) remains descriptive; i.e., it depends on the respective experiment. For the second and third contributions (skeleton slip and monolayer slip), the underlying mechanisms on a molecular scale are known, however, the respective frictional coefficients are unknown. Steps to widen the knowledge are suggested. If successful, it is expected to determine the relative share of the three contributions to the energy dissipation in the various deformation modes human red cells are subjected to in vivo as well as in vitro. The values of the characteristic times associated with the second and third contributions will determine to what extent skeleton slip and monolayer slip actually contribute in a particular situation. Finally, it is expected that the improved modeling will facilitate the interpretation of experiments performed to provide information on the reference configuration of the red cell skeleton the last big unknown of mechanical membrane properties.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks Junfeng Zhang, Bharti School of Engineering and Computer Science, Laurentian University, Sudbury, ON, Canada and Zhangli Peng, Department of Biomedical Engineering, University of Illinois at Chicago for critically reading the manuscript.

Conflicts of Interest

The author states no conflict of interest.

Abbreviations

Text
two-dimensional2D
bilayer couple bendingbcb
hemoglobin free red cellghost
single layer bendingslb
tank-treadTT
Quantities
bending stiffness of the membrane in bcb B bc
frictional coefficient opposing a lateral motion between two layers c f
2D shear strain of the membrane skeleton ϵ
Newtonian membrane viscosity η mem ¯
non Newtonian membrane viscosity η mem
viscosity of the suspending phase η 0
shear rate of the undisturbed simple shear flow γ ˙
distance between the neutral planes of the two monolayers of a bilayerh
biaxial modulus of the membraneK
shear modulus of the membrane μ
undisturbed density of Band 3 complex ρ 0
actual density of Band 3 complex ρ
shear stress opposing a lateral motion between skeleton and bilayer σ sb
characteristic time of relaxation from an extension τ
TT frequency T T F
slip velocity between skeleton and bilayer v sb

References

  1. Heinrich, V.; Ritchie, K.; Mohandas, N.; Evans, E. Elastic Thickness Compressibilty of the Red Cell Membrane. Biophys. J. 2001, 81, 1452–1463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Matteoli, P.; Nicoud, F.; Mendez, S. Impact of the membrane viscosity on the tank-treading behavior of red blood cells. Phys. Rev. Fluids 2021, 6, 043602. [Google Scholar] [CrossRef]
  3. Tsubota, K. Elongation deformation of a red blood cell under shear flow as stretch testing. J. Mech. Phys. Solids 2021, 152, 104345. [Google Scholar] [CrossRef]
  4. Guglietta, F.; Behr, M.; Biferale, L.; Falcucci, G.; Sbragaglia, M. On the effects of membrane viscosity on transient red blood cell dynamics. Soft Matter 2020, 16, 6191–6205. [Google Scholar] [CrossRef]
  5. Rezghi, A.; Zhang, J. Tank-treading dynamics of red blood cells in shear flow: On the membrane viscosity rheology. Biophys. J. 2022, 121, 3393–3410. [Google Scholar] [CrossRef] [PubMed]
  6. Evans, E.; La Celle, P. Intrinsic material properties of the erythrocyte membrane indicated by mechanical analysis of deformation. Blood 1975, 45, 29–43. [Google Scholar] [CrossRef] [Green Version]
  7. Evans, E.A.; Hochmuth, R.M. Membrane Viscoelasticity. Biophys. J. 1976, 16, 1–11. [Google Scholar] [CrossRef] [Green Version]
  8. Chien, S.; Sung, K.; Skalak, R.; Usami, S.; Tözeren, A. Theoretical and experimental studies on viscoelastic properties of erythrocyte membrane. Biophys. J. 1978, 24, 463–487. [Google Scholar] [CrossRef] [Green Version]
  9. Tözeren, A.; Skalak, R.; Fedorciw, B.; Sung, K.; Chien, S. Constitutive equations of erythrocyte membrane incorporating evolving preferred configuration. Biophys. J. 1984, 45, 541–549. [Google Scholar] [CrossRef] [Green Version]
  10. Fischer, T.M. Is the surface area of the red cell membrane skeleton locally conserved? Biophys. J. 1992, 61, 298–305. [Google Scholar] [CrossRef]
  11. Hochmuth, R.M.; Worthy, P.R.; Evans, E.A. Red Cell Extensional Recovery and the Determination of Membrane Viscosity. Biophys. J. 1979, 26, 101–114. [Google Scholar] [CrossRef] [Green Version]
  12. Tran-Son-Tay, R.; Sutera, S.P.; Rao, P.R. Determination of Red Blood Cell Membrane Viscosity from Rheoscopic Observations of Tank-Treading Motion. Biophys. J. 1984, 46, 65–72. [Google Scholar] [CrossRef] [Green Version]
  13. Fischer, T.M.; Stöhr-Liesen, M.; Schmid-Schönbein, H. The Red Cell as a Fluid Droplet: Tank Tread-Like Motion of the Human Erythrocyte Membrane in Shear Flow. Science 1978, 202, 894–896. [Google Scholar] [CrossRef]
  14. Schmid-Schönbein, H.; Wells, R. Fluid Drop-Like Transition of Erythrocytes under Shear. Science 1969, 165, 288–291. [Google Scholar] [CrossRef]
  15. Keller, S.R.; Skalak, R. Motion of a tank-treading ellipsoidal particle in a shear flow. J. Fluid Mech. 1982, 120, 27–47. [Google Scholar] [CrossRef]
  16. Sutera, S.P.; Pierre, P.R.; Zahalak, G.I. Deduction of intrinsic mechanical properties of the erythrocyte membrane from observations of tank-treading in the rheoscope. Biorheology 1989, 26, 177–197. [Google Scholar] [CrossRef]
  17. Waugh, R.E. Surface Viscosity Measurements From Large Bilayer Vesicle Tether Formation II. Experiments. Biophys. J. 1982, 38, 29–37. [Google Scholar] [CrossRef] [Green Version]
  18. Fischer, T.M. Tank-Tread Frequency of the Red Cell Membrane: Dependence on the Viscosity of the Suspending Medium. Biophys. J. 2007, 93, 2553–2561. [Google Scholar] [CrossRef] [Green Version]
  19. Mills, J.P.; Qie, L.; Dao, M.; Lim, C.T.; Suresh, S. Nonlinear elastic and viscoelastic deformation of the human red blood cell with optical tweezers. Mech. Chem. Biosys. 2004, 1, 169–180. [Google Scholar]
  20. Fischer, T.M.; Korzeniewski, R. Angle of Inclination of Tank-Treading Red Cells: Dependence on Shear Rate and Suspending Medium. Biophys. J. 2015, 108, 1352–1360. [Google Scholar] [CrossRef] [Green Version]
  21. Fischer, T.M.; Korzeniewski, R. Effects of shear rate and suspending medium viscosity on elongation of red cells tank-treading in shear flow. Cytom. Part A 2011, 79A, 946–951. [Google Scholar] [CrossRef]
  22. Stokke, B.T.; Mikkelsen, A.; Elgsaeter, A. The human erythrocyte membrane skeleton may be an ionic gel. I. Membrane mechanochemical properties. Eur. Biophys. J. 1986, 13, 203–218. [Google Scholar] [CrossRef]
  23. Stokke, B.T.; Mikkelsen, A.; Elgsaeter, A. The human erythrocyte membrane skeleton may be an ionic gel. III. Micropipette aspiration of unswollen erythrocytes. J. Theor. Biol. 1986, 123, 205–211. [Google Scholar] [CrossRef]
  24. Discher, D.E.; Mohandas, N.; Evans, E.A. Molecular Maps of Red Cell Deformation: Hidden Elasticity and in Situ Connectivity. Science 1994, 226, 1032–1035. [Google Scholar] [CrossRef] [Green Version]
  25. Saffman, P.G. Brownian motion in thin sheets of viscous fluid. J. Fluid Mech. 1976, 73, 593–602. [Google Scholar] [CrossRef] [Green Version]
  26. Peng, Z.; Asaro, R.J.; Zhu, Q. Multiscale modelling of erythrocytes in Stokes flow. J. Fluid Mech. 2011, 686, 299–337. [Google Scholar] [CrossRef] [Green Version]
  27. Peng, Z.; Zhu, Q. Deformation of the erythrocyte cytoskeleton in tank treading motions. Soft Matter 2013, 9, 7617–7627. [Google Scholar] [CrossRef]
  28. Kodippili, G.C.; Spector, J.; Sullivan, C.; Kuypers, F.A.; Labotka, R.; Gallagher, P.G.; Ritchie, K.; Low, P.S. Imaging of the diffusion of single band 3 molecules on normal and mutant erythrocytes. Blood 2009, 113, 6237–6245. [Google Scholar] [CrossRef] [Green Version]
  29. Turlier, H.; Fedosov, D.A.; Audoly, B.; Auth, T.; Gov, N.S.; Sykes, C.; Joanny, J.F.; Gompper, G.; Betz, T. Equilibrium physics breakdown reveals the active nature of red blood cell flickering. Nat. Phys. 2016, 12, 513–519. [Google Scholar] [CrossRef] [Green Version]
  30. Saffman, P.G.; Delbrück, M. Brownian motion in biological membranes. Proc. Natl. Acad. Sci. USA 1975, 72, 3111–3113. [Google Scholar] [CrossRef] [Green Version]
  31. Evans, E.A.; Skalak, R. Mechanics and Thermodynamics of Biomembranes. CRC Crit. Rev. Bioeng. 1979, 3, 181–418. [Google Scholar]
  32. Helfrich, W. Elastic properties of lipid bilayers: Theory and possible experiments. Z. Naturforsch. C Biosci. 1973, 28, 693–703. [Google Scholar] [CrossRef]
  33. Servuss, R.; Harbich, V.; Helfrich, W. Measurement of the curvature-elastic modulus of egg lecithin bilayers. Biochim. Biophys. Acta (BBA)—Biomembr. 1976, 436, 900–903. [Google Scholar] [CrossRef]
  34. Duwe, H.P.; Zeman, K.; Sackmann, E. Bending undulations of lipid bilayers and the red blood cell membrane: A comparative study. Prog. Colloid polym. Sci. 1989, 79, 6–10. [Google Scholar]
  35. Waugh, R.E.; Song, J.; Svetina, S.; Zeks, B. Local and nonlocal curvature elasticity in bilayer membranes by tether formation from lecithin vesicles. Biophys. J. 1992, 61, 974–982. [Google Scholar] [CrossRef] [Green Version]
  36. Hwang, W.; Waugh, R. Energy of dissociation of lipid bilayer from the membrane skeleton of red blood cells. Biophys. J. 1997, 72, 2669–2678. [Google Scholar] [CrossRef] [Green Version]
  37. Fischer, T.M. Bending stiffness of lipid bilayers. I. Bilayer couple or single-layer bending? Biophys. J. 1992, 63, 1328–1335. [Google Scholar] [CrossRef] [Green Version]
  38. Drummond, J.; Tahir, M. Laminar viscous flow through regular arrays of parallel solid cylinders. Int. J. Multiph. Flow 1984, 10, 515–540. [Google Scholar] [CrossRef]
  39. Bo, L.; Waugh, R. Determination of bilayer membrane bending stiffness by tether formation from giant, thin-walled vesicles. Biophys. J. 1989, 55, 509–517. [Google Scholar] [CrossRef] [Green Version]
  40. Fischer, T.M. On the Energy Dissipation in a Tank-Treading Human Red Blood Cell. Biophys. J. 1980, 61, 863–868. [Google Scholar] [CrossRef] [Green Version]
  41. Evans, E.; Yeung, A. Hidden dynamics in rapid changes of bilayer shape. Chem. Phys. Lipids 1994, 73, 39–56. [Google Scholar] [CrossRef]
  42. Khoshnood, A.; Noguchi, H.; Gompper, G. Lipid membranes with transmembrane proteins in shear flow. J. Chem. Phys. 2010, 132, 025101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Rief, M.; Pascual, J.; Saraste, M.; Gaub, H.E. Single molecule force spectroscopy of spectrin repeats: Low unfolding forces in helix bundles. J. Mol. Biol. 1999, 286, 553–561. [Google Scholar] [CrossRef] [PubMed]
  44. Waugh, R. Effects of 2,3-diphosphoglycerate on the mechanical properties of erythrocyte membrane. Blood 1986, 68, 231–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nash, G.; Meiselman, H. Red cell and ghost viscoelasticity. Effects of hemoglobin concentration and in vivo aging. Biophys. J. 1983, 43, 63–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Nans, A.; Mohandas, N.; Stokes, D.L. Native Ultrastructure of the Red Cell Cytoskeleton by Cryo-Electron Tomography. Biophys. J. 2011, 101, 2341–2350. [Google Scholar] [CrossRef] [Green Version]
  47. Burton, N.; Bruce, L. Modelling the structure of the red cell membrane. Biochem. Cell Biol. Biochim. Biol. Cell. 2011, 89, 200–215. [Google Scholar] [CrossRef]
  48. Feng, Z.; Waugh, R.E.; Peng, Z. Constitutive Model of Erythrocyte Membranes with Distributions of Spectrin Orientations and Lengths. Biophys. J. 2020, 119, 2190–2204. [Google Scholar] [CrossRef]
  49. Discher, D.E.; Mohandas, N. Kinematics of Red Cell Aspiration by Fluorescence-Imaged Microdeformation. Biophys. J. 1996, 71, 1680–1694. [Google Scholar] [CrossRef] [Green Version]
  50. Skalak, R.; Tözeren, A.; Zarda, P.; Chien, S. Strain Energy Function of Red Blood Cell Membranes. Biophys. J. 1973, 13, 245–264. [Google Scholar] [CrossRef] [Green Version]
  51. Dimitrakopoulos, P. Analysis of the variation in the determination of the shear modulus of the erythrocyte membrane: Effects of the constitutive law and membrane modeling. Phys. Rev. E 2012, 85, 041917. [Google Scholar] [CrossRef] [Green Version]
  52. Lenormand, G.; Hénon, S.; Richert, A.; Siméon, J.; Gallet, F. Direct Measurement of the Area Expansion and Shear Moduli of the Human Red Blood Cell Membrane Skeleton. Biophys. J. 2001, 81, 43–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Fung, Y.C.; Tsang, W.C.O.; Patitucci, P. High-Resolution Data on the Geometry of Red Blood Cells. Biorheology 1981, 18, 369–385. [Google Scholar] [CrossRef] [PubMed]
  54. Fischer, T.M. The Shape of Human Red Blood Cells Suspended in Autologous Plasma and Serum. Cells 2022, 11, 1941. [Google Scholar] [CrossRef]
  55. Fischer, T.M. Creep and stress relaxation of human red cell membrane. Biomech. Model. Mechanobiol. 2017, 16, 239–247. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic of the kinematics of the tank-tread (TT) motion. (A) View along the vorticity of a simple shear flow (middle). Left, a red cell at rest in edge-on orientation. Right, TTing under the imposed shear flow. The orientation of the elongated cell is almost stable due to a compensation of two turning moments [15]. The (black) arrow indicates the motion of the membrane. The dimple region is marked by a membrane-attached sphere (dark green). (B) View along the gradient of the simple shear flow. In an experiment, the TT motion is visualized by a back-and-forth motion of membrane-attached spheres [13].
Figure 1. Schematic of the kinematics of the tank-tread (TT) motion. (A) View along the vorticity of a simple shear flow (middle). Left, a red cell at rest in edge-on orientation. Right, TTing under the imposed shear flow. The orientation of the elongated cell is almost stable due to a compensation of two turning moments [15]. The (black) arrow indicates the motion of the membrane. The dimple region is marked by a membrane-attached sphere (dark green). (B) View along the gradient of the simple shear flow. In an experiment, the TT motion is visualized by a back-and-forth motion of membrane-attached spheres [13].
Biomolecules 13 00130 g001
Figure 2. (A) Schematic of a triangle in the idealized network of the skeleton. The circles indicate the cross sections of the intrinsic portions of Band 3 complexes and Band 3 dimers drawn to scale. The diameters are 18.5, 15.8, 6.2 nm for the junction complex (blue), suspension complex (green), and Band 3 dimer (brown), respectively. In vivo, the six hetero dimers of spectrin connecting junction complexes and suspension complexes are strongly convoluted and rise above the bilayer plane [46]. Here, they are indicated by straight lines. (B) Schematic of two hexagons each composed of six triangles as in A. The left hexagon is elongated as in a static elongation at the start of a relaxation process. The right one is relaxed indicating the situation at the end of the relaxation process. The red arrows in the left hexagon indicate the direction of motion of the Band 3 complexes during the relaxation. (C) Schematic of the fate of a hexagon during the aspiration of a membrane tongue into a small bore pipette. The straight arrows indicate the direction of motion of the Band 3 complexes during the aspiration. The curved arrow indicates the displacement of the hexagon to its final position at static equilibrium.
Figure 2. (A) Schematic of a triangle in the idealized network of the skeleton. The circles indicate the cross sections of the intrinsic portions of Band 3 complexes and Band 3 dimers drawn to scale. The diameters are 18.5, 15.8, 6.2 nm for the junction complex (blue), suspension complex (green), and Band 3 dimer (brown), respectively. In vivo, the six hetero dimers of spectrin connecting junction complexes and suspension complexes are strongly convoluted and rise above the bilayer plane [46]. Here, they are indicated by straight lines. (B) Schematic of two hexagons each composed of six triangles as in A. The left hexagon is elongated as in a static elongation at the start of a relaxation process. The right one is relaxed indicating the situation at the end of the relaxation process. The red arrows in the left hexagon indicate the direction of motion of the Band 3 complexes during the relaxation. (C) Schematic of the fate of a hexagon during the aspiration of a membrane tongue into a small bore pipette. The straight arrows indicate the direction of motion of the Band 3 complexes during the aspiration. The curved arrow indicates the displacement of the hexagon to its final position at static equilibrium.
Biomolecules 13 00130 g002
Table 1. Conflicting results in η mem ¯ upon fitting the same experimental data.
Table 1. Conflicting results in η mem ¯ upon fitting the same experimental data.
η 0 (mPa s)5454Axis Ratio of the
γ ˙ (1/s)8130Stress Free ShapeReference
η mem ¯ ( 10 8  Pa s m)244.50.9Matteoli et al. [2]
η mem ¯ ( 10 8  Pa s m)4019.50.84Tsubota [3]
Table 2. Comparison of the η mem ¯ values from different studies.
Table 2. Comparison of the η mem ¯ values from different studies.
Fitted ParameterReference η mem ¯ ( 10 8  Pa s m)
τ Hochmuth et al. [11]60–80
T T F Tran-Son-Tay et al. [12]10–8.6
τ Guglietta et al. [4]42
T T F Matteoli et al. [2], Tsubota [3]4.5–40
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fischer, T.M. Energy Dissipation in the Human Red Cell Membrane. Biomolecules 2023, 13, 130. https://doi.org/10.3390/biom13010130

AMA Style

Fischer TM. Energy Dissipation in the Human Red Cell Membrane. Biomolecules. 2023; 13(1):130. https://doi.org/10.3390/biom13010130

Chicago/Turabian Style

Fischer, Thomas M. 2023. "Energy Dissipation in the Human Red Cell Membrane" Biomolecules 13, no. 1: 130. https://doi.org/10.3390/biom13010130

APA Style

Fischer, T. M. (2023). Energy Dissipation in the Human Red Cell Membrane. Biomolecules, 13(1), 130. https://doi.org/10.3390/biom13010130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop