Molecular dynamics simulations of the mononuclear zinc-beta-lactamase from Bacillus cereus.

Herein, we report molecular dynamics simulations of the mononuclear form of the Bacillus cereuszinc-beta-lactamase. We studied two different configurations which differ in the presence of a zinc-bound hydroxide or a zinc-bound water and in the protonation state of the essential His210 residue. Contacts of the catalytically important residues (Asp90, His210, Cys168, etc.) with the zinc center are characterized by the MD analyses. The nature of the Zn-OH(2) --> His210 proton transfer pathway connecting the two configurations was studied by means of QM calculations on cluster models while the relative stability of the two configurations was estimated from QM/MM calculations in the enzyme. From these results, a theoretical model for the kinetically active form of the B. cereus metalloenzyme is proposed. Some mechanistic implications and the influence of mutating the Cys168 residue are also discussed.     

CO(2) fixation by Rubisco: computational dissection of the key steps of carboxylation, hydration, and C-C bond cleavage.

Despite intensive experimental and computational studies, some important features of the mechanism of the photosynthetic CO(2)-fixing enzyme, Rubisco, are still not understood. To complement our previous investigation of the first catalytic step, the enolization of D-ribulose-1,5-bisphosphate (King et al., Biochemistry 1998, 44, 15414-15422), we present the first complete computational dissection of subsequent steps of the carboxylation reaction that includes the roles of the central magnesium ion and modeled residues of the active site. We investigated carboxylation, hydration, and C-C bond cleavage using the density functional method and the B3LYP/6-31G(d) level to perform geometry optimizations. The energies were determined by B3LYP/6-311+G(2d,p) single-point calculations. We modeled a fragment of the active site and substrate, taking into account experimental findings that the residues coordinated to the Mg ion, especially the carbamylated Lys-201, play critical roles in this reaction sequence. The carbamate appears to act as a general base, not only for enolization but also for hydration of the beta ketoacid formed by addition of CO(2) and, as well, cleavage of the C2-C3 bond of the hydrate. We show that CO(2) is added directly, without assistance of a Michaelis complex, and that hydration of the resultant beta ketoacid occurs in a separate subsequent step with a discrete transition state. We suggest that two conformations of the hydrate (gem-diol), with different metal coordination, are possible. The step with the highest activation energy during the carboxylation cycle is the C-C bond cleavage. Depending on the conformations of the gem-diol, different pathways are possible for this step. In either case, special arrangements of the metal coordination result in bond breaking occurring at remarkably low activation energies (between 28 and 37 kcal mol(-1)) which might be reduced further in the enzyme environment.     

A DFT study of the mechanism of Ni superoxide dismutase (NiSOD): role of the active site cysteine-6 residue in the oxidative half-reaction

In the present DFT study, the catalytic mechanism of H2O2 formation in the oxidative half-reaction of NiSOD, E-Ni(II) + O2- + 2H+ --> E-Ni(III) + H2O2, has been investigated. The main objective of this study is to investigate the source of two protons required in this half-reaction. The proposed mechanism consists of two steps: superoxide coordination and H2O2 formation. The effect of protonation of Cys6 and the proton donating roles of side chains (S) and backbones (B) of His1, Asp3, Cys6, and Tyr9 residues in these two steps have been studied in detail. For protonated Cys6, superoxide binding generates a Ni(III)-O2H species in a process that is exothermic by 17.4 kcal/mol (in protein environment using the continuum model). From the Ni(III)-O2H species, H2O2 formation occurs through a proton donation by His1 via Tyr9, which relative to the resting position of the enzyme is exothermic by 4.9 kcal/mol. In this pathway, a proton donating role of His1 residue is proposed. However, for unprotonated Cys6, a Ni(II)-O2- species is generated in a process that is exothermic by 11.3 kcal/mol. From the Ni(II)-O2- species, the only feasible pathway for H2O2 formation is through donation of protons by the Tyr9(S)-Asp3(S) pair. The results discussed in this study elucidate the role of the active site residues in the catalytic cycle and provide intricate details of the complex functioning of this enzyme.     

An unusual reaction of the natural compound benaphthamycin B: theoretical study of a model system

