Next Article in Journal
Membrane Separation Processes in Wastewater and Water Purification, Volume II
Previous Article in Journal
In Situ-Grown Al2O3 Nanoflowers and Hydrophobic Modification Enable Superhydrophobic SiC Ceramic Membranes for Membrane Distillation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study

Departamento de Bioquímica y Biología Molecular-A, Facultad de Veterinaria, Campus de Espinardo, Universidad de Murcia, E-30100 Murcia, Spain
*
Author to whom correspondence should be addressed.
Membranes 2024, 14(6), 118; https://doi.org/10.3390/membranes14060118
Submission received: 26 April 2024 / Revised: 14 May 2024 / Accepted: 16 May 2024 / Published: 21 May 2024
(This article belongs to the Section Biological Membrane Composition and Structures)

Abstract

:
The Danshen terpenoid cryptotanshinone (CPT) is gaining enormous interest in light of its various outstanding biological activities. Among those, CPT has been shown to interact with cell membranes and, for instance, to have antibacterial activity. Several works have shown that CPT alone, or in combination with other drugs, can effectively act as an antibiotic against various infectious bacteria. Some authors have related the mechanism underlying this action to CPT–membrane interaction. This work shows that CPT readily partitions into phosphatidylcholine membranes, but there is a limiting capacity of accommodation of ca. 1 mol CPT to 3 mol phospholipid. The addition of CPT to unilamellar liposomes composed of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) causes membrane permeabilization, as shown by fluorescent probe leakage. This process has been kinetically studied, as well as its modulation by incorporation of phosphatidylethanolamine or phosphatidylglycerol, as a model for pathogenic cell membranes. The thermotropic behavior of 1,2-dimyristoylphosphatidylcholine (DMPC) model membranes is weakly affected by CPT, but the terpenoid causes significant dehydration of the polar region of the bilayer and weak disordering of the acyl chain palisade, as observed in Fourier-transform infrared spectroscopy (FTIR) results. Small-angle X-ray scattering (SAXS) shows that CPT increases DMPC bilayer thickness, which could be due to localization near the phospholipid/water interface. Molecular dynamics (MD) simulations show that the lateral diffusion coefficient of the phospholipid increases with the presence of CPT. CPT extends from the polar head region to the center of the bilayer, being centered between the carbonyl groups and the unsaturated region of the POPC, where there is greater overlap. Interestingly, the free energy profiles of a water molecule crossing the lipid membrane show that the POPC membrane becomes more permeable in the presence of CPT. In summary, our results show that CPT perturbs the physicochemical properties of the phospholipid membrane and compromises its barrier function, which could be of relevance to explain part of its antimicrobial or anticancer activities.

1. Introduction

Tanshinones constitute one of the main components, in terms of their biological and therapeutic activities, of Salvia miltiorrhiza, whose dried roots are known as Danshen [1]. Danshen has been used in traditional Chinese medicine for centuries, and it contains biologically active hydrophilic and hydrophobic compounds, tanshinones, constituting the major lipophilic components [2]. Tanshinones are diterpene compounds involving more than forty specimens to date [3], among which cryptotanshinone (CPT) can be cited as one of the principals, based on its biological activities. CPT is a diterpene quinone (Figure 1) that can permeate cells and present numerous pharmacological activities, particularly for the treatment of cancer and cardiovascular diseases [4]. The relevance of this biologically active compound has been highlighted in a number of excellent reviews [4,5,6,7]. There are multiple pieces of evidence on the efficacy of CPT for the treatment of cancer [5,8], among which we can recently cite its use in the targeting of mutations in colorectal cancer [9], or as a codrug for the treatment of osteosarcoma [10]. Other interesting activities include anti-inflammatory, anti-fibrotic, or neuroprotective effects [6].
The potential molecular mechanisms underlying this vast number of pharmacological and biological actions of CPT often involve different cell signaling pathways that have been recently reviewed [4,6]. In addition to the above-mentioned activities, CPT has been shown to present antibacterial activity [11,12,13,14] whose molecular mechanism could be completely different, being the result of its interaction with the cell membrane, resulting in membrane damage and subsequent cell death [11,12]. However, to date, there are no published results on the interaction of CPT with biological membranes, and the modulation of the structure and function of phospholipid bilayers. Hence, to investigate the effect of the terpenoid on the structure and function of phospholipid membrane bilayers, this work presents an experimental and MD study on the interaction of CPT with model phospholipid vesicles of various compositions, using appropriate physical techniques and permeability assays. The provided results show the potential preferent location of CPT within the bilayer, and how it can lead to membrane alteration and permeabilization, resulting in aqueous contents leakage.

2. Materials and Methods

2.1. Materials

1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-glycerol (POPG), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) were purchased from Avanti Polar Lipids Inc. (Birmingham, AL, USA). Cryptotanshinone (CPT) was purchased from Cymit-Química (Barcelona, Spain). 5(6)-carboxyfluorescein (CF) (99% by HPLC) was from Sigma-Aldrich (Madrid, Spain). All the other reagents were of the highest purity available. The inorganic salts and buffers were of analytical grade. Purified water was deionized in Milli-Q equipment from Millipore (Millipore, Bedford, MA, USA) up to a resistivity of ca. 18 MΩ. Stock solutions of the various phospholipids and CPT were prepared in chloroform and stored at −80 °C. Phospholipid phosphorus was determined according to the method of Böttcher [15]. The buffer used throughout the work was 150 mM NaCl and 5 mM Hepes pH 7.4 unless otherwise stated. Water and all the buffer solutions used in this work were filtered through 0.2 μm filters before use. The osmolarity of all the buffers and solutions was checked using an Osmomat 030 osmometer (Gonotec, Berlin, Germany).

2.2. Partition of CPT into Phospholipid Bilayers

To investigate the partitioning of CPT into phospholipid membranes, POPC was conveniently used since this phospholipid is in the fluid state at room temperature. Mixtures of 2 µmol of POPC and various increasing amounts of CPT were prepared essentially as described below for DSC. Multilamellar liposomes were obtained by adding 1 mL of the above-mentioned buffer and vortexing. These liposomes were left to rest at room temperature for 3 h to allow the sedimentation of any CPT insoluble crystals. A small aliquot of the supernatants (free of CPT crystals, as observed by light microscopy) was carefully collected, and the concentration of phospholipid and CPT was determined. CPT was determined by absorbance at 450 nm after the addition of absolute ethanol (4:1 by volume, ethanol/sample ratio). The molar extinction coefficient of CPT at 450 nm, determined in that same solvent mixture, was ɛ = 2.86 mM−1 cm−1.

