Crystal structure of Schistosoma purine nucleoside phosphorylase complexed with a novel monocyclic inhibitor

A novel inhibitor of Schistosoma PNP was identified using an ""in silico"" approach allied to enzyme inhibition assays. The compound has a monocyclic structure which has not been previously described for PNP inhibitors The crystallographic structure of the complex was determined and used to elucidate the binding mode within the active site Furthermore, the predicted pose was very similar to that determined crystallographically, validating the methodology The compound Sm_VS1, despite its low molecular weight, possesses an IC(50) of 1 3 mu M, surprisingly low when compared with purine analogues This is presumably due to the formation of eight hydrogen bonds with key residues in the active site E203, N245 and T244. The results of this study highlight the importance of the use of multiple conformations for the target during virtual screening. Indeed the Sm_VS1 compound was only identified after flipping the N245 side chain It is expected that the structure will be of use in the development of new highly active non-purine based compounds against the Sclustosoma enzyme. (c) 2010 Elsevier B V. All rights reserved

Capillary bioreactors based on human purine nucleoside phosphorylase: A new approach for ligands identification and characterization

The enzyme purine nucleoside phosphorylase (PNP) is a target for the discovery of new lead compounds employed on the treatment severe T-cell mediated disorders. Within this context, the development of new, direct, and reliable methods for ligands screening is an important task. This paper describes the preparation of fused silica capillaries human PNP (HsPNP) immobilized enzyme reactor (IMER). The activity of the obtained IMER is monitored on line in a multidimensional liquid chromatography system, by the quantification of the product formed throughout the enzymatic reaction. The Km value for the immobilized enzyme was about twofold higher than that measured for the enzyme in solution (255 +/- 29.2 mu M and 133 +/- 114.9 mu M, respectively). A new fourth-generation immucillin derivative (DI4G: IC50 = 40.6 +/- 0.36 nM), previously identified and characterized in HsPNP free enzyme assays, was used to validate the IMER as a screening method for HsPNP ligands. The validated method was also used for mechanistic studies with this inhibitor. This new approach is a valuable tool to PNP ligand screening, since it directly measures the hypoxanthine released by inosine phosphorolysis, thus furnishing more reliable results than those one used in a coupled enzymatic spectrophotometric assay. (C) 2011 Elsevier B.V. All rights reserved.

Insights into Phosphate Cooperativity and Influence of Substrate Modifications on Binding and Catalysis of Hexameric Purine Nucleoside Phosphorylases

The hexameric purine nucleoside phosphorylase from Bacillus subtilis (BsPNP233) displays great potential to produce nucleoside analogues in industry and can be exploited in the development of new anti-tumor gene therapies. In order to provide structural basis for enzyme and substrates rational optimization, aiming at those applications, the present work shows a thorough and detailed structural description of the binding mode of substrates and nucleoside analogues to the active site of the hexameric BsPNP233. Here we report the crystal structure of BsPNP233 in the apo form and in complex with 11 ligands, including clinically relevant compounds. The crystal structure of six ligands (adenine, 2'deoxyguanosine, aciclovir, ganciclovir, 8-bromoguanosine, 6-chloroguanosine) in complex with a hexameric PNP are presented for the first time. Our data showed that free bases adopt alternative conformations in the BsPNP233 active site and indicated that binding of the co-substrate (2'deoxy) ribose 1-phosphate might contribute for stabilizing the bases in a favorable orientation for catalysis. The BsPNP233-adenosine complex revealed that a hydrogen bond between the 5' hydroxyl group of adenosine and Arg(43*) side chain contributes for the ribosyl radical to adopt an unusual C3'-endo conformation. The structures with 6-chloroguanosine and 8-bromoguanosine pointed out that the Cl-6 and Br-8 substrate modifications seem to be detrimental for catalysis and can be explored in the design of inhibitors for hexameric PNPs from pathogens. Our data also corroborated the competitive inhibition mechanism of hexameric PNPs by tubercidin and suggested that the acyclic nucleoside ganciclovir is a better inhibitor for hexameric PNPs than aciclovir. Furthermore, comparative structural analyses indicated that the replacement of Ser(90) by a threonine in the B. cereus hexameric adenosine phosphorylase (Thr(91)) is responsible for the lack of negative cooperativity of phosphate binding in this enzyme.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Buckled nano rod - a two state system: quantum effects on its dynamics

We consider a suspended elastic rod under longitudinal compression. The compression can be used to adjust potential energy for transverse displacements from the harmonic to the double well regime. The two minima in potential energy curve describe two possible buckled states. Using transition state theory (TST) we have calculated the rate of conversion from one state to other. If the strain epsilon = 4 epsilon c the simple TST rate diverges. We suggest a method to correct this divergence for quantum calculations. We also find that zero point energy contributions can be quite large so that single mode calculations can lead to large errors in the rate.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

