Oxidation and electronic state dependence of proton transfer in the enzymatic cycle of cytochrome P450eryF

Hydrogen bond networks, consisting of hydrogen bonded waters anchored by polar/acidic amino acid sidechains, are often present in the vicinity of the oxygen binding clefts of P450s. Density functional and quantum dynamics calculations of a O(2) binding cleft network model of cytochrome P450eryF(CYP107A1) indicate that such structural motifs facilitate ultrafast proton transfer from network waters to the dioxygen of the reduced oxyferrous species via a multiple proton translocation mechanism with barriers of 7-10 kcal/mol on its doublet ground state, and that the energies of the proton transfer reactant and constrained proton transfer products have an electronic and oxidation state dependence [J. Am. Chem. Soc. 124 (2002) 1430]. In the present study, the origin of the oxidation state dependence is shown to have its roots in differential proton affinities while the electronic state dependence of the reduced oxyferrous heme has its origins in subtle differences in network topologies near the transition state of the initial proton transfer event. Relaxed potential surface scans and unconstrained proton transfer product optimizations indicate that the proton transfer product in both the singlet oxyferrous heme and the reduced oxyferrous heme species in a quartet state are not viable stable (bound) states relative to the reactant form. While the proton affinity of H(3)O(+) is sufficient for it to protonate both the oxyferrous and the reduced oxyferrous heme species, hydrogen bond network stabilized water is only capable of protonating the reduced oxyferrous form. This interpretation is substantiated by study of the NO bound reduced ferrous heme of P450nor, which is isoelectronic with the oxyferrous heme and has a similar proton affinity. Density functional calculations on a more extensive O(2) binding cleft model support the multiple proton translocation mechanism of transfer but indicates that the significant negative charge density on the bound dioxygen of the reduced oxyferrous heme species, in its doublet ground state, polarizes the associated hydrogen bond network sufficiently so as to result in short, strong, low-barrier hydrogen bonds. The computed O-H-O bond distances are less than 2.55 A and have a near degeneracy of the proton transfer reactant and initial (sudden) proton transfer products. These low-barrier hydrogen bond features, in addition to the finding of a (zero point uncorrected) barrier of 1.3 kcal/mol, indicate that proton transfer from water to the distal oxygen should be rapid, facile and may not require large curvature tunneling as originally suggested by use of a smaller model. An initial assessment of protonation of the reduced oxyferrous heme distal oxygen by a model of 6-deoxyerythronolide B (6-DEB) indicates it to be low barrier (3.8 kcal/mol) and exothermic (-2.9 kcal/mol). The combined results indicate the plausibility of simultaneous diprotonation of the distal oxygen of the reduced oxyferrous heme, leading to O-O bond scission, using the combined water network and 6-DEB substrate protonation agents.     

Quantum chemistry analysis of the mechanism of action of proteolytic enzymes. III. Proton transport in serine proteases

The model system for the proton transfer on the amide atom of the substrate leaving group based on the existence of "charge relay system" in the serine type proteases was analysed by the CNDO/2 method. The unfitness of this model to explain the action mechanism of serine proteases was shown. The model system for proton transfer with the water molecule as the intermediate acceptor of the Ser-195 proton was suggested and analysed by the same method. The acylation activation barrier of this system was shown to localize on the stage of synchronous transfer of the Ser-195 alcoholic proton and the water molecule proton hydrogen bound to the His-57 N epsilon 2-atom on the water molecule oxygen atom and the N epsilon 2-atom, respectively. The protonation of substrate in the case of the model system with the water molecule as the intermediate acceptor of proton was demonstrated to begin before the completion of the tetrahedral intermediate substance and the protonated from of the tetrahedral intermediate was shown to form only. A hypothesis considering the role of this water molecule as the nucleophilic reagent on the deacylation stage is presented.     

Theoretical investigations on the hydrogen bonding of nitrile isomers with H2O,HF, NH3 and H2S

