Adenosine binding to low-molecular-weight purine nucleoside phosphorylase: the structural basis for recognition based on its complex with the enzyme from Schistosoma mansoni

Schistosomes are unable to synthesize purines de novo and depend exclusively on the salvage pathway for their purine requirements. It has been suggested that blockage of this pathway could lead to parasite death. The enzyme purine nucleoside phosphorylase (PNP) is one of its key components and molecules designed to inhibit the low-molecular-weight (LMW) PNPs, which include both the human and schistosome enzymes, are typically analogues of the natural substrates inosine and guanosine. Here, it is shown that adenosine both binds to Schistosoma mansoni PNP and behaves as a weak micromolar inhibitor of inosine phosphorolysis. Furthermore, the first crystal structures of complexes of an LMW PNP with adenosine and adenine are reported, together with those with inosine and hypoxanthine. These are used to propose a structural explanation for the selective binding of adenosine to some LMW PNPs but not to others. The results indicate that transition-state analogues based on adenosine or other 6-amino nucleosides should not be discounted as potential starting points for alternative inhibitors.

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

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

Characterization of the interdependency between residues that bind the substrate in a β-glycosidase

The manner by which effects of simultaneous mutations combine to change enzymatic activity is not easily predictable because these effects are not always additive in a linear manner. Hence, the characterization of the effects of simultaneous mutations of amino acid residues that bind the substrate can make a significant contribution to the understanding of the substrate specificity of enzymes. In the β-glycosidase from Spodoptera frugiperda (Sfβgly), both residues Q39 and E451 interact with the substrate and this is essential for defining substrate specificity. Double mutants of Sfβgly (A451E39, S451E39 and S451N39) were prepared by site-directed mutagenesis, expressed in bacteria and purified using affinity chromatography. These enzymes were characterized using p-nitrophenyl β-galactoside and p-nitrophenyl β-fucoside as substrates. The k cat/Km ratio for single and double mutants of Sfβgly containing site-directed mutations at positions Q39 and E451 was used to demonstrate that the effect on the free energy of ES‡ (enzyme-transition state complex) of the double mutations (∆∆G‡xy) is not the sum of the effects resulting from the single mutations (∆∆G‡x and ∆∆G‡y). This difference in ∆∆G‡ indicates that the effects of the single mutations partially overlap. Hence, this common effect counts only once in ∆∆G‡xy. Crystallographic data on β-glycosidases reveal the presence of a bidentate hydrogen bond involving residues Q39 and E451 and the same hydroxyl group of the substrate. Therefore, both thermodynamic and crystallographic data suggest that residues Q39 and E451 exert a mutual influence on their respective interactions with the substrate.

Uma perspectiva computacional sobre catálise enzimática

Enzymes are extremely efficient catalysts. Here, part of the mechanisms proposed to explain this catalytic power will be compared to quantitative experimental results and computer simulations. Influence of the enzymatic environment over species along the reaction coordinate will be analysed. Concepts of transition state stabilisation and reactant destabilisation will be confronted. Divided site model and near-attack conformation hypotheses will also be discussed. Molecular interactions such as covalent catalysis, general acid-base catalysis, electrostatics, entropic effects, steric hindrance, quantum and dynamical effects will also be analysed as sources of catalysis. Reaction mechanisms, in particular that catalysed by protein tyrosine phosphatases, illustrate the concepts.

Gas-phase solvolysis type reactions of SiCl3+ cations

Gas-phase SiCl3+ ions undergo sequential solvolysis type reactions with water, methanol, ammonia, methylamine and propylene. Studies carried out in a Fourier Transform mass spectrometer reveal that these reactions are facile at 10-8 Torr and give rise to substituted chlorosilyl cations. Ab initio and DFT calculations reveal that these reactions proceed by addition of the silyl cation to the oxygen or nitrogen lone pair followed by a 1,3-H migration in the transition state. These transition states are calculated to lie below the energy of the reactants. By comparison, hydrolysis of gaseous CCl3+ is calculated to involve a substantial positive energy barrier.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Structural complexes of human adenine phosphoribosyltransferase reveal novel features of the APRT catalytic mechanism

Adenine phosphoribosyltransferase (APRT) is an important enzyme component of the purine recycling pathway. Parasitic protozoa of the order Kinetoplastida are unable to synthesize purines de novo and use the salvage pathway for the synthesis of purine bases rendering this biosynthetic pathway an attractive target for antiparasitic drug design. The recombinant human adenine phosphoribosyltransferase (hAPRT) structure was resolved in the presence of AMP in the active site to 1.76 angstrom resolution and with the substrates PRPP and adenine simultaneously bound to the catalytic site to 1.83 angstrom resolution. An additional structure was solved containing one subunit of the dimer in the apo-form to 2.10 angstrom resolution. Comparisons of these three hAPRT structures with other `type I` PRTases revealed several important features of this class of enzymes. Our data indicate that the flexible loop structure adopts an open conformation before and after binding of both substrates adenine and PRPR Comparative analyses presented here provide structural evidence to propose the role of Glu 104 as the residue that abstracts the proton of adenine N9 atom before its nucleophilic attack on the PRPP anomeric carbon. This work leads to new insights to the understanding of the APRT catalytic mechanism.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Evidence for magnetic phase separation in La(0.86)Sr(0.14)Mn(1-x)Cu(x)O(3+delta) manganites from NMR and magnetic measurements

