Variational calculations of rovibrational states: a precise high-energy potential surface for HCN [and discussion]

We report the results of variational calculations of the rovibrational energy levels of HCN for J = 0, 1 and 2, where we reproduce all the ca. 100 observed vibrational states for all observed isotopic species, with energies up to 18000 cm$^{-1}$, to about $\pm $1 cm$^{-1}$, and the corresponding rotational constants to about $\pm $0.001 cm$^{-1}$. We use a hamiltonian expressed in internal coordinates r$_{1}$, r$_{2}$ and $\theta $, using the exact expression for the kinetic energy operator T obtained by direct transformation from the cartesian representation. The potential energy V is expressed as a polynomial expansion in the Morse coordinates y$_{i}$ for the bond stretches and the interbond angle $\theta $. The basis functions are built as products of appropriately scaled Morse functions in the bond-stretches and Legendre or associated Legendre polynomials of cos $\theta $ in the angle bend, and we evaluate matrix elements by Gauss quadrature. The hamiltonian matripx is factorized using the full rovibrational symmetry, and the basis is contracted to an optimized form; the dimensions of the final hamiltonian matrix vary from 240 $\times $ 240 to 1000 $\times $ 1000.We believe that our calculation is converged to better than 1 cm$^{-1}$ at 18 000 cm$^{-1}$. Our potential surface is expressed in terms of 31 parameters, about half of which have been refined by least squares to optimize the fit to the experimental data. The advantages and disadvantages and the future potential of calculations of this type are discussed.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Kinetic and crystallographic studies of glucopyranose spirohydantoin and glucopyranosylamine analogs inhibitors of glycogen phosphorylase

Glycogen phosphorylase (GP) is currently exploited as a target for inhibition of hepatic glycogenolysis under high glucose conditions. Spirohydantoin of glucopyranose and N-acetyl-beta-D-glucopyranosylamine have been identified as the most potent inhibitors of GP that bind at the catalytic site. Four spirohydantoin and three beta-D-glucopyranosylamine analogs have been designed, synthesized and tested for inhibition of GP in kinetic experiments. Depending on the functional group introduced, the K(i) values varied from 16.5 microM to 1200 microM. In order to rationalize the kinetic results, we determined the crystal structures of the analogs in complex with GP. All the inhibitors bound at the catalytic site of the enzyme, by making direct and water-mediated hydrogen bonds with the protein and by inducing minor movements of the side chains of Asp283 and Asn284, of the 280s loop that blocks access of the substrate glycogen to the catalytic site, and changes in the water structure in the vicinity of the site. The differences observed in the Ki values of the analogs can be interpreted in terms of variations in hydrogen bonding and van der Waals interactions, desolvation effects, ligand conformational entropy, and displacement of water molecules on ligand binding to the catalytic site.

A comparison of the reactivity of germylene and dimethylgermylene with some methylgermanes. Direct kinetic and quantum chemical studies

