Molecular dynamics simulations of the intramolecular proton transfer and carbanion stabilization in the pyridoxal 5'-phosphate dependent enzymes L-dopa decarboxylase and alanine racemase

Molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential have been carried out to investigate the internal proton transfer equilibrium of the external aldimine species in l-dopa decarboxylase, and carbanion stabilization by the enzyme cofactor in the active site of alanine racemase. Solvent effects lower the free energy of the O-protonated PLP tautomer both in aqueous solution and in the active site, resulting a free energy difference of about -1 kcal/mol relative to the N-protonated Schiff base in the enzyme. The external aldimine provides the dominant contribution to lowering the free energy barrier for the spontaneous decarboxylation of l-dopa in water, by a remarkable 16 kcal/mol, while the enzyme l-dopa decarboxylase further lowers the barrier by 8 kcal/mol. Kinetic isotope effects were also determined using a path integral free energy perturbation theory on the primary (13)C and the secondary (2)H substitutions. In the case of alanine racemase, if the pyridine ring is unprotonated as that in the active site, there is destabilizing contribution to the formation of the α-carbanion in the gas phase, although when the pyridine ring is protonated the contribution is stabilizing. In aqueous solution and in alanine racemase, the α-carbanion is stabilized both when the pyridine ring is protonated and unprotonated. The computational studies illustrated in this article show that combined QM/MM simulations can help provide a deeper understanding of the mechanisms of PLP-dependent enzymes. This article is part of a Special Issue entitled: Pyridoxal Phosphate Enzymology.     

Molecular modeling of the mechanochemical triggering mechanism for catalysis of carbon-cobalt bond homolysis in coenzyme B12

