Showing 1–50 of 96 results for author: Netz, R R

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

Authors: Alexander P. Fellows, Álvaro Díaz Duque, Vasileios Balos, Louis Lehmann, Roland R. Netz, Martin Wolf, Martin Thämer

Abstract: The air-water interface is a highly prevalent phase boundary with a far-reaching impact on natural and industrial processes. Water molecules behave differently at the interface compared to the bulk, exhibiting anisotropic orientational distributions, reduced intermolecular connectivity in the hydrogen bond network, and significantly slower dynamics. Despite many decades of research, the thickness… ▽ More The air-water interface is a highly prevalent phase boundary with a far-reaching impact on natural and industrial processes. Water molecules behave differently at the interface compared to the bulk, exhibiting anisotropic orientational distributions, reduced intermolecular connectivity in the hydrogen bond network, and significantly slower dynamics. Despite many decades of research, the thickness of the structural anisotropy in the interfacial layer remains controversial, with a direct experimental measurement being absent. In this study, we utilise an advancement in non-linear vibrational spectroscopy to gain access to this important parameter. Combining phase-resolved sum- and difference-frequency generation (SFG and DFG) responses, we directly measure the decay in structural anisotropy of the air-water interface. We find a decay length of ~6-8Å, in excellent agreement with depth-resolved SFG spectra calculated from ab initio parameterised molecular dynamics (MD) simulations. The result reveals surprisingly short anisotropic orientational correlations from the interfacial layer that are even shorter than in the bulk. Furthermore, the recorded SFG and DFG responses are decomposed into a vibrationally resonant and non-resonant contribution through isotopic exchange measurements. Through their separate analysis, we show that the resonant response is a sensitive probe of the structural anisotropy at the interface whereas the non-resonant contribution contains a significant isotropic contribution from the bulk and therefore only partially reports on the interfacial structure. This finding places stringent restrictions on the insight available through both purely non-resonant and second-order intensity studies. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 18 pages, 4 figures, 1 table, supplementary information

Quantifying the efficiency of principal signal transmission modes in proteins

Authors: Anil Kumar Sahoo, Richard Schwarzl, Markus S. Miettinen, Roland R. Netz

Abstract: On the microscopic level, biological signal transmission relies on coordinated structural changes in allosteric proteins that involve sensor and effector modules. The timescales and microscopic details of signal transmission in proteins are often unclear, despite a plethora of structural information on signaling proteins. Based on linear-response theory, we develop a theoretical framework to defin… ▽ More On the microscopic level, biological signal transmission relies on coordinated structural changes in allosteric proteins that involve sensor and effector modules. The timescales and microscopic details of signal transmission in proteins are often unclear, despite a plethora of structural information on signaling proteins. Based on linear-response theory, we develop a theoretical framework to define frequency-dependent force and displacement transmit functions through proteins and, more generally, viscoelastic media. Transmit functions quantify the fraction of a local time-dependent perturbation at one site, be it a deformation, a force, or a combination thereof, that survives at a second site. They are defined in terms of equilibrium fluctuations from simulations or experimental observations. We apply the framework to our all-atom molecular dynamics simulation data of a parallel, homodimeric coiled-coil (CC) motif that connects sensor and effector modules of a blue-light-regulated histidine kinase from bacterial signaling systems extensively studied in experiments. Our analysis reveals that signal transmission through the CC is possible via shift, splay, and twist deformation modes. Based on the results of mutation experiments, we infer that the most relevant mode for the biological function of the histidine kinase protein is the splay deformation. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Rapid state-recrossing kinetics in non-Markovian systems

Authors: Qingyuan Zhou, Roland R. Netz, Benjamin A. Dalton

Abstract: The mean first-passage time (MFPT) is one standard measure for the reaction time in thermally activated barrier-crossing processes. While the relationship between MFPTs and phenomenological rate coefficients is known for systems that satisfy Markovian dynamics, it is not clear how to interpret MFPTs for experimental and simulation time-series data generated by non-Markovian systems. Here, we simul… ▽ More The mean first-passage time (MFPT) is one standard measure for the reaction time in thermally activated barrier-crossing processes. While the relationship between MFPTs and phenomenological rate coefficients is known for systems that satisfy Markovian dynamics, it is not clear how to interpret MFPTs for experimental and simulation time-series data generated by non-Markovian systems. Here, we simulate a one-dimensional generalized Langevin equation (GLE) in a bistable potential and compare two related numerical methods for evaluating MFPTs: one that only incorporates information about first arrivals between subsequent states and is equivalent to calculating the waiting time, or dwell time, and one that incorporates information about all first-passages associated with a given barrier-crossing event and is therefore typically employed to enhance numerical sampling. In the Markovian limit, the two methods are equivalent. However, for significant memory times, the two methods suggest dramatically different reaction kinetics. By focusing on first-passage distributions, we systematically reveal the influence of memory-induced rapid state-recrossing on the MFPTs, which we compare to various other numerical or theoretical descriptions of reaction times. Overall, we demonstrate that it is necessary to consider full first-passage distributions, rather than just the mean barrier-crossing kinetics when analyzing non-Markovian time series data. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 7 pages, 5 figures

pH modulates friction memory effects in protein folding

Authors: Benjamin A. Dalton, Roland R. Netz

Abstract: We study the non-Markovian folding dynamics of the $α$3D protein under low- and neutral-pH conditions. Recently published all-atom simulations of $α$3D by the Shaw group reveal that lowering the pH significantly reduces both native and non-native salt-bridge interactions, which dominate the folding dynamics. Here, we demonstrate that this physiochemical modulation directly perturbs the folding fri… ▽ More We study the non-Markovian folding dynamics of the $α$3D protein under low- and neutral-pH conditions. Recently published all-atom simulations of $α$3D by the Shaw group reveal that lowering the pH significantly reduces both native and non-native salt-bridge interactions, which dominate the folding dynamics. Here, we demonstrate that this physiochemical modulation directly perturbs the folding friction, which we evaluate using non-Markovian memory-kernel-extraction techniques. In doing so, we find that the reduction in pH not only decreases the magnitude of the time-dependent friction acting on the protein but also more dramatically shortens the time scale of the friction memory effects. As a result, the folding dynamics in the low pH system are well described by a purely Markovian model. In the neutral pH system, however, the memory time scale is of the same order as the folding time and is accelerated by a factor of 6 compared to a Markovian model prediction. We demonstrate that this memory-induced barrier-crossing speed-up is predicted by non-Markovian reaction-kinetic theories, confirming that non-Markovian models are, in general, necessary for a quantitative description of protein folding dynamics. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 4 pages, 4 figures

Accurate Memory Kernel Extraction from Discretized Time Series Data

Authors: Lucas Tepper, Benjamin Dalton, Roland R. Netz

Abstract: Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate and insightful, particularly when working w… ▽ More Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate and insightful, particularly when working with high-resolution time series data. However, in experimental systems, high-resolution data is often unavailable, raising questions about the impact of the data resolution on the estimated GLE parameters. This study demonstrates that direct memory extraction remains accurate when the discretization time is below the memory time. To obtain memory functions reliably even when the discretization time exceeds the memory time, we introduce a Gaussian Process Optimization (GPO) scheme. This scheme minimizes the deviation of discretized two-point correlation functions between MD and GLE simulations and is able to estimate accurate memory kernels as long as the discretization time stays below the longest time scale in the data, typically the barrier crossing time. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Journal ref: J. Chem. Theory Comput. 2024, 20, 8

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

Authors: Robert F. Schmidt, Henrik Kiefer, Robert Dalgliesh, Michael Gradzielski, Roland R. Netz

Abstract: Deviations between macrorheological and particle-based microrheological measurements are often considered 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 in fact reveal the nanoscopic viscoelastic properties of the particle-hydrogel interface. Ba… ▽ More Deviations between macrorheological and particle-based microrheological measurements are often considered 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 in fact 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 a nanoscopic interfacial shell around tracers with viscoelastic properties that significantly deviate from bulk, consistently explains our macrorheological and microrheological measurements. The extracted shell diameter is comparable with the PEO end-to-end distance, indicating the importance of dangling chain ends. Our methodology reveals the nanoscopic interfacial rheology of hydrogels and is generally applicable to different kinds of viscoelastic fluids and particles. △ Less

Submitted 16 February, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: 25 pages, 4 figures, 52 pages of Supporting Information

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

Authors: Anton Klimek, Debasmita Mondal, Stephan Block, Prerna Sharma, Roland R. Netz

Abstract: 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 (CR). The GLE is the most general model for active or passive motion of organisms and particles and in particular includes non-Markovian effects, i.e., the… ▽ More 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 (CR). The GLE is the most general model for active or passive motion of organisms and particles 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 that is based on motion trajectories offers wide-ranging applications in biology and medicine. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 15 Pages main text including 5 figures, then 22 Pages Supplementary Information (SI) including 8 figures and 2 video files