2.3. Vesicle Contents Release

CPT-induced vesicle contents release was monitored with the CF assay, for which CF was entrapped within the phospholipid unilamellar vesicles of the indicated composition, and its leakage was followed by the increase in fluorescence due to dilution of the probe to the external medium [16]. Multilamellar vesicles were prepared by the dry film hydration method, by vortexing 5 μmol of POPC with 0.5 mL of a buffer containing 50 mM CF and 10 mM Hepes, pH 7.4, at room temperature, and large unilamellar vesicles were obtained by 21 times extrusion of the multilamellar vesicles through two stacked polycarbonate filters (0.1 μm pore diameter). Sephadex G-50 gel filtration, with 100 mM NaCl and 10 mM Hepes pH 7.4 as elution buffer, was used to separate the vesicles from non-encapsulated CF. CPT was added at the indicated concentrations from stock solutions in DMSO, and it was confirmed that the same volumes of pure DMSO did not induce any leakage. Maximum leakage (100%) was established by disrupting the liposomes with 5% Triton X100, which provoked the complete release of CF to the external medium. The percentage of CF leakage was calculated as:
%   C F   L e a k a g e = F t F i · 100 F d F i
where Ft is the fluorescence at a given time after CPT injection, Fi is the initial fluorescence, and Fd is the maximum fluorescence obtained after addition of the detergent Triton X100.

2.4. Differential Scanning Calorimetry

Multilamellar vesicles for DSC measurements were prepared by the dry film hydration method. Briefly, 2 µmol of DMPC or DPPC and CPT as indicated were mixed in chloroform, and the organic solvent was gently evaporated using dry N2 to obtain a thin film at the bottom of a glass tube. The last traces of solvent were removed by a further 3 h desiccation under a high vacuum. To the dry samples, 2 mL of a buffer containing 100 mM NaCl, 0.1 mM EDTA, and 10 mM Hepes pH 7.4 was added, and vesicles were formed by vortexing the mixture at 50 °C. The experiments were performed using a MicroCal MC2 calorimeter (MicroCal, Northampton, NC, USA) at 1 mM phospholipid concentration and 60 °C h−1 heating scan rate. Three consecutive heating scans were carried out for each sample, the last one being taken for analysis.

2.5. Fourier-Transform Infrared Spectroscopy

Multilamellar vesicles for Fourier-transform infrared spectroscopy (FTIR) were prepared by the dry film hydration method, in 40 μL of the same buffer described above prepared with D2O. The samples were placed between two CaF2 windows (25 × 2 mm) separated by 25 μm Teflon spacers and mounted onto a Symta cell mount. Infrared spectra were acquired in a Nicolet 6700 Fourier-transform infrared spectrometer (Madison, WI, USA). Each spectrum was obtained by collecting 64 interferograms with a nominal resolution of 2 cm−1. The equipment was continuously purged with dry air to minimize the contribution peaks of atmospheric water vapor, and the sample holder was thermostated using a Peltier device (Nicolet Proteus System). Spectra were collected at 2 °C intervals, allowing 5 min equilibration between temperatures. The D2O buffer spectra taken at the same temperatures were subtracted interactively using either Omnic (8.0.342) or Grams (7.02) (Galactic Industries, Salem, NH, USA) software.

2.6. X-ray Scattering

Samples for X-ray scattering analysis were prepared essentially as described above for FTIR, in 1 mL of 100 mM NaCl, 0.1 mM EDTA, and 10 mM Hepes pH 7.4 at 40 °C. The samples were centrifuged in a bench microfuge, and the pellets were placed in the sample holder of the diffractometer with the aid of a spatula. A steel holder with cellophane windows was used, providing good thermal contact to the Peltier heating unit. Typical exposure times were 5 min, leaving a 10 min equilibration period prior to each measurement. Small angle (SAXS) X-ray scattering data were collected using a Kratky compact camera (M. Braun-Graz Optical Systems, Graz, Austria), equipped with a linear position sensitive detector (PSD; M. Braun, Garching, Germany), monitoring the s range (s = 2 sin θ/λ, 2θ = scattering angle, λ = 1.54 Å) between 0.0075 and 0.07 Å−1. Nickel-filtered Cu Kα X-rays were generated by a Philips PW3830 X-ray generator (Eindhoven, The Netherlands) at 50 kV and 30 mA. The calibration of the detector position was performed by using silver stearate (d-spacing at 48.8 Å) as the reference material.

2.7. Molecular Dynamics Simulations

The 3D molecular structure of CPT was obtained from the PubChem Substance and Compound databases [17] through the unique chemical structure identifier CID 160254. All the MD simulations were conducted using GROMACS 5.0.7 and 2018.1 [18]. The CHARMM36 force field parameters for DMPC, POPC, CPT, water, Cl, and Na+ were obtained from CHARMM-GUI [19,20,21]. The membrane bilayers were formed by 2 leaflets oriented normal to the z-axis with a total of 128 molecules of DMPC or POPC with and without 14 molecules of CPT, and a water layer containing a total of 6400 water molecules (TIP3 model), 12 Na+, and 12 Cl. The initial membrane structures were built with the aid of Packmol software (v 20.14.4) [22].
All the systems were simulated using the NpT-ensemble at 312 K for DMPC membranes and at 298 K for POPC membranes. The pressure was controlled semi-isotropically at a pressure of 1 bar and compressibility of 4.5 × 10−5 bar−1. The cutoffs for van der Waals and short-range electrostatic interactions were 1.2 nm, and a force-switching function was applied between 1.0 and 1.2 nm [23]. The simulations were initiated by a 20 ns run, using the V-rescale thermostat and the Berendsen barostat [24], followed by a 200 ns run using the Nosé–Hoover thermostat [25] and the Parrinello–Rahman barostat [26]. Graphical representations were conducted with PyMOL 2.3.0 [27]. Analyses of the trajectories were conducted over the last 60 ns using the Gromacs tools.
Area per lipid was calculated from the lateral dimensions of the Malone simulator box (the area of the xy plane) divided by the number of lipids in each leaflet. The diffusion constants (D) were calculated from the slope of the mean square displacement (MSD) versus time using the Einstein relation:
D = lim Δ t Δ r ¯ t 0 2 2 d Δ t
where MSD = | r ¯ | 2 t 0 , and d is the dimensions (e.g., d = 2 for lateral diffusion). The deuterium order parameter across the acyl chain, SCD, is defined as:
S CD = 3 2 cos 2 θ 1 2
where θ is the angle between a CD bond and the bilayer normal. A value of −0.5 indicates a perfectly ordered acyl chain in an all-trans conformation [28]. The membrane thickness was computed by calculating the phosphorous atom’s distance between both leaflets. The cluster size distribution of CPT in the membranes was calculated as the number of CPT molecules found within 0.30 nm of the analyzed trajectory.
The free energy profile of an H2O molecule crossing the POPC bilayer, in the absence and the presence of CPT molecules, was calculated using the potential of mean force (PMF) procedure [29]. The starting configurations were the final configurations obtained in the MD simulations. An H2O molecule from the water phase was pulled to the lipid bilayer center along the z-axis using the umbrella method [30]. A harmonic restraint of 1000 kJ mol−1 nm−2 and a pulling rate of 0.007 nm ps−1 were applied to distance z between the center of mass of the H2O molecule and the POPC bilayer. This procedure was repeated three times with different water molecules without altering the membrane structure. The H2O molecule configuration was then sampled and constrained at different z distances, allowing free motion in the xy plane. The free energy profile was calculated by using the weighted histogram analysis method [31] included in GROMACS tools. The calculated free energy profile was considered symmetric across the POPC bilayer center.

