Density Functional Theory-Based Conformational Analysis of a Phospholipid Molecule (Dimyristoyl Phosphatidylcholine)

Multi-stimuli responsive tetra-PPO60-PEO20 ethylene diamine block copolymer enables pH, temperature, and solvent regulation of Au nanoparticle composite plasmonic response

Dynamic plasmonic tuning of Au NP organization through pH and temperature mediated polymer architecture transformation from lamellar to bicontinuous cubic.

Altering the edge chemistry of bicelles with peptoids

Cell function is tied to the interactions that occur within and across the cell membrane. Therefore, understanding membrane-affiliated interactions is important to many biomedical applications. Advancing the body of knowledge about these interactions will lead to discoveries in biomarker detection and therapeutic targets for disease detection and treatment. Model membrane systems are an effective way to study membrane proteins for such discoveries, allowing for stable protein structure and maintaining native activity. Bicelles, disc-shaped lipid bilayers created by combining long- and short-chain phospholipids, are the model membrane system of focus in this study. Bicelles are accessible from both sides and have a wide size range, which makes them attractive for studying membrane interactions without affecting function. In this work, bicelles were functionalized with peptoids to alter the edge chemistry. Peptoids are suitable for this application because of the large diversity of available side chain chemistries that can be easily incorporated in a sequence-specific manner. The peptoid sequence consists of three functional regions to promote insertion into the edge of bicelles. The insertion sequence at the C-terminus contains two alkyl chains and two hydrophobic, chiral aromatic groups that anchor into the bicelle edge. The facially amphipathic helix contains chiral aromatic groups on one side that interact with the lipid tails and positively charged groups on the other side, which interact with the lipid head groups. Thiol groups are included at the N-terminus to allow for visualization of peptoid location in the bicelle. Bicelle morphology and size were assessed by transmission electron microscopy (TEM) and dynamic light scattering (DLS). Peptoid location in the bicelle was determined by attachment of gold nanoparticles, which confirmed preferential incorporation of the peptoid into the bicelle edge with 82% specificity. Additionally, the peptoid-functionalized bicelles are of similar size and morphology to non-functionalized bicelles. Results from this study show that peptoid-functionalized bicelles are a promising model membrane system with potential applications in biosensors or bioseparations.

Remarkable functions of sn-3 hydroxy and phosphocholine groups in 1,2-diacyl-sn-glycerolipids to induce clockwise (+)-helicity around the 1,2-diacyl moiety: Evidence from conformation analysis by 1H NMR spectroscopy

Cell-membrane glycerolipids exhibit a common structural backbone of asymmetric 1,2-diacyl-sn-glycerol bearing polar head groups in the sn-3 position. In this study, the possible effects of sn-3 head groups on the helical conformational property around the 1,2-diacyl moiety in the solution state were examined. ¹H NMR Karplus relation studies were carried out using a series of 1,2-dipalmitoyl-sn-glycerols bearing different sn-3 substituents (namely palmitoyl, benzyl, hydrogen, and phosphates). The ¹H NMR analysis indicated that the helical property around the 1,2-diacyl moiety is considerably affected by these sn-3 substituents. The sn-3 hydroxy group induced a unique helical property, which was considerably dependent on the solvents used. In CDCl₃ solution, three staggered conformers, namely gt(+), gg(−) and tg, were randomized, while in more polar solvents, the gt(+) conformer with (+)-helicity was amplified at the expense of gg(−) and tg conformers. The sn-3 phosphocholine in phosphatidylcholine exhibited a greater effect on the gt(+) conformer, which was independent of the solvents used. From the ¹H NMR analysis, the helical conformational properties around the 1,2-diacyl moiety conformed to a simple empirical rule, which permitted the proposal of a conformational diagram for 1,2-dipalmitoyl-sn-glycerols in the solution states.

Nonenzymatic Reactions above Phospholipid Surfaces of Biological Membranes: Reactivity of Phospholipids and Their Oxidation Derivatives

