The structure of 5a,6–anhydrotetracycline and its Mg2+ complexes in aqueous solution

Semiempirical molecular orbital theory has been used for a systematic scan of the binding positions for a Mg²⁺ ion with 5a,6–anhydrotetracycline taking both conformational flexibility and possible different tautomeric forms into account. The magnesium ion has been calculated alone and with four or five complexed water molecules in order to simulate the experimental situation more closely. The results are analyzed by comparing the behavior of the title compound with that of tetracycline itself and possible causes for the stronger induction of the Tetracycline Receptor (TetR) by 5a,6–anhydrotetracycline than by tetracycline are considered. Energetically favored 3D -structure of the zwitteranionic 5a,6-anhydrotetracycline magnesium complex in solution Electronic Supplementary Material Supplementary material is available for this article at     

An experimental and theoretical study of the enantioselective deprotonation of cyclohexene oxide with isopinocampheyl-based chiral lithium amides

The mechanism of the enantioselective deprotonation of cyclohexene oxide with isopinocampheyl-based chiral lithium amide was studied by quantum chemical calculations. The transition states of eight molecules were fully optimized at the ab initio HF/3-21G and density functional B3LYP/3-21G levels with Gaussian 98. The activation energies were calculated at the B3LYP/6-31+G(3df,2p)//B3LYP/3-21G level. We found the theoretical evaluation to be consistent with the experimental data. At the best case, an enantiomeric excess of up to 95% for (R)-2-scyclohexen-1-ol was achieved with (−)-N, N-diisopinocampheyl lithium amide. Enantioselective deprotonation of cyclohexene oxide Electronic Supplementary Material Supplementary material is available for this article at Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Evaluation of the chemical reactivity in lignin precursors using the Fukui function

The hydroxycinnamyl alcohols: p-coumarol, coniferol and sinapol are considered the basic units and precursors of lignins models. In this work, the specific reactivity of these molecules was studied. We investigate their intrinsic chemical reactivity in terms of the Fukui function, applying the principle of hard and soft acids and bases (HSAB) in the framework of the density functional theory (DFT). Comparisons of their nucleophilic, electrophilic and free radical reactivity show their most probably sites to form linkages among them. It is found that the most reactive sites, for reactions involving free radicals, are the carbons at the β-position in the p-coumarol and sinapol molecules, whilst the regions around the carbon-oxygen bond of the phenoxyl group are the most reactive in coniferol. Figure Isocontour plots for the free radical form of the Fukui function f ⁰ (r), showing the reactive sites toward electron-rich/poor reactants     

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

Molecular model for host–guest interaction of tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylthiacalix[4]arene and tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylcalix[4]arene receptors with carboxylate and dicarboxylate guests: an ONIOM study

Geometry optimizations of tetraamino-tert-butylthiacalix[4]arene (tatbtc4a) and tetraamino-tert-butylcalix[4]arene (tatbc4a) complexes with acetate, oxalate, malonate, succinate, glutarate, adipate, and pimelate were carried out using the integrated MO:MO method. Thermodynamic quantities, preorganization energies and complexation energies of these complexes were obtained at the ONIOM(B3LYP/6-31G(d):AM1) level of theory. The relative stabilities of the tatbtc4a and tatbc4a complexes with carboxylate guests are reported. The complexes tatbtc4a/malonate and tatbc4a/oxalate were found to be the most stable species. The selectivity of the tatbtc4a receptor toward to malonate with respect to oxalate, in terms of selectivity coefficient, is 9.90×10². Figure Atom labeling of tatbtc4a/oxalate complex as a representative of host-guest system.     

Interactions between the photosystem II subunit PsbS and xanthophylls studied in vivo and in vitro

