Tension fluctuations in contracting myofibrils and their interpretation

Repulsive pressure has been measured as a function of lattice spacing in gels of tobacco mosaic virus (TMV) and in the filament lattice of vertebrate striated muscle. External pressures up to ten atm have been applied to these lattices by an osmotic stress method. Numerical solutions to the Poisson-Boltzmann equation in hexagonal lattices have been obtained and compared to the TMV and muscle data. The theoretical curves using values for k calculated from the ionic strength give a good fit to experimental data from TMV gels, and an approximate fit to that from the muscle lattice, provided that a charge radius for the muscle thick filaments of approximately 16 nm is assumed. Variations in ionic strength, sarcomere length and state of the muscle give results which agree qualitatively with the theory, though a good fit between experiment and theory in the muscle case will clearly require consideration of other types of forces. We conclude that Poisson-Boltzmann theory can provide a good first approximation to the long-range electrostatic forces operating in such biological gel systems.     

The electrostatic potential of B-DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.     

Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

Membrane surface-charge titration probed by gramicidin A channel conductance.

We manipulate lipid bilayer surface charge and gauge its influence on gramicidin A channel conductance by two strategies: titration of the lipid charge through bulk solution pH and dilution of a charged lipid by neutral. Using diphytanoyl phosphatidylserine (PS) bilayers with CsCl aqueous solutions, we show that the effects of lipid charge titration on channel conductance are masked 1) by conductance saturation with Cs+ ions in the neutral pH range and 2) by increased proton concentration when the bathing solution pH is less than 3. A smeared charge model permits us to separate different contributions to the channel conductance and to introduce a new method for "bilayer pKa" determination. We use the Gouy-Chapman expression for the charged surface potential to obtain equilibria of protons and cations with lipid charges. To calculate cation concentration at the channel mouth, we compare different models for the ion distribution, exact and linearized forms of the planar Poisson-Boltzmann equation, as well as the construction of a "Gibbs dividing surface" between salt bath and charged membrane. All approximations yield the intrinsic pKain of PS lipid in 0.1 M CsCl to be in the range 2.5-3.0. By diluting PS surface charge at a fixed pH with admixed neutral diphytanoyl phosphatidylcholine (PC), we obtain a conductance decrease in magnitude greater than expected from the electrostatic model. This observation is in accord with the different conductance saturation values for PS and PC lipids reported earlier (, Biochim. Biophys. Acta. 552:369-378) and verified in the present work for solvent-free membranes. In addition to electrostatic effects of surface charge, gramicidin A channel conductance is also influenced by lipid-dependent structural factors.     

Modeling the electrophoresis of lysozyme. II. Inclusion of ion relaxation.

In this work, boundary element methods are used to model the electrophoretic mobility of lysozyme over the pH range 2-6. The model treats the protein as a rigid body of arbitrary shape and charge distribution derived from the crystal structure. Extending earlier studies, the present work treats the equilibrium electrostatic potential at the level of the full Poisson-Boltzmann (PB) equation and accounts for ion relaxation. This is achieved by solving simultaneously the Poisson, ion transport, and Navier-Stokes equations by an iterative boundary element procedure. Treating the equilibrium electrostatics at the level of the full rather than the linear PB equation, but leaving relaxation out, does improve agreement between experimental and simulated mobilities, including ion relaxation improves it even more. The effects of nonlinear electrostatics and ion relaxation are greatest at low pH, where the net charge on lysozyme is greatest. In the absence of relaxation, a linear dependence of mobility and average polyion surface potential, (lambda zero)s, is observed, and the mobility is well described by the equation [formula: see text] where epsilon 0 is the dielectric constant of the solvent, and eta is the solvent viscosity. This breaks down, however, when ion relaxation is included and the mobility is less than predicted by the above equation. Whether or not ion relaxation is included, the mobility is found to be fairly insensitive to the charge distribution within the lysozyme model or the internal dielectric constant.     

