| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Center for Theoretical Biological Physics and Department of Physics, University of California at San Diego, La Jolla, California 92093–0374
Correspondence: Address reprint requests to José Nelson Onuchic, E-mail: jonuchic{at}ucsd.edu.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
; Chan and Li, 1990; Berg et al., 2002; Brzezinski et al., 1995; Winkler et al., 1995), a key bioenergetic reaction in aerobic organisms. Fig. 1 shows a schematic view of the key groups involved in ET. A copper dimer (CuA) functions as a primary donor, receiving electrons from cytochrome c molecules. The ET reaction is believed to occur in two subsequent steps: from CuA to heme a, and further to the complex of heme a3 and the molecular oxygen. A proton transfer reaction concomitant to the ET reduces the molecular oxygen to water, releasing the redox energy. This reaction creates and maintains a transmembrane pH gradient, which provides energy for adenosine triphosphate (ATP) synthesis. The produced ATP is then delivered to different locations throughout the cell, providing the necessary energy supply to a broad variety of biochemical reactions (Babcock and Wikstrom, 1992; Chan and Li, 1990; Berg et al., 2002; Brzezinski et al., 1995; Winkler et al., 1995).
|
, Tsukihara et al. (1995, 1996), and Yoshikawa et al. (1998)). Although a number of ET pathways has been suggested, there still is an ongoing discussion over which pathways dominate the ET reaction (Iwata et al., 1995; Ramirez et al., 1995; Williams et al., 1997; Regan et al., 1998; Gamelin et al., 1998; Paula et al., 1999; Farver et al, 2000; Medvedev et al., 2000). In this manuscript, we elucidate the ET pathways in CCO and investigate how thermal motion of the protein atoms affects the interference among the dominant pathway tubes (Regan and Onuchic, 1999). In the conventional theories of biological ET reactions (the Marcus theory (Marcus and Sutin, 1985) and its modifications), nuclear dynamics controls the Franck-Condon factor; however, it was recently shown that nuclear dynamics may also control the effective electronic coupling (Balabin and Onuchic, 2000; Beratan et al, 1987). |
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
, we calculate the effective electronic couplings between the donors and the acceptors using the Green's function technique with an extended Huckel-like electronic Hamiltonian. To study effects of nuclear motion on the electronic coupling, we perform molecular dynamics (MD) simulations of the protein in a native (membrane and solvent) environment.
Molecular dynamics simulations
A CCO structure from bovine heart in the fully oxidized state (Yoshikawa et al., 1998; Protein Data Bank code 2OCC) was simulated. The protein was embedded into a rectangular patch of phosphatidylethanolamine membrane. (See http://www.ks.uiuc.edu/~ilya/Membranes for details of the membrane structure.) The lipid molecules that overlapped with the proteins were removed. The protein-lipid complex was then solvated in a rectangular volume of preequilibrated TIP3 water molecules. Sodium and chlorine ions were added, corresponding to an ionic strength of 0.025 M/L. The final system measured 116 x 85 x 133 Å in size and included 122,084 atoms.
The simulations were performed using a molecular dynamics program NAMD (Kale et al., 1999) with the all-atom Charmm27 force field (MacKerell et al., 1998), particle mesh Ewald (PME) full electrostatics, periodic boundaries, constant pressure (NpT) ensemble, and a Langevin thermostat. In all simulations, the integrator time step was 1 fs. First, with all protein backbone atoms restrained, the system underwent 1,000 conjugate gradient minimization steps followed by 200 ps of equilibration at 310 K. Then the system was equilibrated for 30 ps with no restraints. The root mean square deviation (RMSD) of the protein C
atoms from the crystallographic conformation, along with the system volume, energies, temperature, and membrane surface area were monitored and did not show any significant changes after 30 ps. After that, the system was simulated for another 150 ps, and a trajectory frame was saved each 1 ps, providing 150 "snapshots" that described the protein dynamics. Thus only local equilibrium fluctuations are considered in this study. An equilibrium conformation was obtained by aligning and averaging positions of each atom in these "snapshots". The RMSD between the crystallographic structure and the equilibrium conformation, which describes the overall change in the protein conformation caused by the thermal motion, is 
1.5 Å. The RMSD between the equilibrium conformation and any of the "snapshots", which characterizes the average magnitude of each atom's thermal vibrations around its equilibrium position, is 0.6 Å. To further investigate conformational changes in protein associated with the thermal motion, RMSDs for each protein residue involved in mediating ET were computed. These RMSDs vary between 1 and 3 Å, depending on the side chain size rather than on the residue location in the protein; therefore, no major conformational changes in the protein are indicated.
Electron transfer reaction rate
Because of the typically large donor to acceptor distance (10–25 Å), most biological ET reactions occur in the weak coupling limit, and their rates are described by the Fermi golden rule (Marcus and Sutin, 1985; Hopfield, 1974; Jortner, 1980; Regan and Onuchic, 1999; Gray and Winkler, 1996; Page et al., 1999; Bertini et al., 1994; Sigel and Sigel., 1991)
 | (1) |
):
 | (2) |
