Alper Tunga Celebi

Ion-specific effects at aqueous interfaces are fundamental to chemistry, biology, and environmental science, yet unraveling their thermodynamic origins remains challenging. Here, we quantitatively resolve interfacial thermodynamics at molecular resolution by combining concentration-dependent adsorption isotherms obtained from atomic force microscopy imaging, molecular dynamics simulations, and sum-frequency generation spectroscopy. We show that strongly hydrated ions disrupt and compact interfacial water structure, leading to entropic penalties and increased dielectric screening that weakens their interaction with the surface. Conversely, weakly hydrated ions incur lower entropy costs and substantially decrease screening via water depletion, promoting inner-sphere overadsorption on a charged surface. This balance between entropic penalties and screening governs inner double-layer formation. Our findings establish a thermodynamically grounded molecular framework for ion-specific adsorption, quantitatively linking interfacial entropy, water structure, and electrostatics. This approach generalizes the Hofmeister concept and provides a foundational approach for understanding interfacial structure in complex electrolytes.

Protein-based underwater adhesives of marine organisms exhibit extraordinary binding strength in high salinity based on utilizing a variety of molecular interaction mechanisms. These include acid-base interactions, bidentate bindings or complex hydrogen bonding interactions, and electrochemical manipulation of interfacial bonding. In this Perspective, we briefly review recent progress in the field, and we discuss how interfacial electrochemistry can vary interfacial forces by concerted tuning of surface charging, hydration forces, and tuning of the interfacial ion concentration. We further discuss open questions, controversial findings, and new paths into understanding and utilizing redox-proteins and derived polymers for enhancing underwater adhesion in a complex salt environment.

The human angiotensin-converting enzyme 2 (hACE2) receptor is the primary receptor for SARS-CoV-2 infection. However, the presence of alternative receptors such as the transmembrane glycoprotein CD147 has been proposed as a potential route for SARS-CoV-2 infection. The outcomes of SARS-CoV-2 spike protein binding to receptors have been shown to vary among individuals. Additionally, some patients infected with SARS-CoV-2 develop autoimmune responses. Given that CD147 is involved in the hyperactivation of memory T cells resulting in autoimmunity, we investigated the interaction of the SARS-CoV-2 viral spike protein with CD147 receptor and retinal specific CD147 Ig0 domain in silico using molecular docking and molecular dynamics (MD) simulations. The results indicated that binding involves two critical residues Lys63 and Asp65 in a ubiquitous CD147 isoform, potentially leading to the hyperactivation of T cells for only SARS-CoV-2, but not for SARS-CoV or MERS-CoV. Overall binding was confirmed by docking simulations. Next, MD analyses were completed to verify the docking poses. Polar interactions suggested that the interaction via Lys63 and Asp65 might be one of the determinants associated with severe COVID-19 outcomes. Neither did SARS-CoV nor MERS-CoV bind to these two critical residues when molecular docking analyses were performed. Interestingly, SARS-CoV was able to bind to CD147 with a lower affinity (-4.5 kcal/mol) than SARS-CoV-2 (-5.6 kcal/mol). Furthermore, Delta and Omicron variants of SARS-CoV-2 did not affect the polar interactions with Lys63 and Asp65 in CD147. This study further strengthens the link between SARS-CoV-2 infection and autoimmune responses and provides novel insights for prudent antiviral drug designs for COVID-19 treatment that have implications in the prevention of T cell hyperactivation.

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation–anion and ion–urea decreases, while the number of HBs between urea–urea increases. Radial distribution functions (RDFs) for ChCl–urea and ChCl–ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood–Buff Integrals (KBIs) are computed. KBIs show that interactions of urea–urea become stronger, while interactions of urea–ChCl and ChCl–ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

