The coronavirus is always idealized as a spherical capsid with radially protruding spikes. However, histologically, in the tissues of infected patients, capsids in cross section are elliptical, and only sometimes spherical [Neuman et al., “Supramolecular architecture of severe acute respiratory syndrome coronavirus revealed by electron cryomicroscopy,” J Virol, 80, 7918 (2006)]. This capsid ellipticity implies that coronaviruses are oblate or prolate or both. We call this diversity of shapes, pleomorphism. Recently, the rotational diffusivity of the spherical coronavirus in suspension was calculated, from first principles, using general rigid bead-rod theory [Kanso et al., “Coronavirus rotational diffusivity,” Phys Fluids 32, 113101 (2020)]. We did so by beading the spherical capsid and then also by replacing each of its bulbous spikes with a single bead. In this paper, we use energy minimization for the spreading of the spikes, charged identically, over the oblate or prolate capsids. We use general rigid bead-rod theory to explore the role of such coronavirus cross-sectional ellipticity on its rotational diffusivity, the transport property around which its cell attachment revolves. We learn that coronavirus ellipticity drastically decreases its rotational diffusivity, be it oblate or prolate.

The coronavirus is thought of as a spherical capsid with radially protruding spikes. However, histologically, in the tissues of infected patients, capsids in cross section are aspherical, roughly elliptical, and only sometimes spherical. For instance, its capsid may be elliptical when microtomed in preparation for microscopy (see all 35 panels of Fig. 2 of Ref. 1, which we reproduce here in Fig. 1). This capsid ellipticity implies that coronaviruses are oblate or prolate or both. We call this diversity of shapes, pleomorphism. We know of no live microscopy of the coronavirus, in suspension or otherwise. We thus arrive at our understanding of pleomorphism through the colored lenses of electron microscopy, be it cryogenic fractography1 or particle staining (Fig. 3. of Ref. 2). This is true whether the coronavirus capsid is axisymmetric or not. In cross section, as in microscopic imagery of microtomed coronavirus-infected tissue, a suspension of spherical capsids presents as circles, just circles. Further, when axisymmetric ellipsoidal coronavirus is sectioned, normal to its major axis (prolate), or to its minor axis (oblate), these cuts are also circular. Hence, in cross section, a suspension of aspherical capsids presents as both circular and acircular loops. More specifically, a suspension of ellipsoidal capsids presents as both circles and ellipses. There is thus still much to be learned experimentally about the shape of the coronavirus. This paper is devoted to how its capsid shapes affect its transport properties.

FIG. 1.

Cryo-electron micrographs of pleomorphic SARS-CoV virions (Fig 2 of Ref. 1) Reproduced with permission from Neuman et al., J. Virol. 80, 7918 (2006).

FIG. 1.

Cryo-electron micrographs of pleomorphic SARS-CoV virions (Fig 2 of Ref. 1) Reproduced with permission from Neuman et al., J. Virol. 80, 7918 (2006).

Close modal

General rigid bead-rod theory relies exclusively on macromolecular orientation to explain the rheological properties.3–7 This distinguishes general rigid bead-rod theory from its competing approaches, which include reptation or disentanglement. We refer the reader to Ref. 8 for the detailed derivation of the general rigid bead-rod theory, and specifically, to Sec. III of Ref. 8. We define the characteristic length, L, of our coronavirus bead-rod models as the separation of nearest bead centers. Our general rigid bead-rod theory symbols, dimensional and non-dimensional respectively, are listed in Tables I and II, which follow those of the corresponding textbook treatments (EXAMPLE 16.7-1 of Ref. 9 or EXAMPLE 13.6-1 of Ref. 10). We are attracted to general rigid bead-rod theory first, for its flexibility. We are attracted to general rigid bead-rod theory second, for the accuracy of its simplest special case, the rigid dumbbell, for which many rheological material functions are properly predicted [see romanettes (i)–(xvi) of Sec. I of Ref. 11].

TABLE I.

Dimensional variables. Legend: M—mass, L—length, and t—time.

Name Unit Symbol
Angular frequency  t 1   ω  
Augmented energy functional  M L 2 / t 2   E  
Bead friction coefficient  M / t   ζ  
Capsid radius  L   r c  
Complex viscosity  M / L t   η  
Dielectric permittivity  T 4 I 2 / M L 3   ϵ  
Kinetic molecular energy per molecule  M L 2 / t 2   k T  
Length of the spike of each peplomer  L    
Minus the imaginary part of the complex viscosity  M / L t   η  
Number of dumbbells per unit volume  1 / L 3   n  
Point charge  A s   Q  
Real part of the complex viscosity  M / L t   η  
Relaxation time of rigid dumbbell  t   λ 0  
Relaxation time of solution  t   λ  
Rotational diffusivity  t 1   D r  
Rotatory diffusivity  L 2 / t   D rot  
Shear rate amplitude  t 1   γ ̇ 0  
Solvent viscosity  M / L t   η s  
Total electrostatic energy  M L 2 / t 2   E  
Translational diffusivity  L 2 / t   D t r  
Virus radius  L   r v  
Viscosity, zero-shear  M / L t   η 0  
Zero-shear first normal stress difference  M / L   Ψ 1 , 0  
Name Unit Symbol
Angular frequency  t 1   ω  
Augmented energy functional  M L 2 / t 2   E  
Bead friction coefficient  M / t   ζ  
Capsid radius  L   r c  
Complex viscosity  M / L t   η  
Dielectric permittivity  T 4 I 2 / M L 3   ϵ  
Kinetic molecular energy per molecule  M L 2 / t 2   k T  
Length of the spike of each peplomer  L    
Minus the imaginary part of the complex viscosity  M / L t   η  
Number of dumbbells per unit volume  1 / L 3   n  
Point charge  A s   Q  
Real part of the complex viscosity  M / L t   η  
Relaxation time of rigid dumbbell  t   λ 0  
Relaxation time of solution  t   λ  
Rotational diffusivity  t 1   D r  
Rotatory diffusivity  L 2 / t   D rot  
Shear rate amplitude  t 1   γ ̇ 0  
Solvent viscosity  M / L t   η s  
Total electrostatic energy  M L 2 / t 2   E  
Translational diffusivity  L 2 / t   D t r  
Virus radius  L   r v  
Viscosity, zero-shear  M / L t   η 0  
Zero-shear first normal stress difference  M / L   Ψ 1 , 0  
TABLE II.

