Driving Forces of Protein Diffusion

A single-particle energy-conserving dissipative particle dynamics approach for simulating thermophoresis of nanoparticles in polymer networks

The transport of nanoparticles in polymer networks has critical implications in biology and medicine, especially through thermophoresis in response to temperature gradients. This study presents a single-particle energy-conserving dissipative particle dynamics (seDPD) method by integrating a single-particle model into the energy-conserving DPD model to simulate the mesoscopic thermophoretic behavior of nanoparticles in polymer matrices. We first validate the newly developed seDPD model through comparisons with analytical solutions for nanoparticle viscosity, thermal diffusivity, and hydrodynamic drag and then demonstrate the effectiveness of the seDPD model in capturing thermophoretic forces induced by temperature gradients. The results show that nanoparticles driven by the Soret forces exhibit unique transport characteristics, such as drift velocity and diffusivity, leading to a significant acceleration of nanoparticle diffusion in the polymer network, which has been known as the giant acceleration of diffusion. Quantifying how nanoparticles move in flexible polymer networks sheds light on the interaction dynamics of nanoparticles within polymer networks, providing insight into nanoparticle behavior in complex environments that could be leveraged in various applications from drug delivery to material design.

Remarkable Insensitivity of Protein Diffusion to Protein Charge

Friction to translational diffusion of ionic particles in polar liquids should scale linearly with the squared ion charge, according to standard theories. Substantial slowing of translational diffusion is expected for proteins in water. In contrast, our simulations of charge mutants of green fluorescent proteins in water show remarkable insensitivity of the translational diffusion constant to protein's charge in the range of charges between -29 and +35. The friction coefficient is given as a product of the force variance and the memory function relaxation time. We find remarkably accurate equality between the variance of the electrostatic force and the negative cross-correlation of the electrostatic and van der Waals forces. The charge invariance of the diffusion constant is a combined effect of the force variance and relaxation time invariances with the protein charge. The temperature dependence of the protein diffusion constant is highly non-Arrhenius, with a fragile-to-strong crossover at the glass transition.

Do Electrostatics Control the Diffusive Dynamics of Solitary Water? NMR and MD Studies of Water Translation and Rotation in Dipolar and Ionic Solvents

NMR-based measurements of the diffusion coefficients and rotation times of solitary water and benzene at 300 K are reported in a diverse collection of 13 conventional organic solvents and 10 imidazolium ionic liquids. Proton chemical shifts of water are found to be correlated to water OH-stretching frequencies, confirming the importance of electrostatic interactions in these shifts. However, the influence of magnetic interactions in aromatic solvents renders chemical shifts a less reliable indicator of electrostatics. Diffusion coefficients (DB) and rotational correlation times (τB) of benzene in the solvents examined are accurately described as functions of viscosity (η) by DB ∝ η-0.81 and τB ∝ η0.64. Literature values of DB and τB in alkane and normal alcohols, which were not included among the solvents studied here, are systematically faster than predicted by these correlations, indicating that factors beyond solvent viscosity play a role in determining the friction on benzene. In contrast to benzene, water diffusion and rotation are poorly described in terms of viscosity alone, even in the dipolar and ionic solvents measured here. The present data and the substantial literature data already available on dilute water diffusion show a systematic dependence of DW on solvent polarity among isoviscous solvents. The aspect of solvent polarity most relevant to water dynamics is the ability of a solvent to accept hydrogen bonds from water, as conveniently quantified by the frequency of water's OH stretching band, ΔνOH. The friction on translation, ζtr = kBT/DW, and rotation, ζrot = kBTτW, are both well correlated by functions of the form ζ(η, ΔνOH) = a1ηa2 exp (a3ΔνOH), where the ai are adjustable parameters. Molecular dynamics simulations reveal a strong coupling between electrostatic and nonelectrostatic water-solvent interactions, which makes it impossible to dissect the friction on water into additive dielectric and hydrodynamic components. Simulations also provide a tentative explanation for the unusual form of the correlating function ζ(η, ΔνOH), at least in the case of ζrot.

