Apparent NAC effect in chorismate mutase reflects electrostatic transition state stabilization

The catalytic reaction of chorismate mutase (CM) has been the subject of major current attention. Nevertheless, the origin of the catalytic power of CM remains an open question. In particular, it has not been clear whether the enzyme works by providing electrostatic transition state stabilization (TSS), by applying steric strain, or by populating near attack conformation (NAC). The present work explores this issue by a systematic quantitative analysis. The overall catalytic effect is reproduced by the empirical valence bond (EVB) method. In addition, the binding free energy of the ground state and the transition state is evaluated, demonstrating that the enzyme works by TSS. Furthermore, the evaluation of the electrostatic contribution to the reduction of the activation energy establishes that the TSS results from electrostatic effects. It is also found that the apparent NAC effect is not the reason for the catalytic effect but the result of the TSS. It is concluded that in CM as in other enzymes the key catalytic effect is electrostatic TSS. However, since the charge distribution of the transition state and the reactant state is similar, the stabilization of the transition state leads to reduction in the distance between the reacting atoms in the reactant state.     

Contributions of conformational compression and preferential transition state stabilization to the rate enhancement by chorismate mutase

The rate enhancement provided by the chorismate mutase (CM) enzyme for the Claisen rearrangement of chorismate to prephenate has been investigated by application of the concept of near attack conformations (NACs). Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82% and 100% of chorismate conformers were found to be NAC structures in water and in the CM active site, respectively. Consequently, the conversion of non-NACs to NACs does not contribute to the free energy of activation from preorganization of the substrate into NACs. The FEP calculations yielded differences in free energies of activation that well reproduce the experimental data. Additional calculations indicate that the rate enhancement by CM over the aqueous phase results primarily from conformational compression of NACs by the enzyme and that this process is enthalpically controlled. This suggests that preferential stabilization of the transition state in the enzyme environment relative to water plays a secondary role in the catalysis by CM.     

A transition path sampling study of the reaction catalyzed by the enzyme chorismate mutase

The study of the chemical steps in enzyme-catalyzed reactions represents a challenge for molecular simulation techniques. One concern is how to calculate paths for the reaction. Common techniques include the definition of a reaction coordinate in terms of a small set of (normally) geometrical variables or the determination of minimum energy paths on the potential energy surface of the reacting system. Both have disadvantages, the former because it presupposes knowledge of which variables are likely to be important for reaction and the latter because it provides a static picture and dynamical effects are ignored. In this paper, we employ the transition path sampling method developed by Chandler and co-workers, which overcomes some of these limitations. The reaction that we have chosen is the chorismate-mutase-catalyzed conversion of chorismate into prephenate, which has become something of a test case for simulation studies of enzyme mechanisms. We generated an ensemble of approximately 1000 independent transition paths for the reaction in the enzyme and another approximately 500 for the corresponding reaction in solution. A large variety of analyses of these paths was performed, but we have concentrated on characterizing the transition state ensemble, particularly the flexibility of its structures with respect to other ligands of the enzyme and the time evolution of various geometrical and energetic properties as the reaction proceeds. We have also devised an approximate technique for locating transition state structures along the paths.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

C-H Group as proton donor by formation of a weak hydrogen bond

A systematic CNDO/2 study has been carried out on the weak C-H...F, C-H...O, and C-H...N hydrogen bonds. An H-bond formation has been revealed for proton donors with negative excess charge on the hydrogen. The concept of the positive character of the hydrogen as well as some essential features of the hydrogen bond have been discussed. Some useful correlations have been found and the importance of the charge shifts has been emphasized.     

Co-C bond activation in methylmalonyl-CoA mutase by stabilization of the post-homolysis product Co2+ cobalamin

Despite decades of research, the mechanism by which coenzyme B12 (adenosylcobalamin, AdoCbl)-dependent enzymes promote homolytic cleavage of the cofactor's Co-C bond to initiate catalysis has continued to elude researchers. In this work, we utilized magnetic circular dichroism spectroscopy to explore how the electronic structure of the reduced B12 cofactor (i.e., the post-homolysis product Co2+ Cbl) is modulated by the enzyme methylmalonyl-CoA mutase. Our data reveal a fairly uniform stabilization of the Co 3d orbitals relative to the corrin pi/pi*-based molecular orbitals when Co2+ Cbl is bound to the enzyme active site, particularly in the presence of substrate. Contrastingly, our previous studies (Brooks, A. J.; Vlasie, M.; Banerjee, R.; Brunold, T. C. J. Am. Chem. Soc. 2004, 126, 8167-8180.) showed that when AdoCbl is bound to the MMCM active site, no enzymatic perturbation of the Co3+ Cbl electronic structure occurs, even in the presence of substrate (analogues). Collectively, these observations provide direct evidence that enzymatic Co-C bond activation involves stabilization of the post-homolysis product, Co2+ Cbl, rather than destabilization of the Co3+ Cbl "ground" state.     