Dimensionless variables and groups.

Name Symbol
Aspect ratio  ε c a 1  
Capsid sphere  C  
Coefficient in Eqs. (9) and (10)  a  
Coefficient in Eqs. (9) and (10)(9)  b  
Coefficient in Eqs. (9) and (10)  ν  
Deborah number, oscillatory shear  De λ ω  
Sphere  S  
Total number of beads  N  
Total number of capsid beads  N c  
Total number of peplomers  N p  
Volume fraction  Φ  
Weissenberg number  Wi λ γ ̇ 0  
Name Symbol
Aspect ratio  ε c a 1  
Capsid sphere  C  
Coefficient in Eqs. (9) and (10)  a  
Coefficient in Eqs. (9) and (10)(9)  b  
Coefficient in Eqs. (9) and (10)  ν  
Deborah number, oscillatory shear  De λ ω  
Sphere  S  
Total number of beads  N  
Total number of capsid beads  N c  
Total number of peplomers  N p  
Volume fraction  Φ  
Weissenberg number  Wi λ γ ̇ 0  

Recently, we calculated the rotational diffusivity of the spherical coronavirus in suspension, from first principles, using general rigid bead-rod theory.12–14 We did so by beading the spherical capsid and then also by replacing each of its bulbous spikes with a single bead (see Fig. 5 of Ref. 12). In this paper, we use energy minimization for the spreading of the spikes, charged identically, over the oblate or prolate capsids (Sec. IV). We use general rigid bead-rod theory to explore the role of such coronavirus cross-sectional ellipticity on its rotational diffusivity, the transport property around which its cell attachment revolves. We learn that coronavirus ellipticity decreases its rotational diffusivity for both oblate and prolate capsids.

In the tradition of the transport sciences, we define the rotatory diffusivity as (see Footnote 2 of p. 62 of Ref. 9)

D rot 2 k T ζ ,
(1)

which, for any axisymmetric macromolecule, from general rigid bead-rod theory, gives

D rot 12 L 2 ν D r ,
(2)

which has the dimensions of diffusivity and which is four times the translational diffusivity

D rot 4 D t r
(3)

or

D r ν 3 L 2 D t r .
(4)

In this paper, we depart from said transport tradition of using the rotatory diffusivity, D rot, and frame our results in terms of the rotational diffusivities, D r, of pleomorphic coronavirus particles.

In general rigid bead-rod theory, we construct macromolecules from sets of beads whose positions, relative to one another, are fixed. Our macromolecular bead-rod models of our pleomorphic coronavirus particles are suspended in a Newtonian solvent. In this work, we neglect interactions of the solvent velocity fields, be they between nearest beads,15,16 or nearest macromolecules. With general rigid bead-rod theory, we thus locate beads to sculpt an approximation of the pleomorphic coronavirus particle shapes. In this way, using general rigid bead-rod theory, we can model any virus macromolecular architecture (see Fig. 9 of Ref. 11).

We use Eqs. (3)–(13) in Ref. 12 for the method of computing the rotational diffusivity (see Footnote 2 of p. 62 of Ref. 9)

D r 1 6 λ
(5)

or [Eq. (23) of Ref. 12]

λ 0 D r = ν 72 ,
(6)

which we will use for our results below.

In this paper, we focus on small-amplitude oscillatory shear flow (SAOS). For this flow field, for the molecular definition of small amplitude, general rigid bead-rod theory yields [Eq. (32) of Ref. 12)

λ γ ̇ 0 1 ν 2 ,
(7)

whose left side is the macromolecular Weissenberg number.

The polymer contributions to the complex viscosity17,18

η * η i η
(8)

are [Eqs. (40) and (41) of Ref. 11]

η η s η 0 η s = 1 2 b / a ν + 1 1 1 2 b / a ν + 1 1 + λ ω 2
(9)

and

η η 0 η s = 1 2 b / a ν + 1 1 λ ω 1 + λ ω 2 ,
(10)

where λ ω is the Deborah number. Equations (9) and (10) each capture non-Newtonian behavior: (i) Eq. (9) captures the descent of η ω and (ii) Eq. (10) captures the ascent of η ω from the origin. In this paper, we plot the real and minus the imaginary parts of the shear stress responses to small-amplitude oscillatory shear flow as functions of frequency, following Ferry (Secs. 2.A.4–2.A.6 of Ref. 19) or Bird et al. (Sec. 4.4 of Ref. 20):

As ω 0, for the polymer contribution to the zero-shear viscosity, we get

η 0 η s nkT λ = a ν 2 + b = b 1 + 2 b a ν 2 b a ν 1
(11)

and for the zero-shear first normal stress difference coefficient

Ψ 1 , 0 λ η 0 η s = 2 1 2 b / a ν + 1 1 ,
(12)

which we use in the table of Sec. V. below.

Coronavirus peplomers are charged identically and anchored into a lipid bilayer, and are thus displaced by the repulsions of their nearest neighbors. In this paper, we represent peplomers with single beads. We then locate them by applying the energy minimization scheme of Sec. IV. of Ref. 13 to singly beaded peplomers repelling one another over ellipsoidal surfaces of the pleomorphic virus.

From the literature, we learn that if oblate, the coronavirus shape range is [Figs. 1(a) and 2 of Ref. 1]

2 3 c a 4 5 ,
(13)

and if prolate [Figs. 1(a) and 2 of Ref. 1]

5 4 c a 3 2 .
(14)

For this work, we therefore cover [Eqs. (13)plus(14)]

5 4 c a 4 5 ,
(15)

where, c is the ellipsoidal whole-particle length along the δ3 molecular axis, and a along the δ1. By whole-particle, we mean capsid plus peplomer. Specifically, for this paper, we straddle both ranges Eqs. (13) and (14) with the following set of pleomorphic coronavirus aspect ratios

