Electrostatics of a simple membrane model using Green's functions formalism

Energy landscapes and dynamics of ion translocation through membrane transporters: a meeting ground for physics, chemistry, and biology

The dynamics of ion translocation through membrane transporters is visualized from a comprehensive point of view by a Gibbs energy landscape approach. The ΔG calculations have been performed with the Kirkwood-Tanford-Warshel (KTW) electrostatic theory that properly takes into account the self-energies of the ions. The Gibbs energy landscapes for translocation of a single charge and an ion pair are calculated, compared, and contrasted as a function of the order parameter, and the characteristics of the frustrated system with bistability for the ion pair are described and quantified in considerable detail. These calculations have been compared with experimental data on the ΔG of ion pairs in proteins. It is shown that, under suitable conditions, the adverse Gibbs energy barrier can be almost completely compensated by the sum of the electrostatic energy of the charge-charge interactions and the solvation energy of the ion pair. The maxima in ΔG_KTW with interionic distance in the bound H⁺ - A^- charge pair on the enzyme is interpreted in thermodynamic and molecular mechanistic terms, and biological implications for molecular mechanisms of ATP synthesis are discussed. The timescale at which the order parameter moves between two stable states has been estimated by solving the dynamical equations of motion, and a wealth of novel insights into energy transduction during ATP synthesis by the membrane-bound F_OF₁-ATP synthase transporter is offered. In summary, a unifying analytical framework that integrates physics, chemistry, and biology has been developed for ion translocation by membrane transporters for the first time by means of a Gibbs energy landscape approach.

Coupling mechanisms in ATP synthesis: Rejoinder to "Response to molecular-level understanding of biological energy coupling and transduction"

Recently, an exchange of views on key fundamental aspects of biological energy coupling and ATP synthesis in the vital process of oxidative phosphorylation appeared in the pages of this journal. The very difficult scientific problems are analyzed and clarified. Errors in the mathematical/thermodynamic equations of a previous analysis have been identified that invalidate previous assertions, and the correct equations are derived. The major differences between the two competing models - localized versus delocalized - for biological energy coupling and transduction are discussed from physical, chemical, and mathematical perspectives. The opposing views are summarized, so that the reader can assess for himself or herself the merits of the two coupling mechanisms. A fresh attempt has been made to go to the root of bioenergetics by calculating the desolvation free energy barrier, ∆G_desolvation for ion transport across biomembranes. Several constructive suggestions are made that have the power to resolve the basic contradictions and the areas of fundamental conflict, and reach a consensus by catalyzing the progress of future research in this interdisciplinary field.

Electric dipole image forces in three-layer systems: The classical electrostatic model

General exact analytical expressions have been derived for the image force energy W_i(Z, φ) of a point dipole in a classical three-layer system composed of dispersionless media with arbitrary constant dielectric permittivities ε_i. Here, i = 1-3 is the layer number, and Z and φ are the dipole coordinate and orientation angle, respectively. It was found that the long-range asymptotics W_i(Z→∞,φ) in both covers (i = 1, 3) are reached unexpectedly far from the interlayer (i = 2). Another specific feature of the solution consists in that the interference of the fields created by polarization charges emerging at both interfaces leads to the appearance of a constant contribution inside the interlayer with a non-standard dependence on the dipole orientation angle φ. It was shown that by changing the dielectric constants of the structure components, one can realize two peculiar regimes of the W_i(Z, φ) behavior in the covers; namely, there arises either a potential barrier preventing adsorption or a well far from the interface, both being of a totally electrostatic origin, i.e., without involving the Pauli exchange repulsion, which is taken into account in the conventional theories of physical adsorption. The results obtained provide a fresh insight into the physics of adsorption in physical electronics, chemical physics, and electrochemistry.

Attraction between neutral dielectrics mediated by multivalent ions in an asymmetric ionic fluid

We study the interaction between two neutral plane-parallel dielectric bodies in the presence of a highly asymmetric ionic fluid, containing multivalent as well as monovalent (salt) ions. Image charge interactions, due to dielectric discontinuities at the boundaries, as well as effects from ion confinement in the slit region between the surfaces are taken fully into account, leading to image-generated depletion attraction, ion correlation attraction, and steric-like repulsive interactions. We investigate these effects by employing a combination of Monte Carlo simulation methods, including explicit-ion simulations (where all electrostatic interactions are simulated explicitly) and implicit-ion simulations (where monovalent ions are replaced by an effective screened electrostatic potential between multivalent ions), as well as an approximate analytical theory. The latter incorporates strong ion-image charge correlations, which develop in the presence of high valency ions in the mixture. We show that the implicit-ion simulations and the proposed analytical theory can describe the explicit simulation results on a qualitative level, while excellent quantitative agreement can be obtained for sufficiently large monovalent salt concentrations. The resultant attractive interaction between the neutral surfaces is shown to be significant, as compared with the usual van der Waals interactions between semi-infinite dielectrics, and can thus play an important role at the nano scale.

