Molecular electrostatic potentials from crystal diffraction: the neurotransmitter gamma-aminobutyric acid.

Given the electronic charge parameters obtained from a diffraction study of the charge density distribution in a crystal, a mathematical procedure is presented for deriving the electrostatic potential. The procedure allows the mapping of electrostatic potential for a molecule or group of molecules removed from the crystal structure but with each molecule retaining the effects of polarization owing to the original crystal environment. The method is applied for the neurotransmitter gamma-aminobutyric acid. The potential for a gamma-aminobutyric acid molecule is analyzed in terms of a simple model that is suitable for rapid computations concerned with Coulombic molecular interactions. Outside the molecular envelope at 1.2 A from the atomic nuclei, the total potential is well represented by a sum of spherical atom contributions with V(r) = (q/r) + exp(-beta r2). The most important aspherical component in the potential is the dipole contribution from the hydrogen atoms. This can be represented as V(r, phi) = (0.162 cos phi)/(r2 + 0.615). Here, V is in e/A, r is the distance from each nucleus in A, q is the experimentally determined net atomic charge in electron units, and phi is the angle between r and the bond X-H. We obtain beta = 1.47, 1.66, 1.83 A-2 for C, N, and O respectively. For H, no term in beta is needed.     

Membrane dipole potentials, hydration forces, and the ordering of water at membrane surfaces.

We have compared hydration forces, electrical dipole potentials, and structural parameters of dispersions of dipalmitoylphosphatidylcholine (DPPC) and dihexadecylphosphatidylcholine (DHPC) to evaluate the influence of fatty acid carbonyl groups on phospholipid bilayers. NMR and x-ray investigations performed over a wide range of water concentrations in the samples show, that in the liquid crystalline lamellar phase, the presence of carbonyl groups is not essential for lipid structure and hydration. Within experimental error, the two lipids have identical repulsive hydration forces between their bilayers. The higher transport rate of the negatively charged tetraphenylboron over the positively charged tetraphenylarsonium indicates that the dipole potential is positive inside the membranes of both lipids. However, the lack of fatty acid carbonyl groups in the ether lipid DHPC decreased the potential by (118 +/- 15) mV. By considering the sign of the potential and the orientation of carbonyl groups and headgroups, we conclude that the first layer of water molecules at the lipid water interface makes a major contribution to the dipole potential.     

Dynamical properties of phospholipid bilayers from computer simulation

We present the results of a 10-ns molecular dynamics simulation of a dipalmitoylphosphatidylcholine/water system. The main emphasis of the present study is on the investigation of the stability over a long time and the dynamic properties of the water/membrane system. The motion of the lipid molecules is characterized by the center of mass movement and the displacement of individual atom groups. Because of the slow movement of the headgroup atoms, their contributions to the dipole potential vary slowly and with a large amplitude. Nevertheless, the water molecules compensate the strong fluctuations and maintain an almost constant total dipole potential. From the lateral displacement of the center of masses, we calculate the lateral diffusion coefficient to be Dlat = (3 +/- 0.6) x 10(-7) cm2/s, in agreement with neutron scattering results. The rotational motion is also investigated in our simulations. The calculated value for the rotational diffusion coefficient parallel to the molecular long axis, D = (1.6 +/- 0.1) x 10(8) s-1, is in good agreement with the experiment.     

Fullerene modified supported lipid membrane as sensitive element of sensor for odorants

A supported lipid bilayer membrane (s-BLMs) formed on a freshly cleaved metallic surface by the Tien method was applied for the design of an electrochemical sensor for detection of neutral odorant molecules. The lipid bilayer was modified by saturation with fullerene C₆₀, which possesses electron mediator properties and facilitates a redox reaction occurring at the border of the lipid membrane and metal surface. I₂/I⁻ and ferrocenyl trimethyl bromide were used as electroactive marker ions. The smell compounds adsorb on the surface of the lipid layer and change its structure. As a consequence the ratio of marker ion penetration to the lipid membrane is altered. The magnitude of these changes depends on the amount and chemical structure of adsorbed molecules. The research presented was carried out by cyclic voltammetry. The magnitude of the electrochemical signal generated by smell compounds was correlated with other parameters describing their molecular properties such as: octanol/water partition coefficients and dipole moments.     

Effect of trehalose and sucrose on the hydration and dipole potential of lipid bilayers

