Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Molecular simulation of rapid translocation of cholesterol, diacylglycerol, and ceramide in model raft and nonraft membranes

The translocation of lipids across membranes (flip-flop) is an important biological process. Slow exchange on a physiological timescale allows the creation of asymmetric distributions of lipids across cellular membranes. The location of lipids and their rate of exchange have important biological consequences, especially for lipids involved in cellular signaling. We investigated the translocation of cholesterol, ceramide, and diacylglycerol in two model bilayers using molecular dynamics simulations. We estimate half times for flip-flop for cholesterol, diacylglycerol, and ceramide of 20 μs, 30 μs, and 10 ms in a POPC bilayer, compared with approximately 30 min, 30 ms, and 30 s in a model raft bilayer (1:1:1 PSM, POPC, and cholesterol). Cholesterol has a large (54 kJ/mol) free energy of exchange between the POPC and raft bilayer, and therefore, it strongly prefers the more ordered and rigid raft bilayer over the more liquid POPC bilayer. Ceramide and diacylglycerol have relatively small free energies of exchange, suggesting nearly equal preference for both bilayers. This unexpected result may have implications for ceramide and diacylglycerol signaling and membrane localization.     

Molecular dynamics simulation of d-Benzedrine transmitting through molecular channels within D3R

Dex-Benzedrine (known as d-Benzedrine or SAT) acts in dopamine receptors of central nerve cell system. In clinic, SAT is used to treat a variety of diseases; meanwhile, it has dependence and addiction. In order to investigate the pharmacology and addiction mechanisms of SAT as a medicine, in this paper, we have studied the structure of D₃R complex protein with SAT, and based on which, using potential mean force with umbrella samplings and the simulated phospholipid bilayer membrane (or POPC bilayer membrane), the molecular dynamics simulation was performed to obtain free energy changes upon the trajectories for SAT moving along the molecular channels within D₃R. The free energy change for SAT transmitting toward the outside of cell along the functional molecular channel within D₃R is 83.5 kJ mol^?1. The change of free energy for SAT to permeate into the POPC bilayer membrane along the protective molecular channel within D₃R is 87.7 kJ mol^?1. Our previous work gave that the free energy for Levo-Benzedrine (RAT) transmitting toward the outside of cell along the functional molecular channel within D₃R is 91.4 kJ mol^?1, while it is 117.7 kJ mol^?1 for RAT to permeate into the POPC bilayer membrane along the protective molecular channel within D₃R. The values of free energy suggest that SAT relatively prefers likely to pass through the functional molecular channel within D₃R for increasing the release of dopamine molecules resulting in a variety of functional effects for SAT. The obtained results show that the pharmacology and addiction mechanisms of SAT as a drug are closely related to the molecular dynamics and mechanism for SAT transmitting along molecular channels within D₃R.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

The energetics of transmembrane helix insertion into a lipid bilayer

Free energy profiles for insertion of a hydrophobic transmembrane protein α-helix (M2 from CFTR) into a lipid bilayer have been calculated using coarse-grained molecular dynamics simulations and umbrella sampling to yield potentials of mean force along a reaction path corresponding to translation of a helix across a lipid bilayer. The calculated free energy of insertion is smaller when a bilayer with a thinner hydrophobic region is used. The free energies of insertion from the potentials of mean force are compared with those derived from a number of hydrophobicity scales and with those derived from translocon-mediated insertion. This comparison supports recent models of translocon-mediated insertion and in particular suggests that: 1), helices in an about-to-be-inserted state may be located in a hydrophobic region somewhat thinner than the core of a lipid bilayer; and/or 2), helices in a not-to-be-inserted state may experience an environment more akin (e.g., in polarity/hydrophobicity) to the bilayer/water interface than to bulk water.     

