Next Article in Journal
Al2O3-Based Hollow Fiber Membranes Functionalized by Nitrogen-Doped Titanium Dioxide for Photocatalytic Degradation of Ammonia Gas
Next Article in Special Issue
Self-Assembly of Lipid Mixtures in Solutions: Structures, Dynamics Processes and Mechanical Properties
Previous Article in Journal
Membrane Fouling Prediction Based on Tent-SSA-BP
Previous Article in Special Issue
Notes on the Treatment of Charged Particles for Studying Cyclotide/Membrane Interactions with Dissipative Particle Dynamics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Membranes
Volume 12
Issue 7
10.3390/membranes12070692
membranes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes

by
José Villalaín
Institute of Research, Development, and Innovation in Healthcare Biotechnology (IDiBE), Universidad Miguel Hernández, E-03202 Elche, Spain
Membranes 2022, 12(7), 692; https://doi.org/10.3390/membranes12070692
Submission received: 21 June 2022 / Revised: 29 June 2022 / Accepted: 4 July 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Molecular Simulations of Biomembranes: From Biophysics Fundamentals to Biological Function)
Browse Figures
Versions Notes

Abstract

:
Procyanidins are known for their many benefits to human health and show a plethora of biological effects. One of the most important procyanidin is the procyanidin trimer C1 (PC1). Due to its relatively high lipid–water partition coefficient, the properties of PC1 could be attributed to its capability to interact with the biomembrane, to modulate its structure and dynamics, and to interact with lipids and proteins, however, its biological mechanism is not known. We have used all-atom molecular dynamics in order to determine the position of PC1 in complex membranes and the presence of its specific interactions with membrane lipids, having simulated a membrane mimicking the plasma membrane and another mimicking the mitochondrial membrane. PC1 has a tendency to be located at the membrane interphase, with part of the molecule exposed to the water solvent and part of it reaching the first carbons of the hydrocarbon chains. It has no preferred orientation, and it completely excludes the CHOL molecule. Remarkably, PC1 has a tendency to spontaneously aggregate, forming high-order oligomers. These data suggest that its bioactive properties could be attributed to its membranotropic effects, which therefore supports the development of these molecules as therapeutic molecules, which would open new opportunities for future medical advances.

1. Introduction

Procyanidins are a group of bioactive molecules that are known for their benefits to human health. The use of these molecules is increasing in the treatment of different chronic diseases, such as cancer, diabetes, and cardiovascular diseases, since they supposedly prevent cell damage that is possibly related to oxidative stress. Procyanidins are included within the condensed tannins group and, together with the hydrolysable tannins groups, form the tannin group of heterogeneous phytochemicals with a high molecular weight [1,2]. Tannins, in general, form a defense strategy in plants against both biotic and abiotic stresses [1]. Procyanidins belong to the polyphenols group, which comprise one of the most dominant classes of antioxidant bioactive compounds in nature [3]. Interestingly, the quantity of procyanidins is very high in vegetables, fruits, cereals, and the seeds of numerous species of plants [2,4] and they are very heterogeneous compounds with a diverse chemical structure [5,6]. Most importantly, they have been shown to possess a plethora of valuable properties, which include anticarcinogenic, antimicrobial, antiviral, antiinflammatory, antiallergic, antimutagenic, and antihyperglycaemic effects [3,7,8,9,10,11,12,13,14,15].
The procyanidin trimer C1 (PC1) is one of the most important procyanidins with many biological effects (Figure 1A). For example, PC1 has antiproliferative properties as it can induce apoptosis in breast cancer cells [16]; in addition, it also prevents melanoma cell growth [17]. PC1 has also been shown to inhibit the Na+,K+-ATPase [18] and the H+,K+-ATPase [19] and it has been demonstrated that PC1 also promotes preadipocyte differentiation and enhances insulin sensitivity [20,21]. Moreover, PC1 has senotherapeutic activity and increases the health and the life span of mice, possibly through the impairment of the functional integrity of the mitochondria [22]. Due to its relatively high lipid–water partition coefficient, the possible effects of PC1 on the biological systems could be attributed to its capability to locate the biological membrane and to modulate the membrane structure and dynamics. The interaction of flavonoids in general (procyanidins in particular) with the membrane components, in both proteins and lipids, might be responsible, at least in part, for their effects on human and animal health [23]. The XLogP3 value [24] of PC1 is 3.3, but its hydrophobic properties do not give us any indication about its location, orientation, or its possible interactions with the membrane. Significantly, it has been shown that procyanidins can interact with membrane lipids, which alters the biophysical properties of biomembranes [25,26]. Specifically, hexameric procyanidins can interact with membrane lipid rafts [27], which could imply that other procyanidins (shorter and larger oligomers) could also interact with lipid rafts, either sphingomyelin, cholesterol, or both, as well as with other types of membrane lipids. PC1 can distress the biological membrane, interact with the cellular proteins, or both, however, due to its high number of beneficial and different biological effects, a common mechanism should be at play. It could, therefore, be thought that its mechanism of action could be ascribed to its capability to interact with the biological membrane, modulate its fluidity, morphology, and permeability, and, finally, to interact with its components. The knowledge of the PC1 mechanism of action might open up new avenues for future therapeutic developments of this molecule or other related molecules. It is important to note that, under gastric conditions, procyanidins dimers and trimers do not transform into their monomer constituents, and they are found intact in the plasma after ingestion [28,29]. Moreover, high-concentration and high-frequency treatment of procyanidin C1 has no systemic toxicities [22,30].
Molecular dynamics (MD) is perfectly suitable for learning about the dynamics, the location, the interaction, and the structure of many different types of molecules interacting with biological model membranes [31,32]. We have used atomistic MD to define the location and the orientation of PC1 in the membrane and, at the same time, to determine the presence of any interactions with the membrane lipids. It is difficult to attempt to simulate a system as complex as that of a biological membrane. Furthermore, one has to consider the existence of different biological membranes, due to their different lipid and protein composition. Since the biological effects of PC1 could be related to the cellular oxidative stress, we have studied two different model biomembrane systems; one of them was similar to the lipid composition of the plasma membrane (PM) and the other was similar to the mitochondrial membrane (MIT). We believe that this approximation is a valid one and, therefore, the results obtained could be extrapolated to the situation that occurs in a biological cell; that is, to explain the possible mechanism of action of the PC1 molecule at the membrane level. For that goal, we have studied five different systems, three PM-derived and two MIT-derived, containing different numbers and locations of PC1 molecules (Table 1). Our results suggest that PC1 tends to be primarily placed at the membrane interface, between the phosphate atoms of the phospholipids and the oxygen atom of CHOL, which have a high propensity to aggregate through hydrogen bonding with each other which therefore affects the biophysical properties of the membrane lipids. Our data suggest that the bioactive properties of PC1 could be attributed to its membranotropic effects and accordingly through the modulation of the biophysical membrane properties.

2. Materials and Methods