The water activity in dimyristoylphosphatidylcholine (DMPC) decreases by 60% when the lipid is dehydrated in the presence of trehalose concentrations higher than 0.02 M. In contrast, sucrose in concentrations 10 times higher produced only a 20% decrease in the water activity in the sample. Titrations of a DMPC solution in chloroform yielded 14 water molecules per lipid when pure water was added and seven water molecules per lipid when the titration was done with 0.025 M trehalose. The same concentrations of sucrose produced a turbid solution, which made it impossible to quantify the number of water molecules per lipid. Lipid monolayers spread on an air/water interface showed a decrease from 480 mV in pure water to 425 mV in 0.1 M trehalose. However, the same concentrations of sucrose produced an increase of less than 100 mV. Results obtained with Fourier transform infrared spectroscopy (FTIR) under the same conditions denoted that trehalose binds to the carbonyl groups, while sucrose showed no specific binding. It is concluded that per lipid molecule, 11 of 14 water molecules can be replaced by three trehalose molecules. About four are displaced by changes in the water activity of the bulk solution, and seven by specific interactions with the phospholipids. In this last case, at least two of them are linked to the carbonyls, and this appears to be the cause of the decrease in the dipole potential of the membrane. In contrast, four sucrose molecules displace only three water molecules per lipid, with no effect on the dipole potential or the carbonyl groups.     

Dominant solvation effects from the primary shell of hydration: Approximation for molecular dynamics simulations

A simple approximation is developed to account for the dominant effects of solvation in molecular dynamics simulations of biopolymers. A small number of water molecules are included explicitly in the primary hydration shell around the biopolymer. A nonspherical confining potential responding dynamically to the conformational changes of the biopolymer is applied to prevent evaporation and to approximate the conditions of constant pressure of a bulk solution. Simulations of a spherical system of 25 water molecules are lined to adjust the empirical restraining potential to yield a uniform density distribution close to that in the bulk liquid. The primary hydration shell approach is tested with molecular dynamics simulations of simple hydrated peptides. The conformational equilibrium of alanine dipeptide and alanine tripeptide is examined using umbrella sampling calculations. The relative free energies of the C_7ax (? = 60, ψ = ?80) and α_L (? = 60, ψ = 60) conformations of the alanine dipeptide and the opened and closed conformations of a reversed β-turn modeled with the alanine tripeptide were calculated. The results indicate that the primary hydration shell can reproduce the influence of solvent on small peptides that was observed in simulations involving a much larger number of water molecules. © 1995 John Wiley & Sons, Inc.     

Distribution of halothane in a dipalmitoylphosphatidylcholine bilayer from molecular dynamics calculations

We report a 2-ns constant pressure molecular dynamics simulation of halothane, at a mol fraction of 50%, in the hydrated liquid crystal bilayer phase of dipalmitoylphosphatidylcholine. Halothane molecules are found to preferentially segregate to the upper part of the lipid acyl chains, with a maximum probability near the C(5) methylene groups. However, a finite probability is also observed along the tail region and across the methyl trough. Over 95% of the halothane molecules are located below the lipid carbonyl carbons, in agreement with photolabeling experiments. Halothane induces lateral expansion and a concomitant contraction in the bilayer thickness. A decrease in the acyl chain segment order parameters, S(CD), for the tail portion, and a slight increase for the upper portion compared to neat bilayers, are in agreement with several NMR studies on related systems. The decrease in S(CD) is attributed to a larger accessible volume per lipid in the tail region. Significant changes in the electric properties of the lipid bilayer result from the structural changes, which include a shift and broadening of the choline headgroup dipole (P-N) orientation distribution. Our findings reconcile apparent controversial conclusions from experiments on diverse lipid systems.     

Molecular origin of the internal dipole potential in lipid bilayers: calculation of the electrostatic potential.

The finite difference linearized Poisson-Boltzmann equation was solved for a segment of bilayer for two lipids (phosphatidylcholine dihydrate and phosphatidylethanolamine-acetic acid) in order to obtain the transbilayer electrostatic potential. Atomic coordinates derived from the crystal structures of these lipids were used, and partial changes were assigned to all atoms in the polar parts of the molecules. These calculations confirmed that a dipole potential exists in the uncharged hydrophobic interior of a bilayer. The phosphocholine and phosphoethanolamine groups make negative contributions to the internal potential, and the glycerol acyl esters make positive contributions, but the sum of these terms is negative. The water of hydration in phosphatidylcholine, and the acetic acid which is present in the phosphatidylethanolamine crystal structure, make positive contributions to the internal potential. It is concluded that the water of hydration in fully hydrated lipid bilayers is mainly responsible for the experimentally inferred positive sign of the internal potential.     