Polycrystalline La(0.86)Sr(0.14)Mn(1-x)Cu(x)O(3+delta) (x = 0, 0.05, 0.10, 0.15, 0.20) manganites were investigated by means of magnetic measurements and zero-field (139)La and (55)Mn nuclear magnetic resonance (NMR) spectroscopy. Magnetization versus temperature measurements revealed a paramagnetic to ferromagnetic transition in most samples, with lower Curie temperatures and broader transitions for samples with higher Cu contents. The details of the magnetization measurements suggested a phase-separated scenario, with ferromagnetic clusters embedded in an antiferromagnetic matrix, especially for the samples with large Cu contents (x = 0.15 and 0.20). Zero-field (139)La NMR measurements confirmed this finding, since the spectral features remained almost unchanged for all Cu-doped samples, whereas the bulk magnetization was drastically reduced with increasing Cu content. (55)Mn NMR spectra were again typical of ferromagnetic regions, with a broadening of the resonance line caused by the disorder introduced by the Cu doping. The results indicate a coexistence of different magnetic phases in the manganites studied, with the addition of Cu contributing to the weakening of the double-exchange interaction in most parts of the material.

Subsite substrate specificity of midgut insect chymotrypsins

Insect chymotrypsins are distinctively sensitive to plant protein inhibitors, suggesting that they differ in subsite architecture and hence in substrate specificities. Purified digestive chymotrypsins from insects of three different orders were assayed with internally quenched fluorescent oligopeptides with three different amino acids at P1 (Tyr, Phe, and Leu) and 13 amino acid replacements in positions P1`, P2, and P3. The binding energy (Delta G(s), calculated from Km values) and the activation energy (Delta G(T)(double dagger), determined from k(cat)/K-m values) were calculated. The hydrophobicities of each subsite were calculated from the efficiency of hydrolysis of the different amino acid replacements at that subsite. The results showed that except for S1, the other subsites (S2, S3, and S1`) vary among chymotrypsins. This result contrasts with insect trypsin data that revealed a trend along evolution, putatively associated with resistance to plant inhibitors. In spite of those differences, the data suggested that in lepidopteran chymotrypsins S2 and S1` bind the substrate ground state, whereas only S1` binds the transition state, supporting aspects of the present accepted mechanism of catalysis. 2008 Elsevier Ltd. All rights reserved.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

2-Aminopurine non-radiative decay and emission in aqueous solution: A theoretical study

The minimum energy path along the lowest-lying pi pi* excited state of 2-aminopurine was calculated to elucidate the mechanisms of radiationless decay and emission in water. The sequential Monte Carlo quantum mechanics approach with a multiconfigurational and perturbative description of the wave function was employed to compute the minimum, transition state, and conical intersection. It was found that the barrier in the potential energy surface to access the conical intersection funnel increases in aqueous environment, making the system prone to enlarge the emission yield. These results rationalize the observed enhancement of emission in 2-aminopurine upon increasing of the solvent polarity. (c) 2008 Elsevier B.V. All rights reserved.

Gas-Phase Acylium Ion Transfer Reactions Mediated by a Proton Shuttle Mechanism

The mechanism and the energy profile of the gas-phase reaction that mimics esterification under acidic conditions have been investigated at different levels of theory. These reactions are known to proceed with rate constants close to the collision limit in the gas-phase and questions have been raised as to whether the typical addition-elimination mechanism via a tetrahedral intermediate can explain the ease of these processes. Because these reactions are common to many organic and biochemical processes it is important to understand the intrinsic reactivity of these systems. Our calculations at different levels of theory reveal that a stepwise mechanism via a tetrahedral species is characterized by energy barriers that are inconsistent with the experimental results. For the thermoneutral exchange between protonated acetic acid and water and the exothermic reaction of protonated acetic acid and methanol our calculations show that these reactions proceed initially by a proton shuttle between the carbonyl oxygen and the hydroxy oxygen of acetic acid mediated by water, or methanol, followed by displacement at the acylium ion center. These findings suggest that the reactions in the gas-phase should be viewed as an acylium ion transfer reaction. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1596-1606, 2011

Solvent assisted decomposition of the tetrahedral intermediate of the transesterification reaction to biodiesel production. A density functional study

B3LYP/6-31 + G(d) calculations were employed to investigate the mechanism of the transesterification reaction between a model monoglyceride and the methoxide and ethoxide anions. The gas-phase results reveal that both reactions have essentially the same activation energy (5.9 kcal mol(-1)) for decomposition of the key tetrahedral intermediate. Solvent effects were included by means of both microsolvation and the polarizable continuum solvation model CPCM. Both solvent approaches reduce the activation energy, however, only the microsolvation model is able to introduce some differentiation between methanol and ethanol, yielding a lower activation energy for decomposition of the tetrahedral intermediate in the reaction with methanol (1.1 kcal mol(-1)) than for the corresponding reaction with ethanol (2.8 kcal mol(-1)), in line with experimental evidences. Analysis of the individual energy components within the CPCM approach reveals that electrostatic interactions are the main contribution to stabilization of the transition state. (C) 2009 Elsevier Ltd. All rights reserved.