Detail Analysis for Energy Transfer and Pigment Assembling in C − Phycocyanin through Time − Resolved Spectroscopy and AlphaFold2

: The time-resolved ﬂuorescence spectroscopy of the C − phycocyanin (C − PC) complex from the cyanobacterium Spirulina platensis ( S. platensis ) hexamer with the sub-picosecond resolution was detected, and the energy transfer pathways and related transfer rates were identiﬁed through the multiexponential analysis based on Monte-Carlo method. With the purpose to construct the relationship between the functions and the structure in vivo, the three − dimension (3D) protein structure was predicted via AlphaFold2, and the arrangement of chromophores treated as the energy transfer nodes were obtained. The experiment results have been matched well with the structure prediction. This work suggests a new way to investigate structure prediction in vivo and the corresponding functions.


Introduction
As the most important chemical reaction on earth, photosynthesis allows photosynthetic organisms including plants, algae, cyanobacteria and photosynthetic bacteria to capture and transfer the light's energy into chemical energy [1]. Among these photosynthetic organisms, algae represent is regarded as the intermediate stage in the evolution from photosynthetic bacteria in the ocean to green plants on the land. In algae, the function of light harvesting is fulfilled by the pigment−protein complex called phycobiliproteins (PBPs), including phycocyanin (PC), phycoerythrin (PE) and allophycocyanin (APC) [2,3]. Experimental results show that different PBPs may contain different kinds of chromophores and the chromophores are classified by structure as phycoerythrobilin (PEB), phycocyanobilin (PCB), phycoviolobilin (PVB) and phycourobilin (PUB) [4,5].
In the antenna system, it was clear that energy is transferred in and between the phycobiliproteins at an extremely fast speed, which occurs with a remarkably high conversion efficiency of 90% or higher. The efficient energy transfer process is closely related to the unique structure of PBPs. The early results from the absorption and emission spectra of PBPs indicate that the direction of excitation energy transfer is usually from the short−wavelength band to the long-wavelength band in the antenna system [6]. Corresponding to the structure, energy is transferred from the phycoerythrin, to the phycocyanin, the allophycocyanin, and reaches the reaction centre in sequence [7][8][9][10][11]. However, if one ever wants to mimic the natural process of light−harvesting in an artificial system, it is not enough to know the general direction of energy transfer, the key elements of the process should be fully understood in the view of physics, chemistry and biology. It is also necessary to understand the precise structure to identify the specific path of energy transfer. Commonly, the assignment of energy transfer paths relies on hyperfine protein structures obtained by X−ray crystallography or Cryo−EM and time-domain sequential fluorescence events obtained by time-resolved spectroscopy [12,13]. Techniques such as Cryo−EM and crystallization have accumulated enough key structural information for protein structure. In addition, with the development of computation-based protein structure prediction technology, using more rapid and informative means to efficiently predict protein in vivo may provide more detailed and realistic information for the study of the efficient energy transfer mechanism.
Recent advances in deep learning, combined with the availability of genomic data for inferring co−evolutionary patterns enlightening a new way for protein structure determination to complete the longstanding physics−based approaches. The outstanding performance of AlphaFold2 demonstrates the remarkable power of deep learning in protein structure prediction, which can provide a structural basis for understanding the energy transfer dynamics in PBPs [14]. C−PC is a typical phycobiliprotein in cyanobacteria that thrives around hot springs, with three PCB chromophores attached through thioester linkages at the α84, β84 and β155 positions. Previous research has discovered the energy transfer dynamics among the three PCB chromophores in C−phycocyanin, via theory calculation based on the Förster resonance mechanism and experimental measurement using ultra−fast spectroscopy technology [15][16][17][18][19][20][21][22]. To explore the accuracy and applicability of the AlphaFold2 calculation results in the energy transfer process, we selected the C−PC as the experimental sample.
In order to achieve our research, we first selected a C−PC hexamer from S. platensis to investigate the structure-based energy transfer mechanism. Secondly, AlphaFold2 was applied to simulate the 3D protein structure of the C−PC complex for predicting the plausible structure basis energy transfer routes among the chromophores. Thirdly, sub-picosecond time-resolved fluorescence spectra were detected for obtaining the overall direct experimental evidence for energy transfer routes in the C−PC hexamer. Finally, the energy transfer kinetics between different chromophores were analysed by comparing the structural prediction results and deconvolution of time-resolved fluorescence spectra. Therefore, we obtained the meaningful energy transfer kinetics between different chromophores in C−phycocyanin in the view of the microscopic mechanism of the original process of photosynthesis and found that it may be instructive to design and construct the semi-artificial and artificial solar-energy systems modelled on the photosynthetic processes utilized by photosynthetic organisms at nearly 100−percent efficiency of light harvesting.

