Partially deuterated phospholipids as IR structure probes of conformational order in bulk and monolayer phases

Role of Hydration in Magnesium versus Calcium Ion Pairing with Carboxylate: Solution and the Aqueous Interface

The binding of group II metal cations such as Ca²⁺ and Mg²⁺ has been largely categorized as electrostatic or ionic using carboxylate symmetric and asymmetric stretching frequency assignments that have been historically used with little regard for the solvation environment of aqueous solutions. However, given the importance of these cations and their binding mechanisms related to biological function and in revealing surface enrichment factors for ocean to marine aerosol transfer, it is imperative that a deeper understanding be sought to include hydration effects. Here, infrared reflection-absorption and Raman spectra for surface and solution phase carboxylate binding information, respectively, are compared against bare (unbound) carboxylate and bidentate Zn²⁺:carboxylate spectral signatures. Spectral non-coincidence effect analysis, temperature studies, and spectral and potential of mean force calculations result in a concise interpretation of binding motifs that include the role of mediating water molecules, that is, contact and solvent-shared ion pairs. Calcium directly binds to the carboxylate group in contact ion pairs where magnesium rarely does. Moreover, we reveal the dominance of the solvent-shared ion pair of magnesium with carboxylate at the air-water interface and in solution.

The coarse-grained models of poly(ethylene oxide) and poly(propylene oxide) homopolymers and poloxamers in big multipole water (BMW) and MARTINI frameworks

Polyethylene oxide (PEO) and poly(propylene oxide) (PPO), especially their tri-block copolymers PEO-PPO-PEO (poloxamers), have a broad range of applications in biotechnology and medical science. Understanding their specific interactions with biomembranes is the key to unveil the unique features of poloxamers either as membrane-healing or membrane pore-forming agents. Based on the coarse-graining convention of the MARTINI force field and the big multipole water (BMW) model, which has a three charged site topology and can reproduce the correct dipole moment of four-water clusters, we generated coarse-grained (CG) models with analytical and numerical potentials for PEO and PPO homopolymers and poloxamers in dilute solution. The effective bonded interaction potentials between CG beads were determined from the probability distributions of bond lengths, angles and dihedrals that are determined from atomistic simulations. The nonbonded interaction parameters were fine-tuned to reproduce the conformational properties of atomistic PEO and PPO homopolymers and poloxamers via extensive CG simulations of PEO and PPO homopolymers and poloxamers in a BMW water environment. The reported CG models provide a promising framework for a comprehensive understanding of the microstructural, conformational, and dynamic properties of poloxamers and their delicate interactions with other species in an explicit water environment.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

Interaction of DNA with Cationic Lipid Mixtures-Investigation at Langmuir Lipid Monolayers

Four different binary lipid mixtures composed of a cationic lipid and the zwitterionic colipids DOPE or DPPC, which show different DNA transfer activities in cell culture models, were investigated at the soft air/water interface to identify transfection efficiency determining characteristics. Langmuir films are useful models to investigate the interaction between DNA and lipid mixtures in a two-dimensional model system by using different surface sensitive techniques, namely, epifluorescence microscopy and infrared reflection-absorption spectroscopy. Especially, the effect of adsorbed DNA on the properties of the lipid mixtures has been examined. Distinct differences between the lipid composites were found which are caused by the different colipids of the mixtures. DOPE containing lipid mixtures form fluid monolayers with a uniform distribution of the fluorescent probe in the presence and absence of DNA at physiologically relevant surface pressures. Only at high nonphysiological pressures, the lipid monolayer collapses and phase separation was observed if DNA was present in the subphase. In contrast, DPPC containing lipid mixtures show domains in the liquid condensed phase state in the presence and absence of DNA in the subphase. The adsorption of DNA at the positively charged mixed lipid monolayer induces phase separation which is expressed in the morphology and the point of appearance of these domains.

Coarse-graining poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO) block copolymers using the MARTINI force field