The possible contributions of the mechanochemical triggering effect to the enzymatic activation of the carbon-cobalt bond of coenzyme B12 (5'-deoxyadenosylcobalamin, AdoCbl) for homolytic cleavage have been studied by molecular modeling and semiempirical molecular orbital calculations. Classically, this effect has envisioned enzymatic compression of the axial Co-N bond in the ground state to cause upward folding of the corrin ring and subsequent sterically induced distortion of the Co-C bond leading to its destabilization. The models of this process show that in both methylcobalamin (CH3Cbl) and AdoCbl, compression of the axial Co-N bond does engender upward folding of the corrin ring, and that the extent of such upward folding is smaller in an analog in which the normal 5,6-dimethylbenzimidazole axial ligand is replaced by the sterically smaller ligand, imidazole (CH3(lm)Cbl and Ado(lm)Cbl). Furthermore, in AdoCbl, this upward folding of the corrin is accompanied by increases in the carbon-cobalt bond length and in the Co-C-C bond angle (which are also less pronounced in Ado(Im)Cbl), and which indicate that the Co-C bond is indeed destabilized by this mechanism. However, these effects on the Co-C bond are small, and destabilization of this bond by this mechanism is unlikely to contribute more than ca. 3 kcal mol(-1) towards the enzymatic catalysis of Co-C bond homolysis, far short of the observed ca. 14 kcal mol(-1). A second version of mechanochemical triggering, in which compression of the axial Co-N bond in the transition state for Co-C bond homolysis stabilizes the transition state by increased Co-N orbital overlap, has also been investigated. Stretching the Co-C bond to simulate the approach to the transition state was found to result in an upward folding of the corrin ring, a slight decrease in the axial Co-N bond length, a slight displacement of the metal atom from the plane of the equatorial nitrogens towards the "lower" axial ligand, and a decrease in strain energy amounting to about 8 kcal mol(-1) for both AdoCbl and Ado(Im)Cbl. In such modeled transition states, compression of the axial Co-N bond to just below 2.0 A (the distance subsequently found to provide maximal stabilization of the transition state by increased orbital overlap) required about 4 kcal mol(-1) for AdoCbl, and about 2.5 kcal mol(-1) for Ado(Im)Cbl. ZINDO/1 calculations on slightly simplified structures showed that maximal electronic stabilization of the transition state by about 10 kcal mol(-1) occurred at an axial Co-N bond distance of 1.96 A for both AdoCbl and Ado(Im)Cbl. The net result is that this type of transition state mechanochemical triggering can provide 14 kcal mol(-1) of transition state stabilization for AdoCbl, and about 15.5 kcal mol(-1) for the Ado(Im)Cbl, enough to completely explain the observed enzymatic catalysis. These results are discussed in the light of current knowledge about class I AdoCbl-dependent enzymes, in which the coenzyme is bound in its "base-off" conformation, with the lower axial ligand position occupied by the imidazole moiety of an active site histidine residue, and the class II enzymes, in which AdoCbl binds to the enzyme in its "base-on" conformation, and the pendent 5,6-dimethylbenzimidazole base remains coordinated to the metal during Co-C bond activation.     

A definitive mechanism for chorismate mutase

In previous research presentations, we have described the important features of the chorismate --> prephenate reaction using molecular dynamics (MD) and thermodynamic integration studies. This investigation of the reaction in Escherichia coli and water involves QM/MM procedures (SCCDFTB/MM two-dimensional reaction coordinates to identify transition state structures in the water, enzyme, and gas phase followed by B3LYP/6-31+G* single-point computations which allow the determination of activation energies in water and in the E. coli enzyme). Computed activation energies of 11.3 kcal/mol in enzyme and 20.3 kcal/mol in water may be compared to the experimental values of 12.7 and 20.7 kcal/mol, respectively. The transition state structures in the gas phase, water, and enzyme are much the same. The transition states are characteristic of a concerted pericyclic rearrangement. The very small differences in the partial charges of O13 in NAC and TS support only a small preferential (10%) electrostatic stabilization of TS. The free energy of NAC formation in water exceeds that in enzyme by 8.5 kcal/mol, and it is this favored formation of NAC that provides the major kinetic advantage to the enzymatic reaction. These findings compare most favorably with those previous observations of this laboratory employing molecular dynamics and thermodynamic integrations. A definitive mechanism for the chorismate mutase enzymes is provided.     

Nitrogen isotope effects on glutamate decarboxylase from Escherichia coli

The nitrogen isotope effect on the decarboxylation of glutamic acid by glutamate decarboxylase from Escherichia coli has been measured by comparison of the isotopic composition of the amino nitrogen of the product gamma-aminobutyric acid isolated after 10-20% reaction with that of the starting glutamic acid. At pH 4.7, 37 degrees C, the isotope effect is k14/k15 = 0.9855 +/- 0.0006 when compared to unprotonated glutamic acid. Interpretation of this result requires knowledge of the equilibrium nitrogen isotope effect for Schiff base formation. This equilibrium isotope effect is k14/k15 = 0.9824 for the formation of the unprotonated Schiff base between unprotonated valine and salicylaldehyde. Analysis of the nitrogen isotope effect on decarboxylation of glutamic acid and of the previously measured carbon isotope effect on this same reaction [O'Leary, M.H., Yamada, H., & Yapp, C.J. (1981) Biochemistry 20, 1476] shows that decarboxylation and Schiff base formation are jointly rate limiting. The enzyme-bound Schiff base between glutamate and pyridoxal 5'-phosphate partitions approximately 2:1 between decarboxylation and return to the starting state. The nitrogen isotope effect also reveals that the Schiff base nitrogen is protonated in this intermediate.     

How does GAP catalyze the GTPase reaction of Ras? A computer simulation study

The formation of a complex between p21(ras) and GAP accelerates the GTPase reaction of p21(ras) and terminates the signal for cell proliferation. The understanding of this rate acceleration is important for the elucidation of the role of Ras mutants in tumor formation. In principle there are two main options for the origin of the effect of GAP. One is a direct electrostatic interaction between the residues of GAP and the transition state of the Ras-GAP complex and the other is a GAP-induced shift of the structure of Ras to a configuration that increases the stabilization of the transition state. This work examines the relative importance of these options by computer simulations of the catalytic effect of Ras. The simulations use the empirical valence bond (EVB) method to study the GTPase reaction along the alternative associative and dissociative paths. This approach reproduces the trend in the overall experimentally observed catalytic effect of GAP: the calculated effect is 7 +/- 3 kcal/mol as compared to the observed effect of approximately 6.6 kcal/mol. Furthermore, the calculated effect of mutating Arg789 to a nonpolar residue is 3-4 kcal/mol as compared to the observed effect of 4.5 kcal/mol for the Arg789Ala mutation. It is concluded, in agreement with previous proposals, that the effect of Arg789 is associated with its direct interaction with the transition state charge distribution. However, calculations that use the coordinates of Ras from the Ras-GAP complex (referred to here as Ras') reproduce a significant catalytic effect relative to the Ras coordinates. This indicates that part of the effect of GAP involves a stabilization of a catalytic configuration of Ras. This configuration increases the positive electrostatic potential on the beta-phosphate (relative to the corresponding situation in the free Ras). In other words, GAP stabilizes the GDP bound configuration of Ras relative to that of the GTP-bound conformation. The elusive oncogenic effect of mutating Gln61 is also explored. The calculated effect of such mutations in the Ras-GAP complex are found to be small, while the observed effect is very large (8.7 kcal/mol). Since the Ras is locked in its Ras-GAP configuration in our simulations, we conclude that the oncogenic effect of mutation of Gln61 is indirect and is associated most probably with the structural changes of Ras upon forming the Ras-GAP complex. In view of these and the results for the Ras' we conclude that GAP activates Ras by both direct electrostatic stabilization of the transition state and an indirect allosteric effect that stabilizes the GDP-bound form. The present study also explored the feasibility of the associative and dissociative mechanism in the GTPase reaction of Ras. It is concluded that the reaction is most likely to involve an associative mechanism.     

Initial step of B12-dependent enzymatic catalysis: energetic implications regarding involvement of the one-electron-reduced form of adenosylcobalamin cofactor

Density functional theory analysis was performed to elucidate the impact of one-electron reduction upon the initial step of adenosylcobalamin-dependent enzymatic catalysis. The transition state (TS) corresponding to the Co–C bond cleavage and subsequent hydrogen abstraction from the substrate was located. The intrinsic reaction coordinate calculations predicted that the reaction consisting of Co–C5′ bond cleavage in [Co^III(corrin^•)]–Rib (where Rib is ribosyl) and hydrogen-atom abstraction from the CH₃–CH₂–CHO substrate occurs in a concerted fashion. The computed activation energy barrier of the reaction (15.0 kcal/mol) was lowered by approximately 54.5% in comparison with the reaction involving the positively charged cofactor model (Im–[Co^III(corrin)]–Rib⁺, where Im is imidazole; energy barrier = 33.0 kcal/mol). The Im base was detached during the TS search in the reaction involving the one-electron-reduced analogue. Thus, to compare the energetics of the two reactions, the axial Im ligand detachment energy for the Im–[Co^III(corrin^•)]–Rib model was computed [7.6 kcal/mol (gas phase); 4.6 kcal/mol (water)]. Consequently, the effective activation energy barrier for the reaction mediated by the Im-off [Co^III(corrin^•)]–Rib was estimated to be 22.6 kcal/mol, which implied an overall 31.5% reduction in the energetic demands of the reaction. Considering that the lengthened Co–N_axial bond has been observed in X-ray crystal structure studies of B₁₂-dependent mutases, the catalytic impact induced by one-electron reduction of the cofactor is expected to be higher in the presence of the enzymatic environment.     

Powering DNA repair through substrate electrostatic interactions

The reaction catalyzed by the DNA repair enzyme uracil DNA glycosylase (UDG) proceeds through an unprecedented stepwise mechanism involving a positively charged oxacarbenium ion sugar and uracil anion leaving group. Here we use a novel approach to evaluate the catalytic contribution of electrostatic interactions between four essential phosphodiester groups of the DNA substrate and the cationic transition state. Our strategy was to substitute each of these phosphate groups with an uncharged (R)- or (S)-methylphosphonate linkage (MeP). We then compared the damaging effects of these methylphosphonate substitutions on catalysis with their damaging effects on binding of a cationic 1-azadeoxyribose (1-aza-dR(+)) oxacarbenium ion analogue to the UDG-uracil anion binary complex. A plot of log k(cat)/K(m) for the series of MeP-substituted substrates against log K(D) for binding of the 1-aza-dR(+) inhibitors gives a linear correlation of unit slope, confirming that the electronic features of the transition state resemble that of the 1-aza-dR(+), and that the anionic backbone of DNA is used in transition state stabilization. We estimate that all of the combined phosphodiester interactions with the substrate contribute 6-8 kcal/mol toward lowering the activation barrier, a stabilization that is significant compared to the 16 kcal/mol catalytic power of UDG. However, unlike groups of the enzyme that selectively stabilize the charged transition state by an estimated 7 kcal/mol, these phosphodiester groups also interact strongly in the ground state. To our knowledge, these results provide the first experimental evidence for electrostatic stabilization of a charged enzymatic transition state and intermediate using the anionic backbone of DNA.     

Closing the gate to the active site: effect of the inhibitor methoxyarachidonyl fluorophosphonate on the conformation and membrane binding of fatty acid amide hydrolase

Fatty acid amide hydrolase (FAAH) is a dimeric, membranebound enzyme that degrades neuromodulatory fatty acid amides and esters and is expressed in mammalian brain and peripheral tissues. The cleavage of approximately 30 amino acids from each subunit creates an FAAH variant that is soluble and homogeneous in detergent-containing buffers, opening the avenue to the in vitro mechanistic and structural studies. Here we have studied the stability of FAAH as a function of guanidinium hydrochloride concentration and of hydrostatic pressure. The unfolding transition was observed to be complex and required a fitting procedure based on a three-state process with a monomeric intermediate. The first transition was characterized by dimer dissociation, with a free energy change of approximately 11 kcal/mol that accounted for approximately 80% of the total stabilization energy. This process was also paralleled by a large change in the solvent-accessible surface area, because of the hydration occurring both at the dimeric interface and within the monomers. As a consequence, the isolated subunits were found to be much less stable (DeltaG approximately 3 kcal/mol). The addition of methoxyarachidonyl fluorophosphonate, an irreversible inhibitor of FAAH activity, enhanced the stability of the dimer by approximately 2 kcal/mol, toward denaturant- and pressure-induced unfolding. FAAH inhibition by methoxyarachidonyl fluorophosphonate also reduced the ability of the protein to bind to the membranes. These findings suggest that local conformational changes at the level of the active site might induce a tighter interaction between the subunits of FAAH, affecting the enzymatic activity and the interaction with membranes.     

Photochemical properties of a bacteriorhodopsin analog with 13-demethyl-13-trifluoromethyl retinal

It was shown that the substitution of the CF₃ group in the structure of retinal for the methyl group at C13 causes not only a decrease in the affinity of the proton for the nitrogen in the Schiff base (pK ～ 8.4) but also considerably changes the photochemical properties of the bacteriorhodopsin analog. At pH > 6.5, the rate of the Schiff base reprotonation during M decay depends on the proton concentration in the medium. In the photocycle of the yellow M-like form with the deprotonated Schiff base, a long-wavelength product absorbing at 625 nm is formed, which has a similar pH dependence of decay kinetics. The two processes also have similar activation energies (about 15 ± 1 kcal/mol). It is concluded that both cases involve proton transfer from an aqueous medium through the donor part of the channel to the Schiff base and Asp96, respectively. In the analog, however, the structure of water molecules necessary for the stabilization of the proton on the Schiff base is broken. As a result, dehydration of the preparation gives rise to a fraction of M-like form of bacteriorhodopsin with the deprotonated Schiff base.     

An investigation of the pyranose ring interconversion path of alpha-L-idose calculated using density functional theory

The interconversion pathways of the pyranose ring conformation of alpha-L-idose from a (4)C1 chair to other conformations were investigated using density functional calculations. From these calculations, four different ring interconversion paths and their transition state structures from the (4)C1 chair to other conformations, such as B(3,O), and (1)S3, were obtained. These four transition-state conformations cover four possible combinations of the network patterns of the hydroxyl group hydrogen bonds (clockwise and counterclockwise) and the conformations of the primary alcohol group (tg and gg). The optimized conformations, transition states, and their intrinsic reaction coordinates (IRC) were all calculated at the B3LYP/6-31G** level. The energy differences among the structures obtained were evaluated at the B3LYP/6-311++G** level. The optimized conformations indicate that the conformers of (4)C1, (2)S(O), and B(3,O) have similar energies, while (1)S3 has a higher energy than the others. The comparison of the four transition states and their ring interconversion paths, which were confirmed using the IRC calculation, suggests that the most plausible ring interconversion of the alpha-L-idopyranose ring occurs between (4)C1 and B(3,O) through the E3 envelope, which involves a 5.21 kcal/mol energy barrier.     

Thermodynamic and base-pairing studies of matched and mismatched DNA dodecamer duplexes containing cis-syn, (6-4) and Dewar photoproducts of TT.

Cis-syn dimers, (6-4) products and their Dewar valence isomers are the major photoproducts of DNA and have different mutagenic properties and rates of repair. To begin to understand the physical basis for these differences, the thermal stability and base pairing properties of the corresponding photoproducts of the TT site in d(GAGTATTATGAG) were investigated. The (6-4) and Dewar products destabilize the duplex form by approximately 6 kcal/mol of free energy at 37 degreesC relative to the parent, whereas a cis-syn dimer only destabilizes the duplex form by 1.5 kcal/mol. Duplexes with G opposite the 3'-T of the (6-4) and Dewar products are more stable than those with A by approximately 0.4 kcal/mol, whereas the cis-syn dimer prefers A over G by 0.7 kcal/mol. Proton NMR suggests that wobble base pairing takes place between the 3'-T of the cis-syn dimer and an opposed G, whereas there is no evidence of significant H-bonding between these two bases in the (6-4) product. The thermodynamic and H-bonding data for the (6-4) product are consistent with a 4 nt interior loop structure which may facilitate flipping of the photoproduct in and out of the helix.     

Arginine coordination in enzymatic phosphoryl transfer: evaluation of the effect of Arg166 mutations in Escherichia coli alkaline phosphatase

Arginine residues are commonly found in the active sites of enzymes catalyzing phosphoryl transfer reactions. Numerous site-directed mutagenesis experiments establish the importance of these residues for efficient catalysis, but their role in catalysis is not clear. To examine the role of arginine residues in the phosphoryl transfer reaction, we have measured the consequences of mutations to arginine 166 in Escherichia coli alkaline phosphatase on hydrolysis of ethyl phosphate, on individual reaction steps in the hydrolysis of the covalent enzyme-phosphoryl intermediate, and on thio substitution effects. The results show that the role of the arginine side chain extends beyond its positive charge, as the Arg166Lys mutant is as compromised in activity as Arg166Ser. Through measurement of individual reaction steps, we construct a free energy profile for the hydrolysis of the enzyme-phosphate intermediate. This analysis indicates that the arginine side chain strengthens binding by approximately 3 kcal/mol and provides an additional 1-2 kcal/mol stabilization of the chemical transition state. A 2.1 A X-ray diffraction structure of Arg166Ser AP is presented, which shows little difference in enzyme structure compared to the wild-type enzyme but shows a significant reorientation of the bound phosphate. Altogether, these results support a model in which the arginine contributes to catalysis through binding interactions and through additional transition state stabilization that may arise from complementarity of the guanidinum group to the geometry of the trigonal bipyramidal transition state.     

Cis-trans isomerization and puckering of proline residue

We report here the results on N-acetyl-L-proline-N'-methylamide (Ac-Pro-NHMe) calculated at the HF/6-31+G(d) level with the conductor-like polarizable continuum model (CPCM) of self-consistent reaction field methods to investigate the changes of backbone and prolyl ring along the cis-trans isomerization of the prolyl peptide bond. From the potential energy surface, the barrier to ring flip from the down-puckered conformation to the up-puckered one is estimated to be 2.5 and 3.2 kcal/mol for trans and cis conformers of Ac-Pro-NHMe, respectively. In particular, the ring flip seems to be inaccessible in the intermediate regions between trans and cis conformations, because of higher barriers (approximately 13-19 kcal/mol) to rotation of the prolyl peptide bond. The torsion angles for backbone and prolyl ring vary largely around the transition states at omega' approximately 120 degrees and -70 degrees for the prolyl peptide bond. Three kinds of puckering amplitudes show the same trend of puckering along the cis-trans isomerization although their absolute values are different. In particular, trans and cis conformations have the almost same degree of puckering. The cis populations and barriers to rotation of the prolyl peptide bond for Ac-Pro-NHMe are increased with the increase of solvent polarity, which is mainly ascribed to the decreases of relative free energies for cis conformations and the increase of relative free energies for transition states.     

The PLP cofactor: lessons from studies on model reactions

Experimental probes of the acidity of weak carbon acids have been developed and used to determine the carbon acid pK(a)s of glycine, glycine derivatives and iminium ion adducts of glycine to the carbonyl group, including 5'-deoxypyridoxal (DPL). The high reactivity of the DPL-stabilized glycyl carbanion towards nucleophilic addition to both DPL and the glycine-DPL iminium ion favors the formation of Claisen condensation products at enzyme active sites. The formation of the iminium ion between glycine and DPL is accompanied by a 12-unit decrease in the pK(a) of 29 for glycine. The complicated effects of formation of glycine iminium ions to DPL and other aromatic and aliphatic aldehydes and ketones on carbon acid pK(a) are discussed. These data provide insight into the contribution of the individual pyridine ring substituents to the catalytic efficiency of DPL. It is suggested that the 5'-phosphodianion group of PLP may play an important role in enzymatic catalysis of carbon deprotonation by providing up to 12 kcal/mol of binding energy that is utilized to stabilize the transition state for the enzymatic reaction. This article is part of a Special Issue entitled: Pyridoxal Phospate Enzymology.     

The Effect of Protonation and Electrical Interactions on the Stereochemistry of Retinal Schiff Bases

Based on quantumchemical MNDOC calculations it is shown that the ground-state properties of a retinal Schiff base depend sensitively on its protonation state and charge environment. This is exemplified for the equilibrium geometry, for the distribution of partial charges and, in particular, for the thermal isomerization barriers around the π-bonds. It is demonstrated that a protein, by protonating the retinal Schiff base and by providing one or two negative ions in its environment, can reduce double-bond isomerization barriers from 50 kcal/mol for the unprotonated compound to ~ 5 kcal/mol and can increase single bond barriers from 5 kcal/mol to ~ 20 kcal/mol. Thereby, the specific location of the ions relative to the polyene chain of the protonated retinal Schiff base determines the barrier heights. The results explain the ground-state isomerization reactions of retinal observed in bacteriorhodopsin and in squid retinochrome.     

Evolutionary conservation of a functionally important backbone phosphate group critical for aminoacylation of histidine tRNAs

All histidine tRNA molecules have an extra nucleotide, G-1, at the 5' end of the acceptor stem. In bacteria, archaea, and eukaryotic organelles, G-1 base pairs with C73, while in eukaryotic cytoplasmic tRNAHis, G-1 is opposite A73. Previous studies of Escherichia coli histidyl-tRNA synthetase (HisRS) have demonstrated the importance of the G-1:C73 base pair to tRNAHis identity. Specifically, the 5'-monophosphate of G-1 and the major groove amine of C73 are recognized by E. coli HisRS; these individual atomic groups each contribute approximately 4 kcal/mol to transition state stabilization. In this study, two chemically synthesized 24-nucleotide RNA microhelices, each of which recapitulates the acceptor stem of either E. coli or Saccharomyces cervisiae tRNAHis, were used to facilitate an atomic group "mutagenesis" study of the -1:73 base pair recognition by S. cerevisiae HisRS. Compared with E. coli HisRS, microhelixHis is a much poorer substrate relative to full-length tRNAHis for the yeast enzyme. However, the data presented here suggest that, similar to the E. coli system, the 5' monophosphate of yeast tRNA(His) is critical for aminoacylation by yeast HisRS and contributes approximately 3 kcal/mol to transition state stability. The primary role of the unique -1:73 base pair of yeast tRNAHis appears to be to properly position the critical 5' monophosphate for interaction with the yeast enzyme. Our data also suggest that the eukaryotic HisRS/tRNAHis interaction has coevolved to rely less on specific major groove interactions with base atomic groups than the bacterial system.     

Computational Study of a Model System of Enzyme-Mediated [4+2] Cycloaddition Reaction

A possible mechanistic pathway related to an enzyme-catalyzed [4+2] cycloaddition reac-tion was studied by theoretical calculations at density functional (B3LYP, O3LYP, M062X) and semiempirical levels (PM6-DH2, PM6) performed on a model system. The calculations were carried out for the key [4+2] cycloaddition step considering enzyme-catalyzed biosynthesis of Spinosyn A in a model reaction, where a reliable example of a biological Diels-Alder reaction was reported experimentally. In the present study it was demonstrated that the [4+2] cycloaddition reaction may benefit from moving along the energetically balanced reaction coordinate, which enabled the catalytic rate enhancement of the [4+2] cycloaddition pathway involving a single transition state. Modeling of such a system with coordination of three amino acids indicated a reliable decrease of activation energy by ~18.0 kcal/mol as compared to a non-catalytic transformation.     

A DFT study on the mechanism of the organocatalytic synthesis of a benzoxazine-substituted indolizine derivative

A series of theoretical computations were conducted via density functional theory at the B3LYP/6–31 + G(d,p) level to study the mechanism of the organocatalytic synthesis of a benzoxazine-substituted indolizine derivative. Four possible alternative pathways were considered in this work. The calculated results show that the formation of an N-ylide precursor from 4-dimethylaminopyridine (DMAP) is a key step as it provides the necessary nucleophilic centre for the subsequent H-migration and H-elimination processes. The precursor N-ylide and Schiff base isomers with the most favourable activities in the preliminary work were identified theoretically by analysing the reaction mechanism. The synthetic mechanism to obtain the indolizine derivative was found to be a two-step reaction, with the rate-determining step being the first H migration to form a transition state with a four-membered ring. The catalytic activity of DMAP in the first H-migration step in the overall synthetic process greatly reduces the reaction barrier height. The chiral selectivity of the synthesis is dominated by the spatial geometry of the Schiff base functional group.     

Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin

In dark-adapted bacteriorhodopsin (bR) the retinal moiety populates two conformers: all-trans and (13,15)cis. Here we examine factors influencing the thermodynamic equilibrium and conformational transition between the two forms, using molecular mechanics and dynamics calculations. Adiabatic potential energy mapping indicates that whereas the twofold intrinsic torsional potentials of the C13==C14 and C15==N16 double bonds favor a sequential torsional pathway, the protein environment favors a concerted, bicycle-pedal mechanism. Which of these two pathways will actually occur in bR depends on the as yet unknown relative weight of the intrinsic and environmental effects. The free energy difference between the conformers was computed for wild-type and modified bR, using molecular dynamics simulation. In the wild-type protein the free energy of the (13,15)cis retinal form is calculated to be 1.1 kcal/mol lower than the all-trans retinal form, a value within approximately kBT of experiment. In contrast, in isolated retinal the free energy of the all-trans state is calculated to be 2.1 kcal/mol lower than (13,15)cis. The free energy differences are similar to the adiabatic potential energy differences in the various systems examined, consistent with an essentially enthalpic origin. The stabilization of the (13,15)cis form in bR relative to the isolated retinal molecule is found to originate from improved protein-protein interactions. Removing internal water molecules near the Schiff base strongly stabilizes the (13,15)cis form, whereas a double mutation that removes negative charges in the retinal pocket (Asp85 to Ala; Asp212 to Ala) has the opposite effect.