Growth of Red Algae and Isolation of PBSs
Spirulina platensis was selected as the experiment material and cultured at 30 • C under constant low light illumination (40 W white fluorescent tubes, 16 h light and 8 h dark as the ideal photoperiod) in Zarrouk medium. Contamination of the medium was inhibited by high-pressure steam sterilization and the air supply with 0.22 µm paper filters. Spirulina platensis was harvested by centrifugation after 2 weeks. C-PC hexamer from S. platensis was prepared according to the methods described in our previous studies [23,24]. The harvested algae were added to pH 7.0, 0.002 mol/L phosphate buffer (containing 0.004 mol/L sodium azide) at a ratio of 1:15, stirred well, put into liquid nitrogen to completely frozen, then thawed at 4 • , repeated 3 times, centrifuged at 20,000 r/min, 4 • for 10 min, and the supernatant was taken out. To the upper (NH 4 ) 2 SO 4 , a final concentration of 1.25 mo/L was added to the upper layer of the clear solution, mixed well, left to stand at 4 • for 12 h, centrifuged and the upper layer of the clear solution was taken. The supernatant is the crude extract of PC. The Macro−Prep Methyl hydrophobic column was equilibrated with the above phosphate buffer, and then the crude extract was loaded onto the column and eluted with a gradient of 1.25, 0.9, 0.5, and 0.2 mol of (NH 4 ) 2 SO 4 phosphate buffer, respectively, and the eluate with a purity (A 620 /A 230 ) greater than 4.0 was collected and concentrated by desalting with a 30 kD ultrafiltration tube and then stored at 4 • C.

Determination and Analysis of Time-Resolved Spectra
Time-resolved fluorescence spectra (TRFS) were measured at 77 K by a synchroscan streak camera (Hamamatsu C6860, time−resolution 700 fs) coupled to a polychromator. The silica cuvette containing samples was placed in liquid nitrogen and frozen in the dark. The applied Ti: sapphire laser was selected as the excitation at the wavelength was 570 nm. Spectral data were extracted with a 0.2 nm interval, and TRFS were reconstructed afterwards. Fluorescence delay-associated spectra (FDAS) were calculated based on a global optimization method [25]. The overall time resolution was less than 2 ps. All the decays were re−analysed using the same lifetime parameters. The amplitudes of these exponential components as a function of emission wavelength provided FDAS [26]. In the process of deconvolving the fluorescence spectrum detected by the streak camera, the detected fluorescence intensity F exp can be described as follows: where f pump represents the pump laser pulse. Here we considered the theoretical fluorescence intensity F theo as a sum of multi-exponential form: We applied a deconvolution procedure based on the global optimization method described as follows to fit the isotropic fluorescence. Letter A denotes the amplitude of the fluorescence isotropic decay constant of ε i .

Protein Structure Prediction of C−PC Complex Based on AlphaFold2
Two different subunits form the monomer of C−PC, which serves as the basic unit for the PBP assembly. The amino acid sequence of the two types of C−PC subunits from S. platensis has been obtained (Protein Data Bank in Europe No. p72509 and p72508, Protein ID: 1gh0). The structure of the C−PC complex was predicted with the latest 2.1 version of AlphaFold2. AlphaFold2 runs on the supercomputing cluster provided by the Beikun cloud supercomputing platform with hardware support of 8 vCPUs, 64G of memory and two V100 graphics cards. In addition, since the model does not contain phycobilin pigments, it is necessary to use the molecular docking software, Autodock, to attach the pigments to the protein model [27,28]. MD analysis and NumPy Python packages were applied for the calculation of the mean pLDDT value of each residue extracted from the B-factor column of AlphaFold2 structure files. Molecular visualization and RMSD calculation were performed using PyMOL(The PyMol Molecular Graphics System, Version 2.4.0, Schrödinger Inc, New York, USA) [29].

