Multiscale Modeling of Aqueous Electric Double Layers.

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Fast protein folding is governed by memory-dependent friction.

When described by a low-dimensional reaction coordinate, the folding rates of most proteins are determined by a subtle interplay between free-energy barriers, which separate folded and unfolded states, and friction. While it is commonplace to extract free-energy profiles from molecular trajectories, a direct evaluation of friction is far more elusive and typically relies on fits of measured reaction rates to memoryless reaction-rate theories. Here, using memory-kernel extraction methods founded on a generalized Langevin equation (GLE) formalism, we directly calculate the time-dependent friction acting on the fraction of native contacts reaction coordinate Q, evaluated for eight fast-folding proteins, taken from a published set of large-scale molecular dynamics protein simulations. Our results reveal that, across the diverse range of proteins represented in this dataset, friction is more influential than free-energy barriers in determining protein folding rates. We also show that proteins fold in a regime where the finite decay time of friction significantly reduces the folding times, in some instances by as much as a factor of 10, compared to predictions based on memoryless friction.

Nanoscopic Interfacial Hydrogel Viscoelasticity Revealed from Comparison of Macroscopic and Microscopic Rheology.

Deviations between macrorheological and particle-based microrheological measurements are often considered to be a nuisance and neglected. We study aqueous poly(ethylene oxide) (PEO) hydrogels for varying PEO concentrations and chain lengths that contain microscopic tracer particles and show that these deviations reveal the nanoscopic viscoelastic properties of the particle-hydrogel interface. Based on the transient Stokes equation, we first demonstrate that the deviations are not due to finite particle radius, compressibility, or surface-slip effects. Small-angle neutron scattering rules out hydrogel heterogeneities. Instead, we show that a generalized Stokes-Einstein relation, accounting for an interfacial shell around tracers with viscoelastic properties that deviate from bulk, consistently explains our macrorheological and microrheological measurements. The extracted shell diameter is comparable to the PEO end-to-end distance, indicating the importance of dangling chain ends. Our methodology reveals the nanoscopic interfacial rheology of hydrogels and is applicable to different kinds of viscoelastic fluids and particles.

Data-driven classification of individual cells by their non-Markovian motion.

We present a method to differentiate organisms solely by their motion based on the generalized Langevin equation (GLE) and use it to distinguish two different swimming modes of strongly confined unicellular microalgae Chlamydomonas reinhardtii. The GLE is a general model for active or passive motion of organisms and particles that can be derived from a time-dependent general many-body Hamiltonian and in particular includes non-Markovian effects (i.e., the trajectory memory of its past). We extract all GLE parameters from individual cell trajectories and perform an unbiased cluster analysis to group them into different classes. For the specific cell population employed in the experiments, the GLE-based assignment into the two different swimming modes works perfectly, as checked by control experiments. The classification and sorting of single cells and organisms is important in different areas; our method, which is based on motion trajectories, offers wide-ranging applications in biology and medicine.

Water Structuring Induces Nonuniversal Hydration Repulsion between Polar Surfaces: Quantitative Comparison between Molecular Simulations, Theory, and Experiments.

Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.

How Thick is the Air-Water Interface?─A Direct Experimental Measurement of the Decay Length of the Interfacial Structural Anisotropy.

The air-water interface is a highly prevalent phase boundary impacting many natural and artificial processes. The significance of this interface arises from the unique properties of water molecules within the interfacial region, with a crucial parameter being the thickness of its structural anisotropy, or "healing depth". This quantity has been extensively assessed by various simulations which have converged to a prediction of a remarkably short length of ∼6 Å. Despite the absence of any direct experimental measurement of this quantity, this predicted value has surprisingly become widely accepted as fact. Using an advancement in nonlinear vibrational spectroscopy, we provide the first measurement of this thickness and, indeed, find it to be ∼6-8 Å, finally confirming the prior predictions. Lastly, by combining the experimental results with depth-dependent second-order spectra calculated from ab initio parametrized molecular dynamics simulations, which are also in excellent agreement with this experimental result, we shed light on this surprisingly short correlation length of molecular orientations at the interface.

Collective modes and quantum effects in two-dimensional nanofluidic channels.

