Modulatory Effects of Acidic pH and Membrane Potential on the Adsorption of pH-Sensitive Peptides to Anionic Lipid Membrane.

Anionic lipid membrane electrostatic potential and solution pH can influence cationic peptide adsorption to these bilayers, especially those containing simultaneously acid and basic residues. Here, we investigate the effects of the pH solution on MP1 (IDWKKLLDAAKQIL-NH2) adsorption to anionic (7POPC:3POPG) lipid vesicles in comparison to its analog H-MP1, with histidines substituting lysines. We used the association of adsorption isotherms and constant pH molecular dynamic simulations (CpHMD) to explore the effects of membrane potential and pH on peptides' adsorption on this lipid membrane. We analyzed the fluorescence and zeta potential adsorption isotherms using the Gouy-Chapman theory. In CpHMD simulations for the peptides in solution and adsorbed on the lipid bilayer, we used the conformations obtained by conventional MD simulations at a µs timescale. Non-equilibrium Monte Carlo simulations provided the protonation states of acidic and basic residues. CpHMD showed average pKa shifts of two to three units, resulting in a higher net charge for the analog than for MP1, strongly modulating the peptide adsorption. The fractions of the protonation of acidic and basic residues and the peptides' net charges obtained from the analysis of the adsorption isotherms were in reasonable agreement with those from CpHMD. MP1 adsorption was almost insensitive to solution pH. H-MP1 was much more sensitive to partitioning, at acidic pH, with an affinity ten times higher than in neutral ones.

Adsorption of polyelectrolytes onto oppositely charged cylindrical macroions.

In this work we investigate the adsorption of polyelectrolyte chains onto uniformly charged cylindrical macroions by means of the Metropolis Monte Carlo simulations and weighted histogram analysis method. Adopting a simplified model for macromolecules and treating the electrolytic solution in the Debye-Hückel level, conformational properties of the adsorbed chain, such as the radius of gyration and the thickness of the adsorbed layer, are provided as a function of ionic strength and macroion charge density. By analysis of the free energy profile as a function of the radius of gyration it was possible to identify first-order-like transitions between adsorbed and desorbed states and obtain a macroion charge density dependence of the critical ionic strength in good agreement with experiments.

Protein-ion binding process on finite macromolecular concentration. A Poisson-Boltzmann and Monte Carlo study.

Electrostatic interactions are one of the key driving forces for protein-ligands complexation. Different levels for the theoretical modeling of such processes are available on the literature. Most of the studies on the Molecular Biology field are performed within numerical solutions of the Poisson-Boltzmann Equation and the dielectric continuum models framework. In such dielectric continuum models, there are two pivotal questions: (a) how the protein dielectric medium should be modeled, and (b) what protocol should be used when solving this effective Hamiltonian. By means of Monte Carlo (MC) and Poisson-Boltzmann (PB) calculations, we define the applicability of the PB approach with linear and nonlinear responses for macromolecular electrostatic interactions in electrolyte solution, revealing some physical mechanisms and limitations behind it especially due the raise of both macromolecular charge and concentration out of the strong coupling regime. A discrepancy between PB and MC for binding constant shifts is shown and explained in terms of the manner PB approximates the excess chemical potentials of the ligand, and not as a consequence of the nonlinear thermal treatment and/or explicit ion-ion interactions as it could be argued. Our findings also show that the nonlinear PB predictions with a low dielectric response well reproduce the pK shifts calculations carried out with an uniform dielectric model. This confirms and completes previous results obtained by both MC and linear PB calculations.

Monte Carlo and modified Tanford-Kirkwood results for macromolecular electrostatics calculations.

The understanding of electrostatic interactions is an essential aspect of the complex correlation between structure and function of biological macromolecules. It is also important in protein engineering and design. Theoretical studies of such interactions are predominantly done within the framework of Debye-Hückel theory. A classical example is the Tanford-Kirkwood (TK) model. Besides other limitations, this model assumes an infinitesimally small macromolecule concentration. By comparison to Monte Carlo (MC) simulations, it is shown that TK predictions for the shifts in ion binding constants upon addition of salt become less reliable even at moderately macromolecular concentrations. A simple modification based on colloidal literature is suggested to the TK scheme. The modified TK models suggested here satisfactorily predict MC and experimental shifts in the calcium binding constant as a function of protein concentration for the calbindin D(9k) mutant and calmodulin.