Proton-pumping mechanism of cytochrome c oxidase: a kinetic master-equation approach

Cytochrome c oxidase is an efficient energy transducer that reduces oxygen to water and converts the released chemical energy into an electrochemical membrane potential. As a true proton pump, cytochrome c oxidase translocates protons across the membrane against this potential. Based on a wealth of experiments and calculations, an increasingly detailed picture of the reaction intermediates in the redox cycle has emerged. However, the fundamental mechanism of proton pumping coupled to redox chemistry remains largely unresolved. Here we examine and extend a kinetic master-equation approach to gain insight into redox-coupled proton pumping in cytochrome c oxidase. Basic principles of the cytochrome c oxidase proton pump emerge from an analysis of the simplest kinetic models that retain essential elements of the experimentally determined structure, energetics, and kinetics, and that satisfy fundamental physical principles. The master-equation models allow us to address the question of how pumping can be achieved in a system in which all reaction steps are reversible. Whereas proton pumping does not require the direct modulation of microscopic reaction barriers, such kinetic gating greatly increases the pumping efficiency. Further efficiency gains can be achieved by partially decoupling the proton uptake pathway from the active-site region. Such a mechanism is consistent with the proposed Glu valve, in which the side chain of a key glutamic acid shuttles between the D channel and the active-site region. We also show that the models predict only small proton leaks even in the absence of turnover. The design principles identified here for cytochrome c oxidase provide a blueprint for novel biology-inspired fuel cells, and the master-equation formulation should prove useful also for other molecular machines. .

Influence of a nanorod molecular layer on the biological activity of neuronal cells. A semiclassical model for complex solid/liquid interfaces with carbon nanotubes

A general account of electric effects is given for a biological phase interacting with a nanorod molecular layer by means of the formed hard-soft and solid-liquid interfaces. In particular, the frequency enhancement previously detected for the spontaneous activity of neuronal cultures interfaced with carbon nanotubes is quantitatively explained upon a quantum/semiclassical description, where the duration of a biological signal is viewed as the (average) lifetime of a decaying state (or population of states), and the effect of the carbon phase as a linewidth broadening. Four contributions were principally accounted for, one biological, for the synaptic strength, one electrochemical, for the overall capacitance increase implied by the nanotube double layers, one geometric, for the typical scales ruling the electron and ion conduction mechanisms, and one electromagnetic-like, translating the membrane polarization changes. These calculations predict an enhancement factor equal on average to ≃6.39, against a former experimental value ≃6.08.

Arginine in membranes: the connection between molecular dynamics simulations and translocon-mediated insertion experiments

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

Effects on capacitance by overexpression of membrane proteins

Functional Channelrhodopsin-2 (ChR2) overexpression of about 10(4)channels/mum(2) in the plasma membrane of HEK293 cells was studied by patch-clamp and freeze-fracture electron microscopy. Simultaneous electrorotation measurements revealed that ChR2 expression was accompanied by a marked increase of the area-specific membrane capacitance (C(m)). The C(m) increase apparently resulted partly from an enlargement of the size and/or number of microvilli. This is suggested by a relatively large C(m) of 1.15+/-0.08 microF/cm(2) in ChR2-expressing cells measured under isotonic conditions. This value was much higher than that of the control HEK293 cells (0.79+/-0.02 microF/cm(2)). However, even after complete loss of microvilli under strong hypoosmolar conditions (100 mOsm), the ChR2-expressing cells still exhibited a significantly larger C(m) (0.85+/-0.07 microF/cm(2)) as compared to non-expressing control cells (0.70+/-0.03 microF/cm(2)). Therefore, a second mechanism of capacitance increase may involve changes in the membrane permittivity and/or thickness due to the embedded ChR2 proteins.

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.

Dipole lattice membrane model for protein calculations