The hydrogen bonds formed by the interaction of nitriles with water, hydrogen fluoride, ammonia and hydrogen sulphide have been studied using B3LYP and second-order Møller–Plesset perturbation (MP2) methods and 6-311+ + G(d,p) basis set. The energies and structures of 80 hydrogen-bonded complexes between nitriles and small molecules were examined systematically using B3LYP and MP2 procedure. Categorisation of the hydrogen bonds involved in the various complexes led to an ordering of hydrogen bond donor and acceptor abilities for some functional groups. The interaction energies have been corrected for the basis set superposition error using Boy's counterpoise correction method. The Morokuma energy decomposition analysis reveals that the strong interactions are due to the attractive contributions from the electrostatic (ES), polarisation (PL) and charge transfer (CT) components. The topological parameters, electron density and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital (NBO) analysis has also been performed to study the CT from proton acceptor to the antibonding orbital of the H–Y bond in the proton donor part of complexes. The frequency analysis of C–H…Y bond in the complexes indicates the blue-shifting nature largely in case of sp² hybridised carbon atom.     

Model of serine proteases charge relay system - pcilo study

The PCILO (Perturbative Configuration Interaction Using Localized Orbitals) method has been used to determine the electronic structure of the active center of serine proteases. The results show that the carboxyl group of the aspartic acid residue is the ultimate proton acceptor of the catalytic triad (Asp, His, Ser)^?. In the absence of a substrate the negative charge of the active centre is delocalized, causing polarization of the Ser O^γ-H bond and an increase of the nucleophilicity of the O^γ atom. The proton of the O^γ-H bond of the Ser residue is, however, only partially transferred to the N^?2 atom of imidazole His. The hydration of the model charge relay system is also investigated.     

Structures of prolyl oligopeptidase substrate/inhibitor complexes. Use of inhibitor binding for titration of the catalytic histidine residue

Structure determination of the inactive S554A variant of prolyl oligopeptidase complexed with an octapeptide has shown that substrate binding is restricted to the P4-P2' region. In addition, it has revealed a hydrogen bond network of potential catalytic importance not detected in other serine peptidases. This involves a unique intramolecular hydrogen bond between the P1' amide and P2 carbonyl groups and another between the P2' amide and Nepsilon2 of the catalytic histidine 680 residue. It is argued that both hydrogen bonds promote proton transfer from the imidazolium ion to the leaving group. Another complex formed with the product-like inhibitor benzyloxycarbonyl-glycyl-proline, indicating that the carboxyl group of the inhibitor forms a hydrogen bond with the Nepsilon2 of His(680). Because a protonated histidine makes a stronger interaction with the carboxyl group, it offers a possibility of the determination of the real pK(a) of the catalytic histidine residue. This was found to be 6.25, lower than that of the well studied serine proteases. The new titration method gave a single pK(a) for prolyl oligopeptidase, whose reaction exhibited a complex pH dependence for k(cat)/K(m), and indicated that the observed pK(a) values are apparent. The procedure presented may be applicable for other serine peptidases.     

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.     

Reactions of ketones with aromatics in acid media. The effect of trifluoromethyl groups and the acidity media. A theoretical study

The reactions of acetone, 2,2,2-trifluoroacetone and hexafluoroacetone in methanesulfonic (MSA) and triflic acids (TFSA) with benzene have been studied at M06-2X/6-311+G(d,p) level using cluster-continuum model, where the carbonyl group is explicitly solvated by acid molecules. The introduction of a trifluoromethyl group into the ketone structure reduces the activation energy of the tetrahedral intermediates formation due to an increase of the electrophilicity of the carbonyl group and raises the activation and the reaction energies of the C-O bond cleavage in formed carbinol due to the destabilization of the corresponding carbocation. The introduction of the second trifluoromethyl group inhibits the hydroxyalkylation reaction due to a very strong increase of the reaction and activation energies of the C-O bond cleavage which becomes the rate determining step. The most important catalytic effect of TFSA compared to MSA is not the protonation of the ketone carbonyl, but the reduction of the activation and reaction energies of the carbinol C-O bond cleavage due to better protosolvation properties. Even for TFSA no complete proton transfer to carbonyl oxygen has been observed for free ketones. Therefore, the protonation energies of free ketones cannot be considered as a measure of ketone reactivity in the hydroxyalkylation reaction.

