Incorporating Ion-Specific van der Waals and Soft Repulsive Interactions in the Poisson-Boltzmann Theory of Electrical Double Layers

Electrical double layers (EDLs) arise when an electrolyte is in contact with a charged surface, and are encountered in several application areas including batteries, supercapacitors, electrocatalytic reactors, and colloids. In the modeling of EDLs, a prominent knowledge gap has been the exclusion of van der Waals (vdW) and soft repulsive interactions in modified Poisson-Boltzmann (PB) theories. Although more short-ranged as compared to electrostatic interactions, we show here that vdW interactions can play an important role in determining the structure of the EDL via the formation of a Stern layer and in modulating the differential capacitance of an electrode in solution. To this end, we incorporate ion-ion and wall-ion vdW attraction and soft repulsion via a 12-6 Lennard-Jones (LJ) potential, resulting in a modified PB-LJ approach. The wall-ion LJ interactions were found to have a significant effect on the electrical potential and concentration profiles, especially close to the wall. However, ion-ion LJ interactions do not affect the EDL structure at low bulk ion concentrations (< 1 M). We also derive dimensionless numbers to quantify the impact of ion-ion and wall-ion LJ interactions on the EDL. Furthermore, in the pursuit of capturing ion-specific effects, we apply our model by considering various combinations of ions. We observe how varying parameters such as the electrolyte concentration and electrode potential affect the structure of the EDL due to the competition between ion-specific LJ and electrostatic interactions. Lastly, we show that the inclusion of vdW and soft repulsion interactions as well as hydration effects lead to a better qualitative agreement of the PB models with experimental double-layer differential capacitance data. Overall, the modified PB-LJ approach presented herein will lead to more accurate theoretical descriptions of EDLs in various application areas

How Does Arbidol Inhibit the Novel Coronavirus SARS-CoV-2? Atomistic Insights from Molecular Dynamics Simulations

The COVID-19 pandemic is spreading at an alarming rate, posing an unprecedented threat to the global economy and human health. Broad-spectrum antivirals are currently being administered for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) treatment. China\u27s prevention and treatment guidelines suggest the use of an anti-influenza drug, Arbidol, for the clinical treatment of COVID-19. Reports indicate that Arbidol could neutralize the SARS-CoV-2. Monotherapy with Arbidol is found superior to Lopinavir-Ritonavir or Favipiravir in the treatment of COVID-19. In the SARS-CoV-2, Arbidol acts upon interfering in virus binding to host cells. However, the detailed understanding of Arbidol induced inhibition of SARS-CoV-2 is not known. Here, we present atomistic insights into the Arbidol-induced SARS-CoV-2 membrane fusion inhibition and propose a model of inhibition. 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 RBD than ACE2. Analyses of the detailed decomposition of energy components and binding affinities revealed a substantial increase in the affinity between 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 induced structural rigidity in the virus glycoprotein resulting in restriction of the conformational rearrangements associated with membrane attachment and virus entry.Further, key residues of RBD and ACE2 that interacted with Arbidol were identified, opening the doors for the development of therapeutic strategies and higher efficacy Arbidol derivatives or lead drug candidates.</p

E-Z Isomerization in Guanidine: Second-order Saddle Dynamics, Non-statisticality, and Time-frequency Analysis

Our recent work on the E-Z isomerization reaction of guanidine using ab initio chemical dynamics simulations [Rashmi et al, Regul. Chaotic Dyn. 2021, 26, 119] emphasized the role of second-order saddle (SOS) in the isomerization reaction; however could not unequivocally establish the non-statistical nature of the dynamics followed in the reaction. In the present study, we performed thousands on-the-fly trajectories using forces computed at the MNDO level to investigate the influence of second-order saddle in the E-Z isomerization reaction of guanidine and the role of intramolecular vibrational energy redistribution (IVR) on the reaction dynamics. The simulations reveal that while majority of the trajectories follow the traditional transition state pathways, 15% of the trajectories follow the SOS path. The dynamics was found to be highly non-statistical with the survival probabilities of the reactants showing large deviations from those obtained within the RRKM assumptions. In addition, a detailed analysis of the dynamics using time-dependent frequencies and the frequency ratio spaces reveal the existence of multiple resonance junctions that indicate the existence of regular dynamics and long-lived quasi-periodic trajectories in the phase space associated with non-RRKM behavior