The Plasmodium falciparum Malaria M1 Alanyl Aminopeptidase (PfA-M1): Insights of Catalytic Mechanism and Function from MD Simulations

Malaria caused by several species of Plasmodium is major parasitic disease of humans, causing 1-3 million deaths worldwide annually. The widespread resistance of the human parasite to current drug therapies is of major concern making the identification of new drug targets urgent. While the parasite grows and multiplies inside the host erythrocyte it degrades the host cell hemoglobin and utilizes the released amino acids to synthesize its own proteins. The P. falciparum malarial M1 alanyl-aminopeptidase (PfA-M1) is an enzyme involved in the terminal stages of hemoglobin digestion and the generation of an amino acid pool within the parasite. The enzyme has been validated as a potential drug target since inhibitors of the enzyme block parasite growth in vitro and in vivo. In order to gain further understanding of this enzyme, molecular dynamics simulations using data from a recent crystal structure of PfA-M1 were performed. The results elucidate the pentahedral coordination of the catalytic Zn in these metallo-proteases and provide new insights into the roles of this cation and important active site residues in ligand binding and in the hydrolysis of the peptide bond. Based on the data, we propose a two-step catalytic mechanism, in which the conformation of the active site is altered between the Michaelis complex and the transition state. In addition, the simulations identify global changes in the protein in which conformational transitions in the catalytic domain are transmitted at the opening of the N-terminal 8 angstrom-long channel and at the opening of the 30 angstrom-long C-terminal internal chamber that facilitates entry of peptides to the active site and exit of released amino acids. The possible implications of these global changes with regard to enzyme function are discussed.

Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

Probing myo-inositol 1-phosphate synthase with multisubstrate adducts

The synthesis of a series of carbohydrate-nucleotide hybrids, designed to be multisubstrate adducts mimicking myo-inositol 1-phosphate synthase first oxidative transition state, is reported. Their ability to inhibit the synthase has been assessed and results have been rationalised computationally to estimate their likely binding mode.

GHMP Kinases - Structures, Mechanisms and Potential for Therapeutically Relevant Inhibition