Charge optimization increases the potency and selectivity of a chorismate mutase inhibitor

The highest affinity inhibitor for chorismate mutases, a conformationally constrained oxabicyclic dicarboxylate transition state analogue, was modified as suggested by computational charge optimization methods. As predicted, replacement of the C10 carboxylate in this molecule with a nitro group yields an even more potent inhibitor of a chorismate mutase from Bacillus subtilis (BsCM), but the magnitude of the improvement (roughly 3-fold, corresponding to a DeltaDeltaG of -0.7 kcal/mol) is substantially lower than the gain of 2-3 kcal/mol binding free energy anticipated for the reduced desolvation penalty upon binding. Experiments with a truncated version of the enzyme show that the flexible C terminus, which was only partially resolved in the crystal structure and hence omitted from the calculations, provides favorable interactions with the C10 group that partially compensate for its desolvation. Although truncation diminishes the affinity of the enzyme for both inhibitors, the nitro derivative binds 1.7 kcal/mol more tightly than the dicarboxylate, in reasonable agreement with the calculations. Significantly, substitution of the C10 carboxylate with a nitro group also enhances the selectivity of inhibition of BsCM relative to a chorismate mutase from Escherichia coli (EcCM), which has a completely different fold and binding pocket, by 10-fold. These results experimentally verify the utility of charge optimization methods for improving interactions between proteins and low-molecular weight ligands.     

Multiple high-level QM/MM reaction paths demonstrate transition-state stabilization in chorismate mutase: correlation of barrier height with transition-state stabilization

Multiple profiles for the reaction from chorismate to prephenate in the enzyme chorismate mutase calculated with hybrid density functional combined quantum mechanics/molecular mechanics methods (B3LYP/6-31G(d)-CHARMM27) agree well with experiment, and provide direct evidence of transition-state stabilization by this important enzyme, which is at the centre of current debates about the nature of enzyme catalysis.     

Silanediols: a new class of hydrogen bond donor catalysts

Silanediols are introduced as a new class of hydrogen bond donor catalysts for the activation of nitroalkenes toward nucleophilic attack. Excellent yields of product are obtained for the conjugate addition of indole to β-nitrostyrene catalyzed with a stable, storable dinaphthyl-derived silanediol. The preparation and structural characterization of a C(2)-symmetric chiral silanediol is also reported along with its ability to catalyze the conjugate addition reaction.     

Fluorescence detection of guanine-adenine transition by a hydrogen bond forming small compound

In combination with abasic site-containing oligodeoxynucleotides, 2-amino-4-oxopteridine (pterin) can selectively recognize guanine base over other nucleobases accompanied by fluorescence quenching, which allows clear detection of a guanine-adenine transition with the naked eye.     

Selective imaging of positively charged polar and nonpolar lipids by optimizing matrix solution composition

Previous studies have shown that matrix‐assisted laser desorption/ionization–imaging mass spectrometry (MALDI‐IMS) is useful for studying the distribution of various small metabolites, particularly lipids. However, in this technique, selective ionization of the target molecules is imperative, particularly when analyzing small molecules. Since the sample clean‐up procedures available for the MALDI‐IMS of small metabolites are limited, the tissue sample will contain numerous molecular species other than the target molecules. These molecules will compete for ionization resulting in severe ion suppression. Hence, it is necessary to develop and optimize a sample preparation protocol for the target molecules. In this study, through model experiments using reference compounds, we optimized the composition of the matrix solution used for positively charged lipids in terms of the concentration of the organic solvent and presence/absence of alkali metal salts. We demonstrated that a high concentration of organic solvent in the matrix solution favors the preferential detection of lipids over peptides. The presence of alkali metal salts in the matrix solution was favorable for the detection of polar lipids, while a salt‐free matrix solution was suitable for the detection of nonpolar lipids. Furthermore, potassium salts added to the matrix solution caused merging of various lipid adducts (adducts with proton, sodium, and potassium) into one single potassiated species. Using the optimized protocols, we selectively analyzed phosphatidylcholine (PC) and triacylglycerol (TG) with different fatty acid compositions in a rat kidney section. Copyright © 2009 John Wiley & Sons, Ltd.     

The hydrogen bond in the solid state