Unrestrained all-atom MD simulations were carried out using NAMD 2.14 [33] with the CHARMM36 protein and lipid force fields [34,35,36]. All MD parameters used in this work have been described previously [37,38]. The whole systems, comprising water, ions, membrane, and protein/s, were equilibrated before the simulation for 5 ns after 100,000 steps of minimization, so as to remove unfavorable atomic contacts. The production trajectory for each one of the five protein/membrane systems studied in this work was calculated for a total of 675 ns (Table 1).
We have studied two different model biomembrane systems, one of them was similar to the lipid composition of the plasma membrane (PM) and the other was similar to the mitochondrial membrane (MIT). We have studied five systems, three PM-derived and two MIT-derived, each one containing a different number of PC1 molecules and locations (Table 1). The membrane systems were assembled using the Charmm-Gui web server [39]. The systems contained excess water [40]. Each one of the systems was composed of PC1 molecules, a membrane bilayer, water, and NaCl at physiological conditions, i.e., a concentration of 0.15 M, enclosed in a rectangular box, and in a neutral setting (Table 1) [41,42,43]. The PM system contained 80 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 48 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), 18 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine (POPS), 16 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoinositol-3-phosphate (PI-3P), 34 molecules of N-stearoyl-D-erythro-sphingosylphosphorylcholine (PSM), and 84 molecules of cholesterol (CHOL). The MIT system contained 90 molecules of POPC, 68 molecules of POPE, 2 molecules of POPS, 10 molecules of PI-3P, 2 molecules of PSM, 2 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphate (POPA), 28 molecules of 1′,3′-bis [1-palmitoyl-2-oleoyl-sn-glycero-3-phospho]-glycerol (cardiolipin, CL), and 20 molecules of CHOL [44,45,46] (Table 1). The chemical structures of the lipid molecules used in this study are shown in Supplementary Figure S1. The presence of one oleoyl chain in the phospholipids increase its fluidity in the membrane, makes the membrane more similar to the real one, and reduces the simulation time. The total number of lipids was 280, 140 per leaflet, for the PM system and 222, 111 per leaflet for the MIT system. The normal bilayer was parallel to the z-axis of the membrane and its surface constituted the x–y plane. The height of the simulation box and the cross-sectional area were permitted to fluctuate independently of each other with no constraints. The molecule of PC1 was created and minimized using Discovery Studio 4.0 (Accelrys Inc., San Diego, USA). The CHARMM General Force Field (CGenFF) compatible stream files of PC1 were obtained using the Charmm-Gui web server [39]. The initial arrangements at t = 0 ns are depicted in Figure 1B–F for all five different systems. Systems 1 to 3 contained the PM membrane and systems 4 and 5 contained the MIT membrane. System 1 contained 16 molecules of PC1 at the external part of the membrane, with 8 at each side of the membrane (Figure 1B). System 2 contained 8 molecules of PC1 at the external part of the membrane, with 4 at each side of the membrane (Figure 1C). System 3 contained 4 molecules of PC1 at the middle of the palisade part of the membrane (Figure 1D). System 4 contained 8 molecules of PC1 at the external part of the membrane, with 4 at each side of the membrane (Figure 1E). System 5 contained 16 molecules of PC1 at the external part of the membrane, with 8 at each side of the membrane (Figure 1F). The distance between the PC1 molecules in water and the membrane surface (represented by the phosphate atoms of the phospholipids) was 16 Å.
The average layer of the lipid head-group phosphate atoms defined the membrane surface and were parallel to the x–y plane. VMD software was used for analysis and visualization, as before [37,38,47,48]. The center-of-masses and number of molecular contacts were obtained using standard VMD plugins. SCD order parameters, surface area per lipid, membrane thickness, molecular areas, and molecule tilt were calculated as previously described [49] using VMD “Membplugin” [48]. Mass density profiles were obtained by means of the VMD “Density Profile Tool” plugin [50]. The presence of hydrogen bonds were defined by a distance of less than 3 Å between the acceptor and donor atoms and an acceptor-H-donor angle of at least 150° [51]. The complete simulations were used for analysis, i.e., 675 ns, unless otherwise stated.

3. Results

We have used five different membrane/PC1 systems, each one of them divergent from the others by the number and the location of the PC1 molecules in the system, as well as in the model biomembrane used (Table 1). System one contained the PM membrane and 16 molecules of PC1 at the external part (Figure 1B), system two contained the PM membrane and eight molecules of PC1 at the external part (Figure 1C), system three contained the PM membrane and four molecules of PC1 at the middle part of the bilayer (Figure 1D), system four contained the MIT membrane and eight molecules of PC1 at the external part of the membrane, with four at each side of the membrane (Figure 1E), and system five contained the MIT membrane and 16 molecules of PC1 at the external part of the membrane (Figure 1F). The time variation of the lipid area, as well the time variation of the membrane thickness, were used in order to assess the membrane system equilibration both during and at the end of the simulation [52,53]. The lipids areas remained constant after ~35 ns (Supplementary Figure S2). The mean area of all of the lipids in the different systems for the last 30 ns of simulation are shown in Supplementary Table S1. The obtained data show that average molecular areas were comparable to those that have been previously reported [49,54,55]. The membrane thickness also remained practically constant after ~35 ns (Supplementary Figure S3A). The average membrane thickness for all of the systems for the last 30 ns was 43–44 Å (rounded to the first integer, Supplementary Table S2 and Supplementary Figure S3B), which is comparable to those that have been described for systems containing CHOL [54]. This data indicate that the systems were equilibrated early on in the progress of the simulation and have revealed that the membrane systems attained a steady state after ~35 ns of the simulation.
At first, the PC1 molecules were placed in both the upper and lower water layers and, in the case of system three, in the middle of the membrane (see Figure 1). The final, t = 675 ns, snapshots for all of the systems are shown in Figure 1B–F. At the end of the simulation, the PC1 molecules, which at the beginning of the simulation were located in the middle of the upper and lower water layers, had moved to a position near the membrane interface (systems one and two for the PM membrane system, Figure 1B,C, and systems four and five for the MIT membrane system, Figure 1E,F). In a similar way, the PC1 molecules, which at the beginning of the simulation were inserted into the middle of the membrane, moved to a place near to the membrane interface, except for one system (system three, Figure 1D). However, interestingly, many of the PC1 molecules, depending on the system, aggregated with time. Not only did they form dimers, trimers, or tetramers, but also higher-order oligomers (Table 2). In system one, which had 16 PC1 molecules, only two of them were in the monomer state, four of them were in the dimer state, and ten of the molecules formed an oligomer (Table 2). This is in contrast to system two, which had half of the concentration of system one, i.e., eight PC1 molecules, which had five PC1 molecules in the monomer state and only three forming a trimer (Table 2). System four, which had eight PC1 molecules, had six PC1 molecules in the monomer state and two forming a dimer (Table 2). On the other hand, system five, which had 16 PC1 molecules, had four monomers, one tetramer, and an oligomer that was formed by eight molecules of PC1 (Table 2). System three was different from the others in the sense that it only had four molecules of PC1 inserted into the middle part of the membrane. In this case, all of the molecules were in the monomer state at the termination of the MD simulation (Table 2). From these data, it is clear that PC1 aggregates in water in a concentration-dependent mode and it is independent of the model membrane system. The total number of monomers, taking into account all of the five systems, was 21 and the number of PC1 molecules in aggregated form was 31 (Table 2). We have obtained the number of contacts and hydrogen bonds between the pairs of PC1 molecules for the last 30 ns of the MS simulation, either forming dimers or dimers in higher-number oligomers. The average number of contacts between each pair was 15.7 ± 5.9 and the average number of hydrogen bonds was 0.96 ± 0.66. This data highlights the predisposition of PC1 molecules to aggregate both in solution and at the membrane surface.
In order to check the conduct of all of the PC1 molecules in the membrane throughout the whole simulation, we have obtained the time variation of their center-of-mass (COM) compared with the COMs of the phosphate atoms in both leaflets, which make up the membrane surface (Figure 2). For system one, eight molecules of PC1 were located at the middle of the upper water layer and eight molecules of PC1 were located at the middle of the lower water layer (Figure 1B). At the beginning of the MD simulation, some of them displayed a significant fluctuation in their COM, whereas others did not (Figure 2A). Interestingly, and at around 400–450 ns, several of the freely moving PC1 molecules in the water layer dramatically reduced their motion, coinciding with the time when those PC1 molecules formed a massive aggregate (Figure 2A). The histograms corresponding to the COM of all of the PC1 molecules in system one and for the last 30 ns of the simulation are shown in Supplementary Figure S4, where it can be observed that all of them were close to the membrane surface and some had even crossed it. Interestingly, of the sixteen different isolated molecules at the beginning of the simulation, only two remained in the monomer state (Table 2). In the case of system two, four molecules of PC1 were located at the middle of the upper water layer and four molecules of PC1 were located at the middle of the lower water layer (Figure 1C). As observed in Figure 2B, and at the beginning of the MD simulation, some of them displayed a significant oscillation in their COM. However, at an early stage of less than 200 ns, they moved to a position near to the membrane surface and their erratic movement stopped (Figure 2B). In this case, only one trimer was encountered, the other PC1 molecules remained in the monomer state throughout the MD simulation (Table 2). The histograms corresponding to the COM of all of the PC1 molecules and for the last 30 ns of the simulation is shown in Supplementary Figure S5, where it can be observed that all of them were close to the membrane surface and some had even crossed it. System three was very different to the previous systems, since four molecules of PC1 were located at the middle of the membrane bilayer (Figure 1D). At the beginning of the MD simulation, all of them displayed a significant fluctuation in their COM, but after about 100 ns of the simulation time they did not move significantly from the position until the end (Figure 2C). All four of the molecules remained in the monomer state for all of the MD simulation (Table 2), three of them located to a position near to the membrane surface, but one PC1 molecule remained at the hydrocarbon region of the membrane. This tendency can be observed in the histograms corresponding to the COMs of the last 30 ns of the simulation (Supplementary Figure S6). Systems one to three comprised the PM model membrane. However, systems four and five comprised the MIT model membrane. System four had eight molecules of PC1, of which four were in the middle of the upper water layer, and four were in the middle of the lower water layer (Figure 1E). At the beginning of the MD simulation, and similar to the other systems, many of them displayed a significant fluctuation in their COM, whereas other did not (Figure 2C). However, after about 300 ns, the movement on the z-axis of all of the PC1 molecules stopped since all of them were located near to the surface of the membrane (Figure 2C). For system four, only one dimer was formed, the other six molecules remained in the monomer state (Table 2). The histograms for the last 30 ns of the MD simulation show the location of the PC1 molecules in this system (Supplementary Figure S7). System five had 16 molecules of PC1 in total, of which eight molecules were located at the middle of the upper water layer and eight molecules were located at the middle of the lower water layer (Figure 1F). Similar to the other systems, and at the beginning, all of them had an erratic system at the water layer; however, at about 150 ns they were approaching the membrane surface and remained there, their movement restricted to the layer near to it (Figure 2E). For this system, only four molecules remained in the monomer state, while the others formed a tetramer and an octamer (Table 2). Supplementary Figure S8 shows the histograms of the COMs for the last 30 ns of the simulation. It is possible to observe that all of the PC1 molecules, either forming monomers or oligomers, were located near to the membrane surface.
As observed in Figure 2 and Supplementary Figures S4–S8, nearly all of the PC1 molecules at the conclusion of the simulation time were located near to or at the membrane surface, either with the PM model membrane (systems one to three) or the MIT model membrane (systems four and five). It is very interesting to compare the location of all of the PC1 molecules in all of the systems studied, as shown in Figure 3, where the mean COMs for the last 30 ns of the simulation are shown. All of them, even the aggregated ones, are all less than 10 Å from the surface, except for one PC1 molecule in system three, which is near to the center of the bilayer (Figure 3). If we take only the location data of the PC1 molecules in the monomer state, the mean COM for the positive layer was 19.5 ± 2.8 Å, whereas the mean COM for the negative layer was −19.7 ± 5.8 Å (Figure 3F), i.e., an absolute global mean COM of 19.6 Å. Taking into account that the positive membrane surface was 22.1 ± 0.5 Å and the negative membrane surface was −22 ± 0.6 Å, the preferred location of the PC1 molecules is one in which the COM is about 2.5 Å below the membrane surface (Figure 3). However, we are talking about the COM of the molecule, but we have to take into account the real space that the PC1 molecule occupies. The molecule is, approximately, a rectangle of dimensions of about 15 Å·12 Å·12 Å, so that the molecule tends to be positioned in the membrane, with part of it extending 4–5 Å outside and 8–10 Å inside of it. The average mass density of all of the components for the last 30 ns of the biomembrane model systems are presented in Supplementary Figure S9. All of the profiles are practically symmetric between the two leaflets of the membrane, which implies a comparable behavior for all of the lipids inside of them. By comparing the mass density profiles between both the PM and the MIT systems (Supplementary Figure S9A–E,G–I, respectively), it is interesting to point out that the profiles are very similar, indicating a similar behavior of all of the lipids in both of the systems. As observed in the Figure 3, the global mass density of all of the PC1 molecules extends to the external part of the membrane, as well as to the internal part, passing the phosphate atom layer but never extending beyond the molecules of cholesterol. If we observe the mass density profile of each individual PC1 molecule compared to the phosphate mass density, this pattern looks much better in the sense that their location is more defined. Comparing the location of the PC1 molecules in the monomeric state with those that form aggregates, it can be clearly seen that all of them, whether in the monomeric or the aggregated form, have a tendency to integrate at the membrane interface.
We have measured the number and type of lipid molecules that are located within two Å of the PC1 molecules, i.e., in close contact, and the results are shown in Figure 4. In the case of the PM-containing systems (systems one, two, and three) the most significant result is that CHOL is excluded from proximity to the PC1 molecules, since its global percentage composition is 30%, but reduces to about 18% when in contact with the PC1 molecules. The same can be said about POPS, since its global percentage is about 6.5%, which reduces to about 1.5% in the PM systems. On the contrary, PI-3P slightly increases, since its global percentage is 5.5% and in the surrounding area of the PC1 molecules it is about 9%. As for the other lipids, there is no clear preference for them to be near to the PC1 molecules, since the percentage of molecules near to the PC1 molecules is relatively similar to the global percentage in the system (Figure 4). For the MIT-containing systems, and similar to the PM-containing systems, the CHOL is completely excluded in the surrounding area of the PC1 molecules, since its global percentage composition is 9.5% but reduces to 0% in the surrounding area of the PC1 molecules (Figure 4). Interestingly, the quantity of the POPC molecules surrounding the PC1 molecules increased significantly, since its global percentage is 40% but increased to 51% around the PC1 molecules. POPA increased slightly, since its global percentage is about 0.5%, which increased to about 3% in the surrounding area of PC1 (Figure 4). For CL, the characteristic lipid of the mitochondrial membrane, there is no preference/aversion to being around the PC1 molecules, since there was no difference in the global and the local percentages. As for the other lipids, there is no clear preference for them to be near to the PC1 molecules, since the percentage of molecules near to PC1 molecules is relatively similar to the global percentage in the system (Figure 4). For all of the systems studied, the most significant result is that CHOL is the molecule that tends to be excluded from the surrounding area of the PC1 molecules.
It is known that molecules interacting with the membrane can affect the hydrocarbon chain order of the phospholipid acyl chains. Therefore, we have studied the effect of the PC1 molecules on the hydrocarbon chain order of the phospholipid acyl chains by obtaining their deuterium order parameter, SCD (Supplementary Figures S10 and S11 for the PM and MIT systems, respectively). When there is complete disorder, the SCD value is 0, and a value of 0.5 specifies full order along the normal bilayer [56]. The average –SCD values of the acyl chains of all of the phospholipids in the PM model system, i.e., POPC, POPE, POPS, PI-3P, and PSM, were in agreement with the profiles that were observed earlier for the experimental and simulated data [36,57,58] (Supplementary Figure S10). However, for some of the phospholipids that were near to the PC1 molecules, there were significant, although not dramatic, changes in the SCD profiles. For POPC, POPE, and PI-3P, a general decrease in the SCD values was observed, whereas an increase was observed for POPS (Supplementary Figure S10). In contrast, no effect was observed for PSM. The decrease in SCD that was observed for the phospholipids POPC, POPE, and PI-3P indicates that the PC1 molecules increase the fluidity of the hydrocarbon chains of these phospholipids, whereas the PC1 molecules increase the rigidity of the hydrocarbon chains of POPS. Similar to the PM system discussed above, the average –SCD values of the acyl chains of all of the phospholipids in the MIT model system, i.e., POPC, POPE, POPS, POPA, PI-3P, CL, and PSM, were all in agreement with the profiles for the previous experimental and simulated data [36,57,58] (Supplementary Figure S11). However, for those phospholipids near to the PC1 molecules, there were significant changes in the SCD profiles. In this case, and for all of the phospholipids, the presence of the PC1 molecules decreased the SCD values, indicating an increase in the fluidity of the hydrocarbon chains. It can be inferred from these data that the PC1 molecules insert relatively well in between the hydrocarbon chains of the phospholipids. They do not show a dramatic effect on the anisotropy of the hydrocarbon chains, but the general trend is that the PC1 molecules increase the fluidity of the membrane, both in the PM and in the MIT systems.

4. Discussion

The bioactive procyanidins are known for their ample benefits to human health and have been shown to have a plethora of biological effects [3,7,8,9,10,11,12,13,14,15]. The bioactive effects of procyanidins in general, PC1 in particular, could be ascribed to their capability to localize in the biological membrane and to modulate the membrane structure and/or the dynamics, as well as their interaction with the lipids and proteins [23]. The numerous biological properties of these molecules, including procyanidins, demonstrate that diverse mechanisms are at play, but also show the presence of a common point of action. That common point is that the biological membrane and the bioactive properties of these molecules could be associated to a membrane modulation mechanism [59,60]. However, the biological mechanism of the procyanidins has been not identified, despite the knowledge that it could help in the development of new therapeutic molecules. Their biological, pharmacological, and medicinal properties should be linked to this membrane’s modulation effect, which itself rests on their structure, their interactions, their location, and their orientation in the membrane [59,60]. Therefore, the biological properties of PC1 could therefore be ascribed, not only to its possible interaction with proteins, but also to its ability to modulate the membrane’s biophysical properties.
We have aimed to locate the molecule of PC1 and to determine the possible specific interactions of this molecule with the membrane lipids using MD simulations. For that objective, we have used two different model membranes, one mimicking the plasma membrane and the other mimicking the mitochondrial membrane. Our data reveal that PC1 tends to be located at the membrane interphase, with part of the molecule exposed to the water solvent and part of it reaching the first carbons of the hydrocarbon chains, but never extending beyond that so that they do not reach the middle of the bilayer leaflet (Figure 5, compare with Figure 3F). Furthermore, due to the specific three-dimensional structure of the PC1 molecule, it has no preferred orientation with all of its hydroxyl groups in the outer part of the molecule. One of the most interesting facts that we have found in this study is that the PC1 molecules completely exclude the molecule of CHOL, both in the PM and in the MIT membrane systems. In the case of the PM system, there were differences with the other lipids present, such as the reduction in POPS and an increase in PI-3P; however, no significant differences were found for the other lipids in the system. In the case of the MIT system, both POPC and POPA seem to increase in the surrounding area of the PC1 molecules. Similar to the PM system, no differences were found for the other lipids in the MIT system. Yet, one of the most remarkable facts about PC1 is that, in the solution, it has a tendency to combine forming dimers, trimers, and higher-order aggregates between the different PC1 molecules. These groups of PC1 molecules form spontaneously through the formation of hydrogen bonds; the formation of the hydrogen bonds being independent of the systems studied here, be it either the PM or the MIT system. Interestingly, the formation of the aggregates did not prevent the PC1 molecules from interacting with the membrane, either in the PM or the MIT systems. PC1 is known to be stable under gastric conditions, it is not fragmented into its monomer constituents, it is found in the plasma after ingestion and, fundamentally, relatively high concentrations and high frequency treatments of PC1 have no systemic toxicities. However, the formation of the PC1 aggregates could hamper its bioactive properties and, consequently, this should be taken into account when planning its use in clinical trials in order to choose an appropriate vehicle for its preparation. Our work should help to advance these molecules as therapeutic molecules by opening up new avenues for future medical advances.

5. Conclusions

PC1 tends to be located at the membrane interphase, with a part of the molecule exposed to the solvent and part of it reaching the carbonyl region of the hydrocarbon chains. PC1 molecules, when inserted into the membrane, have no preferred orientation, having all of the hydroxyl groups in the outer part of the molecule. Significantly, PC1 completely excludes the molecule of CHOL in both the PM and the MIT systems, so that the existence of domains with and without cholesterol cannot be ruled out. Outstandingly, PC1 molecules have a tendency to spontaneously associate, forming aggregates, which does not prevent them from interacting with the membrane.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/membranes12070692/s1, Figure S1: Chemical structures of the lipid molecules used in this study.(A) POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine), (B) POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine), (C) POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine), (D) POPA (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphate), (E) PI-3P (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoinositol-3-phosphate), (F) CL (1′,3′-bis[1-palmitoyl-2-oleoyl-sn-glycero-3-phospho]-glycerol), (G) PSM (N-stearoyl-D-erythro-sphingosylphosphorylcholine) and (H) CHOL (cholesterol); Figure S2: Time variation of the molecular areas for the whole simulation time for each lipid species in the (A) system 1, (B) system 2, (C) system 3, (D) system 4 and (E) system 5 membrane model systems. The following lipids are depicted: POPC (black), POPE (blue), POPS (red), PI-3P (magenta), POPA (orange), PSM (olive), CL (violet) and CHOL (wine); Figure S3: Time variation of (A) membrane thickness for the whole simulation time and (B) average membrane thickness for the last 30 n of simulation for system 1 (black), system 2 (blue), system 3 (red), system 4 (magenta) and system 5 (olive); Figure S4: Average COM z-axis for the last 30 n of simulation for all the PC1 molecules in system 1. Blue and black histograms correspond to monomer and oligomer states, respectively. The average COM of the phosphate atoms of the phospholipids is depicted in red color; Figure S5: Average COM z-axis for the last 30 n of simulation for all the PC1 molecules in system 2. Blue and black histograms correspond to monomer and oligomer states, respectively. The average COM of the phosphate atoms of the phospholipids is depicted in red col-our; Figure S6: Average COM z-axis for the last 30 n of simulation for all the PC1 molecules in system 3. Blue and black histograms correspond to monomer and oligomer states, respectively. The average COM of the phosphate atoms of the phospholipids is depicted in red color; Figure S7: Average COM z-axis for the last 30 n of simulation for all the PC1 molecules in system 4. Blue and black histograms correspond to monomer and oligomer states, respectively. The average COM of the phosphate atoms of the phospholipids is depicted in red color; Figure S8: Average COM z-axis for the last 30 n of simulation for all the PC1 molecules in system 5. Blue and black histograms correspond to monomer and oligomer states, respectively. The average COM of the phosphate atoms of the phospholipids is depicted in red color; Figure S9: Mass density profiles for the last 30 ns of the MD simulation for (A, B) system 1, (C, D) system 2, (E, F) system 3, (G, H) system 4 and (I, J) system 5. The global mass density of all the lipids and all the PC1 molecules are shown in (A, C, E, G and I), whereas the mass density of the phosphate atoms of the phospholipids and each one of the individual PC1 molecules are shown in (B, D, F, H and J). The mass density profiles correspond to POPC (blue), POPE (orange), POPS (red), PI-3P (magenta), POPA (wine), PSM (olive), CL (navy), CHOL (violet) and PC1 (black). In (B, D, F, H and J) the PC1 molecules in the monomer state are depicted as dotted red lines; Figure S10: Average deuterium order parameter –SCD calculated for the hydrocarbon chains of the phospholipids in the PM system. (A, C, E, G) oleoyl and (B, D, F, H) palmitoyl acyl chains of (A, B) POPC, (C, D) POPE, (E, F) POPS and (G, H) PI-3P as well as the palmitoyl (I) and sphyngosyl (J) acyl chains of PSM. The data correspond to the bulk phospholipid acyl chains (-■-) and the phospholipid acyl chains within 5 Å of the PC1 molecules (- Membranes 12 00692 i003-). The analysis was carried out for the last 30 ns of simulation; Figure S11: Average deuterium order parameter –SCD calculated for the hydrocarbon chains of the phospholipids in the MIT system. (A, C, E, G, I) oleoyl and (B, D, F, H, J) palmitoyl acyl chains of (A, B) POPC, (C, D) POPE, (E, F) POPS, (G, H) PI-3P and (K, L) POPA, the palmitoyl (K) and vacenoyl (L) acyl chains of CL and the palmitoyl (M) and sphyngosyl (N) acyl chains of PSM. The data corre-spond to the bulk phospholipid acyl chains (-■-) and the phospholipid acyl chains within 5 Å of the PC1 molecules (- Membranes 12 00692 i003-). The analysis was carried out for the last 30 ns of simulation; Table S1: Average molecular area (Å2) and membrane thickness (Å) for the last 30 ns of simulation of all the lipids in the systems studied in this work.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

NAMD was developed by the Theoretical and Computational Biophysics Group in the Beckman Institute for Advanced Science and Technology at the University of Illinois at Urbana-Champaign. We are very grateful to SIPT, Universitas Miguel Hernández (UMH), for the generous use of the UMH Computer Cluster. This work was not funded by any external or internal funding agencies.

Conflicts of Interest

The author declares that there are no conflict of interest.

Abbreviations

CHOLCholesterol
CL(1′,3′-bis[1-palmitoyl-2-oleoyl-sn-glycero-3-phospho]-glycerol)
COMCenter of mass
MITMitochondrial membrane
PC1Procyanidin C1 [(-)-epicatechin–(4β→8)-(-)-epicatechin (4β→8)-(-)-epicatechin]
PI-3P1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoinositol-3-phosphate
PMPlasma membrane
POPA1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphate
POPC1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine
POPE1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine
POPS1-Palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine
PSMN-Palmitoyl-D-erythro-sphingosylphosphorylcholine

References

  1. Smeriglio, A.; Barreca, D.; Bellocco, E.; Trombetta, D. Proanthocyanidins and hydrolysable tannins: Occurrence, dietary intake and pharmacological effects. Br. J. Pharmacol. 2017, 174, 1244–1262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Valencia-Hernandez, L.J.; Wong-Paz, J.E.; Ascacio-Valdes, J.A.; Chavez-Gonzalez, M.L.; Contreras-Esquivel, J.C.; Aguilar, C.N. Procyanidins: From Agro-Industrial Waste to Food as Bioactive Molecules. Foods 2021, 10, 3152. [Google Scholar] [CrossRef] [PubMed]
  3. Kiokias, S.; Oreopoulou, V. A Review of the Health Protective Effects of Phenolic Acids against a Range of Severe Pathologic Conditions (Including Coronavirus-Based Infections). Molecules 2021, 26, 5405. [Google Scholar] [CrossRef] [PubMed]
  4. Aron, P.M.; Kennedy, J.A. Flavan-3-ols: Nature, occurrence and biological activity. Mol. Nutr. Food Res. 2008, 52, 79–104. [Google Scholar] [CrossRef]
  5. Dixon, R.A.; Xie, D.Y.; Sharma, S.B. Proanthocyanidins—A final frontier in flavonoid research? New Phytol. 2005, 165, 9–28. [Google Scholar] [CrossRef] [Green Version]
  6. Heim, K.E.; Tagliaferro, A.R.; Bobilya, D.J. Flavonoid antioxidants: Chemistry, metabolism and structure-activity relationships. J. Nutr. Biochem. 2002, 13, 572–584. [Google Scholar] [CrossRef]
  7. Dorenkott, M.R.; Griffin, L.E.; Goodrich, K.M.; Thompson-Witrick, K.A.; Fundaro, G.; Ye, L.; Stevens, J.R.; Ali, M.; O’Keefe, S.F.; Hulver, M.W.; et al. Oligomeric cocoa procyanidins possess enhanced bioactivity compared to monomeric and polymeric cocoa procyanidins for preventing the development of obesity, insulin resistance, and impaired glucose tolerance during high-fat feeding. J. Agric. Food Chem. 2014, 62, 2216–2227. [Google Scholar] [CrossRef] [Green Version]
  8. Mao, J.T.; Smoake, J.; Park, H.K.; Lu, Q.Y.; Xue, B. Grape Seed Procyanidin Extract Mediates Antineoplastic Effects against Lung Cancer via Modulations of Prostacyclin and 15-HETE Eicosanoid Pathways. Cancer Prev. Res. 2016, 9, 925–932. [Google Scholar] [CrossRef] [Green Version]
  9. David, I.M.B.; de Souza Fernandes, F.; Dos Santos Silva Ferreira, J.B.; Ludtke, D.D.; Martins, D.F.; Bobinski, F.; da Silva, T.; Buffon, L.D.; Kopper, M.B.R.; da Silva, G.S.; et al. Dietary supplementation with procyanidin-rich Pinus pinaster extract is associated with attenuated Ehrlich tumor development in mice. Nutr. Res. 2019, 62, 41–50. [Google Scholar] [CrossRef]
  10. Connell, B.J.; Chang, S.Y.; Prakash, E.; Yousfi, R.; Mohan, V.; Posch, W.; Wilflingseder, D.; Moog, C.; Kodama, E.N.; Clayette, P.; et al. A Cinnamon-Derived Procyanidin Compound Displays Anti-HIV-1 Activity by Blocking Heparan Sulfate- and Co-Receptor- Binding Sites on gp120 and Reverses T Cell Exhaustion via Impeding Tim-3 and PD-1 Upregulation. PLoS ONE 2016, 11, e0165386. [Google Scholar] [CrossRef]
  11. Song, X.; Yang, Y.; Li, J.; He, M.; Zou, Y.; Jia, R.; Li, L.; Hang, J.; Cui, M.; Bai, L.; et al. Tannins extract from Galla Chinensis can protect mice from infection by Enterotoxigenic Escherichia coli O101. BMC Complement. Med. Ther. 2021, 21, 84. [Google Scholar] [CrossRef] [PubMed]
  12. Tang, Y.; Xiong, R.; Wu, A.G.; Yu, C.L.; Zhao, Y.; Qiu, W.Q.; Wang, X.L.; Teng, J.F.; Liu, J.; Chen, H.X.; et al. Polyphenols Derived from Lychee Seed Suppress Abeta (1-42)-Induced Neuroinflammation. Int. J. Mol. Sci. 2018, 19, 2109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Huang, Q.; Liu, X.; Zhao, G.; Hu, T.; Wang, Y. Potential and challenges of tannins as an alternative to in-feed antibiotics for farm animal production. Anim. Nutr. 2018, 4, 137–150. [Google Scholar] [CrossRef] [PubMed]
  14. Kiokias, S.; Proestos, C.; Oreopoulou, V. Phenolic Acids of Plant Origin—A Review on Their Antioxidant Activity In Vitro (O/W Emulsion Systems) Along with Their in Vivo Health Biochemical Properties. Foods 2020, 9, 534. [Google Scholar] [CrossRef]
  15. Tsang, N.Y.; Li, W.F.; Varhegyi, E.; Rong, L.; Zhang, H.J. Ebola Entry Inhibitors Discovered from Maesa perlarius. Int. J. Mol. Sci. 2022, 23, 2620. [Google Scholar] [CrossRef]
  16. Koteswari, L.L.; Kumari, S.; Kumar, A.B.; Malla, R.R. A comparative anticancer study on procyanidin C1 against receptor positive and receptor negative breast cancer. Nat. Prod. Res. 2020, 34, 3267–3274. [Google Scholar] [CrossRef]
  17. Bae, J.; Kumazoe, M.; Murata, K.; Fujimura, Y.; Tachibana, H. Procyanidin C1 Inhibits Melanoma Cell Growth by Activating 67-kDa Laminin Receptor Signaling. Mol. Nutr. Food Res. 2020, 64, e1900986. [Google Scholar] [CrossRef]
  18. Heger, T.; Zatloukal, M.; Kubala, M.; Strnad, M.; Gruz, J. Procyanidin C1 from Viola odorata L. inhibits Na(+),K(+)-ATPase. Sci. Rep. 2022, 12, 7011. [Google Scholar] [CrossRef]
  19. Souccar, C.; Cysneiros, R.M.; Tanae, M.M.; Torres, L.M.; Lima-Landman, M.T.; Lapa, A.J. Inhibition of gastric acid secretion by a standardized aqueous extract of Cecropia glaziovii Sneth and underlying mechanism. Phytomedicine 2008, 15, 462–469. [Google Scholar] [CrossRef]
  20. Sun, P.; Li, K.; Wang, T.; Ji, J.; Wang, Y.; Chen, K.X.; Jia, Q.; Li, Y.M.; Wang, H.Y. Procyanidin C1, a Component of Cinnamon Extracts, Is a Potential Insulin Sensitizer That Targets Adipocytes. J. Agric. Food Chem. 2019, 67, 8839–8846. [Google Scholar] [CrossRef]
  21. Nakagawa, Y.; Ishimura, K.; Oya, S.; Kamino, M.; Fujii, Y.; Nanba, F.; Toda, T.; Ishii, T.; Adachi, T.; Suhara, Y.; et al. Comparison of the sympathetic stimulatory abilities of B-type procyanidins based on induction of uncoupling protein-1 in brown adipose tissue (BAT) and increased plasma catecholamine (CA) in mice. PLoS ONE 2018, 13, e0201203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Xu, Q.; Fu, Q.; Li, Z.; Liu, H.; Wang, Y.; Lin, X.; He, R.; Zhang, X.; Ju, Z.; Campisi, J.; et al. The flavonoid procyanidin C1 has senotherapeutic activity and increases lifespan in mice. Nat. Metab. 2021, 3, 1706–1726. [Google Scholar] [CrossRef] [PubMed]
  23. Verstraeten, S.V.; Fraga, C.G.; Oteiza, P.I. Interactions of flavan-3-ols and procyanidins with membranes: Mechanisms and the physiological relevance. Food Funct. 2015, 6, 32–41. [Google Scholar] [CrossRef] [PubMed]
  24. Cheng, T.; Zhao, Y.; Li, X.; Lin, F.; Xu, Y.; Zhang, X.; Li, Y.; Wang, R.; Lai, L. Computation of octanol-water partition coefficients by guiding an additive model with knowledge. J. Chem. Inf. Model. 2007, 47, 2140–2148. [Google Scholar] [CrossRef]
  25. Verstraeten, S.V.; Keen, C.L.; Schmitz, H.H.; Fraga, C.G.; Oteiza, P.I. Flavan-3-ols and procyanidins protect liposomes against lipid oxidation and disruption of the bilayer structure. Free Radic. Biol. Med. 2003, 34, 84–92. [Google Scholar] [CrossRef]
  26. Verstraeten, S.V.; Hammerstone, J.F.; Keen, C.L.; Fraga, C.G.; Oteiza, P.I. Antioxidant and membrane effects of procyanidin dimers and trimers isolated from peanut and cocoa. J. Agric. Food Chem. 2005, 53, 5041–5048. [Google Scholar] [CrossRef]
  27. Verstraeten, S.V.; Jaggers, G.K.; Fraga, C.G.; Oteiza, P.I. Procyanidins can interact with Caco-2 cell membrane lipid rafts: Involvement of cholesterol. Biochim. Biophys. Acta 2013, 1828, 2646–2653. [Google Scholar] [CrossRef] [Green Version]
  28. Rios, L.Y.; Bennett, R.N.; Lazarus, S.A.; Remesy, C.; Scalbert, A.; Williamson, G. Cocoa procyanidins are stable during gastric transit in humans. Am. J. Clin. Nutr. 2002, 76, 1106–1110. [Google Scholar] [CrossRef] [Green Version]
  29. Serra, A.; Macia, A.; Rubio, L.; Angles, N.; Ortega, N.; Morello, J.R.; Romero, M.P.; Motilva, M.J. Distribution of procyanidins and their metabolites in rat plasma and tissues in relation to ingestion of procyanidin-enriched or procyanidin-rich cocoa creams. Eur. J. Nutr. 2013, 52, 1029–1038. [Google Scholar] [CrossRef]
  30. Sun, P.; Wang, T.; Chen, L.; Yu, B.W.; Jia, Q.; Chen, K.X.; Fan, H.M.; Li, Y.M.; Wang, H.Y. Trimer procyanidin oligomers contribute to the protective effects of cinnamon extracts on pancreatic beta-cells in vitro. Acta Pharmacol. Sin. 2016, 37, 1083–1090. [Google Scholar] [CrossRef] [Green Version]
  31. Ingolfsson, H.I.; Melo, M.N.; van Eerden, F.J.; Arnarez, C.; Lopez, C.A.; Wassenaar, T.A.; Periole, X.; de Vries, A.H.; Tieleman, D.P.; Marrink, S.J. Lipid organization of the plasma membrane. J. Am. Chem. Soc. 2014, 136, 14554–14559. [Google Scholar] [CrossRef] [PubMed]
  32. Perricone, U.; Gulotta, M.R.; Lombino, J.; Parrino, B.; Cascioferro, S.; Diana, P.; Cirrincione, G.; Padova, A. An overview of recent molecular dynamics applications as medicinal chemistry tools for the undruggable site challenge. Medchemcomm 2018, 9, 920–936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Phillips, J.C.; Hardy, D.J.; Maia, J.D.C.; Stone, J.E.; Ribeiro, J.V.; Bernardi, R.C.; Buch, R.; Fiorin, G.; Henin, J.; Jiang, W.; et al. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 2020, 153, 044130. [Google Scholar] [CrossRef] [PubMed]
  34. Best, R.B.; Zhu, X.; Shim, J.; Lopes, P.E.; Mittal, J.; Feig, M.; Mackerell, A.D., Jr. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and chi(2) dihedral angles. J. Chem. Theory Comput. 2012, 8, 3257–3273. [Google Scholar] [CrossRef] [Green Version]
  35. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar] [CrossRef] [Green Version]
  36. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.W.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D., Jr.; Pastor, R.W. Update of the CHARMM all-atom additive force field for lipids: Validation on six lipid types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef] [Green Version]
  37. Villalain, J. Envelope E protein of dengue virus and phospholipid binding to the late endosomal membrane. Biochim. Biophys. Acta Biomembr. 2022, 1864, 183889. [Google Scholar] [CrossRef]
  38. Galiano, V.; Villalain, J. Aggregation of 25-hydroxycholesterol in a complex biomembrane. Differences with cholesterol. Biochim. Biophys. Acta Biomembr. 2020, 1862, 183413. [Google Scholar] [CrossRef]
  39. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.C.; Davila-Contreras, E.M.; Qi, Y.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J. Comput. Chem. 2014, 35, 1997–2004. [Google Scholar] [CrossRef] [Green Version]
  40. Murzyn, K.; Rog, T.; Jezierski, G.; Takaoka, Y.; Pasenkiewicz-Gierula, M. Effects of phospholipid unsaturation on the membrane/water interface: A molecular simulation study. Biophys. J. 2001, 81, 170–183. [Google Scholar] [CrossRef] [Green Version]
  41. Kosinova, P.; Berka, K.; Wykes, M.; Otyepka, M.; Trouillas, P. Positioning of antioxidant quercetin and its metabolites in lipid bilayer membranes: Implication for their lipid-peroxidation inhibition. J. Phys. Chem. B 2012, 116, 1309–1318. [Google Scholar] [CrossRef] [PubMed]
  42. Galiano, V.; Villalain, J. The Location of the Protonated and Unprotonated Forms of Arbidol in the Membrane: A Molecular Dynamics Study. J. Membr. Biol. 2016, 249, 381–391. [Google Scholar] [CrossRef] [PubMed]
  43. Galiano, V.; Villalain, J. Oleuropein aglycone in lipid bilayer membranes. A molecular dynamics study. Biochim. Biophys. Acta 2015, 1848, 2849–2858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Esposito, D.L.; Nguyen, J.B.; DeWitt, D.C.; Rhoades, E.; Modis, Y. Physico-chemical requirements and kinetics of membrane fusion of flavivirus-like particles. J. Gen. Virol. 2015, 96, 1702–1711. [Google Scholar] [CrossRef] [Green Version]
  45. Kobayashi, T.; Beuchat, M.H.; Chevallier, J.; Makino, A.; Mayran, N.; Escola, J.M.; Lebrand, C.; Cosson, P.; Kobayashi, T.; Gruenberg, J. Separation and characterization of late endosomal membrane domains. J. Biol. Chem. 2002, 277, 32157–32164. [Google Scholar] [CrossRef] [Green Version]
  46. van Meer, G.; Voelker, D.R.; Feigenson, G.W. Membrane lipids: Where they are and how they behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112–124. [Google Scholar] [CrossRef]
  47. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 27–38. [Google Scholar] [CrossRef]
  48. Guixa-Gonzalez, R.; Rodriguez-Espigares, I.; Ramirez-Anguita, J.M.; Carrio-Gaspar, P.; Martinez-Seara, H.; Giorgino, T.; Selent, J. MEMBPLUGIN: Studying membrane complexity in VMD. Bioinformatics 2014, 30, 1478–1480. [Google Scholar] [CrossRef] [Green Version]
  49. Villalain, J. Epigallocatechin-3-gallate location and interaction with late endosomal and plasma membrane model membranes by molecular dynamics. J. Biomol. Struct. Dyn. 2019, 37, 3122–3134. [Google Scholar] [CrossRef]
  50. Giorgino, T. Computing 1-D atomic densities in macromolecular simulations: The Density Profile Tool for VMD. Comput. Phys. Commun. 2014, 185, 317–322. [Google Scholar] [CrossRef] [Green Version]
  51. Baylon, J.L.; Tajkhorshid, E. Capturing Spontaneous Membrane Insertion of the Influenza Virus Hemagglutinin Fusion Peptide. J. Phys. Chem. B 2015, 119, 7882–7893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Kandt, C.; Ash, W.L.; Tieleman, D.P. Setting up and running molecular dynamics simulations of membrane proteins. Methods 2007, 41, 475–488. [Google Scholar] [CrossRef] [PubMed]
  53. Anézo, C.; de Vries, A.H.; Höltje, H.-D.; Tieleman, D.P.; Marrink, S.-J. Methodological issues in lipid bilayers simulations. J. Phys. Chem. B 2003, 107, 9424–9433. [Google Scholar] [CrossRef] [Green Version]
  54. Bera, I.; Klauda, J.B. Molecular Simulations of Mixed Lipid Bilayers with Sphingomyelin, Glycerophospholipids, and Cholesterol. J. Phys. Chem. B 2017, 121, 5197–5208. [Google Scholar] [CrossRef]
  55. Mukhopadhyay, P.; Monticelli, L.; Tieleman, D.P. Molecular dynamics simulation of a palmitoyl-oleoyl phosphatidylserine bilayer with Na+ counterions and NaCl. Biophys. J. 2004, 86, 1601–1609. [Google Scholar] [CrossRef] [Green Version]
  56. Tieleman, D.P.; Marrink, S.J.; Berendsen, H.J. A computer perspective of membranes: Molecular dynamics studies of lipid bilayer systems. Biochim. Biophys. Acta 1997, 1331, 235–270. [Google Scholar] [CrossRef]
  57. Tsai, H.H.; Lee, J.B.; Li, H.S.; Hou, T.Y.; Chu, W.Y.; Shen, P.C.; Chen, Y.Y.; Tan, C.J.; Hu, J.C.; Chiu, C.C. Geometrical effects of phospholipid olefinic bonds on the structure and dynamics of membranes: A molecular dynamics study. Biochim. Biophys. Acta 2015, 1848, 1234–1247. [Google Scholar] [CrossRef] [Green Version]
  58. Bockmann, R.A.; Hac, A.; Heimburg, T.; Grubmuller, H. Effect of sodium chloride on a lipid bilayer. Biophys. J. 2003, 85, 1647–1655. [Google Scholar] [CrossRef] [Green Version]
  59. Tsuchiya, H. Membrane Interactions of Phytochemicals as Their Molecular Mechanism Applicable to the Discovery of Drug Leads from Plants. Molecules 2015, 20, 18923–18966. [Google Scholar] [CrossRef]
  60. Margina, D.; Ilie, M.; Manda, G.; Neagoe, I.; Mocanu, M.; Ionescu, D.; Gradinaru, D.; Ganea, C. Quercetin and epigallocatechin gallate effects on the cell membranes biophysical properties correlate with their antioxidant potential. Gen. Physiol. Biophys. 2012, 31, 47–55. [Google Scholar] [CrossRef] [Green Version]
