Theory of electrostatic effects in soft biological interfaces using atomic force microscopy.

We calculated the electrostatic force between a planar interface, such as a planar-supported lipid bilayer membrane, and the tip of a stylus on which another lipid bilayer or some other biomacromolecular system might be deposited. We considered styli with rounded tips as well as conical tips. To take into account the effect of dynamical hydrogen-bonded structures in the aqueous phase, we used a theory of nonlocal electrostatics. We used the Derjaguin approximation and identified the systems for which its use is valid. We pointed out where our approach differs from previous calculations and to what extent the latter are inadequate. We found that 1) the nonlocal interactions have significant effects over distances of 10-15 A from the polar zone and that, at the surface of this zone, the effect on the calculated force can be some orders of magnitude; 2) the lipid dipoles and charges are located a distance L from the hydrophobic layer in the aqueous medium and this can have consequences that may not be appreciated if it is ignored; 3) dipoles, located in the aqueous region, can give rise to forces even though the polar layer is unchanged, and if this is ignored the interpretation of force data can be erroneous if an attempt is made to rationalize an observed force with a knowledge of an uncharged surface; 4) the shape of the stylus tip can be very important, and a failure to take this into account can result in incorrect conclusions, a point made by other workers; and 5) when L is nonzero, the presence of charges and dipoles can yield a force that can be nonmonotonic as a function of ionic concentration.     

Reproduction of correct electrostatic field by charges and dipoles on a closed surface

A general algorithm based on the Green function theorem has been developed to correctly reproduce electrostatic fields inside a closed space by point charges and point dipoles on the surface surrounding the space. For actual computations, limited numbers of point charges, including charge pairs replacing point dipoles, are enough to approximate the inner fields. As examples, reaction fields were reproduced by the current surface charges and dipoles for the dielectric models, where a monopole, dipole, or quadrupole was individually set at the center in a vacuum sphere surrounded by high dielectric continuum. The potentials due to those reaction fields agree well with the analytical ones. As an application of this method to the analysis of the electronic structure of the active site of a protein, a combination of the continuum dielectric model and ab initio molecular orbital calculation was carried out. Other applications to molecular dynamics and quantum mechanical calculations are also discussed.     

The intraglobular electrostatic field of an enzyme. II. Effect of environment polarization

The influence of polarization of surrounding medium on the intraglobular electric field of the alpha-chymotrypsin molecule is considered. The polarization is taken into account by the image charges method, the proper approximations for calculation of the fields due to intraglobular and surface charges are suggested. The polarization of surroundings does not change the qualitative picture of the electric field in the active center of the alpha-chymotrypsin molecule set up by protein dipoles, but reduces almost to zero the intraglobular field set up by surface ions.     

Effects of proteins on thermotropic phase transitions of phospholipid membranes.

A variety of proteins have been studied for their ability to interact and alter the thermotropic properties of phospholipid bilayer membranes as detected by differential scanning calorimeter. The proteins studied included: basic myelin protein (A1 protein), cytochrome c, major apoprotein of myelin proteolipid (N-2 apoprotein), gramicidin A, polylysine, ribonuclease and hemoglobin. The lipids used for the interactions were dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylglycerol. The interactions were grouped in three catagories each having very different effects on the phospholipid phase transition from solid to liquid crystalline. The calorimetric studies were also correlated with data from vesicle permeability and monolayer expansion. Ribonuclease and polylysine which exemplify group 1 interactions, show strong dependence on electrostatic binding. Their effects on lipid bilayers include an increase in the enthalpy of transition (deltaH) accompanied by either an increase or no change in the temperature of transition (Tc). In addition, they show minimal effects on vesicle permeability and monolayer expansion. It was concluded that these interactions represent simple surface binding of the protein on the lipid bilayer without penetration into the hydrocarbon region. Cytochrome c and A1 protein, which exemplify group 2 interactions, also show a strong dependence on the presence of net negative charges on the lipid bilayers for their binding. In contrast to the first group, however, they induce a drastic decrease in both Tc and deltaH of the lipid phase transition. Furthermore, they induce a large increase in the permeability of vesicles and a substantial expansion in area of closely packed monolayers at the air-water interface. It was concluded that group 2 interactions represent surface binding followed by partial penetration and/or deformation of the bilayer. Group 3 interactions, shown by proteolipid apoprotein and gramicidin A, were primarily non-polar in character, not requiring electrostatic charges and not inhibited by salt and pH changes. They had no appreciable effect on the Tc but did induce a linear decrease in the magnitude of the deltaH, proportional to the percentage of protein by weight. Membranes containing 50% proteolipid protein still exhibited a thermotropic transition with a deltaH one half that of the pure lipid, and only a small diminution of the size of the cooperative unit. It was concluded that in this case the protein was embedded within the bilayer, associating with a limited number of molecules via non-polar interactions, while the rest of the bilayer was largely unperturbed.     