Time-resolved studies of germylene, GeH2, and dimethygermylene, GeMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to try to obtain rate coefficients for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas-phase. GeH2 + Me2GeH2 was studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 296-553 K. Only slight pressure dependences were found (at 386, 447 and 553 K). RRKM modelling was carried out to fit these pressure dependences. The high pressure rate coefficients gave the Arrhenius parameters: log(A/cm(3) molecule(-1)s(-1)) = -10.99 +/- 0.07 and E-a = -(7.35 +/- 0.48) kJ mol(-1). No reaction could be found between GeMe2 + Me2GeH2 at any temperature up to 549 K, and upper limits of ca. 10(-14) cm(3) molecule(-1)s(-1) were set for the rate coefficients. A rate coefficient of (1.33 +/- 0.04) x 10(-11)cm(3) molecule(-1)s(-1) was also obtained for GeH2 + MeGeH3 at 296 K. No reaction was found between GeMe2 and MeGeH3. Rate coefficient comparisons showed, inter alia, that in the substrate germane Me-for-H substitution increased the magnitudes of rate coefficients significantly, while in the germylene Me-for-H substitution decreased the magnitudes of rate coefficients by at least four orders of magnitude. Quantum chemical calculations (G2(MP2,SVP)// B3LYP level) supported these findings and showed that the lack of reactivity of GeMe2 is caused by a positive energy barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper.

Direct time-resolved study of the gas-phase reactions of germylene with ethyl- and diethylgermane: absolute rate constants, temperature dependences, and mechanism

Time-resolved studies of germylene, GeH2, generated by the 193 nm laser flash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reactions with ethyl- and diethylgermanes in the gas phase. The reactions were studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 297-564 K. Only slight pressure dependences were found for GeH2 + EtGeH3 (399, 486, and 564 K). The high pressure rate constants gave the following Arrhenius parameters: for GeH2 + EtGeH3, log A = -10.75 +/- 0.08 and E-a = -6.7 +/- 0.6 kJ mol(-1); for GeH2 + Et2GeH2, log A = -10.68 +/- 0.11 and E-a = -6.95 +/- 0.80 kJ mol(-1). These are consistent with fast, near collision-controlled, association processes at 298 K. RRKM modeling calculations are, for the most part, consistent with the observed pressure dependence of GeH2 + EtGeH3. The ethyl substituent effects have been extracted from these results and are much larger than the analogous methyl substituent effects in the SiH2 + methylsilane reaction series. This is consistent with a mechanistic model for Ge-H insertion in which the intermediate complex has a sizable secondary barrier to rearrangement.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

First detection of methylgermylene in the gas phase and time-resolved study of some of its reactions

A new transient species has been produced and detected by the gas-phase, 193 nm laser flash photolysis of 1,3,4-trimethylgermacyclopent-3-ene, TMGCP. The species has strong visible absorptions in the wavelength region 450−520 nm (maximum at 485 nm) and is attributed to the germylene, MeGeH. Time-resolved kinetic studies have led to the first rate constants for its reactions with GeH4, Me2GeH2, C2H2, C2H4, C3H6, i-C4H8, TMGCP, MeOH, HCl, and SO2. The reactivity of MeGeH is compared to those of GeH2 and GeMe2. The Me-for-H substituent effect varies according to reaction type and is not constant from GeH2 to MeGeH to GeMe2.

Gas-phase kinetic and theoretical studies of reactions of germylenes and dimethylstannylene

The results of time-resolved gas phase studies of labile germylenes (GeH2 and GeMe2) and dimethylstannylene (SnMe2) reactions reported to date are considered together with data of quantum-chemical investigations of the potential energy surfaces of these systems. Reaction mechanisms are discussed. A comparison of reactivity in the series of carbene analogs, ER2 (E = Si, Ge, Sn, R = H, Me), is made.

Asymmetric synthesis of (1R,2S,3R)-3-methylcispentacin and (1S,2S,3R)-3-methyltranspentacin by kinetic resolution of tert-butyl (+/-)-3-methylcyclopentene-1-carboxylate

Conjugate addition of lithium dibenzylamide to tert-butyl (+/-)-3-methylcyclopentene-1-carboxylate occurs with high levels of stereocontrol, with preferential addition of lithium dibenzylamide to the face of the cyclic alpha,beta-unsaturated acceptor anti- to the 3-methyl substituent. High levels of enantiorecognition are observed between tert-butyl (+/-)-3-methylcyclopentene-1-carboxylate and an excess of lithium (+/-)-N-benzyl-N-alpha-methylbenzylamide (10 eq.) (E > 140) in their mutual kinetic resolution, while the kinetic resolution of tert-butyl (+/-)-3-methylcyclopentene-1-carboxylate with lithium (S)-N-benzyl-N-alpha-methylbenzylamide proceeds to give, at 51% conversion, tert-butyl (1R, 2S, 3R,alphaS)-3-methyl-2-N-benzyl-N-alpha-methylbenzylaminocyclopentane-1-c arboxylate consistent with E > 130, and in 39% yield and 99 +/- 0.5% de after purification. Subsequent deprotection by hydrogenolysis and ester hydrolysis gives (1R, 2S, 3R)-3-methylcispentacin in > 98% de and 98 +/- 1% ee. Selective epimerisation of tert-butyl (1R, 2S, 3R, alphaS)-3-methyl-2-N- benzyl-N-alpha-methylbenzylaminocyclopentane-1-carboxylate by treatment with (KOBu)-Bu-t in (BuOH)-Bu-t gives tert-butyl (1S, 2S, 3R, alphaS)-3-methyl-2-N-benzyl-N-alpha-methylbenzylaminocyclopentane-1-carb oxylate in quantitative yield and in > 98% de, with subsequent deprotection by hydrogenolysis and ester hydrolysis giving (1S, 2S, 3R)-3-methyltranspentacin hydrochloride in > 98% de and 97 +/- 1% ee.

Metal complexes of a dipyridine octaazamacrocycle: stability constants, structural and modelling studies

Two 28-membered octaazamacrocycles, [28]py(2)N(6) and Me-2[28]py(2)N(6), have been synthesized. The protonation constants of the N-methyl. derivative and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 25degreesC in 0.10 mol dm(-3) KNO3. The high overall basicity of Me-2[28]py(2)N(6) is ascribed to the weaker repulsion between protonated contiguous charged ammonium sites separated by propyl chains. These studies together with NMR, UV-vis and EPR spectroscopies indicated the presence of mono- and di-nuclear species, The single crystal structure of the complex [Ni-2([28]py(2)N(6))(H2O)(4)]Cl-4.3H(2)O was determined, and showed each nickel centre in a distorted octahedral co-ordination environment. The nickel centres are held within the macrocycle at a large distance of 6.991(g) Angstrom from each other. The formation of mononuclear complexes was evaluated theoretically via molecular mechanics (MM) and molecular dynamics (MD) calculations and showed that these large macrocycles have sufficient flexibility to encapsulate metal ions with different stereo-electronic sizes. Structures for small and large metal ions are proposed.

Carboxylate anions binding and sensing by a novel tetraazamacrocycle containing ferrocene as receptor

A tetraazamacrocycle containing ferrocene moieties has been synthesized and characterized. The tetraprotonated form of this compound was evaluated as a receptor (R) for anion recognition of several substrates (S), Cl-, PF6-, HSO4-, H2PO4- and carboxylates, such as p-nitrobenzoate (p-nbz(-)), phthalate (ph(2-)), isophthalate (iph(2-)) and dipicolinate (dipic(2-)). H-1 NMR titrations in CD3OD indicated that this receptor is not suitable for recognizing HSO4- and H2PO4-, but weakly binds p-nbz(-), and strongly interacts with ph(2-), dipic(2-), and iph(2-) anions forming 1 : 2 assembled species. The largest beta(2) binding constant was determined for ph(2-), followed by dipic(2-) and finally iph(2-). The effect of the anionic substrates on the electron-transfer process of the ferrocene units of R was evaluated using cyclic voltammetry (CV) and square wave voltammetry (SWV) in methanol solution and 0.1 mol dm(-3) (CH3)(4)NCl as the supporting electrolyte. Titrations of the receptor were undertaken by addition of anion solutions in their tetrabutylammonium or tetramethylammonium forms. The protonated ligand exhibits a reversible voltammogram, which shifts cathodically in the presence of the substrates. The data revealed kinetic constraints in the formation of the receptor/substrate entity for dipic(2-), ph(2-) and iph(2-) anions, but not for p-nbz(-). In spite of the slow kinetics of assembled species formation with the ph(2-) substrate, this anion provides the largest redox-response when the supramolecular entity is formed, followed by dipic(2-), iph(2-) and finally p-nbz(-) anions. This trend is in agreement with the H-1 NMR results and the values of the binding constants. Single crystal X-ray structures of the receptor with PF6-, ph(2-), iph(2-) and p-nbz(-) were carried out and showed that supermolecules with a RS2 stoichiometry are formed with the first three anions, but RS4 with p-nbz(-). In all cases the binding occurs outside the macrocyclic cavity via N-H center dot center dot center dot O=C hydrogen bonds for carboxylate anions and N - H center dot center dot center dot F hydrogen bonds for the PF6- anion, which is in agreement with the solution results. The macrocyclic framework adopts different conformations in order to interact with each substrate having Fe center dot center dot center dot Fe intramolecular distances ranging from 10.125(14) to 12.783(15) angstrom.

Metal complexes of cyclen and cyclam derivatives useful for medical applications: a discussion based on thermodynamic stability constants and structural data

A series of the most common chelators used in magnetic resonance imaging ( MRI) and in radiopharmaceuticals for medical diagnosis and tumour therapy, H(4)dota, H(4)teta, H(8)dotp and H(8)tetp, is examined from a chemical point of view. Differences between 12- and 14-membered tetraazamacrocyclic derivatives with methylcarboxylate and methylphosphonate pendant arms and their chelates with divalent first-series transition metal and trivalent lanthanide ions are discussed on the basis of their thermodynamic stability constants, X- ray structures and theoretical studies.

Oxidation of biogenic and water-soluble compounds in aqueous and organic aerosol droplets by ozone: a kinetic and product analysis approach using laser Raman tweezers

The results of an experimental study into the oxidative degradation of proxies for atmospheric aerosol are presented. We demonstrate that the laser Raman tweezers method can be used successfully to obtain uptake coeffcients for gaseous oxidants on individual aqueous and organic droplets, whilst the size and composition of the droplets is simultaneously followed. A laser tweezers system was used to trap individual droplets containing an unsaturated organic compound in either an aqueous or organic ( alkane) solvent. The droplet was exposed to gas- phase ozone and the reaction kinetics and products followed using Raman spectroscopy. The reactions of three different organic compounds with ozone were studied: fumarate anions, benzoate anions and alpha pinene. The fumarate and benzoate anions in aqueous solution were used to represent components of humic- like substances, HULIS; a alpha- pinene in an alkane solvent was studied as a proxy for biogenic aerosol. The kinetic analysis shows that for these systems the diffusive transport and mass accommodation of ozone is relatively fast, and that liquid- phase di. ffusion and reaction are the rate determining steps. Uptake coe. ffcients, g, were found to be ( 1.1 +/- 0.7) x 10(-5), ( 1.5 +/- 0.7) x 10 (-5) and ( 3.0 - 7.5) x 10 (-3) for the reactions of ozone with the fumarate, benzoate and a- pinene containing droplets, respectively. Liquid- phase bimolecular rate coe. cients for reactions of dissolved ozone molecules with fumarate, benzoate and a- pinene were also obtained: k(fumarate) = ( 2.7 +/- 2) x 10 (5), k(benzoate) = ( 3.5 +/- 3) x 10 (5) and k(alpha-pinene) = ( 1-3) x 10(7) dm(3) mol (-1) s (- 1). The droplet size was found to remain stable over the course of the oxidation process for the HULIS- proxies and for the oxidation of a- pinene in pentadecane. The study of the alpha- pinene/ ozone system is the first using organic seed particles to show that the hygroscopicity of the particle does not increase dramatically over the course of the oxidation. No products were detected by Raman spectroscopy for the reaction of benzoate ions with ozone. One product peak, consistent with aqueous carbonate anions, was observed when following the oxidation of fumarate ions by ozone. Product peaks observed in the reaction of ozone with alpha- pinene suggest the formation of new species containing carbonyl groups.

Synthesis, structure and properties of cobalt(II/III) seleno-bisphenolate complexes containing NN bidentate co-ligands: kinetic studies on the reduction of a cobalt(III) complex by ascorbic acid in 80% MeCN-20% water

Two cobalt complexes, [Co(L-Se)(phen)]center dot CH2Cl2 (1) and [Co(L-Se)(N,N-Me(2)en)(CH3COO-)] (2) have been synthesized and characterized by elemental analyses, magnetic measurements, i.r. studies etc. Single crystal X- ray studies reveal that in complex (1) cobalt atom is in +2 oxidation state with trigonal bipyramidal geometry, while in complex (2) it is in +3 oxidation state and surrounded octahedrally. The asymmetric unit of complex (2) contains two crystallographically independent discrete molecules. Complex (1) was found to be paramagnetic with mu(eff) = 2.19 BM indicating a low spin cobalt(II) d(7) system, whereas complex (2) is found to be diamagnetic with cobalt(III) in low spin d(6) state. The kinetic studies on the reduction of (2) by ascorbic acid in 80% MeCN-20% H2O (v/v) at 25 degrees C reveal that the reaction proceeds through the rapid formation of inner-sphere adduct, probably by replacing the loosely coordinated AcO- group, followed by electron transfer in a slow step and is supported by a large Q (formation constant) value.