Membranes 12 00692 g001 550
Figure 1. (A) Chemical and molecular structures of procyanidin C1, PC1 [(-)-epicatechin –(4β→8)-(-)-epicatechin-(4β→8)-(-)-epicatechin]. The initial, t = 0 ns, and final, t = 675 ns, systems for (B) system 1, (C) system 2, (D) system 3, (E) system 4, and (F) system 5. The disposition of the PC1 molecules (numbered) in each one of the systems are also displayed. The PC1 molecules are depicted in licorice representation and the phosphate atoms of the phospholipids, defining the upper and lower boundaries of the membrane, are depicted in VDW representation and an orange colour. The lipid and water molecules and the chloride and sodium ions have been removed for clarity.
Figure 1. (A) Chemical and molecular structures of procyanidin C1, PC1 [(-)-epicatechin –(4β→8)-(-)-epicatechin-(4β→8)-(-)-epicatechin]. The initial, t = 0 ns, and final, t = 675 ns, systems for (B) system 1, (C) system 2, (D) system 3, (E) system 4, and (F) system 5. The disposition of the PC1 molecules (numbered) in each one of the systems are also displayed. The PC1 molecules are depicted in licorice representation and the phosphate atoms of the phospholipids, defining the upper and lower boundaries of the membrane, are depicted in VDW representation and an orange colour. The lipid and water molecules and the chloride and sodium ions have been removed for clarity.
Membranes 12 00692 g001
Membranes 12 00692 g002 550
Figure 2. Time variation of the z-axis distance (middle of the membrane as a reference) for the different PC1 molecules in (A) system 1 (16 molecules), (B) system 2 (8 molecules), (C) system 3 (4 molecules), (D) system 4 (8 molecules), and (E) system 5 (16 molecules). PC1 molecules in the monomer state are depicted in blue, whereas the PC1 molecules forming dimers or higher-state oligomerization are depicted in black. The z-axis distance of the phosphate atoms of the phospholipids, defining the membrane surface, is depicted in red (upper, center, and lower boundaries). See text for details.
Figure 2. Time variation of the z-axis distance (middle of the membrane as a reference) for the different PC1 molecules in (A) system 1 (16 molecules), (B) system 2 (8 molecules), (C) system 3 (4 molecules), (D) system 4 (8 molecules), and (E) system 5 (16 molecules). PC1 molecules in the monomer state are depicted in blue, whereas the PC1 molecules forming dimers or higher-state oligomerization are depicted in black. The z-axis distance of the phosphate atoms of the phospholipids, defining the membrane surface, is depicted in red (upper, center, and lower boundaries). See text for details.
Membranes 12 00692 g002
Membranes 12 00692 g003 550
Figure 3. Average z-axis distance (middle of the membrane as a reference) for the different PC1 molecules and for the last 30 ns of the MD simulation in (A) system 1, (B) system 2, (C) system 3, (D) system 4, and (E) system 5. PC1 molecules in monomer ( Membranes 12 00692 i002) and oligomerized (■) states are shown. The mean plus standard deviation for all the PC1 molecules in the monomer state ( Membranes 12 00692 i003) is shown in (F). The dotted lines define the membrane surfaces and the center of the bilayer (phosphate atoms of the phospholipids).
Figure 3. Average z-axis distance (middle of the membrane as a reference) for the different PC1 molecules and for the last 30 ns of the MD simulation in (A) system 1, (B) system 2, (C) system 3, (D) system 4, and (E) system 5. PC1 molecules in monomer ( Membranes 12 00692 i002) and oligomerized (■) states are shown. The mean plus standard deviation for all the PC1 molecules in the monomer state ( Membranes 12 00692 i003) is shown in (F). The dotted lines define the membrane surfaces and the center of the bilayer (phosphate atoms of the phospholipids).
Membranes 12 00692 g003
Membranes 12 00692 g004 550
Figure 4. Relative percentage of the lipid molecules surrounding the PC1 molecules for the PM-containing systems (■) and the MIT-containing systems ( Membranes 12 00692 i001) for the last 30 ns of the MD simulation. Lines correspond to the global percentage of each one of the lipids in the system (— for PM- and—for MIT-containing systems).
Figure 4. Relative percentage of the lipid molecules surrounding the PC1 molecules for the PM-containing systems (■) and the MIT-containing systems ( Membranes 12 00692 i001) for the last 30 ns of the MD simulation. Lines correspond to the global percentage of each one of the lipids in the system (— for PM- and—for MIT-containing systems).
Membranes 12 00692 g004
Membranes 12 00692 g005 550
Figure 5. A view of all PC1 molecules in the monomer state together with the surrounding lipids (less than 7Å from PC1) for the five different systems (centered on PC1, 20 different superposed representations). (A) PC1 molecules in VDW representation, lipids in both licorice and white transparent surface representation and the phosphate atoms of the phospholipids in blue transparent surface representation. (B) PC1 molecules in VDW representation and the phosphate atoms of the phospholipids in blue opaque surface representation.
Figure 5. A view of all PC1 molecules in the monomer state together with the surrounding lipids (less than 7Å from PC1) for the five different systems (centered on PC1, 20 different superposed representations). (A) PC1 molecules in VDW representation, lipids in both licorice and white transparent surface representation and the phosphate atoms of the phospholipids in blue transparent surface representation. (B) PC1 molecules in VDW representation and the phosphate atoms of the phospholipids in blue opaque surface representation.
Membranes 12 00692 g005
Table
Table 1. Systems and number of components used in this study. The NaCl concentration was 0.15 M. The production trajectories for each one of the five different systems were calculated for 675 ns.
Table 1. Systems and number of components used in this study. The NaCl concentration was 0.15 M. The production trajectories for each one of the five different systems were calculated for 675 ns.
Plasma
Membrane Model		Mitochondrial Membrane Model
SYSTEM 1SYSTEM 2SYSTEM 3SYSTEM 4SYSTEM 5
POPC8080809090
POPE4848486868
PI-3P1616161010
POPS18181822
POPA---22
PSM34343422
CL---2828
CHOL8484842020
Total280280280222222
PC11684816
H2O27,91428,42028,92233,59833,136
Na+145146148176172
Cl−7980829694
Dimensions
x/y/z (Å)		104/104/120112/112/124
Table
Table 2. Number of oligomer types and number of molecules forming them (in parenthesis) for the five systems studied in this work *.
Table 2. Number of oligomer types and number of molecules forming them (in parenthesis) for the five systems studied in this work *.
Number of
MonomerDimerTrimerTetramerOligomerPC1 (Total)
SYSTEM 1
PM		2 (2)2 (4)--1 (10)16
SYSTEM 2
PM		5 (5)-1 (3)---8
SYSTEM 3
PM		4 (4)----4
SYSTEM 4
MIT		6 (6)1 (2)---8
SYSTEM 5
MIT		4 (4)--1 (4)1 (8)16
216341852
* PM: plasma membrane, MIT: mitochondrial membrane.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Villalaín, J. Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes. Membranes 2022, 12, 692. https://doi.org/10.3390/membranes12070692

AMA Style

Villalaín J. Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes. Membranes. 2022; 12(7):692. https://doi.org/10.3390/membranes12070692

Chicago/Turabian Style

Villalaín, José. 2022. "Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes" Membranes 12, no. 7: 692. https://doi.org/10.3390/membranes12070692

APA Style

Villalaín, J. (2022). Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes. Membranes, 12(7), 692. https://doi.org/10.3390/membranes12070692

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