Effects of a carane derivative local anesthetic on a phospholipid bilayer studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of two hydrated palmitoyloleoylphosphatidylcholine (POPC) bilayers each containing eight carane derivative (KP-23) local anesthetic (LA) molecules in neutral (POPC-LA) or protonated (POPC-LAH) forms were carried out to investigate the effect of KP-23 and its protonation on the bilayer. 3-ns trajectories were used for analyses. A pure POPC bilayer was employed as a reference system. In both POPC-LA and POPC-LAH systems a few KP-23 molecules intercalated into the bilayer and moved near the bilayer/water interface. They were located on the hydrophobic core side of the interface in the POPC-LA bilayer, but on the water phase side in the POPC-LAH bilayer. The order of the POPC chains was higher in the POPC-LA bilayer than in the pure POPC bilayer and was lower in the POPC-LAH bilayer. Interactions between polar groups of KP-23 and POPC or water were responsible for a lower hydration of POPC headgroups in POPC bilayers containing KP-23 than in the pure POPC bilayer. KP-23 molecules were found to form aggregates both in POPC-LA and POPC-LAH bilayers. Due to higher amphiphilicity of LAH, the LAH aggregate was more micelle-like and larger than the LA one. The results demonstrate the rapid timescales of the initial processes that take place at and near the bilayer interface as well as details of the atomic level interactions between local anesthetic and the lipid matrix of a cell membrane.     

Molecular dynamics simulation of structural changes of lipid bilayers induced by shock waves: Effects of incident angles

Unsteady and nonequilibrium molecular dynamics simulations of the response of dipalmitoylphosphatidylcholine (DPPC) bilayers to the shock waves of various incident angles are presented. The action of an incident shock wave is modeled by adding a momentum in an oblique direction to water molecules adjacent to a bilayer. We thereby elucidate the effects of incident shock angles on (i) collapse and rebound of the bilayer, (ii) lateral displacement of headgroups, (iii) tilts of lipid molecules, (iv) water penetration into the hydrophobic region of the bilayer, and (v) momentum transfer across the bilayer. The number of water molecules delivered into the hydrophobic region is found to be insensitive to incident shock angles. The most important structural changes are the lateral displacement of headgroups and tilts of lipid molecules, which are observed only in the half of the bilayer directly exposed to a shock wave for all incident shock angles studied here. As a result, only the normal component of the added oblique momentum is substantially transferred across the bilayer. This also suggests that the irradiation by shock waves may induce a jet-like streaming of the cytoplasm toward the nucleus.     

Molecular dynamics simulation of structural changes of lipid bilayers induced by shock waves: Effects of incident angles

Unsteady and nonequilibrium molecular dynamics simulations of the response of dipalmitoylphosphatidylcholine (DPPC) bilayers to the shock waves of various incident angles are presented. The action of an incident shock wave is modeled by adding a momentum in an oblique direction to water molecules adjacent to a bilayer. We thereby elucidate the effects of incident shock angles on (i) collapse and rebound of the bilayer, (ii) lateral displacement of headgroups, (iii) tilts of lipid molecules, (iv) water penetration into the hydrophobic region of the bilayer, and (v) momentum transfer across the bilayer. The number of water molecules delivered into the hydrophobic region is found to be insensitive to incident shock angles. The most important structural changes are the lateral displacement of headgroups and tilts of lipid molecules, which are observed only in the half of the bilayer directly exposed to a shock wave for all incident shock angles studied here. As a result, only the normal component of the added oblique momentum is substantially transferred across the bilayer. This also suggests that the irradiation by shock waves may induce a jet-like streaming of the cytoplasm toward the nucleus.     

Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Rhodopsin is the only member of the pharmacologically important superfamily of G-protein-coupled receptors with a known structure at atomic resolution. A molecular dynamics model of rhodopsin in a POPC phospholipid bilayer was simulated for 15 ns, revealing a conformation significantly different from the recent crystal structures. The structure of the bilayer compared with a protein-free POPC control indicated hydrophobic matching with the nonpolar interface of the receptor, in agreement with deuterium NMR experiments. A new generalized molecular surface method, based on a three-dimensional Voronoi cell construction for atoms with different radii, was developed to quantify cross-sectional area profiles for the protein, lipid acyl chains and headgroups, and water. Thus, it was possible to investigate the bilayer deformation due to curvature of the individual lipid monolayers. Moreover, the generalized molecular surface derived hydrophobic interface allowed benchmarking of the hydropathy sequence analysis, an important structural genomics tool. Five water molecules diffused into internal hydration sites during the simulation, yielding a total of 12 internal waters. The cytoplasmic loops and the C-terminal tail, containing the G-protein recognition and protein sorting sequences, exhibited a high mobility, in marked contrast to the extracellular and transmembrane domains. The proposed functional coupling of the highly conserved ERY motif to the lipid-water interface via the cytoplasmic loops provides insight into lipid effects on G-protein-coupled receptor activation in terms of a flexible surface model, involving the spontaneous monolayer curvature.     

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).     