Steady State and Time−Resolved Spectroscopy Results
The steady-state spectra of the C−PC hexamer were measured at room temperature and 77 K, respectively. The results are shown in Figure 1a. The black line represents the absorption spectra of C−PC hexamer at room temperature, which exhibited a peak at 619 nm and a shoulder at 570 nm, with a large half-width owing to the absorption of PCB chromophores. At room temperature, the maximum fluorescence peak occurred at 643 nm, as shown in the blue curve in Figure 1a. It is inferred that the C−PC hexamer was a completely independent unit both in structure and function. The steady-state spectra at room temperature are consistent with the result from Zhao et al. [19]. files represent the fluorescence intensity recorded at various delay times after excitation. Such representation helps us to understand the energy flow more conceptually. On excitation at 570 nm that preferentially excited PCB chromophores, the time-resolved fluorescence spectra showed an initial spike of the PC fluorescence band around 650 nm, followed by the red shift of the fluorescence band. At 280 ps after excitation, the time-resolved fluorescence spectra showed the maximum intensity at the wavelength of around 560 nm, indicating the energy has been transferred from the PC fluorescence chromophore of 550 nm to the chromophore of 560 nm.

Protein Structure Prediction Results
In the study of photosynthesis, the resonance energy transfer theory is usually used to describe the energy transfer between phycobilin pigments. Förster theory gives the relationship between energy transfer efficiency and orientation factor, the distance of pigment-pigment and spectral overlap integral. In addition, the reciprocal of the sixth power The fluorescence spectra of the C−PC hexamer measured at the temperature of 77 K upon excitation at 570 nm are also shown as the red line in Figure 1a, which exhibits a narrower spectra band and 15 nm peak red-shift compared with the steady-state fluorescence at room temperature. This is due to the low temperature at 77 K that can lead to the re−conformation of phycobiliprotein, resulting in the change of the electronic state of the chromophore group [30]. The emission at the wavelength region longer than 700 nm might relate to the vibrational bands.
Fluorescence spectra in a steady state are a combination of emissions from fluorescence components. Energy transfer among them shall be described on the basis of individual kinetics. Thus, the time-resolved fluorescence spectra of C−PC were analysed in detail, as shown in Figure 1b. The time-resolved fluorescence spectrum was conducted at 77 K, at which the molecular vibration effect and solvent effect that affect the precise reflection of the energy transfer process could be neglected, allowing large signal-to-noise ratios to guarantee that the time-resolved fluorescence spectrum is close to the necessary assumption of the energy transfer in theoretical models. Figure 1c shows the reconstructed three-dimensional time-resolved fluorescence spectra of C−PC from S. platensis at 77 K, represented in colour gradation. The colour profiles represent the fluorescence intensity recorded at various delay times after excitation. Such representation helps us to understand the energy flow more conceptually. On excitation at 570 nm that preferentially excited PCB chromophores, the time-resolved fluorescence spectra showed an initial spike of the PC fluorescence band around 650 nm, followed by the red shift of the fluorescence band. At 280 ps after excitation, the time-resolved fluorescence spectra showed the maximum intensity at the wavelength of around 560 nm, indicating the energy has been transferred from the PC fluorescence chromophore of 550 nm to the chromophore of 560 nm.

