A QM/MM Investigation of the Catalytic Mechanism of Metal‐Ion‐Independent Core 2 β1,6‐N‐Acetylglucosaminyltransferase

β1,6‐GlcNAc‐transferase (C2GnT) is an important controlling factor of biological functions for many glycoproteins and its activity has been found to be altered in breast, colon, and lung cancer cells, in leukemia cells, in the lymhomonocytes of multiple sclerosis patients, leukocytes from diabetes patients, and in conditions causing an immune deficiency. The result of the action of C2GnT is the core 2 structure that is essential for the further elongation of the carbohydrate chains of O‐glycans. The catalytic mechanism of this metal‐ion‐independent glycosyltransferase is of paramount importance and is investigated here by using quantum mechanical (QM) (density functional theory (DFT))/molecular modeling (MM) methods with different levels of theory. The structural model of the reaction site used in this report is based on the crystal structures of C2GnT. The entire enzyme–substrate system was subdivided into two different subsystems: the QM subsystem containing 206 atoms and the MM region containing 5914 atoms. Three predefined reaction coordinates were employed to investigate the catalytic mechanism. The calculated potential energy surfaces discovered the existence of a concerted S_N2‐like mechanism. In this mechanism, a nucleophilic attack by O6 facilitated by proton transfer to the catalytic base and the separation of the leaving group all occur almost simultaneously. The transition state for the proposed reaction mechanism at the M06‐2X/6‐31G** (with diffuse functions on the O1′, O5′, O_Glu, and O6 atoms) level was located at C1? O6=1.74 Å and C1? O1=2.86 Å. The activation energy for this mechanism was estimated to be between 20 and 29 kcal mol^?1, depending on the method used. These calculations also identified a low‐barrier hydrogen bond between the nucleophile O6H and the catalytic base Glu320, and a hydrogen bond between the N‐acetamino group and the glycosidic oxygen of the donor in the TS. It is proposed that these interactions contribute to a stabilization of TS and participate in the catalytic mechanism.     

Understanding Bacterial Bioluminescence: A Theoretical Study of the Entire Process,from Reduced Flavin to Light Emission

Bacterial bioluminescence (BL) has been successfully applied in water‐quality monitoring and in vivo imaging. The attention of researchers has been attracted for several decades, but the mechanism of bacterial BL is still largely unknown due to the complexity of the multistep reaction process. Debates mainly focus on three key questions: How is the bioluminophore produced? What is the exact chemical form of the bioluminophore? How does the protein environment affect the light emission? Using quantum mechanics (QM), combined QM and molecular mechanics (QM/MM) and molecular dynamic (MD) calculations in gas‐phase, solvent and protein environments, the entire process of bacterial BL was investigated, from flavin reduction to light emission. This investigation revealed that: 1) the chemiluminescent decomposition of flavin peroxyhemiacetal does not occur through the intramolecular chemical initiated electron exchange luminescence (CIEEL) or the “dioxirane” mechanism, as suggested in the literature. Instead, the decomposition occurs according to the charge‐transfer initiated luminescence (CTIL) mechanism for the thermolysis of dioxetanone. 2) The first excited state of 4a‐hydroxy‐4a,5‐dihydroFMN (HFOH) was affirmed to be the bioluminophore of bacterial BL. This study provides details regarding the mechanism by which bacterial BL is produced and is helpful in understanding bacterial BL in general.     

A QM/MM Investigation of Thymine Dimer Radical Anion Splitting Catalyzed by DNA Photolyase

On the mend : The repair reaction of the thymine dimer by DNA photolyase (see picture) is studied by hybrid quantum mechanical/molecular mechanical dynamics simulations based on the X‐ray structure of the enzyme–DNA complex. The dynamics of splitting of the thymine dimer radical anion within the DNA photolyase active site is characterized. The model includes the protein environment.

     

Revisiting the Origin of Bacterial Bioluminescence: QM/MM Study on Oxygenation Reaction of Reduced Flavin in Protein

