New model for a theoretical density functional theory investigation of the mechanism of the carbonic anhydrase: how does the internal bicarbonate rearrangement occur?

A theoretical density functional theory (DFT, B3LYP) investigation has been carried out on the catalytic cycle of the carbonic anhydrase. A model system including the Glu106 and Thr199 residues and the "deep" water molecule has been used. It has been found that the nucleophilic attack of the zinc-bound OH on the CO(2) molecule has a negligible barrier (only 1.2 kcal mol(-1)). This small value is due to a hydrogen-bond network involving Glu106, Thr199, and the deep water molecule. The two usually proposed mechanisms for the internal bicarbonate rearrangement have been carefully examined. In the presence of the two Glu106 and Thr199 residues, the direct proton transfer (Lipscomb mechanism) is a two-step process, which proceeds via a proton relay network characterized by two activation barriers of 4.4 and 9.0 kcal mol(-1). This pathway can effectively compete with a rotational mechanism (Lindskog mechanism), which has a barrier of 13.2 kcal mol(-1). The fast proton transfer found here is basically due to the effect of the Glu106 residue, which stabilizes an intermediate situation where the Glu106 fragment is protonated. In the absence of Glu106, the barrier for the proton transfer is much larger (32.3 kcal mol(-1)) and the Lindskog mechanism becomes favored.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.     

Mechanism of the tyrosine ammonia lyase reaction-tandem nucleophilic and electrophilic enhancement by a proton transfer

Quantum mechanics/molecular mechanics calculations in tyrosine ammonia lyase (TAL) ruled out the hypothetical Friedel-Crafts (FC) route for ammonia elimination from L-tyrosine due to the high energy of FC intermediates. The calculated pathway from the zwitterionic L-tyrosine-binding state (0.0?kcal mol(-1)) to the product-binding state ((E)-coumarate+H(2)N-MIO; -24.0?kcal mol(-1); MIO = 3,5-dihydro-5-methylidene-4H-imidazol-4-one) involves an intermediate (IS, -19.9?kcal mol(-1)), which has a covalent bond between the N atom of the substrate and MIO, as well as two transition states (TS1 and TS2). TS1 (14.4?kcal mol(-1)) corresponds to a proton transfer from the substrate to the N1 atom of MIO by Tyr300-OH. Thus, a tandem nucleophilic activation of the substrate and electrophilic activation of MIO happens. TS2 (5.2?kcal mol(-1)) indicates a concerted C-N bond breaking of the N-MIO intermediate and deprotonation of the pro-S β position by Tyr60. Calculations elucidate the role of enzymic bases (Tyr60 and Tyr300) and other catalytically relevant residues (Asn203, Arg303, and Asn333, Asn435), which are fully conserved in the amino acid sequences and in 3D structures of all known MIO-containing ammonia lyases and 2,3-aminomutases.     

New computational evidence for the catalytic mechanism of carbonic anhydrase

Some aspects of the catalytic mechanism of HCA have been investigated. Either a zinc-bound water or a zinc-bound hydroxide has been considered as a nucleophile attacking CO ₂. No reaction path exists in the former case, while a transition state for the nucleophilic attack has been located in the latter (barrier of 7.6 kcal mol⁻¹). This activation energy is determined by the breaking of the hydrogen-bond network that shields the zinc-bound hydroxide when the CO ₂ molecule approaches the reaction center. No ambiguity exists about the mechanism for the internal rearrangement of the zinc–bicarbonate complex. The rotation pathway (Lindskog mechanism) proposed by many authors is too energy demanding since it causes the breaking of the hydrogen-bond network around the bicarbonate. The only possible rearrangement mechanism is a proton transfer (Lipscomb) that occurs in two steps (each step corresponding to a double proton transfer) and involves the Thr199 residue as a proton shuttle. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Contribution to the Fernando Bernardi Memorial Issue.     

DFT study of the glutathione peroxidase-like activity of phenylselenol incorporating solvent-assisted proton exchange