is the reorganization energy, and 
G0 is the driving force of the reaction. In this article, we focus on qualitative effects of the protein structure and dynamics on the effective tunneling coupling, and we do not aim at obtaining quantitative estimates for the ET reaction rate. As such, we hereafter compute the maximum ET reaction rate, which occurs at the activationless limit ( + G0 
0), and the Franck-Condon factor is optimized. Although this assumption may lead to overestimating the reaction rate by up to three orders of magnitude, it does not interfere with the analysis of the effective coupling. Previous studies reported reorganization energies of about 0.3 eV for the reaction step from CuA to heme a (Regan et al., 1998; Randall et al., 2000; Brzezinski, 1996) and about 0.76 eV for the reaction step from heme a to heme a3 (Regan et al., 1998; Brzezinski, 1996; Adelroth et al., 1995). Assuming that the driving force is small, the maximum ET reaction rate is about 2 x 1016 (TDA)2 s-1, whereas the actual rate (computed with nonoptimized Franck-Condon factor) is 2 x 1015 (TDA)2 s-1 for the first step and 2 x 1013 (TDA)2 s-1 for the second step, respectively.
Effective tunneling matrix element calculations
Following the perturbation approach by Marcus (Marcus and Sutin, 1985), the effective tunneling matrix element for a given protein conformation can be calculated as follows:
 | (3) |
is the transformed Green's function of the bridge (Balabin and Onuchic, 1996)
 | (4) |
The overlap matrix was computed using a minimum basis set (four valence orbitals per heavy atom and one orbital per hydrogen atom) with the Slater type atomic orbitals. The bridge orbitals included all protein residues that participate in mediating the relevant electron transfer pathways. The relevant protein residues were identified in a way similar to the method described in Siddarth and Marcus (1990). Calculations of TDA between CuA and heme a included residues ArgA438, ArgA439, CysB196, CysB200, GluB198, HisA61, HisB161, HisB204, HisA378, IleB199, MetB207, SerB197, and TyrA440. Calculations of TDA between heme a and heme a3 included residues AlaA375, ArgA438, ArgA439, HisA61, HisA240, HisA290, HisA291, HisA376, HisA378, PheA377, TrpA126, TyrA440, and ValA374. Calculations of TDA between CuA and heme a3 included all of the above residues.
The Hamiltonian matrix was computed using a standard extended Huckel parametrization (Yates, 1978): the diagonal matrix elements were the ionization potentials of the respective atom, Hii = i, and the off-diagonal elements were assigned as Hij = (K/2) (i + j) Sij, where K = 1.5 was the standard Huckel constant (Yates, 1978). The overlap and Hamiltonian matrices were converted into the bonding-antibonding orbital basis using the standard transformations (Naray-Szabo et al, 1987). Although calculations of the effective electronic coupling with the above Hamiltonian cannot provide quantitative accuracy, they have proved, in the past experience of us (Balabin and Onuchic, 1996, 2000) and others (Regan and Onuchic, 1999; Medvedev et al, 2000), to provide reasonable estimates and to capture effects of the protein structure and dynamics.
Nuclear dynamics effects on the tunneling matrix element
Thermal motion of the protein atoms influences the overlap and the Hamiltonian matrices, therefore modulating the effective electronic coupling TDA. To investigate how this modulation affects the ET reaction rate, one needs to calculate the average square effective coupling (see article (Balabin and Onuchic, 2000) and the supporting information therein):
 | (5) |