Journal ref: Biophysical Journal 123, 1173-1183 (2024)

Multi-point distribution for Gaussian non-equilibrium non-Markovian observables

Authors: Roland R Netz

Abstract: When analyzing experimental or simulation time-series data, the question arises whether it is possible to tell from a one-dimensional time-dependent trajectory whether the system is in equilibrium or not. We here consider the non-equilibrium version of the generalized Langevin equation for a Gaussian observable and show that i) the multi-point joint distribution solely depends on the two-point cor… ▽ More When analyzing experimental or simulation time-series data, the question arises whether it is possible to tell from a one-dimensional time-dependent trajectory whether the system is in equilibrium or not. We here consider the non-equilibrium version of the generalized Langevin equation for a Gaussian observable and show that i) the multi-point joint distribution solely depends on the two-point correlation function and that ii) the two-point correlation function for a non-equilibrium process is identical to an equilibrium process with uniquely determined parameters. Since the multi-point joint distribution completely characterizes the dynamics of an observable, this means that the non-equilibrium character of a system, in contrast to its non-Markovianity, cannot be read off from a one-dimensional trajectory. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: arXiv admin note: text overlap with arXiv:2310.00748

Derivation of the non-equilibrium generalized Langevin equation from a generic time-dependent Hamiltonian

Authors: Roland R. Netz

Abstract: It has been become standard practice to describe steady-state non-equilibrium phenomena by Langevin equations with colored noise and time-dependent friction kernels that do not obey the fluctuation-dissipation theorem, but since these Langevin equations are typically not derived from first-principle Hamiltonian dynamics it is not clear whether they correspond to physically realizable scenarios. By… ▽ More It has been become standard practice to describe steady-state non-equilibrium phenomena by Langevin equations with colored noise and time-dependent friction kernels that do not obey the fluctuation-dissipation theorem, but since these Langevin equations are typically not derived from first-principle Hamiltonian dynamics it is not clear whether they correspond to physically realizable scenarios. By exact Mori projection in phase space we derive the non-equilibrium generalized Langevin equation (GLE) from a generic many-body Hamiltonian with a time-dependent force h(t) acting on an arbitrary phase-space dependent observable $A$. The GLE is obtained in explicit form to all orders in $h(t)$. For non-equilibrium observables that correspond to a Gaussian process, the resultant GLE has the same form as the equilibrium Mori GLE, in particular the memory kernel is proportional to the total force autocorrelation function. This means that the extraction and simulation methods developed for equilibrium GLEs can be used also for non-equilibrium Gaussian variables. This is a non-trivial and very useful result, as many observables that characterize non-equilibrium systems display Gaussian statistics. For non-Gaussian non-equilibrium variables correction terms appear in the GLE and in the relation between the complementary force autocorrelations and the memory kernels, which are explicitly given in term of cubic correlation functions of $A$. Interpreting the time-dependent force h(t) as a stochastic process, we derive non-equilibrium corrections to the fluctuation-dissipation theorem and methods to extract all GLE parameters from experimental or simulation data, thus making our non-equilibrium GLE a practical tool to study and model general non-equilibrium systems. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Conformational isomerization dynamics in solvent violates both the Stokes-Einstein relation and Kramers' theory

Authors: Benjamin A. Dalton, Henrik Kiefer, Roland R. Netz

Abstract: Molecular isomerization kinetics in liquid solvents are determined by a complex interplay between the friction acting on a rotating dihedral due to interactions with the solvent, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the free energy profile over which a dihedral rotates. Currently, it is not understood how these quantities are related at… ▽ More Molecular isomerization kinetics in liquid solvents are determined by a complex interplay between the friction acting on a rotating dihedral due to interactions with the solvent, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the free energy profile over which a dihedral rotates. Currently, it is not understood how these quantities are related at the molecular scale. Here, we combine molecular dynamics simulations of isomerizing n-alkane chains and dipeptide molecules in mixed water-glycerol solvents with memory-kernel extraction techniques to directly evaluate the frequency-dependent friction acting on a rotating dihedral. We extract the friction and isomerization times over a range of glycerol concentrations and accurately evaluate the relationships between solvent viscosity, isomerization kinetics, and dihedral friction. We show that the total friction acting on a rotating dihedral does not scale linearly with solvent viscosity, thus violating the Stokes-Einstein relation. Additionally, we demonstrate that the kinetics of isomerization are significantly faster compared to the Kramers prediction in the overdamped limit. We suggest that isomerization kinetics are determined by the multi-time-scale friction coupling between a rotating dihedral and its solvent environment, which results in non-Markovian kinetic speed-up effects. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 9 pages, 4 figures

Collective modes and quantum effects in two-dimensional nanofluidic channels

Authors: Baptiste Coquinot, Maximilian Becker, Roland R. Netz, Lydéric Bocquet, Nikita Kavokine

Abstract: 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… ▽ More 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 A. 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 in nanoscale confinement. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Insights into the Structure and Ion Transport of Pectin-[BMIM][PF$_6$] Electrolytes

Authors: Sipra Mohapatra, Sougata Halder, Sachin R. Chaudhary, Roland R. Netz, Santosh Mogurampelly

Abstract: 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, w… ▽ More 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-IL electrolytes offer improved properties for battery applications, including ionic conductivity, mechanical stability, and biodegradability. △ Less

Submitted 30 September, 2023; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Journal of Chemical Physics

Ion Transport Mechanisms in Pectin-containing EC-LiTFSI Electrolytes

Authors: Sipra Mohapatra, Hema Teherpuria, Sapta Sindhu Paul Chowdhury, Suleman Jalilahmad Ansari, Prabhat K Jaiswal, Roland R. Netz, Santosh Mogurampelly

Abstract: Using all-atom molecular dynamics simulations, we report the structure and ion transport characteristics of a new class of solid polymer electrolytes that contain biodegradable and mechanically stable biopolymer pectin. We simulate a highly conducting ethylene carbonate (EC) as a solvent for lithium-trifluoromethanesulfonimide (LiTFSI) salt containing different weight percentages of pectin. Our si… ▽ More Using all-atom molecular dynamics simulations, we report the structure and ion transport characteristics of a new class of solid polymer electrolytes that contain biodegradable and mechanically stable biopolymer pectin. We simulate a highly conducting ethylene carbonate (EC) as a solvent for lithium-trifluoromethanesulfonimide (LiTFSI) salt containing different weight percentages of pectin. Our simulations reveal that the pectin chains reduce the coordination numbers of lithium ions around the counterions (and vice-versa) because of stronger lithium-pectin interactions compared to lithium-TFSI interactions. Further, the pectin is found to promote smaller ionic aggregates over larger ones, contrary to results typically reported for liquid and polymer electrolytes. We observe that the loading of pectin in EC-LiTFSI electrolytes increases the viscosity ($η$) and relaxation timescales ($τ_c$), indicating higher mechanical stability and, consequently, the mean squared displacements, diffusion coefficients (D), and the Nernst-Einstein conductivity ($σ_{NE}$) decrease. Interestingly, while the lithium diffusivities are related to the ion-pair relaxation timescales as $D_+\sim τ_c^{-3.1}$, the TFSI$^-$ diffusivities exhibit excellent correlations with ion-pair relaxation timescales as $D_-\sim τ_c^{-0.95}$. On the other hand, the NE conductivities are dictated by distinct transport mechanisms and scales with ion-pair relaxation timescale as $σ_{NE}\sim τ_c^{-1.85}$. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Dielectric properties of aqueous electrolytes at the nanoscale

Authors: Maximilian R. Becker, Philip Loche, Douwe Jan Bonthuis, Dominique Mouhanna, Roland R. Netz, Hélène Berthoumieux

Abstract: Despite the ubiquity of nanoconfined aqueous electrolytes, a theoretical framework that accounts for the nonlinear coupling of water and ion polarization is still missing. We introduce a nonlocal and nonlinear field theory for the nanoscale polarization of ions and water and derive the electrolyte dielectric properties as a function of salt concentration to first order in a loop expansion. Classic… ▽ More Despite the ubiquity of nanoconfined aqueous electrolytes, a theoretical framework that accounts for the nonlinear coupling of water and ion polarization is still missing. We introduce a nonlocal and nonlinear field theory for the nanoscale polarization of ions and water and derive the electrolyte dielectric properties as a function of salt concentration to first order in a loop expansion. Classical molecular dynamics simulations are favorably compared with the calculated dielectric response functions. The theory correctly predicts the dielectric permittivity decrement with rising salt concentration and furthermore shows that salt induces a Debye screening in the longitudinal susceptibility but leaves the short-range water organization remarkably unchanged. △ Less