c a = 2 3 , 4 5 , 1 , 5 4 , 3 2
(16)

for which, after energy minimization, we construct the bead-rod models of Figures 2–6 (Multimedia views) for N c = 256 and N p = 74. For the N c = 256, we rely on the capsid beading study (see subsection VII A of Ref. 12). For the average number of peplomers, N p = 74, we rely on our previous literature review (see Table X of Ref. 12).

FIG. 2.

General rigid bead-rod model of prolate coronavirus, N c = 256, N p = 74, and c / a = 5 / 4. Multimedia view: https://doi.org/10.1063/5.0094771.1

FIG. 2.

General rigid bead-rod model of prolate coronavirus, N c = 256, N p = 74, and c / a = 5 / 4. Multimedia view: https://doi.org/10.1063/5.0094771.1

Close modal
FIG. 3.

General rigid bead-rod model of prolate coronavirus, N c = 256, N p = 74, and c / a = 3 / 2. Multimedia view: https://doi.org/10.1063/5.0094771.2

FIG. 3.

General rigid bead-rod model of prolate coronavirus, N c = 256, N p = 74, and c / a = 3 / 2. Multimedia view: https://doi.org/10.1063/5.0094771.2

Close modal
FIG. 4.

General rigid bead-rod model of spherical coronavirus, N c = 256, N p = 74, and c / a = 1. Multimedia view: https://doi.org/10.1063/5.0094771.3

FIG. 4.

General rigid bead-rod model of spherical coronavirus, N c = 256, N p = 74, and c / a = 1. Multimedia view: https://doi.org/10.1063/5.0094771.3

Close modal
FIG. 5.

General rigid bead-rod model of oblate coronavirus, N c = 256, N p = 74, and c / a = 4 / 5. Multimedia view: https://doi.org/10.1063/5.0094771.4

FIG. 5.

General rigid bead-rod model of oblate coronavirus, N c = 256, N p = 74, and c / a = 4 / 5. Multimedia view: https://doi.org/10.1063/5.0094771.4

Close modal
FIG. 6.

General rigid bead-rod model of oblate coronavirus, N c = 256, N p = 74, and c / a = 2 / 3. Multimedia view: https://doi.org/10.1063/5.0094771.5

FIG. 6.

General rigid bead-rod model of oblate coronavirus, N c = 256, N p = 74, and c / a = 2 / 3. Multimedia view: https://doi.org/10.1063/5.0094771.5

Close modal

We next employ the framework developed by Chaurasia et al.21 (see also Chaurasia22) to find equilibrium solutions of a system consisting of flexible structures, specifically charged elastic loops constrained to a sphere. We do so for identical point charges spreading over the surfaces of oblate and prolate ellipsoids.

Let C be the ellipsoidal capsid with axes lengths a, b, and c. Let o be the origin of a Cartesian coordinate system with an orthonormal basis i , j , k, such that the center of C coincides with the origin o and the axes of C are along the vectors i, j, and k, respectively. Let N p be the number of single bead peplomers of identical spike length attached along the normal to the ellipsoid C at the point of contact. The point of contact of the spike of the i th bead to C is defined by

x i i + y i j + z i k ,
(17)

where i =1,2,…,N. Then, the quantities x i, y i, and z i must satisfy

x i 2 a 2 + y i 2 b 2 + z i 2 c 2 = 1.
(18)

We assume that each peplomer spike is attached to the ellipsoidal capsid C along the normal to C. The normal vector at the point of contact x i i + y i j + z i k of the spike of the ith peplomer to C is given by

n i x i a 2 i + y i b 2 j + z i c 2 k = 1.
(19)

Thus, introducing a scalar quantity t i, i = 1 , 2 , , N, the position vector r i of the i th bead is given by

r i x i i + y i j + z i k + t i n i ,
(20)

where, by using Eq. (19), Eq. (20) simplifies to

r i x i 1 + t i a 2 i + y i 1 + t i b 2 i + z i 1 + t i c 2 i ,
(21)

where i = 1 , 2 , , N. We note from Eq. (21) that the length of the spike of the ith peplomer is given by r i x i i + y i j + z i k = t i n i, i = 1 , 2 , , N. For simplification, we assume that the length of the spike of each peplomer is equal to ; thus, the quantities x i, y i, z i, and t i and Eq. (19) of n i in t i n i 2 = 2 must satisfy

t i 2 x i 2 a 4 + y i 2 b 4 + z i 2 c 4 = 2 ,
(22)

where i = 1 , 2 , , N. Thus, the position vector r i, defined by Eq. (21), of the ith peplomer bead of a given spike length is entirely determined in terms of four scalar quantities x i, y i, z i, and t i satisfying Eqs. (18) and (22). For an axisymmetric capsid, a = b.

Let each single bead peplomer be endowed with a point charge Q. The total electrostatic energy of N p peplomers, constrained to the ellipsoidal capsid C, is given by

E = Q 4 π ϵ i = 1 N j = 1 N 1 r i r j ,
(23)

where ϵ is dielectric permittivity and r i, defined in Eq. (21), is the position vector of the i th peplomer.23 

We have assumed that each bead is endowed with an identical point charge Q for simplification. Therefore, the beads repel each other and would prefer to distribute themselves as far as possible from each other to minimize the electrostatic energy E, defined in Eq. (23). Using a constrained minimization approach, we find an equilibrium distribution of the beads, defined in Eq. (20), that locally minimizes the energy in Eq. (23) while satisfying the kinematic constraints in Eqs. (18) and (22), for given values of N p. Since the charge Q appears only as a prefactor in Eq. (23), its value plays no role in determining equilibrium solutions. Dropping that prefactor, we define the augmented energy functional

E ̂ r i , Λ i , λ i i = 1 , N p = i = 1 N p j = 1 j i N p 1 r i r j + 1 2 i = 1 N p Λ i x i 2 a 2 + y i 2 b 2 + z i 2 c 2 1 + 1 2 i = 1 N p Λ i t i 2 x i 2 a 4 + y i 2 b 4 + z i 2 c 4 2 ,
(24)

