BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Joachim Dzubiella

Joachim Dzubiella appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.chem-phBiological Physicscond-mat.stat-mechBiomoleculescond-mat.mtrl-scicond-mat.mes-hallphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

36works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Nonequilibrium free energy during polymer chain growth

During fast diffusion-influenced polymerization, nonequilibrium behavior of the polymer chains and the surrounding reactive monomers has been reported recently. Based on the laws of thermodynamics, the emerging nonequilibrium structures should be characterisable by some "extra free energy" (excess over the equilibrium Helmholtz free energy). Here, we study the nonequilibrium thermodynamics of chain-growth polymerization of ideal chains in a dispersion of free reactive monomers, using off-lattice, reactive Brownian Dynamics (R-BD) computer simulations in conjunction with approximative statistical mechanics and relative entropy (Gibbs-Shannon and Kullback-Leibler) concepts. In case of fast growing polymers, we indeed report increased nonequilibrium free energies of several kBT compared to equilibrium and near-equilibrium, slowly growing chains. Interestingly, the extra free energy is a non-monotonic function of the degree of polymerization and thus also of time. Our decomposition of the thermodynamic contributions shows that the initial dominant extra free energy is stored in the nonequilibrium inhomogeneous density profiles of the free monomer gas (showing density depletion and wakes) in the vicinity of the active center at the propagating polymer end. At later stages of the polymerization process, we report significant extra contributions stored in the nonequilibrium polymer conformations. Finally, our study implies a nontrivial relaxation kinetics and restoring of the extra free energy during the equilibration process after polymerization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Active interaction switching controls the dynamic heterogeneity of soft colloidal dispersions

We employ Reactive Dynamical Density Functional Theory (R-DDFT) and Reactive Brownian Dynamics (R-BD) simulations to investigate the dynamics of a suspension of active soft Gaussian colloids with binary interaction switching, i.e., a one-component colloidal system in which every particle stochastically switches at predefined rates between two interaction states with different mobility. Using R-DDFT we extend a theory previously developed to access the dynamics of inhomogeneous liquids to study the influence of the switching activity on the self and distinct part of the Van Hove function in bulk solution, and determine the corresponding mean squared displacement of the switching particles. Our results demonstrate that, even though the average diffusion coefficient is not affected by the switching activity, it significantly modifies the non-equilibrium dynamics and diffusion coefficients of the individual particles, leading to a crossover from short to long times, with a regime for intermediate times showing anomalous diffusion. In addition, the self-part of the van Hove function has a Gaussian form at short and long times, but becomes non-Gaussian at intermediates ones, having a crossover between short and large displacements. The corresponding self-intermediate scattering function shows the two-step relaxation patters typically observed in soft materials with heterogeneous dynamics such as glasses and gels. We also introduce a phenomenological Continuous Time Random Walk (CTRW) theory to understand the heterogeneous diffusion of this system. R-DDFT results are in excellent agreement with R-BD simulations and the analytical predictions of CTRW theory, thus confirming that R-DDFT constitutes a powerful method to investigate not only the structure and phase behavior, but also the dynamical properties of non-equilibrium active switching colloidal suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nanochannels and nanodroplets in polymer membranes controlling ionic transport

Polymer materials with low water uptake exhibit a highly heterogeneous interior, characterized by water clusters in the form of nanodroplets and nanochannels. Here, based on our recent insights from computer simulations, we argue that water cluster structure has large implications for ionic transport and selective permeability in polymer membranes. Importantly, we demonstrate that the two key quantities for transport, the ion diffusion and the solvation free energy inside the polymer, are extremely sensitive to molecular details of the water clusters. In particular, we highlight the significance of water droplet interface potentials and the nature of hopping diffusion through transient water channels. These mechanisms can be harvested and fine-tuned to optimize selectivity in ionic transport in a wide range of applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Transport in polymer membranes beyond linear response: Controlling permselectivity by the driving force