The photosystem II subunit PsbS is essential for excess energy dissipation (qE); however, both lutein and zeaxanthin are needed for its full activation. Based on previous work, two models can be proposed in which PsbS is either 1) the gene product where the quenching activity is located or 2) a proton-sensing trigger that activates the quencher molecules. The first hypothesis requires xanthophyll binding to two PsbS-binding sites, each activated by the protonation of a dicyclohexylcarbodiimide-binding lumen-exposed glutamic acid residue. To assess the existence and properties of these xanthophyll-binding sites, PsbS point mutants on each of the two Glu residues PsbS E122Q and PsbS E226Q were crossed with the npq1/npq4 and lut2/npq4 mutants lacking zeaxanthin and lutein, respectively. Double mutants E122Q/npq1 and E226Q/npq1 had no qE, whereas E122Q/lut2 and E226Q/lut2 showed a strong qE reduction with respect to both lut2 and single glutamate mutants. These findings exclude a specific interaction between lutein or zeaxanthin and a dicyclohexylcarbodiimide-binding site and suggest that the dependence of nonphotochemical quenching on xanthophyll composition is not due to pigment binding to PsbS. To verify, in vitro, the capacity of xanthophylls to bind PsbS, we have produced recombinant PsbS refolded with purified pigments and shown that Raman signals, previously attributed to PsbS-zeaxanthin interactions, are in fact due to xanthophyll aggregation. We conclude that the xanthophyll dependence of qE is not due to PsbS but to other pigment-binding proteins, probably of the Lhcb type.     

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

The design of organic catalysis for epoxidation by hydrogen peroxide

The potential of various organic species to catalyze epoxidation of ethene by hydrogen peroxide is explored with B3LYP/6-31G* DFT calculations. Electronic Supplementary Materials Supplementary material is available for this article at http://dx.doi.org/10.1007/s00894-005-0044-4.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

DFT investigation of the 1-octene metathesis reaction mechanism with the Phobcat precatalyst

The productive self-metathesis reaction of 1-octene in the presence of the Phobcat precatalyst [RuCl₂(Phoban-Cy)₂(=CHPh)] using density functional theory was investigated and compared to the Grubbs 1 precatalyst [RuCl₂(PCy₃)₂(=CHPh)]. At the GGA-PW91/DNP level, the geometry optimization of all the participating species and the PES scans of the various activation and catalytic cycles in the dissociative mechanism were performed. The formation of the catalytically active heptylidene species is kinetically and thermodynamically favored, while the formation of trans-tetradecene is thermodynamically favored.      

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

CHIH-DFT determination of the molecular structure and IR and UV spectra of solanidine

Solanidine is the steroidal aglycon of some potato glycoalkaloids and a very important precursor for the synthesis of hormones and some pharmacologically active compounds. In this work, we make use of a new chemistry model within Density Functional Theory, called CHIH-DFT, to calculate the molecular structure of solanidine, as well to predict its infrared and ultraviolet spectra. The calculated values are compared with the experimental data available for this molecule as a means of validation of our proposed chemistry model. Figure Molecular structure of solanidine calculated with the CHIH(small) model chemistry     

Three-dimensional structure of the catalytic domain of the yeast β-(1,3)-glucan transferase Gas1: a molecular modeling investigation

The three-dimensional (3D) structure of the catalytic domain of Gas1p, a protein belonging to the only family of β-(1,3)-glucan transferases so far identified in yeasts and some pathogenic fungi (family GH-72), has been predicted by combining results derived from threading methods, multiple sequence alignments and secondary-structure predictions. The 3D model has allowed the identification of several residues that are predicted to play a crucial role in structural integrity, substrate recognition and catalysis. In particular, the model of the catalytic domain can be useful for designing site-directed mutagenesis experiments and for developing inhibitors of Gas1p enzymatic activity. Figure Three-dimensional models of the Gas1p catalytic domain as predicted using as template 7A3H (PDB code) protein Electronic Supplementary Material Supplementary material is available for this article at     

Periodic trends and easy estimation of relative stabilities in 11-vertex nido-p-block-heteroboranes and -borates

