POPC Bilayers Supported on Nanoporous Substrates: Specific Effects of Silica-Type Surface Hydroxylation and Charge Density

Recent advances in nanotechnology bring to the forefront a new class of extrinsic constraints for remodeling lipid bilayers. In this next-generation technology, membranes are supported over nanoporous substrates. The nanometer-sized pores in the substrate are too small for bilayers to follow the substrate topology; consequently, the bilayers hang over the pores. Experiments demonstrate that nanoporous substrates remodel lipid bilayers differently from continuous substrates. The underlying molecular mechanisms, however, remain largely undetermined. Here we use molecular dynamics (MD) simulations to probe the effects of silica-type hydroxylation and charge densities on adsorbed palmitoyl-oleoylphosphatidylcholine (POPC) bilayers. We find that a 50% porous substrate decorated with a surface density of 4.6 hydroxyls/nm(2) adsorbs a POPC bilayer at a distance of 4.5 Å, a result consistent with neutron reflectivity experiments conducted on topologically similar silica constructs under highly acidic conditions. Although such an adsorption distance suggests that the interaction between the bilayer and the substrate will be buffered by water molecules, we find that the substrate does interact directly with the bilayer. The substrate modifies several properties of the bilayer-it dampens transverse lipid fluctuations, reduces lipid diffusion rates, and modifies transverse charge densities significantly. Additionally, it affects lipid properties differently in the two leaflets. Compared to substrates functionalized with sparser surface hydroxylation densities, this substrate adheres to bilayers at smaller distances and also remodels POPC more extensively, suggesting a direct correspondence between substrate hydrophilicity and membrane properties. A partial deprotonation of surface hydroxyls, as expected of a silica substrate under mildly acidic conditions, however, produces an inverse effect: it increases the substrate-bilayer distance, which we attribute to the formation of an electric double layer over the negatively charged substrate, and restores, at least partially, leaflet asymmetry and headgroup orientations. Overall, this study highlights the intrinsic complexity of lipid-substrate interactions and suggests the prospect of making two surface attributes-dipole densities and charge densities-work antagonistically toward remodeling lipid bilayer properties.

Effect of intrinsic and extrinsic factors on the simulated D-band length of type I collagen

A signature feature of collagen is its axial periodicity visible in TEM as alternating dark and light bands. In mature, type I collagen, this repeating unit, D, is 67 nm long. This periodicity reflects an underlying packing of constituent triple-helix polypeptide monomers wherein the dark bands represent gaps between axially adjacent monomers. This organization is visible distinctly in the microfibrillar model of collagen obtained from fiber diffraction. However, to date, no atomistic simulations of this diffraction model under zero-stress conditions have reported a preservation of this structural feature. Such a demonstration is important as it provides the baseline to infer response functions of physiological stimuli. In contrast, simulations predict a considerable shrinkage of the D-band (11-19%). Here we evaluate systemically the effect of several factors on D-band shrinkage. Using force fields employed in previous studies we find that irrespective of the temperature/pressure coupling algorithms, assumed salt concentration or hydration level, and whether or not the monomers are cross-linked, the D-band shrinks considerably. This shrinkage is associated with the bending and widening of individual monomers, but employing a force field whose backbone dihedral energy landscape matches more closely with our computed CCSD(T) values produces a small D-band shrinkage of < 3%. Since this force field also performs better against other experimental data, it appears that the large shrinkage observed in earlier simulations is a force-field artifact. The residual shrinkage could be due to the absence of certain atomic-level details, such as glycosylation sites, for which we do not yet have suitable data.

Nonintercalating nanosubstrates create asymmetry between bilayer leaflets

The physical properties of lipid bilayers can be remodeled by a variety of environmental factors. Here we investigate using molecular dynamics simulations the specific effects of nanoscopic substrates or external contact points on lipid membranes. We expose palmitoyl-oleoyl phosphatidylcholine bilayers unilaterally and separately to various model nanosized substrates differing in surface hydroxyl densities. We find that a surface hydroxyl density as low as 10% is sufficient to keep the bilayer juxtaposed to the substrate. The bilayer interacts with the substrate indirectly through multiple layers of water molecules; however, despite such buffered interaction, the bilayers exhibit certain properties different from unsupported bilayers. The substrates modify transverse lipid fluctuations, charge density profiles, and lipid diffusion rates, although differently in the two leaflets, which creates an asymmetry between bilayer leaflets. Other properties that include lipid cross-sectional areas, component volumes, and order parameters are minimally affected. The extent of asymmetry that we observe between bilayer leaflets is well beyond what has been reported for bilayers adsorbed on infinite solid supports. This is perhaps because the bilayers are much closer to our nanosized finite supports than to infinite solid supports, resulting in a stronger support-bilayer electrostatic coupling. The exposure of membranes to nanoscopic contact points, therefore, cannot be considered as a simple linear interpolation between unsupported membranes and membranes supported on infinite supports. In the biological context, this suggests that the exposure of membranes to nonintercalating proteins, such as those belonging to the cytoskeleton, should not always be considered as passive nonconsequential interactions.

