Non-Brownian diffusion in lipid membranes: Experiments and simulations.

The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Nonergodicity, fluctuations, and criticality in heterogeneous diffusion processes.

We study the stochastic behavior of heterogeneous diffusion processes with the power-law dependence D(x) ∼ |x|(α) of the generalized diffusion coefficient encompassing sub- and superdiffusive anomalous diffusion. Based on statistical measures such as the amplitude scatter of the time-averaged mean-squared displacement of individual realizations, the ergodicity breaking and non-Gaussianity parameters, as well as the probability density function P(x,t), we analyze the weakly nonergodic character of the heterogeneous diffusion process and, particularly, the degree of irreproducibility of individual realizations. As we show, the fluctuations between individual realizations increase with growing modulus |α| of the scaling exponent. The fluctuations appear to diverge when the critical value α = 2 is approached, while for even larger α the fluctuations decrease, again. At criticality, the power-law behavior of the mean-squared displacement changes to an exponentially fast growth, and the fluctuations of the time-averaged mean-squared displacement do not converge for increasing number of realizations. From a systematic comparison we observe some striking similarities of the heterogeneous diffusion process with the familiar subdiffusive continuous time random walk process with power-law waiting time distribution and diverging characteristic waiting time.

Structure-driven homology pairing of chromatin fibers: the role of electrostatics and protein-induced bridging.

Chromatin domains formed in vivo are characterized by different types of 3D organization of interconnected nucleosomes and architectural proteins. Here, we quantitatively test a hypothesis that the similarities in the structure of chromatin fibers (which we call "structural homology") can affect their mutual electrostatic and protein-mediated bridging interactions. For example, highly repetitive DNA sequences in heterochromatic regions can position nucleosomes so that preferred inter-nucleosomal distances are preserved on the surfaces of neighboring fibers. On the contrary, the segments of chromatin fiber formed on unrelated DNA sequences have different geometrical parameters and lack structural complementarity pivotal for stable association and cohesion. Furthermore, specific functional elements such as insulator regions, transcription start and termination sites, and replication origins are characterized by strong nucleosome ordering that might induce structure-driven iterations of chromatin fibers. We propose that shape-specific protein-bridging interactions facilitate long-range pairing of chromatin fragments, while for closely-juxtaposed fibers electrostatic forces can in addition yield fine-tuned structure-specific recognition and pairing. These pairing effects can account for some features observed for mitotic and inter-phase chromatins.

Electrical monitoring of polyelectrolyte multilayer formation by means of capacitive field-effect devices.

The semiconductor field-effect platform represents a powerful tool for detecting the adsorption and binding of charged macromolecules with direct electrical readout. In this work, a capacitive electrolyte-insulator-semiconductor (EIS) field-effect sensor consisting of an Al-p-Si-SiO2 structure has been applied for real-time in situ electrical monitoring of the layer-by-layer formation of polyelectrolyte (PE) multilayers (PEM). The PEMs were deposited directly onto the SiO2 surface without any precursor layer or drying procedures. Anionic poly(sodium 4-styrene sulfonate) and cationic weak polyelectrolyte poly(allylamine hydrochloride) have been chosen as a model system. The effect of the ionic strength of the solution, polyelectrolyte concentration, number and polarity of the PE layers on the characteristics of the PEM-modified EIS sensors have been studied by means of capacitance-voltage and constant-capacitance methods. In addition, the thickness, surface morphology, roughness and wettabilityof the PE mono- and multilayers have been characterised by ellipsometry, atomic force microscopy and water contact-angle methods, respectively. To explain potential oscillations on the gate surface and signal behaviour of the capacitive field-effect EIS sensor modified with a PEM, a simplified electrostatic model that takes into account the reduced electrostatic screening of PE charges by mobile ions within the PEM has been proposed and discussed.

Detection of DNA hybridization by field-effect DNA-based biosensors: mechanisms of signal generation and open questions.

