Membrane stability under electrical stress: A nonlocal electroelastic treatment

Electroelastic coupling between membrane surface fluctuations and membrane-embedded charges: continuum multidielectric treatment

The coupling of electric fields and charges with membrane-water interfacial fluctuations affects membrane electroporation, ionic conductance, and voltage gating. A modified continuum model is introduced to study charge interaction with membrane-water interfacial fluctuations in multidielectric environments. By surrounding a point charge with a low dielectric sphere, the linear Poisson-Boltzmann equation is directly solved by calculating the reaction field potential via a method that eliminates singularity contributions. This allows treatment of charges located at dielectric boundaries. Two complementary mechanisms governing charge-fluctuation interactions are considered: (1) electroelastic deformation (EED), treating the membrane as an elastic slab (smectic bilayer model), and (2) electrohydrophobic solvation (EHS), accounting for water penetration into the membrane's hydrophobic core. EED often leads to large membrane thickness perturbations, far larger than those consistent with elastic model descriptions [M. B. Partenskii, G. V. Miloshevsky, and P. C. Jordan, Isr. J. Chem. 47, 385 (2007)]. We argue that a switch from EED to EHS can be energetically advantageous at intermediate perturbation amplitudes. Both perturbation mechanisms are simulated by introducing adjustable shapes optimized by the kinetic Monte Carlo reaction path following approach [G. V. Miloshevsky and P. C. Jordan, J. Chem. Phys. 122, 214901 (2005)]. The resulting energy profiles agree with those of recent atomistic molecular dynamics studies on translating a charged residue across a lipid bilayer [S. Dorairaj and T. W. Allen, Proc. Natl. Acad. Sci. U.S.A. 104, 4943 (2007)].

"Squishy capacitor" model for electrical double layers and the stability of charged interfaces

Negative capacitance (NC), predicted by various electrical double layer (EDL) theories, is critically reviewed. Physically possible for individual components of the EDL, the compact or diffuse layer, it is strictly prohibited for the whole EDL or for an electrochemical cell with two electrodes. However, NC is allowed for the artificial conditions of sigma control, where an EDL is described by the equilibrium electric response of electrolyte to a field of fixed, and typically uniform, surface charge-density distributions, sigma. The contradiction is only apparent; in fact local sigma cannot be set independently, but is established by the equilibrium response to physically controllable variables, i.e., applied voltage phi (phi control) or total surface charge q (q control). NC predictions in studies based on sigma control signify potential instabilities and phase transitions for physically realizable conditions. Building on our previous study of phi control [M. B. Partenskii and P. C. Jordan, Phys. Rev. E 77, 061117 (2008)], here we analyze critical behavior under q control, clarifying the basic picture using an exactly solvable "squishy capacitor" toy model. We find that phi can change discontinuously in the presence of a lateral transition, specify stability conditions for an electrochemical cell, analyze the origin of the EDL's critical point in terms of compact and diffuse serial contributions, and discuss perspectives and challenges for theoretical studies not limited by sigma control.

Photosynthetic electron transport and proton flux under moderate heat stress

Moderate heat stress has been reported to increase PSI cyclic electron flow (CEF). We subjected leaves of Arabidopsis (Arabidopsis thaliana) mutants disrupted in the regulation of one or the other pathway of CEF flow-crr2 (chlororespiratory reduction, deficient in regulation of chloroplast NAD(P)H dehydrogenase-dependent CEF) and pgr5 (proton gradient regulation, proposed to have reduced efficiency of antimycin-A-sensitive-CEF regulation) to moderate heat stress. Light-adapted leaves were switched from 23 to 40 degrees C in 2 min. Gas exchange, chlorophyll fluorescence, the electrochromic shift (ECS), and P700 were measured. Photosynthesis of crr2 and pgr5 was more sensitive to heat and had less ability to recover than the genetic background gl. The proton conductance in light was increased by heat and it was twice as much in pgr5, which had much smaller light-induced proton motive force. We confirmed that P700 becomes more reduced at high temperature and show that, in contrast, the proportion of PSII open centers (with Q (A) oxidized) increases. The two mutants had much slower P700(+) reduction rate during and after heat than gl. The proportion of light absorbed by PSI versus PSII was increased in gl and crr2 during and after heat treatment, but not in pgr5. We propose that heat alters the redox balance away from PSII and toward PSI and that the regulation of CEF helps photosynthesis tolerate heat stress.

Transbilayer Pores Induced by Thickness Fluctuations