Submitted 26 March, 2023; originally announced March 2023.

Comments: 6 pages, 3 figures

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

Authors: Laura Scalfi, Domenico Vitali, Henrik Kiefer, Roland R. Netz

Abstract: 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… ▽ More 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 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 rather small simulation boxes. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: 6 pages, 4 figures, 9 pages of Appendices

Fast protein folding is governed by memory-dependent friction

Authors: Benjamin A. Dalton, Cihan Ayaz, Lucas Tepper, Roland R. Netz

Abstract: When described by a low-dimensional reaction coordinate, the rates of protein folding are determined by a subtle interplay between free-energy barriers and friction. While it is commonplace to extract free-energy profiles from molecular trajectories, a direct evaluation of friction is far more elusive, and one typically evaluates it indirectly via memoryless reaction rate theories. Here, using mem… ▽ More When described by a low-dimensional reaction coordinate, the rates of protein folding are determined by a subtle interplay between free-energy barriers and friction. While it is commonplace to extract free-energy profiles from molecular trajectories, a direct evaluation of friction is far more elusive, and one typically evaluates it indirectly via memoryless reaction rate theories. Here, using memory-kernel extraction methods founded on a generalised Langevin equation (GLE) formalism, we directly calculate the memory-dependent friction for eight fast-folding proteins, taken from a published set of large-scale molecular dynamics protein simulations. Our results reveal that, contrary to common expectation, friction is more important than free energy barriers in determining protein folding rates, particularly for larger proteins. 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. △ Less

Submitted 21 October, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 10 pages, 4 figures

Non-Markovian modeling of non-equilibrium fluctuations and dissipation in active viscoelastic biomatter

Authors: Amir Abbasi, Roland R. Netz, Ali Naji

Abstract: Based on a Hamiltonian that incorporates the elastic coupling between a tracer and active particles, we derive a generalized Langevin model for the non-equilibrium mechanical response of active viscoelastic biomatter. Our model accounts for the power-law viscoelastic response of the embedding polymeric network as well as for the non-equilibrium energy transfer between active and tracer particles.… ▽ More Based on a Hamiltonian that incorporates the elastic coupling between a tracer and active particles, we derive a generalized Langevin model for the non-equilibrium mechanical response of active viscoelastic biomatter. Our model accounts for the power-law viscoelastic response of the embedding polymeric network as well as for the non-equilibrium energy transfer between active and tracer particles. Our analytical expressions for the frequency-dependent response function and the positional autocorrelation function agree nicely with experimental data for red blood cells and actomyosin networks with and without ATP. The fitted effective active-particle temperature, elastic constants and effective friction coefficients of our model allow straightforward physical interpretation. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 6 pages, 2 figures; supplemental material (16 pages, 5 figures) is appended

Barrier-crossing times for different non-Markovian friction in well and barrier $-$ A numerical study

Authors: Florian N. Brünig, Roland R. Netz, Julian Kappler

Abstract: We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the dependence of the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier m… ▽ More We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the dependence of the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier memory time can differ by several orders of magnitude. Our model furthermore allows to test an analytic rate theory for space-inhomogeneous memory, which disagrees with our numerical results in the long well-memory regime. In this regime, we find that short barrier memory decreases the barrier-crossing time as compared to long barrier memory. This is in contrast with the short well-memory regime, where both our numerical results and GH theory predict an acceleration of the barrier crossing time with increasing barrier memory time. Both effects, the `Markovian-barrier acceleration' and GH `non-Markovian-barrier acceleration' can be understood from a committor analysis. Our model combines finite relaxation times of orthogonal degrees of freedom with a space-inhomogeneous coupling to such degrees, and represents a step towards more realistic modeling of physical reaction coordinates. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

Authors: Florian N. Brünig, Manuel Rammler, Ellen M. Adams, Martina Havenith, Roland R. Netz

Abstract: Signatures of solvated excess protons in infrared difference absorption spectra, such as the continuum band between the water bend and stretch bands, have been experimentally known for a long time, but the theoretical basis for linking spectral signatures with the microscopic proton-transfer mechanism so far relied on normal-mode analysis. We analyze the excess-proton dynamics in ab initio molecul… ▽ More Signatures of solvated excess protons in infrared difference absorption spectra, such as the continuum band between the water bend and stretch bands, have been experimentally known for a long time, but the theoretical basis for linking spectral signatures with the microscopic proton-transfer mechanism so far relied on normal-mode analysis. We analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric acid solutions by trajectory-decomposition techniques. The continuum band in the 2000 - 3000 cm$^{-1 }$ range is shown to be due to normal-mode oscillations of temporary H$_3$O$^+$ complexes. An additional prominent peak at 400 cm$^{-1}$ reports on the coupling of excess-proton motion to the relative vibrations of the two flanking water molecules. The actual proton transfer between two water molecules, which for large water separations involves crossing of a barrier and thus is not a normal mode, is characterized by two characteristic time scales: Firstly, the waiting time for transfer to occur in the range of 200 - 300 fs, which leads to a broad weak shoulder around ~100 cm$^{-1}$, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm$^{-1}$ and agrees in location and width with previous experimental mid-infrared spectra. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Journal ref: Nature Communications 13, 4210 (2022)

Generalized Langevin Equation with a Non-Linear Potential of Mean Force and Non-Linear Memory Friction From a Hybrid Projection Scheme

Authors: Cihan Ayaz, Benjamin A. Dalton, Roland R. Netz

Abstract: We introduce a hybrid projection scheme that combines linear Mori projection and conditional Zwanzig projection techniques and use it to derive a Generalized Langevin Equation (GLE) for a general interacting many-body system. The resulting GLE includes i) explicitly the potential of mean force (PMF) that describes the equilibrium distribution of the system in the chosen space of reaction coordinat… ▽ More We introduce a hybrid projection scheme that combines linear Mori projection and conditional Zwanzig projection techniques and use it to derive a Generalized Langevin Equation (GLE) for a general interacting many-body system. The resulting GLE includes i) explicitly the potential of mean force (PMF) that describes the equilibrium distribution of the system in the chosen space of reaction coordinates, ii) a random force term that explicitly depends on the initial state of the system, and iii) a memory friction contribution that splits into two parts: a part that is linear in the past reaction-coordinate velocity and a part that is in general non-linear in the past reaction coordinates but does not depend on velocities. Our hybrid scheme thus combines all desirable properties of the Zwanzig and Mori projection schemes. The non-linear memory friction contribution is shown to be related to correlations between the reaction-coordinate velocity and the random force. We present a numerical method to compute all parameters of our GLE, in particular, the non-linear memory friction function and the random force distribution, from a trajectory in reaction coordinate space. We apply our method to the dihedral-angle dynamics of a butane molecule in water obtained from atomistic molecular dynamics simulations. For this example, we demonstrate that non-linear memory friction is present and that the random force exhibits significant non-Gaussian corrections. We also present the derivation of the GLE for multidimensional reaction coordinates that are general functions of all positions in the phase space of the underlying many-body system; this corresponds to a systematic coarse-graining procedure that preserves not only the correct equilibrium behavior but also the correct dynamics of the coarse-grained system. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Journal ref: Phys. Rev. E 105 (2022)

Time-dependent friction effects on vibrational infrared frequencies and line shapes of liquid water

Authors: Florian N. Brünig, Otto Geburtig, Alexander von Canal, Julian Kappler, Roland R. Netz

Abstract: From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged non-linear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits adiabatic contributions at and below the vibrational time scales, but also much slower non-adiabatic contributions, reflecting homogeneous and inhomogeneous line broadening m… ▽ More From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged non-linear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits adiabatic contributions at and below the vibrational time scales, but also much slower non-adiabatic contributions, reflecting homogeneous and inhomogeneous line broadening mechanisms, respectively. Compared to the gas phase, hydration softens both stretch and bend potentials, which by itself would lead to a red-shift of the corresponding vibrational bands. In contrast, non-adiabatic friction contributions cause a spectral blue shift. For the stretch mode, the potential effect dominates and thus a significant red shift when going from gas to the liquid phase results. For the bend mode, potential and non-adiabatic friction effects are of comparable magnitude, so that a slight blue shift results, in agreement with well-known but puzzling experimental findings. The observed line broadening is shown to be roughly equally caused by adiabatic and non-adiabatic friction contributions for both, the stretch and bend modes in liquid water. Thus, the understanding of infrared vibrational frequencies and line shapes is considerably advanced by the quantitative analysis of the time-dependent friction that acts on vibrational modes in liquids. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Journal ref: J. Phys. Chem. B 2022, 126, 7, 1579-1589

Role of Entropy in Determining the Phase Behavior of Protein Solutions Induced by Multivalent Ions