Effects of dielectric saturation on planar electric double layers in salt solutions

The effects of dielectric saturation on planar electric double layers in salt solutions are examined by solving the Poisson-Boltzmann equation analytically where the dielectric constant is given as a function of the electric displacement. The activity and the distribution of small ions, the surface potential and the Donnan potential are calculated. The salt exclusion parameter and the Donnan potential decrease while the surface potential increases as a result of the dielectric saturation. The electrostatic entropy is affected considerably by the dielectric saturation while the electrostatic energy is little influenced. Generally, the effects of dielectric saturation on the distribution of small ions and the thermodynamic properties are enhanced by the addition of salt.     

Electrostatic contribution to the bending of DNA

A model is derived that accounts for the short-range electrostatic contribution to the bending of DNA molecule in solution and in complexes with proteins in terms of the non-linear Poisson-Boltzmann equation. We defined that the short-range electrostatic interactions depend on the changes of the polyion surface charge density under deformation, while the long-range interactions depend on the bending-induced changes in distances between each two points along the polyion axis. After an appropriate simplification of the Poisson-Boltzmann equation, the short-range term is calculated separately giving the lower limit for the electrostatic contribution to the DNA persistence length. The result is compared with the theoretical approaches developed earlier [M. Fixman, J. Chem. Phys. 76 (1982) 6346; M. Le Bret, J. Chem. Phys. 76 (1982) 6243] and with the experimental data. The conclusion is made that the results of Fixman-Le Bret, which took into account both types of the electrostatic interactions for a uniformly bent polyion, give the upper limit for the electrostatic persistence length at low ionic strength, and the actual behavior of the DNA persistence length lies between two theoretical limits. Only the short-range term is significant at moderate-to-high ionic strength where our results coincide with the predictions of Fixman-Le Bret. The bending of DNA on the protein surface that is accompanied by an asymmetric neutralization of the DNA charge is also analyzed. In this case, the electrostatic bending energy gives a significant favorite contribution to the total bending energy of DNA. Important implications to the mechanisms of DNA-protein interactions, particularly in the nucleosome particle, are discussed.     

Specificity of desmin to avian and mammalian muscle cells.

The extent of invariance and heterogeneity in desmin, the major component of the muscle form of 100 Å filaments, has been investigated in avian and mammalian muscle and nonmuscle cells with two-dimensional gel electrophoresis and indirect immunofluorescence. Desmin from chick, duck and quail, smooth, skeletal and cardiac muscle cells is resolved into two isoelectric variants, α and β, with each possessing the same charge and electrophoretic mobility in all three avian species irrespective of muscle type. Guinea pig and rat muscle desmin resolves into only one variant; it also possesses the same charge and electrophoretic mobility in the two mammalian species, but it is more acidic and slower in electrophoretic mobility than the two avian variants.In immunofluorescence, desmin is localized together with α-actinin along myofibril Z lines. Antibodies to chick smooth muscle desmin, prepared against the protein purified by preparative SDS gel electrophoresis prior to immunization, cross-react with myofibril Z lines in all three avian species. These antibodies do not cross-react with either rat or guinea pig myofibril Z lines. Similarly, they do not cross-react with avian or mammalian nonmuscle cells grown in tissue culture and known to contain cytoplasmic 100 Å filaments.These results demonstrate that desmin is highly conserved within avian muscle cells and within mammalian muscle cells. It is, however, both biochemically and immunologically distinguishable between avian and mammalian muscle cells, and between muscle and nonmuscle cells. We conclude that there are biochemically and immunologically specific forms of desmin for avian and mammalian muscle cells. Furthermore, within a particular vertebrate species, there are at least two separate classes of 100 Å filaments: the muscle class whose major component is desmin, and the nonmuscle class whose major component is distinct from desmin. Taking into consideration the immunological specificity reported by other laboratories for the 100 Å filaments in glial cells, for neurofilaments and for the epidermal 80 Å keratin filaments, we propose that a given vertebrate species contains at least four major distinguishable classes of 100 Å filaments: muscle 100 Å filaments (desmin filaments), glial filaments, neurofilaments and epidermal keratin filaments.     