We model theoretically the electrostatic effects taking place upon DNA hybridization in dense DNA arrays immobilized on a layer of Au nano-particles deposited on the surface of a field-effect-based DNA capacitive biosensor. We consider the influence of separation of a charged analyte from the sensor surface and the salinity of electrolyte solution, in the framework of both linear and nonlinear Poisson-Boltzmann theories. The latter predicts a substantially weaker sensor signals due to electrostatic saturation effects that is the main conclusion of this paper. We analyze how different physical parameters of dense DNA brushes affect the magnitude of hybridization signals. The list includes the fraction of DNA charge neutralization, the length and spatial conformations of adsorbed DNA molecules, as well as the discreteness of DNA charges. We also examine the effect of Donnan ionic equilibrium in DNA lattices on the sensor response. The validity of theoretical models is contrasted against recent experimental observations on detection of DNA hybridization via its intrinsic electric charge. The sensitivity of such biochemical sensing devices, their detection limit, and DNA hybridization efficiency are briefly discussed in the end.

Polyelectrolyte adsorption onto oppositely charged interfaces: image-charge repulsion and surface curvature.

We analyze theoretically the influence of low-dielectric boundaries on the adsorption of flexible polyelectrolytes onto planar and spherical oppositely charged surfaces in electrolyte solutions. We rationalize to what extent polymer chains are depleted from adsorbing interfaces by repulsive image forces. We employ the WKB (Wentzel-Kramers-Brillouin) quantum mechanical method for the Green function of the Edwards equation to determine the adsorption equilibrium. Scaling relations are determined for the critical adsorption strength required to initiate polymer adsorption onto these low-dielectric supports. Image-force repulsion shifts the equilibrium toward the desorbed state, demanding larger surface charge densities and polyelectrolyte linear charge densities for the adsorption to take place. The effect is particularly pronounced for a planar interface in a low-salt regime, where a dramatic change in the scaling behavior for the adsorption-desorption transition is predicted. For the adsorbed state, polymers with higher charge densities are displaced further from the interface by image-charge repulsions. We discuss relevant experimental evidence and argue about possible biological applications of the results.

Critical polyelectrolyte adsorption under confinement: planar slit, cylindrical pore, and spherical cavity.

We explore the properties of adsorption of flexible polyelectrolyte chains in confined spaces between the oppositely charged surfaces in three basic geometries. A method of approximate uniformly valid solutions for the Green function equation for the eigenfunctions of polymer density distributions is developed to rationalize the critical adsorption conditions. The same approach was implemented in our recent study for the "inverse" problem of polyelectrolyte adsorption onto a planar surface, and on the outer surface of rod-like and spherical obstacles. For the three adsorption geometries investigated, the theory yields simple scaling relations for the minimal surface charge density that triggers the chain adsorption, as a function of the Debye screening length and surface curvature. The encapsulation of polyelectrolytes is governed by interplay of the electrostatic attraction energy toward the adsorbing surface and entropic repulsion of the chain squeezed into a thin slit or small cavities. Under the conditions of surface-mediated confinement, substantially larger polymer linear charge densities are required to adsorb a polyelectrolyte inside a charged spherical cavity, relative to a cylindrical pore and to a planar slit (at the same interfacial surface charge density). Possible biological implications are discussed briefly in the end.

Polyelectrolyte adsorption onto oppositely charged interfaces: unified approach for plane, cylinder, and sphere.

A universal description is presented for weak adsorption of flexible polyelectrolyte chains onto oppositely charged planar and curved surfaces. It is based on the WKB (Wentzel-Kramers-Brillouin) quantum mechanical method for the Green function equation in the ground state dominance limit. The approach provides a unified picture for the scaling behavior of the critical characteristics of polyelectrolyte adsorption and the thickness of the adsorbed polymer layer formed adjacent to the interface. We find, particularly at low-salt conditions, that curved convex surfaces necessitate much larger surface charge densities to trigger polyelectrolyte adsorption, as compared to a planar interface in the same solution. In addition, we demonstrate that the different surface geometries yield very distinct scaling laws for the critical surface charge density required to initiate chain adsorption. Namely, in the low-salt limit, the surface charge density scales cubical with the inverse Debye screening length for a plane, quadratic for an adsorbing cylinder, and linear for a sphere. As the radius of surface curvature grows, the parameter of critical chain adsorption onto a rod and a sphere turns asymptotically into that of a planar interface. The transition occurs when the radius of surface curvature becomes comparable to the Debye screening length. The general scaling trends derived appear to be consistent with the complex-formation experiments of polyelectrolyte chains with oppositely charged spherical and cylindrical micelles. Finally, the WKB results are compared with the existing theories of polyelectrolyte adsorption and future perspectives are outlined.