The MARTINI coarse-grain (CG) force field is extended for a class of triblock block copolymers known as Pluronics. Existing MARTINI bead types are used to model the non-bonded part of the potential while single chain properties for both homopolymers, poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO), are used to develop the bonded interactions. The new set of force field parameters reproduces structural and dynamical properties of high molecular weight homo- and copolymers. The CG model is moderately transferable in solvents of different polarity and concentration; however, the PEO homopolymer model presents a reduced thermodynamic transferability especially in water probably due to the lack of hydrogen bonds with the solvent. Our simulations of a monolayer of Pluronic L44 show polymer-brush-like characteristics for the PEO segments which protrude into the aqueous phase. Other membrane properties not easily accessible using experimental techniques such as its membrane thickness are also calculated.

Physicochemical studies on the interaction of serum albumin with pulmonary surfactant extract in films and bulk bilayer phase

Functionality, structure and composition of the adsorbed films of bovine lipid extract surfactant (BLES), in the absence and presence of bovine serum albumin (BSA), at the air-buffer interface was characterized through surface tension, atomic force microscopy and time of flight secondary ion mass spectrometric methods. Gel and fluid domains of BLES films were found to be altered significantly in the presence of BSA. Differential scanning calorimetric studies on BLES dispersions in presence of BSA revealed that the perturbations of the lipid bilayer structures were significant only at higher amount of BSA. FTIR studies on the BLES dispersions in buffer solution revealed that BSA could affect the lipid head-group hydrations in bilayer as well as the methylene and methyl vibration modes of fatty acyl chains of the phospholipids present in BLES. Serum albumin could perturb the film structure at pathophysiological concentration while higher amount of BSA was required in perturbing the bilayer structures. The studies suggest a connected perturbed bilayer to monolayer transition model for surfactant inactivation at the alveolar-air interface in dysfunctional surfactants.

Electrochemical properties of phospholipid monolayers at liquid-liquid interfaces

Biomembrane models built at the interface between two immiscible electrolytes (ITIES) are useful systems to study phenomena of biological relevance by means of their electrochemical processes. The unique properties of ITIES allow one either to control or measure the potential difference across the biomimetic membranes. Herein we focus on phospholipid monolayers adsorbed at liquid-liquid interfaces, and besides discussing recent developments on the subject, we describe electrochemical techniques that can be used to get insight on the interfacial processes and electrostatic properties of phospholipid membranes at the ITIES. In particular, we examine the electrochemical and physicochemical properties of (modified) phospholipid monolayers and their interaction with other biologically relevant compounds. The use of liquid-liquid electrochemistry as a powerful tool to characterize drug properties is outlined. Although this review is not a survey of all the work in the field, it provides a comprehensive referencing to current research.

Molecular dynamics simulation study of a pulmonary surfactant film interacting with a carbonaceous nanoparticle

This article reports an all-atom molecular dynamics simulation to study a model pulmonary surfactant film interacting with a carbonaceous nanoparticle. The pulmonary surfactant is modeled as a dipalmitoylphosphatidylcholine monolayer with a peptide consisting of the first 25 residues from surfactant protein B. The nanoparticle model with a chemical formula C188H53 was generated using a computational code for combustion conditions. The nanoparticle has a carbon cage structure reminiscent of the buckyballs with open ends. A series of molecular-scale structural and dynamical properties of the surfactant film in the absence and presence of nanoparticle are analyzed, including radial distribution functions, mean-square displacements of lipids and nanoparticle, chain tilt angle, and the surfactant protein B peptide helix tilt angle. The results show that the nanoparticle affects the structure and packing of the lipids and peptide in the film, and it appears that the nanoparticle and peptide repel each other. The ability of the nanoparticle to translocate the surfactant film is one of the most important predictions of this study. The potential of mean force for dragging the particle through the film provides such information. The reported potential of mean force suggests that the nanoparticle can easily penetrate the monolayer but further translocation to the water phase is energetically prohibitive. The implication is that nanoparticles can interact with the lung surfactant, as supported by recent experimental data by Bakshi et al.

Different properties of H2O- and D2O-containing phospholipid-based reverse micelles near a critical temperature

