Prerpints.org logo
Preprint
Review

This version is not peer-reviewed.

Dynamic and Energetic Aspects of Carotenoids In-and-Around Model Lipid Membranes Revealed in Molecular Modelling

A peer-reviewed article of this preprint also exists.

Submitted:

25 June 2024

Posted:

26 June 2024

You are already at the latest version

Abstract
In the introductory part of this review, the relevant results of mostly the latest experimental studies on carotenoid absorption and availability in the human body are briefly presented. The subsequent parts summarise such topics as carotenoid carriers, body transport and tissue delivery, before finally reporting on molecular-level studies of carotenoid binding by membrane receptors. The main message of the review is contained in the section describing computational investigations of carotenoid intercalation and dynamic behaviour in lipid bilayers. The relevance of these computational studies lies in showing the direct link between the microscopic behaviour of molecules and the characteristics of their macroscopic ensembles. Furthermore, studying the interactions between carotenoids and lipid bilayers, and certainly proteins, on the molecular- and atomic-level using computational methods facilitates interpretation and explanation of the macroscopic properties of these systems and, hopefully, helps to better understand the biological functions of carotenoids.
Keywords: 
spontaneous membrane intercalation
;  
preferred membrane orientation
;  
free energy gain
;  
timescales
;  
lipid-soluble antioxidant
;  
trans-cis photoisomerisation
Subject: 
Biology and Life Sciences  -   Biophysics

Carotenoids: Basic Information

Carotenoids are pigments synthesised predominantly by photosynthetic organisms, i.e., plants, algae and cyanobacteria, but some non-photosynthetic bacteria and fungi can also do so, e.g., [1,2]. In contrast, animals and humans are unable to synthesise carotenoids de novo even though they are critically important to their nutrition and health. In plants, carotenoids play fundamental roles in photosynthesis and photoprotection, e.g., [3,4,5]. In humans, carotenoids are vital in maintaining good health, e.g., as antioxidants they reduce the risk of various chronic diseases such as cancer and cardiovascular diseases; β- and α-carotenes are precursors of vitamin A, whose deficiency can cause blindness; xanthophylls lutein and zeaxanthin as macular pigments decrease the onset of age-related macular degeneration and cataract, e.g., [2,6,7,8,9]. Lutein, but also zeaxanthin, accumulated in the human brain improve its cognitive functions [10,11,12].
At present, 1204 natural carotenoids have been identified from 722 source organisms. Their chemical structures and other information are collected in the Carotenoids Database (http://carotenoiddb.jp) established by Yabuzaki [13]. Among the known carotenoids, lycopene, β-carotene (β,β-carotene), α-carotene ((6’R)-β,ε-carotene), lutein and zeaxanthin (Figure 1) are indicated as the most beneficial to humans, e.g., [6,14,15]. The key structural element of all of them is a nearly linear [16] nonpolar polyene chain consisting of eight isoprene units (40-carbon carotenoids) with conjugated double bonds in the trans conformation. With the exception of lycopene, these carotenoids have two terminal ionone rings in C6 and C6′ positions of the polyene chain (Figure 1). In the β- and α-carotene the rings are unsubstituted, in lutein and zeaxanthin they are monohydroxylated at positions C3 and C3’. Thus, lycopene, β- and α-carotene are nonpolar and belong to the group of carotenes, whereas lutein and zeaxanthin are oxygenated derivatives of α- and β-carotene, respectively, and belong to the group of xanthophylls. It is worth noting that both halves (C6–C15 and C15′– C6′, Figure 1) of the polyene chain of the carotenoids are related by 180° rotation about a vertical axis placed at the chain centre (C2h symmetry). Details of carotenoid biosynthesis are given in Ref. [17]. Apart from the polyene chain symmetry, lycopene, β-carotene and zeaxanthin are symmetric molecules (C2h symmetry), as both ionone rings of β-carotene and zeaxanthin are the same (β) and lycopene has no terminal rings. In contrast, both α-carotene and lutein have two different rings (β and ε), thus, the molecules are asymmetric. Each of the small differences in the chemical structures of the molecules has its biological relevance and determines carotenoid-protein and carotenoid-lipid bilayer interactions.

The Journey of Carotenoids through the Human Body

As was mentioned above, humans are unable to synthesise carotenoids de novo, so they have to obtain them from their diet. The main source of dietary carotenoids are coloured fruits and green leafy vegetables, e.g., [6]. To become available to the human organism, carotenoids must first be released from the cellular matrix of fruit- and vegetable-based food products. This takes place in the gastrointestinal tract during digestion. As carotenoids are highly nonpolar, to facilitate their absorption in the intestine they are dispersed into the lipid droplets in the stomach. Then, they are incorporated into bile salt and other lipid mixed micelles in the small intestine [18,19,20,21]. Next, free carotenoids are absorbed by intestinal absorptive cells via scavenger receptors or passive diffusion, incorporated into ultra-low-density lipoproteins (chylomicrons) and secreted into the lymph, e.g., [18,19,20,22,23,24,25]. In general, chylomicrons are taken up by the liver where carotenoids are integrated into either very low-density lipoproteins (VLDL), low-density lipoproteins (LDL) or high-density lipoproteins (HDL), e.g., [26,27,28,29] or stored. Roughly, xanthophylls are associated mainly with HDL and carotenes mainly with LDL and VLDL [24,26,27,29,30], although the exact distribution of carotenoids among lipoproteins is not known for certain [28]. Lipoproteins transport carotenoids in the blood and deliver them to specific tissues for more or less selective uptake.
During their digestion, secretion and delivery, carotenoids exchange, in the free form, between tissues and supramolecular ensembles which they are temporarily associated with, several times before they reach the sites of their main biological activity. There is a large number of such sites in the human organism thus, carotenoids are widely distributed among various organs and tissues. However, different tissues accumulate different amounts of carotenoids and for different purposes. For example, the liver collects virtually all carotenoids mainly to transfer them to different lipoproteins and release them into the bloodstream. Adipose tissue accumulates carotenoids rather indiscriminately mainly to store them. In contrast, lutein and zeaxanthin are the only carotenoids that incorporate into the macula lutea, a specialised tissue at the centre of the retina of the human eyes. Additionally, lutein and zeaxanthin are predominant carotenoids in human brain tissue. Carotenoids also accumulate in other tissues, but those mentioned above are the most representative. There are details concerning the tissue distribution of carotenoids and their functions, e.g., in Refs [10,11,24,31,32,33], but as pointed out by Landrum [34] carotenoids are rather non-specifically accumulated in tissues, with important exceptions for lycopene in the prostate and lutein and zeaxanthin in the retina.
From a practical perspective, information on how carotenoids are transferred from the plasma lipoproteins to specific tissues is, undoubtedly, the most useful. It would be natural to expect that there are several transfer mechanisms, specific to the tissue and the type of carrier. Extensive molecular biology research has revealed that the selective transfer of carotenoid molecules from their carriers to tissues occurs via cell membrane receptors, e.g., [35]. In the case of macular pigments, intake of zeaxanthin into the retinal pigment epithelia (RPE) cells proceeds via scavenger receptor class B type 1 (abbreviated as SR-B1, SRB1 or SCARB1), which is a multifunctional receptor that binds a broad range of lipoproteins, with high affinity for HDL [36,37,38]. The structure and the detailed structure-function analysis of the full-length homology model of human SR-B1 predicted by the molecular modelling software package trRosetta [39] are presented in Ref. [40]. However, intake of lutein into the RPE cells proceeds via the LDL receptor (LDLR) [36,37]. The receptor for LDL-bound β-carotene, let alone that for LDL-bound α-carotene, has not yet been identified conclusively, although it is likely that it is LDLR [33,36,41] but SR-B1 is also considered [33,42].
β- and α-carotene are provitamin A carotenoids. In the human intestine about half of them are converted to vitamin A and the other half is absorbed intact. In β-carotene conversion, the molecule is cleaved by β-carotene 15,15′-oxygenase 1 at the C15,C15′ double bond into two retinal molecules that are subsequently converted into two retinol (vitamin A) molecules [43]. Cleavage of α-carotene gives eventually one retinol and one α-retinol molecule. In the liver, retinol associates with soluble retinol binding protein 4 (RBP4). The RBP4-retinol complex is released into the bloodstream and delivered to tissues via stimulated by retinoic acid 6 (STRA6) receptor which transports retinol across the cell membrane [44,45,46]. In peripheral cells, all-trans-retinol (vitamin A) is converted into 11-cis-retinal, a light-sensitive molecule, and into all-trans retinoic acid, a signalling molecule [45]. In the eye, 11-cis-retinal is the chromophore of rod and cone opsins. After binding retinal, apo-protein opsin becomes rod (rhodopsin) or cone (cone opsin) photoreceptor protein. As a signalling molecule, all-trans retinoic acid binds to the ligand-binding domains of two retinoic acid nuclear receptors, RAR and RXR, while RAR/RXR heterodimers bind specifically to DNA [45,47].
RBP4 that binds and delivers retinol to specific tissue does not bind α-retinol. Nevertheless, animals fed with α-retinol maintained about 50% of the growth rate of those fed with retinol [48]. This indicated that even though α-retinol has no vitamin A activity, it probably shares the ability of retinoic acid to act as a transcription factor [48]. This also indicates that α-retinol is transported as α-retinyl esters in lipoproteins, as does retinyl esters, which is a storage form of retinol.

Carotenoid-Receptor Binding: Molecular-Level Studies

The transfer of carotenoids from plasma lipoproteins to their receptors on the cells of the targeted tissues is a molecular-level event, thus, to obtain better knowledge of this process, it has to be studied using methods applicable to this level. The methods employed in such studies are both experimental and computational. Naturally, molecular-level experiments can provide information at different levels of resolution. Thus, experiments carried out at the macroscopic molecular level have provided physicochemical and biochemical data on carotenoid-protein interactions; the latest publications on the subject include, e.g., Refs [28,29,46,49,50,51,52,53,54]. Experiments carried out at the microscopic molecular level, in combination with a computational approach, have provided high-resolution data of the atomic structure of some carotenoid-binding proteins and the location and structure of the carotenoid molecules bound by them, e.g., [35,51,55,56,57,58,59,60,61]. The results of these and earlier studies have greatly enhanced knowledge and deepened understanding of carotenoid absorption, transport and tissue delivery. However, there are still many topics concerning carotenoid binding, transfer and insertion to and between proteins and lipoproteins that have yet to be revealed. It is thus important to study the processes and elucidate their molecular and atomic-level details, as a deeper understanding of them may help in improving the bioavailability of carotenoids in the body [24,62].

Carotenoids: Transfer from Receptor to Membrane

In the body, carotenoids are transported in lipoproteins and from them are transferred to the cell through more or less specific membrane-bound receptors, but the ultimate location of biologically active carotenoids in the cell is predominantly the cell membrane. Unfortunately, there is not much well-confirmed information regarding how they get there. This is mainly because determining the three-dimensional structures of large multi-domain membrane proteins is a complex and difficult process and not always successful; consequently the structure of none of the human carotenoid receptors has been solved so far. Knowledge of the structure will certainly make predicting how a carotenoid molecule inserts into the membrane feasible. Fortunately, important progress has been made recently as the full-length human SR-B1 receptor has been expressed and purified [54] and, moreover, its 3D structure predicted by homology modelling has been published [40]. Actually, both in the experiment and the computation SR-B1 has been identified as a receptor mediating HDL-cholesterol transport, but it is justified to anticipate that in the case of HDL-carotenoid transport the basic structure of SR-B1 is the same [28,51,63]. The homology model structure of SR-B1 has also been predicted with SWISS-MODEL [64] in Ref. [35]} and with MODELLER [65] in Ref. [61]. All three predicted structures show that the extracellular domain of SR-B1 has a large, predominantly hydrophobic [61] cavity that serves as a tunnel through which a lipid (carotenoid) passes from the lipoprotein to the plasma membrane. The putative sites where the lipid (carotenoid) enters and exits the tunnel [35] are marked in Figure 2.
The situation concerning the membrane receptor, STRA6, for retinol binding protein RBP is more fortunate as its dimeric structure has been solved by cryo-electron microscopy [46]. Knowledge of the STRA6 structure has made it possible to indicate the probable routes of retinol transfer from RBP through STRA6 to the membrane [46]. The routes are schematically shown in Figure 3B.
Experimental studies with atomic resolution provide detailed structures of biomolecules but are not able to record their motion simultaneously. This is because they do not possess sufficiently high concurrent spatial and temporal resolution. At present, only computational molecular modelling methods possess such resolution. Unfortunately, computer simulations of the lipid translocation from the donor lipoprotein to the lipid bilayer through the receptor tunnel have not yet been published, so this issue still remains elusive.