Investigations of group IVA transition metal mediated carbon-carbon bond forming reactions

Zirconocene aldehyde and ketone complexes were synthesized in high yield by treatment of zirconocene acyl complexes with trimethylaluminum or diisobutylaluminum hydride. These complexes, which are activated by dialkylaluminum chloride ligands, inserted unsaturated substrates such as alkynes, allenes, ethylene, nitriles, ketenes, aldehydes, ketones, lactones, and acid chlorides with moderate to high conversion. Insertion of aldehyde substrates yielded zirconocene diolate complexes with up to 20:1 (anti:syn) diastereoselectivity. The zirconocene diolates were hydrolyzed to afford unsymmetrical 1,2-diols in 40-80% isolated yield. Unsymmetrical ketones gave similar insertion yields with little or no diastereoselectivity. A high yielding one-pot method was developed that coupled carbonyl substrates with zirconocene aldehyde complexes that were derived from olefins by hydrozirconation and carbonylation. The zirconocene aldehyde complexes also inserted carbon monoxide and gave acyloins in 50% yield after hydrolysis.

The insertion reaction of aryl epoxides with the trimethylphoshine adduct of titanocene methylidene was examined. The resulting oxytitanacyclopentanes were carbonylated and oxidatively cleaved with dioxygen to afford y-lactones in moderate yields. Due to the instability and difficult isolation of titanocene methylidene trimethylphoshine adducts, a one-pot method involving the addition of catalytic amounts of trimethylphosphine to β,β-dimethyltitanacyclobutane was developed. A series of disubstituted aryl epoxides were examined which gave mixtures of diastereomeric insertion products. Based on these results, as well as earlier Hammett studies and labeling experiments, a biradical transition state intermediate is proposed. The method is limited to aryl substituted epoxide substrates with aliphatic examples showing no insertion reactivity.

The third study involved the use of magnesium chloride supported titanium catalysts for the Lewis acid catalyzed silyl group transfer condensation of enol silanes with aldehydes. The reaction resulted in silylated aldol products with as many as 140 catalytic turnovers before catalyst inactivation. Low diastereoselectivities favoring the anti-isomer were consistent with an open transition state involving a titanium atom bound to the catalyst surface. The catalysts were also used for the aldol group transfer polymerization of t-butyldimethylsilyloxy-1-ethene resulting in polymers with molecular weights of 5000-31,000 and molar mass dispersities of 1.5-2.8. Attempts to polymerize methylmethacrylate using GTP proved unsuccessful with these catalysts.

State-to-state dynamics of elementary chemical reactions using Rydberg H-atom translational spectroscopy

In this review, a few examples of state-to-state dynamics studies of both unimolecular and bimolecular reactions using the H-atom Rydberg tagging TOF technique were presented. From the H2O photodissociation at 157 nm, a direction dissociation example is provided, while photodissociation of H2O at 121.6 has provided an excellent dynamical case of complicated, yet direct dissociation process through conical intersections. The studies of the O(D-1) + H-2 --> OH+H reaction has also been reviewed here. A prototype example of state-to-state dynamics of pure insertion chemical reaction is provided. Effect of the reagent rotational excitation and the isotope effect on the dynamics of this reaction have also been investigated. The detailed mechanism for abstraction channel in this reaction has also been closely studied. The experimental investigations of the simplest chemical reaction, the H-3 system, have also been described here. Through extensive collaborations between theory and experiment, the mechanism for forward scattering product at high collision energies for the H+HD reaction was clarified, which is attributed to a slow down mechanism on the top of a quantized barrier transition state. Oscillations in the product quantum state resolved different cross sections have also been observed in the H+D-2 reaction, and were attributed to the interference of adiabatic transition state pathways from detailed theoretical analysis. The results reviewed here clearly show the significant advances we have made in the studies of the state-to-state molecular reaction dynamics.

The near-edge structure in energy-loss spectroscopy: many-electron and magnetic effects in transition metal nitrides and carbides

We investigate the ability of the local density approximation (LDA) in density functional theory to predict the near-edge structure in electron energy-loss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion K-edge threshold energies are systematically overestimated by 4.22 +/- 0.44 eV in twelve transition metal carbides and nitrides in the rock-salt (B1) structure. When we apply this 'universal' many-electron correction to energy-loss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the near-edge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.

Trends in C-O and C-N bond formations over transition metal surfaces: An insight into kinetic sensitivity in catalytic reactions