Nanoscale fluid transport is typically pictured in terms of atomic-scale dynamics, as is natural in the real-space framework of molecular simulations. An alternative Fourier-space picture, that involves the collective charge fluctuation modes of both the liquid and the confining wall, has recently been successful at predicting new nanofluidic phenomena such as quantum friction and near-field heat transfer, that rely on the coupling of those fluctuations. Here, we study the charge fluctuation modes of a two-dimensional (planar) nanofluidic channel. Introducing confined response functions that generalize the notion of surface response function, we show that the channel walls exhibit coupled plasmon modes as soon as the confinement is comparable to the plasmon wavelength. Conversely, the water fluctuations remain remarkably bulk-like, with significant confinement effects arising only when the wall spacing is reduced to 7 Å. We apply the confined response formalism to predict the dependence of the solid-water quantum friction and thermal boundary conductance on channel width for model channel wall materials. Our results provide a general framework for Coulomb interactions of fluctuating matter under nanoscale confinement.

Understanding the "Berg limit": the 65° contact angle as the universal adhesion threshold of biomatter.

Surface phenomena in aqueous environments such as long-range hydrophobic attraction, macromolecular adhesion, and even biofouling are predominantly influenced by a fundamental parameter-the water contact angle. The minimal contact angle required for these and related phenomena to occur has been repeatedly reported to be around 65° and is commonly referred to as the "Berg limit." However, the universality of this specific threshold across diverse contexts has remained puzzling. In this perspective article, we aim to rationalize the reoccurrence of this enigmatic contact angle. We show that the relevant scenarios can be effectively conceptualized as three-phase problems involving the surface of interest, water, and a generic oil-like material that is representative of the nonpolar constituents within interacting entities. Our analysis reveals that attraction and adhesion emerge when substrates display an underwater oleophilic character, corresponding to a "hydrophobicity under oil", which occurs for contact angles above approximately 65°. This streamlined view provides valuable insights into macromolecular interactions and holds implications for technological applications.

Force field for halide and alkali ions in water based on single-ion and ion-pair thermodynamic properties for a wide range of concentrations.

A classical non-polarizable force field for the common halide (F-, Cl-, Br-, and I-) and alkali (Li+, Na+, K+, and Cs+) ions in SPC/E water is presented. This is an extension of the force field developed by Loche et al. for Na+, K+, Cl-, and Br- (JPCB 125, 8581-8587, 2021): in the present work, we additionally optimize Lennard-Jones parameters for Li+, I-, Cs+, and F- ions. Li+ and F- are particularly challenging ions to model due to their small size. The force field is optimized with respect to experimental solvation free energies and activity coefficients, which are the necessary and sufficient quantities to accurately reproduce the electrolyte thermodynamics. Good agreement with experimental reference data is achieved for a wide range of concentrations (up to 4 mol/l). We find that standard Lorentz-Berthelot combination rules are sufficient for all ions except F-, for which modified combination rules are necessary. With the optimized parameters, we show that, although the force field is only optimized based on thermodynamic properties, structural properties are reproduced quantitatively, while ion diffusion coefficients are in qualitative agreement with experimental values.

Propensity of hydroxide and hydronium ions for the air-water and graphene-water interfaces from ab initio and force field simulations.

Understanding acids and bases at interfaces is relevant for a range of applications from environmental chemistry to energy storage. We present combined ab initio and force-field molecular dynamics simulations of hydrochloric acid and sodium hydroxide highly concentrated electrolytes at the interface with air and graphene. In agreement with surface tension measurements at the air-water interface, we find that HCl presents an ionic surface excess, while NaOH displays an ionic surface depletion, for both interfaces. We further show that graphene becomes less hydrophilic as the water ions concentration increases, with a transition to being hydrophobic for highly basic solutions. For HCl, we observe that hydronium adsorbs to both interfaces and orients strongly toward the water phase, due to the hydrogen bonding behavior of hydronium ions, which donate three hydrogen bonds to bulk water molecules when adsorbed at the interface. For NaOH, we observe density peaks of strongly oriented hydroxide ions at the interface with air and graphene. To extrapolate our results from concentrated electrolytes to dilute solutions, we perform single ion-pair ab initio simulations, as well as develop force-field parameters for ions and graphene that reproduce the density profiles at high concentrations. We find the behavior of hydronium ions to be rather independent of concentration. For NaOH electrolytes, the force-field simulations of dilute NaOH solutions suggest no hydroxide adsorption but some adsorption at high concentrations. For both interfaces, we predict that the surface potential is positive for HCl and close to neutral for NaOH.