Bacterial bioluminescence is initiated by the oxygenation reaction of reduced flavin mononucleotide in luciferase. This enzymatic oxygenation occurs in a wide range of biological processes including cellular redox metabolism, biocatalysis, biosynthesis and homeostasis. However, little is known about the mechanism of the enzymatic reaction between singlet reduced flavin and triplet oxygen. To explore the enigmatic oxygenation, for the first time, the reaction of reduced flavin anion with oxygen was studied in bacterial luciferase by a combined quantum mechanics and molecular mechanics method as well as molecular dynamics simulation. The calculated results demonstrate that the reaction proceeds via a proton-coupled electron transfer (PCET) pathway, and the essential αHis44 acts as a catalytic acid to provide the proton. The currently proposed PCET mechanism clearly describes the initial steps of bacterial bioluminescence, and could be suitable for the other flavin oxygenation reactions in enzymes.     

Binding affinities by alchemical perturbation using QM/MM with a large QM system and polarizable MM model

The most general way to improve the accuracy of binding‐affinity calculations for protein–ligand systems is to use quantum‐mechanical (QM) methods together with rigorous alchemical‐perturbation (AP) methods. We explore this approach by calculating the relative binding free energy of two synthetic disaccharides binding to galectin‐3 at a reasonably high QM level (dispersion‐corrected density functional theory with a triple‐zeta basis set) and with a sufficiently large QM system to include all short‐range interactions with the ligand (744–748 atoms). The rest of the protein is treated as a collection of atomic multipoles (up to quadrupoles) and polarizabilities. Several methods for evaluating the binding free energy from the 3600 QM calculations are investigated in terms of stability and accuracy. In particular, methods using QM calculations only at the endpoints of the transformation are compared with the recently proposed non‐Boltzmann Bennett acceptance ratio (NBB) method that uses QM calculations at several stages of the transformation. Unfortunately, none of the rigorous approaches give sufficient statistical precision. However, a novel approximate method, involving the direct use of QM energies in the Bennett acceptance ratio method, gives similar results as NBB but with better precision, ～3 kJ/mol. The statistical error can be further reduced by performing a greater number of QM calculations. © 2015 Wiley Periodicals, Inc.     

A QM/MM study of cisplatin-DNA oligonucleotides: from simple models to realistic systems

QM/MM calculations were employed to investigate the role of hydrogen bonding and pi stacking in several single- and double-stranded cisplatin-DNA structures. Computed geometrical parameters reproduce experimental structures of cisplatin and its complex with guanine-phosphate-guanine. Following QM/MM optimisation, single-point DFT calculations allowed estimation of intermolecular forces through atoms in molecules (AIM) analysis. Binding energies of platinated single-strand DNA qualitatively agree with myriad experimental and theoretical studies showing that complexes of guanine are stronger than those of adenine. The topology of all studied complexes confirms that platination strongly affects the stability of both single- and double-stranded DNAs: Pt-N-H...X (X = N or O) interactions are ubiquitous in these complexes and account for over 70 % of all H-bonding interactions. The pi stacking is greatly reduced by both mono- and bifunctional complexation: the former causes a loss of about 3-4 kcal mol(-1), whereas the latter leads to more drastic disruption. The effect of platination on Watson-Crick GC is similar to that found in previous studies: major redistribution of energy occurs, but the overall stability is barely affected. The BH&H/AMBER/AIM approach was also used to study platination of a double-stranded DNA octamer d(CCTG*G*TCC)d(GGACCAGG), for which an experimental structure is available. Comparison between theory and experiment is satisfactory, and also reproduces previous DFT-based studies of analogous structures. The effect of platination is similar to that seen in model systems, although the effect on GC pairing was more pronounced. These calculations also reveal weaker, secondary interactions of the form Pt...O and Pt...N, detected in several single- and double-stranded DNA.     

On the Regio- and Stereospecificity of Arachidonic Acid Peroxidation Catalyzed by Mammalian 15-Lypoxygenases: A Combined Molecular Dynamics and QM/MM Study