Protein Structure Prediction Results
In the study of photosynthesis, the resonance energy transfer theory is usually used to describe the energy transfer between phycobilin pigments. Förster theory gives the relationship between energy transfer efficiency and orientation factor, the distance of pigment-pigment and spectral overlap integral. In addition, the reciprocal of the sixth power of the distance between the pigments has a linear relationship with the energy transfer efficiency. Distance strongly affects energy transfer efficiency. As a novel neural network architecture, AlphaFold2 has greatly improved in the field of the accuracy of structure prediction. Here, its structure prediction has been treated as the model for the energy transfer function of C−PC from S. platensis. Predicting protein structure using Alphafold2 requires the amino acid sequence of the protein. Therefore, the FASTA file with amino acid information was proved, which contains monomers, trimers and hexamers of C−PC. Alphafold2 will build protein models from those FASTA files and reorders the sequence numbers of these predicted structures with the confidence of the model. Typically, the model with the highest ordinal number corresponds to the highest confidence. As shown in Figure 2, two different subunits (α and β subunits) form the αβ monomer, which serves as the basic unit for the PBP assembly. Three αβ monomers further assemble into a ring−shaped (αβ) 3 trimer, and two trimers then stack to make a rodlike (αβ) 6 hexamer. The α−chain is 162 residues long and the β−chain is 172 residues long, as shown in Figure  S1 which has been added in Supplementary Materials. of the distance between the pigments has a linear relationship with the energy transfer efficiency. Distance strongly affects energy transfer efficiency. As a novel neural network architecture, AlphaFold2 has greatly improved in the field of the accuracy of structure prediction. Here, its structure prediction has been treated as the model for the energy transfer function of C−PC from S. platensis. Predicting protein structure using Alphafold2 requires the amino acid sequence of the protein. Therefore, the FASTA file with amino acid information was proved, which contains monomers, trimers and hexamers of C−PC. Alphafold2 will build protein models from those FASTA files and reorders the sequence numbers of these predicted structures with the confidence of the model. Typically, the model with the highest ordinal number corresponds to the highest confidence. As shown in Figure 2, two different subunits (α and β subunits) form the αβ monomer, which serves as the basic unit for the PBP assembly. Three αβ monomers further assemble into a ring−shaped (αβ)3 trimer, and two trimers then stack to make a rodlike (αβ)6 hexamer. The α−chain is 162 residues long and the β−chain is 172 residues long, as shown in Figure S1 which has been added in Supplement Materials. Although a C−PC hexamer model was obtained by AlphaFold2, the confidence of the model requires further characterisation. With regards to this, the output file of Al-phaFold2 was processed, the b−factor of the model was extracted, and the pLDDT curve of the model was drawn. The average confidence for the monomer was 0.9642, while the averages for the trimer and hexamer were 0.9705 and 0.9401. Subsequently, the Al-phaFold2 model was aligned with the crystallographic model, and the RMSD values of the two models were calculated by the PyMol script (RMSD = 0.963), as shown in Supplementary Materials Figure S2a. The two results indicated that the computational model of AlphaFold2 has a high match with the crystallographic model. However, it can be seen from the value of the b-factor for the single amino acid site that the AlphaFold2 model has a small structure with low confidence, as shown in Supplementary Materials Figure S2b. Any variation in protein structure may lead to significant changes in the distance between pigments, which in turn alters the efficiency of energy transfer. For this reason, the cysteines of the dangling pigment on the crystallographic model were aligned with the cysteines at the same position on AlphaFold2. From the comparison results, the lower confidence amino acid site was not near the cysteine. It was found that its position was farther from the cysteine bound to phycocyanin, as shown in Supplementary Materials Figure  S2c.

Chromophore Arrangement Results
As described above, the possible energy transfer routes could be predicted from the distances among the chromophores and the energy transfer efficiency is inversely related to the distance between the chromophores. However, since phycobilin pigments are not Although a C−PC hexamer model was obtained by AlphaFold2, the confidence of the model requires further characterisation. With regards to this, the output file of AlphaFold2 was processed, the b−factor of the model was extracted, and the pLDDT curve of the model was drawn. The average confidence for the monomer was 0.9642, while the averages for the trimer and hexamer were 0.9705 and 0.9401. Subsequently, the AlphaFold2 model was aligned with the crystallographic model, and the RMSD values of the two models were calculated by the PyMol script (RMSD = 0.963), as shown in Supplementary Materials Figure S2a. The two results indicated that the computational model of AlphaFold2 has a high match with the crystallographic model. However, it can be seen from the value of the b-factor for the single amino acid site that the AlphaFold2 model has a small structure with low confidence, as shown in Supplementary Materials Figure S2b. Any variation in protein structure may lead to significant changes in the distance between pigments, which in turn alters the efficiency of energy transfer. For this reason, the cysteines of the dangling pigment on the crystallographic model were aligned with the cysteines at the same position on AlphaFold2. From the comparison results, the lower confidence amino acid site was not near the cysteine. It was found that its position was farther from the cysteine bound to phycocyanin, as shown in Supplementary Materials Figure S2c

