Isopentenyl diphosphate isomerase catalyzed reactions in D2O: product release limits the rate of this sluggish enzyme-catalyzed reaction

The E. coli isopentenyl diphosphate isomerase (IDI) catalyzed reaction of isopentenyl diphosphate (IPP) in D(2)O gives a 66% yield of dimethylallyl diphosphate labeled with deuterium at the (E)-methyl group (d-DMAPP) and a 34% yield of IPP labeled with 1 mol of deuterium at C-2 (d-IPP). This shows that the release to D(2)O of the initial product of the IDI-catalyzed reaction (d-DMAPP) is slower than its conversion to d-IPP. Product dissociation is therefore rate determining for isomerization of IPP with a rate constant k(dis) ≈ k(cat) = 0.08 s(-1). The data provide an estimated rate constant of k(as) = 6 × 10(3) M(-1) s(-1) for binding of DMAPP to E. coli IDI that is similar to rate constants determined for the binding of N-protonated 2-amino ethyl diphosphate intermediate analogs to IDI from yeast [Reardon, J. E.; Abeles, R. H. Biochemistry1986, 25, 5609-5616]. We propose that ligand binding to IDI is relatively slow because there is a significant kinetic barrier to reorganization of the initial encounter complex between enzyme, substrate, and an essential Mg(2+) to form the Michaelis complex where the metal cation bridges the protein and the substrate diphosphate group.     

Direct measurement of Criegee intermediate (CH2OO) reactions with acetone, acetaldehyde, and hexafluoroacetone

Criegee biradicals, i.e., carbonyl oxides, are critical intermediates in ozonolysis and have been implicated in autoignition chemistry and other hydrocarbon oxidation systems, but until recently the direct measurement of their gas-phase kinetics has not been feasible. Indirect determinations of Criegee intermediate kinetics often rely on the introduction of a scavenger molecule into an ozonolysis system and analysis of the effects of the scavenger on yields of products associated with Criegee intermediate reactions. Carbonyl species, in particular hexafluoroacetone (CF(3)COCF(3)), have often been used as scavengers. In this work, the reactions of the simplest Criegee intermediate, CH(2)OO (formaldehyde oxide), with three carbonyl species have been measured by laser photolysis/tunable synchrotron photoionization mass spectrometry. Diiodomethane photolysis produces CH(2)I radicals, which react with O(2) to yield CH(2)OO + I. The formaldehyde oxide is reacted with a large excess of a carbonyl reactant and both the disappearance of CH(2)OO and the formation of reaction products are monitored. The rate coefficient for CH(2)OO + hexafluoroacetone is k(1) = (3.0 ± 0.3) × 10(-11) cm(3) molecule(-1) s(-1), supporting the use of hexafluoroacetone as a Criegee-intermediate scavenger. The reactions with acetaldehyde, k(2) = (9.5 ± 0.7) × 10(-13) cm(3) molecule(-1) s(-1), and with acetone, k(3) = (2.3 ± 0.3) × 10(-13) cm(3) molecule(-1) s(-1), are substantially slower. Secondary ozonides and products of ozonide isomerization are observed from the reactions of CH(2)OO with acetone and hexafluoroacetone. Their photoionization spectra are interpreted with the aid of quantum-chemical and Franck-Condon-factor calculations. No secondary ozonide was observable in the reaction of CH(2)OO with acetaldehyde, but acetic acid was identified as a product under the conditions used (4 Torr and 293 K).     

A theoretical investigation of the relative stability of hydrated glycine and methylcarbamic acid--from water clusters to interstellar ices

We have theoretically investigated how the low-energy conformers of the neutral and the zwitterionic forms of glycine as well as methylcarbamic acid are stabilized by the presence water. The MP2/6-311++G(d,p) method was utilized to conduct calculations on glycine and methylcarbamic acid in both isolated clusters and in clusters embedded in the conductor-like polarizable continuum model (C-PCM), where the clusters explicitly contain between one and ten water molecules. The neutral forms of glycine and methylcarbamic acid were found to have similar hydration energies, whereas the neutral methylcarbamic acid was determined to be approximately 32 kJ mol(-1) more stable than the neutral glycine in the isolated clusters and 30 kJ mol(-1) more stable in the C-PCM embedded clusters. Both the number and strength of the hydrogen bonding interactions between water and the zwitterions drive the stability. This lowers the relative energy of the glycine zwitterion from 50 kJ mol(-1) above neutral glycine, when there are two water molecules in the clusters to 11 kJ mol(-1) below for the clusters containing ten water molecules. For the methylcarbamic acid clusters with two water molecules, the zwitterion is 51 kJ mol(-1) higher in energy than the neutral form, but it remains 13 kJ mol(-1) above the neutral methylcarbamic acid in the clusters containing ten water molecules. When the bulk water environment is simulated by the C-PCM calculations, we find both the methylcarbamic acid and glycine zwitterionic forms have similar energies at 20 kJ mol(-1) above the neutral methylcarbamic acid energy and 10 kJ mol(-1) lower than the neutral glycine energy. Although neither methylcarbamic acid nor glycine have been detected in the interstellar medium yet, our findings indicate that methylcarbamic acid is the more stable product from methylamine and carbon dioxide reactions in a water ice. This suggests that methylcarbamic acid likely plays a role in the intermediate steps if glycine is formed in the interstellar medium.     

Hybrid modeling and simulation of stochastic effects on progression through the eukaryotic cell cycle