The GHMP kinases are a structurally related family of small molecule kinases named after four of its members - galactokinase, homoserine kinase, mevalonate kinase and phosphomevalonate kinase. The group also includes the enzymes N-acetylgalactosamine kinase, arabinose kinase, mevalonate 5-diphosphate decarboxylase, archeal shikimate kinase and 4-(cytidine 5'-diphospho)-2-c-methyl-D-erythritol kinase. In addition the group includes two members not known to be catalytically active, the Caenorhabditis elegans sex-fate determining protein XOL-1 and the Saccharomyces cerevisiae transcriptional activator Gal3p. Two catalytic mechanisms have been proposed for GHMP kinases. The structure of mevalonate kinase suggests that an aspartate residue acts as an active site base, removing a proton from the substrate to facilitate attack on the ? phosphate of MgATP. In contrast, in homoserine kinase there is no potential catalytic base and it is proposed that catalysis is driven by transition state stabilisation. Potential chemotherapeutic interventions against GHMP kinases fall into three main categories: inhibition of galactokinase to assist suffers of galactosemia, inhibition of mevalonate kinase or mevalonate 5-diphosphate decarboxylase to reduce flux through the cholesterol biosynthesis pathway and inhibition of bacterial GHMP kinases for novel anti-microbial therapies. These are in the early stages of development, but the accumulation of structural and mechanistic data will assist future progress.

Linear Optics Simulation of Quantum Non-Markovian Dynamics

The simulation of open quantum dynamics has recently allowed the direct investigation of the features of system-environment interaction and of their consequences on the evolution of a quantum system. Such interaction threatens the quantum properties of the system, spoiling them and causing the phenomenon of decoherence. Sometimes however a coherent exchange of information takes place between system and environment, memory effects arise and the dynamics of the system becomes non-Markovian. Here we report the experimental realisation of a non-Markovian process where system and environment are coupled through a simulated transverse Ising model. By engineering the evolution in a photonic quantum simulator, we demonstrate the role played by system-environment correlations in the emergence of memory effects.

Kinetics of levoglucosan and formaldehyde formation during cellulose pyrolysis process

The mechanisms and kinetics studies of the formation of levoglucosan and formaldehyde from anhydroglucose radical have been carried out theoretically in this paper. The geometries and frequencies of all the stationary points are calculated at the B3LYP/6-31+G(D,P) level based on quantum mechanics, Six elementary reactions are found, and three global reactions are involved. The variational transition-state rate constants for the elementary reactions are calculated within 450-1500 K. The global rate constants for every pathway are evaluated from the sum of the individual elementary reaction rate constants. The first-order Arrhenius expressions for these six elementary reactions and the three pathways are suggested. By comparing with the experimental data, computational methods without tunneling correction give good description for Path1 (the formation of levoglucosan); while methods with tunneling correction (zero-curvature tunneling and small-curvature tunneling correction) give good results for Path2 (the first possibility for the formation of formaldehyde), all the test methods give similar results for Path3 (the second possibility for the formation of formaldehyde), all the modeling results for Path3 are in good agreement with the experimental data, verifying that it is the most possible way for the formation of formaldehyde during cellulose pyrolysis. © 2012 Elsevier Ltd. All rights reserved.

Kinetics study on thermal dissociation of levoglucosan during cellulose pyrolysis

(Chemical Equation Presented) The mechanisms and kinetics studies of the levoglucosan (LG) primary decomposition during cellulose pyrolysis have been carried out theoretically in this paper. Three decomposition mechanisms (C-O bond scission, C-C bond scission, and LG dehydration) including nine pathways and 16 elementary reactions were studied at the B3LYP/6-31 + G(D,P) level based on quantum mechanics. The variational transi-tion- state rate constants for every elementary reaction and every pathway were calculated within 298-1550 K. The first-order Arrhenius expressions for these 16 elementary reactions and nine pathways were suggested. It was concluded that computational method using transition state theory (TST) without tunneling correction gives good description for LG decomposition by comparing with the experimental result. With the temperature range of 667-1327 K, one dehydration pathway, with one water molecule composed of a hydrogen atom from C3 and a hydroxyl group from C2, is a preferred LG decomposition pathway by fitting well with the experimental results. The calculated Arrhenius plot of C-O bond scission mechanism is better agreed with the experimental Arrhenius plot than that of C-C bond scission. This C-O bond scission mechanism starts with breaking of C1-O5 and C6-O1 bonds with formation of CO molecule (C1-O1) simultaneously. C-C bond scission mechanism is the highest energetic barrier pathway for LG decomposition. © 2013 Elsevier Ltd. All rights reserved.

Aza-annulation of enaminones with itaconic anhydride:kinetic preference for exocyclic enamide products

Enaminones react with itaconic anhydride in methylene chloride at room temperature to give exocyclic enamides sis the major products. These can be readily equilibrated to the thermodynamically more stable endocyclic enamides. In substrates where the exocyclic isomer could not be formed only intractable materials were produced from this reaction. An intermediate in this two step process was detected and identified by proton and C-13 NMR spectroscopy. In two cases chiral enaminones were employed and the relative stereochemistry at the new chiral centre in the product was established by a crystal structure of compound 27.

Density Functional Theory Studies of Ethanol Decomposition on Rh(211)

Density functional theory calculations were carried out to examine the mechanism of ethanol decomposition on the Rh(211) surface. We found that there are two possible decomposition pathways: (1) CH(3)CH(2)OH -> CH(3)CHOH -> CH(3)COH -> CH(3)CO -> CH(3) + CO -> CH(2) + CO -> CH + CO -> C + CO and (2) CH(3)CH(2)OH -> CH(3)CHOH -> CH(3)COH -> CH(2)COH -> CHCOH -> CHCO -> CH + CO -> C + CO. Both pathways have a common intermediate of CH(3)COH, and the key step is the formation of CH(3)CHOH species. According to our calculations, the mechanism of ethanol decomposition on Rh(211) is totally different from that on Rh(111): the reaction proceeds via CH(3)COH rather than an oxametallacycle species (-CH(2)CH(2)O- for Rh( 111)), which implies that the decomposition process is structure sensitive. Further analyses on electronic structures revealed that the preference of the initial C(alpha)-H path is mainly due to the significant reduction of d-electron energy in the presence of the transition state (TS) complex, which may stabilize the TS-surface system. The present work first provides a clear picture for ethanol decomposition on stepped Rh(211), which is an important first step to completely understand the more complicated reactions, like ethanol steam reforming and electrooxidation.

Oxygen reduction reaction mechanism on nitrogen-doped graphene:A density functional theory study

Nitrogen-doped graphene (N-graphene) was reported to exhibit a good activity experimentally as an electrocatalyst of oxygen reduction reaction (ORR) on the cathode of fuel cells under the condition of electropotential of similar to 0.04 V (vs. NNE) and pH of 14. This material is promising to replace or partially replace the conventionally used Pt. In order to understand the experimental results. ORR catalyzed by N-graphene is studied using density functional theory (DFT) calculations under experimental conditions taking the solvent, surface adsorbates, and coverages into consideration. Two mechanisms, i.e., dissociative and associative mechanisms, over different N-doping configurations are investigated. The results show that N-graphene surface is covered by O with 1/6 monolayer, which is used for reactions in this work. The transition state of each elementary step was identified using four different approaches, which give rise to a similar chemistry. A full energy profile including all the reaction barriers shows that the associative mechanism is more energetically favored than the dissociative one and the removal of O species from the surface is the rate-determining step. (C) 2011 Elsevier Inc. All rights reserved.

A density functional theory study of CO oxidation on Ru(0001) at low coverage

We have performed ab initio density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on Ru(0001). Several reaction pathways and transition states are identified. A much higher reaction barrier compared to that on Pt(111) is determined, confirming that the Ru is very inactive for CO oxidation under UHV conditions. The origin of the reaction barrier was analyzed. It is found that in the transition state the chemisorbed O atom sits in an unfavorable bonding site and a significant competition for bonding with the same substrate atoms occurs between the CO and the chemisorbed O, resulting in the high barrier. Ab initio molecular dynamics calculations show that the activation of the chemisorbed O atom from the initial hcp hollow site (the most stable site) to the bridge site is the crucial step for the reaction. The CO oxidation on Ru(0001) via the Eley-Rideal mechanism has also been investigated. A comparison with previous theoretical work has been made. (C) 2000 American Institute of Physics. [S0021-9606(00)31223-5].

A first principles study of CH3 dehydrogenation on Ni(111)

Density functional theory with gradient corrections and spin polarization has been used to study the dehydrogenation of CH3 on Ni(111), a crucial step in many important catalytic reactions. The reaction, CH3(ads)--> CH2(ads)+H-(ads), is about 0.5 eV endothermic with an activation energy of more than 1 eV. The overall reaction pathway is rather intriguing. The C moiety translates from a hcp to a fcc site during the course of the reaction. The transition state of the reaction has been identified. The CH3 species is highly distorted, and both C and the active H are centered nearly on top of a row of Ni atoms with a long C-H bond length of 1.80 Angstrom. The local density of states coupled with examination of the real space distribution of individual quantum states has been used to analyze the reaction pathway. (C) 2000 American Institute of Physics. [S0021-9606(00)30218-5].

CO oxidation on Pt(111):An ab initio density functional theory study

CO oxidation on Pt(111) is studied with ab initio density functional theory. The low energy pathway and transition state for the reaction are identified. The key event is the breaking of an O-metal bond prior to the formation of a chemisorbed CO2 molecule. The pathway can be rationalized in terms of competition of the O and C atoms for bonding with the underlying surface, and the predominant energetic barrier is the strength of the O-metal bond.

Towards the enantioselective synthesis of lycoricidine alkaloids /

Two efficient, regio- and stereo controlled synthetic approaches to the synthesis of racemic analogs of pancratistatin have been accomplished and they serve as the model systems for the total synthesis of optically active 7-deoxy-pancratistatin. In the Diels-Alder approach, an efficient [4+2] cycloaddition of 3,4-methylenedioxyco- nitrostyrene with Danishefsky's diene to selectively form an exo-nitro adduct has been developed as the key step in the construction of the C-ring of the target molecule. In the Michael addition approach, the key step was a conjugate addition of an organic zinc-cuprate to the 3,4-methylenedioxy-(B-nitrostyrene, followed by a diastereocontroUed closure to form the cyclohexane C-ring of the target molecule via an intramolecular nitro-aldol cyclization on a neutral alumina surface. A chair-like transition state for such a cyclization has been established and such a chelation controlled transition state can be useful in the prediction of diastereoselectivity in other related 6-exo-trig nitroaldol reactions. Cyclization of the above products fi^om both approaches by using a Bischler-Napieralski type reaction afforded two lycoricidine derivatives 38 and 50 in good yields. The initial results from the above modeling studies as well as the analysis of the synthetic strategy were directed to a chiral pool approach to the total synthesis of optically active 7-deoxy-pancratistatin. Selective monsilylation and iodination of Ltartaric acid provided a chiral precursor for the proposed key Michael transformation. The outlook for the total synthesis of 7-deoxy-pancratistatin by this approach is very promising.A concise synthesis of novel designed, optically pure, Cz-symmetrical disulfonylamide chiral ligands starting from L-tartaric acid has also been achieved. This sequence employs the metallation of indole followed by Sfj2 replacement of a dimesylate as the key step. The activity for this Cz-symmetric chiral disulfonamide ligand in the catalytic enantioselective reaction has been confirmed by nucleophilic addition to benzaldehyde in the disulfonamide-Ti (0-i-Pr)4-diethylzinc system with a 48% yield and a 33% e.e. value. Such a ligand tethered with a suitable metal complex should be also applicable towards the total synthesis of 7-deoxy-pancratistatin.

Attemped asymmetric dehydration with chiral carbodiimides /|nMin-jen Shiao. -- 260 0 St. Catharines [Ont. : s. n.],

This research work has been planned with the intention of synthesizing optically active bicyclo[3,l,0]-hexan-2-one using chiral carbodiimides. Several carbodiimides have been prepared for practice and for attempts at asymmetric induction. The total synthesis of dibenzo[e,g]- (l:3)diazonine and the partial synthesis of l:13-dimethyldibenzo[e,g]- (l:3)diazonine are reported. Attempts to resolve 6,6f-dimethyl-2,2t-diphenic acid were not successful. The NMR spectra of carbodiimides and the related thioureas are compared. The reaction transition state of the 4-hydroxycyclohexanone with optically pure R,R(+)-di(a-phenylethyl)-carbodiimide has been considered. The ORD application to chiral cyclohexanones is discussed.