where using Eq. (21)

r i r j = x i 1 + t i a 2 x j 1 + t j a 2 2 + y i 1 + t i b 2 y j 1 + t j b 2 2 + z i 1 + t i c 2 z j 1 + t j c 2 2 ,
(25)

where Λ i and λ i are the Lagrange multipliers introduced to satisfy the kinematic constraints in Eqs. (18) and (22). For finding the distribution of the beads locally minimizing the electrostatic energy E, defined in Eq. (23), we differentiate the augmented energy functional, E ̂, with scalar quantities x i, y i, z i, and t i, where i = 1 , , N p, resulting in 4 N p equilibrium equations

j = 1 j i N p r i r j i 1 + t i a 2 r i r j 3 / 2 Λ i x i a 2 λ i t i 2 x i a 4 = 0 ,
(26)
j = 1 j i N p r i r j j 1 + t i b 2 r i r j 3 / 2 Λ i y i b 2 λ i t i 2 y i b 4 = 0 ,
(27)
j = 1 j i N p r i r j k 1 + t i c 2 r i r j 3 / 2 Λ i z i c 2 λ i t i 2 z i c 4 = 0 ,
(28)

and

j = 1 j i N p r i r j n i r i r j 3 / 2 λ i t i x i 2 a 4 + y i 2 b 4 + z i 2 c 4 = 0 ,
(29)

respectively, where n i and r i are defined by Eqs. (19) and (20). In total, we solve 6 N p equations, 4 N p equilibrium equations in Eqs. (26)–(29), and 2 N p constraints in Eqs. (18) and (22), simultaneously to determine 6 N p unknowns in x i, y i, z i, t i, λ i, and Λ i, where i = 1 , , N p. We use the Levenberg–Marquardt algorithm from the fsolve package of MATLAB to solve the system of equations with 10 16 error tolerance.

General rigid bead-rod theory can be used either structure-by-structure (see TABLES V–XIII of Ref. 11) or analytically (see TABLE XV of Ref. 11). For large values of N p, such as N p = 74, Eq. (23) is applied numerically. Thus, the bead positions R i, are not derived, but rather we arrive at their floating-point approximations. Our exploration of pleomorphism is thus structure-by-structure.

Our model, which has an ellipsoidal core, amounts to a modified version of the Thomson problem,24 wherein one seeks to find a state that distributes N p electrons over a unit sphere as evenly as possible, with minimum electrostatic energy. Wales et al.,25,26 solved this problem, providing solutions for a large set of values of N p. Our energy minimization recovers accurately the results of Wales25,26 for the Thomson solution (energy minimization over a spherical surface), as it should (see Sec. VII. of Ref. 12). As far as we know, we are the first to perform such an energy minimization over the surface of an ellipsoid.

Table III summarizes our results for the characteristics of ellipsoidal coronavirus particles, be they oblate or prolate, arrived at from general rigid bead-rod theory for N c = 256 and N p = 74. Figure 7 maps columns 2 and 3 of Table III onto the I 3 I 1 plane, showing the balance of moments of each pleomorphic coronavirus bead-rod model. The oblate ones lie above the diagonal, and the prolate, below. The spherically symmetric coronavirus lies on the diagonal, near the origin. Figure 8 shows that pleomorphism, be it oblate or prolate, causes the real part of the complex viscosity, η ω, to descend with frequency. This is not seen for the spherical capsid (black horizontal line). In other words, pleomorphism introduces non-Newtonian behavior. Further, when compared with Fig. 10 of Ref. 12, Fig. 9 shows that pleomorphism, be it oblate or prolate, provokes an imaginary part to the complex viscosity, η ω. In other words, pleomorphism introduces elasticity to coronavirus suspensions. From Fig. 8, we also learn that when prolate coronavirus is compared to its oblate counterpart (where one ellipticity is the reciprocal of the other), the prolate decreases η ω more than the oblate. Further, from Fig. 9, we learn that when prolate coronavirus is compared to its oblate counterpart (where one ellipticity is the reciprocal of the other), the oblate increases η ω less than the prolate. Finally, Fig. 8 shows the order of descent for η ω to be

c a = 1 , 4 5 , 5 4 , 2 3 , 3 2 ,
(30)

and Fig. 9 shows the order of ascent for η ω to be

c a = 3 2 , 2 3 , 5 4 , 4 5 , 1 ,
(31)

which is the reverse of Eq. (30).

TABLE III.

Ellipsoidal coronavirus particle characteristics from general rigid bead-rod theory with Nc=256 and Np=74.

Coronavirus ellipsoid, c/aI1mL2I2mL2I3mL2abν2baνη0ηsnkTλλλ0λ0DrΨ1,0λη0ηs
23  1.06×106  1.06×106  1.46×106  5.66×105  8.70×102  5.66×106  5.43×102  1.69  2.12×106  7.86×108  2.10×101  
45  0.86×106  0.86×106  1.03×106  4.45×105  2.79×102  7.10×106  1.77×102  1.61  1.69×106  9.86×108  6.94×102  
32  4.98×105  4.98×105  3.15×105  2.05×105  8.11×102  1.20×105  6.56×102  1.32  9.97×105  1.67×107  2.46×101  
54  5.58×105  5.58×105  4.41×105  2.55×105  2.67×102  1.07×105  1.95×102  1.39  1.12×106  1.49×107  7.67×102  
1  2.48×102  2.48×102  2.48×102  1.24×102  1.19×108  2.24×102  7.96×109  1.5  4.96×102  3.36×104  1.60×108  
Coronavirus ellipsoid, c/aI1mL2I2mL2I3mL2abν2baνη0ηsnkTλλλ0λ0DrΨ1,0λη0ηs
23  1.06×106  1.06×106  1.46×106  5.66×105  8.70×102  5.66×106  5.43×102  1.69  2.12×106  7.86×108  2.10×101  
45  0.86×106  0.86×106  1.03×106  4.45×105  2.79×102  7.10×106  1.77×102  1.61  1.69×106  9.86×108  6.94×102  
32  4.98×105  4.98×105  3.15×105  2.05×105  8.11×102  1.20×105  6.56×102  1.32  9.97×105  1.67×107  2.46×101  
54  5.58×105  5.58×105  4.41×105  2.55×105  2.67×102  1.07×105  1.95×102  1.39  1.12×106  1.49×107  7.67×102  
1  2.48×102  2.48×102  2.48×102  1.24×102  1.19×108  2.24×102  7.96×109  1.5  4.96×102  3.36×104  1.60×108  
FIG. 7.

