Aqueous solvation of As(OH)3: a Monte Carlo study with flexible polarizable classical interaction potentials.

A theoretical study of the hydration of arsenious acid is presented. This study included ab initio calculations and Monte Carlo simulations. The model potentials used for the simulations were ab initio derived and they include polarizability, nonadditivity, and molecular relaxation. It is shown that with these refined potentials it is possible to reproduce the available experimental evidence and therefore permit the study of clusters, as well as of the hydration process in solution. From the study of stepwise hydration and the Monte Carlo simulation of the condensed phase it is concluded that As(OH)(3) presents a hydration scheme similar to an amphipathic molecule. This phenomenon is explained as due to the existence of both a positive electrostatic potential and a localized lone pair in the vicinity of As. These results are used to rationalize the known passage of As(OH)(3) through aqua-glyceroporines.

Effect of membrane structure on the action of polyenes II: nystatin activity along the phase diagram of ergosterol- and cholesterol-containing POPC membranes.

Pores formed by the polyene antibiotic nystatin were studied in solvent-free lipid membranes. The membranes were formed by the tip-dip technique using 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) with different mol fractions (0-50%) of cholesterol or ergosterol. The effects of the mol fraction of sterol and of temperature variation (15-35°C) on the activity of the pores, their unitary conductances, lifetimes and time average conductances were studied. The results were used to analyze the behavior of nystatin channels along the phase diagrams previously reported for these lipid mixtures and to propose that membrane structure is the determinant factor for the known ergosterol/cholesterol selectivity.

Effect of membrane structure on the action of polyenes: I. Nystatin action in cholesterol- and ergosterol-containing membranes.

A detailed and thorough characterization of nystatin-induced permeability on lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-containing ergosterol or cholesterol is presented. The results show that the same collection of transmembrane pores appears in membranes containing either sterol. The concentration range for the appearance of these pores is sterol-dependent. Another mechanism of action, membrane disruption, is also observed in ergosterol-POPC membranes. The greater potency of nystatin present in ergosterol-containing membranes cannot be explained simply by the longer opening times of its pores, as has been suggested; it is also due to an increased number of events in these membranes. The present results and those of a companion paper lead us to propose that membrane structure is the determining factor for drug selectivity in membranes with different sterols.

On the role of sterol in the formation of the amphotericin B channel.

Amphotericin B is an antimycotic agent that has been studied for a long time, both because of its pharmacological action and the interest in understanding how this ionic channel works. It has been proposed that the channel is formed by a barrel of monomers, and that the presence of sterol is needed for the formation of such a barrel. As a matter of fact this need of a sterol has been used as a guiding idea in attempts to design derivatives more efficient in the discrimination of the cholesterol containing membranes, as compared to the ergosterol containing ones, henceforth diminishing the unwanted side effects in its pharmacological use. In this work we show that unitary channels that appear in a cholesterol containing membrane also appear when this membrane is free of cholesterol. We prove this to be the case for two membranes, a biological one, asolectin, and a synthetic one, DMPC. We then advance the idea that the role of sterols in the formation of the amphotericin B channel is related to the effects they have on the structure of the membrane itself, rather than to a direct involvement in the channel formation. We further look into the effect that different cholesterol concentrations in the membrane produce on the single channel properties.

The role of hydration in the hydrolysis of pyrophosphate. A Monte Carlo simulation with polarizable-type interaction potentials.

The exchange of energy in biochemical reactions involves, in a majority of cases, the hydrolysis of phosphoanhydrides (P-O-P). This discovery has lead to a long discussion about the origin of the high energy of such bonds, and to a proposal that hydration plays a major role in the energetics of the hydrolysis. This idea was supported by recent ab initio quantum mechanical calculations (Saint-Martin et al. (1991) Biochim. Biophys. Acta 1080, 205-214) that predicted the hydrolysis of pyrophosphate is exothermic in the gas phase. This exothermicity can account for only a half of the total energy release that one measures in aqueous solutions. Here we address the problem of hydration of the reactants and products of the pyrophosphate hydrolysis by means of Monte Carlo simulations, employing polarizable potentials whose parameters are fitted to energy surfaces computed at the SCF/6-31G** level of the theory. The present results show that the hydration enthalpies of the reactants and products contribute significantly to the total energy output of the pyrophosphate hydrolysis. The study predicts that both, the orthophosphate and the pyrophosphate, have hydration spheres with the water molecules acting as proton acceptors in the P-OH ... O(water) hydrogen bonds. These water molecules weakly repel the water molecules in the further hydration spheres. The perturbation of the structure of the solvent caused by the presence of the solute molecules is short ranged: after ca. 5 A from the P atoms, the energy and the structure of water correspond to bulk water. Due mainly to nonadditive effects, the molecular structure of the hydrated pyrophosphate is quite different from two fused structures of the hydrated orthophosphates. The hydration sphere of pyrophosphate is very loose and has a limited effect on the water network, whereas for orthophosphate it has a well developed shell structure. Hence, upon hydration there will be both a gain in hydration enthalpy and a gain in entropy because of distortion of the water molecular network.

Ab initio calculations of the pyrophosphate hydrolysis reaction.

Ab initio quantum mechanical calculations were used to study the hydrolysis reaction H4P2O7 + H2O in equilibrium with 2H3PO4, as well as some molecular properties of the reactants and products. SCF calculations with several basis sets ranging from minimal to extended with polarization functions were used to look at the basis dependency of the reaction enthalpies and optimized geometries. Although the minimal basis sets yield erratic predictions of the enthalpy, when a more extended basis (3-21G*) was used for the geometry optimization, and the total energies of the reactants and products were computed with this and larger basis sets, we obtained more consistent predictions of the structural properties of the P-O-P bridge and of the heat of the hydrolysis reaction (delta E = -7.39 kcal/mol at the SCF/6-31G** level). A comparison is made with previous estimates performed with smaller basis sets and without taking into account the electron correlation effects, which are calculated in the present work. The inclusion of the zero point energy calculated using the harmonic approximation, and of the electronic correlation energy determined at the MBPT(2) level, raised the computed heat of the reaction to -3.83 kcal/mol, and when an estimate for the thermal energy was added, the value obtained was of -3.38 kcal/mol. In conclusion, we found that the hydrolysis of pyrophosphate should be exothermic in the gas phase. The implications of this result in relation to some recent theories about enzyme catalysis are discussed.