In the popular solution-diffusion picture, the membrane permeability is defined as the product of the partition ratio and the diffusivity of penetrating solutes inside the membrane in the linear response regime, i.e., in equilibrium. However, of practical importance is the penetrants' flux across the membrane driven by external forces. Here, we study nonequilibrium membrane permeation orchestrated by a uniform external driving field using molecular computer simulations and continuum (Smoluchowski) theory in the stationary state. In the simulations, we explicitly resolve the penetrants' transport across a finite monomer-resolved polymer network, addressing one-component penetrant systems and mixtures. We introduce and discuss possible definitions of nonequilibrium, force-dependent permeability, representing `system' and `membrane' permeability. In particular, we present for the first time a definition of the differential permeability response to the force. We demonstrate that the latter turns out to be significantly nonlinear for low-permeable systems, leading to a high amount of selectiveness in permeability, called `permselectivity', and is tunable by the driving force. Our continuum-level analytical solutions exhibit remarkable qualitative agreement with the penetrant- and polymer-resolved simulations, thereby allowing us to characterize the underlying mechanism of permeabilities and steady-state transport beyond the linear response level.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Affinity, Kinetics, and Pathways of Anisotropic Ligands Binding to Hydrophobic Model Pockets

Using explicit-water molecular dynamics (MD) simulations of a generic pocket-ligand model we investigate how chemical and shape anisotropy of small ligands influences the affinities, kinetic rates and pathways for their association to hydrophobic binding sites. In particular, we investigate aromatic compounds, all of similar molecular size, but distinct by various hydrophilic or hydrophobic residues. We demonstrate that the most hydrophobic sections are in general desolvated primarily upon binding to the cavity, suggesting that specific hydration of the different chemical units can steer the orientation pathways via a `hydrophobic torque'. Moreover, we find that ligands with bimodal orientation fluctuations have significantly increased kinetic barriers for binding compared to the kinetic barriers previously observed for spherical ligands due to translational fluctuations. We exemplify that these kinetic barriers, which are ligand specific, impact both binding and unbinding times for which we observe considerable differences between our studied ligands.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here we investigate the partitioning of simple monoatomic (Na$^+$, K$^+$, Cs$^+$, Cl$^-$, I$^-$) and one molecular ion (4-nitrophenolate; NP$^-$) within a dense, electroneutral poly($N$-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios $K\simeq\>$10$^{-1}$ than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes we observe that the water clusters possess a universal negative electrostatic potential with respect to their surrounding, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monoatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP$^-$ can have a partition ratio much larger than unity, $K\simeq\>$10-30 (depending on the cation type) or even $10^3$ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP$^-$ reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Competitive sorption of mono- versus di-valent ions by highly charged globular macromolecules

When a highly charged globular macromolecule, such as a dendritic polyelectrolyte or charged nanogel, is immersed into a physiological electrolyte solution, monovalent and divalent counterions from the solution bind to the macromolecule in a certain ratio and thereby almost completely electroneutralize it. For charged macromolecules in biological media, the number ratio of bound mono- versus di-valent ions is decisive for the desired function. A theoretical prediction of such a sorption ratio is challenging because of the competition of electrostatic (valency), ion-specific, and binding saturation effects. Here, we devise and discuss a few approximate models to predict such an equilibrium sorption ratio by extending and combining established electrostatic binding theories such as Donnan, Langmuir, Manning as well as Poisson-Boltzmann approaches, to systematically study the competitive uptake of mono- and di-valent counterions by the macromolecule. We compare and fit our models to coarse-grained (implicit-solvent) computer simulation data of the globular polyelectrolyte dendritic polyglycerol sulfate (dPGS) in salt solutions of mixed valencies. The dPGS has high potential to serve in macromolecular carrier applications in biological systems and at the same time constitutes a good model system for a highly charged macromolecule. We finally use the simulation-informed models to extrapolate and predict electrostatic features such as the effective charge as a function of the divalent ion concentration for a wide range of dPGS generations (sizes).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Controlling the microstructure and phase behavior of confined soft colloids by active interaction switching

We explore the microstructure and phase behavior of confined soft colloids which can actively switch their interactions at a predefined kinetic rate. For this, we employ a reaction-diffusion approach based on a reactive dynamical density-functional theory (R-DDFT) and study the effect of a binary (two-state) switching of the size of colloids interacting with a Gaussian pair potential. The rate of switching interpolates between a near-equilibrium binary Gaussian mixture at low rates and an effective monodisperse liquid for large rates, hence strongly affecting the one-body density profiles, adsorption, and pressure at confining walls. Importantly, we demonstrate that sufficiently fast switching impedes the phase separation of an (in equilibrium) unstable liquid, allowing the dynamic control of the degree of mixing/condensation and local micro-structuring in a cellular confinement by tuning the switching rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enhanced Catalytic Activity of Gold@Polydopamine Nanoreactors with Multi-compartment Structure Under NIR Irradiation

Photothermal conversion (PTC) nanostructures have great potential for applications in many fields, and therefore, they have attracted tremendous attention. However, the construction of a PTC nanoreactor with multi-compartment structure to achieve the combination of unique chemical properties and structural feature is still challenging due to the synthetic difficulties. Herein, we designed and synthesized a catalytically active, PTC gold (Au)@polydopamine (PDA) nanoreactor driven by infrared irradiation using assembled PS-b-P2VP nanosphere as soft template. The particles exhibit multi-compartment structure which is revealed by 3D electron tomography characterization technique. They feature permeable shells with tunable shell thickness. Full kinetics for the reduction reaction of 4-nitrophenol has been investigated using these particles as nanoreactors and compared with other reported systems. Notably, a remarkable acceleration of the catalytic reaction upon near-infrared irradiation is demonstrated, which reveals for the first time the importance of the synergistic effect of photothermal conversion and complex inner structure to the kinetics of the catalytic reduction. The ease of synthesis and fresh insights into catalysis will promote a new platform for novel nanoreactor studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling of stimuli-responsive nanoreactors: rational rate control towards the design of colloidal enzymes

In modern applications of heterogeneous liquid-phase nanocatalysis, the catalysts (e.g., metal nanoparticles) need to be typically affixed to a colloidal carrier system for stability and easy handling. "Passive carriers" (e.g., simple polyelectrolytes) serve for a controlled synthesis of the nanoparticles and prevent coagulation during catalysis. Recently, however, hybrid conjugates of nanoparticles and synthetic thermosensitive polymers have been developed that enable to change the catalytic activity of the nanoparticles by external triggers. In particular, nanoparticles embedded in a stimuli-responsive network made from poly(N-isopropylacrylamide) (PNIPAM) have become the most-studied examples of such hybrids. It has been demonstrated that the permeability of the polymer network and thus the reactant flux can be switched and controlled by external stimuli. Such "active carriers" may thus be viewed as true nanoreactors that open up new design routes in nano-catalysis and elevate synthesis to create highly selective, programmable "colloidal enzymes". However, only a comprehensive understanding of these materials on all time and length scales can lead to a rational design of future, highly functional materials. Here we review the current state of the theoretical and multi-scale simulation approaches, aiming at a fundamental understanding of these nanoreactors. In particular, we summarize a theoretical approach for reaction rates of surface-catalyzed bimolecular reactions in responsive nanoreactors in terms of the key material parameters, the polymer shell permeability P and the reactant partition ratio K. We discuss recent computer simulation studies of both atomistic and coarse-grained polymer models in which these quantities have been characterized in some detail. We conclude with an outlook on selected open questions and future theoretical challenges in nanoreactor modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural and Transport Properties of Li/S Battery Electrolytes: Role of the Polysulfide Species

Lithium--sulfur (Li/S) batteries are regarded as one of the most promising energy storage devices beyond lithium-ion batteries because of their high energy density of 2600 Wh/kg and an affordable cost of sulfur. Meanwhile, some challenges inherent to Li/S batteries remain to be tackled, for instance, the polysulfide (PS) shuttle effect, the irreversible solidification of Li$_2$S, and the volume expansion of the cathode material during discharge. On the molecular level, these issues originate from the structural and solubility behavior of the PS species in bulk and in the electrode confinement. In this study, we use classical molecular dynamics (MD) simulations to develop a working model for PS of different chain lengths in applied electrolyte solutions of lithium bistriflimide (LiTFSI) in 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) mixtures. We investigate conductivities, diffusion coefficients, solvation structures, and clustering behavior and verify our simulation model with experimental measurements available in literature and newly performed by us. Our results show that diffusion coefficients and conductivities are significantly influenced by the chain length of PS. The conductivity contribution of the short chains, like Li$_2$S$_4$, is lower than of longer PS chains, such as Li$_2$S$_6$ or Li$_2$S$_8$, despite the fact that the diffusion coefficient of Li$_2$S$_4$ is higher than for longer PS chains. The low conductivity of Li$_2$S$_4$ can be attributed to its low degree of dissociation and even to a formation of large clusters in the solution. It is also found that an addition of 1 M LiTFSI into PS solutions considerably reduces the clustering behavior. Our simulation model enables future systematic studies in various solvating and confining systems for the rational design of Li/S electrolytes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal compaction of disordered and elastin-like polypeptides: a temperature-dependent, sequence-specific coarse-grained simulation model

Elastin-like polypeptides (ELPs) undergo a sharp solubility transition from low temperature solvated phases to coacervates at elevated temperatures, driven by the increased strength of hydrophobic interactions at higher temperatures. The transition temperature, or 'cloud point', critically depends on sequence composition, sequence length, and concentration of the ELPs. In this work, we present a temperature-dependent, implicit solvent, sequence-specific coarse-grained (CG) simulation model that reproduces the transition temperatures as a function of sequence length and guest residue identity of various experimentally probed ELPs to appreciable accuracy. Our model builds upon the self-organized polymer model introduced recently for intrinsically disordered polypeptides (SOP-IDP), and introduces a semi-empirical functional form for the temperature-dependence of hydrophobic interactions. In addition to the fine performance for various ELPs, we demonstrate the ability of our model to capture the thermal compactions in dominantly hydrophobic intrinsically disordered polypeptides (IDPs), consistent with experimental scattering data. With the high computational efficiency afforded by the CG representation, we envisage that the model will be ideally suited for simulations of large-scale structures such as ELP networks and hydrogels, as well as agglomerates of IDPs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tuning the permeability of dense membranes by shaping nanoscale potentials

The permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed non-monotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anti-correlation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics

Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multi-method approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable ('dry' and 'wet') hydration states. The resulting dry-wet transition rates are then used in a spatially-dependent multi-state continuous-time Markov chain Brownian dynamics simulations, and the related Fokker-Planck equation calculations, of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semi-quantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of non-equilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding versus unbinding. Our study thus provides a significant step forward towards efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Angle-resolved effective potentials for disk-shaped molecules

We present an approach for calculating coarse-grained angle-resolved effective pair potentials for uniaxial molecules. For integrating out the intramolecular degrees of freedom we apply umbrella sampling and steered dynamics techniques in atomistically-resolved molecular dynamics (MD) computer simulations. Throughout this study we focus on disk-like molecules such as coronene. To develop the methods we focus on integrating out the van-der-Waals and intramolecular interactions, while electrostatic charge contributions are neglected. The resulting coarse-grained pair potential reveals a strong temperature and angle dependence. In the next step we fit the numerical data with various Gay-Berne-like potentials to be used in more efficient simulations on larger scales. The quality of the resulting coarse-grained results is evaluated by comparing their pair and many-body structure as well as some thermodynamic quantities self-consistently to the outcome of atomistic MD simulations of many-particle systems. We find that angle-resolved potentials are essential not only to accurately describe crystal structures but also for fluid systems where simple isotropic potentials start to fail already for low to moderate packing fractions. Further, in describing these states it is crucial to take into account the pronounced temperature dependence arising in selected pair configurations due to bending fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coarse-grained electrostatic interactions of coronene: Towards the crystalline phase

In this article we present and compare two different, coarse-grained approaches to model electrostatic interactions of disc-shaped aromatic molecules, specifically coronene. Our study builds on previous work [J. Chem. Phys. 141, 214110 (2014)] where we proposed, based on a systematic coarse-graining procedure starting from the atomistic level, an anisotropic effective (Gay-Berne-like) potential capable of describing van-der-Waals contributions to the interaction energy. To take into account electrostatics, we introduce, first, a linear quadrupole moment along the symmetry axis of the coronene disc. The second approach takes into account the fact that the partial charges within the molecules are distributed in a ring-like fashion. We then reparametrize the effective Gay-Berne-like potential such that it matches, at short distances, the ring-ring potential. To investigate the validity of these two approaches, we perform many-particle Molecular Dynamics (MD) simulations, focusing on the crystalline phase (karpatite) where electrostatic interaction effects are expected to be particularly relevant for the formation of tilted stacked columns. Specifically, we investigate various structural parameters as well as the melting transition. We find that the second approach yields consistent results with those from experiments despite the fact that the underlying potential decays with the wrong distance dependence at large molecule separations. Our strategy can be transferred to a broader class of molecules, such as benzene or hexabenzocoronene.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reaction rate of a composite core-shell nanoreactor with multiple, spatially distributed embedded nano-catalysts

We present a detailed theory for the total reaction rate constant of a composite core-shell nanoreactor, consisting of a central solid core surrounded by a hydrogel layer of variable thickness, where a given number of small catalytic nanoparticles are embedded at prescribed positions and are endowed with a prescribed surface reaction rate constant. Besides the precise geometry of the assembly, our theory accounts explicitly for the diffusion coefficients of the reactants in the hydrogel and in the bulk as well as for their transfer free energy jump upon entering the hydrogel shell. Moreover, we work out an approximate analytical formula for the overall rate constant, which is valid in the physically relevant range of geometrical and chemical parameters. We discuss in depth how the diffusion-controlled part of the rate depends on the essential variables, including the size of the central core. In particular, we derive some simple rules for estimating the number of nanocatalysts per nanoreactor for an efficient catalytic performance in the case of small to intermediate core sizes. Our theoretical treatment promises to provide a very useful and flexible tool for the design of superior performing nanoreactor geometries and with optimized nanoparticle load.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resonance Reaction in Diffusion-Influenced Bimolecular Reactions

We investigate the influence of a stochastically fluctuating step-barrier potential on bimolecular reaction rates by exact analytical theory and stochastic simulations. We demonstrate that the system exhibits a new resonant reaction behavior with rate enhancement if an appropriately defined fluctuation decay length is of the order of the system size. Importantly, we find that in the proximity of resonance the standard reciprocal additivity law for diffusion and surface reaction rates is violated due to the dynamical coupling of multiple kinetic processes. Together, these findings may have important repercussions on the correct interpretation of various kinetic reaction problems in complex systems, as, e.g., in biomolecular association or catalysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sensing solvents with ultra-sensitive porous poly(ionic liquid) actuators

We introduced a new concept for fabricating solvent stimulus polymer actuators with unprecedented sensitivity and accuracy. This was accomplished by integrating porous architectures and electrostatic complexation gradients in a poly(ionic liquid) membrane that bears ionic liquid species for solvent sorption. In contact with 1.5 mol% of acetone molecules in water, the actuator membrane (1 mm x 20 mm x 30 um) bent into a closed loop. While the interaction between solvents and the polymer drives the actuation, the continuous gradient in complexation degree combined with the porous architecture optimizes the actuation, giving it a high sensitivity and even the ability to discriminate butanol solvent isomers. The membrane is also capable of cooperative actuation. The design concept is easy to implement and applicable to other polyelectrolyte systems, which substantially underpins their potentials in smart and sensitive signaling microrobotics/devices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solvent fluctuations induce non-Markovian kinetics in hydrophobic pocket-ligand binding

We investigate the impact of water fluctuations on the key-lock association kinetics of a hydrophobic ligand (key) binding to a hydrophobic pocket (lock) by means of a minimalistic stochastic model system. It describes the collective hydration behavior of the pocket by bimodal fluctuations of a water-pocket interface that dynamically couples to the diffusive motion of the approaching ligand via the hydrophobic interaction. This leads to a set of overdamped Langevin equations in 2D-coordinate-space, that is Markovian in each dimension. Numerical simulations demonstrate locally increased friction of the ligand, decelerated binding kinetics, and local non-Markovian (memory) effects in the ligand's reaction coordinate as found previously in explicit-water molecular dynamics studies of model hydrophobic pocket-ligand binding [1,2]. Our minimalistic model elucidates the origin of effectively enhanced friction in the process that can be traced back to long-time decays in the force-autocorrelation function induced by the effective, spatially fluctuating pocket-ligand interaction. Furthermore, we construct a generalized 1D-Langevin description including a spatially local memory function that enables further interpretation and a semi-analytical quantification of the results of the coupled 2D-system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermosensitive Cu2O-PNIPAM core-shell nanoreactors with tunable photocatalytic activity

We report a facile and novel method for the fabrication of Cu2O@PNIPAM core-shell nanoreactors using Cu2O nanocubes as the core. The PNIPAM shell not only effectively protects the Cu2O nanocubes from oxidation, but also improves the colloidal stability of the system. The Cu2O@PNIPAM core-shell microgels can work efficiently as photocatalyst for the decomposition of methyl orange under visible light. A significant enhancement in the catalytic activity has been observed for the core-shell microgels compared with the pure Cu2O nanocubes. Most importantly, the photocatalytic activity of the Cu2O nanocubes can be further tuned by the thermosensitive PNIPAM shell, as rationalized by our recent theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Curvature Dependence of Hydrophobic Hydration Dynamics

We investigate the curvature-dependence of water dynamics in the vicinity of hydrophobic spherical solutes using molecular dynamics simulations. For both, the lateral and perpendicular diffusivity as well as for H-bond kinetics of water in the first hydration shell, we find a non-monotonic solute-size dependence, exhibiting extrema close to the well-known structural crossover length scale for hydrophobic hydration. Additionally, we find an apparently anomalous diffusion for water moving parallel to the surface of small solutes, which, however, can be explained by topology effects. The intimate connection between solute curvature, water structure and dynamics has implications for our understanding of hydration dynamics at heterogeneous biomolecular surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Equilibrium structure and fluctuations of suspensions of colloidal dumbbells

We investigate the structure and equilibrium linear-response dynamics of suspensions of hard colloidal dumbbells using Brownian Dynamics computer simulations. The focus lies on the dense fluid and plastic crystal states of the colloids with investigated aspect (elongation-to-diameter) ratios varying from the hard sphere limit up to 0.39, which is roughly the stability limit of the plastic crystal phase. We find expected structural changes with larger elongation with respect to the hard sphere reference case and very localized orientational correlations, typically just involving next-neighbor couplings. These relatively weak correlations are also reflected in only minor effects on the translational and rotational diffusion coefficients for most of the investigated elongations. However, the linear response shear viscosity exhibits a dramatic increase at high packing fractions ($ϕ\gtrsim 0.5$) beyond a critical anisotropy factor of about $L^* \simeq 0.15$ which is surprising in view of the relatively weak changes found before on the level of colloidal self-dynamics. We suspect that even for the small investigated anisotropies, newly occurring, collective rotational-translational couplings must be made responsible for the slow time scales appearing in the plastic crystal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interaction of Human Serum Albumin with short Polyelectrolytes: A study by Calorimetry and Computer Simulation

We present a comprehensive study of the interaction of human serum albumin (HSA) with poly(acrylic acid) (PAA; number average degree of polymerization: 25) in aqueous solution. The interaction of HSA with PAA is studied in dilute solution as the function of the concentration of added salt (20 - 100 mM) and temperature (25 - 37$^{\circ}$C). Isothermal titration calorimetry (ITC) is used to analyze the interaction and to determine the binding constant and related thermodynamic data. It is found that only one PAA chain is bound per HSA molecule. The free energy of binding $ΔG_b$ increases with temperature significantly. $ΔG_b$ decreases with increasing salt concentration and is dominated by entropic contributions due to the release of bound counterions. Coarse-grained Langevin computer simulations treating the counterions in an explicit manner are used study the process of binding in detail. These simulations demonstrate that the PAA chains are bound in the Sudlow II site of the HSA. Moreover, $ΔG_b$ is calculated from the simulations and found to be in very good agreement with the measured data. The simulations demonstrate clearly that the driving force of binding is the release of counterions in full agreement with the ITC-data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonequilibrium Structure of Colloidal Dumbbells under Oscillatory Shear

We investigate the nonequilibrium behavior of dense, plastic-crystalline suspensions of mildly anisotropic colloidal hard dumbbells under the action of an oscillatory shear field by employing Brownian dynamics computer simulations. In particular, we extend previous investigations, where we uncovered novel nonequilibrium phase transitions, to other aspect ratios and to a larger nonequilibrium parameter space, that is, a wider range of strains and shear frequencies. We compare and discuss selected results in the context of novel scattering and rheological experiments. Both simulations and experiments demonstrate that the previously found transitions from the plastic crystal phase with increasing shear strain also occur at other aspect ratios. We explore the transition behavior in the strain-frequency phase and summarize it in a nonequilibrium phase diagram. Additionally, the experimental rheology results hint at a slowing down of the colloidal dynamics with higher aspect ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theory of Solvation-Controlled Reactions in Stimuli-Responsive Nanoreactors

Metallic nanoparticles embedded in stimuli-responsive polymers can be regarded as nanoreactors since their catalytic activity can be changed within wide limits: the physicochemical properties of the polymer network can be tuned and switched by external parameters, e.g. temperature or pH, and thus allows a selective control of reactant mobility and concentration close to the reaction site. Based on a combination of Debye's model of diffusion through an energy landscape and a two-state model for the polymer, here we develop an analytical expression for the observed reaction rate constant $k_{\rm obs}$. Our formula shows an exponential dependence of this rate on the solvation free enthalpy change $Δ\bar{G}_{\rm sol}$, a quantity which describes the partitioning of the reactant in the network versus bulk. Thus, changes in $Δ\bar{G}_{\rm sol}$, and not in the diffusion coefficient, will be the decisive factor affecting the reaction rate in most cases. A comparison with recent experimental data on switchable, thermosensitive nanoreactors demonstrates the general validity of the concept.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tuning the critical solution temperature of polymers by copolymerization

We study statistical copolymerization effects on the upper critical solution temperature (CST) of generic homopolymers by means of coarse-grained Langevin dynamics computer simulations and mean-field theory. Our systematic investigation reveals that the CST can change monotonically or non-monotonically with copolymerization, as observed in experimental studies, depending on the degree of non-additivity of the monomer (A-B) cross-interactions. The simulation findings are confirmed and qualitatively explained by a combination of a two-component Flory-de Gennes model for polymer collapse and a simple thermodynamic expansion approach. Our findings provide some rationale behind the effects of copolymerization and may be helpful for tuning CST behavior of polymers in soft material design.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anisotropic Electrostatic Friction of Organic Molecules on ZnO Surfaces

We study the long-time self-diffusion of a single conjugated organic para-sexiphenyl (p-6P) molecule physisorbed on the inorganic ZnO $\left(10\overline{1}0\right)$ surface by means of all-atom molecular dynamics computer simulations. We find strongly anisotropic diffusion processes in which the diffusive motion along the polar [0001] direction of the surface can be many orders of magnitudes slower at relevant experimental temperatures than in the perpendicular direction. The observation can be rationalized by the underlying charge pattern of the electrostatically heterogeneous surface which imposes direction-dependent energy barriers to the motion of the molecule. Furthermore, the diffusive behavior is found to be normal and Arrhenius-like, governed by thermally activated energy barrier crossings. The detailed analysis of the underlying potential energy landscape shows, however, that in general the activation barriers cannot be estimated from idealized zero-temperature trajectories but must include the conformational and positional excursion of the molecule along its pathway. Furthermore, the corresponding (Helmholtz) free energy barriers are significantly smaller than the pure energetic barriers with implications on absolute rate prediction at experimentally relevant temperatures. Our findings suggest that adequately engineered substrate charge patterns could be possibly harvested to select desired growth modes of hybrid interfaces for optoelectronic device engineering.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic density functional theory of protein adsorption on polymer-coated nanoparticles

We present a theoretical model for the description of the adsorption kinetics of globular proteins onto charged core-shell microgel particles based on Dynamic Density Functional Theory (DDFT). This model builds on a previous description of protein adsorption thermodynamics [Yigit \textit{et al}, Langmuir 28 (2012)], shown to well interpret the available calorimetric experimental data of binding isotherms. In practice, a spatially-dependent free-energy functional including the same physical interactions is built, and used to study the kinetics via a generalised diffusion equation. To test this model, we apply it to the case study of Lysozyme adsorption on PNIPAM coated nanoparticles, and show that the dynamics obtained within DDFT is consistent with that extrapolated from experiments. We also perform a systematic study of the effect of various parameters in our model, and investigate the loading dynamics as a function of proteins' valence and hydrophobic adsorption energy, as well as their concentration and that of the nanoparticles. Although we concentrated here on the case of adsorption for a single protein type, the model's generality allows to study multi-component system, providing a reliable instrument for future studies of competitive and cooperative adsorption effects often encountered in protein adsorption experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamic Description of the LCST of Charged Thermoresponsive Copolymers

The dependence of the lower critical solution temperature (LCST) of charged, thermosensitive copolymers on their charge fraction and the salt concentration is investigated by employing systematic cloud-point experiments and analytical theory. The latter is based on the concept of the Donnan equilibrium incorporated into a thermodynamic expansion of a two-state free energy around a charge-neutral reference homopolymer and should be applicable for weakly charged (or highly salted) polymer systems. Very good agreement is found between the theoretical description and the experiments for aqueous solutions of the responsive copolymer poly(NIPAM-co-EVImBr) for a wide range of salt concentrations and charge fractions up to 8%, using only two global, physical fitting parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rationalizing polymer swelling and collapse under repulsive to highly attractive cosolvent conditions

The collapse and swelling behavior of a generic homopolymer is studied using implicit-solvent, explicit-cosolvent Langevin dynamics computer simulations for varying interaction strengths. The systematic investigation reveals that polymer swelling is maximal if both monomer-monomer and monomer-cosolute interactions are weakly attractive. In the most swollen state the cosolute density inside the coil is remarkably bulk-like and homogenous. Highly attractive monomer-cosolute interactions, however, are found to induce a collapse of the chain which, in contrast to the collapsed case induced by purely repulsive cosolvents, exhibits a considerably enhanced cosolute density within the globule. Thus, collapsed states, although appearing similar on a first glance, may result from very different mechanisms with distinct final structural and thermodynamics properties. Two theoretical models, one based on an effective one-component description where the cosolutes have been integrated out, and a fully two-component Flory -- de Gennes like model, support the simulation findings above and serve for interpretation. In particular, the picture is supported that collapse in highly attractive cosolvents is driven by crosslinking-like bridging effects, while the ratio of attraction width to cosolute size plays a critical role behind this mechanism. Only if polymer-cosolute interactions are not too short-ranged swelling effects should be observable. Our findings may be important for the interpretation of the effects of cosolutes on polymer and protein conformational structure, in particular for highly attractive interaction combinations, such as provided by urea, GdmCl, NaI, or NaClO4 near peptide-like moieties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sorption of proteins to charged microgels: characterizing binding isotherms and driving forces

We present a set of Langmuir binding models in which electrostatic cooperativity effects to protein sorption is incorporated in the spirit of Guoy-Chapman-Stern models, where the global substrate (microgel) charge state is modified by bound reactants (charged proteins). Application of this approach to lysozyme sorption to oppositely charged core-shell microgels allows us to extract the intrinsic, binding affinity of the protein to the gel, which is salt-concentration independent and mostly hydrophobic in nature. The total binding affinity is found to be mainly electrostatic in nature, changes many orders of magnitude during the sorption process, and is significantly influenced by osmotic deswelling effects. The intrinsic binding affinity is determined to be about 7 kBT for our system. We additionally show that Langmuir binding models and those based on excludedvolume interactions are formally equivalent for low to moderate protein packing, if the nature of the bound state is consistently defined. Having appreciated this, a more quantitative interpretation of binding isotherms in terms of separate physical interactions is possible in future for a wide variety of experimental approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How interface geometry dictates water's thermodynamic signature in hydrophobic association

As a common view the hydrophobic association between molecular-scale binding partners is supposed to be dominantly driven by entropy. Recent calorimetric experiments and computer simulations heavily challenge this established paradigm by reporting that water's thermodynamic signature in the binding of small hydrophobic ligands to similar-sized apolar pockets is enthalpy-driven. Here we show with purely geometric considerations that this controversy can be resolved if the antagonistic effects of concave and convex bending on water interface thermodynamics are properly taken into account. A key prediction of this continuum view is that for fully complementary binding of the convex ligand to the concave counterpart, water shows a thermodynamic signature very similar to planar (large-scale) hydrophobic association, that is, enthalpy-dominated, and hardly depends on the particular pocket/ligand geometry. A detailed comparison to recent simulation data qualitatively supports the validity of our perspective down to subnanometer scales. Our findings have important implications for the interpretation of thermodynamic signatures found in molecular recognition and association processes. Furthermore, traditional implicit solvent models may benefit from our view with respect to their ability to predict binding free energies and entropies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ion-specificity in α-helical folding kinetics

The influence of the salts KCl, NaCl, and NaI at molar concentrations on the α-helical folding kinetics of the alanine-based oligopeptide Ace-AEAAAKEAAAKA-Nme is investigated by means of (explicit-water) molecular dynamics simulations and a diffusional analysis. The mean first passage times for folding and unfolding are found to be highly salt-specific. In particular, the folding times increase about one order of magnitude for the sodium salts. The drastic slowing down can be traced back to long-lived, compact configurations of the partially folded peptide, in which sodium ions are tightly bound by several carbonyl and carboxylate groups. This multiple trapping is found to lead to a non-exponential residence time distribution of the cations in the first solvation shell of the peptide. The analysis of α-helical folding in the framework of diffusion in a reduced (one-dimensional) free energy landscape further shows that the salt not only specifically modifies equilibrium properties, but also induces kinetic barriers due to individual ion binding. In the sodium salts, for instance, the peptide's configurational mobility (or "diffusivity") can decrease about one order of magnitude. This study demonstrates the highly specific action of ions and highlights the intimate coupling of intramolecular friction and solvent effects in protein folding.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Ionic force field optimization based on single-ion and ion-pair solvation properties

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional $\{σ,ε\}$ parameter space of the Lennard Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair-potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs$^+$ and K$^+$, while for the anions I$^-$ and F$^-$ the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules.

0 citations0 reviewsNo code linkedOpen access