3. Results and Discussion

CPT, a diterpene quinone, can be considered a phenanthraquinone derivative (Figure 1); it is an essentially water-insoluble compound but is soluble in most organic solvents. This highly hydrophobic character should make CPT very prompt to partition into phospholipid membranes and, in fact, CPT has been described as a cell-permeable compound [8]. In this work, we carried out an experimental MD approach to study the incorporation of CPT into model phospholipid membranes, and its effect on the structure and function of model phospholipid membranes of various compositions, to obtain information on CPT–phospholipid interactions at the molecular level. Up to five different phospholipid species were used. To study thermotropic transitions, multilamellar vesicles made of DMPC or DPPC were used, since these phospholipids present thermal transitions in a temperature range appropriate for this type of study. These systems were analyzed using DSC, FTIR, and SARS techniques. On the other hand, for the study of CPT-induced membrane permeabilization, through CF fluorescence, POPC, POPG, and POPE were used, since these phospholipids, containing an oleic acid at the sn-2 position of glycerol, form fluid bilayers above 0 °C.

3.1. Partition of CPT into POPC Membranes

Over the course of the study, it was observed that upon the formation of CPT/POPC liposomes, a fraction of the diterpene did not incorporate into the vesicles and appeared in the form of orange crystals, as observed by light microscopy. Since CPT is essentially water-insoluble, it was concluded that the fraction of terpenoid non-incorporated into the membrane readily crystallized.
This finding led to the need to determine the actual concentration of CPT incorporated into the membrane. Mixtures of POPC and CPT at various initial proportions were prepared, and the incorporated CPT was determined (Figure 2). The actual CPT/POPC membrane molar ratio showed a hyperbolic dependence with the total initial concentration of CPT. Thus, the data shown in Figure 2 were adjusted to the following equation:
R m e m = R s a t · C P T t o t P O P C R 50 + C P T t o t P O P C
where Rmem is the molar ratio of CPT in the membrane to total POPC, Rsat is the saturation membrane molar ratio of CPT, R50 is the CPT concentration necessary to reach 50% of Rsat, and CPTtot is the total CPT concentration. This fitting yielded a value for Rsat of 0.30 ± 0.04 and R50 of 1.2 ± 0.3. These interesting results indicated that the fluid membrane of POPC had a limiting capacity on incorporating CPT, corresponding to a CPTmem/POPC molar ratio of ca. 0.33 (the non-incorporated diterpene crystallized in the form of insoluble crystals, as commented above). According to these results, it was clear that POPC membranes had a limited capacity to accommodate CPT, which was around 1 mol CPT to 3 mol phospholipid. This conclusion can be extrapolated to any phosphatidylcholine membrane, since POPC is widely accepted and used as an average representative for fluid phosphatidylcholine membranes [32] and, in particular, for membrane permeability studies [33]. Therefore, it is particularly important to take this fact into consideration when the membrane-related effects of CPT are evaluated, since the actual concentration of this bioactive compound in the membrane is significantly smaller than that added.

3.2. CPT-Induced Membrane Permeabilization

Figure 3 shows the time course curves of CPT-induced content leakage of POPC large unilamellar vesicles (LUV). Increasing concentrations of CPT were added from outside, as a DMSO solution, to pre-formed vesicles at a constant concentration, and the leakage of CF was monitored by the concomitant increase in probe fluorescence. Concentrations of the diterpenoid as low as 4 µM (CPTtot/phospholipid molar ratio = 0.2) already gave rise to a slow but measurable liposome content leakage to the external medium, which progressed to reach 100% if allowed to continue.
The rate and extent of CPT-induced leakage continued progressively increasing as its concentration was raised above that ratio, indicating that the Rsat value was not reached. These results indicate that a fraction of added CPT rapidly incorporated into the POPC membranes and compromised their barrier function, resulting in increased permeability to water-soluble compounds, like CF.
The influence of phospholipid membrane composition on CPT-induced leakage was checked in liposomes containing an addition of either POPG or POPE (Figure 4). POPE was selected because of its small headgroup and rich lipid polymorphism [34], whereas POPG has an anionic character and a voluminous headgroup [35].
It was observed that, irrespective of CPT concentration, the presence of POPE protected the membrane against permeabilization, whereas POPG significantly enhanced CPT-induced vesicle content leakage. This effect was observed both at low and high CPT/phospholipid ratios, suggesting that it is the nature of the bilayer itself, and not specific phospholipid–CPT interactions that account for it. This is a very interesting result since it showed that target membrane composition determined the effect of CPT. It means that CPT membrane action is not indiscriminate, and it is expected to present different effects on different cells, as commented below.
The protective effect of phosphatidylethanolamine is not new, since this phospholipid has been reported to inhibit leakage in other systems [36,37,38,39,40]. Due to its small headgroup and the interlipid hydrogen bonding initiated by the amine group [41], phosphatidylethanolamine can increase membrane compactness and stabilize the surface state of the phospholipid membrane [36], which would impede the proper insertion of CPT necessary for permeabilization, thus reducing leakage.
The significant enhancement in CF leakage by POPG is of great importance, given the fact that both bacterial and eukaryotic cancer cells have negatively charged phospholipids in the outer monolayers [42,43,44]. Thus, it seems that negatively charged membranes, as a model of pathogenic membranes [45], are more prompt to permeabilization, which might contribute to explaining part of the anticancer and antimicrobial actions of CPT.

3.3. Modulation of DMPC Thermotropic Transitions by CPT