Electrostatic interactions in biological DNA-related systems.

In this perspective article, we focus on recent developments in the theory of charge effects in biological DNA-related systems. The electrostatic effects on different levels of DNA organization are considered, including the DNA-DNA interactions, DNA complexation with cationic lipid membranes, DNA condensates and DNA-dense cholesteric phases, protein-DNA recognition, DNA wrapping in nucleosomes, and inter-nucleosomal interactions. For these systems, we develop a theoretical framework to describe the physical-chemical mechanisms of structure formation and anticipate some biological consequences. General biophysical principles of DNA compaction in chromatin fibers and DNA spooling inside viral capsids are discussed in the end, with emphasis on electrostatic aspects.

DNA cyclization: suppression or enhancement by electrostatic repulsions?

First, we develop a model of counterion condensation on highly charged polyelectrolyte rings. Using the known analytical results for the electrostatic energy of ring formation, a stronger counterion adsorption is anticipated onto a cyclized polyelectrolyte, as compared to the Manning prediction for a straight rod-like polyelectrolyte. This fact ensures a lower energetic cost of polyelectrolyte bending into a ring. In the main part of the work, we investigate the impact of charges on cyclization of short DNA fragments, both theoretically and by computer simulations. An approximate expression for the electrostatically renormalized DNA cyclization probability is proposed that incorporates the electrostatic energies of polyelectrolyte cyclization and dimerization reactions. Depending on concentration of simple salt and chain length, the probability of formation of ideal polyelectrolyte rings can be either electrostatically inhibited or enhanced. The latter effect is quite counterintuitive. Afterward, simple computer simulations are performed to enumerate the effects of DNA thermal fluctuations onto the electrostatic energies of cyclized and dimerized DNA fragments in solution. Their outcomes support the possibility of electrostatically enhanced polyelectrolyte ring formation reaction in solution. In the end, we discuss some implications of the results obtained for the future DNA cyclization experiments and provide a short analysis of possible DNA-related features neglected in the modeling.

DNA-DNA sequence homology recognition: physical mechanisms and open questions.

First, we summarize recent experimental facts on homologous DNA pairing in vitro and discuss possible mechanisms of DNA-DNA sequence recognition. Then, we overview the mechanism of DNA-DNA recognition based on complementarity of DNA charge patterns. The theory predicts the recognition energy up to 10 k(B) T for close parallel homologous DNA fragments of gene-relevant lengths. We argue why this estimate cannot be directly applied to pairing of homologous DNA loci in experiments on yeast chromosomes. Namely, DNA-DNA distances assessed from experiments are much larger than those typically used in theory. Finally, we suggest some experiments that could help to judge whether short-range electrostatic forces indeed govern DNA pairing. This viewpoint paper introduces recently developed theoretical concepts to molecular biologists, with a hope to generate a junction between theory and future experiments on DNA recognition.

Looping charged elastic rods: applications to protein-induced DNA loop formation.

We analyze looping of thin charged elastic filaments under applied torques and end forces, using the solution of linear elasticity theory equations. In application to DNA, we account for its polyelectrolyte character and charge renormalization, calculating electrostatic energies stored in the loops. We argue that the standard theory of electrostatic persistence is only valid when the loop's radius of curvature and close-contact distance are much larger than the Debye screening length. We predict that larger twist rates are required to trigger looping of charged rods as compared with neutral ones. We then analyze loop shapes formed on charged filaments of finite length, mimicking DNA looping by proteins with two DNA-binding domains. We find optimal loop shapes at different salt amounts, minimizing the sum of DNA elastic, DNA electrostatic, and protein elastic energies. We implement a simple model where intercharge repulsions do not affect the loop shape directly but can choose the energy-optimized shape from the allowed loop types. At low salt concentrations more open loops are favored due to enhanced repulsion of DNA charges, consistent with the results of computer simulations on formation of DNA loops by lac repressor. Then, we model the precise geometry of DNA binding by the lac tetramer and explore loop shapes, varying the confined DNA length and protein opening angle. The characteristics of complexes obtained, such as the total loop energy, stretching forces required to maintain its shape, and the reduction of electrostatic energy with increment of salt, are in good agreement with the outcomes of more elaborate numerical calculations for lac-repressor-induced DNA looping.

