Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “electrostatics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources

Particle deposition onto a human head: influence of electrostatic and wind fields.

This study investigates electrostatic fields surrounding the human head and particle deposition onto facial skin and eyes caused by the combined effect of electrostatic and wind fields. The electrostatic fields are calculated by a three-dimensional numerical model calculating the field strength between a field source and a human head. The deposition velocity can be viewed as determined by the sum of two contributions: that of an electrostatic field and that of a wind field. Deposition velocities are calculated by a semiempirical particle deposition model that considers particle transport from the free stream to the human face. The particle deposition model uses the electrostatic field model results as input parameters and is applied to the forehead and eyes of two facial shapes for two different turbulence conditions and aerosol charge distributions. The results of different practical working conditions, under which the potential difference between head (person) and source ranges from 5.6 to 15.0 kV, indicates that the presence of electrostatic fields always increases particle deposition for industrial aerosols. For aged aerosols an effect is only present for submicron particles.

Air↗

Classification of protein sequences by homology modeling and quantitative analysis of electrostatic similarity.

Protein electrostatics plays a key role in ligand binding and protein-protein interactions. Therefore, similarities or dissimilarities in electrostatic potentials can be used as indicators of similarities or dissimilarities in protein function. We here describe a method to compare the electrostatic properties within protein families objectively and quantitatively. Three-dimensional structures are built from database sequences by comparative modeling. Molecular potentials are then computed for these with a continuum solvation model by finite difference solution of the Poisson-Boltzmann equation or analytically as a multipole expansion that permits rapid comparison of very large datasets. This approach is applied to 104 members of the Pleckstrin homology (PH) domain family. The deviation of the potentials of the homology models from those of the corresponding experimental structures is comparable to the variation of the potential in an ensemble of structures from nuclear magnetic resonance data or between snapshots from a molecular dynamics simulation. For this dataset, the results for analysis of the full electrostatic potential and the analysis using only monopole and dipole terms are very similar. The electrostatic properties of the PH domains are generally conserved despite the extreme sequence divergence in this family. Notable exceptions from this conservation are seen for PH domains linked to a Db1 homology (DH) domain and in proteins with internal PH domain repeats.

Data Interpretation, Statistical↗

Hydration effects on the electrostatic potential around tuftsin.

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr-Lys-Pro-Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson-Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed.

Computer Simulation↗

Characterization of electrostatic binding sites of extracellular polymers by linear programming analysis of titration data.

Electrostatic binding sites of extracellular polymeric substances (EPS) were characterized from titration data using linear programming analysis. Test results for three synthetic solutions of given solutes comprising amino, carboxyl, and phenolic groups indicated that this method was able to identify the electrostatic binding sites. For the six sites with pK(a) between 3 and 10, the estimated pK(a) deviated 0.11 +/- 0.09 from the theoretical values, and the estimated concentrations deviated 3.0% +/- 0.9% from the actual concentrations. Two EPS samples were then extracted from a hydrogen-producing sludge (HPS) and a sulfate-reducing biofilm (SRB). Analysis of charge excess data in titration from pH 3 to 11 indicated that the EPS of HPS comprised of five electrostatic binding sites with pK(a) ranging from 3 to 11. The pK(a) values of these binding sites and the possible corresponding functional groups were pK(a) 4.8 (carboxyl), pK(a) 6.0 (carboxyl/phosphoric), pK(a) 7.0 (phosphoric), pK(a) 9.8 (amine/phenolic), and pK(a) 11.0 (hydroxyl). EPS of the SRB comprised five of similar binding sites (with corresponding pK(a) values of 4.4, 6.0, 7.4, 9.4, and 11.0), plus one extra site at pK(a) 8.2, which was likely corresponding to the sulfhydryl group. The total electrostatic binding site concentration of EPS extracted from HPS were 10.88 mmol/g-EPS, of which the highest concentration was from the site of pK(a) 11.0. The corresponding values for the EPS extracted from SRB were 16.44 mmol/g-EPS and pK(a) 4.4. The total concentrations of electrostatic binding sites found in this study were 20- to 30-fold of those reported for bacterial cell surface, implying that EPS might be more crucial in biosorption of metals than bacterial cell surface in wastewater treatment and in bioremediation.

Bacteria↗

Optimization of the electrostatic interactions in proteins of different functional and folding type.

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed alpha beta folding type are electrostatically better optimized than pure alpha-helical or beta-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.

Computer Simulation↗

Electrostatic contribution to the binding stability of protein-protein complexes.

To investigate roles of electrostatic interactions in protein binding stability, electrostatic calculations were carried out on a set of 64 mutations over six protein-protein complexes. These mutations alter polar interactions across the interface and were selected for putative dominance of electrostatic contributions to the binding stability. Three protocols of implementing the Poisson-Boltzmann model were tested. In vdW4 the dielectric boundary between the protein low dielectric and the solvent high dielectric is defined as the protein van der Waals surface and the protein dielectric constant is set to 4. In SE4 and SE20, the dielectric boundary is defined as the surface of the protein interior inaccessible to a 1.4-A solvent probe, and the protein dielectric constant is set to 4 and 20, respectively. In line with earlier studies on the barnase-barstar complex, the vdW4 results on the large set of mutations showed the closest agreement with experimental data. The agreement between vdW4 and experiment supports the contention of dominant electrostatic contributions for the mutations, but their differences also suggest van der Waals and hydrophobic contributions. The results presented here will serve as a guide for future refinement in electrostatic calculation and inclusion of nonelectrostatic effects.

Acetylcholinesterase↗

Electrostatic fields at the active site of ribulose-1,5-bisphosphate carboxylase.

A macroscopic approach has been employed to calculate the electrostatic potential field of nonactivated ribulose-1,5-bisphosphate carboxylase and of some complexes of the enzyme with activator and substrate. The overall electrostatic field of the L2-type enzyme from the photosynthetic bacterium Rhodospirillum rubrum shows that the core of the dimer, consisting of the two C-terminal domains, has a predominantly positive potential. These domains provide the binding sites for the negatively charged phosphate groups of the substrate. The two N-terminal domains have mainly negative potential. At the active site situated between the C-terminal domain of one subunit and the N-terminal domain of the second subunit, a large potential gradient at the substrate binding site is found. This might be important for polarization of chemical bonds of the substrate and the movement of protons during catalysis. The immediate surroundings of the activator lysine, K191, provide a positive potential area which might cause the pK value for this residue to be lowered. This observation suggests that the electrostatic field at the active site is responsible for the specific carbamylation of the epsilon-amino group of this lysine side chain during activation. Activation causes a shift in the electrostatic potential at the position of K166 to more positive values, which is reflected in the unusually low pK of K166 in the activated enzyme species. The overall shape of the electrostatic potential field in the L2 building block of the L8S8-type Rubisco from spinach is, despite only 30% amino acid homology for the L-chains, strikingly similar to that of the L2-type Rubisco from Rhodospirillum rubrum. A significant difference between the two species is that the potential is in general more positive in the higher plant Rubisco. In particular, the second phosphate binding site has a considerably more positive potential, which might be responsible for the higher affinity for the substrate of L8S8-type enzymes. The higher potential at this site might be due to two remote histidine residues, which are conserved in the plant enzymes.

Binding Sites↗

The role of arginine 143 in the electrostatics and mechanism of Cu,Zn superoxide dismutase: computational and experimental evaluation by mutational analysis.