Molecular dynamic simulation study of cholesterol and conjugated double bonds in lipid bilayers

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is believed to have anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be associated with anti-obesity effects. In this paper we extend earlier molecular dynamics (MD) simulations of pure CLA-phosphatidylcholine bilayers to investigate the comparative effects of cholesterol on bilayers composed of the two respective isomers. Simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed in which, for each isomer, the simulated bilayers contained 10% and 30% cholesterol (Chol). From MD trajectories we calculate and compare structural properties of the bilayers, including areas per molecule, thickness of bilayers, tilt angle of cholesterols, order parameter profiles, and one and two-dimensional radial distribution function (RDF), as functions of Chol concentration. While the structural effect of cholesterol is approximately the same for both isomers, we find differences at an atomistic level in order parameter profiles and in two-dimensional radial distribution functions.

Self-consistent mean-field model for palmitoyloleoylphosphatidylcholine-palmitoyl sphingomyelin-cholesterol lipid bilayers

The connection between membrane inhomogeneity and the structural basis of lipid rafts has sparked interest in the lateral organization of model lipid bilayers of two and three components. In an effort to investigate anisotropic lipid distribution in mixed bilayers, a self-consistent mean-field theoretical model is applied to palmitoyloleoylphosphatidylcholine (POPC)--palmitoyl sphingomyelin (PSM)--cholesterol mixtures. The compositional dependence of lateral organization in these mixtures is mapped onto a ternary plot. The model utilizes molecular dynamics simulations to estimate interaction parameters and to construct chain conformation libraries. We find that at some concentration ratios the bilayers separate spatially into regions of higher and lower chain order coinciding with areas enriched with PSM and POPC, respectively. To examine the effect of the asymmetric chain structure of POPC on bilayer lateral inhomogeneity, we consider POPC-lipid interactions with and without angular dependence. Results are compared with experimental data and with results from a similar model for mixtures of dioleoylphosphatidylcholine, steroyl sphingomyelin, and cholesterol.

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.

A Coarse-Grained Model Based on Morse Potential for Water and n-Alkanes

In order to extend the time and distance scales of molecular dynamics simulations, it is essential to create accurate coarse-grained force fields, in which each particle contains several atoms. Coarse-grained force fields that utilize the Lennard-Jones potential form for pairwise nonbonded interactions have been shown to suffer from serious inaccuracy, notably with respect to describing the behavior of water. In this paper, we describe a coarse-grained force field for water, in which each particle contains four water molecules, based on the Morse potential form. By molecular dynamics simulations, we show that our force field closely replicates important water properties. We also describe a Morse potential force field for alkanes and a simulation method for alkanes in which individual particles may have variable size, providing flexibility in constructing complex molecules comprised partly or solely of alkane groups. We find that, in addition to being more accurate, the Morse potential also provides the ability to take larger time steps than the Lennard-Jones, because the short distance repulsion potential profile is less steep. We suggest that the Morse potential form should be considered as an alternative for the Lennard-Jones form for coarse-grained molecular dynamics simulations.

Lateral organization of complex lipid mixtures from multiscale modeling

The organizational properties of complex lipid mixtures can give rise to functionally important structures in cell membranes. In model membranes, ternary lipid-cholesterol (CHOL) mixtures are often used as representative systems to investigate the formation and stabilization of localized structural domains ("rafts"). In this work, we describe a self-consistent mean-field model that builds on molecular dynamics simulations to incorporate multiple lipid components and to investigate the lateral organization of such mixtures. The model predictions reveal regions of bimodal order on ternary plots that are in good agreement with experiment. Specifically, we have applied the model to ternary mixtures composed of dioleoylphosphatidylcholine:18:0 sphingomyelin:CHOL. This work provides insight into the specific intermolecular interactions that drive the formation of localized domains in these mixtures. The model makes use of molecular dynamics simulations to extract interaction parameters and to provide chain configuration order parameter libraries.

An improved united atom force field for simulation of mixed lipid bilayers

