Theoretical modelling of biomolecular systems I. Large-scale QM/MM calculations of hydrogen-bonding networks of the oxygen evolving complex of photosystem II

Mapping the Oxygens in the Oxygen-Evolving Complex of Photosystem II by Their Nucleophilicity Using Quantum Descriptors

The oxygen-evolving complex (OEC) of Photosystem II catalyzes the water-splitting reaction using solar energy. Thus, understanding the reaction mechanism will inspire the design of biomimetic artificial catalysts that convert solar energy to chemical energy. Conceptual Density Functional Theory (CDFT) focuses on understanding the reactivity of molecules and the atomic contribution to the overall nucleophilicity and electrophilicity of the molecule using quantum descriptors. However, this method has not been applied to the OEC before. Here, we use Fukui functions and the dual descriptor to provide quantitative measures of the nucleophilicity and electrophilicity of oxygens in the OEC for different models in different S states. Our results show that the μ-oxo bridges connected to terminal Mn4 are nucleophilic, and those in the cube formed by Mn1, Mn2, and Mn3 are mostly electrophilic. The dual descriptors of the bridging oxygens in the OEC showed a similar reactivity to that of bridging oxygens in Mn model compounds. However, the terminal water W1, which is bound to Mn4, showed very strong reactivity in some of the S3 models. Thus, our calculations support the model that proposes the formation of the O2 molecule through nucleophilic attack by a terminal water.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox 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

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.

Electronic and spin structures of CaMn4Ox clusters in the S0 state of the oxygen evolving complex of photosystem II. Domain-based local pair natural orbital (DLPNO) coupled-cluster (CC) calculations using optimized geometries and natural orbitals (UNO) by hybrid density functional theory (HDFT) calculations

Domain-based local pair natural orbital (DLPNO) coupled cluster single and double (CCSD) with triple perturbation (T) correction methods were performed to elucidate the relative stabilities of ten different intermediate structures of the CaMn4Ox cluster in the S0 state of the oxygen evolving complex (OEC) of photosystem II (PSII). Full geometry optimizations of all the S0 intermediates were performed by the UB3LYP-D3/Def2-TZVP methods, providing the assumed geometrical structures and starting natural orbitals (UNO) for DLPNO-CCSD(T)/Def2TZVP calculations. The effective exchange integrals (J) for the spin Hamiltonian models for the ten intermediates were obtained by the UB3LYP/Def2-TZVP calculations followed by the general spin projections. DLPNO-CCSD(T) calculations followed by the CBS extrapolation procedure elucidated that the (II, III, IV, IV) and (III, III, III, IV) valence states in the CaMn4O5 cluster of the OEC of the PS II were nearly degenerated in energy in the S0 state, indicating an important role of dynamical electron correlation effects for the valence and spin fluctuations in strongly correlated electron systems (SCESs) consisting of 3d transition metals.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

The S3 State of the Oxygen-Evolving Complex: Overview of Spectroscopy and XFEL Crystallography with a Critical Evaluation of Early-Onset Models for O–O Bond Formation

The catalytic cycle of the oxygen-evolving complex (OEC) of photosystem II (PSII) comprises five intermediate states Si (i = 0–4), from the most reduced S0 state to the most oxidized S4, which spontaneously evolves dioxygen. The precise geometric and electronic structure of the Si states, and hence the mechanism of O–O bond formation in the OEC, remain under investigation, particularly for the final steps of the catalytic cycle. Recent advances in protein crystallography based on X-ray free-electron lasers (XFELs) have produced new structural models for the S3 state, which indicate that two of the oxygen atoms of the inorganic Mn4CaO6 core of the OEC are in very close proximity. This has been interpreted as possible evidence for “early-onset” O–O bond formation in the S3 state, as opposed to the more widely accepted view that the O–O bond is formed in the final state of the cycle, S4. Peroxo or superoxo formation in S3 has received partial support from computational studies. Here, a brief overview is provided of spectroscopic information, recent crystallographic results, and computational models for the S3 state. Emphasis is placed on computational S3 models that involve O–O formation, which are discussed with respect to their agreement with structural information, experimental evidence from various spectroscopic studies, and substrate exchange kinetics. Despite seemingly better agreement with some of the available crystallographic interpretations for the S3 state, models that implicate early-onset O–O bond formation are hard to reconcile with the complete line of experimental evidence, especially with X-ray absorption, X-ray emission, and magnetic resonance spectroscopic observations. Specifically with respect to quantum chemical studies, the inconclusive energetics for the possible isoforms of S3 is an acute problem that is probably beyond the capabilities of standard density functional theory.

Geometric and electronic structures of the synthetic Mn₄CaO₄ model compound mimicking the photosynthetic oxygen-evolving complex