To obtain information on CPT–phospholipid interactions, a saturated species of phosphatidylcholine, namely DMPC, was chosen. This phospholipid is appropriate for these types of studies since, upon heating, it displays the so-called pre-transition, from the gel Lβ phase to the gel ripple Pβ′ phase, at around a Tc of 12.1 °C, and a main transition from the ripple Pβ′ phase to the liquid crystalline fluid, Lα, phase around a Tc of 23.7 °C (Figure 5).
CPT was incorporated into DMPC liposomes by comixing prior to vesicle formation, at two different concentrations. It can be observed that the presence of CPTtot/DMPC 0.1 and 0.2, respectively, of the quinone had a significant effect on the pre-transition, which was shifted toward lower temperatures, from 12.1 to 9.3 and 8.0 °C, respectively, progressively widened, and decreased in area. However, the effect on the main gel-to-liquid crystalline phase transition was weak, the Tc shifting from 23.7 to 22.8 and 22.7 °C, 0.1 and 0.2 CPTtot/DMPC, respectively. In addition, the main transition peak was slightly widened. It was clear that the membrane-incorporated CPT did not establish sufficiently strong interactions with DMPC to significantly perturb its thermotropic behavior. Similar measurements were carried out using DPPC, and the results were qualitatively like those described above for DMPC. Since the pre-transition is mostly caused by the rearrangement of lipid headgroups, and the main transition by the melting of lipid acyl chains, our results showed a stronger effect of CPT at the level of the phospholipid polar region and a much smaller effect on lipid packing.

3.4. Effect of CPT on Acyl Chain and Polar Headgroup Regions of the DMPC Bilayer

Using FTIR, it was possible to elucidate the potential effects of CPT on the two main moieties of the DMPC molecule: the acyl chain palisade and the lipid/water interface. Thus, the acyl chain region was studied by monitoring the frequency of the CH2 symmetric stretching band, νCH2, and the polar region from the frequency of the C=O stretching band, νCO (Figure 6). The incorporation of 10 mol% of CPT into the DMPC resulted in a ~1 cm−1 shift of the maximum of the νCH2 band towards higher values, in the whole temperature range. Further increasing the concentration of CPT to 20 mol% gave rise to an additional small increase of barely 0.3 cm−1, particularly in the fluid phase. These weak effects were in good agreement with the weak effect of CPT on the main gel-to-liquid crystalline phase transition of DMPC described above. In the case of the νCO band, the inclusion of 10 mol% of CPT shifted the maximum of the band about 2.4 cm−1 towards higher values, and, again, 20 mol% of the quinone just slightly enhanced this shift. This indicates that CPT was causing dehydration of the polar region of the membrane, in line with the effect of the terpenoid on the DMPC pre-transition shown above.
A similar study was conducted with liposomes composed of various phosphatidylcholines, and the main results are presented in Table 1. It was observed that in the case of DPPC, with longer acyl chains than DMPC, the effects were similar but significantly weaker. Furthermore, in the case of POPC, which is fluid at 25 °C, CPT also caused an additional shift in the maximum of both bands, as in the case of DMPC.
Taken together, the DSC and FTIR results clearly indicated that CPT did not significantly affect lipid packing, its effects being concentrated at the polar part of the bilayer.

3.5. Effect of CPT on DMPC Structural Parameters

SAXS was carried out to obtain information on the influence of CPT on the DMPC bilayers’ structural parameters. Measurements were carried out at three different temperatures corresponding to the Lβ phase (7 °C), the gel ripple Pβ′ phase (18 °C), and the liquid crystalline Lα phase (34 °C) (Figure 7). Pure DMPC showed three reflections with relative distances at 1:1/1:2/1:3, corresponding to a multilamellar organization. The largest peak (the first-order reflection) gives the interlamellar repeat distance (d-spacing), a summation of the bilayer thickness, and the thickness of the water layer between bilayers [46]. The first-order reflection of pure DMPC appeared at a d-spacing of 59.6 Å at 7 °C, 64.6 Å at 18 °C, and 60.0 Å at 34 °C, in agreement with previous data [47]. The incorporation of 0.1 CPTtot/DMPC resulted in an increase in d-spacing below the main gel-to-liquid crystalline phase transition: 64.6 Å at 7 °C and 66.3 Å at 18 °C, whereas it did not change in the Lα phase: 60.5 Å at 34 °C. This widening effect could occur by an increase in the water layer between the DMPC membranes, or by a real increase in the bilayer thickness, due to the insertion of CPT.

3.6. Molecular Dynamics Simulations

The area per lipid parameter is frequently used as a property of a lipid bilayer for validating MD simulations [48]. Figure 8 shows the time course of the area per lipid in different simulations.
It can be observed that the area per lipid remained mainly constant during the MD simulation time. The area per lipid of pure POPC and DMPC bilayers above the phase transition was 0.63 nm2 and 0.61 nm2, respectively (Table 2). Thus, our results are in good agreement with the reported data for POPC (0.61–0.64 nm2) [49,50] and for DMPC (0.61 nm2) [50,51]. In the presence of CPT, the area per lipid increased by around 5% in both membranes, probably due to the broadening produced solely by the presence of more non-phospholipid molecules in the membranes. Neither the order parameter (SCD) nor the thickness of the membranes changed significantly in the presence of CPT (Table 2), in total agreement with the SAXS data shown above.
In the presence of CPT, the lateral diffusion coefficient of the phospholipid increased by 41% and 57% in POPC and DMPC membranes, respectively (Table 2). This result indicated a higher molecular mobility in the membrane in the presence of CPT, which can result in local membrane destabilization (‘pores’ formation) and leakage. In this case, the term ‘pore’ within should be understood as zones in the membrane where membrane stability is decreased, resulting in membrane permeabilization [52].
Figure 9 shows the mass density profile of the POPC membrane in the presence of CPT. Some important groups of the POPC molecule have been included: the polar head region (P atoms), the carbonyl groups, and the location of the terminal methyl of sn-1 chains of POPC. CPT extended from the polar headgroup region up to near the center of the bilayer (z = 0), but there was more overlapping with the carbonyl groups profile (phospholipid/water interface) than with that of the terminal carbons of POPC acyl chains, indicating that the terpenoid was not deeply inserted within the bilayer.
The study of cluster formation of CPT in the membrane showed that CPT molecules are mostly found in the form of monomers, around 66% in POPC and DMPC membranes (Table 2), a result that might seem in contrast to the high hydrophobicity of this molecule. However, observing its molecular structure, the intermolecular interactions would result in steric hindrance, hence the presence of monomers was favored, and only clusters of a few molecules were found. This result excludes CPT clusters as the possible responsible for membrane permeabilization.
Figure 10 shows a snapshot corresponding to a final configuration of POPC + CPT membrane simulation. It was observed that clustering of CPT was minimum, and most molecules were monomers or dimers, as shown above. This localization of CPT, closer to the POPC carbonyl groups than to the terminal methyl, was not expected to produce large perturbations in the thermotropic behavior of these systems, as shown above by DSC, due to a minimum effect on lipid packing. However, it will explain the stronger effect on the region of the C=O groups, where significant dehydration was observed by FTIR.
To simulate the effect of CPT on the permeability properties of the membrane to polar substances, the PMF of a water molecule crossing the membrane was determined. The free energy profiles of such a process are shown in Figure 11. The free energy was set to zero in the aqueous phase. In the POPC membrane, the free energy for a water molecule to cross the membrane was about 6.8 ± 1.37 kcal/mol, which is in good agreement with previously reported values ranging from 6.2 to 7 kcal/mol [48,53,54].
It can be observed that in the presence of CPT, the free energy decreased to 4.84 ± 1.22 kcal/mol, clearly indicating that the POPC membrane became more permeable to water in the presence of CPT, making it likely that, by extrapolation, it was also more permeable in general to polar water-soluble compounds. These MD results agree with our experimental data on CF leakage shown above since CF is also a highly polar water-soluble compound, which easily permeates upon the addition of CPT.