Is there a weak H‐bond → LBHB transition on tetrahedral complex formation in serine proteases?

The transformation of a weak hydrogen bond in the free enzyme into a low-barrier hydrogen bond (LBHB) in the tetrahedral intermediate has been suggested as an important factor facilitating catalysis in serine proteases. In this work, we examine the structure of the H-bond in the Asp102-His57 diad of serine proteases in the free enzyme and in a covalent tetrahedral complex (TC) with a trifluoromethylketone inhibitor. We apply ab initio quantum mechanical calculations to models consisting of a large molecular fragment of the enzyme active site, and the combined effect of the rest of the protein body and the solvation by surrounding bulk water was simulated by a self-consistent reaction field method in our novel QM/SCRF(VS) approach. Potential profiles of adiabatic proton transfer in the Asp102-His57 diad in these model systems were calculated. We conclude that the hydrogen bond in both the free enzyme and in the enzyme-inhibitor TC is a strong ionic asymmetric one-well hydrogen bond, in contrast to a previous suggestion that it is a weak H-bond in the former and a double-well LBHB in the latter.     

Minimal transition state charge stabilization of the oxyanion during peptide bond formation by the ribosome

Peptide bond formation during ribosomal protein synthesis involves an aminolysis reaction between the aminoacyl α-amino group and the carbonyl ester of the growing peptide via a transition state with a developing negative charge, the oxyanion. Structural and molecular dynamic studies have suggested that the ribosome may stabilize the oxyanion in the transition state of peptide bond formation via a highly ordered water molecule. To biochemically investigate this mechanistic hypothesis, we estimated the energetic contribution to catalytic charge stabilization of the oxyanion using a series of transition state mimics that contain different charge distributions and hydrogen bond potential on the functional group mimicking the oxyanion. Inhibitors containing an oxyanion mimic that carried a neutral charge and a mimic that preserved the negative charge but could not form hydrogen bonds had less than a 3-fold effect on inhibitor binding affinity. These observations argue that the ribosome provides minimal transition state charge stabilization to the oxyanion during peptide bond formation via the water molecule. This is in contrast to the substantial level of oxyanion stabilization provided by serine proteases. This suggests that the oxyanion may be neutralized via a proton shuttle, resulting in an uncharged transition state.     

Mechanistic consequences of charge transfer systems in serine proteases and angiotensin: semiempirical computations

Structures and relative energies for the triads of interacting groups in the serine charge relay system of serine proteases and the proposed tyrosine charge relay system of angiotensin II, respectively, were computed according to the standard MNDOC procedure. The most stable configuration obtained for both systems was one in which the histidine residue was negatively charged. These findings indicate that the histidine ring and not the serine hydroxyl group at the active site of serine proteases would be the nucleophilic center which is acylated by substrate. Similarly, the extreme nucleophilicity of the imidazole anion produced by the proposed triad of interacting groups in angiotensin could provoke the formation of a transient covalent bond (acyl intermediate) between receptor and peptide in the receptor activation mechanism.     

On the inactivity of thiol-subtilisin. The role of the intramolecular electric field