Structure-Property Relationships in Liquid Crystals: Can Modelling do Better than Empiricism?

Abstract

This paper reviews the relationships observed experimentally between the physical properties of liquid crystals and the molecular structures of the constituent molecules, and reports molecular mechanics calculations designed to provide a predictive and interpretative basis for structure/property relationships in liquid crystals. The calculations are of the minimum energy configurations of dimers of interacting liquid crystal molecules, and the geometry of the dimers, relative orientation of molecular dipoles and the extent of parallel correlation are related to liquid crystal properties. A large number of structural types are examined and the results discussed in terms of dipole correlation, apolar angular correlation, chiral twist sense and transition temperatures of liquid crystals.     

A Molecular Dynamics Study of DMPC Lipid Bilayers Interacting with Dimethylsulfoxide–Water Mixtures

The diffusion process of dimethylsulfoxide (DMSO) through zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer was studied by means of molecular dynamics (MD) simulations. To account for the cryoprotectant concentration difference between the inside and the outside of the cell, dual DMPC lipid bilayers which separate two aqueous reservoirs with and without DMSO were modeled. The initial configuration of the simulation model had DMSO molecules present in one of the aqueous phases (outside the cell) at two different concentrations of ~3 and ~6?mol%. MD simulations were performed on the systems for 50?ns at 323?K and 1?bar. Although the simulation time considered in the study was insufficient for the DMSO molecules to reach the other aqueous phase and equilibrium, early stages of the diffusion process indicated that DMSO molecules had a tendency to diffuse towards the other aqueous phase. The effects of DMSO on bilayer structural characteristics during the diffusion process were investigated. Simulations were analyzed to correlate the following properties of lipid bilayers in the presence of two different aqueous phases: area per lipid, lipid thickness, mass density profiles, lipid tail order parameter and water dipole orientation. Area per lipid calculated for the leaflet facing the aqueous DMSO?Cwater mixture did not show any significant difference compared to area per lipid for the DMSO-free pure DMPC bilayer. Mass density profiles revealed that DMSO molecules had a strong tendency to diffuse toward the aqueous phase with pure water. The lipid tail order parameter calculated for the sn-1 tail of the leaflet facing the aqueous DMSO?Cwater mixture showed that the ordering of lipid tails decreased compared to the leaflet exposed to pure water. However, the ordering of lipid tails in a system where a single bilayer is hydrated by an aqueous DMSO?Cwater mixture is far lower.     

Influence of vitamin C on alcohol binding to phospholipid monolayers

The simple model of the biological membrane is provided by well-controlled lipid monolayers at the air-water interface. The Maxwell displacement current technique (MDC) provides novel approach to conformation study of the membrane models. The effect of alcohols is interaction with membrane molecules, mainly with the lipid head group and consequent changes in physical-chemical properties of the membrane. The aim of study is to detect changes in structural, electrical and mechanical properties of dipalmitoyl-phosphatidylcholine (DPPC) monolayer on the subphase of methanol-water and ethanol-water mixtures before and after addition of antioxidant agent, vitamin C. Monolayers properties are investigated by a surface pressure analysis (including mechanical properties evaluation) and the Maxwell displacement current measurement, the dipole moment projection calculation. Surface pressure-area isotherms show similar behaviour of the DPPC monolayer on alcohol-water mixtures independently on presence of vitamin C. Binding/adsorption process induces change of electron density distribution across monolayer and thus the molecular dipole moment. We observe small or negligible binding of methanol molecules on oxygen bonds of DPPC. Thus the antioxidant, vitamin C, has no significant effect. For ethanol-water mixtures is observed recovery of electrical properties in presence of antioxidant agent. We suppose that vitamin C regulates DPPC-ethanol molecules interaction.     

Characterization of the sequence of interactions of the fusion domain of the simian immunodeficiency virus with membranes. Role of the membrane dipole potential.