4. Conclusions

The hydrophobic character of CPT suggests it will exert a strong influence on phospholipid membrane structure and function. We have shown that upon addition to phospholipid liposomes, a fraction of CPT is incorporated into the membrane, resulting in membrane perturbation and content leakage. The finding that CPT presents a limited incorporation into fluid phospholipid membranes, and its tendency to crystallize in water, is an important fact to be taken into consideration relative to its pharmacological use. The effect of CPT mainly occurs at the level of the phospholipid polar headgroups, with little influence on lipid packing, according to FTIR. Thus, CPT increases POPC mobility and dehydrates the polar headgroup region, which could explain membrane permeabilization. The most relevant finding is that CPT membrane effects depend on the membrane phospholipid composition. In particular, the presence of negatively charged phosphatidylglycerol, as a model for cancer and bacterial cells, enhances the activity of CPT, which is important for the explanation of some of its pharmacological actions, since either bacterial or cancer cells normally present negatively charged phospholipids at the outer monolayer. The described effects of CPT on model membranes add to the vast number of biological actions of the terpenoid and enhance its potential, both as an antibacterial agent and for addressing pathogenic tumor cells.

Author Contributions

J.O.: data curation, validation, investigation; F.J.A.: validation; J.A.T.: validation, data curation; A.O.: conceptualization, supervision, resources, formal analysis, validation, writing—original draft, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the Computational Service of the University of Murcia for allowing the use of their facilities for the MD calculations.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jiang, Z.; Gao, W.; Huang, L. Tanshinones, Critical Pharmacological Components in Salvia miltiorrhiza. Front. Pharmacol. 2019, 10, 439237. [Google Scholar] [CrossRef] [PubMed]
  2. Ke, L.; Zhong, C.; Chen, Z.; Zheng, Z.; Li, S.; Chen, B.; Wu, Q.; Yao, H. Tanshinone I: Pharmacological activities, molecular mechanisms against diseases and future perspectives. Phytomedicine 2023, 110, 154632. [Google Scholar] [CrossRef] [PubMed]
  3. Gao, H.; Huang, L.; Ding, F.; Yang, K.; Feng, Y.; Tang, H.; Xu, Q.M.; Feng, J.; Yang, S. Simultaneous purification of dihydrotanshinone, tanshinone I, cryptotanshinone, and tanshinone IIA from Salvia miltiorrhiza and their anti-inflammatory activities investigation. Sci. Rep. 2018, 8, 8460. [Google Scholar] [CrossRef] [PubMed]
  4. Wu, Y.R.Y.H.; Wu, Y.R.Y.H.; Li, B.; Yan, Z.Y. Cryptotanshinone: A review of its pharmacology activities and molecular mechanisms. Fitoterapia 2020, 145, 104633. [Google Scholar] [CrossRef] [PubMed]
  5. Ashrafizadeh, M.; Zarrabi, A.; Orouei, S.; Saberifar, S.; Salami, S.; Hushmandi, K.; Najafi, M. Recent advances and future directions in anti-tumor activity of cryptotanshinone: A mechanistic review. Phyther. Res. 2021, 35, 155–179. [Google Scholar] [CrossRef] [PubMed]
  6. Li, H.; Gao, C.; Liu, C.; Liu, L.; Zhuang, J.; Yang, J.; Zhou, C.; Feng, F.; Sun, C.; Wu, J. A review of the biological activity and pharmacology of cryptotanshinone, an important active constituent in Danshen. Biomed. Pharmacother. 2021, 137, 111332. [Google Scholar] [CrossRef]
  7. Wang, X.; Yang, Y.; Liu, X.; Gao, X. Pharmacological Properties of Tanshinones, the Natural Products from Salvia miltiorrhiza, 1st ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2020; Volume 87. [Google Scholar]
  8. Chen, W.; Lu, Y.; Chen, G.; Huang, S. Molecular Evidence of Cryptotanshinone for Treatment and Prevention of Human Cancer. Anticancer Agents Med. Chem. 2013, 13, 979–987. [Google Scholar] [CrossRef]
  9. Vundavilli, H.; Datta, A.; Sima, C.; Hua, J.; Lopes, R.; Bittner, M. Targeting oncogenic mutations in colorectal cancer using cryptotanshinone. PLoS ONE 2021, 16, e0247190. [Google Scholar] [CrossRef] [PubMed]
  10. Saraf, R.; Datta, A.; Sima, C.; Hua, J.; Lopes, R.; Bittner, M.L.; Miller, T.; Wilson-Robles, H.M. In Silico Modeling of the Induction of Apoptosis by Cryptotanshinone in Osteosarcoma Cell Lines. IEEE/ACM Trans. Comput. Biol. Bioinform. 2022, 19, 1683–1693. [Google Scholar] [CrossRef]
  11. Chen, B.C.; Ding, Z.S.; Dai, J.S.; Chen, N.P.; Gong, X.W.; Ma, L.F.; Qian, C.D. New Insights Into the Antibacterial Mechanism of Cryptotanshinone, a Representative Diterpenoid Quinone from Salvia miltiorrhiza Bunge. Front. Microbiol. 2021, 12, 647289. [Google Scholar] [CrossRef]
  12. Wang, D.; Zhang, W.; Wang, T.; Li, N.; Mu, H.; Zhang, J.; Duan, J. Unveiling the mode of action of two antibacterial tanshinone derivatives. Int. J. Mol. Sci. 2015, 16, 17668–17681. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, D.; Lee, S.; Noh, J.; Hong, S. Antibacterial Activities of Cryptotanshinone and Dihydrotanshinone I from a Medicinal Herb, Salvia miltiorrhiza Bunge. Biosci. Biotechnol. Biochem. 1999, 63, 2236–2239. [Google Scholar] [CrossRef] [PubMed]
  14. Cha, J.D.; Lee, J.H.; Choi, K.M.; Choi, S.M.; Park, J.H. Synergistic effect between cryptotanshinone and antibiotics against clinic methicillin and vancomycin-resistant Staphylococcus aureus. Evid.-Based Complement. Altern. Med. 2014, 2014, 450572. [Google Scholar] [CrossRef] [PubMed]
  15. Böttcher, C.J.F.; van Gent, C.M.; Pries, C. A rapid and sensitive sub-micro phosphorus determination. Anal. Chim. Acta 1961, 24, 203–204. [Google Scholar] [CrossRef]
  16. Coronel, J.R.; Aranda, F.J.; Teruel, J.A.; Marqués, A.; Manresa, Á.; Ortiz, A. Kinetic and Structural Aspects of the Permeabilization of Biological and Model Membranes by Lichenysin. Langmuir 2016, 32, 78–87. [Google Scholar] [CrossRef] [PubMed]
  17. Kim, S.; Thiessen, P.A.; Bolton, E.E.; Chen, J.; Fu, G.; Gindulyte, A.; Han, L.; He, J.; He, S.; Shoemaker, B.A.; et al. PubChem Substance and Compound databases. Nucleic Acids Res. 2016, 44, D1202–D1213. [Google Scholar] [CrossRef] [PubMed]
  18. Abraham, M.J.; van der Spoel, D.; Lindahl, E.; Hess, B.; the GROMACS Developent Team. GROMACS User Manual Version 2018. 2018. Available online: www.gromacs.org (accessed on 26 April 2024).
  19. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  20. Brooks, B.R.; Brooks, C.L.; Mackerell, A.D.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef] [PubMed]
  21. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef]
  22. Martinez, L.; Andrade, R.; Birgin, E.G.; Martínez, J.M. PACKMOL: A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157–2164. [Google Scholar] [CrossRef]
  23. Bjelkmar, P.; Larsson, P.; Cuendet, M.A.; Hess, B.; Lindahl, E. Implementation of the CHARMM Force Field in GROMACS: Analysis of Protein Stability Effects from Correction Maps, Virtual Interaction Sites, and Water Models. J. Chem. Theory Comput. 2010, 6, 459–466. [Google Scholar] [CrossRef] [PubMed]
  24. Berendsen, H.J.C.; Postma, J.P.M.; Van Gunsteren, W.F.; Dinola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  25. Hoover, W.G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [PubMed]
  26. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  27. Schrödinger, L. PyMOL. 2015. Available online: https://github.com/schrodinger/pymol-open-source (accessed on 26 April 2024).
  28. Seelig, J. Deuterium magnetic resonance: Theory and application to lipid membranes. Q. Rev. Biophys. 1977, 10, 353–418. [Google Scholar] [CrossRef] [PubMed]
  29. Woo, H.J.; Roux, B. Calculation of absolute protein-ligand binding free energy from computer simulations. Proc. Natl. Acad. Sci. USA 2005, 102, 6825–6830. [Google Scholar] [CrossRef] [PubMed]
  30. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  31. Hub, J.S.; De Groot, B.L.; Van Der Spoel, D. g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. J. Chem. Theory Comput. 2010, 6, 3713–3720. [Google Scholar] [CrossRef]
  32. Oberholzer, T.; Luisi, P.L. The Use of Liposomes for Constructing Cell Models. J. Biol. Phys. 2002, 28, 733–744. [Google Scholar] [CrossRef]
  33. Nasr, G.; Greige-Gerges, H.; Elaissari, A.; Khreich, N. Liposomal Membrane Permeability Assessment by Fluorescence Techniques: Main Permeabilizing Agents, Applications and Challenges; Elsevier B.V.: Amsterdam, The Netherlands, 2020; Volume 580, p. 119198. [Google Scholar]
  34. van der Veen, J.N.; Kennelly, J.P.; Wan, S.; Vance, J.E.; Vance, D.E.; Jacobs, R.L. The critical role of phosphatidylcholine and phosphatidylethanolamine metabolism in health and disease. Biochim. Biophys. Acta—Biomembr. 2017, 1859, 1558–1572. [Google Scholar] [CrossRef]
  35. Manna, M.; Mukhopadhyay, C. Molecular dynamics simulations of the interactions of kinin peptides with an anionic POPG bilayer. Langmuir 2011, 27, 3713–3722. [Google Scholar] [CrossRef] [PubMed]
  36. Nicol, F.; Nir, S.; Szoka, F.C. Effect of phospholipid composition on an amphipathic peptide-mediated pore formation in bilayer vesicles. Biophys. J. 2000, 78, 818. [Google Scholar] [CrossRef] [PubMed]
  37. Yan, L.; Miao, Q.; Sun, Y.; Yang, F. tBid forms a pore in the liposome membrane. FEBS Lett. 2003, 555, 545–550. [Google Scholar] [CrossRef] [PubMed]
  38. Zaragoza, A.; Aranda, F.J.; Espuny, M.J.; Teruel, J.A.; Marqués, A.; Manresa, A.; Ortiz, A. Mechanism of membrane permeabilization by a bacterial trehalose lipid biosurfactant produced by Rhodococcus sp. Langmuir 2009, 25, 7892–7898. [Google Scholar] [CrossRef]
  39. Sánchez, M.; Aranda, F.J.; Teruel, J.A.; Espuny, M.J.; Marqués, A.; Manresa, Á.; Ortiz, A. Permeabilization of biological and artificial membranes by a bacterial dirhamnolipid produced by Pseudomonas aeruginosa. J. Colloid Interface Sci. 2010, 341, 240–247. [Google Scholar] [CrossRef] [PubMed]
  40. Fiedler, S.; Heerklotz, H. Vesicle Leakage Reflects the Target Selectivity of Antimicrobial Lipopeptides from Bacillus subtilis. Biophys. J. 2015, 109, 2079–2089. [Google Scholar] [CrossRef] [PubMed]
  41. Boggs, J.M. Lipid intermolecular hydrogen bonding: Influence on structural organization and membrane function. Biochim. Biophys. Acta—Rev. Biomembr. 1987, 906, 353–404. [Google Scholar] [CrossRef] [PubMed]
  42. Ran, S.; Downes, A.; Thorpe, P.E. Increased Exposure of Anionic Phospholipids on the Surface of Tumor Blood Vessels. Cancer Res. 2002, 62, 6132–6140. [Google Scholar]
  43. Travkova, O.G.; Moehwald, H.; Brezesinski, G. The interaction of antimicrobial peptides with membranes. Adv. Colloid Interface Sci. 2017, 247, 521–532. [Google Scholar] [CrossRef]
  44. Teixeira, V.; Feio, M.J.; Bastos, M. Role of lipids in the interaction of antimicrobial peptides with membranes. Prog. Lipid Res. 2012, 51, 149–177. [Google Scholar] [CrossRef]
  45. Battista, F.; Oliva, R.; Del Vecchio, P.; Winter, R.; Petraccone, L. Insights into the Action Mechanism of the Antimicrobial Peptide Lasioglossin III. Int. J. Mol. Sci. 2021, 22, 2857. [Google Scholar] [CrossRef] [PubMed]
  46. Rappolt, M.; Hickel, A.; Bringezu, F.; Lohner, K. Mechanism of the Lamellar/Inverse Hexagonal Phase Transition Examined by High Resolution X-ray Diffraction. Biophys. J. 2003, 84, 3111. [Google Scholar] [CrossRef] [PubMed]
  47. Ortiz, A.; Teruel, J.A.; Espuny, M.J.; Marqués, A.; Manresa, A.; Aranda, F.J. Effects of dirhamnolipid on the structural properties of phosphatidylcholine membranes. Int. J. Pharm. 2006, 325, 99–107. [Google Scholar] [CrossRef] [PubMed]
  48. Venable, R.M.; Krämer, A.; Pastor, R.W. Molecular Dynamics Simulations of Membrane Permeability. Chem. Rev. 2019, 119, 5954–5997. [Google Scholar] [CrossRef] [PubMed]
  49. Kučerka, N.; Nieh, M.P.; Katsaras, J. Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. Biochim. Biophys. Acta—Biomembr. 2011, 1808, 2761–2771. [Google Scholar] [CrossRef] [PubMed]
  50. Drabik, D.; Chodaczek, G.; Kraszewski, S.; Langner, M. Mechanical Properties Determination of DMPC, DPPC, DSPC, and HSPC Solid-Ordered Bilayers. Langmuir 2020, 36, 3826–3835. [Google Scholar] [CrossRef] [PubMed]
  51. Kučerka, N.; Liu, Y.; Chu, N.; Petrache, H.I.; Tristram-Nagle, S.; Nagle, J.F. Structure of Fully Hydrated Fluid Phase DMPC and DLPC Lipid Bilayers Using X-ray Scattering from Oriented Multilamellar Arrays and from Unilamellar Vesicles. Biophys. J. 2005, 88, 2626. [Google Scholar] [CrossRef] [PubMed]
  52. Ortiz, J.; Teruel, J.A.; Aranda, F.J.; Ortiz, A. On the Mechanism of Membrane Permeabilization by Tamoxifen and 4-Hydroxytamoxifen. Membranes 2023, 13, 292. [Google Scholar] [CrossRef]
  53. Comer, J.; Schulten, K.; Chipot, C. Calculation of Lipid-Bilayer Permeabilities Using an Average Force. J. Chem. Theory Comput. 2014, 10, 554–564. [Google Scholar] [CrossRef]
  54. Sajadi, F.; Rowley, C.N. Simulations of lipid bilayers using the CHARMM36 force field with the TIP3P-FB and TIP4P-FB water models. PeerJ 2018, 6, e5472. [Google Scholar] [CrossRef]