We introduce a new force field (43A1-S3) for simulation of membranes by the Gromacs simulation package. Construction of the force fields is by standard methods of electronic structure computations for bond parameters and charge distribution and specific volumes and heats of vaporization for small-molecule components of the larger lipid molecules for van der Waals parameters. Some parameters from the earlier 43A1 force field are found to be correct in the context of these calculations, while others are modified. The validity of the force fields is demonstrated by correct replication of X-ray form factors and NMR order parameters over a wide range of membrane compositions in semi-isotropic NTP 1 atm simulations. 43-A1-S3 compares favorably with other force fields used in conjunction with the Gromacs simulation package with respect to the breadth of phenomena that it accurately reproduces.

Multiscale simulations of heterogeneous model membranes

This review will focus on computer modeling aimed at providing insights into the existence, structure, size, and thermodynamic stability of localized domains in membranes of heterogeneous composition. Modeling the lateral organization within a membrane is problematic due to the relatively slow lateral diffusion rate for lipid molecules so that microsecond or longer time scales are needed to fully model the formation and stability of a raft in a membrane. Although atomistic simulations currently are not able to reach this scale, they can provide data on the intermolecular forces and correlations that are involved in lateral organization. These data can be used to define coarse grained models that are capable of predictions of lateral organization in membranes. In this paper, we review modeling efforts that use interaction data from MD simulations to construct coarse grained models for heterogeneous bilayers. In this review we will discuss MD simulations done with the aim of gaining the information needed to build accurate coarse-grained models. We will then review some of the coarse-graining work, emphasizing modeling that has resulted from or has a basis in atomistic simulations.

Cholesterol packing around lipids with saturated and unsaturated chains: a simulation study

The fundamental role of cholesterol in the regulation of eukaryotic membrane structure is well-established. However the manner in which atomic level interactions between cholesterol and lipids, with varying degrees of chain unsaturation and polar groups, affect the overall structure and organization of the bilayer is only beginning to be understood. In this paper we describe a series of Molecular Dynamics simulations designed to provide new insights into lipid-cholesterol interactions as a function of chain unsaturation. We have run simulations of varying concentrations of cholesterol in dipalmitoyl phosphatidylcholine (DPPC), palmitoyl-oleyol phosphatidylcholine (POPC), and dioleyol phosphatidylcholine (DOPC) bilayers. Structural analysis of the simulations reveals both atomistic and systemic details of the interactions and are presented here. In particular, we find that the minimum partial molecular area of cholesterol occurs in POPC-Chol mixtures implying the most favorable packing. Physically, this appears to be related to the fact that the two faces of the cholesterol molecule are different from each other and that the steric cross section of cholesterol molecules drops sharply near the small chain tails.

Atomistic and coarse-grained computer simulations of raft-like lipid mixtures

Computer modeling can provide insights into the existence, structure, size, and thermodynamic stability of localized raft-like regions in membranes. However, the challenges in the construction and simulation of accurate models of heterogeneous membranes are great. The primary obstacle in modeling the lateral organization within a membrane is the relatively slow lateral diffusion rate for lipid molecules. Microsecond or longer time-scales are needed to fully model the formation and stability of a raft in a membra ne. Atomistic simulations currently are not able to reach this scale, but they do provide quantitative information on the intermolecular forces and correlations that are involved in lateral organization. In this chapter, the steps needed to carry out and analyze atomistic simulations of hydrated lipid bilayers having heterogeneous composition are outlined. It is then shown how the data from a molecular dynamics simulation can be used to construct a coarse-grained model for the heterogeneous bilayer that can predict the lateral organization and stability of rafts at up to millisecond time-scales.

Lateral organization in lipid-cholesterol mixed bilayers

Interactions between lipid and cholesterol molecules in membranes play an important role in the structural and functional properties of cell membranes. Although structural properties of lipid-cholesterol mixtures have been extensively studied, an understanding of the role of cholesterol in the lateral organization of bilayers has been elusive. In this article, we propose a simple yet powerful model, based on self-consistent mean-field theory and molecular dynamics simulations, for lipid bilayers containing cholesterol. Properties predicted by our model are shown to be in excellent agreement with experimental data. Our model predicts that cholesterol induces structural changes in the bilayer through the formation of regions of ordered lipids surrounding each cholesterol molecule. We find that the "smooth" and "rough" sides of cholesterol play crucial roles in formation and distribution of the ordered regions. Our model is predictive in that input parameters are obtained from independent atomistic molecular dynamics simulations. The model and method are general enough to describe other heterogeneous lipid bilayers, including lipid rafts.

Cholesterol surrogates: a comparison of cholesterol and 16:0 ceramide in POPC bilayers