Carotenoids in the Membrane: Atomic-Level Motional Studies

Classical molecular dynamics (MD) simulation is one of the computational molecular modelling methods. This widely used method has atomic spatial resolution and sub-picosecond (ps) temporal resolution but has its limitations too. Its main limitations concern the size of the simulated system and the time of its “observation”. Nevertheless, at present, the all-atom MD simulation allows systems to be studied that contain over a million atoms and processes on the µs timescale sampled with 2-4 fs time steps [66]. However, it should be remembered that each model predicted in MD simulation has to be validated through experiment.
The first MD simulation study on the behaviour of carotenoids in lipid bilayers was carried out on a simple model where four β-carotene molecules were inserted into a hydrated palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayer, parallel to the bilayer normal (vertically). During 4-ns MD simulation the molecules tilted by ~20° and located their rings near the positions of the phospholipid carbonyl groups, although the simulation was too short to allow for larger configurational changes in the system [67]. The next MD simulation study focused on spontaneous intercalation into and subsequent orientation of lutein molecules in the POPC bilayer [68]. In each of two bilayer systems, six lutein molecules were placed horizontally on the bilayer surfaces on the side of the hydrating water. In each of them, one molecule intercalated into the bilayer within the first two ns of MD simulation from the β-ring-side and five did not. These ten molecules (five in each system) aggregated in the water and during 200-ns simulation, remained there as aggregates. The remaining two lutein molecules intercalated into the upper leaflet of the bilayer, one vertically and the other horizontally. At ~70-80 ns of MD simulations the molecules changed their positions to horizontal and vertical (transmembrane), respectively, and remained in these positions until the end of the 200-ns simulations. These orientations were mainly stabilised via hydrogen bonds (H-bond) between the lutein OH groups, and, in the case of the vertical orientation, phosphate and carbonyl groups of POPC and water, and in the case of the horizontal orientation, carbonyl groups of POPC and water [68].
The study of lutein intercalation into a POPC bilayer was pursued further in Ref. [69]. In this study, four bilayer systems were built, each containing six lutein molecules. In the initial structures, the luteins were situated almost vertically to the bilayer plane with the “entering ring” of each molecule placed within the regions of the PC phosphate groups and the rest of the molecule immersed in the hydrating water. In two systems, the “entering ring” was the lutein β ring, and in another two, it was the ε ring. The systems were MD simulated for 20–100 ns depending on the time it took for lutein to intercalate fully into the bilayer. Of the 24 lutein molecules present in the systems, twelve either fully or partially spontaneously intercalated in the bilayer, ten from the β ring side (~80%) and two from the ε ring side (~17%). The remaining twelve molecules aggregated in the water phase and stayed there as aggregates during the whole simulation time. Five of the intercalated molecules were analysed in detail. Four of them intercalated vertically, one horizontally. The time needed by a vertically translocating molecule to reach the interface of the other bilayer side, ranged between ~7.5–15 and ~95 ns, both for the β and ε “entering ring”. It is interesting to note that both the β and the ε ring of the lutein that intercalated into the bilayer horizontally were H-bonded with 1-3 water molecules during the whole simulation time. Irrespective of the final orientation of the lutein in the bilayer, both its rings were H-bonded with water and PC polar groups; these interactions stabilised the orientations [69].
To understand the reason for the discrepancy between the number of lutein molecules intercalated with the “entering β ring” and with the “entering ε ring”, the barriers to the β and ε ring lutein intercalation into the bilayer were evaluated from the free energy profiles of transfer of a lutein molecule from the water phase into the bilayer [69]. The profiles were calculated using umbrella sampling simulations [70]. The positions of the barriers to lutein intercalation from the β and the ε ring end coincided with the region of the rigid POPC glycerol groups. The barriers were low, although that encountered by the ε ring was ~2 kcal/mol higher than that encountered by the β ring. The difference in the barrier heights stemmed from different orientations of the β and ε rings relative to the flat polyene chain, which resulted from sp2 hybridisation of the C6 carbon atom and sp3 hybridisation of the C6′ carbon atom. When lutein intercalated into the bilayer from the β ring end, the system gained 19.5 ± 1.0 kcal/mol and from the ε ring end, the gain was 17.9 ± 1.1 kcal/mol [69].
The MD simulation study described in Ref. [71] was first, where reorientation of carotenoid molecules in the lipid bilayer was analysed in detail. In this study, either a single β-carotene or zeaxanthin molecule was inserted into a dimyristoyl-PC (DMPC) bilayer either vertically or horizontally; the latter in the central bilayer region. During 100-ns and 200-ns simulations, the orientation of the β-carotenes varied more than that of the zeaxanthins but the molecules eventually assumed relatively stable orientations and the β-carotenes were more tilted than the zeaxanthins. These results were confirmed in umbrella sampling simulations [70], which provided free energy profiles for the rotation of β-carotene and zeaxanthin from the horizontal to the vertical position in the bilayer. The relatively flat profile for β-carotene between 90 and 35° implied that in this range of angles, β-carotene had no preferential orientation, whereas the preferential orientation of zeaxanthin in the bilayer was ~35°. Moreover, the system gained over 3 kcal/mol when zeaxanthin rotated from 90 to 35°. For both carotenoids, orientations between 35 and 0° were energetically unfavoured. Moreover, on the basis of the model that utilised a detailed conformational analysis of the single-bond torsion angles along the conjugated polyene chain of β-carotene and zeaxanthin, the authors estimated the effective length of the chain of the molecules in the bilayer and in water. The transmembrane position of zeaxanthin made its chain more extended, which might result in its enhanced antioxidant capability. However, this was not the case with β-carotene [71].
Orientations of xanthophylls in the phospholipid bilayer were also analysed in Ref. [72]. In this study, which was carried out both experimentally and computationally, free energy profiles for rotating zeaxanthin and lutein from the horizontal (90°) to the vertical (10°) position in the DMPC bilayer were calculated using umbrella sampling simulations [70]. The profile for the zeaxanthin was similar to that in Ref. [71]. The minor differences concerned the free energy minimum, which was at ~25° in Ref. [72] and the energy gain due to rotation from 90 to 25°, which was ~4 kcal/mol in Ref. [72]. The free energy profile for rotating lutein in the bilayer showed that the energy gain due to rotation was smaller than that for zeaxanthin, of 3 kcal/mol, but the free energy minimum was also at ~25°. As their experiment has no means of determining the molecule’s orientation at the free energy minimum, the authors compared the average values of the xanthophyll tilt obtained from unbiased MD simulations and from experimental measurements; the values turned out to be close to each other, which mutually validated the model and the measurements [72]. The values of the average tilt of lutein derived experimentally and computationally were virtually the same as those, respectively, of zeaxanthin. This result was surprising because, as could be anticipated from the difference in the free energy profiles, the horizontal orientation of lutein in the bilayer was ~10 times more probable than that of zeaxanthin [72].
The study carried out in Ref. [73] was a cutting-edge piece of research performed both experimentally and computationally. Its aim was to identify a regulatory mechanism of controlling the intensity of light reaching photoreceptors in the human eye. This molecular-level mechanism, which took place in the retinal cell membranes, was based on the rotation of xanthophyll molecules induced by trans-cis photoisomerisation. Mutually perpendicular orientations of trans and cis xanthophyll isomers acted as modulators of light absorption. In the computational part of the research, unbiased MD and umbrella sampling [70] simulations were employed to predict the orientation of zeaxanthin and its two isomers, 9-cis and 13-cis, in two bilayers. One of them was composed of dipalmitoyl-PC (DPPC), the other of lipids typical for the retinal cell membrane, distearoyl-PC (DSPC), stearoyl-docosahexaenoyl-phosphatidylethanolamine (SDPE) and SD-phosphatidylserine (SDPS). The simulations confirmed the previous results [72] that trans zeaxanthin in the bilayer was oriented vertically, whereas 9-cis and 13-cis zeaxanthin were in 2 and 45%, respectively, of cases oriented horizontally (Figure 4). Both the simulations and the experimental measurements identified the same molecular mechanisms of controlling the intensity of light reaching photoreceptors in the retina [73].
Orientations of echinenone (β,β-caroten-4-one, 4-keto-β-carotene) and β-carotene in the lipid bilayer were investigated in Ref. [74]. In this study a mixed-lipid bilayer composed of phospholipids present in natural egg yolk, i.e., PC and PE with 18 and 16 carbon atom acyl chains of varying unsaturation, was used. The bilayer containing one molecule of echinenone or β-carotene was MD simulated for 1 µs [74]. The results showed that although β-carotene could reorient in the bilayer rather freely, it resided most often (78% of the simulation time) in the middle of one of the bilayer leaflets oriented almost horizontally, at 80-85° with the bilayer normal. Less often (22% of the simulation time) β-carotene was oriented 40-45° relative to the bilayer normal. Echinenone in the bilayer was in a transmembrane position with a preferential tilt at ~30°. The results for β-carotene were in general agreement with those in Ref. [71] and the orientation and position of echinenone in the bilayer were similar to those of zeaxanthin in Ref. [72], although detailed comparison is not fully justified, as in both studies related, but not the same quantities were analysed and the bilayers used had different lipid compositions.
In Ref. [75] orientations of β-carotene and zeaxanthin in six mixed-lipid bilayers in different phases were compared. The bilayers were composed of POPC, DSPC and cholesterol in different proportions such that some bilayers were in the liquid disordered and some in the liquid ordered phase. In each phase, the distribution of the β-carotene orientations was broader than that of zeaxanthin. Moreover, in each bilayer the carotenoids caused bilayer thinning. This result was at variance with the result in Ref. [76] and with the literature data that both β-carotene and zeaxanthin increased the POPC bilayer thickness. The contradicting results might stem from different lipid compositions of bilayers used in both studies, and possibly from too short an equilibration time of the bilayers in Ref. [75].
The MD simulations described in Ref. [77] were carried out to investigate the reorientational dynamics and preferred orientation of two xanthophylls in four POPC bilayers on the µs timescale. In the initial structures, six molecules of lutein or zeaxanthin were placed in the bilayer either vertically or horizontally. The molecules oriented horizontally were located in the region between POPC glycerol and phosphate groups, thus, no molecules aggregated in the water phase. Each system was MD simulated for 1.1 µs. The molecules could reorient freely, although none of the vertically placed xanthophyll molecules reoriented to the horizontal position. In contrast, all horizontally placed zeaxanthin molecules reoriented to the vertical position during MD simulation times ranging from 10 to 500 ns. Whereas five out of six horizontally placed lutein molecules changed their orientation to vertical; the times of full reorientation were in the same range as those for zeaxanthin. One lutein molecule remained in the horizontal position for the whole MD simulation time [77] (Figure 5).
To answer the question as to why none of the zeaxanthin molecules remained in the horizontal orientation, whereas one out of six lutein molecules did remain so, detailed conformational analyses of the C5-C6-C7-C8 (β ring) and C5′-C6′-C7′-C8′ (ε ring) torsion angles (called β-ring torsion and ε-ring torsion, respectively) were performed. The β-ring torsion of lutein and zeaxanthin in vacuum, water and the bilayer had two low energy conformations that were close to each other, of 30° and −30°, and separated by a relatively low energy barrier; therefore, transitions between them were very frequent. The ε-ring torsion of lutein in vacuum had two low energy conformations, of 130° and −50°. Even though the angular distance between them was quite large and they were separated by a relatively high barrier, transitions between them were moderately frequent. However, for lutein in water and in the bilayer, the torsion occupied only the lower energy conformation of 130°. H-bonds of lutein OH groups with polar groups of water and PC together with a relatively high energy barrier, hindered transitions between the conformational states. Moreover, analyses of the interactions between water molecules and the methyl groups of the polyene chain as well as the OH groups of the ionone rings, of horizontally oriented lutein and zeaxanthin revealed which interactions with water and which orientations of the rings played a role in stabilising the horizontal orientation of the molecules in the bilayer. The result that the ε-ring torsion of lutein in the bilayer occupied only one conformation was contrary to a commonly made claim that “free” rotation of the ε ring about the single C6′–C7′ bond was the possible cause of the horizontal location of lutein in the bilayer. Conversely, the fixed conformation of the ε-ring torsion, the mutually perpendicular orientations of the lutein ε and β rings (the β ring is coplanar with the polyene chain plane, the ε ring is nearly perpendicular to it), the orientation of the β ring OH group of the horizontal lutein, were the key factors that made a horizontal lutein less likely to rotate to the vertical position than a horizontal zeaxanthin [77]. Of note was the observation that the plane of the polyene chain of any molecule that entered the bilayer from the horizontal orientation was perpendicular to the bilayer surface.
The phospholipid bilayers used in the MD simulation studies described above were, broadly speaking, simple models of an unspecified animal cell membrane. In Ref. [78] the behaviour of two carotenoids, lycopene and zeaxanthin, in a computer model of the human stratum corneum (SC) cell membrane was investigated. In human skin, as in other organs, carotenoids play a protective role against photodamage [79]. The main lipid species of the lipid matrix of the SC cell membrane is ceramides with acyl chains of different lengths. In this study, the lipid bilayer was built of ceramide NS24 (sphingolipid with a 24-carbon-atom fatty acid linked via an amide). The bilayer with one inserted molecule of lycopene or zeaxanthin was MD simulated for 100 ns. The preferred orientation of zeaxanthin in the ceramide bilayer was similar to that in the PC bilayer [71,72,73,77], i.e., vertical [78]. In contrast, the preferred orientation of lycopene was horizontal. The results of unbiased MD simulations of bilayers with a single carotenoid molecule were confirmed by umbrella sampling simulations [70]. The free energy profile for rotating zeaxanthin from 90 to 0° indicated that its preferred orientation in the bilayer was ~25° and that the system gained over 3 kcal/mol due to rotation from 90 to 25°; these results conformed well to those in Refs [71,72,73] even though the bilayers used those studies had different lipid compositions. Rotating lycopene from 90 to 0° resulted in an almost monotonic increase in the system free energy, which indicated that the preferred orientation of lycopene in the bilayer was ~90°, i.e., horizontal [78].
MD simulation of β-carotene in a bilayer made of lipids typical for the brain, POPC and DMPS, was carried out in Ref. [80]. Details concerning the simulation system were nevertheless scarce; the study revealed that even though (an unspecified number of) β-carotene molecules located initially in the water hydrating the bilayer aggregated, they, over time, spontaneously intercalated into the bilayer as an aggregate. As can be seen in the video attached to Ref. [80], a β-carotene molecule that intercalated into the bilayer at the beginning of the simulation as a single species, during (probably) 200-ns MD simulation, had considerable rotational freedom there but was in the vertical position quite often. Once the aggregate intercalated into the bilayer, the single molecule joined it.
A combined experimental and theoretical study of the aggregation of carotenoids in the bilayer was carried out in Ref. [81]. Experiments that employed several methods were performed on chiral carotenoids that differed in polarity, i.e., α-carotene (nonpolar), zeaxanthin and fucoxanthin (polar), in the DPPC bilayer. MD simulations and quantum-chemical calculations were made only on an α-carotene monomer and dimer inserted horizontally into the middle of the DPPC bilayer. Both in a monomeric and a dimeric state, the molecules moved about and rotated in the bilayer core rather freely, thus the “dimer” they formed was not firm, and their tilts relative to the bilayer normal were broadly spread, between 30 and 140°, with an average of 100°. Even though the molecules in the “dimer” were on average ~21 Å apart from each other, they somehow affected each other, e.g., they became more bent than as monomers. The difference in the bending between monomeric and dimeric α-carotene expressed in Å was very small, although it was quite apparent in QM calculations. The stabilising effect of bending on the “dimer” structure shown in simulation and confirmed experimentally indicated the propensity of α-carotene to aggregate in the bilayer in contrast to more polar and vertically oriented zeaxanthin and fucoxanthin [81].