The simian immunodeficiency virus fusion peptide constitutes a 12-residue N-terminal segment of the gp32 protein that is involved in the fusion between the viral and cellular membranes, facilitating the penetration of the virus in the host cell. Simian immunodeficiency virus fusion peptide is a hydrophobic peptide that in Me(2)SO forms aggregates that contain beta-sheet pleated structures. When added to aqueous media the peptide forms large colloidal aggregates. In the presence of lipidic membranes, however, the peptide interacts with the membranes and causes small changes of the membrane electrostatic potential as shown by fluorescein phosphatidylethanolamine fluorescence. Thioflavin T fluorescence and Fourier transformed infrared spectroscopy measurements reveal that the interaction of the peptide with the membrane bilayer results in complete disassembly of the aggregates originating from an Me(2)SO stock solution. Above a lipid/peptide ratio of about 5, the membrane disaggregation and water precipitation processes become dependent on the absolute peptide concentration rather than on the lipid/peptide ratio. A schematic mechanism is proposed, which sheds light on how peptide-peptide interactions can be favored with respect to peptide-lipid interactions at various lipid/peptide ratios. These studies are augmented by the use of the fluorescent dye 1-(3-sulfonatopropyl)-4-[beta[2-(di-n-octylamino)-6-naphthyl]vinyl ] pyridinium betaine that shows the interaction of the peptide with the membranes has a clear effect on the magnitude of the so-called dipole potential that arises from dipolar groups located on the lipid molecules and oriented water molecules at the membrane-water interface. It is shown that the variation of the membrane dipole potential affects the extent of the membrane fusion caused by the peptide and implicates the dipolar properties of membranes in their fusion.     

A pulsed NMR study of lipids,bound water and sodium ions in macroscopically-oriented lecithin/water lyotropic liquid crystal model membrane systems

The temperature and orientation dependence of pulsed NMR ‘free induction decay’ signals have been studied in detail for lipid bilayers macroscopically-oriented between glass slides. Results for the lipid molecules (¹H, ³¹P), bound water (²H₂O) and ions dissolved in the aqueous phase (²³Na) are presented. Bilayers of egg-lecithin, dimyristoyl lecithin and potassium oleate have been investigated. In the liquid crystal phase all the signals, including those from bound water and ions exhibit a |3 cos²? ? 1| dependence on orientation of the bilayer normal to the magnetic field. In the case of DML samples, some orientation dependence of both ¹H and ²H signals persists in the gel phase, indicating that the lipid molecules retain a degree of reorientational freedom about their long axes in this phase. At the gel-liquid crystal transition the ²H quadrupole spittings undergo a discontinuous change. Results are interpreted in terms of a model in which water molecules are bound to individual lipid head groups and reorient with them, while sodium ions are located in the aqueous channel between bilayers.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

A laser-T-jump study of the adsorption of dipolar molecules to planar lipid membranes

Phloretin and structurally related neutral molecules adsorb to the interface of lipid membranes and modify the electric dipole potential of the membrane/water interface. The adsorption process has been studied using a laser-T-jump relaxation technique in combination with an analysis of nonactin mediated potassium transport (see part I, Awiszus and Stark 1988).Deviations from the Langmuir isotherm were observed for most of the substance. The discrepancies were most pronounced at large surface densities, whereas good agreement was found at low concentrations in many cases.The partition coefficient in the limit of low concentrations was compared with that of octanol/water bulk phases. No correlation was found. The individual values of the two partition coefficients differed by more than three orders of magnitude. The contribution, b, of a single adsorbed molecule to the dipole potential could not be predicted from the dipole moment, _L , of the molecule measured in the bulk phase. Different values of b were found at identical values of _L .The study shows the limitations of the use of bulk phase data to predict molecular properties in lipid membranes.     

Theory of one-dimension rotational isomerization: A study of thecis-trans isomerization of HS-NS compared to that of HO-NO