A dipole lattice model for lipid membranes and their interactions with peptides is presented. It uses the Langevin dipole method to calculate electrostatic interactions in the heterogeneous membrane environment. A series of test cases are presented, including spherical charges, dipoles, side chain analogs, and helical peptides. The model consistently produces qualitatively correct results.

The interaction between lipid derivatives of colchicine and tubulin: consequences of the interaction of the alkaloid with lipid membranes

Colchicine is a potent antimitotic poison which is well known to prevent microtubule assembly by binding tubulin very tightly. Colchicine also possesses anti-inflammatory properties which are not well understood yet. Here we show that colchicine tightly interacts with lipid layers. The physical and biological properties of three different lipid derivatives of colchicine are investigated parallel to those of membrane lipids in the presence of colchicine. Upon insertion in the fatty alkyl chains, colchicine rigidifies the lipid monolayers in a fluid phase and fluidifies rigid monolayers. Similarly X-ray diffraction data show that lecithin-water phases are destabilized by colchicine. In addition, an unexpectedly drastic enhancement of the photoisomerization rate of colchicine into lumicolchicine in the lipid environment is observed and further supports insertion of the alkaloid in membranes. Finally the interaction of colchicine with lipids makes the drug inaccessible to tubulin. The possible in vivo significance of these results is discussed.

Direct measurement of specific membrane capacitance in neurons

The specific membrane capacitance (C(m)) of a neuron influences synaptic efficacy and determines the speed with which electrical signals propagate along dendrites and unmyelinated axons. The value of this important parameter remains controversial. In this study, C(m) was estimated for the somatic membrane of cortical pyramidal neurons, spinal cord neurons, and hippocampal neurons. A nucleated patch was pulled and a voltage-clamp step was applied. The exponential decay of the capacitative charging current was analyzed to give the total membrane capacitance, which was then divided by the observed surface area of the patch. C(m) was 0.9 microF/cm(2) for each class of neuron. To test the possibility that membrane proteins may alter C(m), embryonic kidney cells (HEK-293) were studied before and after transfection with a plasmid coding for glycine receptor/channels. The value of C(m) was indistinguishable in untransfected cells and in transfected cells expressing a high level of glycine channels, indicating that differences in transmembrane protein content do not significantly affect C(m). Thus, to a first approximation, C(m) may be treated as a "biological constant" across many classes of neuron.

(In)validity of the constant field and constant currents assumptions in theories of ion transport

Constant electric fields and constant ion currents are often considered in theories of ion transport. Therefore, it is important to understand the validity of these helpful concepts. The constant field assumption requires that the charge density of permeant ions and flexible polar groups is virtually voltage independent. We present analytic relations that indicate the conditions under which the constant field approximation applies. Barrier models are frequently fitted to experimental current-voltage curves to describe ion transport. These models are based on three fundamental characteristics: a constant electric field, negligible concerted motions of ions inside the channel (an ion can enter only an empty site), and concentration-independent energy profiles. An analysis of those fundamental assumptions of barrier models shows that those approximations require large barriers because the electrostatic interaction is strong and has a long range. In the constant currents assumption, the current of each permeating ion species is considered to be constant throughout the channel; thus ion pairing is explicitly ignored. In inhomogeneous steady-state systems, the association rate constant determines the strength of ion pairing. Among permeable ions, however, the ion association rate constants are not small, according to modern diffusion-limited reaction rate theories. A mathematical formulation of a constant currents condition indicates that ion pairing very likely has an effect but does not dominate ion transport.

Electrostatics and the ion selectivity of ligand-gated channels

The nicotinic acetylcholine receptor (nAChR) is a cation-selective ion channel that opens in response to acetylcholine binding. The related glycine receptor (GlyR) is anion selective. The pore-lining domain of each protein may be modeled as a bundle of five parallel M2 helices. Models of the pore-lining domains of homopentameric nAChR and GlyR have been used in continuum electrostatics calculations to probe the origins of ion selectivity. Calculated pKA values suggest that "rings" of acidic or basic side chains at the mouths of the nAChR or GlyR M2 helix bundles, respectively, may not be fully ionized. In particular, for the nAChR the ring of glutamate side chains at the extracellular mouth of the pore is predicted to be largely protonated at neutral pH, whereas those glutamate side chains in the intracellular and intermediate rings (at the opposite mouth of the pore) are predicted to be fully ionized. Inclusion of the other domains of each protein represented as an irregular cylindrical tube in which the M2 bundles are embedded suggests that both the M2 helices and the extramembrane domains play significant roles in determining ion selectivity.