Authors: Anil Kumar Sahoo, Frank Schreiber, Roland R. Netz, Prabal K. Maiti

Abstract: Recent experiments have reported lower critical solution temperature (LCST) phase behavior of aqueous solutions of proteins induced by multivalent ions, where the solution phase separates upon heating. This phenomenon is linked to complex hydration effects that result in a net entropy gain upon phase separation. To decipher the underlying molecular mechanism, we use all-atom molecular dynamics sim… ▽ More Recent experiments have reported lower critical solution temperature (LCST) phase behavior of aqueous solutions of proteins induced by multivalent ions, where the solution phase separates upon heating. This phenomenon is linked to complex hydration effects that result in a net entropy gain upon phase separation. To decipher the underlying molecular mechanism, we use all-atom molecular dynamics simulations along with the two-phase thermodynamic method for entropy calculation. Based on simulations of a single BSA protein in various salt solutions (NaCl, CaCl_2, MgCl_2, and YCl_3) at temperatures (T) ranging 283-323 K, we find that the cation-protein binding affinity increases with T, reflecting its thermodynamic driving force to be entropic in origin. We show that in the cation binding process, many tightly bound water molecules from the solvation shells of a cation and the protein are released to the bulk, resulting in entropy gain. To rationalize the LCST behavior, we calculate the ζ-potential that shows charge inversion of the protein for solutions containing multivalent ions. The ζ-potential increases with T. Performing simulations of two BSA proteins, we demonstrate that the protein-protein binding is mediated by multiple cation bridges and involves similar dehydration effects that cause a large entropy gain which more than compensates for rotational and translational entropy losses of the proteins. Thus, the LCST behavior is entropy-driven, but the associated solvation effects are markedly different from hydrophobic hydration. Our findings have direct implications for tuning the phase behavior of biological and soft-matter systems, e.g., protein condensation and crystallization. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Journal ref: Soft Matter, 2022

Proton-transfer spectroscopy beyond the normal-mode scenario

Authors: Florian N. Brünig, Paul Hillmann, Won Kyu Kim, Jan O. Daldrop, Roland R. Netz

Abstract: A stochastic theory is developed to predict the spectral signature of proton-transfer processes and applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single H$_5$O$_2{}^{+}$ cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contr… ▽ More A stochastic theory is developed to predict the spectral signature of proton-transfer processes and applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single H$_5$O$_2{}^{+}$ cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively report on the proton-transfer kinetics. In accordance with reaction kinetic theory, the transfer-waiting-contribution frequency depends inverse exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height. △ Less

Submitted 28 July, 2022; v1 submitted 17 September, 2021; originally announced September 2021.

Water evaporation from solute-containing aerosol droplets: effects of internal concentration and diffusivity profiles and onset of crust formation

Authors: Majid Rezaei, Roland R. Netz

Abstract: Saliva is primarily composed of water, but additionally includes a variety of organic and inorganic substances such as salt, proteins, peptides, mucins, virions, etc. The presence of such solutes affects the evaporation time of respiratory droplets that are sedimenting in air, and thereby the airborne transmission of infections. From solutions of the coupled heat-conduction and water-diffusion equ… ▽ More Saliva is primarily composed of water, but additionally includes a variety of organic and inorganic substances such as salt, proteins, peptides, mucins, virions, etc. The presence of such solutes affects the evaporation time of respiratory droplets that are sedimenting in air, and thereby the airborne transmission of infections. From solutions of the coupled heat-conduction and water-diffusion equations within the droplet and in the ambient vapor phase, we find that the solute-induced water vapor-pressure reduction considerably slows down the evaporation process and dominates the solute-concentration dependence of the droplet evaporation time. The evaporation-induced solute-concentration gradient near the droplet surface, which is accounted for using a two-stage evaporation model, is found to further intensify the slowing down of the drying process. On the other hand, the presence of solutes is found to reduce evaporation cooling of the droplet, which causes a slight decrease in the evaporation time. Overall, the first two effects are dominant, meaning that the droplet evaporation time increases in the presence of solutes. The solute-concentration dependence of the water diffusivity inside the droplet does not significantly change the evaporation time. Finally, crust formation on the droplet surface is found to increase the final equilibrium size of the droplet. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 40 pages, 8 figures

Fluids at the Nanoscale: from continuum to sub-continuum transport

Authors: Nikita Kavokine, Roland R. Netz, Lydéric Bocquet

Abstract: Nanofluidics has firmly established itself as a new field in fluid mechanics, as novel properties have been shown to emerge in fluids at the nanometric scale. Thanks to recent developments in fabrication technology, artificial nanofluidic systems are now being designed at the scale of biological nanopores. This ultimate step in scale reduction has pushed the development of new experimental techniq… ▽ More Nanofluidics has firmly established itself as a new field in fluid mechanics, as novel properties have been shown to emerge in fluids at the nanometric scale. Thanks to recent developments in fabrication technology, artificial nanofluidic systems are now being designed at the scale of biological nanopores. This ultimate step in scale reduction has pushed the development of new experimental techniques and new theoretical tools, bridging fluid mechanics, statistical mechanics and condensed matter physics. This review is intended as a toolbox for fluids at the nanometre scale. After presenting the basic equations that govern fluid behaviour in the continuum limit, we will show how these equations break down and new properties emerge in molecular scale confinement. △ Less

Submitted 28 November, 2020; originally announced November 2020.

Journal ref: Annual Review of Fluid Mechanics, 53 (2021)

Approach to equilibrium and non-equilibrium stationary distributions of interacting many-particle systems that are coupled to different heat baths

Authors: Roland R. Netz

Abstract: A Hamiltonian-based model of many harmonically interacting massive particles that are subject to linear friction and coupled to heat baths at different temperatures is used to study the dynamic approach to equilibrium and non-equilibrium stationary states. Based on the exactly calculated dynamic approach to the stationary distribution, the functional that governs this approach, which is called the… ▽ More A Hamiltonian-based model of many harmonically interacting massive particles that are subject to linear friction and coupled to heat baths at different temperatures is used to study the dynamic approach to equilibrium and non-equilibrium stationary states. Based on the exactly calculated dynamic approach to the stationary distribution, the functional that governs this approach, which is called the free entropy, is constructed. For the stationary distribution the free entropy becomes maximal and its time derivative is minimal and vanishes. Thus, the free entropy characterizes equilibrium as well as non-equilibrium stationary distributions by their extremal and stability properties. For an equilibrium system, i.e. if all heat baths have the same temperature, the free entropy equals the negative free energy divided by temperature. Using a systematic perturbative scheme for calculating velocity and position correlations in the overdamped massless limit, explicit results for few particles are presented: For two particles localization in position and momentum space is demonstrated in the non-equilibrium stationary state, indicative of a tendency to phase separate. For three elastically interacting particles heat flows from a particle coupled to a cold reservoir to a particle coupled to a warm reservoir if the third reservoir is sufficiently hot. Active particle models can be described in the same general framework, which thereby allows to characterize their entropy production not only in the stationary state but also in the approach to the stationary non-equilibrium state. Finally, the connection to non-equilibrium thermodynamics formulations that include the reservoir entropy production is discussed. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Journal ref: Physical Review E 101, 022120 (2020)

Particle Diffusivity and Free-Energy Profiles in Inhomogeneous Hydrogel Systems from Time-Resolved Penetration Profiles

Authors: Amanuel Wolde-Kidan, Anna Herrmann, Albert Prause, Michael Gradzielski, Rainer Haag, Stephan Block, Roland R. Netz

Abstract: A combined experimental/theoretical method to simultaneously determine diffusivity and free-energy profiles of particles that penetrate into inhomogeneous hydrogel systems is presented. As the only input, arbitrarily normalized concentration profiles from fluorescence intensity data of labeled tracer particles for different penetration times are needed. The method is applied to dextran molecules o… ▽ More A combined experimental/theoretical method to simultaneously determine diffusivity and free-energy profiles of particles that penetrate into inhomogeneous hydrogel systems is presented. As the only input, arbitrarily normalized concentration profiles from fluorescence intensity data of labeled tracer particles for different penetration times are needed. The method is applied to dextran molecules of varying size which penetrate into hydrogels of polyethylene-glycol (PEG) chains with different lengths that are covalently cross-linked by hyperbranched polyglycerol (hPG) hubs. Extracted dextran bulk diffusivities agree well with fluorescence correlation spectroscopy data obtained separately. Scaling laws for dextran diffusivities and free energies inside the hydrogel are identified as a function of the dextran mass. An elastic free-volume model that includes dextran as well as PEG linker flexibility describes the repulsive dextran-hydrogel interaction free energy, which is of steric origin, quantitatively and furthermore suggests that the hydrogel mesh-size distribution is rather broad and particle penetration is dominated by large hydrogel pores. Particle penetration into hydrogels is for steric particle-hydrogel interactions thus suggested to be governed by an elastic size-filtering mechanism that involves the tail of the hydrogel pore-size distribution. △ Less