Modeling of the glutathione peroxidase-like activity of phenylselenol has been accomplished using density-functional theory and solvent-assisted proton exchange (SAPE). SAPE is a modeling technique intended to mimic solvent participation in proton transfer associated with chemical reaction. Within this method, explicit water molecules incorporated into the gas-phase model allow relay of a proton through the water molecules from the site of protonation in the reactant to that in the product. The activation barriers obtained by SAPE for the three steps of the GPx-like mechanism of PhSeH fall within the limits expected for a catalytic system at physiological temperatures (DeltaG(1)++ = 19.1 kcal/mol; DeltaG(2)++= 6.6 kcal/mol; G(3)++ = 21.7 kcal/mol) and are significantly lower than studies which require direct proton transfer. The size of the SAPE network is also considered for the model of the reduction of the selenenic acid, step 2 of the GPx-like cycle. Use of a four-water network better accommodates the reaction pathway and reduces the activation barrier by 5 kcal/mol over the two-water model.     

Theoretical studies on pyridoxal 5'-phosphate-dependent transamination of alpha-amino acids

Density functional methods have been applied to investigate the irreversible transamination between glyoxylic acid and pyridoxamine analog and the catalytic mechanism for the critical [1,3] proton transfer step in aspartate aminotransferase (AATase). The results indicate that the catalytic effect of pyridoxal 5'-phosphate (PLP) may be attributed to its ability to stabilize related transition states through structural resonance. Additionally, the PLP hydroxyl group and the carboxylic group of the amino acid can shuttle proton, thereby lowering the barrier. The rate-limiting step is the tautomeric conversion of the aldimine to ketimine by [1,3] proton transfer, with a barrier of 36.3 kcal/mol in water solvent. A quantum chemical model consisting 142 atoms was constructed based on the crystal structure of the native AATase complex with the product L-glutamate. The electron-withdrawing stabilization by various residues, involving Arg386, Tyr225, Asp222, Asn194, and peptide backbone, enhances the carbon acidity of 4'-C of PLP and Calpha of amino acid. The calculations support the proposed proton transfer mechanism in which Lys258 acts as a base to shuttle a proton from the 4'-C of PLP to Calpha of amino acid. The first step (proton transfer from 4'-C to lysine) is shown to be the rate-limiting step. Furthermore, we provided an explanation for the reversibility and specificity of the transamination in AATase.     

QM/MM study of catalytic methyl transfer by the N5-glutamine SAM-dependent methyltransferase and its inhibition by the nitrogen analogue of coenzyme

The combined density functional quantum mechanical/molecular mechanical (QM/MM) approach has been used to investigate methyl-transfer reactions catalyzed by the N(5)-glutamine S-adenosyl-L-methionine (SAM)-dependent methyltransferase (HemK) and the coenzyme-modified HemK with the replacement of SAM by a nitrogen analogue. Calculations reveal that the catalytic methyl transfer by HemK is an energy-favored process with an activation barrier of 15.7 kcal/mol and an exothermicity of 12.0 kcal/mol, while the coenzyme-modified HemK is unable to catalyze the methyl transfer because of a substantial barrier of 20.6 kcal/mol and instability of the product intermediate. The results lend support to the experimental proposal that the nitrogen analogue of the SAM coenzyme should be a practicable inhibitor for the catalytic methyl transfer by HemK. Comparative QM/MM calculations show that the protein environment, especially the residues Asn197 and Pro198 in the active site, plays a pivotal role in stabilizing the transition state and regulating the positioning of reactive groups.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Fundamental reaction pathway and free energy profile for inhibition of proteasome by Epoxomicin

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to provide the first detailed computational study on the possible mechanisms for reaction of proteasome with a representative peptide inhibitor, Epoxomicin (EPX). The calculated results reveal that the most favorable reaction pathway consists of five steps. The first is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is nucleophilic attack on the carbonyl carbon of EPX by the negatively charged Thr1-O(γ) atom, followed by a proton transfer from Thr1-N(z) to the carbonyl oxygen of EPX (third step). Then, Thr1-N(z) attacks on the carbon of the epoxide group of EPX, accompanied by the epoxide ring-opening (S(N)2 nucleophilic substitution) such that a zwitterionic morpholino ring is formed between residue Thr1 and EPX. Finally, the product of morpholino ring is generated via another proton transfer. Noteworthy, Thr1-O(γ) can be activated directly by Thr1-N(z) to form the zwitterionic intermediate (with a free energy barrier of only 9.9 kcal/mol), and water cannot assist the rate-determining step, which is remarkably different from the previous perception that a water molecule should mediate the activation process. The fourth reaction step has the highest free energy barrier (23.6 kcal/mol) which is reasonably close to the activation free energy (～21-22 kcal/mol) derived from experimental kinetic data. The obtained novel mechanistic insights should be valuable for not only future rational design of more efficient proteasome inhibitors but also understanding the general reaction mechanism of proteasome with a peptide or protein.     