Collapse of highly charged polyelectrolytes triggered by attractive dipole-dipole and correlation-induced electrostatic interactions.

In the first part of the paper, we study the collapse of flexible highly charged polyelectrolyte chains induced by attractive dipole-dipole interactions. The latter emerge due to the formation of dipoles between the chain monomers and counterions condensed on the polyelectrolyte from solution. Using the statistics of slightly perturbed Gaussian polymers, we obtain the scaling relations for the chain dimensions as a function of polyelectrolyte linear charge density in the limit of compacting chains. The results are in good agreement with the outcomes of recent molecular dynamics simulations of the collapse of flexible polyelectrolytes in the presence of explicit counterions. In the second part, we analyze the results of molecular dynamics simulations for the complex formation by two highly charged polyelectrolyte chains carrying opposite charges. We use the scaling arguments based on the picture of complexation of electrostatic blobs in order to rationalize the size of the complexes of two polyelectrolyte chains in the collapsed state. Similar scaling relation for the complex size was recently obtained in computer simulations of complexation of diblock polyampholytes and was described theoretically on the basis of similar electrostatic blob concepts. We also analyze the density of complexes formed and polyelectrolyte linear charge densities required for the onset onto the collapse as a function of interaction strength between the monomers. In both parts of the paper, we overview the scaling relationships obtained for similar systems with alternating charges from other theoretical approaches.

Probing DNA-DNA electrostatic friction in tight superhelical DNA plies.

We estimate theoretically the strength of DNA-DNA electrostatic friction forces emerging upon a slow drag of one DNA over another one in a close juxtaposition. For ideally helical DNA duplexes, this friction occurs due to correlations in electrostatic potential near the DNA surface. The latter originate from the intrinsic helicity of DNA phosphates and adsorbed cations on a scale of 3.4 nm. They produce positive-negative charge interlocking along the DNA-DNA contact. For realistic nonideally helical DNAs, where electrostatic potential barriers become decorrelated due to accumulation of mismatches in DNA structure, DNA-DNA frictional forces are strongly impeded. We discuss possibilities of probing the DNA-DNA intermolecular interactions in strongly confined DNA superhelical plies, as obtained in single-molecule experiments.

Positively charged residues in DNA-binding domains of structural proteins follow sequence-specific positions of DNA phosphate groups.

We study electrostatic charge complementarity along interfaces of DNA-protein complexes. We use the Protein Data Bank atomic coordinates of DNA-protein complexes for some DNA-binding proteins to study the distribution of positively charged protein residues in the close contact with DNA. We show that large structural proteins reveal a peculiar nonuniform distribution of Arg, Lys, and His amino acids in the frame of negatively charged DNA phosphate strands. We study the nucleosome core particles, DNA complexes with prokaryotic DNA-bending histone analogues, but also the basic binding motifs of small DNA-binding proteins. For large DNA-protein complexes, where extensive DNA wrapping around protein cores occurs, we show that positive amino acids on the proteins track sequence-specific positions of individual DNA phosphates. This specificity of electrostatic interactions can contribute to DNA recognition by DNA-binding proteins, which is governed for many DNA-protein complexes primarily by the hydrogen bond formation between protein residues and DNA bases.

DNA cholesteric phases: the role of DNA molecular chirality and DNA-DNA electrostatic interactions.

DNA molecules form dense liquid-crystalline twisted phases both in vivo and in vitro. How the microscopic DNA chirality is transferred into intermolecular twist in these mesophases and what is the role of chiral DNA-DNA electrostatic interactions is still not completely clear. In this paper, we first give an extended overview of experimental observations on DNA cholesteric phases and discuss the factors affecting their stability. Then, we consider the effects of steric and electrostatic interactions of grooved helical molecules on the sign of cholesteric twist. We present some theoretical results on the strength of DNA-DNA chiral electrostatic interactions, on DNA-DNA azimuthal correlations in cholesteric phases, on the value of DNA cholesteric pitch, and on the regions of existence of DNA chiral phases stabilized by electrostatic interactions. We suggest for instance that 146 bp long DNA fragments with stronger affinities for the nucleosome formation can form less chiral cholesteric phases, with a larger left-handed cholesteric pitch. Also, the value of left-handed pitch formed in assemblies of homologous DNA fragments is predicted to be smaller than that of randomly sequenced DNAs. We expect also the cholesteric assemblies of several-kbp-long DNAs to require higher external osmotic pressures for their stability than twisted phases of short nucleosomal DNA fragments at the same DNA lattice density.