Electrostatic potential between surfaces bearing ionizable groups in ionic equilibrium with physiologic saline solution

In the past calculations of electrostatic properties of interacting cellular surfaces have been restricted by assumptions of fixed surface charge or surface potential. For the most part these calculations have been confined to a linear approximation and neglect the small but important complement of divalent cations in the cellular environment. In the present paper these limitations are removed. Solutions are obtained to the full non-linear Poisson-Boltzmann equation, treating the fraction of dissociated ionizable surface groups as a self-consistent functional of the electrostatic potential.     

Structure-based prediction of DNA-binding sites on proteins using the empirical preference of electrostatic potential and the shape of molecular surfaces

Protein-DNA interactions play an essential role in the genetic activities of life. Many structures of protein-DNA complexes are already known, but the common rules on how and where proteins bind to DNA have not emerged. Many attempts have been made to predict protein-DNA interactions using structural information, but the success rate is still about 80%. We analyzed 63 protein-DNA complexes by focusing our attention on the shape of the molecular surface of the protein and DNA, along with the electrostatic potential on the surface, and constructed a new statistical evaluation function to make predictions of DNA interaction sites on protein molecular surfaces. The shape of the molecular surface was described by a combination of local and global average curvature, which are intended to describe the small convex and concave and the large-scale concave curvatures of the protein surface preferentially appearing at DNA-binding sites. Using these structural features, along with the electrostatic potential obtained by solving the Poisson-Boltzmann equation numerically, we have developed prediction schemes with 86% and 96% accuracy for DNA-binding and non-DNA-binding proteins, respectively.     

The δ-function counterion condensate around a line charge derived from the Fuoss,Katchalsky, and Lifson polyelectrolyte theory

It has recently been proven that the counterion condensate around an isolated line charge in an electrolyte, as characterized by nonlinear Poisson-Boltzmann theory, is an encapsulating δ-function. Here the identical result is shown to hold in the framework of the polyelectrolyte theory of Fuoss, Katchalsky, and Lifson. The proof fully exploits analytic solutions to the differential equation which are not available for the nonlinear, cylindrical Poisson-Boltzmann equation.     

Glycocalyx electrostatic potential profile analysis: ion, pH, steric, and charge effects

The Poisson-Boltzmann equation is modified to consider charge ionogenicity, steric exclusion, and charge distribution in order to describe the perimembranous electrostatic potential profile in a manner consistent with the known morphology and biochemical composition of the cell's glycocalyx. Exact numerical and approximate analytical solutions are given for various charge distributions and for an extended form of the Donnan potential model. The interrelated effects of ionic conditions, bulk pH, ion binding, local dielectric, steric volume exclusion, and charge distribution on the local potential, pH, and charge density within the glycocalyx are examined. Local charge-induced, potential-mediated pH reductions cause glycocalyx charge neutralization. Under certain conditions, local potentials may be insensitive to ionic strength or may decrease in spite of increasing charge density. The volume exclusion of the glycocalyx reduces the local ion concentration, thereby increasing the local potential. With neutral lipid membranes, the Donnan and surface potential agree if the glycocalyx charge distribution is both uniform and several times thicker than the Debye length (approximately 20 A in thickness under physiological conditions). Model limitations in terms of application to microdomains or protein endo- and ectodomains are discussed.     

Equilibrium surface charge distribution in phospholipid vesicles. I. Method of calculation

This paper presents a method of calculation of the surface charge equilibrium distribution between the two surfaces of a spherically closed phospholipid bilayer suspended in aqueous electrolyte solution. The net surface charge is supposed to be provided by the ionized polar groups of the phospholipid molecules. Its equilibrium distribution is found by minimization of the free electrostatic energy. The procedure of minimization utilizes the solution of the Poisson-Boltzmann equation which describes the double electric layers of the membrane and an expression for the membrane potential derived under the assumption of absence of charges in the membrane phase. An analytical solution of the problem in the range of validity of the linearized Poisson-Boltzman equation is obtained. It is shown that in this case an equilibrium transmembrane potential exists, and the surface charge density is greater at the outer surface of the vesicle.     