15-Lipoxygenases (15-LOs) catalyse the peroxidation reaction of arachidonic acid (AA) in mammals with remarkable regio- and stereospecificity. This positional-specific peroxidation is of paramount importance because it determines the nature and biological functions of the final metabolites generated by each LO as a result of the oxidative metabolism of AA. Although several hypotheses have been formulated concerning the regio- and stereospecificity of LOs, the molecular basis of such behaviour is still unclear. Herein, we combined quantum mechanics/molecular mechanics calculations with molecular dynamics simulations of the complete rabbit 15-LO/AA solvated model to examine the most accepted hypotheses for the regio- and stereospecificity of LOs. We have found that the clue to explain this specificity is the oxygen-targeting hypothesis through steric shielding of specific residues (mainly Leu597, Gln548 and Phe175, as well as the AA tail itself). Our deductions are based primarily on the analysis of the energy barrier heights from the oxygen addition reaction profiles.     

Cover Picture: QM/MM Car‐Parrinello Molecular Dynamics Study of the Solvent Effects on the Ground State and on the First Excited Singlet State of Acetone in Water (ChemPhysChem 11/2003)

The cover picture shows …?‥the solvent effect of water upon the lowest‐lying singlet excitation in acetone. The transition, which involves the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), is blue‐shifted in water (sol) with respect to the gas phase (vac), since the HOMO is stronger stabilized in water than the LUMO. The authors calculate the absorption and fluorescence spectra of acetone in water with a hybrid Car–Parrinello quantum chemical/classical molecular dynamics approach and investigate the influence of the solvent. The lower part of the picture shows the excitation energy during a simulation. One configuration with a very high excitation energy (three hydrogen bonds, short C?O bond) and one configuration with a very low excitation energy (two hydrogen bonds, long C?O bond) are shown in detail. Find out more in the article by Rothlisberger and co‐workers on pages pp. 1177–1182.     

New QM/MM implementation of the DFTB3 method in the gromacs package

The approximate density‐functional tight‐binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle–mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free‐energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron‐induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization. © 2015 Wiley Periodicals, Inc.     

TDDFT Study of the Electronic Structure,Absorption and Emission Spectra of the Light Emitters of the Amazing Firefly Bioluminescence and Solvation Effects on the Spectra

The ground and excited state properties of luciferin (LH₂) and oxyluciferin (OxyLH₂), the bioluminescent chemicals in the firefly, have been characterized using density functional theory (DFT) and time dependent DFT (TDDFT) methods. The effects of solvation on the electronic absorption and emission spectra of luciferin and oxyluciferin were predicted with a self‐consistent isodensity polarized continuum model of the solvent using TDDFT. The S₀→S₁ vertical excitation energies in the gas phase and in water were obtained. Optimizations of the excited state geometries permitted the first predictions of the fluorescence spectra for these biologically important molecules. Shifts in both of the absorption and emission spectra on proceeding from the gas phase to aqueous solution were also predicted.     

Bioluminophore and Flavin Mononucleotide Fluorescence Quenching of Bacterial Bioluminescence—A Theoretical Study

Bacterial bioluminescence with continuous glow has been applied to the fields of environmental toxin monitoring, drug screening, and in vivo imaging. Nonetheless, the chemical form of the bacterial bioluminophore is still a bone of contention. Flavin mononucleotide (FMN), one of the light‐emitting products, and 4a‐hydroxy‐5‐hydro flavin mononucleotide (HFOH), an intermediate of the chemical reactions, have both been assumed candidates for the light emitter because they have similar molecular structures and fluorescence wavelengths. The latter is preferred in experiments and was assigned in our previous density functional study. HFOH displays weak fluorescence in solutions, but exhibits strong bioluminescence in the bacterial luciferase. FMN shows the opposite behavior; its fluorescence is quenched when it is bound to the luciferase. This is the first example of flavin fluorescence quenching observed in bioluminescent systems and is merely an observation, both the quenching mechanism and quencher are still unclear. Based on theoretical analysis of high‐level quantum mechanics (QM), combined QM and molecular mechanics (QM/MM), and molecular dynamics (MD), this paper confirms that HFOH in its first singlet excited state is the bioluminophore of bacterial bioluminescence. More importantly, the computational results indicate that Tyr110 in the luciferase quenches the FMN fluorescence via an electron‐transfer mechanism.     