The surprising and complex transformation of benaphthamycin B to give quinone 2a is investigated theoretically with a model compound, 1,5-dihydroxy-4-methoxy-2,3-dimethylanthraquinone (3). The detailed study is performed using both DFT and perturbation theory under inclusion of solvent effects. Several individual steps (reduction and hydrolysis, water elimination, ether cleavage, and oxidation) of the proposed reaction cascade calculated at the PCM-MP2/6-31G(d)//B3LYP/6-31G(d) level of theory are presented and discussed. It is shown that the key step, the ether cleavage as an S(N)2 reaction leading to the anthrone 12a, possesses a smaller activation barrier compared to the alternative process yielding 12b. Therefore, the formation of the thermodynamically preferred model quinone 13a is also the kinetically favored pathway: The results of the calculated model reaction should also be valid for benaphthamycin B (1).     

Theoretical study of the amazing firefly bioluminescence: the formation and structures of the light emitters

Reaction mechanisms for the formation of the keto-form of oxyluciferin (OxyLH(2)) from the luciferin of fireflies via a dioxetanone intermediate are predicted using the B3LYP/6-31G theoretical method. The ring opening of a model dioxetanone and the decarboxylation proceed in one step via a singlet diradical transition structure with an activation barrier of 18.1 and an exothermicity of 90.8 kcal/mol. The S(0) --> S(1) vertical excitation energies predicted with time dependent density functional theory, TDDFT B3LYP/6-31+G, for the anionic and neutral forms of OxyLH(2) are in the range of 60 to 80 kcal/mol. These energetic results support the generally accepted theory of chemically initiated electron exchange luminescence (CIEEL). The chemical origin of the multicolor bioluminescence from OxyLH(2) is examined theoretically using the TDDFT B3LYP/6-31+G, ZINDO//B3LYP/6-31+G, and CIS/6-31G methods. A change in color of the light emission upon rotation of the two rings in the S(1) excited state of OxyLH(2) is unlikely because both possible emitters, the planar keto- and enol-forms, are minima on the S(1) potential energy surface. The participation of the enol-forms of OxyLH(2) in bioluminescence is plausible but not required to explain the multicolor emission. According to predictions at the TDDFT B3LYP level, the color of the bioluminescence depends on the polarization of the OxyLH(2) in the microenvironment of the enzyme-OxyLH(2) complex.     

Experimental and theoretical investigations of the loss of amino acid side chains in electron capture dissociation of model peptides

Loss of side chains from different amino acid residues in a model peptide framework of RGGGXGGGR under electron capture dissociation conditions were systematically investigated, where X represents one of the twenty common amino acid residues. The alpha-carbon radical cations initially formed by N-Calpha cleavage of peptide ions were shown to undergo secondary dissociation through losses of even-electron and/or odd-electron side-chain moieties. Among the twenty common amino acid residues studied, thirteen of them were found to lose their characteristic side chains in terms of odd-electron neutral fragments, and nine of them were found to lose even-electron neutral side chains. Several generalized dissociation pathways were proposed and were evaluated theoretically with truncated leucine-containing models using ab initio calculations at B3-PMP2/6-311++G(3df,2p)//B3LYP/6-31++G(d,p) level. Elimination of odd-electron side chain was associated with the initial abstraction of the hydrogen from the alpha-carbon bearing the side chain by the N-terminal alpha-carbon radical. Subsequent formation of alpha-beta carbon-carbon double bond leads to the elimination of the odd-electron side chain. The energy barrier for this reaction pathway was 89 kJmol-1. This reaction pathway was 111 kJmol-1 more favorable than the previously proposed pathway involving the formation of cyclic lactam. Elimination of even-electron side chain was associated with the initial abstraction of the gamma-hydrogen from the side chain by the N-terminal alpha-carbon radical. Subsequent formation of beta-gamma carbon-carbon double bond leads to the elimination of the even-electron side chain and the migration of the radical center to the alpha-carbon. The energy barrier for this fragmentation reaction was found to be 50 kJmol-1.     

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.     

Theoretical models of directional proton molecular transport

The important topic of proton transport through molecular wires is usually associated with the Grotthuss mechanism. In this paper we propose an alternative conductor based on chains of lone pairs. B3LYP/6-31+G** and PW91 DFT calculations on model compounds (1,2,3,4-tetrasubstituted benzenes) show that these compounds could play the role of proton conductors.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Quantum chemical study of ester aminolysis catalyzed by a single adenine: a reference reaction for the ribosomal peptide synthesis