Density-functional study of mechanisms for the cofactor-free decarboxylation performed by uroporphyrinogen III decarboxylase

Uroporphyrinogen III decarboxylase catalyzes the fifth step in heme biosynthesis: the elimination of carboxyl groups from the four acetate side chains of uroporphyrinogen III to yield coproporphyrinogen III. The enzyme acts by successively protonating each of the four pyrrole rings present in the substrate, thereby allowing decarboxylation of their side chains, but the identity of the proton donors has not been established yet. Tyr164 has been suggested as a proton donor, and Asp86 has been proposed to act either as a proton donor or as an intermediate-stabilizing residue. We have performed density-functional calculations to study this reaction mechanism, and found that the rate-limiting step is substrate protonation, rather than decarboxylation. Surprisingly, whereas Tyr164 is unable to protonate the substrate, this protonation can be effected by a nearby arginine residue (Arg37), with a free energy barrier of 21.4 kcal.mol(-1), in remarkable agreement with the experimental value of 19.5 kcal.mol(-1). The central positioning of this residue in close proximity to all four pyrrole rings in the substrate may play a key role in the sequential activation of each of these moieties.     

The proton uptake channel of bacteriorhodopsin as studied by a photoelectrochemical method

A series of the mutant proteins (D96N, D96N/D85N, D115N, L93T, T46V, V49A) where the residues are located at the cytoplasmic domain of bacteriorhodopsin (bR) were studied photoelectrochemically and their photocurrent response characteristics at the electrode/electrolyte interface were compared with those of the wild-type bR. While the wild-type bR of normal proton pumping activity yields symmetrical cathodic (positive) and anodic (negative) responses, corresponding to proton release and proton uptake, respectively, these mutants, with the exception of D115N, showed diminished amplitudes in the negative response. This indicates retardation of proton translocation from the cytoplasmic surface to the retinal Schiff base. The mutation that gave the strongest influence on the negative response was D96N while moderate influence was obtained with L93T, T46V, and V49A. These results suggest that residues other than D96 also participate in the cytoplasmic proton uptake channel, either by interacting with D96 directly or by forming a hydrogen-bonded network with water molecules. The D96N/D85N double mutant yielded little response at neutral pH, but the response was partially recovered by addition of azide, while it was fully recovered in the single mutant D96N. The D115N mutant showed the response profile that closely resembles the wild-type, indicating that D115 is not crucially involved in the event of proton transfer relay at the cytoplasmic region. It was also found that every mutant in this study releases protons prior to uptake at the other membrane surface, as does the wild-type.     

Computational study of the phosphoryl transfer catalyzed by a cyclin-dependent kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.     

Theoretical study of the mechanism of peptide ring formation in green fluorescent protein