Transition metal catalyzed bond formation is a fundamental process in catalysis and is of general interest throughout chemistry. To date, however, the knowledge of association reactions is rather limited, relative to what is known about dissociative processes. For example, surprisingly little is known about how the bond-forming ability of a metal, in general, varies across the Periodic Table. In particular, the effect of reactant valency on such trends is poorly understood. Herein, the authors examine these key issues by using density functional theory calculations to study CO and CN formations over the 4d metals. The calculations reveal that the chemistries differ in a fundamental way. In the case of CO formation, the reaction enthalpies span a much greater range than those of CN formation. Moreover, CO formation is found to be kinetically sensitive to the metal; here the reaction barriers (E-a) are found to be influenced by the reaction enthalpy. CN formation, conversely, is found to be relatively kinetically insensitive to the metal, and there is no correlation found between the reaction barriers and the reaction enthalpy. Analysis has shown that at the final adsorbed state, the interaction between N and the surface is relatively greater than that of O. Furthermore, in comparison with O, relatively less bonding between the surface and N is observed to be lost during transition state formation. These greater interactions between N and the surface, which can be related to the larger valency of N, are found to be responsible for the relatively smaller enthalpy range and limited variation in E-a for CN formation. (C) 2007 American Institute of Physics.

Importance of electronegativity differences and surface structure in molecular dissociation reactions at transition metal surfaces

The dissociative adsorption of N-2 has been studied at both monatomic steps and flat regions on the surfaces of the 4d transition metals from Zr to Pd. Using density functional theory (DFT) calculations, we have determined and analyzed the trends in both straight reactivity and structure sensitivity across the periodic table. With regards to reactivity, we find that the trend in activation energy (Ea) is determined mainly by a charge transfer from the surface metal atoms to the N atoms during transition state formation, namely, the degree of ionicity of the N-surface bond at the transition state. Indeed, we find that the strength of the metal-N bond at the transition state (and therefore the trend in Ea) can be predicted by the difference in Mulliken electronegativity between the metal and N. Structure sensitivity is analyzed in terms of geometric and electronic effects. We find that the lowering of Ea due to steps is more pronounced on the right-hand side of the periodic table. It is found that for the early transition metals the geometric and electronic effects work in opposition when going from terrace to step active site. In the case of the late 4d metals, however, these effects work in combination, producing a more marked reduction in Ea.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

Insights into the Staggered Nature of Hydrogenation Reactivity over the 4d Transition Metals

Hydrogenation reactions at transition metal surfaces comprise a key set of reactions in heterogeneous catalysis. In this paper, density functional theory methods are employed to take an in-depth look at this fundamental reaction type. The energetics of hydrogenation of atomic C, N, and O have been studied in some detail over low index Zr, Nb, Mo, Tc, Ru, Rh, and Pd surfaces. Detailed bonding analysis has also been employed to track carefully the chemical changes taking place during reaction. A number of interesting horizontal and vertical trends have been uncovered relating to reactant valency and metal d-band filling. A general correlation has also been found between the reaction barriers and the reaction potential energies. Moreover, when each reaction is considered independently, correlation has been found to improve with decreasing reactant valency. Bonding analysis has pointed to this being related to the relative position of the transition state along the reaction coordinate and has shown that as reactant valency decreases, the transition states become progressively later.

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.

Femtosecond real-time probing of reactions. 12. Vectorial dynamics of transition states

Femtosecond time-resolved techniques with KETOF (kinetic energy time-of-flight) detection in a molecular beam are developed for studies of the vectorial dynamics of transition states. Application to the dissociation reaction of IHgI is presented. For this system, the complex [I---Hg---I](++)* is unstable and, through the symmetric and asymmetric stretch motions, yields different product fragments: [I---Hg---I](++)* -> HgI(X^2/sigma^+) + I(^2P_3/2) [or I*(^2P_l/2)] (1a); [I---Hg---I](++)* -> Hg(^1S_0) + I(^2P_3/2) + I(^2P_3/2) [or I* (^2P_1/2)] (1 b). These two channels, (1a) and (1b), lead to different kinetic energy distributions in the products. It is shown that the motion of the wave packet in the transition-state region can be observed by MPI mass detection; the transient time ranges from 120 to 300 fs depending on the available energy. With polarized pulses, the vectorial properties (transition moments alignment relative to recoil direction) are studied for fragment separations on the femtosecond time scale. The results indicate the nature of the structure (symmetry properties) and the correlation to final products. For 311-nm excitation, no evidence of crossing between the I and I* potentials is found at the internuclear separations studied. (Results for 287-nm excitation are also presented.) Molecular dynamics simulations and studies by laser-induced fluorescence support these findings.