Estimating computational limits on theoretical descriptions of biological cells.

There has been much success recently in theoretically simulating parts of complex biological systems on the molecular level, with the goal of first-principles modeling of whole cells. However, there is the question of whether such simulations can be performed because of the enormous complexity of cells. We establish approximate equations to estimate computation times required to simulate highly simplified models of cells by either molecular dynamics calculations or by solving molecular kinetic equations. Our equations place limits on the complexity of cells that can be theoretically understood with these two methods and provide a first step in developing what can be considered biological uncertainty relations for molecular models of cells. While a molecular kinetics description of the genetically simplest bacterial cell may indeed soon be possible, neither theoretical description for a multicellular system, such as the human brain, will be possible for many decades and may never be possible even with quantum computing.

Non-Markovian modeling of protein folding.

We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide [Formula: see text] for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote-Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects.

Functionalized Fullerene for Inhibition of SARS-CoV-2 Variants.

As virus outbreaks continue to pose a challenge, a nonspecific viral inhibitor can provide significant benefits, especially against respiratory viruses. Polyglycerol sulfates recently emerge as promising agents that mediate interactions between cells and viruses through electrostatics, leading to virus inhibition. Similarly, hydrophobic C60 fullerene can prevent virus infection via interactions with hydrophobic cavities of surface proteins. Here, two strategies are combined to inhibit infection of SARS-CoV-2 variants in vitro. Effective inhibitory concentrations in the millimolar range highlight the significance of bare fullerene's hydrophobic moiety and electrostatic interactions of polysulfates with surface proteins of SARS-CoV-2. Furthermore, microscale thermophoresis measurements support that fullerene linear polyglycerol sulfates interact with the SARS-CoV-2 virus via its spike protein, and highlight importance of electrostatic interactions within it. All-atom molecular dynamics simulations reveal that the fullerene binding site is situated close to the receptor binding domain, within 4 nm of polyglycerol sulfate binding sites, feasibly allowing both portions of the material to interact simultaneously.

Non-Markovian Modeling of Nonequilibrium Fluctuations and Dissipation in Active Viscoelastic Biomatter.

Based on a Hamiltonian that incorporates the elastic coupling between a tracer particle and the embedding active viscoelastic biomatter, we derive a generalized non-Markovian Langevin model for the nonequilibrium mechanical tracer response. Our analytical expressions for the frequency-dependent tracer response function and the tracer positional autocorrelation function agree quantitatively with experimental data for red blood cells and actomyosin networks with and without adenosine triphosphate over the entire frequency range and in particular reproduce the low-frequency violation of the fluctuation-dissipation theorem. The viscoelastic power laws, the elastic constants and effective friction coefficients extracted from the experimental data allow straightforward physical interpretation.

Insights into the structure and Ion transport of pectin-[BMIM][PF6] electrolytes.

We investigate the effect of pectin on the structure and ion transport properties of the room-temperature ionic liquid electrolyte 1-n-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]) using molecular dynamics simulations. We find that pectin induces intriguing structural changes in the electrolyte that disrupt large ionic aggregates and promote the formation of smaller ionic clusters, which is a promising finding for ionic conductivity. Due to pectin in [BMIM][PF6] electrolytes, the diffusion coefficient of cations and anions is observed to decrease by a factor of four for a loading of 25 wt. % of pectin in [BMIM][PF6] electrolyte. A strong correlation between the ionic diffusivities (D) and ion-pair relaxation timescales (τc) is observed such that D ∼ τc-0.75 for cations and D ∼ τc-0.82 for anions. The relaxation timescale exponents indicate that the ion transport mechanisms in pectin-[BMIM][PF6] electrolytes are slightly distinct from those found in neat [BMIM][PF6] electrolytes (D∼τc-1). Since pectin marginally affects ionic diffusivities at the gain of smaller ionic aggregates and viscosity, our results suggest that pectin-ionic liquid electrolytes offer improved properties for battery applications, including ionic conductivity, mechanical stability, and biodegradability.