Figure 1. The chemical structure of cryptotanshinone.
Figure 1. The chemical structure of cryptotanshinone.
Membranes 14 00118 g001
Figure 2. Partition of CPT into POPC membranes. The actual molar ratio of CPT in the membrane (CPTmem/POPC), determined as described in Materials and Methods, is plotted against the total initial CPT to POPC molar ratio (CPTtot/POPC). Data correspond to the mean of three independent experiments ± SD (error bars).
Figure 2. Partition of CPT into POPC membranes. The actual molar ratio of CPT in the membrane (CPTmem/POPC), determined as described in Materials and Methods, is plotted against the total initial CPT to POPC molar ratio (CPTtot/POPC). Data correspond to the mean of three independent experiments ± SD (error bars).
Membranes 14 00118 g002
Figure 3. Content leakage curves for CPT-induced permeabilization of POPC LUV. The concentration of POPC was kept constant at 20 µM. Numbers on the curves indicate the molar ratio of total CPT to phospholipid in the cuvette (CPTtot/POPC).
Figure 3. Content leakage curves for CPT-induced permeabilization of POPC LUV. The concentration of POPC was kept constant at 20 µM. Numbers on the curves indicate the molar ratio of total CPT to phospholipid in the cuvette (CPTtot/POPC).
Membranes 14 00118 g003
Figure 4. The influence of membrane composition on CPT-induced leakage. Experiments were conducted under the same conditions as in Figure 3, for liposomes composed of pure POPC (black), POPC/POPG (5:1.7, mol/mol) (blue), and POPC/POPE (5:1.7, mol/mol) (red), at two different CPT concentrations, as indicated.
Figure 4. The influence of membrane composition on CPT-induced leakage. Experiments were conducted under the same conditions as in Figure 3, for liposomes composed of pure POPC (black), POPC/POPG (5:1.7, mol/mol) (blue), and POPC/POPE (5:1.7, mol/mol) (red), at two different CPT concentrations, as indicated.
Membranes 14 00118 g004
Figure 5. High-sensitivity DSC heating thermograms for mixtures of CPT with DMPC. CPTtot/DMPC molar ratios are indicated on the curves. Scans were carried out at 60 °C h−1.
Figure 5. High-sensitivity DSC heating thermograms for mixtures of CPT with DMPC. CPTtot/DMPC molar ratios are indicated on the curves. Scans were carried out at 60 °C h−1.
Membranes 14 00118 g005
Figure 6. The effect of CPT on the region of the DMPC acyl chains and polar headgroups, determined by FTIR. Top: the effect on the maximum frequency of the νCH2 symmetric stretching band. Bottom: the effect on the maximum frequency of the νCO stretching band. Plots correspond to pure DMPC (blue), and CPTtot/DMPC (mol/mol) 0.1 (red) and 0.25 (green). Spectra were collected from 10 to 34 °C, every 2 °C. Data correspond to the mean of three independent repetitions ± SD (error bars).
Figure 6. The effect of CPT on the region of the DMPC acyl chains and polar headgroups, determined by FTIR. Top: the effect on the maximum frequency of the νCH2 symmetric stretching band. Bottom: the effect on the maximum frequency of the νCO stretching band. Plots correspond to pure DMPC (blue), and CPTtot/DMPC (mol/mol) 0.1 (red) and 0.25 (green). Spectra were collected from 10 to 34 °C, every 2 °C. Data correspond to the mean of three independent repetitions ± SD (error bars).
Membranes 14 00118 g006
Figure 7. Effect of CPT on the SAXS profiles of DMPC. Diffractograms correspond to pure DMPC (red) and CPTtot/DMPC 0.1 (mol/mol) (blue), at 7 °C (A), 18 °C (B), and 34 °C (C).
Figure 7. Effect of CPT on the SAXS profiles of DMPC. Diffractograms correspond to pure DMPC (red) and CPTtot/DMPC 0.1 (mol/mol) (blue), at 7 °C (A), 18 °C (B), and 34 °C (C).
Membranes 14 00118 g007
Figure 8. Area per lipid vs. simulated time for all simulated systems: POPC at 298 K (black), POPC + CPT at 298 K (red), DMPC at 312 K (blue), and DMPC + CPT at 312 K (green).
Figure 8. Area per lipid vs. simulated time for all simulated systems: POPC at 298 K (black), POPC + CPT at 298 K (red), DMPC at 312 K (blue), and DMPC + CPT at 312 K (green).
Membranes 14 00118 g008
Figure 9. Mass density profiles along the z-axis of the simulation box of the simulated system of POPC + CPT. POPC phosphorus atoms are in green, CPT in red, POPC terminal methyl of sn-1 chains in blue, and POPC carbonyl groups in black. Curves are symmetrized around the center of the bilayer.
Figure 9. Mass density profiles along the z-axis of the simulation box of the simulated system of POPC + CPT. POPC phosphorus atoms are in green, CPT in red, POPC terminal methyl of sn-1 chains in blue, and POPC carbonyl groups in black. Curves are symmetrized around the center of the bilayer.
Membranes 14 00118 g009
Figure 10. Final snapshot of the simulation box of POPC + CPT membranes. Water molecules are shown in red lines, CPT in orange sticks, POPC atoms in green lines, POPC carbonyl groups in red spheres, POPC methyl terminals of sn-1 chains in blue spheres, and phosphorous atoms in yellow spheres.
Figure 10. Final snapshot of the simulation box of POPC + CPT membranes. Water molecules are shown in red lines, CPT in orange sticks, POPC atoms in green lines, POPC carbonyl groups in red spheres, POPC methyl terminals of sn-1 chains in blue spheres, and phosphorous atoms in yellow spheres.
Membranes 14 00118 g010
Figure 11. Free energy profiles for a water molecule crossing a POPC membrane (black line) and a POPC + CPT membrane (red line). Profiles are assumed to be symmetric across the bilayer center.
Figure 11. Free energy profiles for a water molecule crossing a POPC membrane (black line) and a POPC + CPT membrane (red line). Profiles are assumed to be symmetric across the bilayer center.
Membranes 14 00118 g011
Table 1. The effect of CPT on the maximum frequency of the FTIR νCH2 symmetric stretching and νCO stretching bands of various phosphatidylcholines. The CPTtot/phospholipid molar ratio was 0.1 in all cases. Data correspond to the mean of three independent repetitions ± SD.
Table 1. The effect of CPT on the maximum frequency of the FTIR νCH2 symmetric stretching and νCO stretching bands of various phosphatidylcholines. The CPTtot/phospholipid molar ratio was 0.1 in all cases. Data correspond to the mean of three independent repetitions ± SD.
ConditionsCH2 Symmetric StretchingC=O Stretching
DMPC 12 °C2848.3 ± 0.121730.9 ± 0.06
DMPC + CPT 12 °C2849.0 ± 0.051733.5 ± 0.09
DMPC 34 °C2850.1 ± 0.081729.8 ± 0.12
DMPC + CPT 34 °C2851.0 ± 0.041731.6 ± 0.06
DPPC 24 °C2848.6 ± 0.111731.4 ± 0.04
DPPC + CPT 24 °C2849.1 ± 0.031733.4 ± 0.03
DPPC 50 °C2851.1 ± 0.031730.7 ± 0.13
DPPC + CPT 50 °C2851.5 ± 0.031732.4 ± 0.10
POPC 25 °C2848.9 ± 0.081726.9 ± 0.05
POPC + CPT 25 °C2850.3 ± 0.061728.7 ± 0.07
Table 2. Parameters obtained from MD simulations for pure POPC and DMPC membranes, and after incorporation of CPT.
Table 2. Parameters obtained from MD simulations for pure POPC and DMPC membranes, and after incorporation of CPT.
POPCPOPC + CPTDMPCDMPC + CPT
Area per lipid (nm2)0.63 ± 0.010.66 ± 0.020.61 ± 0.010.64 ± 0.01
Order parameter (SCD)0.14 ± 0.070.15 ± 0.070.20 ± 0.050.193 ± 0.05
Thickness (nm)4.04 ± 0.073.93 ± 0.083.81 ± 0.063.65 ± 0.07
Diffusion coefficient (108 cm2 s−1)5.45 ± 2.197.67 ± 2.167.3 ± 1.0311.5 ± 3.2
Cluster (% monomers)-66.6 ± 0.9-65.6 ± 0.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ortiz, J.; Aranda, F.J.; Teruel, J.A.; Ortiz, A. Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study. Membranes 2024, 14, 118. https://doi.org/10.3390/membranes14060118

AMA Style

Ortiz J, Aranda FJ, Teruel JA, Ortiz A. Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study. Membranes. 2024; 14(6):118. https://doi.org/10.3390/membranes14060118

Chicago/Turabian Style

Ortiz, Julia, Francisco J. Aranda, José A. Teruel, and Antonio Ortiz. 2024. "Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study" Membranes 14, no. 6: 118. https://doi.org/10.3390/membranes14060118

APA Style

Ortiz, J., Aranda, F. J., Teruel, J. A., & Ortiz, A. (2024). Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study. Membranes, 14(6), 118. https://doi.org/10.3390/membranes14060118

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop