Distribution of pentachlorophenol in phospholipid bilayers: a molecular dynamics study

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Genome-Wide Analysis of Transcriptional Changes and Genes That Contribute to Fitness during Degradation of the Anthropogenic Pollutant Pentachlorophenol by Sphingobium chlorophenolicum

Pentachlorophenol (PCP) is a highly toxic pesticide that was first introduced in the 1930s. The alphaproteobacterium Sphingobium chlorophenolicum, which was isolated from PCP-contaminated sediment, has assembled a metabolic pathway capable of completely degrading PCP. This pathway produces four toxic intermediates, including a chlorinated benzoquinone that is a potent alkylating agent and three chlorinated hydroquinones that react with O2 to produce reactive oxygen species (ROS). RNA-seq analysis revealed that PCP causes a global stress response that resembles responses to proton motive force uncoupling and membrane disruption, while surprisingly, little of the response resembles the responses expected to be produced by the PCP degradation intermediates. Tn-seq was used to identify genes important for fitness in the presence of PCP. By comparing the genes that are important for fitness in wild-type S. chlorophenolicum and a non-PCP-degrading mutant, we identified genes that are important only when the PCP degradation intermediates are produced. These include genes encoding two enzymes that are likely to be involved in protection against ROS. In addition to these enzymes, the endogenous levels of other enzymes that protect cells from oxidative stress appear to mitigate the toxic effects of the chlorinated benzoquinone and hydroquinone metabolites of PCP. The combination of RNA-seq and Tn-seq results identify important mechanisms for defense against the toxicity of PCP. IMPORTANCE Phenolic compounds such as pentachlorophenol (PCP), triclosan, and 2,4-dichlorophenoxyacetic acid (2,4-D) represent a common class of anthropogenic biocides. Despite the novelty of these compounds, many can be degraded by microbes isolated from contaminated sites. However, degradation of this class of chemicals often generates toxic intermediates, which may contribute to their recalcitrance to biodegradation. We have addressed the stresses associated with degradation of PCP by Sphingobium chlorophenolicum by examining the transcriptional response after PCP exposure and identifying genes necessary for growth during both exposure to and degradation of PCP. This work identifies some of the mechanisms that protect cells from this toxic compound and facilitate its degradation. This information could be used to engineer strains capable of improved biodegradation of PCP or similar phenolic pollutants.

Molecular dynamics simulations of T-2410 and T-2429 HIV fusion inhibitors interacting with model membranes: Insight into peptide behavior, structure and dynamics

T-2410 and T-2429 are HIV fusion inhibitor peptides (FI) designed to present a higher efficiency even against HIV strains that developed resistance against other FIs. Similar peptides were shown to interact with model membranes both in the liquid disordered phase and in the liquid ordered state. Those results indicated that such interaction is important to function and could be correlated with their effectiveness. Extensive molecular dynamics simulations were carried out to investigate the interactions between both T-2410 and T-2429 with bilayers of pure 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and a mixture of POPC/cholesterol (Chol) (1:1). It was observed that both peptides interact strongly with both membrane systems, especially with the POPC/Chol systems, where these peptides show the highest number of H-bonds observed so far. T-2410 and T-2429 showed higher extent of interaction with bilayers when compared to T-20 or T-1249 in previous studies. This is most notable in POPC/Chol membranes where, although able to form H-bonds with Chol, they do so to a lesser extent than T-1249 does, the latter being the only FI peptide so far that was observed to form H-bonds with Chol. This behavior suggests that interaction of FI peptides with rigid Chol rich membranes may not be as dependent from peptide/Chol H-bond formation as previous results of T-1249 behavior led to believe. As in other similar peptides, the higher ability to interact with membranes shown by T-2410 and T2429 is probably correlated with its higher inhibitory efficiency.

Kinetic and Structural Aspects of the Permeabilization of Biological and Model Membranes by Lichenysin

The various lichenysins produced by Bacillus licheniformis are anionic surfactants with interesting properties. Here it is shown that lichenysin caused hemolysis of human erythrocytes, which varied with lichenysin concentration in a sigmoidal manner. The release of K(+) from red blood cells induced by lichenysin preceded the leakage of hemoglobin, and in addition, hemolysis could be impeded by the presence of compounds in the external medium having a size larger than that of PEG 3350, indicating a colloid-osmotic mechanism for hemolysis. Lichenysin also caused permeabilization of model phospholipid membranes, which was a slow process with an initial lag period of 10-20 s observed for all lichenysin concentrations. A high cholesterol ratio in the membrane decreased the extent of leakage as compared to that of pure POPC, whereas at lower ratios the effect of cholesterol was the opposite, enhancing the extent of leakage. POPE was found to decrease the extent of leakage at all the concentrations assayed, and inclusion of DPPC resulted in a considerable increase in CF leakage extent. From this scenario it was concluded that lipid membrane composition plays a role in the target membrane selectivity of lichenysin. Molecular dynamics simulations indicated that lichenysin is well distributed along the bilayer, and Na(+) ions can penetrate inside the bilayer through the lichenysin molecules. The presence of lichenysin in the membrane increases the permeability of the membrane to hydrophilic molecules facilitating its flux across the lipid palisade. The results presented in this work contribute to understanding the molecular mechanisms that explain the biological actions of lichenysin related to biomembranes.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Amphiphilic interactions of ionic liquids with lipid biomembranes: a molecular simulation study

Current bottlenecks in the large-scale commercial use of many ionic liquids (ILs) include their high costs, low biodegradability, and often unknown toxicities. As a proactive effort to better understand the molecular mechanisms of ionic liquid toxicities, the work herein presents a comprehensive molecular simulation study on the interactions of 1-n-alkyl-3-methylimidazolium-based ILs with a phosphatidylcholine (PC) lipid bilayer. We explore the effects of increasing alkyl chain length (n = 4, 8, and 12) in the cation and anion hydrophobicity on the interactions with the lipid bilayer. Bulk atomistic molecular dynamics (MD) simulations performed at millimolar (mM) IL concentrations show spontaneous insertion of cations into the lipid bilayer regardless of the alkyl chain length and a favorable orientational preference once a cation is inserted. Cations also exhibit the ability to "flip" inside the lipid bilayer (as is common for amphiphiles) if partially inserted with an unfavorable orientation. Moreover, structural analysis of the lipid bilayer show that cationic insertion induces roughening of the bilayer surface, which may be a precursor to bilayer disruption. To overcome the limitation in the timescale of our simulations, free energies for a single IL cation and anion insertion have been determined based on potential of mean force calculations. These results show a decrease in free energy in response to both short and long alkyl chain IL cation insertion, and likewise for a single hydrophobic anion insertion, but an increase in free energy for the insertion of a hydrophilic chloride anion. Both bulk MD simulations and free energy calculations suggest that toxicity mechanisms toward biological systems are likely caused by ILs behaving as ionic surfactants. [Yoo et al., Soft Matter, 2014].

Interaction of neurotransmitters with a phospholipid bilayer: a molecular dynamics study

We have performed a series of molecular dynamics simulations to study the interactions between the neurotransmitters (NTs) γ-aminobutyrate (GABA), glycine (GLY), acetylcholine (ACH) and glutamate (GLU) as well as the amidated/acetylated γ-aminobutyrate (GABA(neu)) and the osmolyte molecule glycerol (GOL) with a dipalmitoylphosphatidylcholine (DPPC) bilayer. In agreement with previously published experimental data, we found the lowest membrane affinity for the charged molecules and a moderate affinity for zwitterionic and polar molecules. The affinity can be ranked as follows: ACH-GLU<<GABA<GLY<<GABA(neu)<<GOL. The latter three penetrated the bilayer at most with the deepest location being close to the glycerol backbone of the phospholipids. Even at that position, these solutes were noticeably hydrated and carried ∼30-80% of the bulk water along. The mobility of hydration water at the solute is also affected by the penetration into the bilayer. Two time scales of exchanging water molecules could be determined. In the bulk phase, the hydration layer contains ∼20% slow exchanging water molecules which increases 2-3 times as the solutes entered the bilayer. Our results indicate that there is no intermediate exchange of slow moving water molecules from the solutes to the lipid atoms and vice versa. Instead, the exchange relies on the reservoir of unbounded ("free") water molecules in the interfacial bilayer region. Results from the equilibrium simulations are in good agreement with the results from umbrella sampling simulations, which were conducted for the four naturally occurring NTs. Free energy profiles for ACH and GLU show a minimum of ∼2-3 kJ/mol close to the bilayer interface, while for GABA and GLY, a minimum of respectively ∼2 kJ/mol and ∼5 kJ/mol is observed when these NTs are located in the vicinity of the lipid glycerol backbone. The most important interaction of NTs with the bilayer is the charged amino group of NTs with the lipid phosphate group.