I 3 I 1 plane locating bead-rod models of pleomorphic coronavirus (from columns 2 and 3 of Table III), with the oblate ones lying above the diagonal ( I 3 > I 1), the prolate, below ( I 3 < I 1), and the spherical, near the origin ( I 3 = I 1).

FIG. 7.

I 3 I 1 plane locating bead-rod models of pleomorphic coronavirus (from columns 2 and 3 of Table III), with the oblate ones lying above the diagonal ( I 3 > I 1), the prolate, below ( I 3 < I 1), and the spherical, near the origin ( I 3 = I 1).

Close modal
FIG. 8.

Effect of ellipticity on the dimensionless real part of the complex viscosity, η η s / η 0 η s, using Eq. (9), spherical c / a = 1 (black), prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green) ( N c = 256, N p = 74, respectively, rows 1–4 of Table III).

FIG. 8.

Effect of ellipticity on the dimensionless real part of the complex viscosity, η η s / η 0 η s, using Eq. (9), spherical c / a = 1 (black), prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green) ( N c = 256, N p = 74, respectively, rows 1–4 of Table III).

Close modal
FIG. 9.

Effect of ellipticity on the dimensionless minus the imaginary part of the complex viscosity, η / η 0 η s, using Eq. (10), prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green) ( N c = 256, N p = 74, respectively, rows 1–4 of Table III).

FIG. 9.

Effect of ellipticity on the dimensionless minus the imaginary part of the complex viscosity, η / η 0 η s, using Eq. (10), prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green) ( N c = 256, N p = 74, respectively, rows 1–4 of Table III).

Close modal

Figure 11 shows that pleomorphism, be it oblate or prolate, causes the dimensionless rotational diffusivity to decrease. Further, from Fig. 11 we learn that when prolate coronavirus is compared to its oblate counterpart (where one ellipticity is the reciprocal of the other), the oblate decreases the dimensionless rotational diffusivity more than the prolate. Figure 10 recalls the canonical dimensionless rotational diffusivity behavior of spherical coronavirus particles (Fig. 12 of Ref. 12 and Fig. 5 of Ref. 13). Comparing Fig. 10 to Fig. 11, we discover that with the reported pleomorphisms [Eq. (15)], λ o D r lands about 3 orders of magnitude below the canonical rotational diffusivity of spherical coronavirus. Finally, mindful of Fig. 11, Fig. 10 shows the order of descent for λ o D r to be

c a = 1 , 3 2 , 5 4 , 4 5 , 2 3 ,
(32)

which differs from both Eqs. (30) and (31), and is not monotonic.

FIG. 10.

Dimensionless rotational diffusivity λ 0 D r from Eq. (6)vs peplomer population N p ( N c = 256). Spherical c / a = 1.

FIG. 10.

Dimensionless rotational diffusivity λ 0 D r from Eq. (6)vs peplomer population N p ( N c = 256). Spherical c / a = 1.

Close modal
FIG. 11.

Dimensionless rotational diffusivity λ 0 D r from Eq. (6) vs peplomer population N p ( N c = 256). Prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green).

FIG. 11.

Dimensionless rotational diffusivity λ 0 D r from Eq. (6) vs peplomer population N p ( N c = 256). Prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green).

Close modal

From column 8 of Table III and Fig. 12, we find that (i) oblate coronavirus pleomorphism increases the dimensionless zero-shear viscosity, η 0, and (ii) prolate coronavirus pleomorphism decreases the dimensionless zero-shear viscosity. We also find, from column 11 of Table III and Fig. 12, that coronavirus pleomorphism, be it oblate or prolate, increases the zero-shear first normal stress coefficient, Ψ 1 , 0.

Through the lens of general rigid bead-rod theory, we have explored the role coronavirus cross-sectional ellipticity on its rotational diffusivity, the transport property around which its cell attachment revolves. We find that coronavirus pleomorphism, be it oblate or prolate, injects into the complex viscosity (i) non-Newtonian behavior, and (ii) elasticity (Figs. 8 and 9, respectively). We also learn that when prolate coronavirus is compared to its oblate counterpart, (i) the prolate decreases η ω more than the oblate, and (ii) the oblate increases η ω less than the prolate. We find that (i) oblate coronavirus pleomorphism increases the dimensionless zero-shear viscosity, η 0, and (ii) prolate coronavirus pleomorphism decreases the dimensionless zero-shear viscosity (column 8 of Table III). We also find that coronavirus pleomorphism, be it oblate or prolate, increases the zero-shear first normal stress coefficient, Ψ 1 , 0 (column 11 of Table III). We further find that rotational diffusivity decreases with coronavirus pleomorphism, be it oblate or prolate, and specifically that this descent is monotonic with the ratio of major to minor axis lengths, c / a (see Sec. IV A).

In this paper, we explored coronavirus pleomorphism using single beads to approximate the hydrodynamic resistance contributed by its spikes. However, from previous work, we know that the bulbous triangular equidimensional shape of the coronavirus spikes reduces rotational diffusivity (see Figs. 2, 4, and 5 of Ref. 13). We leave the important combination of coronavirus pleomorphism with spike triangularity for another day.

FIG. 12.

Dimensionless zero-shear viscosity, η 0, and zero-shear first normal stress difference coefficient, Ψ 1 , 0 from Eqs. (11) and (13), respectively, vs ellipticity. Spherical c / a = 1 (black), Prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green).

FIG. 12.

Dimensionless zero-shear viscosity, η 0, and zero-shear first normal stress difference coefficient, Ψ 1 , 0 from Eqs. (11) and (13), respectively, vs ellipticity. Spherical c / a = 1 (black), Prolate c / a = 3 / 2 (blue), prolate c / a = 5 / 4 (red), oblate c / a = 2 / 3 (magenta), and oblate c / a = 4 / 5 (green).

