Solution of the Poisson-Boltzmann equation for a system with four ionic species

The analytical solution of the Poisson-Boltzmann equation in an electrolyte with four ionic species (2:2:1:1), in the presence of a charged planar membrane or surface is presented. The function describing the mean electrical potential provides a convenient description that helps the understanding of electrical processes of biological interest.

INTEGRATION OF THE NONLINEAR POISSON BOLTZMANN-EQUATION FOR CHARGED VESICLES IN ELECTROLYTIC SOLUTIONS

The Poisson-Boltzmann equation (PBE), with specific ion-surface interactions and a cell model, was used to calculate the electrostatic properties of aqueous solutions containing vesicles of ionic amphiphiles. Vesicles are assumed to be water- and ion-permeable hollow spheres and specific ion adsorption at the surfaces was calculated using a Volmer isotherm. We solved the PBE numerically for a range of amphiphile and salt concentrations (up to 0.1 M) and calculated co-ion and counterion distributions in the inside and outside of vesicles as well as the fields and electrical potentials. The calculations yield results that are consistent with measured values for vesicles of synthetic amphiphiles.

SSPBE: um programa para solução numérica da equação de Poisson-Boltzmann em simetria esférica com modelo de adsorção

A Fortran77 program, SSPBE, designed to solve the spherically symmetric Poisson-Boltzmann equation using cell model for ionic macromolecular aggregates or macroions is presented. The program includes an adsorption model for ions at the aggregate surface. The working algorithm solves the Poisson-Boltzmann equation in the integral representation using the Picard iteration method. Input parameters are introduced via an ASCII file, sspbe.txt. Output files yield the radial distances versus mean field potentials and average molar ion concentrations, the molar concentration of ions at the cell boundary, the self-consistent degree of ion adsorption from the surface and other related data. Ion binding to ionic, zwitterionic and reverse micelles are presented as representative examples of the applications of the SSPBE program.

Aplicação da equação de Poisson-Boltzmann ao cálculo de propriedades dependentes do pH em proteínas

The ability of biomolecules to catalyze chemical reactions is due chiefly to their sensitivity to variations of the pH in the surrounding environment. The reason for this is that they are made up of chemical groups whose ionization states are modulated by pH changes that are of the order of 0.4 units. The determination of the protonation states of such chemical groups as a function of conformation of the biomolecule and the pH of the environment can be useful in the elucidation of important biological processes from enzymatic catalysis to protein folding and molecular recognition. In the past 15 years, the theory of Poisson-Boltzmann has been successfully used to estimate the pKa of ionizable sites in proteins yielding results, which may differ by 0.1 unit from the experimental values. In this study, we review the theory of Poisson-Boltzmann under the perspective of its application to the calculation of pKa in proteins.

Protein-Ion Binding Process on Finite Macromolecular Concentration. A Poisson-Boltzmann and Monte Carlo Study

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

pK determination. A mean field, Poisson-Boltzmann approach

A model describing dissociation of monoprotonic acid and a method for the determination of its pK value are presented. The model is based on a mean field approximation. The Poisson-Boltzmann equation, adopting spherical symmetry, is numerically solved, and the solution of its linearized form is written. By use of the pH values of a dilution experiment of galacturonic acid as the entry data, the proposed method allowed estimation of the value of pK = 3.25 at a temperature of 25 degrees C. Values for the complex dimensions and dissociation degree are calculated using experimental pH values for solution concentration values ranging from 0.1 to 60 mM. The present analysis leads to the conclusion that the Poisson-Boltzmann equation or its linear form is equally suited for the description of such systems.

pH at the micellar interface: Synthesis of pH probes derived from salicylic acid, acid-base dissociation in sodium dodecyl sulfate micelles, and Poisson-Boltzmann simulation

The study of the H+ concentration at the micellar interface is a convenient system for modeling the distribution of H+ at interfaces. We have synthesized salicylic acid derivatives to analyze the proton dissociation of both the carboxylic and phenol groups of' the probes, determining spectrophotometrically the apparent pK(a)'s (pK(ap)) in sodium dodecyl Sulfate, SDS, micelles with and without added salt. The synthesized probes were 2-hydroxy-5-(2-trimethylammoniumacetyl)benzoate; 2-hydroxy-5-(2-dimethylhexadecylammoniumacetyl)benzoate- 2-hydroxy-5-(2-dimethylhexadecylammoniumhexanoyl)benzoate-, 2-hydroxy-5-(2-diniethylhexadecylammoniumundecanoyl)betizoate; 2-hydroxy-5-acetylbenzoic acids and 2-hydroxy-5-dodecanoylbenzoic acid. Upon incorporation into SDS micelles the pK(ap)'s of both carboxylic and phenol groups increased by ca. 3 pH units and NaCl addition caused a decrease in the probe-incorporated pKap. The experimental results were fitted with a cell model Poisson-Boltzmann (P-B) equation taking in consideration the effect of salt on the aggregation number of SDS and using the distance of' the dissociating group as a parameter. The conformations of the probes were analyzed theoretically using two dielectric constants, e.g., 2 and 78. Both the P-B analysis and conformation calculations can be interpreted by assuming that the acid groups dissociate very close to, or at, the interface. Our results are consistent with the assumption that the intrinsic pK(a)'s of both carboxylic and phenol groups of the salicylic acid probes used here can be taken as those in water. Using this assumption the micellar and salt effects on the pKap's of the (trialkylammonium)benzoate probes were described accurately using a cell model P-B analysis. (c) 2005 Elsevier B.V. All rights reserved.

Contractive metrics for a Boltzmann equation for granular gases: diffusive equilibriaThe Patterson-Sullivan embedding and minimal volume entropy for outer space

We quantify the long-time behavior of a system of (partially) inelastic particles in a stochastic thermostat by means of the contractivity of a suitable metric in the set of probability measures. Existence, uniqueness, boundedness of moments and regularity of a steady state are derived from this basic property. The solutions of the kinetic model are proved to converge exponentially as t→ ∞ to this diffusive equilibrium in this distance metrizing the weak convergence of measures. Then, we prove a uniform bound in time on Sobolev norms of the solution, provided the initial data has a finite norm in the corresponding Sobolev space. These results are then combined, using interpolation inequalities, to obtain exponential convergence to the diffusive equilibrium in the strong L¹-norm, as well as various Sobolev norms.

A equação de Poisson-Boltzmann em regiões com fronteira irregular

Propomos uma idealização da situação em que uma macromolécula é ionizada em um solvente. Neste modelo a área da superfície da molécula é suposta ser grande com respeito a seu diâmetro. A molécula é considerada como um dielétrico com uma distribuição de cargas em sua superfície. Utilizando as condições de transmissão, a distribuição de Boltzmann no solvente e resultados recentes sobre espaços de Sobolev no contexto de espaços métricos, bem como de integração sobre superfícies irregulares, o problema é formulado em forma variacional. Resultados clássicos do cálculo de variações permitem a resolução analítica do problema.

Aplicação da equação de Paisson-Boltzmann com condição de contorno para regulação carga-potencial a sistemas polieletrolíticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Iterative and direct methods for solving Poisson's equation and their adaptability to ILLIAC IV,

Includes bibliography.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.

Influence of fixed electric charges on potential profile across the squid axon membrane

The potential profile for a model of squid axon membrane has been determined for two physiological states: resting and action states. The non-linear Poisson-Boltzmann equation has been solved by considering the volumetric charge densities due to charges dissolved in an electrolytic solution and fixed on both glycocalyx and cytoplasmatic proteins. Results showing the features of the potential profile along the outer electrolytic region are similar for both resting and action states. However, the potential fall along glycocalyx at action state is lower than at resting. A small variation in the Na+ concentration drastically affects the surface membrane potentials and vice versa. We conclude that effects on the potential profile due to surface lipidic bilayer charge and contiguous electric double layers are more relevant than those provoked by fixed charges distributed along the cell cytoplasm. (c) 2007 Elsevier B.V. All rights reserved.