Molecular dynamics simulation of HIV fusion inhibitor T-1249: insights on peptide-lipid interaction

T-1249 is a peptide that inhibits the fusion of HIV envelope with the target cell membrane. Recent results indicate that T-1249, as in the case of related inhibitor peptide T-20 (enfuvirtide), interacts with membranes, more extensively in the bilayer liquid disordered phase than in the liquid ordered state, which could be linked to its effectiveness. Extensive molecular dynamics simulations (100 ns) were carried out to investigate the interaction between T-1249 and bilayers of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and POPC/cholesterol (1 : 1). It was observed that T-1249 interacts to different extents with both membrane systems and that peptide interaction with the bilayer surface has a local effect on membrane structure. Formation of hydrogen bonding between certain peptide residues and several acceptor and donor groups in the bilayer molecules was observed. T-1249 showed higher extent of interaction with bilayers when compared to T-20. This is most notable in POPC/Chol membranes, owing to more peptide residues acting as H bond donors and acceptors between the peptide and the bilayer lipids, including H-bonds formed with cholesterol. This behavior is at variance with that of T-20, which forms no H bonds with cholesterol. This higher ability to interact with membranes is probably correlated with its higher inhibitory efficiency.

Structural effects and translocation of doxorubicin in a DPPC/Chol bilayer: the role of cholesterol

We use molecular dynamics simulations to characterize the influence of cholesterol (Chol) on the interaction between the anticancer drug doxorubicin (DOX) and a dipalmitoyl phosphatidylcholine/Chol lipid bilayer. We calculate the potential of mean force, which gives us an estimate of the free energy barrier for DOX translocation across the membrane. We find free energy barriers of 23.1 ± 3.1 k(B)T, 36.8 ± 5.1 k(B)T, and 54.5 ± 4.7 k(B)T for systems composed of 0%, 15%, and 30% Chol, respectively. Our predictions agree with Arrhenius activation energies from experiments using phospholipid membranes, including 20 k(B)T for 0% Chol and 37.2 k(B)T for 20% Chol. The location of the free energy barrier for translocation across the bilayer is dependent on composition. As Chol concentration increases, this barrier changes from the release of DOX into the water to flip-flop over the membrane center. The drug greatly affects local membrane structure by attracting dipalmitoyl phosphatidylcholine headgroups, curving the membrane, and allowing water penetration. Despite its hydrophobicity, DOX facilitates water transport via its polar groups.

Interaction of 7-nitrobenz-2-oxa-1,3-diazol-4-yl-labeled fatty amines with 1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine bilayers: a molecular dynamics study

A complete homologous series of fluorescent 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-labeled fatty amines of varying alkyl chain length, NBD-C(n), inserted in 1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers, was studied using atomistic molecular dynamics (MD) simulations. For all amphiphiles, the NBD fluorophore locates near the glycerol backbone/carbonyl region of POPC and establishes stable hydrogen bonding with POPC ester oxygen atoms. Small differences observed in the transverse location of the fluorophore correlate with other calculated parameters and with small discrepancies recently measured in the photophysical properties of the molecules. The longer-chained NBD-C(n) amphiphiles show significant mass density near the bilayer midplane, and the chains of these derivatives interdigitate to some extent the opposite bilayer leaflet. This phenomenon leads to a slower lateral diffusion for the longer-chained derivatives (n > 12). Effects of these amphiphiles on the structure and dynamics of the host lipid were found to be relatively mild, in comparison with acyl-chain-labeled NBD probes. The molecular details obtained by this work allow the rationalization of the nonmonotonic behavior, recently obtained experimentally, for the photophysical parameters of the amphiphiles and the kinetic and thermodynamic parameters for their interaction with the POPC membranes.

NMR-NOE and MD simulation study on phospholipid membranes: dependence on membrane diameter and multiple time scale dynamics