Voltage dependence of liver gap-junction channels reconstituted into liposomes and incorporated into planar bilayers.

The voltage dependence of rat liver gap junctions was investigated using non-denaturing solubilization and reconstitution of gap-junction protein into proteoliposomes in controlled conditions of connexon aggregation. The presence of liver connexin 32 in reconstituted proteoliposomes was checked with specific antibodies. The proteoliposomes were inserted into planar lipid bilayers by fusion. The single-channel conductance was voltage independent, and its magnitude was 700-1900 pS in 1 M NaCl, as expected from other reports, assuming that conductance is linear with ion activity. The channels were open at zero voltage and completely closed above 40 mV in either direction. This steep voltage dependence corresponded to an open/closed-state voltage difference of 19 mV and to 3.5 gating charges moving through the field. When several channels were inserted into the bilayer, a large fraction of the membrane conductance became voltage insensitive. These results show that the isolated channel units are highly voltage dependent and are consistent with the assumption that aggregated connexons interact through links which prevent voltage-sensitive conformational changes.     

The intrinsic molecular potential of glyceryl monooleate layers and its effect on the conformation and orientation of an inserted molecule: example of gramicidin A.

It is shown by explicit calculation that the distribution of the atomic charges in the constituent molecules of a lipid monolayer or bilayer of glyceryl monooleate creates an intrinsic potential difference between the head region and the hydrocarbon region which tends to repel positive charges towards the exterior and attract negative charges to the interior. The analogies and differences between a bilayer and a monolayer are analyzed. The possible consequences of the intrinsic potential gradient in a lipid layer on the preferred orientation and conformation of a polar neutral molecule are illustrated on the case of a gramicidin A monomer.     

Modulating dipoles for structure-function correlations in the gramicidin A channel.

Dipoles of the tryptophan indole side chains have a direct impact on ion conductance in the gramicidin channel. Here, fluorination of the indoles (both 5- and 6-fluoro) is used to manipulate both the orientations and the magnitudes of the dipoles. The orientations and positions with respect to the channel axis were determined using (2)H solid state NMR of uniformly aligned lipid bilayer preparations. By exchange of the remaining four protons in the indole ring for deuterium, comparison could be made to d(5)-indole spectra that have previously been recorded for each of the four indoles of gramicidin A. After making the assignments which were aided by the observation of (19)F-(2)H dipolar interactions, we found that fluorination caused only minor changes in side chain conformation. With the high-resolution structural characterization of the fluorinated indoles in position 11, 13, and 15, the electrostatic interactions with a cation at the channel and bilayer center can be predicted and the influence of the modified dipoles on ion conductance estimated. The importance of the long-range electrostatic interaction was recently documented with the observation of alpha-helical dipoles oriented toward the bilayer center on the ion conductance pathway for the Streptomyces K(+) channel. We present direct measurements of the orientation of gramicidin channel F-Trp positions for use in analysis of dipole effects on channel permeation.     

Theory of initial current density in membrane voltage-clamp experiments