is an effective coupling computed for the i-th protein "snapshot" conformation. The dynamical effects can be quantified by means of the following coherence parameter (Balabin and Onuchic, 2000): 
. When the interference among the ET pathways is predominantly constructive, TDA is robust to small changes of the protein conformation, and C is close to 1. In the opposite limit, when there is a pronounced destructive interference among the pathways, TDA is sensitive to details of the protein conformation, and C is close to 0. Therefore, C quantitatively describes the sensitivity of the effective coupling to the thermal motion of the protein atoms, characterizing the regime of the interference among the ET pathways (Balabin and Onuchic, 2000). To separate effects of the overall protein conformation change (the equilibrium conformation at 310 K versus the crystallographic conformation) from those of the thermal vibrations around the equilibrium, we computed the effective couplings for the crystallographic and the equilibrium conformations (Eq. 3) as well as the dynamic coupling (Eq. 5).
Ligand-metal bonds
ET from CuA is mediated by the copper ligands, which form metal bonds with the copper atoms. The coupling provided by these bonds (HDi in Eq. 3) is usually substantially weaker than the one provided by typical covalent bonds. In our calculations, these couplings are estimated using the metal bond covalencies (Gamelin et al., 1998; Williams et al., 1997; Solomon et al., 2000; Randall et al., 2000):
 | (6) |
), and Ci are the respective covalencies. Following the recent experimental data and ab initio calculation results (Gamelin et al., 1998; Williams et al., 1997; Solomon et al., 2000; Randall et al., 2000; George et al., 2001), we used the covalency values of 0.23 for the cysteins, 0.04 for the the histidine, and 0.01 for the glutamine.
Pathways model analysis
The Pathways model (Beratan et al., 1991; Betts et al., 1992; Curry et al., 1995) has provided a simple but efficient way to estimate the effective electronic coupling in proteins. It assumes that the effective coupling between the donor and the acceptor is proportional to a product of all partial decays in the electronic path (Beratan et al., 1991; Betts et al., 1992):
 | (7) |
is a decay per covalent bond mediated step, 
is a decay per hydrogen bond mediated step (R is the distance between the heavy atoms, Å), and 
is the decay per through-space jump (R is the jump distance, Å). Consistently with earlier estimates (Regan et al., 1999; Priyadarshy et al., 1998; Curry et al., 1995; Page et al., 1999), we assume that the maximum rate is 1013 s-1 at the donor to acceptor separation of 3 Å (two covalent bond steps), leading to the following approximation that relates the pathway structure with the maximum (not adjusted for the Franck-Condon factor) ET reaction rate:
 | (8) |