Motional correlation times between the hydrophilic and hydrophobic terminal groups in lipid membranes are studied over a wide range of curvatures using the solution-state (1)H NMR-nuclear Overhauser effect (NOE) and molecular dynamics (MD) simulation. To enable (1)H NMR-NOE measurements for large vesicles, the transient NOE method is combined with the spin-echo method, and is successfully applied to a micelle of 1-palmitoyl-lysophosphatidylcholine (PaLPC) with diameter of 5 nm and to vesicles of dipalmitoylphosphatidylcholine (DPPC) with diameters ranging from 30 to 800 nm. It is found that the NOE intensity increases with the diameter up to ∼100 nm, and the model membrane is considered planar on the molecular level beyond ∼100 nm. While the NOE between the hydrophilic terminal and hydrophobic terminal methyl groups is absent for the micelle, its intensity is comparable to that for the neighboring group for vesicles with larger diameters. The origin of NOE signals between distant sites is analyzed by MD simulations of PaLPC micelles and DPPC planar bilayers. The slow relaxation is shown to yield an observable NOE signal even for the hydrophilic and hydrophobic terminal sites. Since the information on distance and dynamics cannot be separated in the experimental NOE alone, the correlation time in large vesicles is determined by combining the experimental NOE intensity and MD-based distance distribution. For large vesicles, the correlation time is found to vary by 2 orders of magnitude over the proton sites. This study shows that NOE provides dynamic information on large vesicles when combined with MD, which provides structural information.

Simulation of nanoparticle permeation through a lipid membrane

A metric of nanoparticle toxicity is the passive permeability rate through cellular membranes. To assess the influence of nanoparticle morphology on this process, the permeability of buckyball-sized molecules through a representative lipid bilayer was investigated by molecular-dynamics simulation. When C(60) was compared with a prototypical opened C(60) molecule and a representative combustion-generated particle, C(68)H(29), the calculated free-energy profiles along the permeation coordinate revealed a sizable variation in form and depth. The orientation of the anisotropic molecules was determined by monitoring the principal axis corresponding to the largest moment of inertia, and free rotation was shown to be hindered in the bilayer interior. Diffusion constant values of the permeant molecules were calculated from a statistical average of seven to 10 trajectories at five locations along the permeation coordinate. A relatively minor variation of the values was observed in the bilayer interior; however, local resistance values spanned up to 24 orders of magnitude from the water layer to the bilayer center, due primarily to its exponential dependence on free energy. The permeability coefficient values calculated for the three similarly sized but structurally distinct nanoparticles showed a significant variance. The use of C(60) to represent similarly sized carbonaceous nanoparticles for assessments of toxicity is questioned.

Bridging the gap between molecular descriptors and mechanism: cases studies by molecular dynamics simulations

In recent years, both classification models and quantitative structure-activity relationships (QSARs) have been developed to discriminate the acute toxicity of polar narcotics and uncouplers. One of fundamental issues is how to select and interpret the molecular descriptors used in both methods. In this work, we first employed support vector machine on a dataset containing 155 polar narcotics and 19 uncouplers to filter the predictive hydrophobic and hydrogen bonding descriptors. Molecular dynamics simulations were then conducted to investigate the behavior of salicylate and pentachlorophenol molecules in the context of a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. The results demonstrated that their equilibrium properties in the lipid bilayer were closely associated with hydrophobic and hydrogen bonding descriptors. The preferable occupations of these molecules in the lipid bilayer were discussed in terms of their modes of toxic action. The observations from molecular dynamics simulations facilitated to elucidate the mechanism of polar narcotics and uncouplers.

Distribution and dynamics of adamantanes in a lipid bilayer

The adamantanes are a class of compounds that have found use in the treatment of influenza A and Parkinson's disease, among others. The mode of action for influenza A is based on the adamantanes' interaction with the transmembrane M2 channel, whereas the treatment of Parkinson's disease is thought to relate to a channel block of N-methyl-D-aspartate receptors. An understanding of how these compounds interact with the lipid bilayer is thus of great interest. We used molecular-dynamics simulations to calculate the potential of mean force of adamantanes in a lipid bilayer. Our results demonstrate a preference for the interfacial region of the lipid bilayer for both protonated and deprotonated species, with the protonated species proving significantly more favorable. However, the protonated species have a large free-energy barrier in the center of the membrane. In contrast, there is no barrier (compared with aqueous solution) at the center of the bilayer for deprotonated species, suggesting that the permeant species is indeed the neutral form, as commonly assumed. We discuss the results with respect to proposed mechanisms of action and implications for drug-delivery in general.