Electrostatic repulsion of ion penetrable charged membranes: role of Donnan potential

A model is proposed for the electrostatic repulsion between two ion-penetrable charged membranes in which fixed charges are uniformly distributed. This model assumes that the electric potential far inside the membrane is always equal to the Donnan potential, independent of the membrane separation. In this respect the present model completely differs from usual models for the electrostatic interaction of colloidal particles which assume that the surface potential or the surface charge density remains constant during interaction. It is shown that the rise in potential in the interacting membranes caused by their approach is greatly suppressed so that the potential in the membrane does not exceed the Donnan potential. Numerical results of the calculation of the repulsion by the non-linear Poisson-Boltzmann equation are displayed as a function of the membrane separation and an approximate formula is also derived.     

The relationship between the poisson-boltzmann model and the condensation hypothesis: an analysis based on the low salt form of the Donnan coefficient

Two common models for the interaction of counterions with cylindrical polyions are considered in the context of the Donnan membrane equilibrium. General analytic expressions are obtained from the Poisson-Boltzmann equation for the Donnan coefficient in terms of the potential at the surface of the polyion or the local concentration of unbound ions at the surface. Analysis based on these expressions shows that if, and only if, the polyion charge density exceeds a certain critical value a large local concentration of ions will persist near the polyion surface at low ionic strengths. We therefore conclude that this principal hypothesis of the condensation model is consistent with the characteristics of the Poisson-Boltzmann potential at the surface of the polyion.     

Monte Carlo determination of the distribution of ions about a cylindrical polyelectrolyte

We report a calculation of the distribution of small ions around a charged cylinder representing a polyelectrolyte molecule in solution. The Monte Carlo method of Metropolis, Rosenbluth, and Teller was used to avoid the inaccuracies known to be associated with the Poisson-Boltzmann equation. The systems examined contained a long polyelectrolyte cylinder with charge parameter, χ, equal to 4.2, corresponding approximately to a DNA molecule. In one model, the cylinder had charges on its axis and an exclusion radius to the center of the small ions equal to 10 Å, while the small ions had various radii in the range from 1 to 10 Å and one or two protonic charges. Various systems were studied; some had one species of small ion alone, others had mixtures of different types. The results showed good agreement with the solution of the Poisson-Boltzmann equation when only the species with 1-Å radius was present, but considerable discrepancies appeared with larger ions as a result of excluded volume interactions between the latter. Deviations from the Poisson-Boltzmann equation also appeared when both positive and negative small ions were present; the deviations were in the direction of a higher concentration of both counter- and co-ions, but particularly co-ions, close to the polyelectrolyte. In another model, the charges were arranged along two helices on the surface of the cylinder; the resulting radial distribution of small ions was not much different from that found when the charges were situated on the axis. In all cases there was a striking accumulation of counterions in a layer of concentration exceeding 1 mol/L at the surface of the polyion.     

Force Balances in Systems of Cylindrical Polyelectrolytes

A detailed analysis is made of the model system of two parallel cylindrical polyelectrolytes which contain ionizable groups on their surfaces and are immersed in an ionic bathing medium. The interaction between the cylinders is examined by considering the interplay between repulsive electrostatic forces and attractive forces of electrodynamic origin. The repulsive force arises from the screened coulomb interaction between the surface charge distributions on the cylinders and has been treated by developing a solution to the linearized Poisson-Boltzmann equation. The boundary condition at the cylinder surfaces is determined as a self-consistent functional of the potential, with the input consisting of the density of ionizable groups and their dissociation constants. It is suggested that a reasonably accurate representation for the form of the attractive force can be obtained by performing a pairwise summation of the individual interatomic forces. A quantitative estimate is obtained using a Hamaker constant chosen on the basis of rigorous calculations on simpler systems. It is found that a balance exists between these repulsive and attractive forces at separations in good agreement with those observed in arrays of tobacco mosaic virus and in the A band myosin lattice in striated muscle. The behavior of the balance point as a function of the pH and ionic strength of the bathing medium closely parallels that seen experimentally.     