Cu,Zn superoxide dismutase protects cells from oxidative damage by removing superoxide radicals in one of the fastest enzyme reactions known. The redox reaction at the active-site Cu ion is rate-limited by diffusion and enhanced by electrostatic guidance. To quantitatively define the electrostatic and mechanistic contributions of sequence-invariant Arg-143 in human Cu,Zn superoxide dismutase, single-site mutants at this position were investigated experimentally and computationally. Rate constants for several Arg-143 mutants were determined at different pH and ionic strength conditions using pulse radiolytic methods and compared to results from Brownian dynamics simulations. At physiological pH, substitution of Arg-143 by Lys caused a 2-fold drop in rate, neutral substitutions (Ile, Ala) reduced the rate about 10-fold, while charge-reversing substitutions (Asp, Glu) caused a 100-fold decrease. Position 143 mutants showed pH dependencies not seen in other mutants. At low pH, the acidic residue mutations exhibited protonation/deprotonation effects. At high pH, all enzymes showed typical decreases in rate except the Lys mutant in which the rate dropped off at an unusually low pH. Increasing ionic strength at acidic pH decreased the rates of the wild-type enzyme and Lys mutant, while the rate of the Glu mutant was unaffected. Increasing ionic strength at higher pH (> 10) increased the rates of the Lys and Glu mutants while the rate of the wild-type enzyme was unaffected. Reaction simulations with Brownian dynamics incorporating electrostatic effects tested computational predictability of ionic strength dependencies of the wild-type enzyme and the Lys, Ile, and Glu mutants. The calculated and experimental ionic strength profiles gave similar slopes in all but the Glu mutant, indicating that the electrostatic attraction of the substrate is accurately modeled. Differences between the calculated and experimental rates for the Glu and Lys mutants reflect the mechanistic contribution of Arg-143. Results from this joint analysis establish that, aside from the Cu ligands, Arg-143 is the single most important residue in Cu,Zn superoxide dismutase both electrostatically and mechanistically, and provide an explanation for the evolutionary selection of arginine at position 143.

Arginine↗

Proton Binding to Humic Acids: Electrostatic and Intrinsic Interactions.

The proton adsorption behavior of eight humic acids (HA) and a fulvic acid (FA) was studied as a function of pH and KNO(3) concentration. The emphasis is on the comparison of the different humics with respect to their ion binding properties and on the comparison of two different models to describe the electrostatic interactions: the Donnan model and the impermeable sphere (IS) model. Viscosimetric data were used to estimate the hydrodynamic volumes and radii of the HA molecules. These data were incorporated in the electrostatic models and calculations could be carried out without any adjustable parameter. The Donnan model in combination with hydrodynamic volumes obtained by viscometry cannot adequately describe the electrostatic effects related to changes of the electrolyte concentration. This model leads to good prediction of the HA behavior if unrealistically large volumes are used for fulvics and unrealistically large volume-salt concentration dependencies are used for humics. The IS model can successfully reproduce experimental proton adsorption data with physically realistic radii. The good performance of the IS model and the poor performance of the Donnan model is directly related to the fact that the hydrodynamic volumes of the molecules are too small to allow for charge compensation within the molecular limits. The combination of viscometry with the IS model leads to a consistent description of the electrostatic in humics and to a consistent way of positioning the master curves. Therefore, the electrostatic potentials and the intrinsic affinity distributions of the different samples can be compared on an equal basis. The similarities in the intrinsic affinity distributions give faith for the possibility to develop a generic model to describe the ion binding to humics. Copyright 1999 Academic Press.

Journal Article↗

Collision Kinetics and Electrostatic Dispersion of Airborne Submicrometer Fractal Agglomerates.

Collision and electrostatic dispersion rates of airborne submicrometer TiO(2) agglomerates were measured and compared with the classical collision theory for spheres as well as with models accounting for the agglomerate structure in terms of the fractal dimension and electrostatic effects such as Coulomb and van der Waals interactions. According to the authors' knowledge, this is the first time that the agglomerate fractal dimension and electrostatic effects have been considered simultaneously in determining the collision frequency function of agglomerates. The observed enhancement in the collision frequency of agglomerates was found mainly to be a result of electrostatic particle interactions. Nonspherical particle shape has only a comparatively small influence on the collision probability, on the order of 10-20%. Electrostatic dispersion coefficients of agglomerates were found to be similar to those of spheres. Copyright 2001 Academic Press.