Herein we report the results of a HF/6-31+G** and B3LYP/6-31+G** computational investigation of the title reaction for the peptide bond synthesis catalyzed by a single adenine. This system constitutes a reference reaction to study the basic chemical events that have been proposed to occur at the peptidyl transferase active site of ribosomes on the basis of structural determinations (Science 2000, 239, 920--931). Thus, the analysis of the geometry, charge distribution, and energetics of the critical structures involved in this title reaction yields insight into the catalytic mode of action of RNA molecules. Our computational results give further support to the hypotheses that the activated nucleotide A2451 in the ribosome acts as a base catalyst and that this role is similar to that of the His residue in the catalytic triad of serine proteases.     

An exhaustive conformational analysis of N-acetyl-L-cysteine-N-methylamide. Identification of the complete set of interconversion pathways on the ab initio and DFT potential energy hypersurface

The full conformational space of N-acetyl-l-cysteine-N-methylamide was explored by ab initio (RHF/ 6-31G(d)) and DFT (B3LYP/6-31G(d)) computations. Multidimensional conformational analysis predicts 81 structures in N-acetyl-l-cysteine-N-methylamide, but only 47 relaxed structures were previously determined at the RHF/3-21G level of theory. These structures were now optimized using RHF/6-31G(d) and B3LYP/6-31G(d) approaches. Seven conformational migrations were observed when recalculated at higher level of theory. Besides these major changes, only smaller conformational shifts were operative for the remaining stationary points. The exploration of the whole conformational space of N-acetyl-l-cysteine-N-methylamide, including the transition-state structures allowing the conformational interconversion among the low-energy forms, was analyzed in this study. Our results offer new insights into the influence of polar side chains on the conformational preferences of peptide structures.     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

Size Effect on the Raman Spectra and Electronic Structure of the Glycine-alanine Oligopeptide Chains

A theoretical study on oligopeptide chains of glycine-alanine by density functional theory(DFT) is given in this paper. Raman spectra of the oligopeptide chains are examined. The geometric structures, frontier orbital, energy gap, atomic charge distribution, density of states and chemical activity of the side chain are studied at the B3LYP/6-31G(d) level. Results show that, with the number of residues increasing, vibrations of typical functional groups present Raman frequency shift, and the energy gap is gradually reduced. The HOMO and LUMO focus on the amino and carboxyl at the ends of oligopeptides. It is helpful for oligopeptides to self-assemble into chains. In addition, different residues(glycine or alanine) at the ends of chains result in the even-odd effect of orbital energy in the growth process. The size effects of physical and chemical properties only exist when the oligopeptides are shorter, and the phenomenon disappeared as the chain continues to grow.     

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.     

Design of radical-resistant amino acid residues: a combined theoretical and experimental investigation

Ab initio calculations have been used to design radical-resistant amino acid residues. Optimized structures of free and protected amino acids and their corresponding alpha-carbon-centered radicals were determined with B3-LYP/6-31G(d). Single-point RMP2/6-31G(d) calculations on these structures were then used to obtain radical stabilization energies, to examine the effect of steric repulsion between the side chains and amide carbonyl groups on the stability of alpha-carbon-centered peptide radicals. Relative to glycine, the destabilization for alanine and valine residues was found to be approximately 9 and 18 kJ mol(-1), respectively, which correlates with the reactivity of analogous amino acid residues in peptides toward hydrogen atom abstraction in conventional free radical reactions. To design amino acid residues that would resist radical reactions, strategies by which the steric effects could be magnified were considered. This resulted in the identification of tert-leucine and 3,3,3-trifluoroalanine as suitable molecules. With these amino acid residues, the destabilization of the alpha-carbon-centered radicals relative to that of glycine is increased substantially to approximately 36 and 41 kJ mol(-1), respectively. The theoretical predictions have been supported by experimental observations: a tert-leucine derivative was shown to be very slow to react with N-bromosuccinimide, while the corresponding trifluoroalanine derivative was found to be inert.     

硝酸氯冰溶胶水解反应过程的计算模拟

用二级微扰(MP2)和密度泛函理论(B3LYP),辅以不同的基组,对硝酸氯在冰表面上水解反应的机理进行了理论计算研究.根据关键部位化学键的松弛效应和关键原子的电荷分布,对冰表面催化的原因进行了深入分析.水分子一方面作为桥,辅助分子间质子发生迁移;另一方面作为连续介质,通过偶极相互作用加快硝酸氯的水解过程.