Density functional calculations using hybrid functionals (B3LYP) have been performed to study the mechanism of peptide ring formation in green fluorescent protein (GFP). Several different chemical models were used ranging from a minimal model of the ring formation to a model including the full side chains of the groups involved in forming the peptide ring. The surrounding protein was described using a dielectric cavity model. The previously most accepted mechanism was found to lead to an endothermic cyclization of about 10 kcal/mol, independent of chemical model used. The formation of the required dihydro‐imidazolone intermediate was found to be even more endothermic with 16–18 kcal/mol. In contrast, another mechanism where the dehydration of residue 66 precedes cyclization was found to be exothermic by 1.9 kcal/mol and to go over an endothermic intermediate of only 6.7 kcal/mol. Correcting these results using the more accurate G2‐M scheme leads to an intermediate with an energy of only +3.7 kcal/mol and an overall exothermicity of 4.7 kcal/mol. Possible transition states involving proton transfer steps were also investigated. Comparisons are made to the similar and more well‐known deamination reaction of Asn‐Gly sequences in peptides, for which good agreement is obtained with experiments. The results are discussed with respect to available experiments. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 169–186, 2001     

A Marcus treatment of rate constants for protonation of ring-substituted alpha-methoxystyrenes: intrinsic reaction barriers and the shape of the reaction coordinate

Rate and equilibrium constants were determined for protonation of ring-substituted -methoxystyrenes by hydronium ion and by carboxylic acids to form the corresponding ring-substituted alpha-methyl alpha-methoxybenzyl carbocations at 25 degrees C and I = 1.0 (KCl). The thermodynamic barrier to carbocation formation increases by 14.5 kcal/mol as the phenyl ring substituent(s) is changed from 4-MeO- to 3,5-di-NO2-, and as the carboxylic acid is changed from dichloroacetic to acetic acid. The Br?nsted coefficient alpha for protonation by carboxylic acids increases from 0.67 to 0.77 over this range of phenyl ring substituents, and the Br?nsted coefficient beta for proton transfer increases from 0.63 to 0.69 as the carboxylic acid is changed from dichloroacetic to acetic acid. The change in these Br?nsted coefficients with changing reaction driving force, (inverted theta)alpha/ (inverted theta) deltaG(av) degrees=(inverted theta)beta/(inverted theta)delta G(av) degrees= 1/8lambda = 0.011, is used to calculate a Marcus intrinsic reaction barrier of lambda= 11 kcal/mol which is close to the barrier of 13 kcal/mol for thermoneutral proton transfer between this series of acids and bases. The value of alpha= 0.66 for thermoneutral proton transfer is greater than alpha= 0.50 required by a reaction that follows the Marcus equation. This elevated value of beta may be due to an asymmetry in the reaction coordinate that arises from the difference in the intrinsic barriers for proton transfer at the oxygen acid reactant and resonance-stabilized carbon acid product.     

Computational study on nonenzymatic peptide bond cleavage at asparagine and aspartic acid

Nonenzymatic peptide bond cleavage at asparagine (Asn) and glutamine (Gln) residues has been observed during peptide deamidation experiments; cleavage has also been reported at aspartic acid (Asp) and glutamic acid (Glu) residues. Although peptide backbone cleavage at Asn is known to be slower than deamidation, fragmentation products are often observed during peptide deamidation experiments. In this study, mechanisms leading to the cleavage of the carboxyl-side peptide bond of Asn and Asp residues were investigated using computational methods (B3LYP/6-31+G**). Single-point solvent calculations at the B3LYP/6-31++G** level were carried out in water, utilizing the integral equation formalism-polarizable continuum (IEF-PCM) model. Mechanism and energetics of peptide fragmentation at Asn were comparatively analyzed with previous calculations on deamidation of Asn. When deamidation proceeds through direct hydrolysis of the Asn side chain or through cyclic imide formationvia a tautomerization routeit exhibits lower activation barriers than peptide bond cleavage at Asn. The fundamental distinction between the mechanisms leading to deamidationvia a succinimideand backbone cleavage was found to be the difference in nucleophilic entities involved in the cyclization process (backbone versus side-chain amide nitrogen). If deamidation is prevented by protein three-dimensional structure, cleavage may become a competing pathway. Fragmentation of the peptide backbone at Asp was also computationally studied to understand the likelihood of Asn deamidation preceding backbone cleavage. The activation barrier for backbone cleavage at Asp residues is much lower (approximately 10 kcal/mol) than that at Asn. This suggests that peptide bond cleavage at Asn residues is more likely to take place after it has deamidated into Asp.     