Comparison of hard-cylinder and screened coulomb interactions in the modeling of supercoiled DNAs

A 1000 base pair (bp) model supercoiled DNA is simulated using spherical screened Coulomb interactions between subunits on one hand and equivalent hard-cylinder interactions on the other. The amplitudes, or effective charges, of the spherical screened Coulomb electrostatic potentials are chosen so that the electrostatic potential surrounding the middle of a linear array of 2001 subunits (31.8 Å diameter) closely matches the solution of the nonlinear Poisson-Boltzmann equation for a cylinder with 12 Å radius and the full linear charge density of DNA at all distances beyond the 24 Å hard-core diameter. This superposition of spherical screened Coulomb potentials is practically identical to the particular solution of the cylindrical linearized Poisson-Boltzmann equation that matches the solution of the nonlinear Poisson-Boltzmann equation at large distances. The interaction energy between subunits is reckoned from the effective charges according to the standard DLVO expression. The equivalent hard-cylinder diameter is chosen following Stigter's protocol for matching second virial coefficients, but for the full linear charge density of DNA. The electrostatic persistence length of the model with screened Coulomb interactions is extremely sensitive to the (arbitrarily) chosen subunit length at the higher salt concentrations. The persistence length of the hard-cylinder model is adjusted to match that of the screened Coulomb model for each ionic condition. Simulations for a superhelix density σ = -0.05 using a spherical screened Coulomb interaction plus a 24 Å hard-cylinder core (SCPHC) potential indicate that the radius of gyration of this 1000 bp DNA actually undergoes a slight increase as the NaCl concentration is raised from 0.01 to 1.0M. Thus, merely softening the potential from hard-cylinder to screened Coulomb form does not produce a large decrease in radius of gyration with increasing NaCl concentration for DNAs of this size. Radii of gyration, static structure factors, and diffusion coefficients obtained using the equivalent hard-cylinder (EHC) potential agree well with those obtained using the SCPHC potential in 1.0M NaCl, but in 0.1M NaCl the agreement is not as good, and in 0.01M NaCl the agreement is definitely unsatisfactory. These conclusions differ in significant respects from those obtained in previous studies. © 1997 John Wiley & Sons, Inc. Biopoly 42: 455–470, 1997     

Electrostatic effects in proteins: comparison of dielectric and charge models.

Two approaches for calculating electrostatic effects in proteins are compared and ana analysis is presented of the dependence of calculated properties on the model used to define the charge distribution. Changes in electrostatic free energy have been calculated using a screened Coulomb potential (SCP) with a distance-dependent effective dielectric permittivity to model bulk solvent effects and a finite difference approach to solve the Poisson-Boltzmann (FDPB) equation. The properties calculated include shifts in dissociation constants of ionizable groups, the effect of annihilating surface charges on the binding of metals, and shifts in redox potentials due to changes in the charge of ionizable groups. In the proteins considered the charged sites are separated by 3.5-12 A. It is shown that for the systems studied in this distance range the SCP yields calculated values which are at least as accurate as those obtained from solution of the FDPB equation. In addition, in the distance range 3-5 A the SCP gives substantially better results than the FDPB equation. Possible sources of this difference between the two methods are discussed. Shifts in binding constants and redox potentials were calculated with several standard charge sets, and the resulting values show a variation of 20-40% between the 'best' and 'worst' cases. From this study it is concluded that in most applications, changes in electrostatic free energies can be calculated economically and reliably using an SCP approach with a single functional form of the screening function.