Dipalmitoylphosphatidylcholine (DPPC)/water/pyridine reverse micelles have been found to transform from a clear liquid into a glass when the DPPC-to-water volume fraction is in the 0.78-0.89 range at 28 or 26 degrees C depending on whether water is H2O or D2O. Their study by SANS, FT-IR, and 1H NMR for this composition has shown remarkable effects of the isotopic nature of water on their structural and dynamic properties. By SANS, between 38 and 43.5 degrees C, micelles appear as either flexible polymer-like cylinders or short rods depending on whether water is H2O or D2O. On the basis of this dual aspect, micelles have been visualized as branched cylinders whose quasi-spherical branching points would be prone to assemble into short rods. In addition, when water contains more than 40% of D2O, a Bragg reflection emerges at 0.12 A(-1) on SANS spectra, evidencing an organization of micelles. In addition, FT-IR spectra show that DPPC phosphate groups are D bonded only when water is D2O. Consequently, we assumed that forces prone to organize the D2O-containing micelles are D-bonded water bridges between neighboring micelles at the level of their branching points. In fact, ab initio calculations have shown that water dimers are more stable when the bridging atom is D rather than H. These water bridges could be formed due to the fact that branching points, able to slide along micelles, keep close for a longer time when water is D2O than when it is H2O. Indeed, it has been shown experimentally that the lateral diffusion of phospholipid molecules in any layer is slower in the first case. Formation of such bridges triggers a deuteron migration between micelles evidenced by the 1/T1 relaxation rate of deuterons of water in D2O-containing micelles measured at 43 degrees C by 1H NMR.

Molecular Dynamics Simulations of Neutral Chlorpromazine in Zwitterionic Phospholipid Monolayers

Molecular dynamics simulations have been performed to investigate the interactions between chlorpromazine (CPZ), a neuroleptic drug used in the treatment of psychiatric disorders, and dipalmitoylphosphatidylcholine (DPPC), a zwitterionic phospholipid, in Langmuir monolayers. The results from simulations carried out at different monolayer surface densities were able to capture important features of the CPZ-lipid interaction. We find that neutral (unprotonated) CPZ is preferentially located in the lipid tail region of the phospholipids, in little contact with the aqueous phase, and that the orientation of its rigid ring structure and tail conformation vary with lipid surface density. CPZ is found to promote ordering of the lipid tails for all surface densities because of a reduction in the effective surface area per lipid upon addition of the drug. Similar effects have been observed in previous studies of cholesterol in DPPC monolayers, in which lipid tails were seen to order around the solute. This feature, however, is quite distinct from what we observe for the most dense monolayer considered here (area per lipid of 50 A(2)), for which we find that CPZ promotes a local distortion of the lipid tails in its immediate vicinity and a concomitant ordering of lipid tails located further away from the solute. This view is further supported by the results obtained for an approximated nonlinear vibrational sum frequency generation susceptibility, which showed greater tail disorder close to CPZ.

Dynamics of phospholipid monolayers on polarised liquid–liquid interfaces

Quasi-elastic laser light scattering (QELS) is used to investigate dynamics of the polarised water/1,2-dichloroethane (DCE) interface in the presence of adsorbed DL-alpha-dipalmitoylphosphatidylcholine (DPPC) over a range of the interfacial potential differences (+/-0.3 V) and DPPC concentrations (0-20 microM). An analysis of the frequency of thermally excited capillary waves reveals some novel features in the adsorption of DPPC. The effect of the capillary wavenumber on the capillary wave frequency and the damping factor suggest that the dynamic behaviour of ITIES is consistent with the theoretical predictions for a sharp liquid-liquid interface.

Simulations of zwitterionic and anionic phospholipid monolayers