Close modal

Our previous work on the rotational diffusivity of coronavirus suspensions has excluded interferences of Stokes flow velocity fields between nearby spikes.12,13 We did so because no theory for incorporating hydrodynamic interactions in general rigid bead rod theory had been advanced. However, such a method for incorporating hydrodynamic interactions analytically has now been developed (Sec. III. of Ref. 27) and used (Sec. V. of Refs. 27 and 28). For the spherical capsid, this problem has also recently been attacked by means of molecular dynamics simulations.29 We leave the exploration of how spike hydrodynamic interactions affect coronavirus rotational diffusivity for future work.

General rigid bead-rod theory can only be applied to suspensions for whom the suspended particle structure is known. Thus, general rigid bead-rod theory cannot be more accurate than our knowledge of this structure. Of coronaviruses, Bàrcena et al. wisely observed: “As a result of their pleomorphic nature, our structural insight into the coronavirion is still rudimentary, and it is based mainly on 2D electron microscopy.” Table X of Ref. 12 summarizes what we know about the average values of the physical characteristics of coronavirus particles, and specifically, how we arrived at N p = 74 to rely upon in this study. From this study, we can see that a distribution of particle ellipticities will result in a corresponding distribution of particle rotational diffusivities. The virus likelihood of attaching, we would expect, will be determined by this distribution. To handle such mixtures of different species, we can rewrite Eqs. (9) and (10) for a dispersed system (see problem14C.2 of Ref. 9, see also Sec. 26 of Ref. 30). We leave this for another day.

Though coronavirus pleomorphism is well-known, little is understood about its causes. Peplomer repulsions, for instance, might themselves deform the capsid. So might flow. Coronavirus pleomorphism may also reflect different asymmetrical arrangements of its capsid contents.31 Our work is silent on these important considerations.

Pleomorphism causes the peplomer spacing to differ over the oblate or prolate coronavirus surfaces. In as much as attachment requires the alignment of an adjacent peplomer pair with its nominally rectangular dimeric receptor, this pleomorphic peplomer spacing will complicate Eq. (2) of Ref. 12 for the attachment probability. In other words, the probability of finding a peplomer both (i) aligned with said receptor and (ii) matching the dimeric spacing of said receptor differs from place to place over the surface of a pleomorphic coronavirus.

The uninitiated may wonder why our bead-rod model employs an empty capsid. After all, the capsid contents are what causes disease. However, only the parts of the suspended particle resisted hydrodynamically by the surrounding Newtonian fluid can contribute to the resisted rotation and to the complex viscosity of the suspension. Thus, the rotational diffusivity of an empty capsid does not differ from a packed one. Ιn general rigid bead-rod theory, only beads dragged through solvent can contribute to the rheology. Pieces of the macromolecule that are shielded from the surrounding solvent, such as capsid contents, are not to be counted in the general bead-rod model.

In this paper, we just considered axisymmetric pleomorphism, both oblate and prolate. We know of no characterization, microscopic or otherwise, detailed enough to distinguish oblate from prolate, or for that matter, asymmetric from axisymmetric. Moreover, we still know of no rheological characterization of the coronavirus suspension. We consider these experimental measurements to be crucial next steps for those equipped to handle live coronavirus.

Our general rigid bead-rod theory relies entirely on orientation to explain the polymer contribution to the viscosity and elasticity of polymeric liquids. Our theory is not to be confused with the competing theory, for suspensions of ellipsoids of revolution, which yields, for the real and minus the imaginary parts of the complex viscosity [after Eq. (11) of Ref. 34)

η η s η 0 η s = 5 2 + 26 147 ε 2 + 3 5 ε 2 1 + λ 2 ω 2 5 2 + 571 735 ε 2 ,
(33)
η η 0 η s = 3 5 ε 2 λ ω 1 + λ 2 ω 2 5 2 + 571 735 ε 2 ,
(34)

where

η 0 η s = η s Φ 5 2 + 571 735 ε 2 ,
(35)

and where

ε c a 1 ,
(36)

where c / a is the capsid aspect ratio. Equations (33)–(36) are thus for suspensions of otherwise featureless ellipsoids (spikeless), be they prolate or oblate, arrived at elegantly through the competing ellipsoid suspension orientation theory. Unlike general rigid bead-rod theory, however, ellipsoid suspension orientation is silent on the relation between λ and the shapes and sizes of the suspended objects.

We have restricted this work to small-amplitude, namely, where Eq. (7) obtains. For large-amplitude oscillatory shear flow, where Eq. (7) does not obtain, we bridge to the corotational Jeffreys fluid, following the method of Sec. X of Ref. 11. We leave this intriguing problem to future work.

Whereas much prior work on fluid physics related to the virus has attacked transmission,32,33 this paper focuses on the ab initio calculation of coronavirus transport properties. Specifically, we have determined the rotational diffusivity, the property governing the particle alignment for cell attachment (see Sec. I of Ref. 12). Although our work is mainly curiosity driven, it may deepen our understanding of drug, vaccine, and cellular infection mechanisms.

One might prefer to add other comparisons with experiment, and specifically with relevant complex viscosity vs coronavirus observations. However, we know of no experimental measurements on coronavirus suspensions. In other words, to our knowledge, the transport properties explored in this paper have yet to be explored in the laboratory.

Under the microscope, we see agglomeration of coronavirus particles, mechanically interlocked by interdigitation of the bulbous spikes [see Fig. 1(d) of Ref. 35]. The charge repulsion of the interdigitating spikes will, of course, rearrange the spikes. By rearrange, we mean spreading near the interdigitation, and crowding elsewhere. Figure 13 (Multimedia view) shows our general rigid bead-rod model of an interdigitated cluster of coronaviruses. We leave the energy minimization for interdigitated clusters, and the calculation of the coronavirus cluster rotational diffusivity for another day.

FIG. 13.