The Pathways model estimates were additionally verified by computing the extended coupling at the extended Huckel level for isolated individual pathways, one pathway at a time. Only atoms that mediate the pathway were included in the calculations. Effects of protein structural groups on the effective coupling were investigated by computing the effective coupling for the core path alone (no side groups), the core path with its nearest neighbor (NN) atoms (chemically bound to the core path atoms), and the core path with the NN and next to NN (chemically bound to the NN atoms) atoms, as well as with the full pathway structure (Balabin and Onuchic, 1998).
|
RESULTS AND DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
and 
orbitals (from -11 to -9 eV) (Balabin and Onuchic, 2000). It has been shown that the accuracy of extended Huckel level calculations of TDA can be improved by adjusting the tunneling energy based on independent spectral data (Siddarth and Marcus, 1990). In this article, however, we focus on the sensitivity of TDA to the tunneling energy and perform the calculations for a relevant tunneling energy range and not for a single, even though adjusted, tunneling energy.
|
orbital level (about -9 eV) (Balabin and Onuchic, 2000). At these energies, the dynamical coupling (see Eq. 5) is close to the one for the equilibrium conformation (Eq. 3), and neither one significantly depends on Etun. The effective coupling for the crystallographic conformation is smaller than either of the above couplings by at least an order of magnitude, and it does not significantly depend on Etun either. The coherence parameter is considerably large at all tunneling energies, indicating that the interference among the dominant pathways is constructive, and the effective coupling is robust to the thermal motion of the protein atoms. The major effect of the protein dynamics is shortening of some through-space jump distances, which increases the effective coupling (Balabin and Onuchic, 2000). The dynamical coupling translates into the maximum ET rate of 6 x 104 s-1, which is in a close agreement with the experimentally observed rate of 2 x 104 s-1 (Farver et al., 2000; Brzezinski, 1996; Adelroth et al., 1995). Notably, the much smaller maximum ET rate for the crystallographic structure illustrates how calculations based only on the crystallographic conformation can be grossly misleading (Balabin and Onuchic, 2000).
Protein groups between CuA and heme a mediates several possible pathways from CuA to heme a, which start on different copper ligands (Fig. 1). Two pathways that start on HisB204 and CysB196, respectively, have been proposed before (Ramirez et al., 1995; Williams et al., 1997). Using the Pathways analysis, we identified two other potentially relevant pathways that start on CysB200 and GluB198, respectively. (Independently of our work, the CysB200 pathway was suggested in Medvedev et al. (2000)). The effective coupling mediated by each of the above pathways and the corresponding maximum rate computed for the crystallographic and the equilibrium conformations using the Pathways model are presented in Table 1.
|
0.5 Å, increasing the maximum ET rate via the pathway that starts on HisB204 by almost an order of magnitude. Conversely, the distance of the through-space jump between IleB199 and ArgA438 increases in the equilibrium conformation by >2 Å, substantially decreasing the contribution of the pathway that starts on GluB198 to the effective coupling. The contributions of the other two pathways remain small.
Electron transfer from heme a to heme a3
Fig. 3 shows the effective coupling between heme a and heme a3 calculated at the extended Huckel level. At all energies, the dynamic effective coupling is close to that for the equilibrium conformation, neither one significantly depending on Etun. The effective coupling for the crystallographic conformation approaches the above couplings at the most relevant tunneling energies close to -9 eV. Thus, the pathway interference regime is constructive, and the effective coupling is robust to the thermal motions and the tunneling energy.
|
3 x 108 s-1, substantially exceeding the experimentally observed rate of 3 x 105 s-1 (Brzezinski, 1996; Adelroth et al., 1995; Verkhovsky et al., 1992). This difference may be attributed to a nonoptimized Franck-Condon factor (Eq. 2) (Regan et al., 1998). Recently, Verkhovsky and others (Verkhovsky et al, 2001) suggested that ET from heme a to heme a3 includes two kinetically different components: the ultrafast (108 s-1) transfer, which is energetically optimized, and the regular transfer, which is controlled by either heme-protein relaxation or CO dissociation from CuB. Although this conclusion remains controversial (Namslauer et al., 2002), our calculations indirectly support the possibility of the ultrafast ET component.
Protein groups surrounding heme a and heme a3 mediate a variety of different pathways (Fig. 1). Two pathways that have been suggested before include the fully covalent pathway, which is mediated by HisA378, PheA377, and HisA376 (Regan et al., 1998), and the direct pathway that includes a through-space jump between the 
propionates of the two hemes (Medvedev et al., 2000). Using the Pathways analysis, we identified two other relevant pathways: one mediated by ArgA438, which includes two hydrogen bonds, and another mediated by HisA378, AlaA375, and HisA376, which includes one hydrogen bond. The numbers of steps and the corresponding maximum rates for these pathways are shown in Table 2. In the crystallographic structure, the covalent pathway and the direct pathway interfere constructively, dominating the effective coupling. Thermal motion shortens the distance of the through-space jump (CMD-HMD[123] groups in either heme) to 2 Å (Fig. 4) (and possibly less, on longer timescales), increasing the corresponding ET reaction rate by an order of magnitude and making the direct pathway dominant (Fig. 3). Recently, it was proposed (Zheng et al., 2003) that the ET from heme a to heme a3 is facilitated by water molecules inside the protein. Although inner water is certainly a possibility in CCO, its influence on the ET between the two hemes appears unlikely. As shown in Fig. 5, the protein groups surrounding the dominant direct pathway are strongly hydrophobic; in addition, distances between the two hemes do not leave enough space for water molecules.
|
|
|
Direct or two-step reaction?
A possibility of the direct ET reaction from CuA to a3 has been brought on by Tsukihara and others (Tsukihara et al., 1995). The reaction could be mediated by four pathways, which start on the copper ligands CysB196, GluB198, CysB200, and HisB204. The effective coupling and the corresponding maximum rate via any of these pathways were estimated at the Pathways level. The ligand-metal bonds were treated as covalent bonds, and the reorganization energy was set to 0.8 eV. For both the crystallographic and equilibrium conformations, the maximum rate for the direct reaction turned out to be >2 orders of magnitude smaller than the rate of the two-step reaction. Even accounting for possible inner water molecules, which could facilitate the effective coupling for the direct ET reaction, the latter does not appear to be relevant.
|
CONCLUSIONS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
). |
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
This work was supported by the National Science Foundation (grant MCB-0084797 with additional support from PHY-0216576, 0225630, and 0107594). Any opinions, findings, and conclusions or recommendations expressed in this publication are those of the authors and do not necessarily reflect the views of the National Science Foundation.
|
FOOTNOTES |
|---|
Ilya Balabin's present address is P. M. Gross Chemical Laboratory, Duke University, Box 90349, Durham, NC 27708.
Submitted on July 18, 2003; accepted for publication October 21, 2003.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONCONCLUSIONSACKNOWLEDGEMENTSREFERENCES |
|---|
Babcock, G. T., and M. Wikstrom. 1992. Oxygen activation and the conservation of energy in cell respiration. Nature. 356:301–309.[Medline]
Balabin, I. A., and J. N. Onuchic. 1996. Connection between simple models and quantum chemical models for electron-transfer tunneling matrix element calculations: a Dyson's equations-based approach. J. Phys. Chem. 100:11573–11580.
Balabin, I. A., and J. N. Onuchic. 1998. A new framework for electron-transfer calculations: beyond the Pathways-like models. J. Phys. Chem. 102:7497–7505.
Balabin, I. A., and J. N. Onuchic. 2000. Dynamically controlled protein tunneling paths in photosynthetic reaction centers. Science. 290:114–117.
Beratan, D. N., J. N. Onuchic, and J. J. Hopfield. 1987. Electron tunneling through covalent and noncovalent pathways in proteins. J. Chem. Phys. 86:4488–4498.
Beratan, D. N., J. N. Betts, and J. N. Onuchic. 1991. Protein electron transfer rates set by the bridging secondary and tertiary structure. Science. 252:1285–1288.
Berg, J. M., J. L. Tymoczko, and L. Stryer. 2002. Biochemistry. W. H. Freeman, New York, NY.
Bertini, I., H. B. Gray, S. Lippard, and J. S. Valentine. 1994. Bioinorganic Chemistry. University Science Books, Mill Valley, CA .
Betts, J. N., D. N. Beratan, and J. N. Onuchic. 1992. Mapping electron tunneling pathways: an algorithm that finds the minimum length maximum coupling pathway between electron donors and acceptors in proteins. J. Am. Chem. Soc. 114:4043–4046.
Brzezinski, P., M. Sundahl, P. Adelroth, M. T. Wilson,, B. el-Agez, P. Wittung, and B. G. Malmstrom. 1995. Triplet-state quenchng in complexes between Zn-cytochrome-C and cytochrome-oxidase or its Cu-A domain. Biophys. Chem. 54:191–197.[Medline]
Brzezinski, P. 1996. Internal electron-transfer reactions in cytochrome c oxidase. Biochemistry. 35:5611–5615.[Medline]
Chan, S. I., and P. M. Li. 1990. Cytochrome c oxidase: understanding nature's design of a proton pump. Biochemistry. 29:1–12.[Medline]
Curry, W. B., M. D. Grabe, I. V. Kurnikov, S. S. Skourtis, D. N. Beratan, J. J. Regan, A. J. A. Aquino, P. Beroza, and J. N. Onuchic. 1995. Pathways, pathway tubes, pathway docking, and propagators in electron - transfer proteins. J. Bioenerg. Biomembr. 27:285–293.[Medline]
Farver, O., O. Einarsdottir, and I. Pecht. 2000. Electron transfer rates and equilibrium within cytochrome c oxidase. Eur. J. Biochem. 267:950–954.[Medline]
Gamelin, D. R., D. W. Randall, M. T. Hay, R. P. Houser, T. C. Mulder, G. W. Canters, S. de Vries, W. B. Tolman, Y. Lu, and E. I. Solomon. 1998. Spectroscopy of mixed-valence Cu-A-type centers: ligand-field control of ground-state properties related to electron transfer. J. Am. Chem. Soc. 120:5246–5263.
George, S. D., M. Metz, R. K. Szilagyi, H. X. Wang, S. P. Cramer, Y. Lu, W. B. Tolman, B. Hedman, K. O. Hodgson, and E. I. Solomon. 2001. A quantitative description of the ground-state wave function of Cu-A by x-ray absorption spectroscopy: comparison to plastocyanin and relevance to electron transfer. J. Am. Chem. Soc. 123:5757–5767.[Medline]
Gray, H. B., and J. R. Winkler. 1996. Electron transfer in proteins. Annu. Rev. Biochem. 65:537–561.[Medline]
Hopfield, J. J. 1974. Electron transfer between biological molecules by thermally activated tunneling. Proc. Natl. Acad. Sci. USA. 71:3640–3644.
Iwata, S., C. Ostermeier, B. Ludwig, and H. Michel. 1995. Structure at 2.8-angstrom resolution of cytochrome c oxidase from paracoccus denitrificans. Nature. 376:660–669.[Medline]
Jortner, J. 1980. Dynamics of electron transfer in bacterial photosynthesis. Biochim. Biophys. Acta. 594:193–230.[Medline]
Kale, L., R. Skeel, M. Bhandarkar, R. Brunner, A. Gursoy, N. Krawetz, J. Phillips, A. Shinozaki, K. Varadarajan, and K. Schulten. 1999. NAMD2: greater scalability for parallel molecular dynamics. J. Comput. Phys. 151:283–312.
MacKerell, A. D., Jr., D. Bashford, M. Bellott, R. L. Dunbrack, J. D. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph-McCarthy, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, W. E. Reiher, B. Roux, M. Schlenkrich, J. C. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, D. Yin, and M. Karplus. 1998. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B. 102:3586–3616.
Marcus, R. A., and N. Sutin. 1985. Electron transfers in chemistry and biology. Biochim. Biophys. Acta. 811:265–322.
Medvedev, D. M., I. Daizadeh, and A. A. Stuchebrukhov. 2000. Electron transfer tunneling pathways in bovine heart cytochrome c oxidase. J. Am. Chem. Soc. 122:6571–6582.
Namslauer, A., M. Branden, and P. Brzezinski. 2002. The rate of internal heme-heme electron transfer in cytochrome c oxidase. Biochemistry. 41:10369–10374.[Medline]
Naray-Szabo, G., P. R. Surjan, and J. Angyan. 1987. Applied Quantum Chemistry. Dordrecht, Boston, MA.
Page, C. C., C. C. Moser, X. Chen, and P. L. Dutton. 1999. Natural engineering principles of electron tunnelling in biological oxidation-reduction. Nature. 402:47–52.[Medline]
Paula, S., A. Sucheta, I. Szundi, and O. Einarsdottir. 1999. Proton and electron transfer during the reduction of molecular oxygen by fully reduced cytochrome c oxidase: a flow-flash investigation using optical multichannel detection. Biochemistry. 38:3025–3033.[Medline]
Priyadarshy, S., S. M. Risser, and D. N. Beratan. 1998. DNA-mediated electron transfer. JBIC. 3:196–200.
Priyadarshy, S., S. S. Skourtis, S. M. Risser, and D. N. Beratan. 1996. Bridge-mediated electronic interactions: differences between Hamiltonian and Green function partitioning in a non-orthogonal basis. J. Chem. Phys. 104:9473–9481.
Ramirez, B. E., B. G. Malmstrom, J. R. Winkler, and H. B. Gray. 1995. The currents of life: the terminal electron-transfer complex of respiration. Proc. Natl. Acad. Sci. USA. 92:11949–11951.
Randall, D. W., D. R. Gamelin, L. B. LaCroix, and E. I. Solomon. 2000. Electronic structure contributions to electron transfer in blue Cu and Cu-A. JBIC. 5:16–29.[Medline]
Regan, J. J., B. E. Ramirez, J. R. Winkler, H. B. Gray, and B. G. Malmstrom. 1998. Pathways for electron tunneling in cytochrome c oxidase. J. Bioenerg. Biomembr. 30:35–39.[Medline]
Regan, J. J., and J. N. Onuchic. 1999. Electron-transfer tubes. Adv. Chem. Phys. 107:497–553.
Siddarth, P., and R. A. Marcus. 1990. Electron-transfer reactions in proteins: a calculation of electronic coupling. J. Phys. Chem. 94:8430–8434.
Sigel, H., and A. Sigel. 1991. Electron Transfer Reactions in Metalloproteins. H. Sigel and A. Sigel, editors. Marcel Dekker, New York, NY.
Solomon, E. I., D. W. Randall, and T. Glaser. 2000. Electronic structures of active sites in electron transfer metalloproteins: contributions to reactivity. Coord. Chem. Rev. 200:595–632.
Tsukihara, T., H. Aoyama, E. Yamashita, T. Tomizaki, H. Yamaguchi, K. Shinzawa-Itoh, R. Nakashima, R. Yaono, and S. Yoshikawa. 1995. Structures of metal sites of oxidized bovine heart cytochrome c oxidase at 2.8 angstrom. Science. 269:1069–1074.
Tsukihara, T., H. Aoyama, E. Yamashita, T. Tomizaki, H. Yamaguchi, K. Shinzawa-Itoh, R. Nakashima, R. Yaono, and S. Yoshikawa. 1996. The whole structure of the 13-subunit oxidized cytochrome c oxidase at 2.8 angstrom. Science. 272:1136–1144.[Abstract]
Verkhovsky, M. I., J. E. Morgan, and M. Wikstrom. 1992. Intramolecular electron transfer in cytochrome c oxidase: a cascade of equilibria. Biochemistry. 31:11860–11863.[Medline]
Verkhovsky, M. I., A. Jasaitis, and M. Wikstrom. 2001. Ultrafast haem-haem electron transfer in cytochrome c oxidase. Biochim. Biophys. Acta. 1506:143–146.[Medline]
Williams, K. R., D. R. Gamelin, L. B. LaCroix, R. P. Houser, W. B. Tolman, T. C. Mulder, S. de Vries, B. Hedman, K. O. Hodgson, and E. I. Solomon. 1997. Influence of copper-sulfur covalency and copper-copper bonding on valence delocalization and electron transfer in the Cu-A site of cytochrome c oxidase. J. Am. Chem. Soc. 119:613–614.
Winkler, J. R., B. G. Malmstrom, and H. B. Gray. 1995. Rapid electron injection into multisite metalloproteins: intramolecular electron-transfer in cytochrome-oxidase. Biophys. Chem. 54:199–209.[Medline]
Yates, K. 1978. Huckel Molecular Orbital Theory. Academic Press, New York, NY.
Yoshikawa, S., K. Shinzawa-Itoh, R. Nakashima, R. Yaono, E. Yamashita, N. Inoue, M. Yao, M. J. Fei, C. P. Libeu, T. Mizushima, H. Yamaguchi, T. Tomizaki, and T. Tsukihara. 1998. Redox-coupled crystal structural changes in bovine heart cytochrome c oxidase. Science. 280:1723–1729.
Zheng, X. H., D. M. Medvedev, J. Swanson, and A. A. Stuchebrukhov. 2003. Computer simulation of water in cytochrome c oxidase. Biochim. Biophys. Acta Bioenerg. 1557:99–107.
This article has been cited by other articles:
 |