Protein--DNA interactions: reaching and recognizing the targets.

Protein searching and recognizing the targets on DNA was the subject of many experimental and theoretical studies. It is often argued that some proteins are capable of finding their targets 10-100 times faster than predicted by the three-dimensional diffusion rate. However, recent single-molecule experiments showed that the diffusion constants of the protein motion along DNA are usually small. This controversy pushed us to revisit this problem. We present a theoretical approach that describes some physical-chemical aspects of the target search and recognition. We consider the search process as a sequence of cycles, with each cycle consisting of three-dimensional and one-dimensional tracks. It is argued that the search time contains three terms: for the motion on three-dimensional and one-dimensional segments, and the correlation term. Our analysis shows that the acceleration in the search time is achieved at some intermediate strength of the protein-DNA binding energy and it is partially "apparent" because it is in fact reached by parallel scanning for the target by many proteins. We also show how the complementarity of the charge patterns on a target DNA sequence and on the protein may result in electrostatic recognition of a specific track on DNA and subsequent protein pinning. Within the scope of a model, we obtain an analytical expression for the capturing well. We estimate the depth and width of such a potential well and the typical time that a protein spends in it.

Effect of a low-dielectric interior on DNA electrostatic response to twisting and bending.

We study the effects of a low-dielectric core of rod-like macromolecules on their electrostatic persistence lengths. We use the exact solution of the linear Poisson-Boltzmann equation for the potential of a charge on the surface of a low-dielectric cylinder. We apply the results to the B-DNA molecule, modeled as a double helical array of discrete charges wound on the surface of a low-dielectric rod. For this charge geometry, we calculate the change in the electrostatic twist persistence as compared to DNA with a water-permeable core. We also discuss possible effects of the low-dielectric molecular core on DNA bending persistence.

Electrostatics of DNA complexes with cationic lipid membranes.

We present the exact solutions of the linear Poisson-Boltzmann equation for several problems relevant to electrostatics of DNA complexes with cationic lipids. We calculate the electrostatic potential and electrostatic energy for lamellar and inverted hexagonal phases, concentrating on the effects of dielectric boundaries. We compare our results for the complex energy with the known results of numerical solution of the nonlinear Poisson-Boltzmann equation. Using the solution for the lamellar phase, we calculate the compressibility modulus and compare our findings with the experimental data available. Also, we treat charge-charge interactions across, along, and between two low-dielectric membranes. We obtain an estimate for the strength of electrostatic interactions of one-dimensional DNA smectic layers across the lipid membrane. We discuss in the end some aspects of two-dimensional DNA condensation and DNA-DNA attraction in the DNA-lipid lamellar phase in the presence of di- and trivalent cations. We analyze the equilibrium DNA-DNA separations in lamellar complexes using the recently developed theory of electrostatic interactions of DNA helical charge motifs.

Strong and weak adsorptions of polyelectrolyte chains onto oppositely charged spheres.

We investigate the complexation of long thin polyelectrolyte (PE) chains with oppositely charged spheres. In the limit of strong adsorption, when strongly charged PE chains adapt a definite wrapped conformation on the sphere surface, we analytically solve the linear Poisson-Boltzmann equation and calculate the electrostatic potential and the energy of the complex. We discuss some biological applications of the obtained results. For weak adsorption, when a flexible weakly charged PE chain is localized next to the sphere in solution, we solve the Edwards equation for PE conformations in the Hulthen potential, which is used as an approximation for the screened Debye-Huckel potential of the sphere. We predict the critical conditions for PE adsorption. We find that the critical sphere charge density exhibits a distinctively different dependence on the Debye screening length than for PE adsorption onto a flat surface. We compare our findings with experimental measurements on complexation of various PEs with oppositely charged colloidal particles. We also present some numerical results of the coupled Poisson-Boltzmann and self-consistent field equation for PE adsorption in an assembly of oppositely charged spheres.