The hydrogen bond is the most important of all directional intermolecular interactions. It is operative in determining molecular conformation, molecular aggregation, and the function of a vast number of chemical systems ranging from inorganic to biological. Research into hydrogen bonds experienced a stagnant period in the 1980s, but re-opened around 1990, and has been in rapid development since then. In terms of modern concepts, the hydrogen bond is understood as a very broad phenomenon, and it is accepted that there are open borders to other effects. There are dozens of different types of X-H.A hydrogen bonds that occur commonly in the condensed phases, and in addition there are innumerable less common ones. Dissociation energies span more than two orders of magnitude (about 0.2-40 kcal mol(-1)). Within this range, the nature of the interaction is not constant, but its electrostatic, covalent, and dispersion contributions vary in their relative weights. The hydrogen bond has broad transition regions that merge continuously with the covalent bond, the van der Waals interaction, the ionic interaction, and also the cation-pi interaction. All hydrogen bonds can be considered as incipient proton transfer reactions, and for strong hydrogen bonds, this reaction can be in a very advanced state. In this review, a coherent survey is given on all these matters.     

Exploring the quantum mechanical/molecular mechanical replica path method: a pathway optimization of the chorismate to prephenate Claisen rearrangement catalyzed by chorismate mutase

 A replica path method has been developed and extended for use in complex systems involving hybrid quantum/classical (quantum mechanical/molecular mechanical) coupled potentials. This method involves the definition of a reaction path via replication of a set of macromolecular atoms. An “important” subset of these replicated atoms is restrained with a penalty function based on weighted root-mean-square rotation/translation best-fit distances between adjacent (i±1) and next adjacent (i±2) pathway steps. An independent subset of the replicated atoms may be treated quantum mechanically using the computational engine Gamess-UK. This treatment can be performed in a highly parallel manner in which many dozens of processors can be efficiently employed. Computed forces may be projected onto a reference pathway and integrated to yield a potential of mean force (PMF). This PMF, which does not suffer from large errors associated with calculated potential-energy differences, is extremely advantageous. As an example, the QM/MM replica path method is applied to the study of the Claisen rearrangement of chorismate to prephenate which is catalyzed by the Bacillus subtilis isolated, chorismate mutase. Results of the QM/MM pathway minimizations yielded an activation enthalpy ΔH ^†† of 14.9 kcal/mol and a reaction enthalpy of −19.5 kcal/mol at the B3LYP/6-31G(d) level of theory. The resultant pathway was compared and contrasted with one obtained using a forced transition approach based on a reaction coordinate constrained repeated walk procedure (ΔH ^†† =20.1 kcal/mol, ΔH _rxn = −20.1 kcal/mol, RHF/4-31G). The optimized replica path results compare favorably to the experimental activation enthalpy of 12.7±0.4 kcal/mol. Received: 16 December 2001 / Accepted: 6 September 2002 / Published online: 8 April 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 22nd National Meeting of the American Chemical Society, 2001. Correspondence to: H.L. Woodcock e-mail: hlwood@ccqc.uga.edu Acknowledgements. The authors thank Eric Billings, Xiongwu Wu, and Stephen Bogusz for helpful discussions and related work. The authors also show grateful appreciation to The National Institutes of Health and The National Science Foundation for support of the current research.     

Multiple-steering QM-MM calculation of the free energy profile in chorismate mutase

A novel technique for computing free energy profiles in enzymatic reactions using the multiple steering molecular dynamics approach in the context of an efficient QM-MM density functional scheme is presented. The conversion reaction of chorismate to prephenate catalyzed by the Bacillus subtilis enzyme chorismate mutase has been chosen as an illustrative example.     

Selective covalent binding of a positively charged water-soluble benzoheterocycle triosmium cluster to single- and double-stranded DNA

The water-soluble triosmium cluster [Os₃(CO)₉(μ-η²-(4-CHO)C₉H₅N)(μ-H)(P(OCH₂CH₂N(CH₃)₃I)₃)] (4) was tested for its reactivity with plasmid DNA. In contrast to the band retardation previously observed with a related series of positively charged clusters, an intensification and retardation of three discrete bands was observed with increasing cluster concentration. In order to further investigate the apparent modification of DNA by 4, its interaction with a 22-oligomer (sequence 5′-AGT TGT GGT GAC TTT CCC AGG C-3′) was examined. Incubation with this oligonucleotide (pH 7.4 in Tris-HCl buffer and 100 mM NaCl) followed by HPLC analysis revealed the formation of three dose dependent products assigned as covalent modifications at three sites of the oligonucleotide. Incubation of 4 with ³²P-ATP labeled oligonucleotide at the 5′-end followed by treatment with piperidine and comparison with the standard Maxam-Gilbert sequencing protocol products revealed only general background cleavage, indicating that the modification products are piperidine labile and suggesting that the modification involved formation of a Schiff base. An alternative approach was then pursued which involved annealing the 4-oligonucleotide products with their complementary strand and treatment of the resulting duplex DNAwith the exonuclease, Exo III. This assay indicated three exonuclease stops, consistent with the three products observed by HPLC whose electrophoretic mobility approximately matched guanine containing fragments when compared with the Maxam-Gilbert sequencing lanes. Reduction of the 4-oligonucleotide products with borohydride reducing agents, followed by treatment with piperidine, resulted in the formation of one product (by HPLC) with the same electrophoretic mobility as the AGTT fragment based on comparison with the Maxam-Gilbert sequencing lanes. This product most likely results from reduction of an initially formed Schiff base adduct (to the corresponding amine) with the guanine of the TGT fragment of the oligonucleotide, and corresponds to the most stable of the three Schiff base adducts detected by HPLC and by incubation with the exonuclease. The other two products are less stable and competitive reduction of the free aldehyde functionality on the cluster in equilibrium with these adducts precludes their detection after treatment with the reducing agents. The formation of the Schiff base adduct is further corroborated by the model reaction of [Os₃(CO)₁₀(μ-η²-(4-CHO)C₉H₅N)(μ-H)] (4′) with acetylated guanine in nonaqueous solvents where disappearance of the aldehyde resonance and the appearance of several new resonances in the 6-9 ppm region of the ¹H NMR of the reaction mixture is noted.     