Herein we show that the voltage-clamp current density at zero time calculated from electrodiffusion equations is linear in the clamping voltage for a simple membrane (no charge structure) and for a membrae with fixed charges. Such membranes are nonexcitable. Excitable membranes can be represented by a homogeneous membrane with dipole layers at the surface. In this case the initial current density will be linear in the clamping voltage if a critical field for a dipole layer reorientation is not passed through in changing from holding to clamping potential. Otherwise, deviation from nonlinearity may occur. This is in agreement with experimental data for the squid giant axon.     

Outer-membrane protein PhoE from Escherichia coli forms anion-selective pores in lipid-bilayer membranes

Porin PhoE of the outer membrane of Escherichia coli was isolated and purified. Reconstitution experiments with lipid bilayer membranes showed that this protein formed pores which had a single channel conductance of 210 pS at 0.1 M KCl. The PhoE pores were obviously not voltage-controlled or regulated. In contrast to pores formed by the OmpF porin from E. coli the PhoE channel was found to be anion-selective at neutral pH. Chloride is about three to ten times more permeable through the pore than alkali ions. On the basis of the observed pH dependence of the permeability ratio of anions and cations, this anionic selectivity is explained by the assumption that the PhoE pore contains an excess of fixed positive charges.     

The intraglobular electrostatic field of an enzyme. 1. The primary field created by the polypeptide core, functional groups and ions of the alpha-chymotrypsin molecule

The electric field set up by the dipoles of peptide groups ad other dipoles at the atoms of substrate and catalytic groups of alpha-chymotrypsin is considered. It is shown that substantial electric potentials reaching some tenths of volts exist in the active center of the enzyme, the fact which must influence significantly the reactivity of corresponding groups. In contrast to low molecular weight liquids, the contribution to the total potential of dipoles located at different distances from the point under consideration often changes nonmonotonically with the distance, sometimes the predominant influence being exerted not by the nearest polar groups but by the more distant ones. The existence of electric fields having a complex spatial configuration determined by the protein structure can be defined as the effect of the polar medium preorganization. Emphasis is placed on the necessity of taking into account the polarization of the external medium by charges of protein atoms and ions (the difference of primary and secondary electric fields).     

An experimental test of new theoretical models for the electrokinetic properties of biological membranes. The effect of UO2++ and tetracaine on the electrophoretic mobility of bilayer membranes and human erythrocytes

For a large smooth particle with charges at the surface, the electrophoretic mobility is proportional to the zeta potential, which is related to the charge density by the Gouy-Chapman theory of the diffuse double layer. This classical model adequately describes the dependence of the electrophoretic mobility of phospholipid vesicles on charge density and salt concentration, but it is not applicable to most biological cells, for which new theoretical models have been developed. We tested these new models experimentally by measuring the effect of UO2++ on the electrophoretic mobility of model membranes and human erythrocytes in 0.15 M NaCl at pH 5. We used UO2++ for these studies because it should adsorb specifically to the bilayer surface of the erythrocyte and should not change the density of fixed charges in the glycocalyx. Our experiments demonstrate that it forms high-affinity complexes with the phosphate groups of several phospholipids in a bilayer but does not bind significantly to sialic acid residues. As observed previously, UO2++ adsorbs strongly to egg phosphatidylcholine (PC) vesicles: 0.1 mM UO2++ changes the zeta potential of PC vesicles from 0 to +40 mV. It also has a large effect on the electrophoretic mobility of vesicles formed from mixtures of PC and the negative phospholipid phosphatidylserine (PS): 0.1 mM UO2++ changes the zeta potential of PC/PS vesicles (10 mol % PS) from -13 to +37 mV. In contrast, UO2++ has only a small effect on the electrophoretic mobility of either vesicles formed from mixtures of PC and the negative ganglioside GM1 or erythrocytes: 0.1 mM UO2++ changes the apparent zeta potential of PC/GM1 vesicles (17 mol % GM1) from -11 to +5 mV and the apparent zeta potential of erythrocytes from -12 to -4 mV. The new theoretical models suggest why UO2++ has a small effect on PC/GM1 vesicles and erythrocytes. First, large groups (e.g., sugar moieties) protruding from the surface of the PC/GM1 vesicles and erythrocytes exert hydrodynamic drag. Second, charges at the surface of a particle (e.g., adsorbed UO2++) exert a smaller effect on the mobility than charges located some distance from the surface (e.g., sialic acid residues).     