Chromophore Arrangement Results
As described above, the possible energy transfer routes could be predicted from the distances among the chromophores and the energy transfer efficiency is inversely related to the distance between the chromophores. However, since phycobilin pigments are not included in the model, it is necessary to perform molecular docking on the model. Normally, phycobilisomes hang from proteins in the form of thioether bonds, and the monomers that make up the C−PC have three hanging sites. Figure 3 shows the PCB chromophores arrangement in C−PC. Three PCB chromophores α84, β84 and β155 are covalently attached to cysteine residues by thioester bonds in the αβ monomer. α84−PCB is associated with Cys84 in the α-chain, while the β84−PCB and β155−PCB are linked to Cys82 and Cys153 in the β−chain, respectively, as shown in Figure 3a. The PCB chromophores arrangement in the (αβ) 3 trimer was also shown in Figure 3b. Since the C−PC hexamer is assembled from two trimers in a face−to−face manner, there is some overlap in the positions of the pigments, as shown in Figure 3c. included in the model, it is necessary to perform molecular docking on the model. Normally, phycobilisomes hang from proteins in the form of thioether bonds, and the monomers that make up the C−PC have three hanging sites. Figure 3 shows the PCB chromophores arrangement in C−PC. Three PCB chromophores α84, β84 and β155 are covalently attached to cysteine residues by thioester bonds in the αβ monomer. α84−PCB is associated with Cys84 in the α-chain, while the β84−PCB and β155−PCB are linked to Cys82 and Cys153 in the β−chain, respectively, as shown in Figure 3a. The PCB chromophores arrangement in the (αβ)3 trimer was also shown in Figure 3b. Since the C−PC hexamer is assembled from two trimers in a face−to−face manner, there is some overlap in the positions of the pigments, as shown in Figure 3c. According to the AlphaFold2 predicted protein structure and the further covalent bond docking, eighteen phycocyanins were suspended at various active sites of C−PC. The distances between different pigments can be read out by the molecular visualization software, PyMol. Table 1 shows the information on the distance distribution between different pigments read out by PyMol. The results show that the distance distribution between pigments is less than 40 Å, whether in the monomer, trimer or hexamer of C−PC. In addition, the distance between β84 and β155 is about 35 Å in the αβ monomer. However, compared to the monomer, the α84 and β84 in the trimer have a much shorter distance between adjacent monomers and the shortest distance of about 20 Å. The distance between β-84 pairs in the adjacent monomer is about 35 Å, and the α84 and β155 in the adjacent monomer are about 39 Å. Similarly, in the (αβ)6 hexamer, the distance distribution between α84 pairs, β84 pairs and β155 pairs in two face−to−face trimers are 26 Å, 35 Å, and 28 Å, respectively. The distance between α84-β84 pairs and α84-β155 pairs are 35 Å and 37 Å, respectively.  According to the AlphaFold2 predicted protein structure and the further covalent bond docking, eighteen phycocyanins were suspended at various active sites of C−PC. The distances between different pigments can be read out by the molecular visualization software, PyMol. Table 1 shows the information on the distance distribution between different pigments read out by PyMol. The results show that the distance distribution between pigments is less than 40 Å, whether in the monomer, trimer or hexamer of C−PC. In addition, the distance between β84 and β155 is about 35 Å in the αβ monomer. However, compared to the monomer, the α84 and β84 in the trimer have a much shorter distance between adjacent monomers and the shortest distance of about 20 Å. The distance between β-84 pairs in the adjacent monomer is about 35 Å, and the α84 and β155 in the adjacent monomer are about 39 Å. Similarly, in the (αβ) 6 hexamer, the distance distribution between α84 pairs, β84 pairs and β155 pairs in two face−to−face trimers are 26 Å, 35 Å, and 28 Å, respectively. The distance between α84-β84 pairs and α84-β155 pairs are 35 Å and 37 Å, respectively.

Discussion
To determine the constants of fluorescence decay of C−PC, deconvolution was conducted in a global optimization method. The time-resolved fluorescence decay at different detected wavelengths was resolved by multi-exponential deconvolution and the Monte−Carlo method was adopted to process the experimental data, as shown in Figure 4a.  Based on the results of the deconvolution procedure, we can obtain the fluorescence decay-associated spectra (FDAS), which reveal the energy transfer dynamics of the decay constants. Figure 4b shows the FDAS of the C−PC hexamer. The negative and positive bands indicate the rise and decay kinetics in populations in the excited state, respectively. The coupling of negative and positive bands is a clear indication of energy transfer from the chromophore giving the positive band to those giving the negative band.
According to the previous report, the energy transfer process in the monomer and trimer is different due to the change in the C−PC structure. In the αβ monomer, the fluorescence decay constant with the shortest lifetime is 52 ps and was assigned to the energy transfer progress from β155 to β84 [15]. Moreover, from the structural information, the distance between the two chromophores was 35 Å. While in trimers, two new energy transfer paths emerged due to the complexity of the spatial structure: 1 ps decay constant from the α84 to β84 and 40 ps decay constant between identical chromophore pairs. In the adjacent monomer (α1-84 and β2-84), the α84 and β84 have the shortest distance of about 20Å, and the energy transfer time between α1-84 and β2-84 has been identified as less than 1 ps [16]. However, the energy transfer progress which has the fluorescence decay constant of about less than 1 ps cannot be deconvoluted definitely. Moreover, the energy transfer progress between identical chromophore β−84 pairs (β1-84 and β2-84) around the trimer ring has an energy transfer time of 46 ps, with a distance of 35 Å. Therefore, the fluorescence decay constant of 7, 22, 250, and 1560 ps can be identified by combining the data on the distance between the chromophores and the FDAS results. It is reasonable to assign the decay time of 7ps to the energy transfer time between the identical chromophore α-84 pairs (α1-84 to α4-84 and α2-84 to α5-84 and α3-84 to α6-84) of the two trimers in (αβ)6 hexamer, with the chromophore distance of 26 Å. The fluorescence decay constant of 22 ps might be assigned to the energy transfer time between the identical chromophore β-155 pairs (β1-155 to β5-155 and β2-155 to β6-155 and β3-155 to β4-155) of the two trimers in (αβ)6 hexamer, with the chromophore distance of 28 Å. Furthermore, the predominant energy transfer pathways in the C−PC hexamer were likely from β−155(or α−84) in the top trimer to β−155(or α−84) in the adjacent bottom trimer, rather than from β-84 in the top trimer to β-84 in adjacent bottom trimer [31], with the chromophore distance of 35 Å. According to the FDAS results, the fluorescence decay constant of 250 ps might be assigned to the emission of α84 and β84 chromophores, because the two chromophores of α84−PCB and β84−PCB have similar fluorescence spectra, Based on the results of the deconvolution procedure, we can obtain the fluorescence decay-associated spectra (FDAS), which reveal the energy transfer dynamics of the decay constants. Figure 4b shows the FDAS of the C−PC hexamer. The negative and positive bands indicate the rise and decay kinetics in populations in the excited state, respectively. The coupling of negative and positive bands is a clear indication of energy transfer from the chromophore giving the positive band to those giving the negative band.
According to the previous report, the energy transfer process in the monomer and trimer is different due to the change in the C−PC structure. In the αβ monomer, the fluorescence decay constant with the shortest lifetime is 52 ps and was assigned to the energy transfer progress from β155 to β84 [15]. Moreover, from the structural information, the distance between the two chromophores was 35 Å. While in trimers, two new energy transfer paths emerged due to the complexity of the spatial structure: 1 ps decay constant from the α84 to β84 and 40 ps decay constant between identical chromophore pairs. In the adjacent monomer (α1-84 and β2-84), the α84 and β84 have the shortest distance of about 20Å, and the energy transfer time between α1-84 and β2-84 has been identified as less than 1 ps [16]. However, the energy transfer progress which has the fluorescence decay constant of about less than 1 ps cannot be deconvoluted definitely. Moreover, the energy transfer progress between identical chromophore β−84 pairs (β1-84 and β2-84) around the trimer ring has an energy transfer time of 46 ps, with a distance of 35 Å. Therefore, the fluorescence decay constant of 7, 22, 250, and 1560 ps can be identified by combining the data on the distance between the chromophores and the FDAS results. It is reasonable to assign the decay time of 7ps to the energy transfer time between the identical chromophore α-84 pairs (α1-84 to α4-84 and α2-84 to α5-84 and α3-84 to α6-84) of the two trimers in (αβ) 6 hexamer, with the chromophore distance of 26 Å. The fluorescence decay constant of 22 ps might be assigned to the energy transfer time between the identical chromophore β-155 pairs (β1-155 to β5-155 and β2-155 to β6-155 and β3-155 to β4-155) of the two trimers in (αβ) 6 hexamer, with the chromophore distance of 28 Å. Furthermore, the predominant energy transfer pathways in the C−PC hexamer were likely from β−155(or α−84) in the top trimer to β−155(or α−84) in the adjacent bottom trimer, rather than from β-84 in the top trimer to β-84 in adjacent bottom trimer [31], with the chromophore distance of 35 Å. According to the FDAS results, the fluorescence decay constant of 250 ps might be assigned to the emission of α84 and β84 chromophores, because the two chromophores of α84−PCB and β84−PCB have similar fluorescence spectra, and the energy transfer path between α84-PCB and β84−PCB are approximately reversible [32]. The fluorescence decay constant of 1560 ps might reflect the final emission of the anchor polypeptide in the C−PC hexamer.
Accurate identification of the energy transfer process requires precise protein structure and high-time-resolution spectral information. However, the acquisition of both is challenging work. To this end, we used computational methods to assign fluorescence decay constant reasonably between high-resolution protein structures and existing time-resolved fluorescence spectra, and also verified the applicability of AI-based protein structure calculation for the energy transfer process.

