<Emphasis Type="Italic">Ab initio</Emphasis> fragment molecular orbital studies of influenza virus hemagglutinin–sialosaccharide complexes toward chemical clarification about the virus host range determination

If we predict the host range of new or mutant influenza virus in advance, we are able to measure against pandemic human influenza immediately after the new virus emerges somewhere. Influenza viral hemagglutinin(HA)–sialoside receptor interaction is a target event for in silico chemical prediction studies about the virus host range determination. We theoretically studied avian and human influenza A virus HA H3 subtype complexed with avian or human type receptor Neu5Acα(2-3 or 2-6)Gal analogues by ab initio fragment molecular orbital (FMO) method at the second order Møller–Plesset (MP2)/6–31G level, which can evaluate correctly not only electrostatic interactions but also lipophilic interactions based on van der Waals dispersion force. Avian H3 bound to avian α2-3 11.4 kcal/mol stronger than to human α2-6 in the model complexes with taking account of intermolecular lipophilic interaction. A substitution at the position 226 between Gln(avian) and Leu(human) on influenza H3 HA1 has altered its virus host range between avian and human. In the ab initio FMO studies, binding energy of avian Gln226Leu H3–human α2-6 was quite similar to that in the human H3–human α2-6 complex with amino acid sequence differences at nine positions in the models. This similarity indicates that avian Gln226Leu H3 virus can infect human with the same level as human H3 virus. Opposite mutation Leu226Gln in the human H3 gave the moderate binding energies to avian α2-3 with similarity to avian H3–α2-3 complex that supported our previous virus-sialoside binding assay. Ab initio FMO studies have revealed the relationship between influenza H3 virus host range and H3–α(2-3 or 2-6) receptors binding. Our theoretical approach may predict the infectious level of new viruses and point out some unknown dangerous mutation positions on HA in advance.     

Conformational behaviour of the architectural units of peptides and proteins. Assessment of current understanding by ab initio quantum mechanical methods and refinement of the dipeptide energy surface.

The conformational energy surfaces of analogues of the dipeptide unit of polypeptides and proteins are calculated by ab initio methods using extended basis sets.The calculations are not particularly sensitive to the choice of (extended) basis set.The calculations are shown to support a particular empirical method parameterized with respect to crystal data. Non-hydrogen bonded conformations agree to within 3 kcal mol^?1, even for conformations in which quite considerable degrees of atomic overlap occur.Hydrogen bonded conformations, are, however, in less satisfactory agreement and it is the ab initio calculations which appear to be at fault.A simple correction is applied to the ab initio energy for hydrogen bonded conformations, and with the use of the empirical energy surface a full quantum mechanical conformational energy map is interpolated for the alanyl dipeptide.The effect of flexibility in the peptide backbone is taken into account, and supports recent empirical findings that distortions in valence angles must be considered in calculations of the conformational behaviour of peptides.     

Mapping Enzymatic Catalysis Using the Effective Fragment Molecular Orbital Method: Towards all ab initio Biochemistry

We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol⁻¹ for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively.     

Application of ab initio molecular-orbital calculations to the structural moieties of carbohydrates. Part VI

Ab initio RHF/4–31G molecular-orbital calculations have been conducted on methoxymethyl formate and methoxymethyl acetate as models for examining the anomeric effect and stereochemistry of 1-O-acetylglycopyranoses. The results indicate that, as with the methyl glycopyranosides, the α-⁴C₁(D) configurations are more stable than the β-⁴C₁(D), except that the energy difference is more dependent on the disposition about the glycosidic bond. The lowest-energy conformations occur with glycosidic torsion-angles of ?  180°, where the anomeric energy is about 4 kcal/mol. There is a secondary energy-minimum at ?  90°, for which the anomeric energy is less, about 2 kcal/mol. This orientation corresponds to the conformation most commonly observed in the crystal structures of peracetylated glycopyranoses. Small differences in the CO single-bond lengths, which are observed experimentally in both the α and β anomers, are reproduced by the theoretical calculations.     

Molecular recognition mechanism of FK506 binding protein: an all-electron fragment molecular orbital study