Lipid bilayer deformation and the free energy of interaction of a Kv channel gating-modifier toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.     

Binding and reorientation of melittin in a POPC bilayer: Computer simulations

We performed, using an all-atom force field, molecular dynamics computer simulations to study the binding of melittin to the POPC bilayer and its subsequent reorientation in this bilayer. The binding process involves a simultaneous folding and adsorption of the peptide to the bilayer, followed by the creation of a "U shaped" conformation. The reorientation of melittin from the parallel to the perpendicular conformation requires charged residues to cross the hydrophobic core of the bilayer. This is accomplished by a creation of defects in the bilayer that are filled out with water. The defects are caused by peptide charged residues dragging the lipid headgroup atoms along with them, as they reorient. With increased concentration of melittin water defects form stable pores; this makes it easier for the peptide N-terminus to reorient. Our results complement experimental and computational observations of the melittin/lipid bilayer interaction.     

Tuning a polar molecule for selective cytoplasmic delivery by a pH (Low) insertion peptide

Drug molecules are typically hydrophobic and small in order to traverse membranes to reach cytoplasmic targets, but we have discovered that more polar molecules can be delivered across membranes using water-soluble, moderately hydrophobic membrane peptides of the pHLIP (pH low insertion peptide) family. Delivery of polar cargo molecules could expand the chemical landscape for pharmacological agents that have useful activity but are too polar by normal drug criteria. The spontaneous insertion and folding of the pHLIP peptide across a lipid bilayer seeks a free energy minimum, and insertion is accompanied by a release of energy that can be used to translocate cell-impermeable cargo molecules. In this study, we report our first attempt to tune the hydrophobicity of a polar cargo, phallacidin, in a systematic manner. We present the design, synthesis, and characterization of three phallacidin cargoes, where the hydrophobicity of the cargo was tuned by the attachment of diamines of various lengths of hydrophobic chains. The phallacidin cargoes were conjugated to pHLIP and shown to selectively inhibit the proliferation of cancer cells in a concentration-dependent manner at low pH.     

Thermodynamic measurements of bilayer insertion of a single transmembrane helix chaperoned by fluorinated surfactants

Accurate determination of the free energy of transfer of a helical segment from an aqueous into a transmembrane (TM) conformation is essential for understanding and predicting the folding and stability of membrane proteins. Until recently, direct thermodynamically sound measurements of free energy of insertion of hydrophobic TM peptides were impossible due to peptide aggregation outside the lipid bilayer. Here, we overcome this problem by using fluorinated surfactants that are capable of preventing aggregation but, unlike detergents, do not themselves interact with the bilayer. We have applied the fluorescence correlation spectroscopy methodology to study surfactant-chaperoned insertion into preformed POPC (palmitoyloleoylphosphatidylcholine) vesicles of the two well-studied dye-labeled TM peptides of different lengths: WALP23 and WALP27. Extrapolation of the apparent free-energy values measured in the presence of surfactants to a zero surfactant concentration yielded free-energy values of -9.0±0.1 and -10.0±0.1 kcal/mol for insertion of WALP23 and WALP27, respectively. Circular dichroism measurements confirmed helical structure of peptides in lipid bilayer, in the presence of surfactants, and in aqueous mixtures of organic solvents. From a combination of thermodynamic and conformational measurements, we conclude that the partitioning of a four-residue L-A-L-A segment in the context of a continuous helical conformation from an aqueous environment into the hydrocarbon core of the membrane has a favorable free energy of 1 kcal/mol. Our measurements, combined with the predictions of two independent experimental hydrophobicity scales, indicate that the per-residue cost of transfer of the helical backbone from water to the hydrocarbon core of the lipid bilayer is unfavorable and is equal to +2.13±0.17 kcal/mol.     

Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study.