COSMOmic: a mechanistic approach to the calculation of membrane-water partition coefficients and internal distributions within membranes and micelles

A new approach for the modeling of molecules in micellar systems and especially in biomembranes, COSMOmic, is presented, and its performance is validated on the example of the partitioning of molecules between water and biological membranes. Starting from quantum chemical calculations of the surfactant, solvent, and solute molecules, and being based on the COSMO-RS method for fluid-phase thermodynamic properties, COSMOmic is essentially free of additional adjustable parameters. The inclusion of an elastic energy correction into the COSMOmic model did not turn out to yield any significant improvement. The novel COSMOmic method allows for the efficient prediction of the distribution of molecules in micellar systems.

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.

Defects in ordered aggregates of cardiolipin visualized by atomic force microscopy

The formation and the nature of defects in ordered aggregates of cardiolipin (tetra acyl diphosphatidylglycerol) supported on solid substrates have been investigated by atomic force microscopy (AFM). The experiments were performed on two model systems, i.e. three-dimensional liquid crystals dispersed in water and partially de-hydrated on a hydrophilic surface, and two-dimensional films of molecules self-assembled onto an isotropic hydrophobic surface. Defects were induced both by varying the preparation temperature and by treatment with specific chemicals known to modify the order parameters in natural and artificial membranes, specifically: 2,4-dinitro-phenol (DNP) and pentachloro-phenol (PCP). The effect of lipid oxidation on the nanocrystalline order was also investigated. The images obtained by AFM allow to characterize the type of defects and their local density at nanoscale level. They also provide additional information to differentiate the specific role of acyl chains and polar heads in the process of lipid self-organization.

Liposomal drug transport: a molecular perspective from molecular dynamics simulations in lipid bilayers

Computational methods to predict drug permeability across biomembranes prior to synthesis are increasingly desirable to minimize the investment in drug design and development. Significant progress in molecular dynamics (MD) simulation methodologies applied to lipid bilayer membranes, for example, is making it possible to move beyond characterization of the membranes themselves to explore various thermodynamic and kinetic processes governing membrane binding and transport. Such methods are also likely to be directly applicable to the design and optimization of liposomal delivery systems. MD simulations are particularly valuable in addressing issues that are difficult to explore in laboratory experiments due to the heterogeneity of lipid bilayer membranes at the molecular level. Insights emerging from MD simulations are contributing to an understanding of which regions within bilayers are most and least favored by solutes at equilibrium as the solute structure is varied, local diffusivities of permeants, and the origin of the amplified selectivity to permeant size imposed by lipid bilayer membranes, particularly as changes in composition increase acyl chain ordering.

Permeability of psoralen derivatives in lipid membranes

Molecular dynamics simulations have been performed to explore the distribution and translocation of a set of furocoumarins (psoralen derivatives) inside saturated and partially unsaturated lipid membranes. Within the simulations, strong accumulation of the photodynamic drugs is observed near the polar headgroup region, although the populations also extend out into the membrane/water interface as well as to the membrane center. The computed transverse (Dz) diffusion coefficients are in the range 0.01-0.03 x 10(-5) cm2 s(-1)-significantly slower than those reported for small molecules like water, ethane, and ammonia-and are related to the low mobility inside the polar headgroup region. Trimethylpsoralen (TMP) has a very low free energy barrier to transversion, only approximately 10 kJ/mol, whereas 5- and 8-methoxy psoralens (5-MOP, 8-MOP) have the largest barriers of the compounds studied-between 25 and 40 kJ/mol. Upper bounds to the permeation coefficients, obtained by integrating the resistance profiles across the bilayers, range from 5.2 x 10(-8) cm s(-1) for TMP to 4.1 x 10(-12) cm s(-1) for 5-MOP. The current simulations explain the high level of furocoumarin-lipid membrane complexes found in experimental studies of albino Wistar rats exposed to topical application of 8-MOP, and points to the possibility of membrane photodamage as a viable mechanism in psoralen ultraviolet-A treatment.