Effects of proteins on the thermotropic phase transitions of phospholipid membranes

A variety of proteins have been studied for their ability to interact and alter the thermotropic properties of phospholipid bilayer membranes as detected by differential scanning calorimeter. The proteins studied included: basic myelin protein (A1 protein), cytochrome c, major apoprotein of myelin proteolipid (N-2 apoprotein), gramicidin A, polylysine, ribonuclease and hemoglobin. The lipids used for the interactions were dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylglycerol. The interactions were grouped in three categories each having very different effects on the phospholipid phase transition from solid to liquid crystalline. The calorimetric studies were also correlated with data from vesicle permeability and monolayer expansion.Ribonuclease and polylysine which exemplify group 1 interactions, show strong dependence on electrostatic binding. Their effects on lipid bilayers include an increase in the enthalpy of transition (ΔH) accompanied by either an increase or no change in the temperature of transition (T_c). In addition, they show minimal effects on vesicle permeability and monolayer expansion. It was concluded that these interactions represent simple surface binding of the protein on the lipid bilayer without penetration into the hydrocarbon region.Cytochrome c and Al protein, which exemplify group 2 interactions, also show a strong dependence on the presence of net negative charges on the lipid bilayers for their binding. In contrast to the first group, however, they induce a drastic decrease in both T_c and ΔH of the lipid phase transition. Furthermore, they induce a large increase in the permeability of vesicles and a substantial expansion in area of closely packed monolayers at the air-water interface. It was concluded that group 2 interactions represent surface binding followed by partial penetration and/or deformation of the bilayer.Group 3 interactions, shown by proteolipid apoprotein and gramicidin A, were primarily non-polar in character, not requiring electrostatic charges and not inhibited by salt and pH changes. They had no appreciable effect on the T_c but did induce a linear decrease in the magnitude of the ΔH, proportional to the percentage of protein by weight. Membranes containing 50% proteolipid protein still exhibited a thermotropic transition with a ΔH one half that of the pure lipid, and only a small diminution of the size of the cooperative unit. It was concluded that in this case the protein was embedded within the bilayer, associating with a limited number of molecules via non-polar interactions, while the rest of the bilayer was largely unperturbed.     

Control of reduction potential by protein matrix: lesson from a spherical protein model

 Factors that contribute to the control of reduction potential by protein matrix are examined within a spherical protein model. These include the nonpolar nature of protein matrices, solvent accessibility of the redox center, and net charges and dipoles of surrounding amino acids. Simple rules on their effects are established. In particular, surface charges have little effect on the reduction potential, and polar groups may either increase or decrease the reduction potential, depending on their orientations relative to the redox center. The effects of complex formation, proton titration, and ionic strength are also discussed. Received, accepted: 26 November 1996     

Probing the environment of signal-anchor sequences during topogenesis in the endoplasmic reticulum

Signal sequences for insertion of protein into the mammalian endoplasmic reticulum orient themselves in the translocon on the basis of their flanking charges. It has recently been shown that hydrophobic N-terminal signals initially insert head-on before they invert their orientation to translocate the C-terminus. The rate of inversion is reduced with the increasing hydrophobicity of the signal due to an increased affinity for the initial bound state at the translocon. To probe the environment of the signal while its orientation is determined, different hydrophobic residues were inserted at various positions throughout a uniform oligoleucine signal sequence and the constructs were expressed in transfected COS-7 cells. The resulting topologies revealed a strikingly symmetric position dependence specifically for bulky aromatic amino acids, reflecting the structure of a lipid bilayer. Maximal N-translocation was observed when the guest residues were placed at the N- or C-terminus of the hydrophobic sequence or in the very center, corresponding to the positions of highest expected affinity of the signal sequence as a membrane-spanning helix for the bilayer. The results support the model that during topogenesis in vivo the signal sequence is exposed to the lipid membrane.     