Journal Article↗

Effect of electrostatic interaction on the adsorption of globular proteins on octacalcium phosphate crystal film.

The electrostatic effect on the adsorption of globular proteins, such as bovine serum albumin (BSA), hen egg white lysozyme (LZM), and beta-lactoglobulin (beta-Lg), on octacalcium phosphate (OCP)-like crystal thin films was investigated. A poorly crystalline thin film was synthesized on a tissue culture polystyrene (TCP) surface and used as a model surface in this study. The solution pH clearly affected the electrostatic properties of both proteins and surface. The adsorbed amounts obtained at quasi-steady state were readily related to the solution pH for each protein. The adsorption rate is fast during the initial period and levels off gradually. The maximum adsorbed mass occurred at pH 7 for BSA and at pH 9 for LZM. beta-Lg adsorbed similar amounts at pHs lower than 9, but the adsorbed mass decreased at pHs higher than 9 where electrostatic repulsion exists. The pH values where the maximum adsorbed mass occurred may be considered as the conditions where electrostatic attraction is most favorable. The adsorbed mass of beta-Lg was the greatest among the proteins of interest while BSA adsorbed the least despite its greater molecular mass. LZM falls into the intermediate region. According to these observations, BSA has undergone conformational changes that prevent further adsorption to a greater extent than the others. A simple relationship between the adsorption rate and the electrostatic properties was not established. However, the order of magnitude of the adsorption rate at the initial period tends to be the same as that of maximum adsorbed mass for each protein.

Adsorption↗

Electrostatic channeling in the bifunctional enzyme dihydrofolate reductase-thymidylate synthase.

The bifunctional enzyme dihydrofolate reductase-thymidylate synthase (DHFR-TS) carries out two distinct reactions, with the dihydrofolate produced by the TS-catalyzed reaction acting as the substrate for the DHFR-catalyzed reaction. Brownian dynamics simulation techniques were used to investigate the possible role of electrostatics in determining efficient channeling of the substrate, by explicitly simulating substrate diffusion between the two active sites. With a substrate charge of -2, almost all (> 95%) substrate molecules leaving the TS active site reached the DHFR active site at zero ionic strength. Under the same conditions, but in the absence of electrostatic effects, successful channeling was reduced to only around 6%: electrostatic effects therefore appear essential to explain the efficient channeling observed experimentally. The importance of substrate charge, the relative contributions of specific basic residues in the protein, the role played by the second monomer of the dimer in channeling and the effects of changing ionic strength were all investigated. Simulations performed for substrate transfer in the opposite direction suggest that channeling in DHFR-TS is not strongly directional and that the role of electrostatics is perhaps more one of restricting diffusion of the substrate than one of actively guiding it from the TS to the DHFR active site. The results demonstrate that electrostatic channeling can be a highly efficient means of transferring charged substrates between active sites in solvent-exposed environments.

Binding Sites↗

Electrostatic effects in homeodomain-DNA interactions.

We report here an investigation of the role of electrostatics in homeodomain-DNA interactions using techniques based around the use of the Poisson-Boltzmann equation. In the present case such a study is of particular interest, since in contrast to other proteins previously studied with this method, the homeodomain is a small, highly charged protein that forms extensive ion pairs upon binding DNA. We have investigated the salt dependence of the binding constant for specific association and for a variety of models for non-specific association. The results indicate that, in line with the models proposed by Manning and Record, the entropy of counterion release accounts for a significant fraction of the salt dependence of the binding free energy, though this is perhaps due to fortuitous cancellation of other contributing terms. The thermodynamic effects of a number of specific homeodomain mutants were also investigated, and partly rationalized in terms of favorable electrostatic interactions in the major goove of DNA. Investigation of the temperature-dependence of the free energy of association indicates that the electrostatic contributions become increasingly favorable as the temperature rises. For this particular system, however, there appears to be no significant electrostatic contribution to the delta(delta C(p)) of association. Finally, an analysis of the free energy of interaction when the homeodomain is moved ca one Debye length from the DNA suggests that pure electrostatic forces are able to steer the homeodomain into a partially correct orientation for binding to the DNA.