Computer simulation of partitioning of ten pentapeptides Ace-WLXLL at the cyclohexane/water and phospholipid/water interfaces

BACKGROUND

Peptide-membrane interactions play a key role in the binding, partitioning and folding of membrane proteins, the activity of antimicrobial and fusion peptides, and a number of other processes. To gain a better understanding of the thermodynamics of such interactions, White and Wimley created an interfacial hydrophobicity scale based of the transfer free energy from water to octanol or lipid bilayers of a series of synthetic peptapeptides (Ace-WLXLL, with X being any of the twenty natural amino acids) (White and Wimley (1996) Nat. Struct. Biol. 3, 842-848). In this study, we performed molecular dynamics simulations of a representative set of ten of these peptides (X = D, K, R, N, A, T, S, I, F and W) in two membrane mimetic interfaces: water-cyclohexane (10 ns) and a fully solvated dioleoylphosphatidylcholine (DOPC) bilayer (50 ns) using both constant pressure and constant area ensembles. We focus on partitioning of the ten peptides at the cyclohexane/water and lipid/water interfaces.

RESULTS

The peptides rapidly equilibrate (< 2 ns) and partition at the cyclohexane/water interface. The X3 guest residue assumes average orientations that depend on the nature of the side chain. At the DOPC/water interface, dynamics is much slower and convergence is difficult to achieve on a 50 ns timescale. Nonetheless, all peptides partition to the lipid/water interface with distributions with widths of 1-2 nm. The peptides assume a broad range of side chain and backbone orientations and have only a small effect on the area of the unit cell. On average, hydrophobic guest residues partition deeper into the hydrophobic core than hydrophilic residues. In some cases the peptides penetrate sufficiently deep to somewhat affect the distribution of the C=C double bond in DOPC. The relative distribution of the X3 guest residue compared to W1 and L5 is similar in the water/cyclohexane and water/lipid simulations. Snapshots show mostly extended backbone conformations in both environments. There is little difference between simulations at a constant area of 0.66 nm2 and simulations at constant pressure that approximately yield the same average area of 0.66 nm2.

CONCLUSION

These peptides were designed to assume extended conformations, which is confirmed by the simulations. The distribution of the X3 side chain depends on its nature, and can be determined from molecular dynamics simulations. The time scale of peptide motion at a phospholipids-water interface is too long to directly calculate the experimentally measured hydrophobicity scale to test and improve the simulation parameters. This should be possible at the water/cyclohexane interface and likely will become feasible in the future for the phospholipids/water case.

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.

Behaviour of small solutes and large drugs in a lipid bilayer from computer simulations

To reach their biological target, drugs have to cross cell membranes, and understanding passive membrane permeation is therefore crucial for rational drug design. Molecular dynamics simulations offer a powerful way of studying permeation at the single molecule level. Starting from a computer model proven to be able to reproduce the physical properties of a biological membrane, the behaviour of small solutes and large drugs in a lipid bilayer has been studied. Analysis of dihedral angles shows that a few nano seconds are sufficient for the simulations to converge towards common values for those angles, even if the starting structures belong to different conformations. Results clearly show that, despite their difference in size, small solutes and large drugs tend to lie parallel to the bilayer normal and that, when moving from water solution into biomembranes, permeants lose degrees of freedom. This explains the experimental observation that partitioning and permeation are highly affected by entropic effects and are size-dependent. Tilted orientations, however, occur when they make possible the formation of hydrogen bonds. This helps to understand the reason why hydrogen bonding possibilities are an important parameter in cruder approaches which predict drug absorption after administration. Interestingly, hydration is found to occur even in the membrane core, which is usually considered an almost hydrophobic region. Simulations suggest the possibility for highly polar compounds like acetic acid to cross biological membranes while hydrated. These simulations prove useful for drug design in rationalising experimental observations and predicting solute behaviour in biomembranes.

Partitioning of Tetrachlorophenol into Lipid Bilayers and Sarcoplasmic Reticulum: Effect of Length of Acyl Chains, Carbonyl Group of Lipids and Biomembrane Structure

