The COVID-19 pandemic has spread rapidly and posed an unprecedented threat to the global economy and human health. Broad-spectrum antivirals are currently being administered to treat severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). China's prevention and treatment guidelines suggest the use of an antiinfluenza drug, arbidol, for the clinical treatment of COVID-19. Reports indicate that arbidol could neutralize SARS-CoV-2. Monotherapy with arbidol is found to be superior to lopinavir-ritonavir or favipiravir for treating COVID-19. In SARS-CoV-2 infection, arbidol acts by interfering with viral binding to the host cells. However, the detailed mechanism through which arbidol induces the inhibition of SARS-CoV-2 is not known. Here, we present atomistic insights into the mechanism underlying membrane fusion inhibition by arbidol for SARS-CoV-2. Molecular dynamics (MD) simulation-based analyses demonstrate that arbidol binds and stabilizes at the receptor-binding domain (RBD)/ACE2 interface with a high affinity. It forms stronger intermolecular interactions with the RBD than ACE2. Analyses of the detailed decomposition of energy components and binding affinities revealed a substantial increase in the affinity between the RBD and ACE2 in the arbidol-bound RBD/ACE2 complex, suggesting that arbidol could generate favorable interactions between them. Based on our MD simulation results, we propose that the binding of arbidol induces structural rigidity in the viral glycoprotein, resulting in a restriction of the conformational rearrangements associated with membrane fusion and virus entry. Furthermore, key residues of the RBD and ACE2 that interact with arbidol were identified, opening the door for developing therapeutic strategies and higher-efficacy arbidol derivatives or lead drug candidates.

Atomic-scale defects are ubiquitous in nanomaterials, yet their role in modulating fluid flow is inadequately understood. Hexagonal boron nitride (hBN) is an important two-dimensional material with applications in desalination and osmotic power. Although pristine hBN offers higher friction to the flow of water than graphene, we show here that certain defects can enhance water slippage on hBN. Using classical molecular dynamics simulations assisted by quantum-mechanical density functional theory, we compute the friction coefficient of water on hBN containing various vacancies (B, N, BN, B 2 N, and B 3 N) and the Stone−Wales defect. By investigating two defect concentrations, we obtain friction coefficients ranging from 0.4 to 2.6 times that of pristine hBN, leading to a maximum water slip length of 18.1 nm on hBN with a N vacancy or a Stone−Wales defect. Our work informs the use of defects to tune water flow and reveals defective hBN as an alternative high-slip surface to graphene.

Biological membrane is a complex self-assembly of lipids, sterols, and proteins organized as a fluid bilayer of two closely stacked lipid leaflets. Differential molecular interactions among its diverse constituents give rise to heterogeneities in the membrane lateral organization. Under certain conditions, heterogeneities in the two leaflets can be spatially synchronized and exist as registered domains across the bilayer. Several contrasting theories behind mechanisms that induce registration of nanoscale domains have been suggested. Following a recent study showing the effect of position of lipid tail unsaturation on domain registration behavior, we decided to develop an analytical theory to elucidate the driving forces that create and maintain domain registry across leaflets. Towards this, we formulated a Hamiltonian for a stacked lattice system where site variables capture the lipid molecular properties such as the position of unsaturation and various other interactions that could drive phase separation and interleaflet coupling. We solve the Hamiltonian using Monte Carlo simulations and create a complete phase diagram that reports the presence or absence of registered domains as a function of various Hamiltonian parameters. We find that the interleaflet coupling should be described as a competing enthalpic contribution due to interaction of lipid tail termini, primarily due to saturated-saturated interactions, and an interleaflet entropic contribution from overlap of unsaturated tail termini. A higher position of unsaturation is seen to provide weaker interleaflet coupling. Thermodynamically stable nanodomains could also be observed for certain points in the parameter space in our bilayer model, which were further verified by carrying out extended Monte Carlo simulations. These persistent noncoalescing registered nanodomains close to the lower end of the accepted nanodomain size range also point towards a possible "nanoscale" emulsion description of lateral heterogeneities in biological membrane leaflets.