The fragment molecular orbital (FMO) method has enabled electronic structure calculations and geometry optimizations of very large molecules with ab initio quality. We applied the method to four FK506 binding protein (FKBP) complexes (denoted by their PDB codes 1fkb, 1fkf, 1fkg, and 1fki) containing rapamycin, FK506, and two synthetic ligands. The geometries of reduced complex models were optimized at the restricted Hartree–Fock (FMO‐RHF) level using the 3‐21G basis set, and then for a better estimate of binding, the energetics were refined at a higher level of theory (2nd order Møller–Plesset perturbation theory FMO‐MP2 with the 6‐31G* basis set). Thus, obtained binding energies were ?103.9 (?82.0), ?102.2 (?69.2), ?70.1 (?57.7), and ?71.3 (?55.3) kcal/mol for 1fkb, 1fkf, 1fkg, and 1fki, respectively, where the correlation contribution is given in parentheses. The results show that the electron correlation contribution to binding is extremely important, and it accounts for 70–80% of the binding energy. The molecular recognition mechanism of FKBP was analyzed in detail based on the FMO‐pair interactions between protein residues and the ligands. Solvation effects on the protein–ligand binding were estimated using the Poisson–Boltzmann/surface area model. Proteins 2007. © 2007 Wiley‐Liss, Inc.     

Quantum chemical studies of the catalytic mechanism of N-terminal nucleophile hydrolase

Modeling of the catalytic mechanism of penicillin acylase, a member of the N-terminal nucleophile hydrolase superfamily, is for the first time conducted at ab initio quantum chemistry level. The uniqueness of this family of enzymes is that their active site lacks His and Asp (Glu) residues, comprising together with a serine residue the classical catalytic triad. The current investigation confirms that the amino group of the N-terminal serine residue in N-terminal hydrolases is capable of activating its own hydroxyl group. Using the MP2/RHF method with the 6−31+G** basis set, stationary points on the potential energy surface of the considered molecular system were located, corresponding to local minima (complexes of reagents, products, intermediate) and to saddle points (transition states). It turned out that the stage of acyl-serine formation proceeds via two transition states; the first one, which separates reagents from the so-called tetrahedral intermediate, has the highest relative energy (30 kcal/mol). In contrast to recently proposed empiric suggestions, we have found that participation of a bridging water molecule in proton shuttling is not necessary for the catalysis. The quantum chemical calculations showed a crucial role of a specific solvation in decreasing the activation barrier of the reaction by approximately 10 kcal/mol. Published in Russion in Biokhimiya, 2007, Vol. 72, No. 5, pp. 615–621.     

Proton migration in a stack of tyrosine molecules induced by Mg(2+) ion and electric field. Mathematical model

Quantum chemistry methods (ab initio, RHF + MP2(FULL), 6-31G** basis set) were used to study proton migration in tyrosine stacks mimicking the proton channel in tubulin and other proteins. When bound to guanosine-5'-triphosphate, Mg2+ favors the dissociation of water in its first coordination shell, thus initiating subsequent proton shifts in the tyrosine chain composed of spatially remote tyrosine residues of tubulin. The process appears to be thermodynamically allowed, delta G298 < 0, with a potential barrier along the proton shift of no more than 0.75 kcal/mol. The exposure to external electrical field of low intensity, which simulates the electric properties of tubulin, promotes proton migration over long distances.     

A Comparative Semiempirical, Ab initio and DFT Study of Decarbonylation of Cyclic Ketones Part 1: Thermal Fragmentation of Cyclopropanone