Submitted 7 September, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Viral air load due to sedimenting and evaporating droplets produced by speaking

Authors: Roland R. Netz

Abstract: The effect of evaporation on droplet sedimentation times is crucial for estimating the risk of infection from virus-containing airborne droplets. For droplet radii in the range 100 nm < R < 60 μm, evaporation can be described in the stagnant-flow approximation and is diffusion limited. Analytical equations are presented for the droplet evaporation rate, the time-dependent droplet size and the sedi… ▽ More The effect of evaporation on droplet sedimentation times is crucial for estimating the risk of infection from virus-containing airborne droplets. For droplet radii in the range 100 nm < R < 60 μm, evaporation can be described in the stagnant-flow approximation and is diffusion limited. Analytical equations are presented for the droplet evaporation rate, the time-dependent droplet size and the sedimentation time, including the significant effect of evaporation cooling. Evaporation makes the time for large droplets to sediment much longer and thus significantly increases the viral air load. Using recent estimates for SARS-CoV-2 concentrations in sputum and droplet production rates while speaking, a single infected person that constantly speaks without a mouth cover produces a total air load of more than 10^4 virions. In a mid-size closed room, this leads to a viral inhalation frequency of at least 2.5 per minute. Low relative humidity, as encountered inside buildings in winter and in airliners, speeds evaporation and thus keeps initially larger droplets suspended in air. △ Less

Submitted 25 May, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: New version contains estimate of number of virions sedimenting in air in steady state and effect of air-exchange rate in closed room

A rigorous derivation and energetics of a wave equation with fractional damping

Authors: Alexander Mielke, Roland R. Netz, Sina Zendehroud

Abstract: We consider a linear system that consists of a linear wave equation on a horizontal hypersurface and a parabolic equation in the half space below. The model describes longitudinal elastic waves in organic monolayers at the water-air interface, which is an experimental setup that is relevant for understanding wave propagation in biological membranes. We study the scaling regime where the relevant h… ▽ More We consider a linear system that consists of a linear wave equation on a horizontal hypersurface and a parabolic equation in the half space below. The model describes longitudinal elastic waves in organic monolayers at the water-air interface, which is an experimental setup that is relevant for understanding wave propagation in biological membranes. We study the scaling regime where the relevant horizontal length scale is much larger than the vertical length scale and provide a rigorous limit leading to a fractionally-damped wave equation for the membrane. We provide the associated existence results via linear semigroup theory and show convergence of the solutions in the scaling limit. Moreover, based on the energy-dissipation structure for the full model, we derive a natural energy and a natural dissipation function for the fractionally-damped wave equation with a time derivative of order 3/2 △ Less

Submitted 24 April, 2020; originally announced April 2020.

MSC Class: 35Q74

Emergence of molecular friction in liquids: bridging between the atomistic and hydrodynamic pictures

Authors: Arthur V. Straube, Bartosz G. Kowalik, Roland R. Netz, Felix Höfling

Abstract: Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data f… ▽ More Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with rigorous theoretical arguments, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids; as a consequence, its origin is non-local in time. Concomitantly, the molecules experience Brownian forces that display persistent correlations and long-lasting memory. A critical test of the generalised Stokes-Einstein relation, mapping the friction of single molecules to the viscoelastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield novel insights into vitrification mechanisms and the intriguing mechanical properties of soft materials. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Journal ref: Commun. Phys. 3, 126 (2020)

Energy Transfer within the Hydrogen Bonding Network of Water Following Resonant Terahertz Excitation

Authors: Hossam Elgabarty, Tobias Kampfrath, Douwe Jan Bonthuis, Vasileios Balos, Naveen Kumar Kaliannan, Philip Loche, Roland R. Netz, Martin Wolf, Thomas D. Kühne, Mohsen Sajadi

Abstract: Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen-bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle tera… ▽ More Energy dissipation in water is very fast and more efficient than in many other liquids. This behavior is commonly attributed to the intermolecular interactions associated with hydrogen bonding. Here, we investigate the dynamic energy flow in the hydrogen-bond network of liquid water by a pump-probe experiment. We resonantly excite intermolecular degrees of freedom with ultrashort single-cycle terahertz pulses and monitor its Raman response. By using ultrathin sample-cell windows, a background-free bipolar signal whose tail relaxes mono-exponentially is obtained. The relaxation is attributed to the molecular translational motions, using complementary experiments, force-field and ab initio molecular dynamics simulations. They reveal an initial coupling of the terahertz electric field to the molecular rotational degrees of freedom whose energy is rapidly transferred, within the excitation pulse duration, to the restricted-translational motion of neighboring molecules. This rapid energy transfer may be rationalized by the strong anharmonicity of the intermolecular interactions. △ Less

Submitted 24 March, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Molecular interpretation of the non-Newtonian viscoelastic behavior of liquid water at high frequencies

Authors: Julius C. F. Schulz, Alexander Schlaich, Matthias Heyden, Roland R. Netz, Julian Kappler

Abstract: Using classical as well as ab-initio molecular dynamics simulations, we calculate the frequency-dependent shear viscosity of pure water and water-glycerol mixtures. In agreement with recent experiments, we find deviations from Newtonian-fluid behavior in the THz regime. Based on an extension of the Maxwell model, we introduce a viscoelastic model to describe the observed viscosity spectrum of pure… ▽ More Using classical as well as ab-initio molecular dynamics simulations, we calculate the frequency-dependent shear viscosity of pure water and water-glycerol mixtures. In agreement with recent experiments, we find deviations from Newtonian-fluid behavior in the THz regime. Based on an extension of the Maxwell model, we introduce a viscoelastic model to describe the observed viscosity spectrum of pure water. We find four relaxation modes in the spectrum which we attribute to i) hydrogen-bond network topology changes, ii) hydrogen-bond stretch vibrations of water pairs, iii) collective vibrations of water molecule triplets, and iv) librational excitations of individual water molecules. Our model quantitatively describes the viscoelastic response of liquid water on short timescales, where the hydrodynamic description via a Newtonian-fluid model breaks down. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: Supplemental material is provided as ancillary file

Journal ref: Phys. Rev. Fluids 5, 103301 (2020)

Exploring the Absorption Spectrum of Simulated Water from MHz to the Infrared

Authors: Shane Carlson, Florian N. Brünig, Philip Loche, Douwe Jan Bonthuis, Roland R. Netz

Abstract: Absorption spectra of liquid water at 300 K are calculated from both classical and density functional theory molecular dynamics simulation data, which together span from 1 MHz to hundreds of THz, agreeing well with experimental data qualitativley and quantitavely over the entire range, including the IR modes, the microwave peak, and the intermediate THz bands. The spectra are decomposed into singl… ▽ More Absorption spectra of liquid water at 300 K are calculated from both classical and density functional theory molecular dynamics simulation data, which together span from 1 MHz to hundreds of THz, agreeing well with experimental data qualitativley and quantitavely over the entire range, including the IR modes, the microwave peak, and the intermediate THz bands. The spectra are decomposed into single-molecular and collective components as well as into components due to molecular reorientations and changes in induced intramolecular dipole moments. These decompositions shed new light on the motions underlying the librational and translational (hydrogen-bond stretch) bands at 20 and 5 THz respectively: interactions between donor protons and acceptor lone pair electrons are shown to be important for the line shape in both librational and translational regimes, and in- and out-of-phase librational dimer modes are observed and explored. △ Less

Submitted 16 March, 2020; originally announced March 2020.

Comments: Main text: 7 pages with 4 figures. Supplement: 8 pages with 5 figures. Submitted to Nature Communications on February 29, 2020

Non-Markovian barrier crossing with two-time-scale memory is dominated by the faster memory component

Authors: Julian Kappler, Victor B. Hinrichsen, Roland R. Netz

Abstract: We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times $τ_1$ and $τ_2$. Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier crossing time becomes independent of the longer memory… ▽ More We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times $τ_1$ and $τ_2$. Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier crossing time becomes independent of the longer memory time if at least one of the two memory times is larger than the intrinsic diffusion time. When we associate memory effects with coupled degrees of freedom that are orthogonal to a one-dimensional reaction coordinate, this counterintuitive result shows that the faster orthogonal degrees of freedom dominate barrier-crossing kinetics in the non-Markovian limit and that the slower orthogonal degrees become negligible, quite contrary to the standard time-scale separation assumption and with important consequences for the proper setup of coarse-graining procedures in the non-Markovian case. By asymptotic matching and symmetry arguments, we construct a crossover formula for the barrier crossing time that is valid for general multi-exponential memory kernels. This formula can be used to estimate barrier-crossing times for general memory functions for high friction, i.e. in the overdamped regime, as well as for low friction, i.e. in the inertial regime. Typical examples where our results are important include protein folding in the high-friction limit and chemical reactions such as proton-transfer reactions in the low-friction limit. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Cold denaturation of RNA secondary structures with loop entropy and quenched disorder