Experimental evidence indicates that, under some circumstances, "surrogate" molecules may play the same role as cholesterol in ordering membrane lipids. The simplest molecule in this class is Ceramide. In this article, we describe atomic-level molecular dynamics simulations designed to shed light on this phenomenon. We run simulations of hydrated phosphoryl-oleoyl phosphatidylcholine (POPC) bilayers containing cholesterol, and containing ceramide, in concentrations ranging from 5% to 33%. We also perform a simulation of a pure POPC bilayer to verify the simulation force fields against experimental structural data for POPC. Our simulation data are in good agreement with experimental data for the partial molecular volumes, areas, form factors, and order parameters. These simulations suggest that ceramide and cholesterol have a very similar effect on the POPC bilayer, although ceramide is less effective in inducing order in the bilayer compared with cholesterol at the same concentrations.

Molecular-dynamics simulation of a ceramide bilayer

Ceramide is the simplest lipid in the biologically important class of glycosphingolipids. Ceramide is an important signaling molecule and a major component of the strateum corneum layer in the skin. In order to begin to understand the biophysical properties of ceramide, we have carried out a molecular-dynamics simulation of a hydrated 16:0 ceramide lipid bilayer at 368 K (5 degrees above the main phase transition). In this paper we describe the simulation and present the resulting properties of the bilayer. We compare the properties of the simulated ceramide bilayer to an earlier simulation of 18:0 sphingomyelin, and we discuss the results as they relate to experimental data for ceramide and other sphingolipids. The most significant differences arise at the lipid/water interface, where the lack of a large ceramide polar group leads to a different electron density and a different electrostatic potential but, surprisingly, not a different overall "dipole potential," when ceramide is compared to sphingomyelin.

Self-consistent mean-field model based on molecular dynamics: application to lipid-cholesterol bilayers

We have developed a dynamic self-consistent mean-field model, based on molecular-dynamics simulations, to study lipid-cholesterol bilayers. In this model the lipid bilayer is represented as a two-dimensional lattice field in the lipid chain order parameters, while cholesterol molecules are represented by hard rods. The motion of rods in the system is continuous and is not confined to lattice cells. The statistical mechanics of chain ordering is described by a mean field derived from an extension of a model due to Marcelja. The time evolution of the system is governed by stochastic equations. The ensemble of chain configurations required in partition sums, and the energies of interaction, are taken from atomistic level molecular-dynamics simulations of lipid bilayers. The model allows us to simulate systems 500 nm in lateral size for 20 micros time scales, or greater. We have applied the model to dipalmitoyl-phosphatidylcholine-cholesterol (Chol) bilayers at 50 degrees C for Chol concentrations between 2% and 33%. At low concentrations of Chol (2%-4%), the model predicts the formation of isolated clusters of Chol surrounded by relatively ordered lipid chains, randomly dispersed in the disordered bilayer. With increasing Chol composition, regions of Chol-induced order begin to overlap. Starting from about 11% Chol this ordering effect becomes system wide and regions unaffected by Chol are no longer detectable. From the analysis of properties of the model we conclude that the change in lipid chain order with increasing Chol concentration is continuous over the 20-mus scale of the simulations. We also conclude that at 50 degrees C no large-scale Chol-rich and Chol-depleted coexisting phase-separated regions form at any concentration. At no point in any of the simulations do we observe a higher degree of lateral organization, such as Chol-based superlattice structures.

Simulation of the early stages of nano-domain formation in mixed bilayers of sphingomyelin, cholesterol, and dioleylphosphatidylcholine

It is known from experimental studies that lipid bilayers composed of unsaturated phospholipids, sphingomyelin, and cholesterol contain microdomains rich in sphingomyelin and cholesterol. These domains are similar to "rafts" isolated from cell membranes, although the latter are much smaller in lateral size. Such domain formation can be a result of very specific and subtle lipid-lipid interactions. To identify and study these interactions, we have performed two molecular dynamics simulations, of 200-ns duration, of dioleylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol) systems, a 1:1:1 mixture of DOPC/SM/Chol, and a 1:1 mixture of DOPC/SM. The simulations show initial stages of the onset of spontaneous phase-separated domains in the systems. On the simulation timescale cholesterol favors a position at the interface between the ordered SM region and the disordered DOPC region in the ternary system and accelerates the process of domain formation. We find that the smooth alpha-face of Chol preferentially packs next to SM molecules. Based on a comparative analysis of interaction energies, we find that Chol molecules do not show a preference for SM or DOPC. We conclude that Chol molecules assist in the process of domain formation and the process is driven by entropic factors rather than differences in interaction energies.

Sphingomyelin-cholesterol domains in phospholipid membranes: atomistic simulation