We report results of a partitioning study of 2,3,4,6-tetrachlorophenol (TeCP). In the study we explored (1) the effect of the length of acyl chains of lipids (C16:1 - C24:1) and alkanes (C6-C16), (2) the role of the carbonyl group of lipids, and (3) the effect of molecular structure of the sarcoplasmic reticulum membrane on TeCP partitioning. Mole fraction partition coefficients have been measured using equilibrium dialysis for un-ionized (HA), and ionized (A) species, Kp(x) (HA), Kp(x) (A). Their values are concentration-dependent. Partition coefficients were analyzed in terms of a model that accounts for saturation of membrane associated with the finite area of partition site, and electrostatic interactions of (A-) species with charged membrane. Limiting values of partition coefficients, corresponding to infinite dilution of solute, Kp(x0) (HA), Kp(x0) (A) were obtained. Kp(x0) (HA) and Kp(x0) (A ) measure the strength of solute-membrane interactions. Studies were done with single-layered vesicles of lipids with variable chain length: 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (C16:1), 1,2-dioleoyl-sn-glycero-3-phosphocholine (C18:1), 1,2-dierucoyl-sn-glycero-3-phosphocholine (C22:1), and 1 ,2-dinervonoyl-sn-glycero-3-phosphocholine (C24:1), and egg-PC. Kp(x0) for transfer of TeCP from water into lipid membranes was found to be independent of the length of acyl chains, whereas Kp(x0) for transfer from water into alkanes increased with the length of alkane. The effect of the carbonyl CO group of lipids on partitioning was measured using 1,2-di-o-octadecenyl-sn-glycero-3-phosphocholine (CO absent) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (CO present) liposomes. Carbonyl groups, known to change dipolar potential, had no effect on partitioning. Partition coefficients of un-ionized and ionized forms of TeCP were invariant to the presence of proteins and other membrane components of sarcoplasmic reticulum (SR) membrane.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

2H-NMR study and molecular dynamics simulation of the location, alignment, and mobility of pyrene in POPC bilayers

The alignment of pyrene in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer was investigated using two different approaches, namely solid-state (2)H-NMR spectroscopy and molecular dynamics (MD) simulations. Quadrupolar splittings from (2)H-NMR spectra of deuterated pyrene-d(10) in an oriented lipid bilayer give information about the orientation of C-D bonds with respect to the membrane normal. From MD simulations, geometric information is accessible via trajectories. By defining molecular and bond order parameters, the data from MD trajectories and NMR spectra can be compared straightforwardly. To ensure that the results from both methods are comparable, parameters of the experimental and the simulation setup were chosen to be as similar as possible. From simulations, we saw that pyrene prefers a position inside the lipid membrane near the headgroups and has no tendency to diffuse from one monolayer of the membrane to the other. The results from simulation and NMR show that the normal of the molecular plane is aligned nearly perpendicular to the bilayer normal. The long axis of pyrene lies preferentially parallel to the bilayer normal within a range of +/-30 degrees . The results from the two different methods are remarkably consistent. The good agreement can be explained by the fact that the different kind of motions of a pyrene molecule are already averaged within a few nanoseconds, which is the timescale covered by the MD simulation.

Molecular dynamics simulation of a palmitoyl-oleoyl phosphatidylserine bilayer with Na+ counterions and NaCl

Two 40 ns molecular dynamics simulations of a palmitoyl-oleoyl phosphatidylserine (POPS) lipid bilayer in the liquid crystalline phase with Na(+) counterions and NaCl were carried out to investigate the structure of the negatively charged lipid bilayer and the effect of salt (NaCl) on the lipid bilayer structure. Na(+) ions were found to penetrate deep into the ester region of the water/lipid interface of the bilayer. Interaction of the Na(+) ions with the lipid bilayer is accompanied by a loss of water molecules around the ion and a simultaneous increase in the number of ester carbonyl oxygen atoms binding the ion, which define an octahedral and square pyramidal geometry. The amine group of the lipid molecule is involved in the formation of inter- and intramolecular hydrogen bonds with the carboxylate and the phosphodiester groups of the lipid molecule. The area per lipid of the POPS bilayer is unaffected by the presence of 0.15M NaCl. There is a small increase in the order parameter of carbon atoms in the beginning of the alkyl chain in the presence of NaCl. This is due to a greater number of Na(+) ions being coordinated by the ester carbonyl oxygen atoms in the water/lipid interface region of the POPS bilayer.