QM/MM modeling of enantioselective pybox-ruthenium- and box-copper-catalyzed cyclopropanation reactions: scope, performance, and applications to ligand design

An extensive comparison of full-QM (B3LYP) and QM/MM (B3LYP:UFF) levels of theory has been made for two enantioselective catalytic systems, namely, Pybox-Ru and Box-Cu complexes, in the cyclopropanation of alkenes (ethylene and styrene) with methyl diazoacetate. The geometries of the key reaction intermediates and transition structures calculated at the QM/MM level are generally in satisfactory agreement with full-QM calculated geometries. More importantly, the relative energies calculated at the QM/MM level are in good agreement with those calculated at the full-QM level in all cases. Furthermore, the QM/MM energies are often in better agreement with the stereoselectivity experimentally observed, and this suggests that QM/MM calculations can be superior to full-QM calculations when subtle differences in inter- and intramolecular interactions are important in determining the selectivity, as is the case in enantioselective catalysis. The predictive value of the model presented is validated by the explanation of the unusual enantioselectivity behavior exhibited by a new bis-oxazoline ligand, the stereogenic centers of which are quaternary carbon atoms.     

Computational Study of the Interactions between Benzene and Crystalline Ice Ih: Ground and Excited States

Ground‐state geometries of benzene on crystalline ice cluster model surfaces (I_h) are investigated. It is found that the binding energies of benzene‐bound ice complexes are sensitive to the dangling features of the binding sites. We used time‐dependent DFT to study the UV spectroscopy of benzene, ice clusters, and benzene–ice complexes, by employing the M06‐2X functional. It is observed that the size of the ice cluster and the dangling features have minor effects on the UV spectral characteristics. Benzene‐mediated electronic excitations of water towards longer wavelengths (above 170 nm) are noted in benzene‐bound ice clusters, where the cross‐section of photon absorption by water is negligible, in good agreement with recent experimental results (Thrower et al., J. Vac. Sci. Technol. A, 2008, 26 , 919–924). The intensities of peaks associated with water excitations in benzene–ice complexes are found to be higher than in isolated ice clusters. The π→π* electronic transition of benzene in benzene–ice complexes undergoes a small redshift compared with the isolated benzene molecule, and this holds for all benzene‐bound ice complexes.     

Theoretical Investigation on the Origin of Yellow‐Green Firefly Bioluminescence by Time‐Dependent Density Functional Theory

The question whether the emitter of yellow‐green firefly bioluminescence is the enol or keto‐constrained form of oxyluciferin (OxyLH₂) still has no definitive answer from experiment or theory. In this study, Arg220, His247, adenosine monophosphate (AMP), Water324, Phe249, Gly343, and Ser349, which make the dominant contributions to color tuning of the fluorescence, are selected to simulate the luciferase (Luc) environment and thus elucidate the origin of firefly bioluminescence. Their respective and compositive effects on OxyLH₂ are considered and the electronic absorption and emission spectra are investigated with B3LYP, B3PW91, and PBE1KCIS methods. Comparing the respective effects in the gas and aqueous phases revealed that the emission transition is prohibited in the gas phase but allowed in the aqueous phase. For the compositive effects, the optimized geometry shows that OxyLH₂ exists in the keto(?1) form when Arg220, His247, AMP, Water324, Phe249, Gly343, and Ser349 are all included in the model. Furthermore, the emission maximum wavelength of keto(?1)+Arg+His+AMP+H₂O+Phe+Gly+Ser is close to the experimental value (560 nm). We conclude that the keto(?1) form of OxyLH₂ is a possible emitter which can produce yellow‐green bioluminescence because of the compositive effects of Arg220, His247, AMP, Water324, Phe249, Gly343, and Ser349 in the luciferase environment. Moreover, AMP may be involved in enolization of the keto(?1) form of OxyLH₂. Water324 is indispensable with respect to the environmental factors around luciferin (LH₂).