Density functional theory computations were carried out for 11-vertex nido-p-block-hetero(carba)boranes and -borates containing silicon, germanium, tin, arsenic, antimony, sulfur, selenium and tellurium heteroatoms. A set of quantitative values called “estimated energy penalties” was derived by comparing the energies of two reference structures that differ with respect to one structural feature only. These energy penalties behave additively, i.e., they allow us to reproduce the DFT-computed relative stabilities of 11-vertex nido-heteroboranes in general with good accuracy and to predict the thermodynamic stabilities of unknown structures easily. Energy penalties for neighboring heteroatoms (HetHet and HetHet′) decrease down the group and increase along the period (indirectly proportional to covalent radii). Energy penalties for a five- rather than four-coordinate heteroatom, [Het_5k(1) and Het_5k(2)], generally, increase down group 14 but decrease down group 16, while there are mixed trends for group 15 heteroatoms. The sum of HetHet′ energy penalties results in different but easily predictable open-face heteroatom positions in the thermodynamically most stable mixed heterocarbaboranes and -borates with more than two heteroatoms. Figure Correlation of HetHet′ and HetC increments with covalent radii of group 15 heteroatoms Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Acclimation- and mutation-induced enhancement of PsbS levels affects the kinetics of non-photochemical quenching in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

The efficiency of photosystem II antenna complexes (LHCs) in higher plants must be regulated to avoid potentially damaging overexcitation of the reaction centre in excess light. Regulation is achieved via a feedback mechanism known as non-photochemical quenching (NPQ), triggered the proton gradient (ΔpH) causing heat dissipation within the LHC antenna. ΔpH causes protonation of the LHCs, the PsbS protein and triggers the enzymatic de-epoxidation of the xanthophyll, violaxanthin, to zeaxanthin. A key step in understanding the mechanism is to decipher whether PsbS and zeaxanthin cooperate to promote NPQ. To obtain clues about their respective functions we studied the effects of PsbS and zeaxanthin on the rates of NPQ formation and relaxation in wild-type Arabidopsis leaves and those overexpressing PsbS (L17) or lacking zeaxanthin (npq1). Overexpression of PsbS was found to increase the rate of NPQ formation, as previously reported for zeaxanthin. However, PsbS overexpression also increased the rate of NPQ relaxation, unlike zeaxanthin, which is known decrease the rate. The enhancement of PsbS levels in plants lacking zeaxanthin (npq1) by either acclimation to high light or crossing with L17 plants showed that the effect of PsbS was independent of zeaxanthin. PsbS levels also affected the kinetics of the 535 nm absorption change (ΔA535), which monitors the formation of the conformational state of the LHC antenna associated with NPQ, in an identical way. The antagonistic action of PsbS and zeaxanthin with respect to NPQ and ΔA535 relaxation kinetics suggests that the two molecules have distinct regulatory functions.     

Molecular models of protein kinase 6 from Plasmodium falciparum

Cyclin-dependent kinases (CDKs) have been identified as potential targets for development of drugs, mainly against cancer. These studies generated a vast library of chemical inhibitors of CDKs, and some of these molecules can also inhibit kinases identified in the Plasmodium falciparum genome. Here we describe structural models for Protein Kinase 6 from P. falciparum (PfPK6) complexed with Roscovitine and Olomoucine. These models show clear structural evidence for differences observed in the inhibition, and may help designing inhibitors for PfPK6 generating new potential drugs against malaria. Figure Ribbon diagram of PfPK6 complexed with a roscovitine and b olomoucine     

Insights into catalytic activity of industrial enzyme Co-nitrile hydratase. Docking studies of nitriles and amides

Nitrile hydratase (NHase) is an enzyme containing non-corrin Co³⁺ in the non-standard active site. NHases from Pseudonocardia thermophila JCM 3095 catalyse hydration of nitriles to corresponding amides. The efficiency of the enzyme is 100 times higher for aliphatic nitriles then aromatic ones. In order to understand better this selectivity dockings of a series of aliphatic and aromatic nitriles and related amides into a model protein based on an X-ray structure were performed. Substantial differences in binding modes were observed, showing better conformational freedom of aliphatic compounds. Distinct interactions with postranslationally modified cysteines present in the active site of the enzyme were observed. Modeling shows that water molecule activated by a metal ion may easily directly attack the docked acrylonitrile to transform this molecule into acryloamide. Thus docking studies provide support for one of the reaction mechanisms discussed in the literature. Figure Crystalographic structure of Pseudonocardia thermophila JCM 3095 nitrile hydratase (a) and the non-standard active site (b)     

