Computational Investigations of the Chalcogen-Substituted Carboxylic Acids RC(=O)XH and RC(=X)OH and their Dimers [RC(=O)XH]2 and [RC(=X)OH]2 (X=S, Se, Te)

Semiempirical and ab initio theoretical methods have been used to investigate molecular structures of the chalcogen-substituted carboxylic acid isomers RC(=O)XH (chalcogenol acid) and RC(=X)OH (chalcogenon acid). A recent experimental report suggests that the chalcogenon isomers, although less stable at room temperature, predominate at low temperature in polar solvents and that there is only a small barrier to isomerization between the isomers. Theoretical calculations have been used to locate minimum energy structures of chalcogen-substituted carboxylic acid isomers and to calculate energy differences between pairs of isomers. Carboxylic acids are well known to dimerize, especially in the gas phase and in non-polar solvents. We have, therefore, also calculated energies of dimerization of the chalcogen-substituted acids by optimizing the geometries of the symmetric dimers. We note that the PM3 level of theory is only qualitatively correct for sulfur- and selenium-containing species but fails even qualitatively for the tellurium-containing compounds. Ab initio results confirm the experimental observations and provide good estimates of both isomerization and dimerization energies. We conclude that for many functional groups with tautomers RC(=X)YH and RC(=Y)XH, the more acidic tautomer is the one with the acid proton on the smaller, more electronegative atom, although in many cases this may not be the more stable tautomer.Electronic Supplementary Material available.     

Kinetic and molecular orbital studies on the rate of oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II in comparison with the case of horseradish peroxidase

The second-order rate constant (k4) for the oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II was examined by Chance's method [B. Chance, Arch. Biochem. Biophys. 71 (1957), 130–136]. When the electronic states of these substrates were calculated by an ab initio molecular orbital method, it was found that the log k4 value correlates well with the highest occupied molecular orbital (HOMO) energy level but not with the net charge or frontier electron density. These results are essentially similar to those reported previously in the case of horseradish peroxidase [J. Sakurada, R. Sekiguchi, K. Sato, and T. Hosoya, Biochemistry 29 (1990), 4093–4098], showing some dissimilar features which are considered to reflect the structural difference between the two enzymes.Abbreviations HOMO highest occupied molecular orbital - HRP horseradish peroxidase - LPO lactoperoxidase (EC 1.11.1.7) - LUMO lowest unoccupied molecular orbital     

Molecular Parameterisation and the Theoretical Calculation of Electrode Potentials

The difference in reduction potentials between ortho and para-benzoquinones has been calculated. The employs gas phase ab initio and semi-empirical computations in combination with free energy perturbation theory applied to gas and solution phase Monte Carlo simulations. The effects on calculated results of altering solute electrostatic parameterisation in solution phase simulations is examined. Atom centred charges derived from the molecular electrostatic potentials, MEPs, from optimised ab initio wavefunctions and charges generated by consideration of hydrogen bonded complexes are considered. Parameterisation of hydroxyl torsions in hydroquinone molecules is treated in a physically realistic manner. The coupled torsional system of the ortho-hydrobenzoquinone molecule is described by a potential energy surface calculated using gas phase AM1 semi-empirical computations rather than the simple torsional energy functions frequently employed in such calculations. Calculated differences in electrode potentials show that the electrostatic interactions of quinone and hydroquinone molecules in aqueous solution are not well described by atom centred charges derived from ab initio calculated MEPs. Moreover, results in good agreement with the experimental reduction potential difference can be obtained by employing high level ab initio calculations and solution phase electrostatic parameters developed by consideration of hydrogen bonded complexes.     

Chemical shifts and three-dimensional protein structures

Summary During the past three years it has become possible to compute ab initio the ¹³C, ¹⁵N and ¹⁹F NMR chemical shifts of many sites in native proteins. Chemical shifts are beginning to become a useful supplement to more established methods of solution structure determination, and may find utility in solid-state analysis as well. From ¹³C NMR, information on , and torsions can be obtained, permitting both assignment verification, and structure refinement and prediction. For ¹⁵N, both torsional and hydrogen-bonding effects are important, while for ¹⁹F, chemical shifts are primarily indicators of the local charge field. Chemical shift calculations are still slow, but shielding hypersurfaces — the shift as a function of the dihedral angles that define the molecular conformation — are becoming accessible. Over the next few years, theoretical and computer hardware improvements will enable more routine use of chemical shifts in structural studies, including the study of metal-ligand interactions, the analysis of drug and substrate binding and catalysis, the study of folding/unfolding pathways, as well as the characterization of conformational substates. Rather than simply being a necessary prerequisite for multidimensional NMR, chemical shifts and chemical shift non-equivalence due to folding are now beginning to be useful for structural characterization.     

蝉的变音调复合声和发声机制的分析

蛙鸣蝉(Meimuna opalifera (Walk).ab.punctata Kato)的自然鸣声为“ji…guái”的重复单变调复合声.“ji”为主音频约4800Hz的准单音;“guái”的波形和主音频呈明显的演变,优势主频约2100Hz和2800Hz的变音调声.鸟鸣蝉(Meimuna opalifera (Walk)var.formosana Kato)的自然鸣声由重复的“jiū…ruǎ”和“jiū…gū…”合成的双变调复合声.“jiū”为基频和主频分别约625Hz和2100—2300Hz的准单音;“ruǎ”的波形和主音频呈明显的演变,基音和优势主频分别约575—625Hz和1550—1750Hz的变音调声.“gū”为优势主频约625Hz的准单音.变音调复合声不仅与腹部运动有关,而主要取决于发声肌的收缩特性和发声膜肋结构的振动特性.     