Amino Acid Sequence↗

Role of main-chain electrostatics, hydrophobic effect and side-chain conformational entropy in determining the secondary structure of proteins.

The physiochemical bases of amino acid preferences for alpha-helical, beta-strand, and other main-chain conformational states in proteins is controversial. Hydrophobic effect, side-chain conformational entropy, steric factors, and main-chain electrostatic interactions have all been advanced as the dominant physical factors which determine these preferences. Many attempts to resolve the controversy have focused on small model systems. The disadvantage of such systems is that the amino acids in small molecules are largely exposed to the solvent. In proteins, however, the amino acids are in contact with the solvent to a different degree, causing a large variability of strengths of all interactions. The estimates of mean strengths of interactions in the actual protein environment are therefore essential to resolve the controversy. In this work the experimental protein structures are used to estimate the mean strengths of various interactions in proteins. The free energy contributions of the interactions are implemented into the Lifson-Roig theory to calculate the helix and strand free energy profiles. From the profiles the secondary structures of proteins and peptides are predicted using simple rules. The role of hydrophobic effect, side-chain conformational entropy, and main-chain electrostatic interactions in determining the secondary structure of proteins is assessed from the abilities of different models, describing stability of secondary structures, to correctly predict alpha-helices, beta-strands and coil in 130 proteins. The three-state accuracy of the model, which contains only the free energy terms due to the main-chain electrostatics with 40 coefficients, is 68.7%. This accuracy is approaching to the accuracy of currently the best secondary structure prediction algorithm based on neural networks (72%); however, many thousands of parameters have to be optimized during the training of the neural networks to reach this level of accuracy. The correlation coefficient between the calculated and the experimental helix contents of 37 alanine based peptides is 0.91. If the hydrophobic and the side-chain conformational entropy terms are included into the helix-coil transition parameters, the accuracy of the algorithm does not improve significantly. However, if the main-chain electrostatic interactions are excluded from the helix-coil and strand-coil transition parameters, the accuracy of the algorithm reaches only 59.5%. These results support the dominant role of the short-range main-chain electrostatics in determining the secondary structure of proteins and peptides. The role of the hydrophobic effect and the side-chain conformational entropy is small.

Agglutinins↗

Electrostatic contributions to the stability of halophilic proteins.

Examination of the first crystal structures of proteins from a halophilic organism suggests that an abundance of acidic residues distributed over the protein surface is a key determinant of adaptation to high-salt conditions. Although one extant theory suggests that acidic residues are favored because of their superior water-binding capacity, it is clear that extensive repulsive electrostatic interactions will also be present in such proteins at physiological pH. To investigate the magnitude and importance of such electrostatic interactions, we conducted a theoretical analysis of their contributions to the salt and pH-dependence of stability of two halophilic proteins. Our approach centers on use of the Poisson-Boltzmann equation of classical electrostatics, applied at an atomic level of detail to crystal structures of the proteins. We first show that in using the method, it is important to account for the fact that the dielectric constant of water decreases at high salt concentrations, in order to reproduce experimental changes in pKa values of small acids and bases. We then conduct a comparison of salt and pH effects on the stability of 2Fe-2S ferredoxins from the halophile Haloarcula marismortui and the non-halophile anabaena. In both proteins, substantial upward shifts in pKa accompany protein folding, though shifts are considerably larger, on average, in the halophile. Upward shifts for basic residues occur because of favorable salt-bridge interactions, whilst upward shifts for acidic residues result from unfavorable electrostatic interactions with other acidic groups. Our calculations suggest that at pH 7 the stability of the halophilic protein is decreased by 18.2 kcal/mol on lowering the salt concentration from 5 M to 100 mM, a result that is in line with the fact that halophilic proteins generally unfold at low salt concentrations. For comparison, the non-halophilic ferredoxin is calculated to be destabilized by only 5.1 kcal/mol over the same range. Analysis of the pH stability curve suggests that lowering the pH should increase the intrinsic stability of the halophilic protein at low salt concentrations, although in practice this is not observed because of aggregation effects. We report the results of a similar analysis carried out on the tetrameric malate dehydrogenase from H. marismortui. In this case, we investigated the salt and pH dependence of the various monomer-monomer interactions present in the tetramer. All monomer-monomer interactions are found to make substantial contributions to the salt-dependence of stability of the tetramer. Excellent agreement is obtained between our calculated results for the stability of the tetramer and experimental results. In particular, the finding that at 4 M NaCl, the tetramer is stable only between pH 4.8 and 10 is accurately reproduced. Taken together, our results suggest that repulsive electrostatic interactions between acidic residues are a major factor in the destabilization of halophilic proteins in low-salt conditions, and that these interactions remain destabilizing even at high salt concentrations. As a consequence, the role of acidic residues in halophilic proteins may be more to prevent aggregation than to make a positive contribution to intrinsic protein stability.