Thermally-induced fluctuations of individual phospholipids in a bilayer lipid membrane (BLM) are converted into collective motions due to the intermolecular interactions. Here, we demonstrate that transbilayer stochastic pores can be generated via collective thermal movements (CTM). Using the elastic theory of continuous media applied to smectic-A liquid crystals, we estimate the pore radius and the energetic requirements for pore appearance. Three types of thermally-induced transbilayer pores could be formed through BLMs: open and stable, open and unstable, and closed. In most of the situations, two open and stable pores with different radii could be generated. Notably, the two pores have the same generation probability. Unstable pores are possible to appear across thin bilayers that contain phospholipids with a large polar headgroup. Closed pores are present throughout the cases that we have inspected. The effects of hydrophobic thickness, polar headgroup size of phospholipids, temperature, surface tension, and elastic compression on the pore formation and pore stability have been examined as well.

Ionic permeation free energy in gramicidin: a semimicroscopic perspective

Ion permeation through the gramicidin channel is studied using a model that circumvents two major difficulties inherent to standard simulational methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization, accounting for the influence of electronic polarization by embedding the channel in a dielectric milieu that describes this polarization in a mean sense. The explicit mobile moieties are the ion, multipolar waters, and the carbonyls and amides of the peptide backbone. The model treats the influence of aromatic residues and the membrane dipole potential. A new electrical geometry is introduced that treats long-range electrostatics exactly and avoids problems related to periodic boundary conditions. It permits the translocating ion to make a seamless transition from nearby electrolyte to the channel interior. Other degrees of freedom (more distant bulk electrolyte and nonpolar lipid) are treated as dielectric continua. Reasonable permeation free energy profiles are obtained for potassium, rubidium, and cesium; binding wells are shallow and the central barrier is small. Estimated cationic single-channel conductances are smaller than experiment, but only by factors between 2 (rubidium) and 50 (potassium). When applied to chloride the internal barrier is large, with a corresponding miniscule single-channel conductance. The estimated relative single-channel conductances of gramicidin A, B, and C agree well with experiment.

Dynamic domain formation in membranes: thickness-modulation-induced phase separation

A simple model investigates the amplification of fluctuations on membranes constituted of two lipids having different lengths. Van der Waals and electrostatic interactions across the lipid bilayer result in a destabilization favoring thickness variations of the membrane. Close to spontaneous demixing of the two components, the additional gain in free energy due to thickness undulations shifts the stability boundary which promotes phase separation into domains. Interestingly, this effect can be induced by an applied electric field or membrane potential. In biological systems, the dynamic model presented here indicates that electric fields might be important for controlling phase separation and the formation of domains called "rafts".

Simulation of pore formation in lipid bilayers by mechanical stress and electric fields

We describe computer simulations of pore formation and membrane rupture of phospholipid bilayers under mechanical and electrical stress. On the nanosecond simulation time scale, pores are induced by a lateral pressure exceeding -200 bar or by an applied electric field of 0.5 V/nm.

Modeling permeation energetics in the KcsA potassium channel

The thermodynamics of cation permeation through the KcsA K(+) channel selectivity filter is studied from the perspective of a physically transparent semimicroscopic model using Monte Carlo free energy integration. The computational approach chosen permits dissection of the separate contributions to ionic stabilization arising from different parts of the channel (selectivity filter carbonyls, single-file water, cavity water, reaction field of bulk water, inner helices, ionizable residues). All features play important roles; their relative significance varies with the ion's position in the filter. The cavity appears to act as an electrostatic buffer, shielding filter ions from structural changes in the inner pore. The model exhibits K(+) vs. Na(+) selectivity, and roughly isoenergetic profiles for K(+) and Rb(+), and discriminates against Cs(+), all in agreement with experimental data. It also indicates that Ba(2+) and Na(+) compete effectively with permeant ions at a site near the boundary between the filter and the cavity, in the vicinity of the barium blocker site.

Voltage-induced nonconductive pre-pores and metastable single pores in unmodified planar lipid bilayer

Electric fields promote pore formation in both biological and model membranes. We clamped unmodified planar bilayers at 150-550 mV to monitor transient single pores for a long period of time. We observed fast transitions between different conductance levels reflecting opening and closing of metastable lipid pores. Although mean lifetime of the pores was 3 +/- 0.8 ms (250 mV), some pores remained open for up to approximately 1 s. The mean amplitude of conductance fluctuations (approximately 500 pS) was independent of voltage and close for bilayers of different area (40,000 and 10 microm(2)), indicating the local nature of the conductive defects. The distribution of pore conductance was rather broad (dispersion of approximately 250 pS). Based on the conductance value and its dependence of the ion size, the radius of the average pore was estimated as approximately 1 nm. Short bursts of conductance spikes (opening and closing of pores) were often separated by periods of background conductance. Within the same burst the conductance between spikes was indistinguishable from the background. The mean time interval between spikes in the burst was much smaller than that between adjacent bursts. These data indicate that opening and closing of lipidic pores proceed through some electrically invisible (silent) pre-pores. Similar pre-pore defects and metastable conductive pores might be involved in remodeling of cell membranes in different biologically relevant processes.