Molecular Dynamics Simulation of Conformational Flexibility of Alamethicin Fragments in Aqueous and Membranous Environment

Abstract

We present here results on molecular dynamics (MD) simulation on two fragments of channel forming antibiotic peptide Alamethicin, containing isoamino butyric acid (Aib). Simulations are carried out in aqueous and membranous environment in a bilayer of 39 molecules of Dimyristoyl phosphatidyl choline (DMPC). The peptides Boc—;Pro-Aib-Ala-Aib- OBzl (Alam 1) and Boc-Leu-Aib-Pro-OBzl (Alam 2) were simulated from their crystallography coordinates. The bilayers were built from two different conformations (A and B) of DMPC reported in crystal data. The P-N dipoles were arranged hexagonally with surface area per lipid molecule 66.5 A°² and P-P separation across the bilayer 34 A°. They were hydrated by 28.6 and 25.5 water molecules per DMPC molecule. Simulations are done using AMBER 4.0 package in constant number volume temperature (NVT) condition for 100 pico seconds (ps) in aqueous environment and 250 ps of equilibrated bilayer. Geometric parameters of lipids as: bilayer thickness, order parameter of the chains, transfraction of chain torsional angles were monitored. We also monitored geometric parameters of the peptides as backbone torsional angles, distances amongst Ca atoms, angles between Cα atoms, movement of center of gravity (CG) along and perpendicular to bilayer normal. We find that membrane bilayer is slightly disturbed due to the presence of peptides. In case of alam 2 in water angles ψ1 and ψ3 showed larger variation in water compared to same in the bilayer. The peptide conformation is more stable in DMPC bilayer. However the peptides showed movement along and perpendicular to bilayer normal. This we believe is due to hydrophobic nature of these peptides.     

Surface dipoles,surface charges and negative steady-state resistance in biological membranes

A theoretical expression for the steady-state current-voltage characteristic of biological membranes is given, justifying various assumptions which are usually made in the framework of the electrodiffusion theory, and taking into account the effect of surface charges and surface dipoles. It is shown that the orientation of the dipoles can be treated by a two-state model. But a comparison with the experimental curves of Gilbert and Ehrenstein leads to the conclusion that the observed negative resistances cannot be explained by a reorientation of the dipoles alone. An additional modification of the surface charges could be sufficient. Some features can be well explained by an absorption of divalent cations on the “outside” surface.     

Molecular Dynamics Simulations with Interaction Potentials Including Polarization Development of a Noniterative Method and Application to Water

A general method is suggested for the implementation of polarization in molecular dynamics simulations of small molecules. Induced dipole moments are evaluated on selected polarizability centers and represented by separation of charges. The positive polarization charges reside on the selected atoms. The negative polarization charges are treated as additional particles. The positions of these polarization charges are determined from the electrical fields due to the permanent charges of the system. Thus the induction is treated explicitly, while the higher order contributions, the polarization due to induced dipoles, are taken into account in an average way by modification of potential parameters. The forces can be evaluated for the new charge distribution in the conventional way. As an illustration of this approach initial results are reported for the development of a polarizable water model. The higher order polarization is treated in an average way by slight increase of the permanent charges as compared to the values that would give the gas phase dipole moment. The increase in CPU time is comparable to the addition of one atom per polarizable center.     

Nonlinear Voltage Dependence of the Light-Driven Proton Pump Current of Bacteriorhodopsin

The light-driven proton pump current generated by bacteriorhodopsin reconstituted in asymmetric planar bilayer membranes was investigated. The current-voltage dependence was found to be nonlinear and can be approximated by an exponential at least below +50 mV. The current changed e-fold when the membrane potential was changed by 80 mV. The voltage dependence was analyzed in terms of a barrier model. This analysis revealed an effective displacement of 0.63 elementary charges across the membrane during the rate-limiting step. Comparison of this value with the results from flash-induced photovoltage signals suggests that one proton is pumped per cycle.