The proton transfer activity of the light-driven proton pump, bacteriorhodopsin (bR) in the photochemical cycle might imply internal water molecules. The free energy of inserting water molecules in specific sites along the bR transmembrane channel has been calculated using molecular dynamics simulations based on a microscopic model. The existence of internal hydration is related to the free energy change on transfer of a water molecule from bulk solvent into a specific binding site. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each hydrated model from molecular dynamics simulations of the creation of water molecules into specific protein-binding sites. A rigorous statistical mechanical formulation allowing the calculation of the free energy of transfer of water molecules from the bulk to a protein cavity is used to estimate the probabilities of occupancy in the putative bR proton channel. The channel contains a region lined primarily by nonpolar side-chains. Nevertheless, the results indicate that the transfer of four water molecules from bulk water to this apparently hydrophobic region is thermodynamically permitted. The column forms a continuous hydrogen-bonded chain over 12 A between a proton donor, Asp 96, and the retinal Schiff base acceptor. The presence of two water molecules in direct hydrogen-bonding association with the Schiff base is found to be strongly favorable thermodynamically. The implications of these results for the mechanism of proton transfer in bR are discussed.     

Molecular dynamics simulation of six β-blocker drugs passing across POPC bilayer

Six selected β-blocker drugs (alprenolol, atenolol, metoprolol, nadolol, pindolol and propranolol) passing across 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer were studied using all-atom molecular dynamics simulation. The free energy profiles can be divided into two groups, according to their shapes: the free energy curve of group one (atenolol, nadolol and pindolol) has an obvious minimum while that of the other group (propranolol, metoprolol and alprenolol) is flat inside membrane. Energy analysis shows that electrostatic interaction plays an important role for the first group drugs. The hydrogen bond analysis results also certify that the first group drugs form more hydrogen bonds than the other β-blockers. The calculated permeability sequence agrees with the experimental ones. Our calculation suggests that the permeability model using potential of mean force (PMF) method can be also applied to chemically similar compounds besides chemically diverse compounds.     

Effects of salinity on the CO<Subscript>2</Subscript> permeation across lipid bilayer for microalgae biofixation: a molecular dynamics study

The continuous threat of increasing CO₂ concentration in the atmosphere has altered the carbon balance of our planet causing global climate change. Biological fixation of atmospheric CO₂ by unicellular microorganisms such as microalgae is a promising technology pursued extensively by researchers as a means for carbon capture. The study aimed to provide an atomic level of study that will demonstrate the effect of the salinity on the mechanism of CO₂ absorption across microalgae lipid bilayer. Molecular dynamics simulations were utilized to calculate the free energies of CO₂ molecule as it permeates inside the microalgae cell. In thermodynamics, the transport process of a molecule can be demonstrated through its free energy gradient. Thus, calculating the free energies of CO₂ molecule across microalgae lipid bilayer can elucidate the mechanisms of permeation processes. Four microalgae lipid bilayer structures were constructed that contains 128-DPPC (dipalmitoylphosphatidylcholine) lipid bilayer with 3640 water molecules with different NaCl concentrations: 0, 3, 13, and 19 NaCl molecules which correspond to a salinity level of 0, 50, 200, and 300 mM, respectively. The cavity insertion Widom method was used to calculate the free energy of CO₂ molecule along the lipid bilayer. The results demonstrated that the salinity does not affect the free energies significantly, thus, it does not hamper CO₂ transport across microalgae lipid membrane.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.