1
|
Biswas S, Khaing EP, Zhong V, Eaton-Rye JJ. Arg24 and 26 of the D2 protein are important for photosystem II assembly and plastoquinol exchange in Synechocystis sp. PCC 6803. BIOCHIMICA ET BIOPHYSICA ACTA. BIOENERGETICS 2024; 1865:149150. [PMID: 38906313 DOI: 10.1016/j.bbabio.2024.149150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 05/26/2024] [Accepted: 06/13/2024] [Indexed: 06/23/2024]
Abstract
Photosystem II (PS II) assembly is a stepwise process involving preassembly complexes or modules focused around four core PS II proteins. The current model of PS II assembly in cyanobacteria is derived from studies involving the deletion of one or more of these core subunits. Such deletions may destabilize other PS II assembly intermediates, making constructing a clear picture of the intermediate events difficult. Information on plastoquinone exchange pathways operating within PS II is also unclear and relies heavily on computer-aided simulations. Deletion of PsbX in [S. Biswas, J.J. Eaton-Rye, Biochim. Biophys. Acta - Bioenerg. 1863 (2022) 148519] suggested modified QB binding in PS II lacking this subunit. This study has indicated the phenotype of the ∆PsbX mutant arose by disrupting a conserved hydrogen bond between PsbX and the D2 (PsbD) protein. We mutated two conserved arginine residues (D2:Arg24 and D2:Arg26) to further understand the observations made with the ∆PsbX mutant. Mutating Arg24 disrupted the interaction between PsbX and D2, replicating the high-light sensitivity and altered fluorescence decay kinetics observed in the ∆PsbX strain. The Arg26 residue, on the other hand, was more important for either PS II assembly or for stabilizing the fully assembled complex. The effects of mutating both arginine residues to alanine or aspartate were severe enough to render the corresponding double mutants non-photoautotrophic. Our study furthers our knowledge of the amino-acid interactions stabilizing plastoquinone-exchange pathways while providing a platform to study PS II assembly and repair without the actual deletion of any proteins.
Collapse
Affiliation(s)
- Sandeep Biswas
- Department of Biochemistry, University of Otago, Dunedin 9054, New Zealand
| | - Ei Phyo Khaing
- Department of Biochemistry, University of Otago, Dunedin 9054, New Zealand
| | - Victor Zhong
- Department of Biochemistry, University of Otago, Dunedin 9054, New Zealand
| | - Julian J Eaton-Rye
- Department of Biochemistry, University of Otago, Dunedin 9054, New Zealand.
| |
Collapse
|
2
|
Fatima S, Olshansky L. Conformational control over proton-coupled electron transfer in metalloenzymes. Nat Rev Chem 2024; 8:762-775. [PMID: 39223400 DOI: 10.1038/s41570-024-00646-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/29/2024] [Indexed: 09/04/2024]
Abstract
From the reduction of dinitrogen to the oxidation of water, the chemical transformations catalysed by metalloenzymes underlie global geochemical and biochemical cycles. These reactions represent some of the most kinetically and thermodynamically challenging processes known and require the complex choreography of the fundamental building blocks of nature, electrons and protons, to be carried out with utmost precision and accuracy. The rate-determining step of catalysis in many metalloenzymes consists of a protein structural rearrangement, suggesting that nature has evolved to leverage macroscopic changes in protein molecular structure to control subatomic changes in metallocofactor electronic structure. The proton-coupled electron transfer mechanisms operative in nitrogenase, photosystem II and ribonucleotide reductase exemplify this interplay between molecular and electronic structural control. We present the culmination of decades of study on each of these systems and clarify what is known regarding the interplay between structural changes and functional outcomes in these metalloenzyme linchpins.
Collapse
Affiliation(s)
- Saman Fatima
- Department of Chemistry, College of Liberal Arts and Sciences, University of Illinois Urbana-Champaign, Urbana, IL, USA
| | - Lisa Olshansky
- Department of Chemistry, College of Liberal Arts and Sciences, University of Illinois Urbana-Champaign, Urbana, IL, USA.
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, IL, USA.
- Materials Research Laboratory, The Grainger College of Engineering, University of Illinois Urbana-Champaign, Urbana, IL, USA.
- The Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
3
|
Guo Y, He L, Ding Y, Kloo L, Pantazis DA, Messinger J, Sun L. Closing Kok's cycle of nature's water oxidation catalysis. Nat Commun 2024; 15:5982. [PMID: 39013902 PMCID: PMC11252165 DOI: 10.1038/s41467-024-50210-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Accepted: 07/03/2024] [Indexed: 07/18/2024] Open
Abstract
The Mn4CaO5(6) cluster in photosystem II catalyzes water splitting through the Si state cycle (i = 0-4). Molecular O2 is formed and the natural catalyst is reset during the final S3 → (S4) → S0 transition. Only recently experimental breakthroughs have emerged for this transition but without explicit information on the S0-state reconstitution, thus the progression after O2 release remains elusive. In this report, our molecular dynamics simulations combined with density functional calculations suggest a likely missing link for closing the cycle, i.e., restoring the first catalytic state. Specifically, the formation of closed-cubane intermediates with all hexa-coordinate Mn is observed, which would undergo proton release, water dissociation, and ligand transfer to produce the open-cubane structure of the S0 state. Thereby, we theoretically identify the previously unknown structural isomerism in the S0 state that acts as the origin of the proposed structural flexibility prevailing in the cycle, which may be functionally important for nature's water oxidation catalysis.
Collapse
Affiliation(s)
- Yu Guo
- Center of Artificial Photosynthesis for Solar Fuels and Department of Chemistry, School of Science, Westlake University, Hangzhou, 310024, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, 310024, China
| | - Lanlan He
- Center of Artificial Photosynthesis for Solar Fuels and Department of Chemistry, School of Science, Westlake University, Hangzhou, 310024, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, 310024, China
| | - Yunxuan Ding
- Center of Artificial Photosynthesis for Solar Fuels and Department of Chemistry, School of Science, Westlake University, Hangzhou, 310024, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, 310024, China
| | - Lars Kloo
- Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-10044, Stockholm, Sweden
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr, 45470, Germany
| | - Johannes Messinger
- Department of Plant Physiology, Umeå University, Linnaeus väg 6 (KBC huset), SE-90187, Umeå, Sweden
- Molecular Biomimetics, Department of Chemistry - Ångström Laboratory, Uppsala University, SE-75120, Uppsala, Sweden
| | - Licheng Sun
- Center of Artificial Photosynthesis for Solar Fuels and Department of Chemistry, School of Science, Westlake University, Hangzhou, 310024, China.
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, 310024, China.
- Division of Solar Energy Conversion and Catalysis at Westlake University, Zhejiang Baima Lake Laboratory Co., Ltd., Hangzhou, 310000, Zhejiang, China.
| |
Collapse
|
4
|
Krysiak S, Burda K. The Effect of Removal of External Proteins PsbO, PsbP and PsbQ on Flash-Induced Molecular Oxygen Evolution and Its Biphasicity in Tobacco PSII. Curr Issues Mol Biol 2024; 46:7187-7218. [PMID: 39057069 PMCID: PMC11276211 DOI: 10.3390/cimb46070428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2024] [Revised: 06/30/2024] [Accepted: 07/02/2024] [Indexed: 07/28/2024] Open
Abstract
The oxygen evolution within photosystem II (PSII) is one of the most enigmatic processes occurring in nature. It is suggested that external proteins surrounding the oxygen-evolving complex (OEC) not only stabilize it and provide an appropriate ionic environment but also create water channels, which could be involved in triggering the ingress of water and the removal of O2 and protons outside the system. To investigate the influence of these proteins on the rate of oxygen release and the efficiency of OEC function, we developed a measurement protocol for the direct measurement of the kinetics of oxygen release from PSII using a Joliot-type electrode. PSII-enriched tobacco thylakoids were used in the experiments. The results revealed the existence of slow and fast modes of oxygen evolution. This observation is model-independent and requires no specific assumptions about the initial distribution of the OEC states. The gradual removal of exogenous proteins resulted in a slowdown of the rapid phase (~ms) of O2 release and its gradual disappearance while the slow phase (~tens of ms) accelerated. The role of external proteins in regulating the biphasicity and efficiency of oxygen release is discussed based on observed phenomena and current knowledge.
Collapse
Affiliation(s)
| | - Kvetoslava Burda
- Faculty of Physics and Applied Computer Science, AGH University of Krakow, al. Mickiewicza 30, 30-059 Krakow, Poland;
| |
Collapse
|
5
|
Flesher DA, Liu J, Wang J, Gisriel CJ, Yang KR, Batista VS, Debus RJ, Brudvig GW. Mutation-induced shift of the photosystem II active site reveals insight into conserved water channels. J Biol Chem 2024; 300:107475. [PMID: 38879008 PMCID: PMC11294709 DOI: 10.1016/j.jbc.2024.107475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2024] [Revised: 06/02/2024] [Accepted: 06/09/2024] [Indexed: 07/11/2024] Open
Abstract
Photosystem II (PSII) is the water-plastoquinone photo-oxidoreductase central to oxygenic photosynthesis. PSII has been extensively studied for its ability to catalyze light-driven water oxidation at a Mn4CaO5 cluster called the oxygen-evolving complex (OEC). Despite these efforts, the complete reaction mechanism for water oxidation by PSII is still heavily debated. Previous mutagenesis studies have investigated the roles of conserved amino acids, but these studies have lacked a direct structural basis that would allow for a more meaningful interpretation. Here, we report a 2.14-Å resolution cryo-EM structure of a PSII complex containing the substitution Asp170Glu on the D1 subunit. This mutation directly perturbs a bridging carboxylate ligand of the OEC, which alters the spectroscopic properties of the OEC without fully abolishing water oxidation. The structure reveals that the mutation shifts the position of the OEC within the active site without markedly distorting the Mn4CaO5 cluster metal-metal geometry, instead shifting the OEC as a rigid body. This shift disturbs the hydrogen-bonding network of structured waters near the OEC, causing disorder in the conserved water channels. This mutation-induced disorder appears consistent with previous FTIR spectroscopic data. We further show using quantum mechanics/molecular mechanics methods that the mutation-induced structural changes can affect the magnetic properties of the OEC by altering the axes of the Jahn-Teller distortion of the Mn(III) ion coordinated to D1-170. These results offer new perspectives on the conserved water channels, the rigid body property of the OEC, and the role of D1-Asp170 in the enzymatic water oxidation mechanism.
Collapse
Affiliation(s)
- David A Flesher
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut, USA
| | - Jinchan Liu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut, USA
| | - Jimin Wang
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut, USA
| | | | - Ke R Yang
- Department of Chemistry, Yale University, New Haven, Connecticut, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut, USA
| | - Richard J Debus
- Department of Biochemistry, University of California, Riverside, California, USA.
| | - Gary W Brudvig
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut, USA; Department of Chemistry, Yale University, New Haven, Connecticut, USA.
| |
Collapse
|
6
|
Hussein R, Graça A, Forsman J, Aydin AO, Hall M, Gaetcke J, Chernev P, Wendler P, Dobbek H, Messinger J, Zouni A, Schröder WP. Cryo-electron microscopy reveals hydrogen positions and water networks in photosystem II. Science 2024; 384:1349-1355. [PMID: 38900892 DOI: 10.1126/science.adn6541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Accepted: 05/16/2024] [Indexed: 06/22/2024]
Abstract
Photosystem II starts the photosynthetic electron transport chain that converts solar energy into chemical energy and thus sustains life on Earth. It catalyzes two chemical reactions: water oxidation to molecular oxygen and plastoquinone reduction. Coupling of electron and proton transfer is crucial for efficiency; however, the molecular basis of these processes remains speculative owing to uncertain water binding sites and the lack of experimentally determined hydrogen positions. We thus collected high-resolution cryo-electron microscopy data of fully hydrated photosystem II from the thermophilic cyanobacterium Thermosynechococcus vestitus to a final resolution of 1.71 angstroms. The structure reveals several previously undetected partially occupied water binding sites and more than half of the hydrogen and proton positions. This clarifies the pathways of substrate water binding and plastoquinone B protonation.
Collapse
Affiliation(s)
- Rana Hussein
- Humboldt-Universität zu Berlin, Department of Biology, D 10099 Berlin, Germany
| | - André Graça
- Department of Chemistry, Umeå University, SE 90187 Umeå, Sweden
- Molecular Biomimetics, Department of Chemistry- Ångström Laboratory, Uppsala University, SE 75120 Uppsala, Sweden
| | - Jack Forsman
- Department of Chemistry, Umeå University, SE 90187 Umeå, Sweden
| | - A Orkun Aydin
- Molecular Biomimetics, Department of Chemistry- Ångström Laboratory, Uppsala University, SE 75120 Uppsala, Sweden
| | - Michael Hall
- Department of Chemistry, Umeå University, SE 90187 Umeå, Sweden
| | - Julia Gaetcke
- Humboldt-Universität zu Berlin, Department of Biology, D 10099 Berlin, Germany
| | - Petko Chernev
- Molecular Biomimetics, Department of Chemistry- Ångström Laboratory, Uppsala University, SE 75120 Uppsala, Sweden
| | - Petra Wendler
- Institute of Biochemistry and Biology, Department of Biochemistry, University of Potsdam, Karl-Liebknecht Strasse 24-25, D 14476, Potsdam-Golm, Germany
| | - Holger Dobbek
- Humboldt-Universität zu Berlin, Department of Biology, D 10099 Berlin, Germany
| | - Johannes Messinger
- Molecular Biomimetics, Department of Chemistry- Ångström Laboratory, Uppsala University, SE 75120 Uppsala, Sweden
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, Sweden
| | - Athina Zouni
- Humboldt-Universität zu Berlin, Department of Biology, D 10099 Berlin, Germany
| | - Wolfgang P Schröder
- Department of Chemistry, Umeå University, SE 90187 Umeå, Sweden
- Umeå Plant Science Centre, Department of Plant Physiology, Umeå University, Sweden
| |
Collapse
|
7
|
Liu J, Yang KR, Long Z, Armstrong WH, Brudvig GW, Batista VS. Water Ligands Regulate the Redox Leveling Mechanism of the Oxygen-Evolving Complex of the Photosystem II. J Am Chem Soc 2024; 146:15986-15999. [PMID: 38833517 DOI: 10.1021/jacs.4c02926] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024]
Abstract
Understanding how water ligands regulate the conformational changes and functionality of the oxygen-evolving complex (OEC) in photosystem II (PSII) throughout the catalytic cycle of oxygen evolution remains a highly intriguing and unresolved challenge. In this study, we investigate the effect of water insertion (WI) on the redox state of the OEC by using the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) hybrid methods. We find that water binding significantly reduces the free energy change for proton-coupled electron transfer (PCET) from Mn to YZ•, underscoring the important regulatory role of water binding, which is essential for enabling the OEC redox-leveling mechanism along the catalytic cycle. We propose a water binding mechanism in which WI is thermodynamically favored by the closed-cubane form of the OEC, with water delivery mediated by Ca2+ ligand exchange. Isomerization from the closed- to open-cubane conformation at three post-WI states highlights the importance of the location of the MnIII center in the OEC and the orientation of its Jahn-Teller axis to conformational changes of the OEC, which might be critical for the formation of the O-O bond. These findings reveal a complex interplay between conformational changes in the OEC and the ligand environment during the activation of the OEC by YZ•. Analogous regulatory effects due to water ligand binding are expected to be important for a wide range of catalysts activated by redox state transitions in aqueous environments.
Collapse
Affiliation(s)
- Jinchan Liu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Ke R Yang
- Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), College of Chemistry, Nankai University, Tianjin 300071, China
| | - Zhuoran Long
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - William H Armstrong
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Gary W Brudvig
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
8
|
Bhattacharjee S, Arra S, Daidone I, Pantazis DA. Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations. Chem Sci 2024; 15:7269-7284. [PMID: 38756808 PMCID: PMC11095388 DOI: 10.1039/d3sc06714a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Accepted: 04/07/2024] [Indexed: 05/18/2024] Open
Abstract
Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis.
Collapse
Affiliation(s)
- Sinjini Bhattacharjee
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| | - Srilatha Arra
- Department of Physical and Chemical Sciences, University of L'Aquila Via Vetoio (Coppito 1) 67010 L'Aquila Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila Via Vetoio (Coppito 1) 67010 L'Aquila Italy
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| |
Collapse
|
9
|
Debus RJ, Oyala PH. Independent Mutation of Two Bridging Carboxylate Ligands Stabilizes Alternate Conformers of the Photosynthetic O 2-Evolving Mn 4CaO 5 Cluster in Photosystem II. J Phys Chem B 2024; 128:3870-3884. [PMID: 38602496 DOI: 10.1021/acs.jpcb.4c00829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
The O2-evolving Mn4CaO5 cluster in photosystem II is ligated by six carboxylate residues. One of these is D170 of the D1 subunit. This carboxylate bridges between one Mn ion (Mn4) and the Ca ion. A second carboxylate ligand is D342 of the D1 subunit. This carboxylate bridges between two Mn ions (Mn1 and Mn2). D170 and D342 are located on opposite sides of the Mn4CaO5 cluster. Recently, it was shown that the D170E mutation perturbs both the intricate networks of H-bonds that surround the Mn4CaO5 cluster and the equilibrium between different conformers of the cluster in two of its lower oxidation states, S1 and S2, while still supporting O2 evolution at approximately 50% the rate of the wild type. In this study, we show that the D342E mutation produces much the same alterations to the cluster's FTIR and EPR spectra as D170E, while still supporting O2 evolution at approximately 20% the rate of the wild type. Furthermore, the double mutation, D170E + D342E, behaves similarly to the two single mutations. We conclude that D342E alters the equilibrium between different conformers of the cluster in its S1 and S2 states in the same manner as D170E and perturbs the H-bond networks in a similar fashion. This is the second identification of a Mn4CaO5 metal ligand whose mutation influences the equilibrium between the different conformers of the S1 and S2 states without eliminating O2 evolution. This finding has implications for our understanding of the mechanism of O2 formation in terms of catalytically active/inactive conformations of the Mn4CaO5 cluster in its lower oxidation states.
Collapse
Affiliation(s)
- Richard J Debus
- Department of Biochemistry, University of California at Riverside, Riverside, California 92521, United States
| | - Paul H Oyala
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91106, United States
| |
Collapse
|
10
|
Forde A, Maity S, Freixas VM, Fernandez-Alberti S, Neukirch AJ, Kleinekathöfer U, Tretiak S. Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers. J Phys Chem Lett 2024; 15:4142-4150. [PMID: 38593451 DOI: 10.1021/acs.jpclett.4c00362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.
Collapse
Affiliation(s)
- Aaron Forde
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sayan Maity
- School of Science, Constructor University, Campus Ring 1, 28759 Bremen, Germany
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologiia, Univresidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | | | - Amanda J Neukirch
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | | | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
11
|
Noguchi T. Mechanism of Proton Transfer through the D1-E65/D2-E312 Gate during Photosynthetic Water Oxidation. J Phys Chem B 2024; 128:1866-1875. [PMID: 38364371 DOI: 10.1021/acs.jpcb.3c07787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2024]
Abstract
In photosystem II, the D1-E65/D2-E312 dyad in the Cl-1 channel has been proposed to play a pivotal role in proton transfer during water oxidation. However, the precise mechanism remains elusive. Here, the proton transfer mechanism within the Cl-1 channel was investigated using quantum mechanics/molecular mechanics calculations. The molecular vibration of the E65/E312 dyad and its deuteration effect revealed that the recently suggested stepwise proton transfer, i.e., initial proton release from the dyad followed by slow reprotonation, does not occur in the Cl-1 channel. Instead, proton transfer is proposed to take place via a conformational change at the E65/E312 dyad, acting as a gate. In its closed form, a proton is trapped within the dyad, preventing forward proton transfer. This closed form converts into the open form, where protonated D1-E65 provides a hydrogen bond to the water network, thereby facilitating fast Grotthuss-type proton transfer.
Collapse
Affiliation(s)
- Takumi Noguchi
- Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| |
Collapse
|
12
|
Drosou M, Pantazis DA. Comprehensive Evaluation of Models for Ammonia Binding to the Oxygen Evolving Complex of Photosystem II. J Phys Chem B 2024; 128:1333-1349. [PMID: 38299511 PMCID: PMC10875651 DOI: 10.1021/acs.jpcb.3c06304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Revised: 01/08/2024] [Accepted: 01/17/2024] [Indexed: 02/02/2024]
Abstract
The identity and insertion pathway of the substrate oxygen atoms that are coupled to dioxygen by the oxygen-evolving complex (OEC) remains a central question toward understanding Nature's water oxidation mechanism. In several studies, ammonia has been used as a small "water analogue" to elucidate the pathway of substrate access to the OEC and to aid in determining which of the oxygen ligands of the tetramanganese cluster are substrates for O-O bond formation. On the basis of structural and spectroscopic investigations, five first-sphere binding modes of ammonia have been suggested, involving either substitution of an existing H2O/OH-/O2- group or addition as an extra ligand to a metal ion of the Mn4CaO5 cluster. Some of these modes, specifically the ones involving substitution, have already been subject to spectroscopy-oriented quantum chemical investigations, whereas more recent suggestions that postulate the addition of ammonia have not been examined so far with quantum chemistry for their agreement with spectroscopic data. Herein, we use a common structural framework and theoretical methodology to evaluate structural models of the OEC that represent all proposed modes of first-sphere ammonia interaction with the OEC in its S2 state. Criteria include energetic, magnetic, kinetic, and spectroscopic properties compared against available experimental EPR, ENDOR, ESEEM, and EDNMR data. Our results show that models featuring ammonia replacing one of the two terminal water ligands on Mn4 align best with experimental data, while they definitively exclude substitution of a bridging μ-oxo ligand as well as incorporation of ammonia as a sixth ligand on Mn1 or Mn4.
Collapse
Affiliation(s)
- Maria Drosou
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
- Inorganic
Chemistry Laboratory, National and Kapodistrian
University of Athens, Panepistimiopolis, Zografou 15771, Greece
| | - Dimitrios A. Pantazis
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der Ruhr 45470, Germany
| |
Collapse
|
13
|
Yamaguchi K, Miyagawa K, Shoji M, Kawakami T, Isobe H, Yamanaka S, Nakajima T. Theoretical elucidation of the structure, bonding, and reactivity of the CaMn 4O x clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation. PHOTOSYNTHESIS RESEARCH 2023:10.1007/s11120-023-01053-7. [PMID: 37945776 DOI: 10.1007/s11120-023-01053-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Accepted: 09/25/2023] [Indexed: 11/12/2023]
Abstract
This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.
Collapse
Affiliation(s)
- Kizashi Yamaguchi
- Center for Quantum Information and Quantum Biology, Osaka University, Toyonaka, Osaka, 560-0043, Japan.
- RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan.
- SANKEN, Osaka University, Ibaraki, Osaka, 567-0047, Japan.
| | - Koichi Miyagawa
- Center of Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki, 305-8577, Japan
| | - Mitsuo Shoji
- Center of Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki, 305-8577, Japan
| | - Takashi Kawakami
- RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan
- Graduate School of Science, Osaka University, Toyonaka, Osaka, 560-0043, Japan
| | - Hiroshi Isobe
- Research Institute for Interdisciplinary Science, and Graduate School of Natural Science and Technology, Okayama University, Okayama, 700-8530, Japan
| | - Shusuke Yamanaka
- Graduate School of Science, Osaka University, Toyonaka, Osaka, 560-0043, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, Kobe, Hyogo, 650-0047, Japan
| |
Collapse
|
14
|
Sirohiwal A, Pantazis DA. Reaction Center Excitation in Photosystem II: From Multiscale Modeling to Functional Principles. Acc Chem Res 2023; 56:2921-2932. [PMID: 37844298 PMCID: PMC10634305 DOI: 10.1021/acs.accounts.3c00392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Indexed: 10/18/2023]
Abstract
Oxygenic photosynthesis is the fundamental energy-converting process that utilizes sunlight to generate molecular oxygen and the organic compounds that sustain life. Protein-pigment complexes harvest light and transfer excitation energy to specialized pigment assemblies, reaction centers (RC), where electron transfer cascades are initiated. A molecular-level understanding of the primary events is indispensable for elucidating the principles of natural photosynthesis and enabling development of bioinspired technologies. The primary enzyme in oxygenic photosynthesis is Photosystem II (PSII), a membrane-embedded multisubunit complex, that catalyzes the light-driven oxidation of water. The RC of PSII consists of four chlorophyll a and two pheophytin a pigments symmetrically arranged along two core polypeptides; only one branch participates in electron transfer. Despite decades of research, fundamental questions remain, including the origin of this functional asymmetry, the nature of primary charge-transfer states and the identity of the initial electron donor, the origin of the capability of PSII to enact charge separation with far-red photons, i.e., beyond the "red limit" where individual chlorophylls absorb, and the role of protein conformational dynamics in modulating charge-separation pathways.In this Account, we highlight developments in quantum-chemistry based excited-state computations for multipigment assemblies and the refinement of protocols for computing protein-induced electrochromic shifts and charge-transfer excitations calibrated with modern local correlation coupled cluster methods. We emphasize the importance of multiscale atomistic quantum-mechanics/molecular-mechanics and large-scale molecular dynamics simulations, which enabled direct and accurate modeling of primary processes in RC excitation at the quantum mechanical level.Our findings show how differential protein electrostatics enable spectral tuning of RC pigments and generate functional asymmetry in PSII. A chlorophyll pigment on the active branch (ChlD1) has the lowest site energy in PSII and is the primary electron donor. The complete absence of low-lying charge-transfer states within the central pair of chlorophylls excludes a long-held assumption about the initial charge separation. Instead, we identify two primary charge separation pathways, both with the same pheophytin acceptor (PheoD1): a fast pathway with ChlD1 as the primary electron donor (short-range charge-separation) and a slow pathway with PD1PD2 as the initial donor (long-range charge separation). The low-energy spectrum is dominated by two states with significant charge-transfer character, ChlD1δ+PheoD1δ- and PD1δ+PheoD1δ-. The conformational dynamics of PSII allows these charge-transfer states to span wide energy ranges, pushing oxygenic photosynthesis beyond the "red limit". These results provide a quantum mechanical picture of the primary events in the RC of oxygenic photosynthesis, forming a solid basis for interpreting experimental observations and for extending photosynthesis research in new directions.
Collapse
Affiliation(s)
- Abhishek Sirohiwal
- Department
of Biochemistry and Biophysics, Arrhenius Laboratory, Stockholm University, 10691 Stockholm, Sweden
| | - Dimitrios A. Pantazis
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
15
|
Hussein R, Ibrahim M, Bhowmick A, Simon PS, Bogacz I, Doyle MD, Dobbek H, Zouni A, Messinger J, Yachandra VK, Kern JF, Yano J. Evolutionary diversity of proton and water channels on the oxidizing side of photosystem II and their relevance to function. PHOTOSYNTHESIS RESEARCH 2023; 158:91-107. [PMID: 37266800 PMCID: PMC10684718 DOI: 10.1007/s11120-023-01018-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 03/29/2023] [Indexed: 06/03/2023]
Abstract
One of the reasons for the high efficiency and selectivity of biological catalysts arise from their ability to control the pathways of substrates and products using protein channels, and by modulating the transport in the channels using the interaction with the protein residues and the water/hydrogen-bonding network. This process is clearly demonstrated in Photosystem II (PS II), where its light-driven water oxidation reaction catalyzed by the Mn4CaO5 cluster occurs deep inside the protein complex and thus requires the transport of two water molecules to and four protons from the metal center to the bulk water. Based on the recent advances in structural studies of PS II from X-ray crystallography and cryo-electron microscopy, in this review we compare the channels that have been proposed to facilitate this mass transport in cyanobacteria, red and green algae, diatoms, and higher plants. The three major channels (O1, O4, and Cl1 channels) are present in all species investigated; however, some differences exist in the reported structures that arise from the different composition and arrangement of membrane extrinsic subunits between the species. Among the three channels, the Cl1 channel, including the proton gate, is the most conserved among all photosynthetic species. We also found at least one branch for the O1 channel in all organisms, extending all the way from Ca/O1 via the 'water wheel' to the lumen. However, the extending path after the water wheel varies between most species. The O4 channel is, like the Cl1 channel, highly conserved among all species while having different orientations at the end of the path near the bulk. The comparison suggests that the previously proposed functionality of the channels in T. vestitus (Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Hussein et al., Nat Commun 12:6531, 2021) is conserved through the species, i.e. the O1-like channel is used for substrate water intake, and the tighter Cl1 and O4 channels for proton release. The comparison does not eliminate the potential role of O4 channel as a water intake channel. However, the highly ordered hydrogen-bonded water wire connected to the Mn4CaO5 cluster via the O4 may strongly suggest that it functions in proton release, especially during the S0 → S1 transition (Saito et al., Nat Commun 6:8488, 2015; Kern et al., Nature 563:421-425, 2018; Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Sakashita et al., Phys Chem Chem Phys 22:15831-15841, 2020; Hussein et al., Nat Commun 12:6531, 2021).
Collapse
Affiliation(s)
- Rana Hussein
- Department of Biology, Humboldt-Universität Zu Berlin, 10099, Berlin, Germany.
| | - Mohamed Ibrahim
- Department of Biology, Humboldt-Universität Zu Berlin, 10099, Berlin, Germany
| | - Asmit Bhowmick
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Philipp S Simon
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Isabel Bogacz
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Margaret D Doyle
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Holger Dobbek
- Department of Biology, Humboldt-Universität Zu Berlin, 10099, Berlin, Germany
| | - Athina Zouni
- Department of Biology, Humboldt-Universität Zu Berlin, 10099, Berlin, Germany
| | - Johannes Messinger
- Molecular Biomimetics, Department of Chemistry-Ångström, Uppsala University, SE 75120, Uppsala, Sweden
- Department of Chemistry, Umeå University, SE 90187, Umeå, Sweden
| | - Vittal K Yachandra
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Jan F Kern
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| | - Junko Yano
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
| |
Collapse
|
16
|
Khaing EP, Eaton-Rye JJ. Lys264 of the D2 Protein Performs a Dual Role in Photosystem II Modifying Assembly and Electron Transfer through the Quinone-Iron Acceptor Complex. Biochemistry 2023; 62:2738-2750. [PMID: 37606628 DOI: 10.1021/acs.biochem.3c00240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
Bicarbonate (HCO3-) binding regulates electron flow between the primary (QA) and secondary (QB) plastoquinone electron acceptors of Photosystem II (PS II). Lys264 of the D2 subunit of PS II contributes to a hydrogen-bond network that stabilizes HCO3- ligation to the non-heme iron in the QA-Fe-QB complex. Using the cyanobacterium Synechocystis sp. PCC 6803, alanine and glutamate were introduced to create the K264A and K264E mutants. Photoautotrophic growth was slowed in K264E cells but not in the K264A strain. Both mutants accumulated an unassembled CP43 precomplex as well as the CP43-lacking RC47 assembly intermediate, indicating weakened binding of the CP43 precomplex to RC47. Assembly was impeded more in K264E cells than in the K264A strain, but K264A cells were more susceptible to high-light-induced photodamage when assayed using PS II-specific electron acceptors. Furthermore, an impaired repair mechanism was observed in the K264A mutant in protein labeling experiments. Unexpectedly, unlike the K264A strain, the K264E mutant displayed inhibited oxygen evolution following high-light exposure when HCO3- was added to support whole chain electron transport. In both mutants, the decay of chlorophyll fluorescence was slowed, indicating impaired electron transfer between QA and QB. Furthermore, the fluorescence decay kinetics in the K264E strain were insensitive to addition of either formate or HCO3-, whereas HCO3--reversible formate-induced inhibition in the K264A mutant was observed. Exchange of plastoquinol with the membrane plastoquinone pool at the QB-binding site was also retarded in both mutants. Hence, D2-Lys264 possesses key roles in both assembly and activity of PS II.
Collapse
Affiliation(s)
- Ei Phyo Khaing
- Department of Biochemistry, University of Otago, P.O. Box 56, Dunedin 9054, New Zealand
| | - Julian J Eaton-Rye
- Department of Biochemistry, University of Otago, P.O. Box 56, Dunedin 9054, New Zealand
| |
Collapse
|
17
|
Doyle M, Bhowmick A, Wych DC, Lassalle L, Simon PS, Holton J, Sauter NK, Yachandra VK, Kern JF, Yano J, Wall ME. Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography. J Am Chem Soc 2023; 145:14621-14635. [PMID: 37369071 PMCID: PMC10347547 DOI: 10.1021/jacs.3c01412] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Indexed: 06/29/2023]
Abstract
Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.
Collapse
Affiliation(s)
- Margaret
D. Doyle
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Asmit Bhowmick
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - David C. Wych
- Computer,
Computational and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center
for Non-linear Studies, Los Alamos National
Laboratory, Los Alamos, New Mexico 87545, United States
| | - Louise Lassalle
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Philipp S. Simon
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - James Holton
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department
of Biochemistry and Biophysics, University
of California, San Francisco, San
Francisco, California 94158, United States
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Nicholas K. Sauter
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Vittal K. Yachandra
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Jan F. Kern
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Junko Yano
- Molecular
Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Michael E. Wall
- Computer,
Computational and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|