Cryo-electron tomography of pleomorphic septocluster of interdigitated coronavirus (left) [Left panel of Fig. 1(d) of Ref. 35] vs our general rigid bead-rod pleomorphic model of same (right). Left panel reproduced with permission from B.W. Neuman, G. Kiss, G., A.H. Kunding, D. Bhella, M.F. Baksh, S. Connelly, B. Droese, J.P. Klaus, S. Makino, S.G. Sawicki, S.G. Siddell, “A structural analysis of M protein in coronavirus assembly and morphology,” J. Struct. Biol. 174(1), 11–22 (2011). Copyright 2011 Elsevier. Multimedia view: https://doi.org/10.1063/5.0094771.6

FIG. 13.

Cryo-electron tomography of pleomorphic septocluster of interdigitated coronavirus (left) [Left panel of Fig. 1(d) of Ref. 35] vs our general rigid bead-rod pleomorphic model of same (right). Left panel reproduced with permission from B.W. Neuman, G. Kiss, G., A.H. Kunding, D. Bhella, M.F. Baksh, S. Connelly, B. Droese, J.P. Klaus, S. Makino, S.G. Sawicki, S.G. Siddell, “A structural analysis of M protein in coronavirus assembly and morphology,” J. Struct. Biol. 174(1), 11–22 (2011). Copyright 2011 Elsevier. Multimedia view: https://doi.org/10.1063/5.0094771.6

Close modal

This research was undertaken, in part, thanks to support from the Canada Research Chairs program of the Government of Canada for the Natural Sciences and Engineering Research Council of Canada (NSERC) Tier 1 Canada Research Chair in Rheology. This research was also undertaken, in part, thanks to support from the Discovery Grant program of the Natural Sciences and Engineering Research Council of Canada (NSERC) (A. J. Giacomin), Vanier Canada Graduate Scholarship (M. A. Kanso), and the Mitacs Research Training Award (A. J. Giacomin and M. A. Kanso). A. J. Giacomin is indebted to the Faculty of Applied Science and Engineering of Queen's University at Kingston for its support through a Research Initiation Grant (RIG). V. Chaurasia and E. Fried gratefully acknowledge support from the Okinawa Institute of Science and Technology Graduate University with subsidy funding from the Cabinet Office, Government of Japan.

The authors have no conflicts to disclose.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
B. W.
Neuman
,
B. D.
Adair
,
C.
Yoshioka
,
J. D.
Quispe
,
G.
Orca
,
P.
Kuhn
,
R. A.
Milligan
,
M.
Yeager
, and
M. J.
Buchmeier
, “
Supramolecular architecture of severe acute respiratory syndrome coronavirus revealed by electron cryomicroscopy
,”
J. Virol.
80
(
16
),
7918
(
2006
).
2.
N.
Zhu
,
D.
Zhang
,
W.
Wang
,
X.
Li
,
B.
Yang
,
J.
Song
,
X.
Zhao
,
B.
Huang
,
W.
Shi
,
R.
Lu
,
P.
Niu
,
F.
Zhan
 et al, “