Water oxidation by photosystem II (PSII) converts light energy into chemical energy with the concomitant production of molecular oxygen, both of which are indispensable for sustaining life on Earth. This reaction is catalyzed by an oxygen-evolving complex (OEC) embedded in the huge PSII complex, and its mechanism remains elusive in spite of the extensive studies of the geometric and electronic structures. In order to elucidate the water-splitting mechanism, synthetic approaches have been extensively employed to mimic the native OEC. Very recently, a synthetic complex [Mn4CaO4(Bu(t)COO)8(py)(Bu(t)COOH)2] (1) closely mimicking the structure of the native OEC was obtained. In this study, we extensively examined the geometric, electronic and spin structures of 1 using the density functional theory method. Our results showed that the geometric structure of 1 can be accurately reproduced by theoretical calculations, and revealed many similarities in the ground valence and spin states between 1 and the native OEC. We also revealed two different valence states in the one-electron oxidized state of 1 (corresponding to the S2 state), which lie in the lower and higher ground spin states (S = 1/2 and S = 5/2), respectively. One remarkable difference between 1 and the native OEC is the presence of a non-negligible antiferromagnetic interaction between the Mn1 and Mn4 sites, which slightly influenced their ground spin structures (spin alignments). The major reason causing the difference can be attributed to the short Mn1-O5 and Mn1-Mn4 distances in 1. The introduction of the missing O4 atom and the reorientation of the Ca coordinating ligands improved the Mn1-O5 and Mn1-Mn4 distances comparable to the native OEC. These modifications will therefore be important for the synthesis of further advanced model complexes more closely mimicking the native OEC beyond 1.

The MOD-QM/MM Method: Applications to Studies of Photosystem II and DNA G-Quadruplexes

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and catalytic sites embedded in macrobiomolecules (eg, proteins and nucleic acids). QM/MM methodologies are highly efficient since they implement quantum chemistry methods for modeling only the portion of the system involving bond-breaking/forming processes (QM layer), as influenced by the surrounding molecular environment described in terms of molecular mechanics force fields (MM layer). Some of the limitations of QM/MM methods when polarization effects are not explicitly considered include the approximate treatment of electrostatic interactions between QM and MM layers. Here, we review recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in biomacromolecules and structural refinement, beyond the common limitations of QM/MM hybrid methods. We focus on photosystem II (PSII) with emphasis on the description of the oxygen-evolving complex (OEC) and its high-resolution extended X-ray absorption fine structure spectra (EXAFS) in conjunction with Monte Carlo structural refinement. Furthermore, we review QM/MM structural refinement studies of DNA G4 quadruplexes with embedded monovalent cations and direct comparisons to NMR data.

How Accurately Can Extended X-ray Absorption Spectra Be Predicted from First Principles? Implications for Modeling the Oxygen-Evolving Complex in Photosystem II

First principle calculations of extended X-ray absorption fine structure (EXAFS) data have seen widespread use in bioinorganic chemistry, perhaps most notably for modeling the Mn4Ca site in the oxygen evolving complex (OEC) of photosystem II (PSII). The logic implied by the calculations rests on the assumption that it is possible to a priori predict an accurate EXAFS spectrum provided that the underlying geometric structure is correct. The present study investigates the extent to which this is possible using state of the art EXAFS theory. The FEFF program is used to evaluate the ability of a multiple scattering-based approach to directly calculate the EXAFS spectrum of crystallographically defined model complexes. The results of these parameter free predictions are compared with the more traditional approach of fitting FEFF calculated spectra to experimental data. A series of seven crystallographically characterized Mn monomers and dimers is used as a test set. The largest deviations between the FEFF calculated EXAFS spectra and the experimental EXAFS spectra arise from the amplitudes. The amplitude errors result from a combination of errors in calculated S0(2) and Debye-Waller values as well as uncertainties in background subtraction. Additional errors may be attributed to structural parameters, particularly in cases where reliable high-resolution crystal structures are not available. Based on these investigations, the strengths and weaknesses of using first-principle EXAFS calculations as a predictive tool are discussed. We demonstrate that a range of DFT optimized structures of the OEC may all be considered consistent with experimental EXAFS data and that caution must be exercised when using EXAFS data to obtain topological arrangements of complex clusters.

Native structure of photosystem II at 1.95 Å resolution viewed by femtosecond X-ray pulses

Photosynthesis converts light energy into biologically useful chemical energy vital to life on Earth. The initial reaction of photosynthesis takes place in photosystem II (PSII), a 700-kilodalton homodimeric membrane protein complex that catalyses photo-oxidation of water into dioxygen through an S-state cycle of the oxygen evolving complex (OEC). The structure of PSII has been solved by X-ray diffraction (XRD) at 1.9 ångström resolution, which revealed that the OEC is a Mn4CaO5-cluster coordinated by a well defined protein environment. However, extended X-ray absorption fine structure (EXAFS) studies showed that the manganese cations in the OEC are easily reduced by X-ray irradiation, and slight differences were found in the Mn-Mn distances determined by XRD, EXAFS and theoretical studies. Here we report a 'radiation-damage-free' structure of PSII from Thermosynechococcus vulcanus in the S1 state at a resolution of 1.95 ångströms using femtosecond X-ray pulses of the SPring-8 ångström compact free-electron laser (SACLA) and hundreds of large, highly isomorphous PSII crystals. Compared with the structure from XRD, the OEC in the X-ray free electron laser structure has Mn-Mn distances that are shorter by 0.1-0.2 ångströms. The valences of each manganese atom were tentatively assigned as Mn1D(III), Mn2C(IV), Mn3B(IV) and Mn4A(III), based on the average Mn-ligand distances and analysis of the Jahn-Teller axis on Mn(III). One of the oxo-bridged oxygens, O5, has significantly longer distances to Mn than do the other oxo-oxygen atoms, suggesting that O5 is a hydroxide ion instead of a normal oxygen dianion and therefore may serve as one of the substrate oxygen atoms. These findings provide a structural basis for the mechanism of oxygen evolution, and we expect that this structure will provide a blueprint for the design of artificial catalysts for water oxidation.