Authors: Flavio Iannelli, Yevgeni Mamasakhlisov, Roland R. Netz

Abstract: We study the folding of RNA secondary structures with quenched sequence randomness by means of the constrained annealing method. A thermodynamic phase transition is induced by including the conformational weight of loop structures. In addition to the expected melting at high temperature, a cold melting transition appears. Our results suggest that the cold denaturation of RNA found experimentally i… ▽ More We study the folding of RNA secondary structures with quenched sequence randomness by means of the constrained annealing method. A thermodynamic phase transition is induced by including the conformational weight of loop structures. In addition to the expected melting at high temperature, a cold melting transition appears. Our results suggest that the cold denaturation of RNA found experimentally is, in fact, a continuous phase transition triggered by quenched sequence disorder. We calculate both hot and cold melting critical temperatures for the competing energy scenario between favorable and unfavorable base pairs and present a phase diagram as a function of the loop exponent and temperature. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. E 101, 012502 (2020)

Nonlinear fractional waves at elastic interfaces

Authors: Julian Kappler, Shamit Shrivastava, Matthias F. Schneider, Roland R. Netz

Abstract: We derive the nonlinear fractional surface wave equation that governs compression waves at an interface that is coupled to a viscous bulk medium. The fractional character of the differential equation comes from the fact that the effective thickness of the bulk layer that is coupled to the interface is frequency dependent. The nonlinearity arises from the nonlinear dependence of the interface compr… ▽ More We derive the nonlinear fractional surface wave equation that governs compression waves at an interface that is coupled to a viscous bulk medium. The fractional character of the differential equation comes from the fact that the effective thickness of the bulk layer that is coupled to the interface is frequency dependent. The nonlinearity arises from the nonlinear dependence of the interface compressibility on the local compression, which is obtained from experimental measurements and reflects a phase transition at the interface. Numerical solutions of our nonlinear fractional theory reproduce several experimental key features of surface waves in phospholipid monolayers at the air-water interface without freely adjustable fitting parameters. In particular, the propagation length of the surface wave abruptly increases at a threshold excitation amplitude. The wave velocity is found to be of the order of 40 cm/s both in experiments and theory and slightly increases as a function of the excitation amplitude. Nonlinear acoustic switching effects in membranes are thus shown to arise purely based on intrinsic membrane properties, namely the presence of compressibility nonlinearities that accompany phase transitions at the interface. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: 12 pages, 7 figures; supplemental information provided as ancillary file

Journal ref: Phys. Rev. Fluids 2, 114804 (2017)

Shear-induced unfolding and enzymatic cleavage of full-length VWF multimers

Authors: Svenja Lippok, Matthias Radtke, Tobias Obser, Lars Kleemeier, Reinhard Schneppenheim, Ulrich Budde, Roland R. Netz, Joachim O. Rädler

Abstract: Proteolysis of the multimeric blood coagulation protein von Willebrand Factor (VWF) by ADAMTS13 is crucial for prevention of microvascular thrombosis. ADAMTS13 cleaves VWF within the mechanosensitive A2 domain, which is believed to open under shear flow. Here, we combine Fluorescence Correlation Spectroscopy (FCS) and a microfluidic shear cell to monitor real-time kinetics of full-length VWF prote… ▽ More Proteolysis of the multimeric blood coagulation protein von Willebrand Factor (VWF) by ADAMTS13 is crucial for prevention of microvascular thrombosis. ADAMTS13 cleaves VWF within the mechanosensitive A2 domain, which is believed to open under shear flow. Here, we combine Fluorescence Correlation Spectroscopy (FCS) and a microfluidic shear cell to monitor real-time kinetics of full-length VWF proteolysis as a function of shear stress. For comparison, we also measure the Michaelis-Menten kinetics of ADAMTS13 cleavage of wild-type VWF in the absence of shear but partially denaturing conditions. Under shear, ADAMTS13 activity on full-length VWF arises without denaturing agent as evidenced by FCS and gel-based multimer analysis. In agreement with Brownian hydrodynamics simulations, we find a sigmoidal increase of the enzymatic rate as a function of shear at a threshold shear rate 5522/s. The same flow-rate dependence of ADAMTS13 activity we also observe in blood plasma, which is relevant to predict hemostatic dysfunction. △ Less

Submitted 16 December, 2015; originally announced December 2015.

The mean shape of transition and first-passage paths

Authors: Won Kyu Kim, Roland R. Netz

Abstract: We calculate the mean shape of transition paths and first-passage paths based on the one-dimensional Fokker-Planck equation in an arbitrary free energy landscape including a general inhomogeneous diffusivity profile. The transition path ensemble is the collection of all paths that do not revisit the start position $x_A$ and that terminate when first reaching the final position $x_B$. In contrast,… ▽ More We calculate the mean shape of transition paths and first-passage paths based on the one-dimensional Fokker-Planck equation in an arbitrary free energy landscape including a general inhomogeneous diffusivity profile. The transition path ensemble is the collection of all paths that do not revisit the start position $x_A$ and that terminate when first reaching the final position $x_B$. In contrast, a first-passage path can revisit but not cross its start position $x_A$ before it terminates at $x_B$. Our theoretical framework employs the forward and backward Fokker-Planck equations as well as first-passage, passage, last-passage and transition-path time distributions, for which we derive the defining integral equations. We show that the mean time at which the transition path ensemble visits an intermediate position $x$ is equivalent to the mean first-passage time of reaching the starting position $x_A$ from $x$ without ever visiting $x_B$. The mean shape of first-passage paths is related to the mean shape of transition paths by a constant time shift. Since for large barrier height $U$ the mean first-passage time scales exponentially in $U$ while the mean transition path time scales linearly inversely in $U$, the time shift between first-passage and transition path shapes is substantial. We present explicit examples of transition path shapes for linear and harmonic potentials and illustrate our findings by trajectories generated from Brownian dynamics simulations. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Journal ref: J. Chem. Phys. 143, 224108 (2015)

Global analysis of the ground-state wrapping conformation of a charged polymer on an oppositely charged nano-sphere

Authors: Hoda Boroudjerdi, Ali Naji, Roland R. Netz

Abstract: We investigate the wrapping conformations of a strongly adsorbed polymer chain on an oppositely charged nano-sphere by employing a reduced (dimensionless) representation of a primitive chain-sphere model. This enables us to determine the global phase behavior of the chain conformation in a wide range of values for the system parameters including the chain contour length, its linear charge density… ▽ More We investigate the wrapping conformations of a strongly adsorbed polymer chain on an oppositely charged nano-sphere by employing a reduced (dimensionless) representation of a primitive chain-sphere model. This enables us to determine the global phase behavior of the chain conformation in a wide range of values for the system parameters including the chain contour length, its linear charge density and persistence length as well as the nano-sphere charge and radius, and also the salt concentration in the bathing solution. The phase behavior of a charged chain-sphere complex can be described in terms of a few distinct conformational symmetry classes (phases) separated by continuous or discontinuous transition lines which are determined by means of appropriately defined (order) parameters. Our results can be applied to a wide class of strongly coupled polymer-sphere complexes including, for instances, complexes that comprise a mechanically flexible or semiflexible polymer chain or an extremely short or long chain and, as a special case, include the previously studied example of DNA-histone complexes. △ Less

Submitted 21 November, 2013; originally announced November 2013.

Comments: 16 pages, 12 figures

Journal ref: Eur. Phys. J. E 37, 21 (2014)

Elasticity of DNA and the effect of Dendrimer Binding

Authors: Santosh Mogurampelly, Bidisha Nandy, Roland R. Netz, Prabal K. Maiti