Extension of a de novo TIM barrel with a rationally designed secondary structure element

The ability to construct novel enzymes is a major aim in de novo protein design. A popular enzyme fold for design attempts is the TIM barrel. This fold is a common topology for enzymes and can harbor many diverse reactions. The recent de novo design of a four‐fold symmetric TIM barrel provides a well understood minimal scaffold for potential enzyme designs. Here we explore opportunities to extend and diversify this scaffold by adding a short de novo helix on top of the barrel. Due to the size of the protein, we developed a design pipeline based on computational ab initio folding that solves a less complex sub‐problem focused around the helix and its vicinity and adapt it to the entire protein. We provide biochemical characterization and a high‐resolution X‐ray structure for one variant and compare it to our design model. The successful extension of this robust TIM‐barrel scaffold opens opportunities to diversify it towards more pocket like arrangements and as such can be considered a building block for future design of binding or catalytic sites.     

Conformational behaviours of 2-substituted cyclohexanones: a complete basis set,hybrid-DFT study and NBO interpretation

The generalised anomeric effect (GAE) and gauche effect (GE) associated with donor–acceptor delocalisations, dipole–dipole interactions and total steric exchange energies (TSEE) on the conformational properties of 2-methoxy- (1), 2-methylthio- (2), 2-methylseleno- (3), 2-fluoro- (4), 2-chloro- (5) and 2-bromocyclohexanone (6) have been studied by means of ab initio and hybrid density functional theory methods and natural bond orbital (NBO) analysis. All methods used showed that the axial conformation stability increased from 2-methoxy- (1) to 2-methylselenocyclohexanone (3) and also from 2-fluoro- (4) to 2-bromocyclohexanone (6), which is in agreement with reported NMR data. The results obtained by complete basis set 4 (CBS-4), B3LYP/6-311+G** and HF/6-311+G** levels for compounds 1, 5 and 6 are very similar, but the CBS-4 results for compound 4 are not in agreement with the reported experimental data (vapour phase). The NBO analysis showed that the GAE increases from compounds 1 to 3 and also from compounds 4 to 6. The low axial conformer populations of compounds 1 and 4 can be reasonably explained by their small GAE. GE does not have significant impact on the conformational behaviours of compounds 1–6 and GAE succeeds in accounting qualitatively for the increase in the axial preferences in both series of compounds. The results showed that the calculated Δ(TSEE_eq–ax) values decrease from compounds 4 to 6 which contradicts the suggested arguments in the literature about these compounds. On the other hand, the calculated differences between the dipole moment values of the axial and equatorial conformations, Δ(μ_eq ? μ_ax), increase from compounds 1 to 2, but decrease from compounds 2 to 3 and also decrease from compounds 4 to 6. The calculated GAE values are more significant for the explanation of the conformational preferences of compounds 1–6 than the dipole–dipole repulsion effects. The correlations between the GAE, GE, dipole–dipole interactions, Wiberg Bond Index, TSEE, donor and acceptor orbital energies and occupancies, structural parameters and conformational behaviour of compounds 1–6 have been investigated.     

Comparative behavior of Pd and Ag deposition phenomenon on clean α-Al2O3 (001) surface—A first principle study

The nature of bonding at the interface between deposited silver/palladium and clean Al-terminated (001) surface of α-Al₂O₃ has been investigated using a periodic ab initio method. Substantial inter-planar relaxations within the alumina were found at both the interfaces and the bulk. The periodic calculation with both Ag and Pd deposition shows that 10% of loading on alumina results maximum stability. Surface energy and work function calculations were performed to propose the stability for the metals on the studied surfaces. The deposited Ag forms a three-dimensional (3-D) cluster on top of the alumina surface. The Pd cluster formed on the alumina surface is two-dimensional (2-D) and is distorted to accommodate the Ag cluster in its domain. A further low index calculation can explain the reason for a higher stability of the membrane generated over alumina support with silver and palladium. The results are discussed in view of the existing experimental data and models of metal-oxide interface and a reason for the difference of activity of the metal interaction with alumina surface is postulated.     

Ab initio fragment molecular orbital calculations on the specific interactions between human,mouse and rat PPARα and GW409544

To elucidate the specific interactions between the peroxisome proliferator-activated receptor (PPARα) and ligand GW409544 (GW), we obtained the solvated structures of the PPARα+GW complexes for human, mouse and rat by classical molecular mechanics calculations, and investigated their electronic properties by ab initio fragment molecular orbital calculations. The results indicate that the positively charged amino acids (Lys and Arg) of PPARα make a major contribution to the binding between PPARα and GW. In addition, it was clarified that Ser280 and Tyr314 of human and rat PPARα have a large attractive interaction with GW, while Ser280, Tyr314 and His440 of mouse PPARα have large interaction. These results on the difference in specific interactions between human and mouse/rat PPARα will be useful for predicting the effects of new chemicals on the human body based on the biomedical studies for the experimental animals such as mouse and rat.     

FT-IR,FT-Raman spectra and DFT vibrational analysis of 2-aminobiphenyl

In this work, the Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectra of 2-aminobiphenyl (2ABP) were recorded in the solid phase. The optimised geometry, frequency and intensity of the vibrational bands of 2ABP were obtained by the density functional theory (BLYP and B3LYP) methods with complete relaxation in the potential energy surface using 6-31G(d) basis set. The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically constructed spectrograms.