We have carried out an atomic-level molecular dynamics simulation of a system of nanoscopic size containing a domain of 18:0 sphingomyelin and cholesterol embedded in a fully hydrated dioleylposphatidylcholine (DOPC) bilayer. To analyze the interaction between the domain and the surrounding phospholipid, we calculate order parameters and area per molecule as a function of molecule type and proximity to the domain. We propose an algorithm based on Voronoi tessellation for the calculation of the area per molecule of various constituents in this ternary mixture. The calculated areas per sphingomyelin and cholesterol are in agreement with previous simulations. The simulation reveals that the presence of the liquid-ordered domain changes the packing properties of DOPC bilayer at a distance as large as approximately 8 nm. We calculate electron density profiles and also calculate the difference in the thickness between the domain and the surrounding DOPC bilayer. The calculated difference in thickness is consistent with data obtained in atomic force microscopy experiments.

Structure of sphingomyelin bilayers: a simulation study

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50 degrees C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.

Combined Monte Carlo and molecular dynamics simulation of hydrated 18:0 sphingomyelin–cholesterol lipid bilayers

We have carried out atomic level molecular dynamics and Monte Carlo simulations of hydrated 18:0 sphingomyelin (SM)-cholesterol (CHOL) bilayers at temperatures of 20 and 50 degrees C. The simulated systems each contained 266 SM, 122 CHOL, and 11861 water molecules. Each simulation was run for 10 ns under semi-isotropic pressure boundary conditions. The particle-mesh Ewald method was used for long-range electrostatic interactions. Properties of the systems were calculated over the final 3 ns. We compare the properties of 20 and 50 degrees C bilayer systems with each other, with experimental data, and with experimental and simulated properties of pure SM bilayers and dipalmitoyl phospatidyl choline (DPPC)-CHOL bilayers. The simulations reveal an overall similarity of both systems, despite the 30 degrees C temperature difference which brackets the pure SM main phase transition. The area per molecule, lipid chain order parameter profiles, atom distributions, and electron density profiles are all very similar for the two simulated systems. Consistent with simulations from our lab and others, we find strong intramolecular hydrogen bonding in SM molecules between the phosphate ester oxygen and the hydroxyl hydrogen atoms. We also find that cholesterol hydroxyl groups tend to form hydrogen bonds primarily with SM carbonyl, methyl, and amide moieties and to a lesser extent methyl and hydroxyl oxygens.

Synaptic vesicle transport and synaptic membrane transporter sites in excitatory amino acid nerve terminals in Alzheimer disease