Calculations using different quantum mechanical methods including semiempirical (MNDO,AM1 and PM3), ab initio (RHF and MP2 calculations using the 6-311G and 6-311++G^** basis sets), and density functional theory (LSDA, BP, MIXBP and B3LYP, i.e., B3LYP/6-311+G**//B3LYP/6-31G*) have been performed on the thermal fragmentation of cyclopropanone to ethylene and carbon monoxide. All RHF calculations predict a concerted single step mechanism for this conversion. The estimated activation energies vary from 34.4 to 54.6 kcal·mol^-1, mainly localized around 37±2 kcal·mol^-1, depending on the method. Whereas the calculated RHF reaction energies also varied from 14.5 to -33.3 kcal·mol^-1, the B3LYP/6-311+G**//B3LYP/6-31G* method predicts the experimental value (-17.7 kcal·mol^-1) within experimental uncertainties. Remarkably, semiempirical AM1 and PM3 methods and simple DFT calculations, LSDA, predict comparable results to the more advanced methods. UHF ab initio calculations predict the same single step mechanism, whereas a multistep biradical mechanism with an unrealistically low activation energy is favored by the semiempirical methods. Structures of the activated complex of the single step mechanism, estimated by different methods, are very similar and consistent with a nonlinear cheletropic [_2s + _2a] reaction, as predicted by the orbital symmetry rules and earlier EHT calculations.Electronic Supplementary Material available.     

De(side chain) model of epothilone: bioconformer interconversions DFT study

Using ab initio methods, we have studied conformations of the de(sidechain)de(dioxy)difluoroepothilone model to quantify the effect of stability change between the exo and endo conformers of the epoxy ring. The DFT minimization of the macrolactone ring reveals four low energy conformers, although MP2 predicted five stable structures. The model tested with DFT hybride functional (B3LYP/6–31+G(d,p)) exhibits the global minimum for one of the exo forms (C), experimentally observed in the solid state, but unexpectedly with the MP2 electron correlation method for the virtual endo form (W). Using the QST3 technique, several pathways were found for the conversion of the low energy conformers to the other low energy exo representatives, as well as within the endo analog subset. The potential energy relationships obtained for several exo forms suggest a high conformational mobility between three, experimentally observed, conformers. The high rotational barrier, however, excludes direct equilibrium with experimental EC-derived endo form S. The highest calculated transition state for the conversion of the most stable exo M interligand to the endo S form is approximately a 28 kcal/mol above the energy of the former. The two-step interconversion of the exo H conformer to the endo S requires at least 28 kcal/mol. Surprisingly, we found that the transition state energy of the H form to the virtual endo W has the acceptable value of about 9 kcal/mol and the next energy barrier for free interconversion of endo W to endo S is 13 kcal/mol. Figure DFT Model of Epothilone Interconversions Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evaluating High-Throughput Ab Initio Gene Finders to Discover Proteins Encoded in Eukaryotic Pathogen Genomes Missed by Laboratory Techniques

Next generation sequencing technology is advancing genome sequencing at an unprecedented level. By unravelling the code within a pathogen’s genome, every possible protein (prior to post-translational modifications) can theoretically be discovered, irrespective of life cycle stages and environmental stimuli. Now more than ever there is a great need for high-throughput ab initio gene finding. Ab initio gene finders use statistical models to predict genes and their exon-intron structures from the genome sequence alone. This paper evaluates whether existing ab initio gene finders can effectively predict genes to deduce proteins that have presently missed capture by laboratory techniques. An aim here is to identify possible patterns of prediction inaccuracies for gene finders as a whole irrespective of the target pathogen. All currently available ab initio gene finders are considered in the evaluation but only four fulfil high-throughput capability: AUGUSTUS, GeneMark_hmm, GlimmerHMM, and SNAP. These gene finders require training data specific to a target pathogen and consequently the evaluation results are inextricably linked to the availability and quality of the data. The pathogen, Toxoplasma gondii, is used to illustrate the evaluation methods. The results support current opinion that predicted exons by ab initio gene finders are inaccurate in the absence of experimental evidence. However, the results reveal some patterns of inaccuracy that are common to all gene finders and these inaccuracies may provide a focus area for future gene finder developers.     

Ab Initio Calculations of Vibrational Frequencies,Potential Functions of Internal Rotations and Vibrational Infrared and Raman Spectra for 3,3,3-Trifluoropropanal

The conformational behavior and structure of 3,3,3,-trifluoropropal have been investigated by utilizing ab initio calculations with the 6-31G^** basis set (valence double zeta basis with polarization functions on all atoms) at the restricted Hartree Fock (RHF), second-order Møller-Plesset perturbation (MP2), and Density Functional (B3LYP) levels. The molecule is predicted to have a cis Û gauche conformational equilibrium. Full optimization of the transition states was performed and the rotational barriers of both the CHO and CF₃ rotors were calculated. Vibrational frequencies were computed at the three levels and the zero-point corrections were included into the calculated asymmetric CHO rotational barrier. Complete vibrational assignments were made on the basis of normal coordinate calculations for both stable conformers of the molecule.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.     

Ab initio analysis of the Cope rearrangement of germacrane sesquiterpenoids

The energetics of the Cope rearrangement of 17 germacrane sesquiterpenoids to their respective elemane forms have been calculated using both density functional theory (B3LYP/6-31G*) and post Hartee-Fock (MP2/6-31G**) ab initio methods. The calculations are in qualitative agreement with experimentally observed Cope rearrangements, but the two methods give slightly different results. MP2 calculations generally show more favorable elemene energies compared to the respective germacrenes (by around 3–4 kcal mol⁻¹) and smaller activation energies (by 2–3 kcal mol⁻¹). Additionally, neither method is accurate enough to consistently reproduce the germacrene/elemene equilibrium. Apparently, the generally small energy differences between the two forms in these sesquiterpenoids cannot be adequately reproduced at these levels of calculation. Figure The Cope rearrangement of the germacrane sesquiterpenoid bacchascandon to the elemane shyobunone     

Ab initio Models for Six-Center Multiple Proton Exchange and Ion Pair Formation Assisted by Lewis Acids

High level ab initio and density functional calculations, extrapolated to QCISD(T)/6-311+G(3df,2p)//MP2/6-31+G^**+ZPE, reveal that cyclic ion pairs can form in the hydrogen bonded complexes of haloboric acids BH_nX_3-n–HX, X=F, Cl, with Lewis bases HX, H₂O, CH₃OH, and NH₃, even in isolation (e.g., in the gas phase). The intrinsic acidities (deprotonation energies) required for protonation of these bases with formation of gas phase ion pairs are calculated to be <295 kcal/mol for water, <301 kcal/mol for methanol, and <306 kcal/mol for ammonia; such values are common for acidic sites in zeolites. All gas phase ion pairs prefer symmetric bidentate or tridentate structures. In the other cases where hydrogen bonded complexes prevail, symmetric ion pair-like transition structures for multiple hydrogen exchange are computed.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089400060563     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Theoretical study on rate constants for the reactions of CF3CH2NH2 (TFEA) with the hydroxyl radical at 298 K and atmospheric pressure

Theoretical investigations are carried out on reaction mechanism of the reactions of CF₃CH₂NH₂ (TFEA) with the OH radical by means of ab initio and DFT methods. The electronic structure information on the potential energy surface for each reaction is obtained at MPWB1K/6-31+G(d,p) level and energetic information is further refined by calculating the energy of the species with a Gaussian-2 method, G2(MP2). The existence of transition states on the corresponding potential energy surface is ascertained by performing intrinsic reaction coordinate (IRC) calculation. Our calculation indicates that the H abstraction from –NH₂ group is the dominant reaction channel because of lower energy barrier. The rate constants of the reaction calculated using canonical transition state theory (CTST) utilizing the ab initio data. The agreement between the theoretical and experimental rate constants is good at the measured temperature. From the comparison with CH₃CH₂NH₂, it is shown that the fluorine substution decreases the reactivity of the C-H bond.     

A Computational Methodology to Screen Activities of Enzyme Variants

We present a fast computational method to efficiently screen enzyme activity. In the presented method, the effect of mutations on the barrier height of an enzyme-catalysed reaction can be computed within 24 hours on roughly 10 processors. The methodology is based on the PM6 and MOZYME methods as implemented in MOPAC2009, and is tested on the first step of the amide hydrolysis reaction catalyzed by the Candida Antarctica lipase B (CalB) enzyme. The barrier heights are estimated using adiabatic mapping and shown to give barrier heights to within 3 kcal/mol of B3LYP/6-31G(d)//RHF/3-21G results for a small model system. Relatively strict convergence criteria (0.5 kcal/(molÅ)), long NDDO cutoff distances within the MOZYME method (15 Å) and single point evaluations using conventional PM6 are needed for reliable results. The generation of mutant structures and subsequent setup of the semiempirical calculations are automated so that the effect on barrier heights can be estimated for hundreds of mutants in a matter of weeks using high performance computing.     

Ab initio modeling of glycosyl torsions and anomeric effects in a model carbohydrate: 2-ethoxy tetrahydropyran

A range of ab initio calculations were carried out on the axial and equatorial anomers of the model carbohydrate 2-ethoxy tetrahydropyran to evaluate the level of theory required to accurately evaluate the glycosyl dihedral angle and the anomeric ratio. Vacuum CCSD(T)/CBS extrapolations at the global minimum yield DeltaE = E(equatorial) - E(axial) = 1.42 kcal/mol. When corrected for solvent (by the IEFPCM model), zero-point vibrations and entropy, DeltaG(298) = 0.49 kcal/mol, in excellent agreement with the experimental value of 0.47 +/- 0.3 kcal/mol. A new additivity scheme, the layered composite method (LCM), yields DeltaE to within 0.1 kcal/mol of the CCSD(T)/CBS result at a fraction of the computer requirements. Anomeric ratios and one-dimensional torsional surfaces generated by LCM and the even more efficient MP2/cc-pVTZ level of theory are in excellent agreement, indicating that the latter is suitable for force-field parameterization of carbohydrates. Hartree-Fock and density functional theory differ from CCSD(T)/CBS for DeltaE by approximately 1 kcal/mol; they show similar deviations in torsional surfaces evaluated from LCM. A comparison of vacuum and solvent-corrected one- and two-dimensional torsional surfaces indicates the equatorial form of 2-ethoxy tetrahydropyran is more sensitive to solvent than the axial.