Archaeal Proteins↗

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.

Amino Acid Sequence↗

Balance of electrostatic and hydrophobic interactions in the lysis of model membranes by E. coli alpha-haemolysin.

The relative weight of electrostatic interactions and hydrophobic forces in the process of membrane disruption caused by E. coli alpha-haemolysin (HlyA) has been studied with a purified protein preparation and a model system consisting of large unilamellar vesicles loaded with water-soluble fluorescent probes. Vesicles were prepared in buffers of different ionic strengths, or pHs, and the net surface charge of the bilayers was also modified by addition of negatively (e.g., phosphatidylinositol) or positively (e.g., stearylamine) charged lipids. The results can be interpreted in terms of a multiple equilibrium in which alpha-haemolysin may exist: aggregated HlyA <==> monomeric HlyA <==> membrane-bound HlyA. In these equilibria both electrostatic and hydrophobic forces are significant. Electrostatic forces become substantial under certain circumstances, e.g., membrane binding when bilayer and protein have opposite electric charges. Protein adsorption to the bilayer is more sensitive to electrostatic forces than membrane disruption itself. In the latter case, the irreversible nature of protein insertion may overcome electrostatic repulsions. Also of interest is the complex effect of pH on the degree of aggregation of an amphipathic toxin like alpha-haemolysin, since pH changes are not only influencing the net protein charge but may also be inducing protein conformational transitions shown by changes in the protein intrinsic fluorescence and in its susceptibility to protease digestion, that appear to regulate the presence of hydrophobic patches at the surface of the molecule, thus modifying the ability of the toxin to either aggregate or become inserted in membranes.

Bacterial Proteins↗

Electrostatic interactions across a charged lipid bilayer.

We present theoretical work in which the degree of electrostatic coupling across a charged lipid bilayer in aqueous solution is analyzed on the basis of nonlinear Poisson-Boltzmann theory. In particular, we consider the electrostatic interaction of a single, large macroion with the two apposed leaflets of an oppositely charged lipid bilayer where the macroion is allowed to optimize its distance to the membrane. Three regimes are identified: a weak and a high macroion charge regime, separated by a regime of close macroion-membrane contact for intermediate charge densities. The corresponding free energies are used to estimate the degree of electrostatic coupling in a lamellar cationic lipid-DNA complex. That is, we calculate to what extent the one-dimensional DNA arrays in a sandwich-like lipoplex interact across the cationic membranes. We find that, in spite of the low dielectric constant inside a lipid membranes, there can be a significant electrostatic contribution to the experimentally observed cross-bilayer orientational ordering of the DNA arrays. Our approximate analytical model is complemented and supported by numerical calculations of the electrostatic potentials and free energies of the lamellar lipoplex geometry. To this end, we solve the nonlinear Poisson-Boltzmann equation within a unit cell of the lamellar lipoplex using a new lattice Boltzmann method.

Cations↗