Frequency-dependent hydrodynamic finite size correction in molecular simulations reveals the long-time hydrodynamic tail.

Finite-size effects are challenging in molecular dynamics simulations because they have significant effects on computed static and dynamic properties, in particular diffusion constants, friction coefficients, and time- or frequency-dependent response functions. We investigate the influence of periodic boundary conditions on the velocity autocorrelation function and the frequency-dependent friction of a particle in a fluid, and show that the long-time behavior (starting at the picosecond timescale) is significantly affected. We develop an analytical correction allowing us to subtract the periodic boundary condition effects. By this, we unmask the power-law long-time tails of the memory kernel and the velocity autocorrelation function in liquid water and a Lennard-Jones fluid from simulations with rather small box sizes.

Physics of virus transmission by speaking droplets.

To make the physics of person-to-person virus transmission from emitted droplets of oral fluid while speaking easily understood, we present simple and transparent algebraic equations that capture the essential physics of the problem. Calculations with these equations provide a straightforward way of determining whether emitted droplets remain airborne or rapidly fall to the ground, after accounting for the decrease in droplet size from water evaporation. At a relative humidity of 50%, for example, droplets with initial radii larger than about 50 µm rapidly fall to the ground, while smaller, potentially virus-containing droplets shrink in size from water evaporation and remain airborne for many minutes. Estimates of airborne virion emission rates while speaking strongly support the proposal that mouth coverings can help contain the COVID-19 pandemic.

Cavitation in lipid bilayers poses strict negative pressure stability limit in biological liquids.

Biological and technological processes that involve liquids under negative pressure are vulnerable to the formation of cavities. Maximal negative pressures found in plants are around -100 bar, even though cavitation in pure bulk water only occurs at much more negative pressures on the relevant timescales. Here, we investigate the influence of small solutes and lipid bilayers, both constituents of all biological liquids, on the formation of cavities under negative pressures. By combining molecular dynamics simulations with kinetic modeling, we quantify cavitation rates on biologically relevant length scales and timescales. We find that lipid bilayers, in contrast to small solutes, increase the rate of cavitation, which remains unproblematically low at the pressures found in most plants. Only when the negative pressures approach -100 bar does cavitation occur on biologically relevant timescales. Our results suggest that bilayer-based cavitation is what generally limits the magnitude of negative pressures in liquids that contain lipid bilayers.

Charge Matters: Mutations in Omicron Variant Favor Binding to Cells.

Evidence is strengthening to suggest that the novel SARS-CoV-2 mutant Omicron, with its more than 60 mutations, will spread and dominate worldwide. Although the mutations in the spike protein are known, the molecular basis for why the additional mutations in the spike protein that have not previously occurred account for Omicron's higher infection potential, is not understood. We propose, based on chemical rational and molecular dynamics simulations, that the elevated occurrence of positively charged amino acids in certain domains of the spike protein (Delta: +4; Omicron: +5 vs. wild type) increases binding to cellular polyanionic receptors, such as heparan sulfate due to multivalent charge-charge interactions. This observation is a starting point for targeted drug development.

Basement membrane stiffness determines metastases formation.

The basement membrane (BM) is a special type of extracellular matrix and presents the major barrier cancer cells have to overcome multiple times to form metastases. Here we show that BM stiffness is a major determinant of metastases formation in several tissues and identify netrin-4 (Net4) as a key regulator of BM stiffness. Mechanistically, our biophysical and functional analyses in combination with mathematical simulations show that Net4 softens the mechanical properties of native BMs by opening laminin node complexes, decreasing cancer cell potential to transmigrate this barrier despite creating bigger pores. Our results therefore reveal that BM stiffness is dominant over pore size, and that the mechanical properties of 'normal' BMs determine metastases formation and patient survival independent of cancer-mediated alterations. Thus, identifying individual Net4 protein levels within native BMs in major metastatic organs may have the potential to define patient survival even before tumour formation. The ratio of Net4 to laminin molecules determines BM stiffness, such that the more Net4, the softer the BM, thereby decreasing cancer cell invasion activity.