Based on computed proton affinities for several model systems, the energetics of proton transfer and the acidity of the catalytic triads Cys-His-Asn (papain). Cys-His-Asp (thiol-subtilisin) and Ser-His-Asp (subtilisin) are discussed. It is shown that in papain the ion-pair Cys--HisH+ exists owing to the intramolecular electric field, and that a similar situation is found in thiol-subtilisin. but not in subtilisin. Assuming similar reaction mechanisms for papain and thiol-subtilisin - i.e. proton transfer from HisH+ to the NH group of the scissile peptide bond - the inactivity of thil-subtilisin towards proteins is explained by the much greater basicity of His in the complex His-Asp- than in His-Asn. In order for this explanation to be consistent, it is tentatively concluded that the catalytic mechanism of the serine proteases is different from that of the cystein proteases, and involves direct transfer of the serine proton to the leaving group in the acylation step.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Reaction mechanism of the binuclear zinc enzyme glyoxalase II - A theoretical study

The glyoxalase system catalyzes the conversion of toxic methylglyoxal to nontoxic d-lactic acid using glutathione (GSH) as a coenzyme. Glyoxalase II (GlxII) is a binuclear Zn enzyme that catalyzes the second step of this conversion, namely the hydrolysis of S-d-lactoylglutathione, which is the product of the Glyoxalase I (GlxI) reaction. In this paper we use density functional theory method to investigate the reaction mechanism of GlxII. A model of the active site is constructed on the basis of the X-ray crystal structure of the native enzyme. Stationary points along the reaction pathway are optimized and the potential energy surface for the reaction is calculated. The calculations give strong support to the previously proposed mechanism. It is found that the bridging hydroxide is capable of performing nucleophilic attack at the substrate carbonyl to form a tetrahedral intermediate. This step is followed by a proton transfer from the bridging oxygen to Asp58 and finally C-S bond cleavage. The roles of the two zinc ions in the reaction mechanism are analyzed. Zn2 is found to stabilize the charge of tetrahedral intermediate thereby lowering the barrier for the nucleophilic attack, while Zn1 stabilizes the charge of the thiolate product, thereby facilitating the C-S bond cleavage. Finally, the energies involved in the product release and active-site regeneration are estimated and a new possible mechanism is suggested.     

Water as a structural element in a channel: gating in the Kcsa channel, and implications for voltage-gated ion channels