A theoretical approach, in which the potential functions representing rotational isomerization processes are expressed in terms of linear combinations of local potentials, is presented. Partitioning the torsional potential allows identification of specific contributions that are at the origin of the shape of potential curves at different regions along the torsional variable. Key properties, such as barrier heights, may then be expressed parametrically in terms of properties associated to the stable conformations. Simple analytical expressions intended to explore, quantitatively and qualitatively, the main characteristics of the transition states connecting stable isomers are formulated. As a first step towards the study of complex systems, we use this procedure to analyseab initio results concerning thecis-trans isomerization reaction of two simple prototype molecules: HSNS and HONO. We determine the relative stabilities of the different isomers and molecular structures and evaluate the associated potential barriers. It is shown that the mathematical procedure used to get potential functions is quite convenient and may be applied to the study of more complex isomerization reactions. Numerical results concerning molecular structures, potential barriers, ionization potentials and dipole moments are discussed. Comparing the values for barrier heights suggests that S(O)-S(O) bonding through the mechanism of hyperconjugation may be present, to some extent, especially in HSNS.     

Integral equation theories for predicting water structure around molecules

Water plays a crucial role in the structure and function of proteins and other biological macromolecules; thus, theories of aqueous solvation for these molecules are of great importance. However, water is a complex solvent whose properties are still not completely understood. Statistical mechanical integral equation theories predict the density distribution of water molecules around a solute so that all particles are fully represented and thus potentially both molecular and macroscopic properties are included. Here we discuss how several theoretical tools we have developed have been integrated into an integral equation theory designed for globular macromolecular solutes such as proteins. Our approach predicts the three-dimensional spatial and orientational distribution of water molecules around a solute. Beginning with a three-dimensional Ornstein-Zernike equation, a separation is made between a reference part dependent only on the spatial distribution of solvent and a perturbation part dependent also on the orientational distribution of solvent. The spatial part is treated at a molecular level by a modified hypernetted chain closure whereas the orientational part is treated as a Boltzmann prefactor using a quasi-continuum theory we developed for solvation of simple ions. A potential energy function for water molecules is also needed and the sticky dipole models of water, such as our recently developed soft-sticky dipole (SSD) model, are ideal for the proposed separation. Moreover, SSD water is as good as or better than three point models typically used for simulations of biological macromolecules in structural, dielectric and dynamics properties and yet is seven times faster in Monte Carlo and four times faster in molecular dynamics simulations. Since our integral equation theory accurately predicts results from Monte Carlo simulations for solvation of a variety of test cases from a single water or ion to ice-like clusters and ion pairs, the application of this theory to biological macromolecules is promising.     

Differential effect of cholesterol and its biosynthetic precursors on membrane dipole potential

Dipole potential is the potential difference within the membrane bilayer, which originates due to the nonrandom arrangement of lipid dipoles and water molecules at the membrane interface. Cholesterol, a representative sterol in higher eukaryotic membranes, is known to increase membrane dipole potential. In this work, we explored the effects of immediate (7-DHC and desmosterol) and evolutionary (ergosterol) precursors of cholesterol on membrane dipole potential, monitored by the dual wavelength ratiometric approach utilizing the probe di-8-ANEPPS. Our results show that the effect of these precursors on membrane dipole potential is very different from that observed with cholesterol, although the structural differences among them are subtle. These results assume relevance, since accumulation of cholesterol precursors due to defective cholesterol biosynthesis has been reported to result in several inherited metabolic disorders such as the Smith-Lemli-Opitz syndrome. Interestingly, cholesterol (and its precursors) has a negligible effect on dipole potential in polyunsaturated membranes. We interpret these results in terms of noncanonical orientation of cholesterol in these membranes. Our results constitute the first report on the effect of biosynthetic and evolutionary precursors of cholesterol on dipole potential, and imply that a subtle change in sterol structure can significantly alter the dipolar field at the membrane interface.     

Dipole lattice membrane model for protein calculations

A dipole lattice model for lipid membranes and their interactions with peptides is presented. It uses the Langevin dipole method to calculate electrostatic interactions in the heterogeneous membrane environment. A series of test cases are presented, including spherical charges, dipoles, side chain analogs, and helical peptides. The model consistently produces qualitatively correct results.     

Partitioning of lipids at domain boundaries in model membranes

Line-active molecules (“linactants”) that bind to the boundary interface between different fluid lipid domains in membranes have a strong potential as regulators of the lateral heterogeneity that is important for many biological processes. Here, we use molecular dynamics simulations in combination with a coarse-grain model that retains near-atomic resolution to identify lipid species that can act as linactants in a model membrane that is segregated into two lipid domains of different fluidity. Our simulations predict that certain hybrid saturated/unsaturated chain lipids can bind to the interface and lower the line tension, whereas cone-shaped lysolipids have a less pronounced effect.