Phospholipids play multiple and essential roles in cells, as components of biological membranes. Although phospholipid bilayers provide the supporting matrix and surface for many enzymatic reactions, their inherent reactivity and possible catalytic role have not been highlighted. As other biomolecules, phospholipids are frequent targets of nonenzymatic modifications by reactive substances including oxidants and glycating agents which conduct to the formation of advanced lipoxidation end products (ALEs) and advanced glycation end products (AGEs). There are some theoretical studies about the mechanisms of reactions related to these processes on phosphatidylethanolamine surfaces, which hypothesize that cell membrane phospholipids surface environment could enhance some reactions through a catalyst effect. On the other hand, the phospholipid bilayers are susceptible to oxidative damage by oxidant agents as reactive oxygen species (ROS). Molecular dynamics simulations performed on phospholipid bilayers models, which include modified phospholipids by these reactions and subsequent reactions that conduct to formation of ALEs and AGEs, have revealed changes in the molecular interactions and biophysical properties of these bilayers as consequence of these reactions. Then, more studies are desirable which could correlate the biophysics of modified phospholipids with metabolism in processes such as aging and diseases such as diabetes, atherosclerosis, and Alzheimer's disease.

Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes

Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10(-5) cm(2)/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10(-8) cm(2)/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2-8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a "wagging of the tails" frequency around 30 cm(-1), which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.

Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration

We report parametrization of dipalmitoyl-phosphatidylcholine (DPPC) in the framework of the Optimized Parameters for Liquid Simulations all-atom (OPLS-AA) force field. We chose DPPC as it is one of the most studied phospholipid species and thus has plenty of experimental data necessary for model validation, and it is also one of the highly important and abundant lipid types, e.g., in lung surfactant. Overall, PCs have not been previously parametrized in the OPLS-AA force field; thus, there is a need to derive its bonding and nonbonding parameters for both the polar and nonpolar parts of the molecule. In the present study, we determined the parameters for torsion angles in the phosphatidylcholine and glycerol moieties and in the acyl chains, as well the partial atomic charges. In these calculations, we used three methods: (1) Hartree-Fock (HF), (2) second order Møller-Plesset perturbation theory (MP2), and (3) density functional theory (DFT). We also tested the effect of the polar environment by using the polarizable continuum model (PCM), and for acyl chains the van der Waals parameters were also adjusted. In effect, six parameter sets were generated and tested on a DPPC bilayer. Out of these six sets, only one was found to be able to satisfactorily reproduce experimental data for the lipid bilayer. The successful DPPC model was obtained from MP2 calculations in an implicit polar environment (PCM).

Investigation of the interface in silica-encapsulated liposomes by combining solid state NMR and first principles calculations

In the context of nanomedicine, liposils (liposomes and silica) have a strong potential for drug storage and release schemes: such materials combine the intrinsic properties of liposome (encapsulation) and silica (increased rigidity, protective coating, pH degradability). In this work, an original approach combining solid state NMR, molecular dynamics, first principles geometry optimization, and NMR parameters calculation allows the building of a precise representation of the organic/inorganic interface in liposils. {(1)H-(29)Si}(1)H and {(1)H-(31)P}(1)H Double Cross-Polarization (CP) MAS NMR experiments were implemented in order to explore the proton chemical environments around the silica and the phospholipids, respectively. Using VASP (Vienna Ab Initio Simulation Package), DFT calculations including molecular dynamics, and geometry optimization lead to the determination of energetically favorable configurations of a DPPC (dipalmitoylphosphatidylcholine) headgroup adsorbed onto a hydroxylated silica surface that corresponds to a realistic model of an amorphous silica slab. These data combined with first principles NMR parameters calculations by GIPAW (Gauge Included Projected Augmented Wave) show that the phosphate moieties are not directly interacting with silanols. The stabilization of the interface is achieved through the presence of water molecules located in-between the head groups of the phospholipids and the silica surface forming an interfacial H-bonded water layer. A detailed study of the (31)P chemical shift anisotropy (CSA) parameters allows us to interpret the local dynamics of DPPC in liposils. Finally, the VASP/solid state NMR/GIPAW combined approach can be extended to a large variety of organic-inorganic hybrid interfaces.

Headgroup mediated water insertion into the DPPC bilayer: a molecular dynamics study

Molecular dynamics simulation was performed on the 1,2-dipalmitoyl-sn-phosphocholine (DPPC) bilayer-water system using the GROMOS96 53a6 united atom force field. The transferability of force field was tested by reproducing the area per lipid within 3% accuracy from the experimental value. The simulation shows that water can penetrate much deeper inside the bilayer almost up to the starting point of the aliphatic chain. There is significant evidence from experiments that water goes deep in the DPPC bilayer, but it has not been reported from theoretical work. The mechanism of insertion of water deep inside the lipid bilayer is still not clear. In this report, for the first time, the mechanism of water insertion deep into the bilayer has been proposed. Water transport occurs by the headgroup and its first solvation shell. The trimethyl ammonium (NMe(3)) group (headgroup of DPPC) has two stable conformations at the bilayer-water interface, one outside the bilayer and another inside it. The NMe(3) group has a large clustering of water around it and takes the water molecules inside the bilayer with it during its entry into the bilayer. The water molecules penetrate into the bilayer with the help of the NMe(3) group present at the headgroup of DPPC and eventually form hydrogen bonds with carbonyl oxygen present deep inside the bilayer. Structural characteristics at the bilayer-water interface region are also reported.

Reactivity of a phospholipid monolayer model under periodic boundary conditions: a density functional theory study of the Schiff base formation between phosphatidylethanolamine and acetaldehyde

A mechanism for the formation of the Schiff base between an acetaldehyde and an amine-phospholipid monolayer model based on Dmol3/density functional theory calculations under periodic boundary conditions was constructed. This is the first time such a system has been modeled to examine its chemical reactivity at this computation level. Each unit cell contains two phospholipid molecules, one acetaldehyde molecule, and nine water molecules. One of the amine-phospholipid molecules in the cell possesses a neutral amino group that is used to model the nucleophilic attack on the carboxyl group of acetaldehyde, whereas the other has a charged amino group acting as a proton donor. The nine water molecules form a hydrogen bond network along the polar heads of the phospholipids that facilitates very fast proton conduction at the interface. Using periodic boundary conditions afforded proton transfer between different cells. The reaction takes place in two steps, namely, (1) formation of a carbinolamine and (2) its dehydration to the Schiff base. The carbinolamine is the primary reaction intermediate, and dehydration is the rate-determining step of the process, consistent with available experimental evidence for similar reactions. On the basis of the results, the cell membrane surface environment may boost phospholipid glycation via a neighboring catalyst effect.

Structure, dynamics, and energetics of lysobisphosphatidic acid (LBPA) isomers

Lysobisphosphatidic acid (LBPA), or bis(monoacylglycerol)phosphate, is a very interesting lipid, that is mainly found in late endosomes. It has several intriguing characteristics, which differ from those of other animal glycerophospholipids, that may be related to its specific functions, particularly in the metabolism of cholesterol. Its phosphodiester group is bonded at the sn-1 (sn-1') positions of the glycerols rather than at sn-3 (sn-3'); the position of the two fatty acid chains is still under debate but, increasingly, arguments favor the sn-2, sn-2' position in the native molecule, whereas isolation procedures or acidic conditions lead to the thermodynamically more stable sn-3, sn-3' structure. Because of these peculiar features, it can be expected that LBPA shape and interactions with membrane lipids and proteins are related to its structure at the molecular level. We applied quantum mechanical methods to study the structures and stabilities of the 2,2' and 3,3' LBPA isomers, using a step-by-step procedure from glycerol to precursors (in vitro syntheses) and to the final isoforms. The structures of the two positional LBPA isomers are substantially different, showing that the binding positions of the fatty acid chains on the glycerol backbone determine the shape of the LBPA molecule and thus, possibly, its functions. The 3,3' LBPA structures obtained are more stable with respect to the 2,2' form, as expected from experiment. If one argues that the in vivo synthesis starts from the present glycerol conformers and considering the most stable bis(glycero)phosphate structures, the 2,2' isoform should be the most probable isomer.

Conformational study of methylphosphocholine: a prototype for phospholipid headgroups in membranes

Phospholipid bilayers constitute the largest structural component of cell membranes, in which choline phospholipids are abundant. In this study, through a theoretical sampling on a methylphosphocholine (MePC) potential energy surface, a set of conformers was selected as a prototype for the membrane phospholipid head. We performed a detailed conformational study of such a prototype, both as an isolated moiety and in a solvated system. We used the polarizable continuum model (PCM) to account for solvation effects. We used a quantum-mechanical methodology based on density functional theory (DFT) and the 6-31G(d,p) basis set for the calculations. Through this methodology we were able to obtain a set of conformations that presented a mirror-image pattern, in good agreement with the experimental geometric values for the different phosphocholine derivatives. Potential curves for the main parameters of the dihedral space of MePC were obtained and are provided to guide future force-field parameterizations.

Structural aspects of lipid monolayers: computer simulation analyses

Extensive molecular dynamics simulations at room temperature were carried out for model films of two dissimilar lipids (DPPC and dicaprin) at the air/water interface. To study the peculiarities of the organization patterns at different average areas per molecule, surface concentrations corresponding to five almost equally spaced points along the isotherms of the two surfactants were considered. A variable of prime interest was the density distribution in a direction normal to the interface of the monolayer components: interfacial water and surfactant on one hand and the separate moieties of the lipids on the other hand. The packing pattern and cluster size dispersion were studied by means of Voronoi tessellation and radial distribution functions. Speculations regarding structural changes upon phase-state changes during film compression were made. Individual characteristics for surfactant heads and tails as well as for interfacial water were outlined and related to the available experimental data. An analysis of the diffusion coefficients revealed the limiting factors for lipid lateral and normal diffusion. Structural arguments in support of changes in monolayer dielectric properties with the area per molecule were provided.

Structural analysis of phosphatidyl choline lipids and glycerol precursors

The structures and stabilities of the glycerol (G) and glycerol 3-phosphate (G3P) isomers have been calculated in the gas phase, using an ab initio density functional theory (DFT) method. The different conformational structures are shown to be at the origin of the various phospholipid conformers, except, obviously, for the alkyl chain torsions. The G3P conformations have been examined taking into account the experimental structures of the complexation sites in the glycerol kinase (GK) and glycerol 3-phosphate acyl transferase (G3PAT) enzymes, which correspond to the first and second steps of the “de novo” phospholipid biogenesis, respectively. The conformational analysis of the glycerophosphate skeleton is shown to determine most of the structural characteristics of the phosphatidyl choline lipids, which differ by the length of their diacyl chains, i.e., dilauroyl (DL), dimyristoyl (DM), and dipalmitoyl (DP) phosphatidylcholines (PC). Higher energy conformers with kinks in the acyl chains have been found, in preparation for molecular dynamics studies of the chain melting phase transformation.

Conformational changes in stratum corneum lipids by effect of bicellar systems

Attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy was applied to study the effects of the bicelles formed by dimyristoyl-glycero-phosphocholine (DMPC) and dihexanoyl-glycero-phosphocholine (DHPC) in porcine stratum corneum (SC) in vitro. A comparison of skin samples treated and untreated with bicelles at different temperatures was carried out. The analysis of variations after treatment in the position of the symmetric CH2 stretching, CH2 scissoring, and CH2 rocking vibrations reported important information about the effect of bicelles on the skin. Bicellar systems caused a phase transition from the gel or solid state to the liquid crystalline state in the lipid conformation of SC, reflecting the major order-disorder transition from hexagonally packed to disordered chains. Grazing incidence small and wide X-ray scattering (GISAXS and GIWAXS) techniques confirmed this effect of bicelles on the SC. These results are probably related to with the permeabilizing effect previously described for the DMPC/DHPC bicelles.