Water is becoming understood as a structural element in proteins. Here we are concerned with one particular type of protein, ion channels. The S. Lividans KcsA K(+) channel, the X-ray structure of which is known, is gated by protons (i.e, by a drop in pH). Ab initio calculations suggest that an H(5)O(2) group, partially charged, connects the E118 residues in the gating region, when the four residues have a -2 net charge, but that the hydrogen bonding is not strong enough to do this when the charge becomes -1. The H(5)O(2) group would block the channel, in the -2 state, and prevent motion of the four transmembrane (TM) segments of the protein, by binding them. With the weaker bond in the -1 state, the TM segments would be able to separate (as they have been found to do experimentally, opening the channel. Voltage gated channels have four additional TM segments for each of the four domains of the channel protein. These appear to allow motion of protons; in fact there is evidence that the initial step in gating must be the transfer of a proton. We have earlier shown that the transfer of a single proton between two methylamines under the influence of a field is possible, as proton tunneling. Subsequent steps are hypothesized to result from four proton transfer cascades of about three protons each, triggered by the initial proton transfer. We suggest that the extra 4 TM segments of the voltage gated channel act as a voltage to proton-current transducer. Water, held by hydrogen bonds, is also suggested as the source of the accessibility data found with MTS reagents, based largely on simulations, our earlier Monte Carlo simulations as well as molecular dynamics studies reported by others. These waters may also play a structural role in the protein.     

Hydrogen bonds in crystal structures of amino acids, peptides and related molecules

The results of a survey of 439 hydrogen bonds in 95 recently determined crystal structures of amino acids, peptides and related molecules suggest that the following generalizations hold true for linear (angle X-H---Y greater than 150 degrees) hydrogen bonds. (1) The charge on the acceptor group does not influence the length of a hydrogen bond. (2) For a given acceptor group, the hydrogen bond lengths increase in the order imidazolium N--H less than ammonium N-H less than guanidinium N-H; this order holds true for oxygen anion acceptor groups. Cl-ions and the uncharged oxygen of water molecules. (3) The uncharged imidazole N-H group forms shorter hydrogen than the amide N-H GROUP. (4) The carboxyl O-H groups form shorter hydrogen bonds than other hydroxyl groups. (5) The hydrogen bonds involving a halogen ion are longer than hydrogen bonds with other acceptors when corrected for their longer van der Walls radii. The observed differences between the lengths of hydrogen bonds formed by different donor and acceptor groups in amino acids and peptides, imply differences in the energetics of their formation.     

Influence of key residues on the reaction mechanism of the cAMP-dependent protein kinase

The reaction mechanism of the catalytic phosphoryl transfer of cAMP-dependent protein kinase (cAPK) was investigated by semi-empirical AM1 molecular orbital computations of an active site model system derived from the crystal structure of the catalytic subunit of the enzyme. The activation barrier is calculated as 20.7 kcal mol(-1) and the reaction itself to be exothermic by 12.2 kcal mol(-1). The active site residue Asp166, which was often proposed to act as a catalytic base, does not accept a proton in any of the reaction steps. Instead, the hydroxyl hydrogen of serine is shifted to the simultaneously transferred phosphate group of ATP. Although the calculated transition state geometry indicates an associative phosphoryl transfer, no concentration of negative charge is found. To study the influence of protein mutations on the reaction mechanism, we compared two-dimensional energy hypersurfaces of the protein kinase wild-type model and a corresponding mutant in which Asp166 was replaced by alanine. Surprisingly, they show similar energy profiles despite the experimentally known decrease of catalytic activity for corresponding mutants. Furthermore, a model structure was examined, where the charged NH3 group of Lys168 was replaced by a neutral methyl group. The energetic hypersurface of this hypothetical mutant shows two possible pathways for phosphoryl transfer, which both require significantly higher activation energies than the other systems investigated, while the energetic stabilization of the reaction product is similar in all systems. As the position of the amino acid side chains and the substrate peptide is virtually unchanged in all model systems, our results suggest that the exchange of Asp166 by other amino acid is less important to the phosphoryl transfer itself, but crucial to maintain the configuration of the active site in vivo. The positively charged side chain of Lys168, however, is necessary to stabilize the intermediate reaction states, particularly the side chain of the substrate peptide.     

Design of molecular switching and signaling based on proton transfer in 2-hydroxy Schiff bases: a computational study

The present work aims to exploit the possibility of using the tautomerism in 2-hydroxy Schiff bases for molecular switching. The enol imine (E)? enaminone (K) tautomerization in a series of 2-hydroxy Schiff bases have been investigated theoretically at the DFT/B3LYP/6-311G** level of theory. The intramolecular proton transfer processes have been explored, transition structures have been located and characterized. The kinetics and thermodynamics of the proton transfer process, and its time scale have been computed and discussed in the framework of the suitability as molecular switches. Substituent effects have been computed and its effect on the enthalpy changes (?H*) and activation energies (?G*) have been analyzed and discussed. Nonspecific solvent effects have also been taken into account by using the polarized continuum model (IPCM) of two different solvent. The tautomerization energies are decreased and hence the endothermic nature of the enol imine ? enaminone tautomerization. The potential energy barriers, on the other hand, are increased due to the relative destabilization of the transition states. The NBO charge populations show that there is a high positive charge on the hydrogen atom during the process in all cases, which confirms that the proton transfer proceeds through a three-center interaction. The proton transfer processes, in all cases studied are kinetically allowed. The low potential energy barrier suggests that interconversion between the two tautomeric forms is spontaneous and the two forms may coexist.     

Mimicking the photosynthetic system with strong hydrogen bonds to promote proton electron concerted reactions

A Copper(2+) complex with a Cu^II–C bond containing sp³ configuration was used to investigate the role of strong hydrogen bonds in proton coupled electron transfer (PCET) reactions. The only example of a Cu^II–C system realized so far is that using tris-(pyridylthio)methyl (tptm) as a tetradentate tripodal ligand. Using this ligand, [CuF(tptm)] and [Cu(tptm)(OH)] have been prepared. The former complex forms supra-molecular arrays of layers of the complex between which hydroquinone is intercalated in the crystalline phase. This hydroquinone intercalation crystal was prepared via the photochemical conversion of quinone during the crystallization process. This conversion reaction probably involves a proton coupled electron transfer process. The nuclear magnetic resonance spectroscopic analysis of the reaction mixture shows the presence of Cu(III) during the conversion reaction. These results strongly suggest the presence of the molecular aggregate of the [CuF(tptm)] complex, water and quinone in the solution phase during the quinone to hydroquinone conversion. The presence of this type of aggregate requires a strong hydrogen bond between the [CuF(tptm)] complex and water. The presence of this particular hydrogen bond is a unique character of such a complex that has the Cu^II–C bond. This complex is used as a model for photosynthetic water splitting since the photoconversion of quinone to hydroquinone also involves the production of oxygen from water.     

On the stereochemistry of catalysis by serine proteases

From stereochemical considerations and model building the following conclusions were drawn for the stereochemistry of the catalytic steps of chymotrypsin and subtilisin. (1) In contrast to previous stereochemical investigations, rotation of 120° or more of the oxygen atom of the “reactive” serine residue is not possible in the course of the reaction with specific substrates. (2) During catalysis the serine oxygen atom is approximately in the position found in the crystalline enzyme, i.e. at a distance of about 3 Å from the nitrogen atom of the catalytically important histidine residue. (3) The detailed stereochemical mechanism involves the formation of a strained tetrahedral intermediate and a strained acylenzyme. The strain energy is supplied by the formation of a hydrogen bond between the enzyme and a specific substrate. (4) The geometry of proton transfers in the intimate encounter complex of chymotrypsin is slightly but significantly different from that of subtilisin.     

Low-barrier hydrogen bond plays key role in active photosystem II--a new model for photosynthetic water oxidation

The function and mechanism of Tyr(Z) in active photosystem II (PSII) is one of the long-standing issues in the study of photosynthetic water oxidation. Based on recent investigations on active PSII and theoretical studies, a new model is proposed, in which D1-His190 acts as a bridge, to form a low-barrier hydrogen bond (LBHB) with Tyr(Z), and a coordination bond to Mn or Ca ion of the Mn-cluster. Accordingly, this new model differs from previous proposals concerning the mechanism of Tyr(Z) function in two aspects. First, the LBHB plays a key role to decrease the activation energy for Tyr(Z) oxidation and Tyr(Z)(.) reduction during photosynthetic water oxidation. Upon the oxidation of Tyr(Z), the hydrogen bond between Tyr(Z) and His190 changes from a LBHB to a weak hydrogen bond, and vice versa upon Tyr(Z)(.) reduction. In both stages, the electron transfer and proton transfer are coupled. Second, the positive charge formed after Tyr(Z) oxidation may play an important role for water oxidation. It can be delocalized on the Mn-cluster, thus helps to accelerate the proton release from substrate water on Mn-cluster. This model is well reconciled with observations of the S-state dependence of Tyr(Z) oxidation and Tyr(Z)(.) reduction, proton release, isotopic effect and recent EPR experiments. Moreover, the difference between Tyr(Z) and Tyr(D) in active PSII can also be readily rationalized. The His190 binding to the Mn-cluster predicted in this model is contradictious to the recent structure data, however, it has been aware that the crystal structure of the Mn-cluster and its environment are significantly modified by X-ray due to radiation damage and are different from that in active PSII. It is suggested that the His190 may be protonated during the radiation damage, which leads to the loss of its binding to Mn-cluster and the strong hydrogen bond with Tyr(Z). This type of change arising from radiation damage has been confirmed in other enzyme systems.