A novel coronavirus from patients with pneumonia In China, 2019
,”
N Engl. J. Med.
382
(
8
),
727
733
(
2020
).
3.
M. A.
Kanso
, “
Polymeric liquid behavior in oscillatory shear flow
,” M.S. thesis (
Queen's University
,
Kingston, Canada
,
2019
).
4.
M. A.
Kanso
and
A. J.
Giacomin
, “
van Gurp-Palmen relations for long-chain branching from general rigid bead-rod theory
,”
Phys. Fluids
32
(
3
),
033101
(
2020
); Erratum: In Eq. (44) “ I 3 > I 1” should be “ I 3 < I 1.” Editor's pick.
5.
M. A.
Kanso
,
A. J.
Giacomin
, and
C.
Saengow
, “
Large-amplitude oscillatory shear flow loops for long-chain branching from general rigid bead-rod theory
,”
Phys. Fluids
32
(
5
),
053102
(
2020
).
6.
M. A.
Kanso
,
A. J.
Giacomin
,
C.
Saengow
, and
J. H.
Piette
, “
Diblock copolymer architecture and complex viscosity
,”
Int. J. Mod. Phys. B
34
,
2040110
(
2020
).
7.
D.
Singhal
,
M. A.
Kanso
,
S. J.
Coombs
, and
A. J.
Giacomin
, “
Complex viscosity of poly[n]catenanes and olympiadanes
,”
Phys. Fluids
34
(
3
),
033112
(
2022
).
8.
O.
Hassager
, “
Kinetic theory and rheology of bead-rod models for macromolecular solutions. II. Linear unsteady flow properties
,”
J. Chem. Phys.
60
(
10
),
4001
4008
(
1974
); Erratum: in Eq. (2) of “ 1 / 2” should be “ 1 / 2” and “ ” should be “ .”
9.
R. B.
Bird
,
C. F.
Curtiss
,
R. C.
Armstrong
, and
O.
Hassager
,
Dynamics of Polymeric Liquids
2nd ed. (
John Wiley & Sons, Inc
.,
New York
,
1987
), Vol.
2
; Errata: On p. 409 of the first printing, the n + m ! in the denominator should be n m !; In Table 16.4–1, under L entry “length of rod” should be “bead center to center length of a rigid dumbbell”; In Fig. 14.1–2 caption, “Multibead rods of length L” should be “Multibead rods of length L + d.”
10.
R. B.
Bird
,
O.
Hassager
,
R. C.
Armstrong
, and
C. F.
Curtiss
,
Dynamics of Polymeric Liquids
, 1st ed. (
John Wiley and Sons, Inc
.,
New York
,
1977
), Vol.
2
.
11.
M. A.
Kanso
,
A. J.
Giacomin
,
C.
Saengow
, and
J. H.
Piette
, “
Macromolecular architecture and complex viscosity
,”
Phys. Fluids
31
(
8
),
087107
(
2019
); Editor's pick. Errata: Ganged in Ref. 8 of Ref. 13.
12.
M. A.
Kanso
,
J. H.
Piette
,
J. A.
Hanna
, and
A. J.
Giacomin
, “
Coronavirus rotational diffusivity
,”
Phys. Fluids
32
(
11
),
113101
(
2020
). Feature article. Cover article.
13.
M. A.
Kanso
,
V.
Chaurasia
,
E.
Fried
, and
A. J.
Giacomin
, “
Peplomer bulb shape and coronavirus rotational diffusivity
,”
Phys. Fluids
33
(
3
),
033115
(
2021
).
14.
A. J.
Giacomin
and
M. A.
Kanso
, “
General rigid bead-rod macromolecular theory
,” in
Recent Advances in Rheology: Theory, Biorheology, Suspension and Interfacial Rheology
, edited by
D.
De Kee
and
A.
Ramachandran
(
AIP Publishing
,
Melville
,
2022
), Chap. II, pp. 2-1–2-32.
15.
W. E.
Stewart
and
J. P.
Sørensen
, “
Hydrodynamic interaction effects in rigid dumbbell suspensions. II. Computations for steady shear flow
,”
Trans. Soc. Rheol.
16
(
1
),
1
13
(
1972
).
16.
J. H.
Piette
,
L. M.
Jbara
,
C.
Saengow
, and
A. J.
Giacomin
, “
Exact coefficients for rigid dumbbell suspensions for steady shear flow material function expansions
,”
Phys. Fluids
31
(
2
),
021212
(
2019
); Erratum: Above Eq. (83), “one other” should be “one other use.”
17.
R. B.
Bird
and
A. J.
Giacomin
, “
Who conceived the complex viscosity?
,”
Rheol. Acta
51
(
6
),
481
486
(
2012
).
18.
A. J.
Giacomin
and
R. B.
Bird
, “
Erratum: Official nomenclature of The Society of Rheology
,”
J. Rheol.
55
(
4
),
921
923
(
2011
).
19.
J. D.
Ferry
,
Viscoelastic Properties of Polymers
, 3rd ed. (
Wiley
,
New York
,
1980
).
20.
R. B.
Bird
,
R. C.
Armstrong
, and
O.
Hassager
,
Dynamics of Polymeric Liquids
, 1st ed. (
Wiley
,
New York
,
1977
), Vol.
1
.
21.
V.
Chaurasia
,
Y.-C.
Chen
, and
E.
Fried
, “
Interacting charged elastic loops on a sphere
,”
J. Mech. Phys. Solids
134
,
103771
(
2020
).
22.
V.
Chaurasia
, “
Variational formulation of charged curves confined to a sphere
,” Ph.D. thesis (
University of Houston
,
Houston
,
2018
).
23.
J.
Nocedal
and
M. L.
Overton
, “
Projected Hessian updating algorithms for nonlinearly constrained optimization
,”
SIAM J. Numer. Anal.
22
(
5
),
821
850
(
1985
).
24.
J. J.
Thomson
, “
XXIV. On the structure of the atom: An investigation of the stability and periods of oscillation of a number of corpuscles arranged at equal intervals around the circumference of a circle; with application of the results to the theory of atomic structure
,”
London, Edinburgh, Dublin Philos. Mag. J. Sci.
7
(
39
),
237
265
(
1904
).
25.
D. J.
Wales
and
S.
Ulker
, “
Structure and dynamics of spherical crystals characterized for the Thomson problem
,”
Phys. Rev. B
74
(
21
),
212101
(
2006
).
26.
D. J.
Wales
,
H.
McKay
, and
E. L.
Altschuler
, “
Defect motifs for spherical topologies
,”
Phys. Rev. B
79
(
22
),
224115
(
2009
).
27.
M. C.
Pak
,
K.-I.
Kim
,
M. A.
Kanso
, and
A. J.
Giacomin
, “
General rigid bead-rod theory with hydrodynamic interaction for polymer viscoelasticity
,”
Phys. Fluids
34
(
2
),
023106
(
2022
).
28.
M. A.
Kanso
,
M. C.
Pak
,
K.-I.
Kim
,
S. J.
Coombs
, and
A. J.
Giacomin
, “
Hydrodynamic interaction and complex viscosity of multi-bead rods
,”
Phys. Fluids
34
(
4
),
043102
(
2022
). Editor's pick.
29.
N.
Moreno
,
D.
Moreno-Chaparro
,
F. B.
Usabiaga
, and
M.
Ellero
, “
Hydrodynamics of spike proteins dictate a transport-affinity competition for SARS-CoV-2 and other enveloped viruses
,” bioRxiv (
2022
).
30.
R. B.
Bird
,
H. R.
Warner
, and
D. C.
Evans
, “
Kinetic theory and rheology of dumbbell suspensions with Brownian motion
,”
Fortschr. Hochpolym.-Forsch.
8
,
1
90
(
1971
).
31.
P. J.
de Pablo
and
I. A. T.
Schaap
, “
Atomic force microscopy of viruses
,” in
Physical Virology
, edited by
U. F.
Greber
(
Springer
,
Cham
,
2019
) Chap. VIII, pp.
159
179
.
32.
T.
Dbouk
and
D.
Drikakis
, “
On coughing and airborne droplet transmission to humans
,”
Phys. Fluids
32
(
5
),
053310
(
2020
).
33.
H.
Wang
,
Z.
Li
,
X.
Zhang
,
L.
Zhu
,
Y.
Liu
, and
S.
Wang
, “
The motion of respiratory droplets produced by coughing
,”
Phys. Fluids
32
(
12
),
125102
(
2020
).
34.
L. G.
Leal
and
E. J.
Hinch
, “
The rheology of a suspension of nearly spherical particles subject to Brownian rotations
,”
J. Fluid Mech.
55
(
4
),
745
765
(
1972
).
35.
B. W.
Neuman
,
G.
Kiss, G.
,
A. H.
Kunding
,
D.
Bhella
,
M. F.
Baksh
,
S.
Connelly
,
B.
Droese
,
J. P.
Klaus
,
S.
Makino
,
S. G.
Sawicki
, and
S. G.
Siddell
, “
A structural analysis of M protein in coronavirus assembly and morphology
,”
J. Struct. Biol.
174
(
1
),
11
22
(
2011
).