Conclusions
The functions of proteins are determined mainly by structures. Therefore, resolving energy transfer processes in light-harvesting complexes relies on accurate structural information. In addition to the traditional methods of experimentally resolving protein structures by nuclear magnetic resonance (NMR), cryo−electron microscopy (cryo−EM) and X−ray crystallography, artificial intelligence (AI) is a new attempt at protein prediction, among which AlphaFold2 currently has more than 96% confidence in monomers. The AlphaFold2 method provides atomic-scale high-precision protein structural information, it also reduces the dependence on expensive equipment such as synchrotron radiation and cryo−EM. High-precision protein structure information is essential to identify the energy transfer process in light-harvesting complexes. Therefore, we conducted the relevant work and verified the possibility of AlphaFold2 for the energy transfer process. According to the structure information, two energy transfer paths are identified. This computing-based structure prediction method opens up a new avenue for understanding the energy transfer process in photosynthesis.
For the energy transfer process, AlphaFold2 is an exciting and promising technology and also showed satisfactory results from the overall effect. However, AlphaFold2 can be further improved in detail, when more precise energy transfer paths need to be considered. Moreover, the energy transfer process with a time constant of less than 1 ps, cannot be directly observed using the exiting technique. However, with the rapid development of the protein structure prediction technology and the ultrafast technology improvement (streak camera with time revolution of 100 fs or pump-probe method with time resolution less than 10 fs), future research on the physical mechanism of energy transfer in the initial photosynthesis process can dig more deeply into the energy localization kinetics of the coherence state of the excited states, which may bring us a more detailed understanding of the energy transfer pathways and rates.
The functions of proteins result from structures. Until now, nearly all the protein structures were determined through experimental methods such as crystallography, NMR and cryo−EM [11,33,34]. As a supplementary theoretical method, the AlphaFold model with high-quality predictions has great potential in protein structure determination which can be further used as the structural basis for energy transfer study in the light-harvesting complex. With the application of AlphaFold2, we may conduct more environmentally dependent analyses enriching the biological information library based on the cryo−electron microscopy Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/cryst12111595/s1, Figure S1: Alignment of the C−PC amino−acid sequences of the α−subunit and β−subunit from S. platensis., Table S1: The deconvolution results of C−PC hexamer.