Computational investigation of reaction mechanism of FET3 protein in yeast

Abstract

In Saccharomyces cerevisiae, a yeast species, iron uptake into the cell takes place with the Reducing Iron Uptake Model. Ferric chelates (Fe3C-L) are degraded on the cell surface by being reduced from Fe3+ to Fe2+ by the cell surface reductases Fre1p and Fre2p. The free reduced Fe2+ ions are taken up by Fet3p-Ftr1p, a high-affinity oxidase-permease complex, or by Fet4, another metal carrier. In this study, the reaction mechanism and the role of Fet3p in reducing iron uptake are examined. Fet3p is a membrane-bound protein and a member of the multicopper oxidase protein family. It metabolizes iron uptake with a high affinity for Fe2+ and plays a role in iron uptake together with iron-permease Ftr1p. Since Ftr1p can only transport the oxidized form of iron, Fe2+ needs to be oxidized before entering the cell. Fet3p couples the four-electron reduction of O2 to H2O with the one-electron oxidation of four Fe2+. The oxidized iron leaves the iron-binding site in Fet3p and is transferred to Ftr1p. Thus, Fe3+ ions are transported into the cytoplasm by a permease, Ftr1p. The understanding of the mechanism of Fet3p is of great importance to shed light on other multicopper oxidase members such as laccases and human ceruloplasmin, some having wide industrial applications. When the active site structure of Fet3p is examined, it has 4 copper as a cofactor in the active site. These coppers are divided into 3 types according to their characteristics: Type 1 (T1), Type 2 (T2), and binuclear Type 3 (T3a and T3b). T2 and T3 coppers form the trinuclear cluster (TNC). Iron as a substrate is not observed in any of the crystal structures of Fet3p. However, according to the information obtained from mutation studies and comparing them with the crystal structures of other MCOs, especially copper efflux oxidase (CueO), the amino acids in the iron-binding region of Fet3p and the location of iron were determined. Fet3p couples four one-electron oxidations of 4 Fe2+ as a substrate to the four-electron reduction of dioxygen to water by taking four protons from the environment. This process is mediated by oxidation-reduction reactions of copper ions as cofactors and consists of two stages. In the first stage, the O2 molecule, which will be reduced to H2O during the reaction, enters the TNC through the solvent channel and binds to the TNC. The O-O bond is cleaved by taking two electrons from two coppers (T1 and one of T3 coppers). Proton donation of E481 to one of the oxygens bridging T3 coppers facilitates this cleavage. Finally, all coppers are oxidized to Cu2+, and one O2- and two OH- ions are formed. In the second stage, the four reductions from Cu2+ to Cu+ with oxidation of four Fe2+ to Fe3+, and four protonations occur, and OH- and O2- ions are converted to two water molecules. In the literature, most of the first stage of the reaction mechanism of MCOs, especially dioxygen-cleavage and peroxy intermediate structure are known. However, the exact mechanism of the second stage, the order of electron and proton transfer reactions is not known because this part occurs fast. Due to the rapidity of these reactions, they have not been studied before and the order of the reaction is unknown due to the difficulty of following the protonation order experimentally. In addition to examining the reaction scheme, it is known that D283 plays an important role in iron binding to substrate-binding site, and electron transfer (ET) is enhanced by D283. However, in the crystal structure, the loop containing D283 is oriented away from the active site, suggesting that it closes only after the binding of Fe2+. Thus, to find out the role of D283 on ET and reaction pathways, the geometries are separately examined when the loop containing D283 is open and closed. In order to elucidate the unknown parts, computational methods were used in this present study, so the possible reaction mechanism will be determined. Thus, it is aimed to understand the mechanisms of other multicopper oxidase members through Fet3p. The calculations and geometry optimizations were carried out using the Quantum Mechanics/Molecular Mechanics QM/MM approach. The M06-2X method, a Density Functional Theory (DFT) method, was used for quantum mechanical (QM) calculations. B3LYP, TPSS, and M06 methods were also used to investigate whether M06-2X is the most suitable method for energy calculations and geometry optimizations of Fet3p containing copper and iron metals. Although M06-2X is not recommended to be used on metals in the literature, all necessary electronic states and spin densities could be obtained only with M06-2X in this study. For this reason, the results were interpreted over the energies obtained with M06-2X. The determination of the QM region to be calculated during the QM/MM calculations is of great importance for the calculations to obtain more accurate results. While choosing the most ideal QM region, residues that have the potential to affect the reaction, electron transfer, and proton exchange, especially close to the region where the reaction took place, were determined. The proximity of amino acids that will contribute to electron transfer around copper and Fe was investigated; therefore, calculations were made accordingly by choosing different QM regions. Considering the computational costs, the most ideal QM region was determined. In the structure where the loop containing D283 is closed, the first Fe2+ oxidation occurs exothermically without protonation while T1 is reduced. Protonation of OH- or O2- ions are not needed due to the cost of protonation. It is examined whether the first electron transfers from T1 to TNC before the second iron binds; nevertheless, the structure could not be obtained without protonation. With the protonation of the TNC region, electron transfer to the TNC has yielded a stable structure. After the oxidized iron leaves, the second Fe2+ binds. Meanwhile, the electron already transferred from the first iron remains in the protonated TNC. Considering the necessity of a second proton transfer before oxidation, the proton taken from D94 returns back to D94 during the optimization, thus the second proton transfer is not necessary. The electron from Fe2+ transfers to T1 copper, and oxidation of the second iron takes place. The second oxidation, which was endothermic when D283 was open, is exothermic in the structure where the loop is closed. The results draw attention to the importance of the loop containing D283. After the second oxidation, the oxidized Fe3+ is replaced with the third Fe2+. For the third iron, structures with two protons, three protons, and four protons are examined. For the third and fourth Fe2+ the geometries when the loop with D283 is open are also examined. According to the results, even three protonations are not enough for third oxidation, and a fourth protonation is needed. When the loop containing D283 is open, the oxidation of the fourth Fe2+ is endothermic even in the presence of four protons in TNC, which is the maximum number of protons TNC can take. In the oxidation reactions protonation of TNC-O2- decrease the negativity of TNC; thus, electron transfer to TNC is more favorable. The protonation of TNC is important to reduce coppers at TNC (T2 and T3 coppers) and for transferring an electron from the substrate to TNC. Similarly, the transfer of an electron from the substrate to TNC and the reduction of TNC coppers force the TNC-O2- or T3-OH- to take proton.