Results of atomistic molecular dynamics simulations of dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylglycerol monolayers at the air/water interface are presented. Dipalmitoylphosphatidylcholine is zwitterionic and dipalmitoylphosphatidylglycerol is anionic at physiological pH. NaCl and CaCl2 water subphases are simulated. The simulations are carried out at different surface densities, and a simulation cell geometry is chosen that greatly facilitates the investigation of phospholipid monolayer properties. Ensemble average monolayer properties calculated from simulation are in agreement with experimental measurements. The dependence of the properties of the monolayers on the surface density, the type of the headgroup, and the ionic environment are explained in terms of atomistically detailed pair distribution functions and electron density profiles, demonstrating the strength of simulations in investigating complex, multicomponent systems of biological importance.

Fourier transform infrared spectroscopy as a probe for the study of the hydration of lipid self-assemblies. I. Methodology and general phenomena

An algorithm for the study of the gradual hydration of phospholipid assemblies by means of Fourier transform infrared (FTIR) spectroscopy is presented. A complete series of diacyl phosphatidylcholines (PCs) including all possible analogues with palmitoyl and oleoyl residues, namely DPPC, DOPC, POPC, and OPPC, was investigated at room temperature. The lipid samples were prepared as cast films probably consisting of aligned multilamellar bilayers. The range of water activities studied in these films was regulated by adsorption via the gas phase corresponding to relative humidities of between 0 and 100%. Analyses of the IR-spectroscopic data have concentrated mainly on determining the amounts of water incorporated by each lipid as well as the hydration-induced response observed for some absorption bands of the different lipids. The water uptake at high relative humidity (RH) increases with the portion of unsaturated acyl chains in the molecular structure of the PCs. Isothermal phase transitions triggered lyotropically have been detected in demonstrating the occurrence of the main transition in POPC and OPPC films at room temperature. Moreover, it appears that both lamellar phases, the gel as well as the liquid-crystalline phase, are not uniform. They seem to comprise an amazingly large span of order/disorder states of the lipid chains generally depending on the degree of hydration. As exemplified by the significant variation in the onset of wavenumber shifts for the PO2- and C=O stretching-vibration modes, obtained as a function of hydration, a sequence of attachment to polar lipid binding sites by water molecules was established for DPPC.

Domain structure and molecular conformation in annexin V/1,2-dimyristoyl-sn-glycero-3-phosphate/Ca2+ aqueous monolayers: a Brewster angle microscopy/infrared reflection-absorption spectroscopy study

Annexins comprise a family of proteins that exhibit a Ca2+-dependent binding to phospholipid membranes that is possibly relevant to their in vivo function. Although substantial structural information about the ternary (protein/lipid/Ca2+) interaction in bulk phases has been derived from a variety of techniques, little is known about the temporal and spatial organization of ternary monolayer films. The effect of Ca2+ on the interactions between annexin V (AxV) and anionic DMPA monolayers was therefore investigated using three complementary approaches: surface pressure measurements, infrared reflection-absorption spectroscopy (IRRAS), and Brewster angle microscopy (BAM). In the absence of Ca2+, the injection of AxV into an aqueous subphase beneath a DMPA monolayer initially in a liquid expanded phase produced BAM images revealing domains of protein presumably surrounded by liquid-expanded lipid. The protein-rich areas expanded with time, resulting in reduction of the area available to the DMPA and, eventually, in the formation of condensed lipid domains in spatial regions separate from the protein film. There was thus no evidence for a specific binary AxV/lipid interaction. In contrast, injection of AxV/Ca2+ at a total Ca2+ concentration of 10 microM beneath a DMPA monolayer revealed no pure protein domains, but rather the slow formation of pinhead structures. This was followed by slow (>2 h) rigidification of the whole film accompanied by an increase in surface pressure, and connection of solid domains to form a structure resembling strings of pearls. These changes were characteristic of this specific ternary interaction. Acyl chain conformational order of the DMPA, as measured by nu(sym)CH2 near 2850 cm(-1), was increased in both the AxV/DMPA and AxV/DMPA/Ca2+ monolayers compared to either DMPA monolayers alone or in the presence of Ca2+. The utility of the combined structural and temporal information derived from these three complementary techniques for the study of monolayers in situ at the air/water interface is evident from this work.