Effects of Substitutions D73E, D73N, D103N and V106M on Signaling and pH Titration of Sensory Rhodopsin II

Abstract— Several mutations in the repellent phototaxis receptor sensory rhodopsin II (SRII), in residues homologous to residues important in the related proton pump bacteriorhodopsin, were expressed in Pho81Wr^–, a Halobacterium salinarum strain deficient in production of SRII and its transducer protein HtrII. The lack of production of SRII and HtrII is shown to be due to insertion of an ISH2 transposon into the promoter region upstream of the htrII - sopH gene pair. Near wild-type phototaxis responses are rescued in Pho81Wr^– by expression of HtrII with D73E, D103N or V106M receptors. Partial responses are restored by the HtrII-D73N pair. From absorption spectroscopy of his-tag-purified receptor protein from mutants D73N and D73E we conclude that Asp73 is the primary counterion to the protonated Schiff base in SRII, like the corresponding Asp85 in bacteriorhodopsin. The absorption maximum of SRII (487 nm) is shifted to 514 nm in mutant D73N, a 1080 cm⁻¹ shift identical to that caused by D85N in bacteriorhodopsin. Acid titration of SRII also induces the red shift with a pK of 3.0 in wild type. The absorption shift and the pK are nearly the same in V106M and D103N, but the pK is raised to 5.1 in D73E, confirming that Asp73 is the residue responsible for this spectral transition.     

Influence of the distal his in imparting imidazolate character to the proximal his in heme peroxidase: (1)h NMR spectroscopic study of cyanide-inhibited his42-->ala horseradish peroxidase

The functional higher oxidation states of heme peroxidases have been proposed to be stabilized by the significant imidazolate character of the proximal His. This is induced by a "push-pull" combination effect produced by the proximal Asp that abstracts ("pulls") the axial His ring N(delta)H, along with the distal protonated His that contributes ("pushes") a strong hydrogen bond to the distal ligand. The molecular and electronic structure of the distal His mutant of cyanide-inhibited horseradish peroxidase, H42A-HRPCN, has been investigated by NMR. This complex is a valid model for the active site hydrogen-bonding network of HRP compound II. The (1)H and (15)N NMR spectral parameters characterize the relative roles of the distal His42 and proximal Asp247 in imparting imidazolate character to the axial His. 1D/2D spectra reveal a heme pocket molecular structure that is highly conserved in the mutant, except for residues in the immediate proximity of the mutation. This conserved structure, together with the observed dipolar shifts of numerous active site residue protons, allowed a quantitative determination of the orientation and anisotropies of the paramagnetic susceptibility tensor, both of which are only minimally perturbed relative to wild-type HRPCN. The quantitated dipolar shifts allowed the factoring of the hyperfine shifts to reveal that the significant changes in hyperfine shifts for the axial His and ligated (15)N-cyanide result primarily from changes in contact shifts that reflect an approximately one-third reduction in the axial His imidazolate character upon abolishing the distal hydrogen-bond to the ligated cyanide. Significant changes in side chain orientation were found for the distal Arg38, whose terminus reorients to partially fill the void left by the substituted His42 side chain. It is concluded that 1D/2D NMR can quantitate both molecular and electronic structural changes in cyanide-inhibited heme peroxidase and that, while both residues contribute, the proximal Asp247 is more important than the distal His42 in imparting imidazole character to the axial His 170.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.