Selective transition state stabilization via hyperconjugative and conjugative assistance: stereoelectronic concept for copper-free click chemistry

Dissection of stereoelectronic effects in the transition states (TSs) for noncatalyzed azide-alkyne cycloadditions suggests two approaches to selective transition state stabilization in this reaction. First, the formation of both 1,4- and 1,5-isomers is facilitated via hyperconjugative assistance to alkyne bending and C···N bond formation provided by antiperiplanar σ-acceptors at the propargylic carbons. In addition, the 1,5-TS can be stabilized via attractive C-H···F interactions. Although the two effects cannot stabilize the same transition state for the cycloaddition to α,α-difluorocyclooctyne (DIFO), they can act in a complementary, rather than competing, fashion in acyclic alkynes where B3LYP calculations predict up to ～1 million-fold rate increase relative to 2-butyne. This analysis of stereoelectronic effects is complemented by the distortion analysis, which provides another clear evidence of selective TS stabilization. Changes in electrostatic potential along the reaction path revealed that azide polarization may create unfavorable electrostatic interactions (i.e., for the 1,5-regioisomer formation from 1-fluoro-2-butyne and methyl azide). This observation suggests that more reactive azides can be designed via manipulation of charge distribution in the azide moiety. Combination of these effects with the other activation strategies should lead to the rational design of robust acyclic and cyclic alkyne reagents for fast and tunable "click chemistry". Further computational and experimental studies confirmed the generality of the above accelerating effects and compared them with the conjugative TS stabilization by π-acceptors.     

A theoretical study of the hydrogen bond donor capability and co-operative effects in the hydrogen bond complexes of the diaza-aromatic betacarbolines

In this work, we present a quantum mechanical investigation on the hydrogen bond interactions of N(9)-methyl-9H-pyrido[3,4-b]indole, MBC, and N(2)-methyl-9H-pyrido[3,4-b]indole, BCA, with different hydrogen bond donors. Thus, it has been analysed the influence that the hydrogen bond donor strength and the co-operative effect of the increasing number of donor molecules have on the shape of the potential energy surfaces versus the N···H distances, r(N–H). To rationalize the nature of the interactions, the Bader theory has been applied and the characteristics of the bond critical points analysed. The results show that two different hydrogen bond complexes can be formed depending on the donor capabilities or the number of donor molecules included in the calculations. The topological parameters from the Bader theory are used to justify the statement that the analysed interactions can be classified as weak or partially covalent hydrogen bond interactions, respectively. As experimentally observed, weak hydrogen bond donors form weak hydrogen bond complexes, called HBC. Upon the increase of the donor strength the N···H proton is shifted nearest to the nitrogen atom giving rise to the observation of a stronger hydrogen bond complex, the proton transfer complex, PTC. The most outstanding result of these studies is the fact that the formation of the PTC can also be managed just by changing the number of donor molecules, that is, by a co-operative effect of the hydrogen bonds.     

Exclusive transition state stabilization in the supramolecular catalysis of Diels-Alder reaction by a uranyl salophen complex

Whereas the parent uranyl salophen is catalytically inactive, its phenyl derivative effectively catalyses with turnover the reaction of benzoquinone with 1,3-cyclohexadiene, while showing no appreciable affinity towards reactants and product.     

Selective stabilization of transition state structures for cope rearrangements of semibullvalene and barbaralane through interactions with halogens

A new means-discovered through quantum chemical studies (B3LYP/6-31G(d))-of selectively stabilizing the delocalized forms of unsubstituted semibullvalene and barbaralane is described. Calculations indicate that complexation of delocalized semibullvalene and barbaralane Cope rearrangement transition structures by halogen-containing molecules (homo- and heteronuclear diatomics and Lewis acids) can reduce the barrier significantly, in many cases causing it to disappear entirely.