With the emergence of hydrophobic deep eutectic solvents (DESs), the scope of applications of DESs has been expanded to include situations in which miscibility with water is undesirable. Whereas most studies have focused on the applications of hydrophobic DESs from a practical standpoint, few theoretical works exist that investigate the structural and thermodynamic properties at the nanoscale. In this study, Molecular Dynamics (MD) simulations have been performed to model DESs composed of tetraalkylammonium chloride hydrogen bond acceptor and decanoic acid hydrogen bond donor (HBD) at a molar ratio of 1:2, with three different cation chain lengths (4, 7, and 8). After fine-tuning force field parameters, densities, viscosities, self-diffusivities, and ionic conductivities of the DESs were computed over a wide temperature range. The liquid structure was examined using radial distribution functions (RDFs) and hydrogen bond analysis. The MD simulations reproduced the experimental density and viscosity data from the literature reasonably well and were used to predict diffusivities and ionic conductivities, for which experimental data are scarce or unavailable. It was found that although an increase in the cation chain length considerably affected the density and transport properties of the DESs (i.e., yielding smaller densities and slower dynamics), no significant influence was observed on the RDFs and the hydrogen bonds. The self-diffusivities showed the following order for the mobility of the various components: HBD > anion > cation. Strong hydrogen bonds between the hydroxyl and carbonyl groups of decanoic acid and between the hydroxyl group of decanoic acid and chloride were observed to dominate the intermolecular interactions.

Deep Eutectic Solvents (DESs) are a new generation of green solvents, which are considered an environmentally friendly alternative to ionic liquids (ILs) and volatile organic compounds. The addition of controlled amounts of water to DESs has a signi cant e ect on their microscopic structure, and thus on their thermodynamic and transport properties. In this way DESs can be modi ed, leading to solvents with improved characteristics. In this work, molecular dynamics (MD) simulations are performed to obtain a better understanding of the relation between the microscopic structure, molecular interactions and thermophysical properties of aqueous reline and ethaline solutions at temperatures ranging from 303.15 K to 363.15 K. For both reline and ethaline solutions, the hydrogen bond (HB) networks disappear with increasing mass fraction of water, and the intensity of RDF peaks decreases. For a mass fraction of water of 40%, most of the HBs between the compounds of reline and ethaline are broken, and DESs are fully dissolved in water. Consequently, a monotonic decrease in viscosities, and an increase in self-di usion coecients are observed. Electrical conductivities show a non-monotonic behavior with increasing water content. Up to to 60% water mass fraction, the electrical conductivities increase with increasing water content. A further increase in the mass fraction of water decreases conductivities. For all studied systems, the hydrogen bond network and the peaks of RDFs show relatively small changes for water mass fractions below 5% and beyond 40%. The MD results show that viscosities decrease with temperature, while di usivities and electrical conductivities increase. The e ect of temperature on the structure of DES-water mixtures is negligible.

Slip lengths reported from molecular dynamics (MD) simulations of water flow in graphene nanochannels show significant scatter in the literature. These discrepancies are in part due to the used water models. We demonstrate self-consistent comparisons of slip characteristics between the SPC, SPC/E, SPC/Fw, TIP3P, TIP4P, and TIP4P/2005 water models. The slip lengths are inferred using an analytical model that employs the shear viscosity of water and channel average velocities obtained from nonequilibrium MD simulations. First, viscosities for each water model are quantified using MD simulations of counterflowing, force-driven flows in periodic domains in the absence of physical walls. While the TIP4P/2005 model predicts water viscosity at the specified thermodynamic state with 1.7% error, the predictions of SPC/Fw and SPC/E models exhibit 13.9% and 23.1% deviations, respectively. Water viscosities obtained from SPC, TIP4P, and TIP3P models show larger deviations. Next, force-driven water flows in rigid (cold) and thermally vibrating (thermal) graphene nanochannels are simulated, resulting in pluglike velocity profiles. Large differences in the flow velocities are observed depending on the used water model and to a lesser extent on the choice of rigid vs thermal walls. Depending on the water model, the slip length of water on cold graphene walls varied between 34.2 nm and 62.9 nm, while the slip lengths of water on thermal graphene walls varied in the range of 38.1 nm–84.3 nm.