The eukaryotic cell cycle is regulated by a complicated chemical reaction network. Although many deterministic models have been proposed, stochastic models are desired to capture noise in the cell resulting from low numbers of critical species. However, converting a deterministic model into one that accurately captures stochastic effects can result in a complex model that is hard to build and expensive to simulate. In this paper, we first apply a hybrid (mixed deterministic and stochastic) simulation method to such a stochastic model. With proper partitioning of reactions between deterministic and stochastic simulation methods, the hybrid method generates the same primary characteristics and the same level of noise as Gillespie's stochastic simulation algorithm, but with better efficiency. By studying the results generated by various partitionings of reactions, we developed a new strategy for hybrid stochastic modeling of the cell cycle. The new approach is not limited to using mass-action rate laws. Numerical experiments demonstrate that our approach is consistent with characteristics of noisy cell cycle progression, and yields cell cycle statistics in accord with experimental observations.     

Exploring aromatic chemical space with NEAT: novel and electronically equivalent aromatic template

In this paper, we describe a lead transformation tool, NEAT (Novel and Electronically equivalent Aromatic Template), which can help identify novel aromatic rings that are estimated to have similar electrostatic potentials, dipoles, and hydrogen bonding capabilities to a query template; hence, they may offer similar bioactivity profiles. In this work, we built a comprehensive heteroaryl database, and precalculated high-level quantum mechanical (QM) properties, including electrostatic potential charges, hydrogen bonding ability, dipole moments, chemical reactivity, and othe properties. NEAT bioisosteric similarities are based on the electrostatic potential surface calculated by Brood, using the precalculated QM ESP charges and other QM properties. Compared with existing commercial lead transformation software, (1) NEAT is the only one that covers the comprehensive heteroaryl chemical space, and (2) NEAT offers a better characterization of novel aryl cores by using high-evel QM properties that are relevant to molecular interactions. NEAT provides unique value to medicinal chemists quickly exploring the largely uncharted aromatic chemical space, and one successful example of its application is discussed herein.     

The theory of surface-enhanced Raman scattering

By considering the molecule and metal to form a conjoined system, we derive an expression for the observed Raman spectrum in surface-enhanced Raman scattering. The metal levels are considered to consist of a continuum with levels filled up to the Fermi level, and empty above, while the molecule has discrete levels filled up to the highest occupied orbital, and empty above that. It is presumed that the Fermi level of the metal lies between the highest filled and the lowest unfilled level of the molecule. The molecule levels are then coupled to the metal continuum both in the filled and unfilled levels, and using the solutions to this problem provided by Fano, we derive an expression for the transition amplitude between the ground stationary state and some excited stationary state of the molecule-metal system. It is shown that three resonances contribute to the overall enhancement; namely, the surface plasmon resonance, the molecular resonances, as well as charge-transfer resonances between the molecule and metal. Furthermore, these resonances are linked by terms in the numerator, which result in SERS selection rules. These linked resonances cannot be separated, accounting for many of the observed SERS phenomena. The molecule-metal coupling is interpreted in terms of a deformation potential which is compared to the Herzberg-Teller vibronic coupling constant. We show that one term in the sum involves coupling between the surface plasmon transition dipole and the molecular transition dipole. They are coupled through the deformation potential connecting to charge-transfer states. Another term is shown to involve coupling between the charge-transfer transition and the molecular transition dipoles. These are coupled by the deformation potential connecting to plasmon resonance states. By applying the selection rules to the cases of dimer and trimer nanoparticles we show that the SERS spectrum can vary considerably with excitation wavelength, depending on which plasmon and/or charge-transfer resonance is excited.     

Synthesis and incorporation into cyclic peptides of tolan amino acids and their hydrogenated congeners: construction of an array of A-B-loop mimetics of the Cε3 domain of human IgE

The disruption of the human immunolobulin E-high affinity receptor I (IgE-FcεRI) protein-protein interaction (PPI) is a validated strategy for the development of anti asthma therapeutics. Here, we describe the synthesis of an array of conformationally constrained cyclic peptides based on an epitope of the A-B loop within the Cε3 domain of IgE. The peptides contain various tolan (i.e., 1,2-biarylethyne) amino acids and their fully and partially hydrogenated congeners as conformational constraints. Modest antagonist activity (IC(50) ～660 μM) is displayed by the peptide containing a 2,2'-tolan, which is the one predicted by molecular modeling to best mimic the conformation of the native A-B loop epitope in IgE.     

A numerically stable restrained electrostatic potential charge fitting method

Inspired by the idea of charge decomposition in calculation of the dipole preserving and polarization consistent charges (Zhang et al., J. Comput. Chem. 2011, 32, 2127), we have proposed a numerically stable restrained electrostatic potential (ESP)‐based charge fitting method for protein. The atomic charge is composed of two parts. The dominant part is fixed to a predefined value (e.g., AMBER charge), and the residual part is to be determined by restrained fitting to residual ESP on grid points around the molecule. Nonuniform weighting factors as a function of the dominant charge are assigned to the atoms. Because the residual part is several folds to several orders smaller than the dominant part, the impact of ill‐conditioning is alleviated. This charge fitting method can be used in quantum mechanical/molecular mechanical (QM/MM) simulations and similar studies, where QM calculated electronic properties are frequently mapped to partial atomic charges. © 2012 Wiley Periodicals, Inc.