Quantitative trait loci identified for sugar related traits in a sugarcane (<Emphasis Type="Italic">Saccharum</Emphasis> spp.) cultivar × <Emphasis Type="Italic">Saccharum officinarum</Emphasis> population

The identification of markers linked to quantitative trait loci (QTLs) for increased sugar accumulation could improve the effectiveness of current breeding strategies in sugarcane. Progeny from a cross between a high sucrose producing cultivar, (denotes Australian plant breeding rights), and a Saccharum officinarum clone, IJ76-514 were grown in two field experiments in different years, and evaluated in the early and mid-season phases of crop maturity, to identify robust QTLs in affecting sucrose content in cane. Using an extensive genetic map constructed for with over 1,000 AFLP and SSR markers, a total of 37 QTLs were identified for brix and pol of which, 16 were detected in both experiments. Of these 37 QTL, 30 were clustered into 12 genomic regions in six of the eight homo(eo)logous groups. Each QTL explained from 3 to 9% of the phenotypic variation observed. Both positive and negative effects were identified and the location of the QTLs on linkage groups belonging to the same homo(eo)logy group suggested that a number of the QTLs were allelic forms of the same genes. Of the 37 QTLs identified, the majority were significant in both early and mature cane, but 8 were identified as early specific QTLs and 9 as mature cane QTLs. In total, 97 interactions were significant (P<10⁻⁵) and these were localised to 32 genomic regions of which 6 were detected with both years’ data. Models including all the QTLs explained from 37 to 66% of the total phenotypic variation, depending on the trait. The results will be subsequently applied in marker assisted breeding. denotes variety covered by Australian plant breeding rights.     

Theoretical study of [ XN5 ]− (X=O, S, Se, Te) systems

A series of [XN₅]⁻ (X=O, S, Se, Te) compounds has been examined with ab initio and Density Functional Theory (DFT) methods. The five-membered nitrogen ring series of structures are global minima and may exist or be characterized due to their significant dissociation barriers (29.7–32.7 kcal mol⁻¹). Nucleus-independent chemical shifts (NICS) criteria and the presence of (4n+2) π-electrons confirmed that the five-membered nitrogen ring in their structures exhibits characteristics of aromaticity. Thus, the strong stability of the five-membered nitrogen ring structures may be attributed partially to their aromaticity.      

Localized and delocalized perfluorosemibullvalenes

On the basis of the experimental Gibbs free-energy barrier of the degenerate Cope arrangement in semibullvalene, B3P86 shows the best agreement, while B3LYP and MP2 underestimate and CCSD(T) overestimates the barrier. The substituent effect proposal by Hoffmann has been verified. In contrast to semibullvalenes with either localized energy-minimum structures or delocalized transition-state structures, perfluorosemibullvalene has both localized and delocalized energy-minimum structures that are very close in energy. Localized and delocalized perfluorosemibullvalenes Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.     

Pattern recognition based on color-coded quantum mechanical surfaces for molecular alignment

A pattern recognition algorithm for the alignment of drug-like molecules has been implemented. The method is based on the calculation of quantum mechanical derived local properties defined on a molecular surface. This approach has been shown to be very useful in attempting to derive generalized, non-atom based representations of molecular structure. The visualization of these surfaces is described together with details of the methodology developed for their use in molecular overlay and similarity calculations. In addition, this paper also introduces an additional local property, the local curvature (C _L), which can be used together with the quantum mechanical properties to describe the local shape. The method is exemplified using some problems representing common tasks encountered in molecular similarity. Figure Molecular surfaces for Lorazepam (left) and Diazepam (right)