Abstract: Negatively charged DNA can be compacted by positively charged dendrimers and the degree of compaction is a delicate balance between the strength of the electrostatic interaction and the elasticity of DNA. We report various elastic properties of short double stranded DNA (dsDNA) and the effect of dendrimer binding using fully atomistic molecular dynamics and numerical simulations. In equilibrium at… ▽ More Negatively charged DNA can be compacted by positively charged dendrimers and the degree of compaction is a delicate balance between the strength of the electrostatic interaction and the elasticity of DNA. We report various elastic properties of short double stranded DNA (dsDNA) and the effect of dendrimer binding using fully atomistic molecular dynamics and numerical simulations. In equilibrium at room temperature, the contour length distribution P(L) and end-to-end distance distribution P(R) are nearly Gaussian, the former gives an estimate of the stretch modulus γ_1 of dsDNA in quantitative agreement with the literature value. The bend angle distribution P(θ) of the dsDNA also has a Gaussian form and allows to extract a persistence length, L_p of 43 nm. When the dsDNA is compacted by positively charged dendrimer, the stretch modulus stays invariant but the effective bending rigidity estimated from the end-to-end distance distribution decreases dramatically due to backbone charge neutralization of dsDNA by dendrimer. We support our observations with numerical solutions of the worm-like-chain (WLC) model as well as using non-equilibrium dsDNA stretching simulations. These results are helpful in understanding the dsDNA elasticity at short length scales as well as how the elasticity is modulated when dsDNA binds to a charged object such as a dendrimer or protein. △ Less

Submitted 24 July, 2013; originally announced July 2013.

Comments: 21 pages, 5 figures

Journal ref: Eur. Phys. J. E (2013) 36: 68

Hydrodynamic effects in driven soft matter

Authors: Manoel Manghi, Xaver Schlagberger, Yong-Woon Kim, Roland R. Netz

Abstract: Recent theoretical works exploring the hydrodynamics of soft material in non-equilibrium situations are reviewed. We discuss the role of hydrodynamic interactions for three different systems: i) the deformation and orientation of sedimenting semiflexible polymers, ii) the propulsion and force-rectification with a nano-machine realized by a rotating elastic rod, and iii) the deformation of a brush… ▽ More Recent theoretical works exploring the hydrodynamics of soft material in non-equilibrium situations are reviewed. We discuss the role of hydrodynamic interactions for three different systems: i) the deformation and orientation of sedimenting semiflexible polymers, ii) the propulsion and force-rectification with a nano-machine realized by a rotating elastic rod, and iii) the deformation of a brush made of grafted semiflexible polymers in shear flows. In all these examples deformable polymers are subject to various hydrodynamic flows and hydrodynamic interactions. Perfect stiff nano-cylinders are known to show no orientational effects as they sediment through a viscous fluid, but it is the coupling between elasticity and hydrodynamic torques that leads to an orientation perpendicular to the direction of sedimentation. Likewise, a rotating stiff rod does not lead to a net propulsion in the Stokes limit, but if bending is allowed an effective thrust develops whose strength and direction is independent of the sense of rotation and thus acts as a rectification device. Lastly, surface-anchored polymers are deformed by shear flows, which modifies the effective hydrodynamic boundary condition in a non-linear fashion. All these results are obtained with hydrodynamic Brownian dynamics simulation techniques, as appropriate for dilute systems. Scaling analyses are presented when possible. The common theme is the interaction between elasticity of soft matter and hydrodynamics, which can lead to qualitatively new effects. △ Less

Submitted 7 March, 2012; originally announced March 2012.

Comments: review article, 17 pages, 12 figures, pdftex

Journal ref: Soft Matter, 2006, vol. 2, pp. 653-668

Friction contribution to water-bond breakage kinetics

Authors: Yann von Hansen, Felix Sedlmeier, Michael Hinczewski, Roland R. Netz

Abstract: Based on the trajectories of the separation between water molecule pairs from MD simulations, we investigate the bond breakage dynamics in bulk water. From the spectrum of mean first-passage times, the Fokker-Planck equation allows us to derive the diffusivity profile along the separation coordinate and thus to unambiguously disentangle the effects of free-energy and local friction on the separati… ▽ More Based on the trajectories of the separation between water molecule pairs from MD simulations, we investigate the bond breakage dynamics in bulk water. From the spectrum of mean first-passage times, the Fokker-Planck equation allows us to derive the diffusivity profile along the separation coordinate and thus to unambiguously disentangle the effects of free-energy and local friction on the separation kinetics. For tightly coordinated water the friction is six times higher than in bulk, which can be interpreted in terms of a dominant reaction path that involves additional orthogonal coordinates. △ Less

Submitted 2 November, 2011; originally announced November 2011.

Comments: 11 pages, 10 figures

Journal ref: Phys. Rev. E 84, 051501 (2011)

Conformational dynamics and internal friction in homo-polymer globules: equilibrium vs. non-equilibrium simulations

Authors: Thomas R. Einert, Charles E. Sing, Alfredo Alexander-Katz, Roland R. Netz

Abstract: We study the conformational dynamics within homo-polymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength ε and the globule size NG is observed. We find two distinct dynamical regimes: a liquid- like regime (for ε < εs) with fast internal dynamics and a solid-like regime (for ε > εs) with slow inte… ▽ More We study the conformational dynamics within homo-polymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength ε and the globule size NG is observed. We find two distinct dynamical regimes: a liquid- like regime (for ε < εs) with fast internal dynamics and a solid-like regime (for ε > εs) with slow internal dynamics. The cohesion strength εs of this freezing transition depends on NG. Equilibrium simulations, where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompassing low enough velocities so that the linear-response, viscous regime is reached). From both simulation protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal friction increases continuously with ε and scales extensive in NG. This suggests an internal friction scenario where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the globular structure. △ Less

Submitted 25 September, 2011; originally announced September 2011.

Comments: 19 pages, 18 figures

Deconvolution of dynamic mechanical networks

Authors: Michael Hinczewski, Yann von Hansen, Roland R. Netz

Abstract: Time-resolved single-molecule biophysical experiments yield data that contain a wealth of dynamic information, in addition to the equilibrium distributions derived from histograms of the time series. In typical force spectroscopic setups the molecule is connected via linkers to a read-out device, forming a mechanically coupled dynamic network. Deconvolution of equilibrium distributions, filtering… ▽ More Time-resolved single-molecule biophysical experiments yield data that contain a wealth of dynamic information, in addition to the equilibrium distributions derived from histograms of the time series. In typical force spectroscopic setups the molecule is connected via linkers to a read-out device, forming a mechanically coupled dynamic network. Deconvolution of equilibrium distributions, filtering out the influence of the linkers, is a straightforward and common practice. We have developed an analogous dynamic deconvolution theory for the more challenging task of extracting kinetic properties of individual components in networks of arbitrary complexity and topology. Our method determines the intrinsic linear response functions of a given molecule in the network, describing the power spectrum of conformational fluctuations. The practicality of our approach is demonstrated for the particular case of a protein linked via DNA handles to two optically trapped beads at constant stretching force, which we mimic through Brownian dynamics simulations. Each well in the protein free energy landscape (corresponding to folded, unfolded, or possibly intermediate states) will have its own characteristic equilibrium fluctuations. The associated linear response function is rich in physical content, since it depends both on the shape of the well and its diffusivity---a measure of the internal friction arising from such processes like the transient breaking and reformation of bonds in the protein structure. Starting from the autocorrelation functions of the equilibrium bead fluctuations measured in this force clamp setup, we show how an experimentalist can accurately extract the state-dependent protein diffusivity using a straightforward two-step procedure. △ Less

Submitted 13 July, 2011; originally announced July 2011.

Comments: 9 pages, 3 figures + supplementary material 14 pages, 4 figures

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 107, 21493 (2010)

Secondary structure formation of homopolymeric single-stranded nucleic acids including force and loop entropy: implications for DNA hybridization

Authors: Thomas R. Einert, Henri Orland, Roland R. Netz

Abstract: Loops are essential secondary structure elements in folded DNA and RNA molecules and proliferate close to the melting transition. Using a theory for nucleic acid secondary structures that accounts for the logarithmic entropy c ln m for a loop of length m, we study homopolymeric single-stranded nucleic acid chains under external force and varying temperature. In the thermodynamic limit of a long st… ▽ More Loops are essential secondary structure elements in folded DNA and RNA molecules and proliferate close to the melting transition. Using a theory for nucleic acid secondary structures that accounts for the logarithmic entropy c ln m for a loop of length m, we study homopolymeric single-stranded nucleic acid chains under external force and varying temperature. In the thermodynamic limit of a long strand, the chain displays a phase transition between a low temperature / low force compact (folded) structure and a high temperature / high force molten (unfolded) structure. The influence of c on phase diagrams, critical exponents, melting, and force extension curves is derived analytically. For vanishing pulling force, only for the limited range of loop exponents 2 < c < 2.479 a melting transition is possible; for c <= 2 the chain is always in the folded phase and for 2.479 < c always in the unfolded phase. A force induced melting transition with singular behavior is possible for all loop exponents c < 2.479 and can be observed experimentally by single molecule force spectroscopy. These findings have implications for the hybridization or denaturation of double stranded nucleic acids. The Poland-Scheraga model for nucleic acid duplex melting does not allow base pairing between nucleotides on the same strand in denatured regions of the double strand. If the sequence allows these intra-strand base pairs, we show that for a realistic loop exponent c ~ 2.1 pronounced secondary structures appear inside the single strands. This leads to a lower melting temperature of the duplex than predicted by the Poland-Scheraga model. Further, these secondary structures renormalize the effective loop exponent c^, which characterizes the weight of a denatured region of the double strand, and thus affect universal aspects of the duplex melting transition. △ Less