Carotenoids in the Membrane: A Summary of Computational Studies and Their Significance

What have we learnt from computer modelling studies about the behaviour of 40-carbon carotenoids near and in a lipid bilayer? Starting from the beginning: when carotenoid molecules are put into water hydrating the bilayer, some of them aggregate there but some spontaneously intercalate into the bilayer. Both outcomes are due to their long nonpolar polyene chains. Unfortunately, simulations do not conclusively resolve whether the molecules intercalate the bilayer as monomers or as larger aggregates. While in some simulations intercalation is monomeric, some others indicate that multimers can also intercalate. In an MD simulation study [82] lutein intercalated into the POPC bilayer both as a monomer and a dimer, although in the subsequent MD simulations the latter never happened [68,69,77]. A large β-carotene aggregate intercalating into the POPC-DMPS bilayer was shown in Ref. [80]. An experiment, in contrast to a simulation, has no means to tell what the multimeric form of carotenoids intercalating the bilayer is, but is able to discriminate between carotenoid monomers and multimers in the bilayer and also the type of aggregate formed, e.g., [81,83,84]. α-Carotene self-assembles readily in the gel-phase DPPC bilayer and forms “head-to-tail” J-aggregates of uncertain sizes that remain stable when temperature is increased [81]. Carotenoids with polar groups, namely, fucoxanthin, zeaxanthin, lutein, form “side-by-side” H-aggregates of rather small sizes that are unstable at elevated temperatures and in liquid-phase bilayers they are in the monomeric form [81,83,84]. This experimental result conforms to the results of 1.1-μs MD simulations of zeaxanthin and lutein in the POPC bilayer at 310 K where the molecules did not show any sign of aggregation [77].
Spontaneity of xanthophyll intercalation into the lipid bilayer directly stems from the negative value of the change in free energy for the process. The calculated free energy gain in lutein (zeaxanthin) intercalating the bilayer is ~20 kcal/mol [69]. MD simulations and free energy calculations also demonstrate that the intercalation of lutein from the β ring end is more probable than from the ε ring end; the barrier for the ε ring to pass through the PC glycerol groups region of the bilayer is ~2 kcal/mol higher than that for the β ring [69].
MD simulations confirm experimental results that the vertical position of xanthophyll molecules in the bilayer is far more probable than the horizontal one, whereas umbrella sampling simulations give numerical values of the free energy gain due to the rotation of lutein and zeaxanthin from the horizontal to the vertical positions; the values are in the range <3, 5> kcal/mol [71,72,73,78]. The gain, which is altogether not very large, is smaller for lutein than zeaxanthin, which makes the horizontal orientation of lutein in the bilayer ~10 times more probable than that of zeaxanthin [72]. These energetical evaluations are in harmony with the results of MD simulations in Ref. [77], where out of twelve xanthophyll molecules placed horizontally in the bilayer, eleven spontaneously reoriented to vertical position and the molecule that did not reorient was lutein. In contrast to polar, nonpolar carotenoids, β-carotene, α-carotene and lycopene are mainly in the horizontal position in the lipid bilayer; this is evident from experimental measurements, MD simulations [71,74,78,81] and free energy calculations [71,78]. Echinenone with only one polar group is in the PC-PE bilayer [74] orientated similarly to that of zeaxanthin in the PC bilayer [71].
The answer to the question as to why the vertical or horizontal orientation of carotenoids in the membrane is important should be provided by experimental or quantum mechanical studies. Here one can only speculate, bearing in mind that the main function of carotenoids is to protect the cell from the damaging effects of light and oxygen. Let us take macular pigments as an example. The maximum absorption of light by macular pigments is at 460 nm, the wavelength corresponding to blue light, which is damaging to retinal cells. The cornea and lens of the eye absorb 99.34% of incident UV light so in practice it does not reach the retina [85]. Macular pigments absorb ~60% of the incident blue light at 460 nm and 40-46% of the entire damaging wavelength range from 400 to 500 nm [86,87]. The highest concentration of macular pigments is in the photoreceptor axon layer, which precedes the photoreceptors; macular pigments thus act as an effective pre-receptoral light filter [87,88]. The transmembrane (vertical) position of lutein and zeaxanthin makes it possible to fulfil their task. However, the task would presumably be accomplished more effectively if the xanthophylls were in the horizontal position, which would allow them to absorb blue light from all directions [89,90]. MD simulations indicate that this might be the case only for lutein and probably on a very short timescale [69,77]. However, recent results of the Gruszecki group [73] demonstrate that light-induced cis-isomerisation of xanthophylls in human retina samples and in model membranes trigger their reorientation from a vertical to a horizontal position. The consequences of the reorientation are briefly discussed above and in detail in Refs [73,91].
Blue light that reaches the retina may initiate photosensitised production of singlet oxygen. Thus xanthophylls as pre-receptoral light filters act as passive antioxidants [34,88]. They are also active antioxidants as they are able to quench singlet oxygen and other reactive oxygen species [34,88,90,91,92]. However, the ability of macular pigments to protect the retina stems mainly from their selective location in the domains of retinal cell membranes that are enriched in phospholipids with long polyunsaturated acyl chains, which are especially vulnerable to oxygen damage. In these domains, they act as a lipid-soluble antioxidant and their vertical as well as horizontal position there plays a crucial role in phospholipid protection [91,93,94,95,96].
Computer simulations supplement experimental data not only by providing numerical values for the free energy gain on xanthophyll intercalation or rotation in the membrane, but also by predicting the timescales of the processes. Full intercalation of a lutein molecule into the POPC bilayer is very fast and ranges between ~7.5–15 and ~95 ns for either β or ε ring intercalation [69]. The same time range applies to zeaxanthin. The predicted time of lutein and zeaxanthin reorientation from the horizontal to the vertical position is in the range from 10 to 500 ns [77]. One lutein molecule did not reorient during the simulation time of 1.1 µs in Ref. [77]. In classical MD simulation, such an observation time is considered rather long but it is very short compared to that of experimental measurements on macroscopic samples. As lutein in a horizontal position has not been observed experimentally [73], it is quite probable that in a longer MD simulation the horizontal lutein reorients to the vertical position. However, one should then recall that photochemical reactions are very fast. Photoexcitation of lutein populates its second singlet excited state, S2, that decays in a few hundred femtoseconds to the lowest singlet excited state, S1, which then relaxes back to the ground state, S0, on a picosecond time scale [97,98]. Thus, during the lifetime of its horizontal orientation (>1 µs) lutein may participate in many light generated processes. Besides, even though the horizontal position is much less probable than the vertical one, both orientations must be in some statistical equilibrium in the bilayer.
This review clearly shows that a combined experimental and computational approach in studying molecular-level biological processes provides us with the opportunity to know and understand them better. It also shows that our knowledge concerning how carotenoids are transported in the body and get into the cell membrane is rather fragmentary. However, hopefully, this review might help the reader in recognising topics that merit more in-depth study.

