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.
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, S
2, that decays in a few hundred femtoseconds to the lowest singlet excited state, S
1, which then relaxes back to the ground state, S
0, 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.