Protein translational diffusion as a way to detect intermolecular interactions

In this work, we analyze the information on the protein intermolecular interactions obtained from macromolecular diffusion. We have shown that the most hopeful results are given by our approach based on analysis of protein translational self-diffusion and collective diffusion obtained by dynamic light scattering and pulsed-field gradient NMR (PFG NMR) spectroscopy with the help of Vink's approach to analyze diffusion motion of particles by frictional formalism of non-equilibrium thermodynamics and the usage of the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloid particles interactions in electrolyte solutions. Early we have shown that integration of Vink's theory with DLVO provides a reliable basis for uniform interpreting of PFG NMR and DLS experiments on concentration dependence of diffusion coefficients. Basic details of theoretical and mathematical procedures and a broad analysis of experimental attestation of proposed conception on proteins of various structural form, size, and shape are presented. In the present review, the main capabilities of our approach obtain the details of intermolecular interactions of proteins with different shapes, internal structures, and mass. The universality of Vink's approach is experimentally shown, which gives the appropriate description of experimental results for proteins of complicated structure and shape.

Relaxation dynamics measure the aggregation propensity of amyloid-β and its mutants

Atomistic molecular dynamics simulations are employed to investigate the global and segmental relaxation dynamics of the amyloid-β protein and its causative and protective mutants. Amyloid-β exhibits significant global/local dynamics that span a broad range of length and time scales due to its intrinsically disordered nature. The relaxation dynamics of the amyloid-β protein and its mutants is quantitatively correlated with its experimentally measured aggregation propensity. The protective mutant has slower relaxation dynamics, whereas the causative mutants exhibit faster global dynamics compared with that of the wild-type amyloid-β. The local dynamics of the amyloid-β protein or its mutants is governed by a complex interplay of the charge, hydrophobicity, and change in the molecular mass of the mutated residue.

Effect of Protein–Protein Interactions on Translational Diffusion of Spheroidal Proteins

One of the commonly accepted approaches to estimate protein–protein interactions (PPI) in aqueous solutions is the analysis of their translational diffusion. The present review article observes a phenomenological approach to analyze PPI effects via concentration dependencies of self- and collective translational diffusion coefficient for several spheroidal proteins derived from the pulsed field gradient NMR (PFG NMR) and dynamic light scattering (DLS), respectively. These proteins are rigid globular α-chymotrypsin (ChTr) and human serum albumin (HSA), and partly disordered α-casein (α-CN) and β-lactoglobulin (β-Lg). The PPI analysis enabled us to reveal the dominance of intermolecular repulsion at low ionic strength of solution (0.003–0.01 M) for all studied proteins. The increase in the ionic strength to 0.1–1.0 M leads to the screening of protein charges, resulting in the decrease of the protein electrostatic potential. The increase of the van der Waals potential for ChTr and α-CN characterizes their propensity towards unstable weak attractive interactions. The decrease of van der Waals interactions for β-Lg is probably associated with the formation of stable oligomers by this protein. The PPI, estimated with the help of interaction potential and idealized spherical molecular geometry, are in good agreement with experimental data.

War and Peace between Electrostatic and van der Walls Forces Regulate Translational and Rotational Diffusion

Diffusion of a Brownian particle or a molecule in a liquid solvent is caused, in the Stokes-Einstein picture, by unbalanced fluctuations of osmotic forces on different sides of the particle. When the particle carries a charge or a higher multipolar moment, this picture is amended by fluctuations of electrostatic forces producing dielectric friction. Dielectric friction slows both the translational and rotational diffusion down. While this picture is well established and is physically sound, standard theories grossly overestimate the magnitude of dielectric friction for small dipolar solutes and larger colloidal particles, such as proteins. Motivated by recent simulation studies, this perspective discusses the interplay between osmotic (van der Waals) and electrostatic forces in promoting molecular and colloidal diffusion. Much can be learned about microscopic friction mechanisms from statistical and dynamical correlations between osmotic and electrostatic forces.

Ionic mobility driven by correlated van der Waals and electrostatic forces

Classical theories of dielectric friction make two critical assumptions: (i) friction due to van der Waals (vdW) forces is described by hydrodynamic drag and is independent of the ionic charge and (ii) vdW and electrostatic forces are statistically independent. Both assumptions turn out to be incorrect when tested against simulations of anions and cations with varying charge magnitude dissolved in water. Both the vdW and electrostatic components of the force variance scale linearly with the ionic charge squared. The two components are strongly anticorrelated producing simple relations for the total force variance in terms of self-variances. The inverse diffusion constant scales linearly with the charge squared. Solvation asymmetry between cations and anions extends to linear transport coefficients.

Ewald sum corrections in simulations of ion and dipole solvation and electron transfer

Periodic boundary conditions and Ewald sums used in standard simulation protocols require finite-size corrections when the total charge of the simulated system is nonzero. Corrections for ion solvation were introduced by Hummer, Pratt, and García, [J. Chem. Phys. 107, 9275 (1997)]. The latter approach is extended here to derive finite-size correction for the Stokes-shift and reorganization energy applied to electron-transfer reactions. The same correction term, scaling inversely with the box size, adds to the reorganization energy from the energy-gap variance but is subtracted from the reorganization energy calculated from the Stokes shift. Finite-size corrections thus widen the gap between these two quantities, which were recently found to diverge for protein electron transfer. Corrections to the free energy of dipole solvation and the variance of the electric field scale as m²/L³ with the solute dipole m and the box size L.

Enhanced Molecular Diffusivity through Destructive Interference between Electrostatic and Osmotic Forces

Molecular charge asymmetrically distributed in a diffusing tagged particle causes a nonzero electrostatic force balanced by an opposing van der Waals (vdW) force. Fluctuations of electrostatic and vdW forces are highly correlated, and they destructively interfere in the force variance. This phenomenology is caused by the formation of a structurally frozen hydration layer for a particle diffusing in water and is responsible for a substantial speedup of translational diffusion compared to traditional theories of dielectric friction. Diffusion of proteins is insensitive to charge mutations, while smaller particles with asymmetric charge distribution can show a strong dependence of translational and rotational diffusion on molecular charge. Dielectric calculations of the electrostatic force require low values of ≃5 for the effective dielectric constant of interfacial water to be consistent with simulations.

Modulation of Nanoparticle Diffusion by Surface Ligand Length and Charge: Analysis with Molecular Dynamics Simulations

To help better interpret experimental measurement of nanoparticle size, it is important to understand how their diffusion depends on the physical and chemical features of surface ligands. In this study, explicit solvent molecular dynamics simulations are used to probe the effect of ligand charge and flexibility on the diffusion of small gold nanoparticles. The results suggest that despite a high bare charge (+18 e), cationic nanoparticles studied here have reduced diffusion constants compared to a hydrophobic gold nanoparticle by merely a modest amount. Increasing the ligand length by 10 CH₂ units also has a limited impact on the diffusion constant. For the three particles studied here, the difference between estimated hydrodynamic radius and radius of gyration is on the order of one solvent layer (3-5 Å), confirming that the significant discrepancies found in the size of similar nanoparticles by recent transmission electron microscopy and dynamic light scattering measurements were due to aggregation under solution conditions. The limited impact of electrostatic friction on the diffusion of highly charged nanoparticles is found to be due to the strong anticorrelation between electrostatic and van der Waals forces between nanoparticle and environment, supporting the generality of recent observation for proteins by Matyushov and co-workers. Including the first shell of solvent molecules as part of the diffusing particle has a minor impact on the total force autocorrelation function but reduces the disparity in relaxation time between the total force and its electrostatic and van der Waals components.