The apparent l-[3H]glutamate uptake rate (v') was measured in synaptic vesicles isolated from cerebral cortex synaptosomes prepared from autopsied Alzheimer and non-Alzheimer dementia cases, and age-matched controls. The initial synaptosome preparations exhibited similar densities of d-[3H]aspartate membrane binding sites (BMAX values) in the three groups. In control brain the temporal cortex d-[3H]aspartate BMAX was 132% of that in motor cortex, parallel with the l-[3H]glutamate v' values (temporal=139% of motor; NS). Unlike d-[3H]aspartate BMAX values, l-[3H]glutamate v' values were markedly and selectively lower in Alzheimer brain preparations than in controls, particularly in temporal cortex. The difference could not be attributed to differential effects of autopsy interval or age at death. Non-Alzheimer dementia cases resembled controls. The selective loss of vesicular glutamate transport is consistent with a dysfunction in the recycling of transmitter glutamate.

Cholesterol-induced modifications in lipid bilayers: a simulation study

We present analysis of new configurational bias Monte Carlo and molecular dynamics simulation data for bilayers of dipalmitoyl phosphatidyl choline and cholesterol for dipalmitoyl phosphatidyl choline:cholesterol ratios of 24:1, 47:3, 11.5:1, 8:1, 7:1, 4:1, 3:1, 2:1, and 1:1, using long molecular dynamics runs and interspersed configurational bias Monte Carlo to ensure equilibration and enhance sampling. In all cases with cholesterol concentrations above 12.5% the area per molecule of the heterogeneous membrane varied linearly with cholesterol fraction. By extrapolation to pure cholesterol, we find the cross-sectional area of cholesterol in these mixtures is approximately 22.3 A(2). From the slope of the area/molecule relationship, we also find that the phospholipid in these mixtures is in a liquid ordered state with an average cross-sectional area per lipid of 50.7 A(2), slightly above the molecular area of a pure phospholipid gel. For lower concentrations of cholesterol, the molecular area rises above the straight line, indicating the "melting" of at least some of the phospholipid into a fluid state. Analysis of the lateral distribution of cholesterol molecules in the leaflets reveals peaks in radial distributions of cholesterols at multiples of approximately 5 A. These peaks grow in size as the simulation progresses, suggesting a tendency for small subunits of one lipid plus one cholesterol, hydrogen bonded together, to act as one composite particle, and perhaps to aggregate with other composites. Our results are consistent with experimentally observed effects of cholesterol, including the condensation effect of cholesterol in phospholipid monolayers and the tendency of cholesterol-rich domains to form in cholesterol-lipid bilayers. We are continuing to analyze this tendency on longer timescales and for larger bilayer patches.

Modeling the lipid component of membranes

During the past several years, there have been a number of advances in the computational and theoretical modeling of lipid bilayer structural and dynamical properties. Molecular dynamics (MD) simulations have increased in length and time scales by about an order of magnitude. MD simulations continue to be applied to more complex systems, including mixed bilayers and bilayer self-assembly. A critical problem is bridging the gap between the still very small MD simulations and the time and length scales of experimental observations. Several new and promising techniques, which use atomic-level correlation and response functions from simulations as input to coarse-grained modeling, are being pursued.

Molecular simulation of dioleoylphosphatidylcholine lipid bilayers at differing levels of hydration

The structure and dynamics of the lipid and water components of dioleoylphosphatidylcholine bilayers at various levels of hydration were studied using molecular dynamics (MD) simulations. Equilibration of these systems proceeded by use of a hybrid MD and configurational-bias Monte Carlo technique using one atmosphere of pressure normal to the membrane and a set point for the lateral area derived from experimental Bragg spacings, combined with experimentally derived specific volumes for each of the system components. Membrane surface tensions were observed to be of the order of tens of dyn/cm. The transbilayer molecular fragment peak positions at low hydration were found to agree with experimental neutron and x-ray scattering profiles and previously published simulations. For hydration levels of 5.4, 11.4, and 16 waters/lipid, molecular fragment distributions and order parameters for the headgroup, lipid chains, and water were quantified. Spin-lattice relaxation rates and lateral self-diffusion coefficients of water agreed well with results from experimental nuclear magnetic resonance studies. Relaxation rates of the choline segments and chemical shift anisotropies for the phosphate and carbonyls were computed. Headgroup orientation, as measured by the P-N vector, showed enhanced alignment with the membrane surface at low hydration. The sign of the membrane dipole potential reversed at low hydration, with the membrane interior negative relative to the interlamellar region. Calculation of the number of water molecules in the headgroup hydration shell, as a function of hydration level, supports the hypothesis that the break point in the curve of Bragg spacing versus hydration level near 12 waters/lipid, observed experimentally by Hristova and White (1988. Biophys. J. 74:2419-2433), marks the completion of the first hydration shell.

Glutamate(NMDA) receptor NR1 subunit mRNA expression in Alzheimer's disease

We analyzed the expression profile of two NMDAR1 mRNA isoform subsets, NR1(0XX) and NR1(1XX), in discrete regions of human cerebral cortex. The subsets are characterized by the absence or presence of a 21-amino acid N-terminal cassette. Reverse transcription polymerase chain reaction for NR1 isoforms was performed on total RNA preparations from spared and susceptible regions from 10 pathologically confirmed Alzheimer's disease (AD) cases and 10 matched controls. Primers spanning the splice insert yielded two bands, 342 bp (NR1(0XX)) and 405 bp (NR1(1XX)), on agarose gel electrophoresis. The bands were visualized with ethidium and quantified by densitometry. NR1(1XX) transcript expression was calculated as a proportion of the NR1(1XX) + NR1(0XX) total. Values were significantly lower in AD cases than in controls in mid-cingulate cortex, p < 0.01, superior temporal cortex, p < 0.01 and hippocampus, p approximately 0.05. Cortical proportionate NR1(1XX) transcript expression was invariant over the range of ages and areas of controls tested, at approximately 50%. This was also true for AD motor and occipital cortex. Proportionate NR1(1XX) expression in AD cingulate and temporal cortex was lower at younger ages and increased with age: this regression was significantly different from that in the homotropic areas of controls. Variations in NR1 N-terminal cassette expression may underlie the local vulnerability to excitotoxic damage of some areas in the AD brain. Alternatively, changes in NR1 mRNA expression may arise as a consequence of the AD disease process.

Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration

We have applied a hybrid equilibration and sampling procedure for the atomic level simulation of a hydrated lipid bilayer to systems consisting of dipalmitoyl phosphatidylcholine (DPPC) and cholesterol, and palmitoyl-oleyl phosphatidylcholine (POPC) at low (approximately 6%) cholesterol concentration. The procedure is applied to bilayers of 94 molecules of DPPC, 6 molecules of cholesterol, and 3205 water molecules, and to bilayers of 120 molecules of POPC, 8 molecules of cholesterol, and 4268 water molecules, at a temperature of 325 K. After equilibration, three separate 400-ps continuous molecular dynamics runs, separated by 10,000 configurational bias Monte Carlo steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of DPPC and cholesterol, and of pure DPPC, and pure POPC. Certain properties of the bilayers are indistinguishable from cholesterol-free bilayers, including lateral diffusion and electron density. Other properties, most notably the order parameter profile, show the effect of cholesterol even at low concentrations.

Combined Monte Carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.

Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane

In this paper we report on the molecular dynamics simulation of a fluid phase hydrated dimyristoylphosphatidylcholine bilayer. The initial configuration of the lipid was the x-ray crystal structure. A distinctive feature of this simulation is that, upon heating the system, the fluid phase emerged from parameters, initial conditions, and boundary conditions determined independently of the collective properties of the fluid phase. The initial conditions did not include chain disorder characteristic of the fluid phase. The partial charges on the lipids were determined by ab initio self-consistent field calculations and required no adjustment to produce a fluid phase. The boundary conditions were constant pressure and temperature. Thus the membrane was not explicitly required to assume an area/phospholipid molecule thought to be characteristic of the fluid phase, as is the case in constant volume simulations. Normal to the membrane plane, the pressure was 1 atmosphere, corresponding to the normal laboratory situation. Parallel to the membrane plane a negative pressure of -100 atmospheres was applied, derived from the measured surface tension of a monolayer at an air-water interface. The measured features of the computed membrane are generally in close agreement with experiment. Our results confirm the concept that, for appropriately matched temperature and surface pressure, a monolayer is a close approximation to one-half of a bilayer. Our results suggest that the surface area per phospholipid molecule for fluid phosphatidylcholine bilayer membranes is smaller than has generally been assumed in computational studies at constant volume. Our results confirm that the basis of the measured dipole potential is primarily water orientations and also suggest the presence of potential barriers for the movement of positive charges across the water-headgroup interfacial region of the phospholipid.

Variant forms of neuronal glutamate transporter sites in Alzheimer's disease cerebral cortex

The displacement of Na(+)-dependent D-[3H]-aspartate binding by unlabeled D-aspartate or the inhibitors DL-threo-beta-hydroxyaspartate, L-cysteate, L-glutamate, dihydrokainate, DL-alpha-aminoadipate, alpha-methyl-DL-glutamate, and 1-aminocyclobutane-cis-1,3-dicarboxylate was used to characterize the high-affinity glutamate/aspartate uptake site in human cerebral cortex. Synaptosomal membranes were prepared from tissue obtained at autopsy from nondemented control, Alzheimer's disease (AD), and diffuse Lewy body disease (DLBD) cases. Areas that are damaged in AD (midtemporal, frontal, caudal cingulate, and hippocampal cortices) were compared with those that are spared (occipital and motor cortices). Profiles of the affinities (Ka values) of the ligands showed that areas spared from damage in AD cases differed significantly from equivalent areas in control (p < 0.001) and DLBD (p < 0.001) cases and also from areas susceptible to damage in the same AD cases (p < 0.001). Areas susceptible to damage in AD showed comparable profiles across the three case groups (p = 0.980). The glutamate/aspartate uptake site may be regionally variant in AD cases, and this may underlie local excitotoxicity. D-[3H]Aspartate binding site density was significantly lower in both dementia groups (control vs. AD, p < 0.001; control vs. DLBD, p = 0.009; but AD vs. DLBD, p = 0.528); within-group differences were not significant (control, p = 0.874; AD, p = 0.285; DLBD, p = 0.741).

Excitotoxic mechanisms in the pathogenesis of dementia

Alzheimer disease and related dementias, in common with most major neurological diseases, are characterized by localized brain damage. An abundance of senile plaques and neurofibrillary tangles in certain brain areas is pathognomic of the disease: of the two, the density of tangles may correlate more closely with disease severity ante mortem. Clinical manifestation of the disease also results from a locally severe loss of neurones. This might be caused by over-stimulation by excitant amino acid transmitters such as glutamate, which would promote cell death. Mechanisms which might give rise to the localization of Alzheimer pathogenesis include hypersensitivity to damage because a cell carries a particular sub-set of post-synaptic receptors; local variations in the efficiency of excitatory amino acid transport; and, possibly, local exacerbation of toxicity by substances such as beta-amyloid. Elucidation of such mechanisms could lead to new pharmacotherapies of dementia.

Lipid-cholesterol interactions in the P beta' phase. Application of a statistical mechanical model

We describe a statistical mechanical model for lipid-cholesterol mixtures in the P beta' (ripple) phase of lipid bilayers. The model is a simple extension of an earlier model for the ripple phase in pure lipid bilayers. The extension consists of adding a degree of freedom to allow for the occupation of underlying lattice sites by cholesterol molecules, and adding a lipid-cholesterol interaction term to the model Hamiltonian. The interaction term was constructed based on numerical calculations of lipid-cholesterol energies for several different packing juxtapositions of the two molecules. Other than the lipid-cholesterol interactions, the extended model uses the same parameter set as the earlier model, so that comparison of the properties of the extended model with experimental data serves as a test of the validity of the original model. Properties of the model were calculated using the Monte Carlo method. Results are displayed as snapshots of the ripple configurations at different cholesterol concentrations. The spacing of the ripples increases with increasing cholesterol concentration and the rate of increase compares very well with experimental data. The success of this model supports the conclusion drawn earlier that frustration arising from anisotropic packing interactions is responsible for the ripple phase in lipid bilayers. In the extended model these packing interactions are responsible for the selective partitioning of cholesterol in the regions between the ripples.

Monte Carlo studies of a model for lipid-gramicidin A bilayers

This paper presents results of Monte Carlo simulations of a full bilayer of 200 lipid chains and one gramicidin A dimer. Simulations are described for systems with lipid chains of 14, 16, and 18 carbons, respectively. Using accepted potential functions to calculate interactions between all non-hydrogen atoms a Monte Carlo configuration sampling is generated from which order parameter profiles are calculated and specific configurations are displayed. Results are compared with experimental data for lipid-gramicidin bilayers.

Lipid-cholesterol interactions. Monte Carlo simulations and theory

Results of Monte Carlo calculations of order parameter profiles of lipid chains interacting with cholesterol are presented. Cholesterol concentrations in the simulations are sufficiently large that it is possible to analyze profiles for chains which are near neighbors of two or more cholesterol molecules, chains which are neighbors to a single cholesterol, and chains which are not near any cholesterol molecules. The profiles, show that cholesterol acts to significantly decrease the ability of neighboring chains to undergo trans-gauche isomeric rotations, although these chains are not all forced into all-trans conformations. The effect is significantly greater for chains which are neighbors to more than one cholesterol. The Monte Carlo results are next used as a guide to develop a theoretical model for lipid-cholesterol mixtures. The properties of this model and the phase diagram which it predicts are described. The phase diagram is then compared with experimentally determined phase diagrams. The model calculations and the computer simulations upon which they are based yield a molecular mechanism for several of the observed phases exhibited by lipid-cholesterol mixtures. The theoretical model predicts that at low temperatures the system should exhibit solid phase immiscibility.

Monte Carlo studies of lipid chains and gramicidin A in a model membrane

The Monte Carlo method has been used to simulate the equilibrium properties of a planar array of 94 saturated lipid chains and one monomer of Gramicidin A. Chains are free to move laterally in the layer plane and to change conformation via gauche rotations and long axis rotations in a continuum. All non-hydrogen atoms on chains and on the Gramicidin A monomer interact via 6-12 potentials, and periodic boundary conditions are imposed. Calculated results consist of order parameter profiles for C-14 and C-16 chains. Profiles are calculated for chains which are neighbors to the Gramicidin A molecule and for chains which are not neighbors to the peptide. The main conclusion is that the average conformations of the chains neighboring the Gramicidin A monomer are very similar to those of the bulk chains.

Lipid chains and cholesterol in model membranes: a Monte Carlo Study

The Monte Carlo method has been employed to study the equilibrium properties of a planar array of hydrocarbon chains interacting with a cholesterol molecule. The chains are arranged to model one monolayer of a lipid bilayer and within this monolayer are allowed to move laterally and change conformations by gauche rotations. In the simulation cell there are 90 lipid chains and a single cholesterol molecule. Periodic boundary conditions are imposed upon the cell. The primary results of the calculations are order parameter profiles for the C-C bonds. These are calculated for (i) all chains, (ii) the 6 chains which are nearest neighbors to the cholesterol, and (iii) the 12 chains which are next-nearest neighbors to the cholesterol. Calculations are carried out for C-14, C-16, and C-18 chains. The results show that cholesterol strongly affects the upper portions of the chains, leaving them less able to change conformations. For C-16 and C-18 chains, the chain termini of the cholesterol neighbors are more disordered than the bulk chain termini. The magnitude of the effect depends strongly on the chain length. The results suggest that the changes in the lipid phase transition caused by cholesterol are a consequence of each cholesterol hindering the rotameric freedom of five to seven lipid chains.