Submitted 28 April, 2011; originally announced April 2011.

Comments: 19 pages, 14 figures, supplementary material

Theory for RNA folding, stretching, and melting including loops and salt

Authors: Thomas R. Einert, Roland R. Netz

Abstract: Secondary structure formation of nucleic acids strongly depends on salt concentration and temperature. We develop a theory for RNA folding that correctly accounts for sequence effects, the entropic contributions associated with loop formation, and salt effects. Using an iterative expression for the partition function that neglects pseudoknots, we calculate folding free energies and minimum free en… ▽ More Secondary structure formation of nucleic acids strongly depends on salt concentration and temperature. We develop a theory for RNA folding that correctly accounts for sequence effects, the entropic contributions associated with loop formation, and salt effects. Using an iterative expression for the partition function that neglects pseudoknots, we calculate folding free energies and minimum free energy configurations based on the experimentally derived base pairing free energies. The configurational entropy of loop formation is modeled by the asymptotic expression -c ln m, where m is the length of the loop and c the loop exponent, which is an adjustable constant. Salt effects enter in two ways: first, we derive salt induced modifications of the free energy parameters for describing base pairing and, second, we include the electrostatic free energy for loop formation. Both effects are modeled on the Debye-Hueckel level including counterion condensation. We validate our theory for two different RNA sequences: For tRNA-phe, the resultant heat capacity curves for thermal denaturation at various salt concentrations accurately reproduce experimental results. For the P5ab RNA hairpin, we derive the global phase diagram in the three-dimensional space spanned by temperature, stretching force, and salt concentration and obtain good agreement with the experimentally determined critical unfolding force. We show that for a proper description of RNA melting and stretching, both salt and loop entropy effects are needed. △ Less

Submitted 20 April, 2011; originally announced April 2011.

Comments: 12 pages, 9 figures, accepted for publication in Biophysical Journal

Salt Modulated Structure of Polyelectrolyte-Macroion Complex Fibers

Authors: Hoda Boroudjerdi, Ali Naji, Roland R. Netz

Abstract: The structure and stability of strongly charged complex fibers formed by complexation of a single long semi-flexible polyelectrolyte (PE) chain and many oppositely charged spherical macroions are investigated numerically at the ground-state level using a chain-sphere cell model. The model takes into account chain elasticity as well as electrostatic interactions between charged spheres and chain se… ▽ More The structure and stability of strongly charged complex fibers formed by complexation of a single long semi-flexible polyelectrolyte (PE) chain and many oppositely charged spherical macroions are investigated numerically at the ground-state level using a chain-sphere cell model. The model takes into account chain elasticity as well as electrostatic interactions between charged spheres and chain segments. Using a numerical optimization method based on a periodically repeated unit cell, we obtain fiber configurations that minimize the total energy. The optimal configurations exhibit a variety of helical structures for the arrangement of macroions including zig-zag, solenoidal and beads-on-a-string patterns. These structures are determined by a competition between attraction between spheres and the PE chain (which favors chain wrapping around the spheres), chain bending and electrostatic repulsion between chain segments (which favor unwrapping of the chain), and the interactions between neighboring sphere-chain complexes which can be attractive or repulsive depending on the system parameters such as medium salt concentration, macroion charge and chain length per macroion (linker size). At about physiological salt concentration, dense zig-zag patterns are found to be energetically most stable when parameters appropriate for the DNA-histone system in a chromatin fiber are adopted. In fact, the predicted fiber diameter in this regime is found to be around 30nm, which appears to agree with the thickness observed in in vitro experiments on chromatin. We also find a macroion density of 5-6 per 11nm which agrees with the zig-zag or cross-linker models of chromatin. Since our study deals primarily with a generic model, these findings suggest that chromatin-like structures should also be observable for PE-macroion complexes formed in solutions of DNA and synthetic nano-colloids of opposite charge. △ Less

Submitted 30 January, 2011; originally announced January 2011.

Comments: 18 pages, 15 figures

Journal ref: Eur. Phys. J. E 34, 72 (2011)

Single polymer adsorption in shear: flattening versus hydrodynamic lift and corrugation effects

Authors: Andreas Serr, Christian Sendner, Florian Mueller, Thomas R. Einert, Roland R. Netz

Abstract: The adsorption of a single polymer to a flat surface in shear is investigated using Brownian hydrodynamics simulations and scaling arguments. Competing effects are disentangled: in the absence of hydrodynamic interactions, shear drag flattens the chain and thus enhances adsorption. Hydrodynamic lift on the other hand gives rise to long-ranged repulsion from the surface which preempts the surface-a… ▽ More The adsorption of a single polymer to a flat surface in shear is investigated using Brownian hydrodynamics simulations and scaling arguments. Competing effects are disentangled: in the absence of hydrodynamic interactions, shear drag flattens the chain and thus enhances adsorption. Hydrodynamic lift on the other hand gives rise to long-ranged repulsion from the surface which preempts the surface-adsorbed state via a discontinuous desorption transition, in agreement with theoretical arguments. Chain flattening is dominated by hydrodynamic lift, so overall, shear flow weakens the adsorption of flexible polymers. Surface friction due to small-wavelength surface potential corrugations is argued to weaken the surface attraction as well. △ Less

Submitted 25 October, 2010; originally announced October 2010.

Comments: 6 pages, 4 figures

Exotic Electrostatics: Unusual Features of Electrostatic Interactions between Macroions

Authors: Ali Naji, Matej Kanduc, Roland R. Netz, Rudolf Podgornik

Abstract: We present an overview of our understanding of electrostatic interactions between charged macromolecular surfaces mediated by mobile counter- and coions. The dichotomy between the weak and the strong coupling regimes is described in detail and the way they engender repulsive and attractive interactions between nominally equally charged macroions. We also introduce the concept of dressed counterion… ▽ More We present an overview of our understanding of electrostatic interactions between charged macromolecular surfaces mediated by mobile counter- and coions. The dichotomy between the weak and the strong coupling regimes is described in detail and the way they engender repulsive and attractive interactions between nominally equally charged macroions. We also introduce the concept of dressed counterions in the case of many-component Coulomb fluids that are partially weakly and partially strongly coupled to local electrostatic fields leading to non-monotonic interactions between equally charged macroions. The effect of quenched surface charge disorder on the counterion-mediated electrostatic interactions is analyzed within the same conceptual framework and shown to lead to unexpected and extraordinary electrostatic interactions between randomly charged surfaces with equal mean surface charge densities or even between effectively neutral macroion surfaces. As a result, these recent developments challenge some cherished notions of pop culture. △ Less

Submitted 2 August, 2010; originally announced August 2010.

Comments: 18 pages, 5 figures

How the diffusivity profile reduces the arbitrariness of protein folding free energies

Authors: Michael Hinczewski, Yann von Hansen, Joachim Dzubiella, Roland R. Netz

Abstract: The concept of a protein diffusing in its free energy folding landscape has been fruitful for both theory and experiment. Yet the choice of the reaction coordinate (RC) introduces an undesirable degree of arbitrariness into the problem. We analyze extensive simulation data of an alpha-helix in explicit water solvent as it stochastically folds and unfolds. The free energy profiles for different R… ▽ More The concept of a protein diffusing in its free energy folding landscape has been fruitful for both theory and experiment. Yet the choice of the reaction coordinate (RC) introduces an undesirable degree of arbitrariness into the problem. We analyze extensive simulation data of an alpha-helix in explicit water solvent as it stochastically folds and unfolds. The free energy profiles for different RCs exhibit significant variation, some having an activation barrier, others not. We show that this variation has little effect on the predicted folding kinetics if the diffusivity profiles are properly taken into account. This kinetic quasi-universality is rationalized by an RC rescaling, which, due to the reparameterization invariance of the Fokker-Planck equation, allows the combination of free energy and diffusivity effects into a single function, the rescaled free energy profile. This rescaled free energy indeed shows less variation among different RCs than the bare free energy and diffusivity profiles separately do, if we properly distinguish between RCs that contain knowledge of the native state and those that are purely geometric in nature. Our method for extracting diffusivity profiles is easily applied to experimental single molecule time series data and might help to reconcile conflicts that arise when comparing results from different experimental probes for the same protein. △ Less

Submitted 12 May, 2010; originally announced May 2010.

Comments: 22 pages, 9 figures + supplementary material (8 pages, 5 figures); to appear in J. Chem. Phys.