The ongoing COVID-19 pandemic has infected millions of people, claimed hundreds of thousands of lives, and made a worldwide health emergency. Understanding the SARS-CoV-2 mechanism of infection is crucial in the development of potential therapeutics and vaccines. The infection process is triggered by direct binding of the SARS-CoV-2 receptor-binding domain (RBD) to the host cell receptor, angiotensin-converting enzyme 2 (ACE2). Many efforts have been made to design or repurpose therapeutics to deactivate the RBD or ACE2 and prevent the initial binding. In addition to direct inhibition strategies, small chemical compounds might be able to interfere and destabilize the meta-stable, pre-fusion complex of ACE2-RBD. This approach can be employed to prevent the further progress of virus infection at its early stages. In this study, molecular docking was employed to analyze the binding of two chemical compounds, SSAA09E2 and Nilotinib, with the druggable pocket of the ACE2-RBD complex. The structural changes as a result of the interference with the ACE2-RBD complex were analyzed by molecular dynamics simulations. Results show that both Nilotinib and SSAA09E2 can induce significant conformational changes in the ACE2-RBD complex, intervene with the hydrogen bonds, and influence the flexibility of proteins. Moreover, essential dynamics analysis suggests that the presence of small molecules can trigger large-scale conformational changes that may destabilize the ACE2-RBD complex.

Transient pore formation on the membrane of red blood cells (RBCs) under high mechanical tensions is of great importance in many biomedical applications, such as RBC damage (hemolysis) and mechanoporation-based drug delivery. The dynamic process of pore formation, growth, and resealing is hard to visualize in experiments. We developed a mesoscale coarse-grained model to study the characteristics of transient pores on a patch of the lipid bilayer that is strengthened by an elastic meshwork representing the cytoskeleton. Unsteady molecular dynamics was used to study the pore formation and reseal at high strain rates close to the physiological ranges. The critical strain for pore formation, pore characteristics, and cytoskeleton effects were studied. Results show that the presence of the cytoskeleton increases the critical strain of pore formation and confines the pore growth. Moreover, the pore recovery process under negative strain rates (compression) is analyzed. Simulations show that pores can remain open for a long time during the high-speed tank-treading induced stretching and compression process that a patch of the RBC membrane usually experiences under high shear flow. Furthermore, complex loading conditions can affect the pore characteristics and result in denser pores. Finally, the effects of strain rate on pore formation are analyzed. Higher rate stretching of membrane patch can result in a significant increase in the critical areal strain and density of pores. Such a model reveals the dynamic molecular process of RBC damage in biomedical devices and mechanoporation that, to our knowledge, has not been reported before.

This work introduces a new Lagrangian strain-based model to predict the shear-induced hemolysis in biomedical devices. Current computational models for device-induced hemolysis usually utilize empirical fitting of the released free hemoglobin (Hb) in plasma from damaged red blood cells (RBCs). These empirical correlations contain parameters that depend on specific device and operating conditions, thus cannot be used to predict hemolysis in a general device. The proposed algorithm

Intracellular drug delivery by rapid squeezing is one of the most recent and simple cell membrane disruption-mediated drug encapsulation approaches. In this method, cell membranes are perforated in a microfluidic setup due to rapid cell deformation during squeezing through constricted channels. While squeezing-based drug loading has been successful in loading drug molecules into various cell types, such as immune cells, cancer cells, and other primary cells, there is so far no comprehensive understanding of the pore opening mechanism on the cell membrane and the systematic analysis on how different channel geometries and squeezing speed influence drug loading. This article aims to develop a three-dimensional computational model to study the intracellular delivery for compound cells squeezing through microfluidic channels. The Lattice Boltzmann method, as the flow solver, integrated with a spring-connected network via frictional coupling, is employed to capture compound capsule dynamics over fast squeezing. The pore size is proportional to the local areal strain of triangular patches on the compound cell through mathematical correlations derived from molecular dynamics and coarse-grained molecular dynamics simulations. We quantify the drug concentration inside the cell cytoplasm by introducing a new mathematical model for passive diffusion after squeezing. Compared to the existing models, the proposed model does not have any empirical parameters that depend on operating conditions and device geometry. Since the compound cell model is new, it is validated by simulating a nucleated cell under a simple shear flow at different capillary numbers and comparing the results with other numerical models reported in literature. The cell deformation during squeezing is also compared with the pattern found from our compound cell squeezing experiment. Afterward, compound cell squeezing is modeled for different cell squeezing velocities, constriction lengths, and constriction widths. We reported the instantaneous cell center velocity, variations of axial and vertical cell dimensions, cell porosity, and normalized drug concentration to shed light on the underlying physics in fast squeezing-based drug delivery. Consistent with experimental findings in the literature, the numerical results confirm that constriction width reduction, constriction length enlargement, and average cell velocity promote intracellular drug delivery. The results show that the existence of the nucleus increases cell porosity and loaded drug concentration after squeezing. Given geometrical parameters and cell average velocity, the maximum porosity is achieved at three different locations: constriction entrance, constriction middle part, and outside the constriction. Our numerical results provide reasonable justifications for experimental findings on the influences of constriction geometry and cell velocity on the performance of cell-squeezing delivery. We expect this model can help design and optimize squeezing-based cargo delivery.

The constant thickness in the microfluidic channel is used for controlled absorption of red and blue light to measure red blood cell hemoglobin and height mapping. High speed recording of the height mapping provides us the membrane fluctuation.