|
 |
|
  D. N. Beratan and I. A. Balabin Heme copper oxidases use tunneling pathways PNAS, January 15, 2008; 105(2): 403 - 404. [Full Text] [PDF]
|
|
|
|
|
|
A. Jasaitis, M. P. Johansson, M. Wikstrom, M. H. Vos, and M. I. Verkhovsky Nanosecond electron tunneling between the hemes in cytochrome bo3 PNAS, December 26, 2007; 104(52): 20811 - 20814. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. H. M. Olsson and A. Warshel Monte Carlo simulations of proton pumps: On the working principles of the biological valve that controls proton pumping in cytochrome c oxidase PNAS, April 25, 2006; 103(17): 6500 - 6505. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Lin, I. A. Balabin, and D. N. Beratan The Nature of Aqueous Tunneling Pathways Between Electron-Transfer Proteins Science, November 25, 2005; 310(5752): 1311 - 1313. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Jasaitis, F. Rappaport, E. Pilet, U. Liebl, and M. H. Vos Activationless electron transfer through the hydrophobic core of cytochrome c oxidase PNAS, August 2, 2005; 102(31): 10882 - 10886. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
O. Miyashita, M. Y. Okamura, and J. N. Onuchic Long-Range Electron Transfer Special Feature: Interprotein electron transfer from cytochrome c2 to photosynthetic reaction center: Tunneling across an aqueous interface PNAS, March 8, 2005; 102(10): 3558 - 3563. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
E. Pilet, A. Jasaitis, U. Liebl, and M. H. Vos Electron transfer between hemes in mammalian cytochrome c oxidase PNAS, November 16, 2004; 101(46): 16198 - 16203. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