Author Contributions

M.P-G wrote the manuscript. J.H. wrote a draft of the paragraph on carotenoid absorption and availability in the human body; M.M. collected, summarised and critically evaluated the literature on the simulation studies on carotenoid-membrane interactions and prepared figures. All authors contributed to the final version of the manuscript.

Acknowledgements

We gratefully acknowledge the Polish high-performance computing infrastructure PLGrid (ACK Cyfronet AGH) for providing computer facilities and support within computational grant no. PLG/2022/015999. MP-G thanks W. K. Subczynski for constructive discussions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fuentes, P.; Pizarro, L.; Moreno, J.C.; Handford, M.; Rodriguez-Concepcion, M.; Stange, C. Light-dependent changes in plastid differentiation influence carotenoid gene expression and accumulation in carrot roots. Plant Mol Biol 2012, 79, 47-59. [CrossRef]
  2. Sun, T.H.; Rao, S.M.; Zhou, X.S.; Li, L. Plant carotenoids: recent advances and future perspectives. Mol Hortic 2022, 2, ARTN 3. [CrossRef]
  3. Ballottari, M.; Mozzo, M.; Girardon, J.; Hienerwadel, R.; Bassi, R. Chlorophyll Triplet Quenching and Photoprotection in the Higher Plant Monomeric Antenna Protein Lhcb5. J. Phys. Chem. B 2013, 117, 11337-11348. [CrossRef]
  4. Dall’Osto, L.; Lico, C.; Alric, J.; Giuliano, G.; Havaux, M.; Bassi, R. Lutein is needed for efficient chlorophyll triplet quenching in the major LHCII antenna complex of higher plants and effective photoprotection in vivo under strong light. BMC Plant Biol. 2006, 6, 32.
  5. Xiao, F.G.; Shen, L.; Ji, H.F. On photoprotective mechanisms of carotenoids in light harvesting complex. Biochem Bioph Res Co 2011, 414, 1-4. [CrossRef]
  6. Eggersdorfer, M.; Wyss, A. Carotenoids in human nutrition and health. Arch. Biochem. Biophys. 2018, 652, 18-26. [CrossRef]
  7. Milani, A.; Basirnejad, M.; Shahbazi, S.; Bolhassani, A. Carotenoids: biochemistry, pharmacology and treatment. Brit J Pharmacol 2017, 174, 1290-1324. [CrossRef]
  8. Fiedor, J.; Burda, K. Potential Role of Carotenoids as Antioxidants in Human Health and Disease. Nutrients 2014, 6, 466-488. [CrossRef]
  9. Zhao, Z.L.; Chen, J.; Ci, F.F.; Pang, H.; Cheng, N.; Xing, A.J. α-Carotene: a valuable carotenoid in biological and medical research. J Sci Food Agr 2022, 102, 5606-5617. [CrossRef]
  10. Stringham, J.M.; Johnson, E.J.; Hammond, B.R. Lutein across the Lifespan: From Childhood Cognitive Performance to the Aging Eye and Brain. Curr Dev Nutr 2019, 3, ARTN nzz066. [CrossRef]
  11. Lopresti, A.L.; Smith, S.J.; Drummond, P.D. The Effects of Lutein and Zeaxanthin Supplementation on Cognitive Function in Adults With Self-Reported Mild Cognitive Complaints: A Randomized, Double-Blind, Placebo-Controlled Study. Front Nutr 2022, 9, ARTN 843512. [CrossRef]
  12. Parekh, R.; Hammond, B.R.; Chandradhara, D. Lutein and Zeaxanthin Supplementation Improves Dynamic Visual and Cognitive Performance in Children: A Randomized, Double-Blind, Parallel, Placebo-Controlled Study. Adv Ther 2024. [CrossRef]
  13. Yabuzaki, J. Carotenoids Database: structures, chemical fingerprints and distribution among organisms. Database-Oxford 2017, ARTN bax004. [CrossRef]
  14. Amengual, J. Bioactive Properties of Carotenoids in Human Health. Nutrients 2019, 11, ARTN 2388. [CrossRef]
  15. Pitts, S.B.J.; Johnson, N.S.; Wu, Q.; Firnhaber, G.C.; Kaur, A.P.; Obasohan, J. A meta-analysis of studies examining associations between resonance Raman spectroscopy-assessed skin carotenoids and plasma carotenoids among adults and children. Nutr Rev 2022, 80, 230-241. [CrossRef]
  16. Fiedor, L.; Dudkowiak, A.; Pilch, M. The origin of the dark S state in carotenoids: a comprehensive model. J R Soc Interface 2019, 16, ARTN 20190191. [CrossRef]
  17. Moise, A.R.; Al-Babili, S.; Wurtzel, E.T. Mechanistic Aspects of Carotenoid Biosynthesis. Chem Rev 2014, 114, 164-193. [CrossRef]
  18. Harrison, E.H. Mechanisms involved in the intestinal absorption of dietary vitamin A and provitamin A carotenoids. Bba-Mol Cell Biol L 2012, 1821, 70-77. [CrossRef]
  19. Lemmens, L.; Colle, I.; Van Buggenhout, S.; Palmero, P.; Van Loey, A.; Hendrickx, M. Carotenoid bioaccessibility in fruit- and vegetable-based food products as affected by product (micro) structural characteristics and the presence of lipids: A review. Trends Food Sci Tech 2014, 38, 125-135. [CrossRef]
  20. Xavier, A.A.O.; Mercadante, A.Z. The bioaccessibility of carotenoids impacts the design of functional foods. Curr Opin Food Sci 2019, 26, 1-8. [CrossRef]
  21. Parker, R.S. Carotenoids .4. Absorption, metabolism, and transport of carotenoids. FASEB J. 1996, 10, 542-551. [CrossRef]
  22. Chacón-Ordóñez, T.; Carle, R.; Schweiggert, R. Bioaccessibility of carotenoids from plant and animal foods. J Sci Food Agr 2019, 99, 3220-3239. [CrossRef]
  23. Reboul, E. Proteins involved in fat-soluble vitamin and carotenoid transport across the intestinal cells: New insights from the past decade. Prog. Lipid Res. 2023, 89, ARTN 101208. [CrossRef]
  24. Furr, H.C.; Clark, R.M. Transport, uptake, and target. Tissue storage of carotenoids. In Carotenoids in Health and Disease, Krinsky, N.I., Mayne, S.T., Sies, H., Eds.; Marcel Dekker: New York, 2004; pp. 229-278.
  25. Gugliucci, A. The chylomicron saga: time to focus on postprandial metabolism. Front Endocrinol 2024, 14, ARTN 1322869. [CrossRef]
  26. Shmarakov, I.O.; Yuen, J.J.; Blaner, W.S. Carotenoid Metabolism and Enzymology. In Carotenoids and Human Health, Tanumihardjo, S.A., Ed.; Humana Press (Springer): New York, 2013; pp. 29-56.
  27. Canene-Adams, K.; Erdman Jr., J.W. Absorption, transport, distribution in tissues and bioavailability. In Carotenoids, Nutrition and Health, Britton, G., Liaaen-Jensen, S., Pfander, H., Eds.; Birkhäuser Verlag: Basel, 2009; Volume 5, pp. 115-144.
  28. Li, B.X.; Vachali, P.; Chang, F.Y.; Gorusupudi, A.; Arunkumar, R.; Shi, L.J.; Rognon, G.T.; Frederick, J.M.; Bernstein, P.S. HDL is the primary transporter for carotenoids from liver to retinal pigment epithelium in transgenic ApoA-I/Bco2 mice. Arch. Biochem. Biophys. 2022, 716, ARTN 109111. [CrossRef]
  29. Arunkumar, R.; Gorusupudi, A.; Bernstein, P.S. The macular carotenoids: A biochemical overview. Bba-Mol Cell Biol L 2020, 1865, ARTN 158617. [CrossRef]
  30. Krinsky, N.I.; Cornwell, D.G.; Oncley, J.L. The Transport of Vitamin-a and Carotenoids in Human Plasma. Arch. Biochem. Biophys. 1958, 73, 233-246. [CrossRef]
  31. Böhm, V.; Lietz, G.; Olmedilla-Alonso, B.; Phelan, D.; Reboul, E.; Bánati, D.; Borel, P.; Corte-Real, J.; de Lera, A.R.; Desmarchelier, C.; et al. From carotenoid intake to carotenoid blood and tissue concentrations - implications for dietary intake recommendations. Nutr Rev 2021, 79, 544-573. [CrossRef]
  32. Elvira-Torales, L.I.; García-Alonso, J.; Periago-Castón, M.J. Nutritional Importance of Carotenoids and Their Effect on Liver Health: A Review. Antioxidants-Basel 2019, 8, ARTN 229. [CrossRef]
  33. Bohn, T.; Desmarchelier, C.; El, S.N.; Keijer, J.; van Schothorst, E.; Rühl, R.; Borel, P. β-Carotene in the human body: metabolic bioactivation pathways - from digestion to tissue distribution and excretion. P Nutr Soc 2019, 78, 68-87. [CrossRef]
  34. Landrum, J.T. Reactive Oxygen and Nitrogen Species in Biological Systems: Reactions and Regulation by Carotenoids. In Carotenoids and Human Health, Tanumihardjo, S.A., Ed.; Humana Press (Springer) New York, 2013; Volume 1, pp. 57-101.
  35. Bandara, S.; von Lintig, J. Aster la vista: Unraveling the biochemical basis of carotenoid homeostasis in the human retina. Bioessays 2022, 44, ARTN e2200133. [CrossRef]
  36. Thomas, S.E.; Harrison, E.H. Mechanisms of selective delivery of xanthophylls to retinal pigment epithelial cells by human lipoproteins. J Lipid Res 2016, 57, 1865-1878. [CrossRef]
  37. Harrison, E.H. Mechanisms of Transport and Delivery of Vitamin A and Carotenoids to the Retinal Pigment Epithelium. Mol Nutr Food Res 2019, 63, ARTN 1801046. [CrossRef]
  38. Shen, W.J.; Asthana, S.; Kraemer, F.B.; Azhar, S. Scavenger receptor B type 1: expression, molecular regulation, and cholesterol transport function. J Lipid Res 2018, 59, 1114-1131. [CrossRef]
  39. Yang, J.Y.; Anishchenko, I.; Park, H.; Peng, Z.L.; Ovchinnikov, S.; Baker, D. Improved protein structure prediction using predicted interresidue orientations. P Natl Acad Sci USA 2020, 117, 1496-1503. [CrossRef]
  40. Powers, H.R.; Sahoo, D. SR-B1’s Next Top Model: Structural Perspectives on the Functions of the HDL Receptor. Curr Atheroscler Rep 2022, 24, 277-288. [CrossRef]
  41. Quadro, L.; Giordano, E.; Costabile, B.K.; Nargis, T.; Iqbal, J.; Kim, Y.; Wassef, L.; Hussain, M.M. Interplay between β-carotene and lipoprotein metabolism at the maternal-fetal barrier. Bba-Mol Cell Biol L 2020, 1865, ARTN 158591. [CrossRef]
  42. Borel, P.; Moussa, M.; Reboul, E.; Lyan, B.; Defoort, C.; Vincent-Baudry, S.; Maillot, M.; Gastaldi, M.; Darmon, M.; Portugal, H.; et al. Human plasma levels of vitamin E and Carotenoids are associated with genetic polymorphisms in genes involved in lipid metabolism. J Nutr 2007, 137, 2653-2659. [CrossRef]
  43. Reboul, E. Mechanisms of Carotenoid Intestinal Absorption: Where Do We Stand? Nutrients 2019, 11, ARTN 838. [CrossRef]
  44. Steinhoff, J.S.; Lass, A.; Schupp, M. Biological Functions of RBP4 and Its Relevance for Human Diseases. Front Physiol 2021, 12, ARTN 659977. [CrossRef]
  45. Ask, N.M.; Leung, M.; Radhakrishnan, R.; Lobo, G.P. Vitamin A Transporters in Visual Function: A Mini Review on Membrane Receptors for Dietary Vitamin A Uptake, Storage, and Transport to the Eye. Nutrients 2021, 13, ARTN 3987. [CrossRef]
  46. Chen, Y.; Clarke, O.B.; Kim, J.; Stowe, S.; Kim, Y.K.; Assur, Z.; Cavalier, M.; Godoy-Ruiz, R.; von Alpen, D.C.; Manzini, C.; et al. Structure of the STRA6 receptor for retinol uptake. Science 2016, 353, ARTN aad8266. [CrossRef]
  47. Petkovich, M.; Chambon, P. Retinoic acid receptors at 35 years. J Mol Endocrinol 2022, 69, T13-T24. [CrossRef]
  48. Riabroy, N.; Dever, J.T.; Tanumihardjo, S.A. α-Retinol and 3,4-didehydroretinol support growth in rats when fed at equimolar amounts and α-retinol is not toxic after repeated administration of large doses. Br. J. Nutr. 2014, 111, 1373-1381. [CrossRef]
  49. Likkei, K.; Moldenhauer, M.; Tavraz, N.N.; Maksimov, E.G.; Sluchanko, N.N.; Friedrich, T. Lipid composition and properties affect protein-mediated carotenoid uptake efficiency from membranes. Biochim. Biophys. Acta-Biomem. 2024, 1866, ARTN 184241. [CrossRef]
  50. Shyam, R.; Vachali, P.; Gorusupudi, A.; Nelson, K.; Bernstein, P.S. All three human scavenger receptor class B proteins can bind and transport all three macular xanthophyll carotenoids. Arch. Biochem. Biophys. 2017, 634, 21-28. [CrossRef]
  51. Bandara, S.; Ramkumar, S.; Imanishi, S.; Thomas, L.D.; Sawant, O.B.; Imanishi, Y.; von Lintig, J. Aster proteins mediate carotenoid transport in mammalian cells. P Natl Acad Sci USA 2022, 119, ARTN e2200068119. [CrossRef]
  52. Slonimskiy, Y.B.; Egorkin, N.A.; Friedrich, T.; Maksimov, E.G.; Sluchanko, N.N. Microalgal protein AstaP is a potent carotenoid solubilizer and delivery module with a broad carotenoid binding repertoire. Febs J 2022, 289, 999-1022. [CrossRef]
  53. Bandara, S.; Moon, J.; Ramkumar, S.; von Lintig, J. ASTER-B regulates mitochondrial carotenoid transport and homeostasis. J Lipid Res 2023, 64, ARTN 100369. [CrossRef]
  54. Powers, H.R.; Jenjak, S.E.; Volkman, B.F.; Sahoo, D. Development and validation of a purification system for functional full-length human SR-B1 and CD36. Journal of Biological Chemistry 2023, 299, ARTN 105187. [CrossRef]
  55. Sluchanko, N.N.; Slonimskiy, Y.B.; Egorkin, N.A.; Varfolomeeva, L.A.; Kleymenov, S.Y.; Minyaev, M.E.; Faletrov, Y.V.; Moysenovich, A.M.; Parshina, E.Y.; Friedrich, T.; et al. Structural basis for the carotenoid binding and transport function of a START domain. Structure 2022, 30, 1647-+. [CrossRef]
  56. Kornilov, F.D.; Slonimskiy, Y.B.; Lunegova, D.A.; Egorkin, N.A.; Savitskaya, A.G.; Kleymenov, S.Y.; Maksimov, E.G.; Goncharuk, S.A.; Mineev, K.S.; Sluchanko, N.N. Structural basis for the ligand promiscuity of the neofunctionalized, carotenoid-binding fasciclin domain protein AstaP. Commun Biol 2023, 6, ARTN 471. [CrossRef]
  57. Horvath, M.P.; George, E.W.; Tran, Q.T.; Baumgardner, K.; Zharov, G.; Lee, S.; Sharifzadeh, H.; Shihab, S.; Mattinson, T.; Li, B.X.; et al. Structure of the lutein-binding domain of human StARD3 at 1.74 angstrom resolution and model of a complex with lutein. Acta Crystallogr F 2016, 72, 609-618. [CrossRef]
  58. Shyam, R.; Gorusupudi, A.; Nelson, K.; Horvath, M.P.; Bernstein, P.S. RPE65 has an additional function as the lutein to meso-zeaxanthin isomerase in the vertebrate eye. P Natl Acad Sci USA 2017, 114, 10882-10887. [CrossRef]
  59. Varfolomeeva, L.A.; Slonimskiy, Y.B.; Egorkin, N.A.; Minyaev, M.E.; Faletrov, Y.V.; Boyko, K.M.; Sluchanko, N.N. Preparation and Structural Studies of the Silkworm Carotenoid-Binding Protein Complexed with a New Pigment. Crystallogr Rep+ 2022, 67, 909-917. [CrossRef]
  60. Hazai, E.; Bikádi, Z.; Zsila, F.; Lockwood, S.F. Molecular modeling of the non-covalent binding of the dietary tomato carotenoids lycopene and lycophyll, and selected oxidative metabolites with 5-lipoxygenase. Bioorgan Med Chem 2006, 14, 6859-6867. [CrossRef]
  61. Neculai, D.; Schwake, M.; Ravichandran, M.; Zunke, F.; Collins, R.F.; Peters, J.; Neculai, M.; Plumb, J.; Loppnau, P.; Pizarro, J.C.; et al. Structure of LIMP-2 provides functional insights with implications for SR-BI and CD36. Nature 2013, 504, 172-176. [CrossRef]
  62. Molteni, C.; La Motta, C.; Valoppi, F. Improving the Bioaccessibility and Bioavailability of Carotenoids by Means of Nanostructured Delivery Systems: A Comprehensive Review. Antioxidants-Basel 2022, 11, ARTN 1931. [CrossRef]
  63. Widjaja-Adhi, M.A.K.; Lobo, G.P.; Golczak, M.; Von Lintig, J. A genetic dissection of intestinal fat-soluble vitamin and carotenoid absorption. Hum Mol Genet 2015, 24, 3206-3219. [CrossRef]
  64. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res 2018, 46, W296-W303. [CrossRef]
  65. Fiser, A.; Sali, A. MODELLER: Generation and refinement of homology-based protein structure models. Method Enzymol 2003, 374, 461-491. [CrossRef]
  66. Rizzuti, B. Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena. Bba-Proteins Proteom 2022, 1870, ARTN 140757. [CrossRef]
  67. Jemiola-Rzeminska, M.; Pasenkiewicz-Gierula, M.; Strzalka, K. The behaviour of β-carotene in the phosphatidylcholine bilayer as revealed by a molecular simulation study. Chem. Phys. Lipids 2005, 135, 27-37. [CrossRef]
  68. Pasenkiewicz-Gierula, M.; Baczynski, K.; Murzyn, K.; Markiewicz, M. Orientation of lutein in a lipid bilayer - revisited. Acta Biochim. Pol. 2012, 59, 115-118.
  69. Makuch, K.; Markiewicz, M.; Pasenkiewicz-Gierula, M. Asymmetric Spontaneous Intercalation of Lutein into a Phospholipid Bilayer, a Computational Study. Comput Struct Biotec 2019, 17, 516-526. [CrossRef]
  70. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187-199. [CrossRef]
  71. Cerezo, J.; Zuniga, J.; Bastida, A.; Requena, A.; Ceron-Carrasco, J.P. Conformational changes of beta-carotene and zeaxanthin immersed in a model membrane through atomistic molecular dynamics simulations. Phys. Chem. Chem. Phys. 2013, 15, 6527-6538. [CrossRef]
  72. Grudzinski, W.; Nierzwicki, L.; Welc, R.; Reszczynska, E.; Luchowski, R.; Czub, J.; Gruszecki, W.I. Localization and Orientation of Xanthophylls in a Lipid Bilayer. Sci. Rep. 2017, 7, 9619. [CrossRef]
  73. Luchowski, R.; Grudzinski, W.; Welc, R.; Pinto, M.M.M.; Sek, A.; Ostrowski, J.; Nierzwicki, L.; Chodnicki, P.; Wieczor, M.; Sowinski, K.; et al. Light-Modulated Sunscreen Mechanism in the Retina of the Human Eye. J. Phys. Chem. B 2021, 125, 6090-6102. [CrossRef]
  74. Semenov, A.N.; Gvozdev, D.A.; Zlenko, D.V.; Protasova, E.A.; Khashimova, A.R.; Parshina, E.Y.; Baizhumanov, A.A.; Lotosh, N.Y.; Kim, E.E.; Kononevich, Y.N.; et al. Modulation of Membrane Microviscosity by Protein-Mediated Carotenoid Delivery as Revealed by Time-Resolved Fluorescence Anisotropy. Membranes-Basel 2022, 12, ARTN 905. [CrossRef]
  75. Johnson, Q.R.; Mostofian, B.; Gomez, G.F.; Smith, J.C.; Cheng, X.L. Effects of carotenoids on lipid bilayers. Phys. Chem. Chem. Phys. 2018, 20, 3795-3804. [CrossRef]
  76. Mostofian, B.; Johnson, Q.R.; Smith, J.C.; Cheng, X.L. Carotenoids promote lateral packing and condensation of lipid membranes. Phys. Chem. Chem. Phys. 2020, 22, 12281-12293. [CrossRef]
  77. Makuch, K.; Hryc, J.; Markiewicz, M.; Pasenkiewicz-Gierula, M. Lutein and Zeaxanthin in the Lipid Bilayer-Similarities and Differences Revealed by Computational Studies. Front Mol Biosci 2021, 8, ARTN 768449. [CrossRef]
  78. Ri, J.S.; Choe, C.S.; Choe, S.H.; Jong, K.H.; Hong, S.N.; Schleusener, J.; Lademann, J.; Darvin, M.E. Lycopene, but not zeaxanthin, serves as a skeleton for the formation of an orthorhombic organization of intercellular lipids within the lamellae in the stratum corneum: Molecular dynamics simulations of the hydrated ceramide NS bilayer model. Biochim. Biophys. Acta-Biomem. 2023, 1865. [CrossRef]
  79. Flieger, J.; Raszewska-Famielec, M.; Radzikowska-Büchner, E.; Flieger, W. Skin Protection by Carotenoid Pigments. Int J Mol Sci 2024, 25, ARTN 1431. [CrossRef]
  80. Gastaldo, I.P.; Himbert, S.; Ram, U.; Rheinstädter, M.C. The Effects of Resveratrol, Caffeine, β-Carotene, and Epigallocatechin Gallate (EGCG) on Amyloid-β Aggregation in Synthetic Brain Membranes. Mol Nutr Food Res 2020, 64, ARTN 2000632. [CrossRef]
  81. Hachlica, N.; Stefanska, M.; Mach, M.; Kowalska, M.; Wydro, P.; Domagala, A.; Kessler, J.; Zajac, G.; Kaczor, A. Organization of Carotenoid Aggregates in Membranes Studied Selectively using Resonance Raman Optical Activity. Small 2024. [CrossRef]
  82. Pasenkiewicz-Gierula, M.; Markiewicz, M. Orientation of lutein in a lipid bilayer – revisited. In Proceedings of the 16th International Symposium on Carotenoids., Krakow, Poland,, 2012; pp. 105-105.
  83. Sujak, A.; Okulski, W.; Gruszecki, W.I. Organisation of xanthophyll pigments lutein and zeaxanthin in lipid membranes formed with dipalmitoylphosphatidylcholine. Biochim. Biophys. Acta-Biomem. 2000, 1509, 255-263. [CrossRef]
  84. Sujak, A.; Mazurek, P.; Gruszecki, W.I. Xanthophyll pigments lutein and zeaxanthin in lipid multibilayers formed with dimyristoylphosphatidylcholine. J. Photochem. Photobiol. B 2002, 68, 39-44. [CrossRef]
  85. Fukuoka, H.; Gali, H.E.; Bu, J.J.; Sella, R.; Afshari, N.A. Ultraviolet light exposure and its penetrance through the eye in a porcine model. Int J Ophthalmol-Chi 2023, 16, 172-177. [CrossRef]
  86. Snodderly, D.M.; Auran, J.D.; Delori, F.C. The Macular Pigment .2. Spatial-Distribution in Primate Retinas. Invest Ophth Vis Sci 1984, 25, 674-685.
  87. Schalch, W.; Landrum, J.T.; Bone, R.A. The Eye. In Carotenoids, Nutrition and Health, Britton, G., Liaaen-Jensen, S., Pfande, H., Eds.; Birkhäuser Verlag Basel, 2009; Volume 5, pp. 301-330.
  88. O’Connell, E.; Neelam, K.; Nolan, J.; Eong, K.G.A.; Beatty, S. Macular carotenoids and age-related maculopathy. Ann Acad Med Singap 2006, 35, 821-830.
  89. Ahmed, S.; Lott, M.; Marcus, D. The macular xantophyls. Survey of Ophthalmology 2005, 50, 183-193.
  90. Vishwanathan, R.; Johnson, E.J. Lutein and Zeaxanthin and Eye Disease. In Carotenoids and Human Health, Tanumihardjo, S.A., Ed.; Humana Press (Springer): New York, 2013; Volume 1, pp. 215-235.
  91. Welc-Stanowska, R.; Pietras, R.; Mielecki, B.; Sarewicz, M.; Luchowski, R.; Widomska, J.; Grudzinski, W.; Osyczka, A.; Gruszecki, W.I. How Do Xanthophylls Protect Lipid Membranes from Oxidative Damage? J Phys Chem Lett 2023. [CrossRef]
  92. Cantrell, A.; McGarvey, D.J.; Truscott, T.G.; Rancan, F.; Böhm, F. Singlet oxygen quenching by dietary carotenoids in a model membrane environment. Arch. Biochem. Biophys. 2003, 412, 47-54. [CrossRef]
  93. Subczynski, W.K.; Wisniewska, A.; Widomska, J. Location of macular xanthophylls in the most vulnerable regions of photoreceptor outer-segment membranes. Arch. Biochem. Biophys. 2010, 504, 61-66. [CrossRef]
  94. Widomska, J.; Gruszecki, W.I.; Subczynski, W.K. Factors Differentiating the Antioxidant Activity of Macular Xanthophylls in the Human Eye Retina. Antioxidants-Basel 2021, 10, Artn 601. [CrossRef]
  95. Widomska, J.; Zareba, M.; Subczynski, W.K. Can Xanthophyll-Membrane Interactions Explain Their Selective Presence in the Retina and Brain? Foods 2016, 5, ARTN 7. [CrossRef]
  96. Wisniewska, A.; Subczynski, W.K. Distribution of macular xanthophylls between domains in a model of photoreceptor outer segment membranes. Free Radic. Biol. Med. 2006, 41, 1257-1265.
  97. Polívka, T.; Sundström, V. Ultrafast dynamics of carotenoid excited states -: From solution to natural and artificial systems. Chem Rev 2004, 104, 2021-2071. [CrossRef]
  98. Accomasso, D.; Arslancan, S.; Cupellini, L.; Granucci, G.; Mennucci, B. Ultrafast Excited-State Dynamics of Carotenoids and the Role of the S State. J Phys Chem Lett 2022, 13, 6762-6769. [CrossRef]
Preprints 110361 g001
Figure 1. The chemical structures of lutein, zeaxanthin, β-carotene, α-carotene and lycopene. The ionone ring and polyene chain atoms that are used in the text are indicated in lutein only. The atoms are numbered according to the IUPAC convention. The C6′ and C7′ atoms of lutein are indicated with red circles. The chemical symbols for carbon atoms, C, and the hydrogen atoms have been omitted except for the OH groups of lutein and zeaxanthin.
Figure 1. The chemical structures of lutein, zeaxanthin, β-carotene, α-carotene and lycopene. The ionone ring and polyene chain atoms that are used in the text are indicated in lutein only. The atoms are numbered according to the IUPAC convention. The C6′ and C7′ atoms of lutein are indicated with red circles. The chemical symbols for carbon atoms, C, and the hydrogen atoms have been omitted except for the OH groups of lutein and zeaxanthin.
Preprints 110361 g001
Preprints 110361 g002
Figure 2. Homology model of class B scavenger receptor. Extracellular domain model of human SCARB1 (SR-B1) receptor with proposed lipid entering and exiting tunnel highlighted with a grey colour mesh. Figure. From Bandara, S., & von Lintig, J. (2022). Aster la vista: Unraveling the biochemical basis of carotenoid homeostasis in the human retina. BioEssays, 44, e2200133. https://doi.org/10.1002/bies.202200133.
Figure 2. Homology model of class B scavenger receptor. Extracellular domain model of human SCARB1 (SR-B1) receptor with proposed lipid entering and exiting tunnel highlighted with a grey colour mesh. Figure. From Bandara, S., & von Lintig, J. (2022). Aster la vista: Unraveling the biochemical basis of carotenoid homeostasis in the human retina. BioEssays, 44, e2200133. https://doi.org/10.1002/bies.202200133.
Preprints 110361 g002
Preprints 110361 g003
Figure 3. Possible mechanism for STRA6-mediated retinol uptake. (A) Ribbon representation of STRA6 (red) and RBP (blue). Retinol is shown as green spheres. The lipid bilayer is shaded light yellow. (B) Schematic of STRA6-mediated retinol (green spheres) release from RBP into the outer cleft and translocation to the lipid bilayer (shown as purple spheres and wavy lines) through the lateral window. Question marks indicate putative ligand binding sites; green and white arrows show two potential retinol exit pathways. (C) Previous work mapped onto the STRA6 structure, details in Ref. [46]. Figure. From: Yunting Chen et al., Structure of the STRA6 receptor for retinol uptake. Science 353, aad8266 (2016). doi:10.1126/science.aad8266. Reprinted with permission from AAAS.
Figure 3. Possible mechanism for STRA6-mediated retinol uptake. (A) Ribbon representation of STRA6 (red) and RBP (blue). Retinol is shown as green spheres. The lipid bilayer is shaded light yellow. (B) Schematic of STRA6-mediated retinol (green spheres) release from RBP into the outer cleft and translocation to the lipid bilayer (shown as purple spheres and wavy lines) through the lateral window. Question marks indicate putative ligand binding sites; green and white arrows show two potential retinol exit pathways. (C) Previous work mapped onto the STRA6 structure, details in Ref. [46]. Figure. From: Yunting Chen et al., Structure of the STRA6 receptor for retinol uptake. Science 353, aad8266 (2016). doi:10.1126/science.aad8266. Reprinted with permission from AAAS.
Preprints 110361 g003
Preprints 110361 g004
Figure 4. Representative structures for the perpendicular orientation of all-trans zeaxanthin (left) and the horizontal orientations of its 9-cis and 13-cis isomers (right) in the DPPC membrane. Figure. Form J. Phys. Chem. B 2021, 125, 23, 6090-6102. doi: 10.1021/acs.jpcb.1c01198. CC-BY 4.0.
Figure 4. Representative structures for the perpendicular orientation of all-trans zeaxanthin (left) and the horizontal orientations of its 9-cis and 13-cis isomers (right) in the DPPC membrane. Figure. Form J. Phys. Chem. B 2021, 125, 23, 6090-6102. doi: 10.1021/acs.jpcb.1c01198. CC-BY 4.0.
Preprints 110361 g004
Preprints 110361 g005
Figure 5. Initial (left column) and final (right column; after 1,100 ns of MD simulations) structures of the POPC bilayer containing lutein (upper row) and zeaxanthin (lower row). In the initial structure, six xanthophyll molecules are placed parallel to the bilayer surface (horizontally). The lutein molecule which remained in the horizontal position during the whole simulation time is presented in the ball-and-stick model. In the figures, only the POPC head groups and xanthophyll molecules are shown, to better illustrate the details. The atoms are represented in standard colours, except for the xanthophyll carbon atoms, which are blue for lutein and red for zeaxanthin. Water is shaded light blue and the lipid nonpolar region is shaded light yellow. Figure. Form Makuch K, Hryc J, Markiewicz M and Pasenkiewicz-Gierula M (2021) Lutein and Zeaxanthin in the Lipid Bilayer–Similarities and Differences Revealed by Computational Studies. Front. Mol. Biosci. 8:768449. doi: 10.3389/fmolb.2021.768449. CC-BY.
Figure 5. Initial (left column) and final (right column; after 1,100 ns of MD simulations) structures of the POPC bilayer containing lutein (upper row) and zeaxanthin (lower row). In the initial structure, six xanthophyll molecules are placed parallel to the bilayer surface (horizontally). The lutein molecule which remained in the horizontal position during the whole simulation time is presented in the ball-and-stick model. In the figures, only the POPC head groups and xanthophyll molecules are shown, to better illustrate the details. The atoms are represented in standard colours, except for the xanthophyll carbon atoms, which are blue for lutein and red for zeaxanthin. Water is shaded light blue and the lipid nonpolar region is shaded light yellow. Figure. Form Makuch K, Hryc J, Markiewicz M and Pasenkiewicz-Gierula M (2021) Lutein and Zeaxanthin in the Lipid Bilayer–Similarities and Differences Revealed by Computational Studies. Front. Mol. Biosci. 8:768449. doi: 10.3389/fmolb.2021.768449. CC-BY.
Preprints 110361 g005
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
bsky
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated