key: cord-0078565-9h4aatjh authors: Mahata, Paritosh; Vennamneni, Laxminarsimharao; Chattopadhyay, Soham title: A mechanical-thermodynamic model for understanding endocytosis of COVID-19 virus SARS-CoV-2 date: 2022-05-17 journal: Proc Inst Mech Eng C J Mech Eng Sci DOI: 10.1177/09544062221098538 sha: 6299e9861c0f17158cf415ad88dd96e451647739 doc_id: 78565 cord_uid: 9h4aatjh We analyze the endocytosis process of COVID-19 (coronavirus disease 2019) virus SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) using a mechanical-thermodynamic model. The virus particle is designed to interface with the cell membrane as a hard sphere. The role of cytoplasmic BAR (Bin/Amphiphysin/RVs) proteins is considered in the endocytosis. Interestingly, the Endophilin N-BAR cytoplasmic proteins show resistance in participating endocytosis, whereas F-BAR, Arfaptin BAR, Amphiphysin N-BAR, and PX-BAR proteins participate in endocytosis. The increase in membrane tension, concentrated force between the cell membrane receptor, and spike glycoprotein present on the surface of virus particle promote the endocytosis. Also, the increase in the bending modulus of membrane leads to the two-phase solution of BAR protein concentration on the interior of cell membrane surface. We observe an unstable region of protein concentration, which may help one to retard the endocytosis process and thus the viral infection. Though the present study is focused on SARS-CoV-2, it can be extended to understand any other viral infections, involving endocytosis process. The effect of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) or COVID-19 (coronavirus disease 2019) on humankind is one of the pandemics with high infection rate facing the globe. The genome of SARS-CoV-2 is similar to the SARS-CoV. [1] [2] [3] The COVID-19 mainly causes respiratory problems while damaging cells in the respiratory system. The SARS-CoV-2 has a spherical structure with spike proteins on its surface. When a SARS-CoV-2 particle approaches cell, the spike proteins will be attached to the receptors present on cell surface. 3 It enters into the cell through the endocytosis and inside cell it will make few more copies of SARS-CoV-2 and they will come out while damaging cell through exocytosis process. [4] [5] [6] Through these endocytosis and exocytosis processes, they damage cells in the respiratory system and leads to respiratory problems. In the following, we explain qualitatively the endocytosis process of SARS-CoV-2. Endocytosis is the mechanism by which micro/nano elements can enter inside a cell to keep its biological functions running. 7, 8 The physical mechanism depicting the endocytosis process of CoV-2 is represented with the help of a schematic shown in Figure 1 . As the CoV-2 virus approaches to the cell membrane, the spike glycoproteins on the CoV-2 are attracted by the receptors present on the cell membrane's outer surface (see Figure 1 (a)). When a glycoprotein binds with the receptor, the cell membrane deforms locally. It helps to attain the curvature of the membrane equal to the CoV-2 as shown in Figure 1 (a). The modified curvature of the cell membrane is sensed by BAR proteins (like F-BAR, PX-BAR, Amphiphysin and Endophilin N-BAR [9] [10] [11] ]) present in the cytoplasm (an aqueous solution inside a cell containing salts and proteins) and they bind on the membrane's interior surface. Other cytoplasmic proteins like clathrin and caveolae may also sense and generate the membrane curvature in endocytosis. 12, 13 In this work, we analyze the role of BAR proteins as they help in endocytosis by forming vesicle and tubule of the cell membrane. 9, 10, [14] [15] [16] The membrane curvature is increased by the local crowding of BAR proteins on the membrane's interior surface. The increase in membrane curvature increases the accessibility of spikes to the receptors and further increases the membrane curvature as shown in Figure 1 (b). It helps the virus to enter into the cell gradually. Finally, the combined efforts made by the receptors on the membrane's outer surface and cytoplasmic BAR proteins in the cell interior help the CoV-2 to enter into the cell as a coated vesicle as shown in Figure 1 (c). In the earlier works of Inoue et al., 17 Burkard et al., 18 and Wang et al., 19 they have done experiments and reported the endocytosis of SARS-CoV. They discovered that SARS-CoV enters target cells via clathrin-mediated endocytosis. However, their findings may not accurately reflect the endocytosis of the new coronavirus SARS-CoV-2, which caused the COVID-19 pandemic. Shang et al. 20 used biochemical and pseudovirus entry tests to study the receptor binding and protease activation of SARS-CoV-2 spike glycoprotein. They discover that SARS-CoV-2 has more affinity for binding to the hACE2 protein than the SARS-CoV. This is one of the most important mechanisms for SARS-CoV-2 endocytosis that scientists have discovered so far. Glebov 21 analyzed the process of SARS-CoV-2 endocytosis in a recent paper and identified that the clinically authorized medications as the prospective candidates for repurposing as blockers of various potential routes. Yang and Shen 22 focus their review on the role of the endocytic pathway and autophagy process in viral infection, as well as the development of therapeutic medicines that target these processes. All these reports are available on probable mechanism of infection of CoV-2 in human body and possible mode of action for its prevention. Mukherjee and Mahata 23 have recently explored the endocytosis of SARS-CoV-2 using a finite element model that they proposed. The virus is modeled as hard sphere, while the cell membrane is modeled as an anisotropic elastic material. The role of BAR proteins is also examined in their study using a linear curvature-coupling model with the BAR protein concentration. But the appropriate mathematical model is missing in the literature for understanding SARS-CoV-2 endocytosis. Although, several other studies have been performed for understanding endocytosis behavior of nanoparticles. For example, Vácha et al. 24 have studied the receptor-mediated endocytosis of nanoparticles of different shapes with the help of a molecular dynamics (MD) simulation. They have found that the efficiency of endocytosis is higher for spherocylindrical particles. Later on the effects of orientations and rotations of nanoparticles are considered by Dasgupta et al. 25 and Tang et al. 26 It is found that the endocytosis follows different wrapping pathway due to the competition between stretching and bending of membrane. 26 This motivates the rotation of the nanoparticles during endocytosis. The high-aspect ratio or slender particle with round tips, participate in the endocytosis process while orienting the long axis of particle parallel to membrane, whereas the small aspect ratio particles with flat tips enter while orienting the long axis of particle perpendicular to membrane surface. 25 The effect of nanoparticle geometry, stiffness, and biomechanical properties of membrane on the endocytosis process discussed in refs. 27, 28. The endocytosis process while considering the adhesion between the nanoparticle and membrane has been studied by Agudo-Canalejio and Lipowsky. 29, 30 As the SARS-CoV-2 enters into human body, different types of serum proteins have tendency to settle on the surface of SARS-CoV-2 called protein corona. This might affect the infection rate of virus particle. 31 In the latest work of Yin et al., 32 using all-atom molecular dynamic simulations and molecular docking, they analyzed the effect of corona proteins on the infection rate of virus particle. In this work, we analyze the endocytosis of SARS-CoV-2 using a simple mechanical-thermodynamic model. Here, using a proposed mathematical model, we analyze the effect of cytoplasmic BAR protein concentration on the endocytosis. The virus is modelled as hard sphere. The cell membrane deformation due to the insertion of virus is modeled by Helfrich model. 33, 34 The adsorbed phase of proteins on interior surface of the membrane is modeled using Mahata and Das. 35 For different virus particle's depth of insertion, the concentrations of BAR proteins on the membrane surface are evaluated from equilibrium configuration of membrane cap (deformed region on the membrane surface adjacent to virus particle) generated by the virus during endocytosis. Here, we analyze parameters like bending modulus and tension of cell membrane on protein concentration. The paper is organized as follows. In Mathematical model, a mathematical model for SARS-CoV-2 endocytosis is established. Non-dimensionalization discusses the normalized form of the system's total free energy. Equilibrium configuration of spherical membrane cap generated by the virus during endocytosis is analyzed in Equilibrium of membrane cap. We discuss different model parameters in Model parameters. The results of the equilibrium configuration of the membrane cap are discussed in Results and discussion. We conclude in Conclusions. We propose a mathematical model to study the endocytosis behavior of SARS-CoV-2 as shown in Figure 2 . The deformation of virus particle is negligible as compared to the deformation of cell membrane and also the molecular structure of SARS-CoV-2 retains its spherical shape 36 during the endocytosis process. Thus, in this work we treat the virus particle as a rigid spherical particle. In general, the BAR protein will undergo small deformation during the endocytosis process. But, in order to understand the membrane deformation, we use a simplified model while assuming both BAR protein and viral particles to be rigid and analyze the endocytosis process of virus particle. Due to the applied force f, the particle as well as the cell membrane are displaced by h and it follows the curvature of spherical particle as shown in Figure 2 . In recent works of Dasgupta et al. 25 and Yi and Gao, 28 the specific receptorligand interaction has been considered to generate initial curvature on the cell membrane, so that the intrinsic curvature of the cytoplasmic BAR protein sense the membrane curvature and bind on it. In the current work, we use an equivalent simplified model 34, 37, 38, 39, 40, 41 with a concentrated force f applied on the virus particle to initiate the nonspecific adsorption of cytoplasmic BAR proteins by generating the cell membrane curvature. The total free energy of the system can be expressed as where F b is the bending energy of the cell membrane, F s is the stretching energy due to lateral tension in the membrane, F m is the mixing free energy due to adsorption of cytoplasmic proteins on the interior of the membrane surface, and f h represents the work done by virus particle, which will be equivalent to the specific interaction energy. In the following, we determine the contributions coming from all these free energies. The cell membrane behaves like a fluid for in-plane shear, but it resists the out-of-plane bending. [42] [43] [44] [45] [46] [47] [48] Due to this nature and small thickness of the membrane, one can consider curvature dependent bending energy obtained from Helfrich model. 33, 34 The bending energy is expressed as where κ is the bending modulus of the membrane, A m is the surface area of membrane cap generated by the virus particle during endocytosis, H m is the mean curvature of the membrane cap, and C s is the spontaneous curvature of the membrane, respectively. For a membrane cap of height h and radius R, A m = 2πRh and H m = 1/R. We assume that the virus generates membrane cap of shape similar to the virus particle during the endocytosis as the cytoplasmic proteins adsorbed on the membrane surface. The protein concentration is assumed to be linear in the spontaneous curvature of the membrane cap, 49,37,39 so that C s = C p θ, where C p is the intrinsic curvature of BAR protein, and θ is adsorbed BAR protein concentration on the membrane cap surface. It is assumed that the intrinsic curvature of the protein does not depend on the curvature of membrane cap. In general, the cell membrane is not much stretchable. Once the lateral tension exceeds few units, the cell membrane suddenly fails without much premonition. 50 We consider the situation that the membrane cap generated by the virus particle is coupled to a membrane area reservoir at constant membrane tension. This is typical for the biological cells. 51, 52 and mimicked in a situation, when a lipid membrane tether is pulled by the use of micropipette aspiration from a lipid membrane reservoir. 34, 37 If σ be the constant stretching energy per unit area due to membrane tension, then the total stretching energy can be written as The contribution of mixing free energy in equation (1) comes due to the adsorption of cytoplasmic proteins onto the surface of membrane cap generated by the virus during endocytosis. It includes the interaction free energy due to protein-protein interactions and the pressure energy due to binding of proteins on membrane surface. Before starting endocytosis of the virus particle, the chemical potential of the proteins in the cytoplasm and binding on cell membrane's inside surface should be equal in equilibrium. Also, we can assume that the cell membrane is locally flat in this configuration. When a virus particle tries to enter into cell, the cytoplasmic proteins bounded on the inside surface of the flat membrane will also come onto the surface of the membrane cap. We should neglect the effects of all these proteins as they are bounded to the cell membrane surface before starting the endocytosis. So, there are Legendre transformation 53 of mixing free energy from the flat membrane surface to the membrane cap surface. The Legendre transformed mixing free energy is expressed as 38,37 where f m (θ) is the free energy per unit area of the adsorbed proteins on the membrane cap surface due to mixing and protein-protein interactions, μ m is the chemical potential of adsorbed proteins on the flat membrane surface, b is the area of an individual protein molecule, and P m is the energy per unit area applied due to pressure of adsorbed proteins onto the membrane surface. As noted in Introduction, we take into account the binding of cytoplasmic BAR proteins to the membrane surface. The projected shape of a BAR protein can be represented by an elliptical particle approximately. Also, the soft interactions between the adsorbed protein molecules are expected on the membrane surface. To consider shape aspect of the protein molecules and soft interactions between them, we have chosen f m (θ) in equation (4) form Mahata and Das, 35 so that where k B is Boltzmann's constant, T is absolute temperature, a a is attraction parameter, a r repulsion parameter, and α represents shape parameter of the protein molecule. For elliptical particle, α > 1. For more details about different parameters in equation (5) can refer Mahata and Das. 35 Substituting equations (2), (3), and (4) into equation (1) with A m = 2πRh and H m = 1/R, we have where f m (θ) is obtained from equation (5), and μ m and P m are derived in the next section. In the following, we nondimensionalize equation (6). We use h (the displacement of virus), κ (bending modulus of membrane), and κ/h to non-dimensionalize length, energy and force, respectively. We obtain the following nondimensional parameters Using the above non-dimensional quantities in equation (6), we have In equation (8), the non-dimensional free energy per unit area due to protein-protein interactions where b ¼ κ h 2 , and a r ¼ κ k B T a r h 2 are nondimensional area, attraction, and repulsion parameters of protein molecules, respectively. We can rewrite the equation (6) as 2 θ 2 . Now, introducing a function G m ðθÞ ¼ A m g m ðθÞ, with A m is the dimensionless surface area of membrane cap, one can obtain and In equations (11) and (12) , N is the number of protein molecules binding on the membrane cap in its reference configuration and θ m is the corresponding protein concentration on the interior surface of membrane. In this analysis, we assume that the size of the cell is much larger as compared to membrane cap. So, we consider the reference configuration of the membrane cap to flat locally. This leads to the radius of curvature of membrane infinity. In deformed configuration, curvature of virus particle and membrane cap are same. Also, μ m is associated with the flat membrane surface. Then, using R = ∞ in equation (11), we have Now, with the help of equation (9), one can evaluate P m and μ m from equations (12) and (13), respectively. For the equilibrium configuration of membrane cap, the free energy of system should be minimum. To minimize free energy of the system, we differentiate the free energy with respect to R and θ and equate them to zero, that is ∂F ∂R ¼ 0 and ∂F ∂θ ¼ 0: This translates in to the following equations and Using equation (15) in equation (14), we can write We solve equation (16) for the protein concentration (θ), for the specified values of C p ,σ, and R. Here, we obtain f m ðθÞ and P m using equations (9) and (12), respectively. The stability criteria for the equilibrium configuration of the membrane cap are defined as 0 which leads us into We should choose the parameters a a , a r , and b in equation (9), so that they will satisfy the above inequality. The details about different parameters used in this analysis will be discussed in the following section. The approximate diameter of SARS-CoV-2 particle is 100 nm. 54 So, in this work, we use the radius of the spherical particle R = 50 nm. The radius of cytoplasmic BAR proteins varies between 10 and 100 nm. For example, the F-BAR proteins can bind with the membrane tubule of 25-100 nm radius. 55 The PX-BAR and Arfaptin BAR scaffold with membrane of radius of curvature 13 and 10 nm, respectively. 56 Also, the Amphiphysin and Endophilin N-BAR proteins bind with membrane having radius of curvature 11 and 15 nm, respectively. 57, 58 Since, R = 50 nm, we cannot expect that the radius of curvature of protein greater than 50 nm can bind with the membrane. So, the BAR proteins having radius of curvature in between 10 to 50 nm can bind with the membrane cap generated during the endocytosis of virus particle. The corresponding intrinsic curvature C p of the protein varies from 0.02 to 0.1 nm À1 . The typical values of cell membrane tension are extensively measured for different kind of cells. It can vary from 3 pN μm À1 for epithelial cells up to about 300 pN μm À1 for keratocytes. 51, 52 So, we can choose the range σ = 3 -300 μm À1 for the stretching energy per unit area due to membrane tension. Also, the typical values of bending modulus of cell membrane are taken in the range of 10 to 100 k B T. 59 As we pointed out in Mixing free energy, the projected shape of the above-mentioned BAR proteins can be approximated by an elliptical molecule as shown in Figure 3 . Considering interactions between the elliptical molecules, one can determine the mixing free energy per unit area f m (θ) as shown in equation (5). 35 For different BAR proteins, the axial ratio of the elliptical particles can be determined from the structural information of the proteins mentioned in Protein Data Bank (PDB). One can express the corresponding shape parameter where g is the axial ratio of the elliptical protein molecules and f represents the angular position of an arbitrary point located on the boundary of the elliptical molecule. 60 For the given major and minor axis values, one can obtain the individual area of protein molecule b. From the area of protein molecules, we can determine the Lennard-Jones potential parameters of protein molecules and subsequently the attraction and repulsive parameters a a and a r , respectively from Mahata and Das. 35 The interaction parameters as well as curvatures of different BAR proteins are mentioned in Table 1 . In this analysis, we choose the size of protein molecules at room temperature (300 K) and assume this temperature as the Boyle temperature for them. In Boyle temperature, second virial coefficient becomes zero 35, 61 and we obtain equal attraction and repulsive parameter for the protein molecules as mentioned in Table 1 . Using the parameter values in equation (16), we obtain the protein concentration θ on membrane cap for different values of θ m . Here, θ m represents the equilibrium protein concentration in the interior of the cell membrane surface before starting the endocytosis. The stability condition of the membrane cap is given by equation (18) . The results and corresponding discussions are mentioned in the following section. For different depth of insertion of virus particle, we obtain the equilibrium concentration θ of BAR proteins Figure 3 . A picture of a F-BAR protein is taken from Reference. 55. The projected shape of the protein can be represented by an elliptical particle. Other BAR proteins will have similar kind of structure. (cytoplasmic proteins) on the surface of membrane cap against different values of BAR protein concentration θ m in the reference configuration as shown in Figure 4 . We consider the presence of different BAR proteins separately in the cytoplasm during the endocytosis of the virus particle. The Endophilin N-BAR protein does not satisfy the stability criteria of the solution of the problem as mentioned in equation (18) . This suggests that the Endophilin N-BAR protein may not contribute in the endocytosis. We analyze the F-BAR, Arfaptin BAR, Amphiphysin N-BAR, and PX-BAR and they participate up to the half depth of insertion of the virus as shown in Figure 4(a) and (b) . Finally, three BAR proteins-F-BAR, PX-BAR, and Arfaptin BAR actively participate up to full depth of insertion of the virus as shown in Figure 4 (d). If θ < θ m , then it signifies that some of the cytoplasmic BAR proteins that have already attached on the interior surface of cell membrane released into the cytoplasm, before the endocytosis process. After a critical value of θ m , θ approaches to θ m . This suggests that no proteins release into the cytoplasm from the cell membrane surface during endocytosis. When θ>θ m , it is expected that extra proteins have come to bind on the cell membrane surface at the time of endocytosis. As we expect that the proteins F-BAR, PX-BAR, and Arfaptin BAR play important role in endocytosis of SARS-CoV-2, we analyze the effect of membrane tension on the concentration for these proteins as shown in Figure 5 . It is observed that the concentration increases with increase in membrane tension σ. One can expect this, because the increase in membrane tension stretches the membrane to occupy more proteins on its surface. For smaller values of θ m , we observe that θ is greater than θ m . This suggests that the additional proteins can bind on the membrane surface with increase in membrane tension. After a critical θ m , θ approaches θ m and cytoplasmic BAR proteins stop to bind on the membrane surface. This kind of behavior we observe for all the three BAR proteins considered in Figure 5 . The effect of membrane bending modulus κ on concentration (θ) for F-BAR, PX-BAR, and Arfaptin BAR shown in Figure 6 . We observe a region of θ, for which we do not get any stable solution of the membrane. This suggests that the proteins do not contribute in the endocytosis process. Here, we may expect two phase solution of the protein concentration on the cell membrane surface. In Figure 6 (a), for F-BAR protein with κ = 20k B T, we observe single phase solution. As the bending modulus of cell membrane increased, the single phase transforms into a Table 1 . Here, we use σ = 3 pN μm À1 and κ = 10k B T. two-phase solution. For Arfaptin and PX-BAR proteins, we always observe two-phase solutions for the range of κ values considered in Figure 6 (b) and (c). In Figures 4-6 , for smaller values of θ m , the equilibrium concentration of BAR proteins θ is absent. Also, for larger values of θ m , θ is invariant for different BAR proteins, different values of σ and k B T, respectively. If we can manipulate θ and θ m , such that they lie in the unstable region (as observed in Figures 4-6) , the virus can then avoid endocytosis. In Figure 7 , we calculate free energy F of the system from equation (10) against the concentration of F-BAR protein for different values of f. The free energy decreases as we increase the concentrated force f. The concentrated force is basically associated with the binding force acting between the receptor and spike glycoprotein during the endocytosis as shown in Figure 2 . If the driving force is increased, then the free energy of the whole system will be Table 1 . Here, we use κ = 10 k B T and h = 2R (full endocytosis). Table 1 . Here, we use h = 2R (full endocytosis) and σ = 3 pN μm À1 . Table 1 . Here, we use h = 2R (full endocytosis), σ = 3 pN μm À1 , and κ = 10k B T. decreased. For a particular f, the free energy decreases as we increase the concentration of F-BAR protein. For small f and θ, we obtain positive free energy of the system. However, larger values of f give us negative free energy of the system in the whole region of the protein concentrations as shown in Figure 7 . Positive free energy signifies the absorption of energy into the system and opposite is applicable to the negative free energy. We observe that the increase in f and θ decreases the free energy. This means more concentrated force adopts more F-BAR proteins on the membrane surface to achieve the endocytosis of the virus particle. We observe that the increasing membrane tension increases the BAR protein concentration on the membrane surface as shown in Figure 5 . This increasing protein concentration reduces the total free energy of the system depicted in Figure 7 . Thus, increasing membrane tension enhances the endocytosis process of the virus particle. However, it contradicts with the results of Deserno. 62 In his model, he did not consider the role of non-specific binding of BAR proteins on the membrane surface during membrane wrapping of colloidal particle. This might be the reason that Deserno's work 62 did not capture increase in endocytosis with increasing membrane tension. Current work can be extended to analyze endocytosis of BAR proteins while considering the kinetics of receptor-mediated endocytosis. 63 Also, a more accurate model while considering the deformation of virus particle, a detailed interaction and corresponding binding energy between protein and cell membrane can be incorporated in the current model to analyze the endocytosis of proteins. We develop a mechanical-thermodynamic model to explore the endocytosis behavior of SARS-CoV-2. With the help of present model, we observe that Endophilin N-BAR cytoplasmic protein does not participate in the endocytosis. The Amphiphysin N-BAR protein participates up to the half depth of insertion of the virus particle. However, the other cytoplasmic BAR proteins like F-BAR, PX-BAR, and Arfaptin BAR participate up to full depth of insertion and play important role in endocytosis. In general, the BAR proteins binding on cell membrane's inside surface assist virus for endocytosis. But, if the tension of the cell membrane is increased, then the BAR proteins from the cytoplasmic solution also bind on the membrane surface to assist virus for endocytosis. As we increase the bending modulus of the cell membrane, two phase solutions of the BAR proteins are observed for full endocytosis of the virus. We obtain an unstable region of protein concentration, for which the endocytosis cannot be achieved. Increasing concentrated force associated with the binding of cell membrane receptor and spike glycoprotein of virus during the endocytosis, adopts more proteins on the membrane surface to achieve endocytosis. Although the present study is focused on endocytosis of SARS-CoV-2, but this can also help one to understand endocytosis process involved in different viral infections. The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. The author(s) received no financial support for the research, authorship, and/or publication of this article. Paritosh Mahata  https://orcid.org/0000-0003-2882-5840 A pneumonia outbreak associated with a new coronavirus of probable bat origin A sequence homology and bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2 COVID-19 and cardiovascular disease: from basic mechanisms to clinical perspectives SARS-CoV-2 infects cells following viral entry via clathrin-mediated endocytosis Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune crossreactivity with SARS-CoV Endocytic uptake of SARS-CoV-2: the critical roles of pH, Ca 2+ and NAADP Molecular mechanisms of clathrin-independent endocytosis Molecular mechanism and physiological functions of clathrin-mediated endocytosis Curving cells inside and out: roles of BAR domain proteins in membrane shaping and its cellular implications When physics takes over: BAR proteins and membrane curvature Deciphering the BAR code of membrane modulators Mechanisms of clathrin-mediated endocytosis Caveolaedependent endocytosis in viral infection Domain proteins: a role in tubulation, scission and actin assembly in clathrin-mediated endocytosis Membrane remodeling in clathrinmediated endocytosis Small gtpases and BAR domain proteins regulate branched actin polymerisation for clathrin and dynamin-independent endocytosis Clathrin-dependent entry of severe acute respiratory syndrome coronavirus into target cells expressing ACE2 with the cytoplasmic tail deleted Coronavirus cell entry occurs through the endo-/lysosomal pathway in a proteolysis-dependent manner Sars coronavirus entry into host cells through a novel clathrin-and caveolae-independent endocytic pathway Cell entry mechanisms of SARS-CoV-2 Understanding SARS-CoV-2 endocytosis for COVID -19 drug repurposing Targeting the endocytic pathway and autophagy process as a novel therapeutic strategy in COVID-19 Computational investigation for endocytosis of CoVID-19 virus SARS-CoV-2 in cell membrane Receptormediated endocytosis of nanoparticles of various shapes Shape and orientation matter for the cellular uptake of nonspherical particles Receptor-mediated endocytosis of nanoparticles: roles of shapes, orientations, and rotations of nanoparticles Physical principles of nanoparticle cellular endocytosis Kinetics of receptor-mediated endocytosis of elastic nanoparticles Adhesive nanoparticles as local probes of membrane curvature Critical particle sizes for the engulfment of nanoparticles by membranes and vesicles with bilayer asymmetry Tailoring the component of protein corona via simple chemistry Interaction of serum proteins with sars-cov-2 rbd Elastic properties of lipid bilayers: theory and possible experiments Thermodynamics and mechanics of membrane curvature generation and sensing by proteins and lipids Two-dimensional convex-molecule fluid model for surface adsorption of proteins: effect of soft interaction on adsorption equilibria Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in wuhan, china: a descriptive study Curvature sorting of proteins on a cylindrical lipid membrane tether connected to a reservoir Nonlinear sorting, curvature generation, and crowding of endophilin N-BAR on tubular membranes Models for membrane curvature sensing of curvature generating proteins Sorting of proteins with shape and curvature anisotropy on a lipid bilayer tube Recent developments in membrane curvature sensing and induction by proteins Generation of wavy structure on lipid membrane by peripheral proteins: a linear elastic analysis The hydrophobic insertion mechanism of membrane curvature generation by proteins Modeling membrane shaping by proteins: focus on EHD2 and N-BAR domains Sensing membrane stresses by protein insertions Deformation of an elastic membrane interacting electrostatically with a rigid curved domain: implications to biosystems Drag force on a liquid domain moving inside a membrane sheet surrounded by aqueous medium Diffusion mediated coagulation and fragmentation based study of domain formation in lipid bilayer membrane Nature of curvature coupling of amphiphysin with membranes depends on its bound density Fluid lipid membranes -a primer Cell surface area regulation and membrane tension Membrane tension and cytoskeleton organization in cell motility Molecular driving forces: statistical thermodynamics in biology, chemistry, physics, and nanoscience A novel coronavirus from patients with pneumonia in china Membrane sculpting by F-BAR domains studied by molecular dynamics simulations Structural basis for membrane binding specificity of the Bin/Amphiphysin/Rvs (BAR) domain of Arfaptin-2 determined by Arl1 GTPase BAR domains as sensors of membrane curvature: the amphiphysin BAR structure Crystal structure of the endophilin-A1 BAR domain Pattern formation by curvature-inducing proteins on spherical membranes Adsorption of globular proteins on locally planar surfaces: models for the effect of excluded surface area and aggregation of adsorbed protein on adsorption equilibria Statistical mechanical theory of a new analytical equation of state Elastic deformation of a fluid membrane upon colloid binding Accurate evaluation on the interactions of SARS-CoV-2 with its receptor ACE2 and antibodies CR3022/CB6