Estudo do coeficiente de difusão no enovelamento de proteína

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Visualização do funil de enovelamento de proteínas

Pós-graduação em Biofísica Molecular - IBILCE

Visualization of Protein Folding Funnels in Lattice Models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo da interação não-nativa no enovelamento de proteínas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudos computacionais de esfingomielinases D: docking, dinâmica molecular e métodos híbridos QM/MM

Pós-graduação em Biofísica Molecular - IBILCE

Estudos computacionais de esfingomielinases D: docking, dinâmica molecular e métodos híbridos QM/MM

Pós-graduação em Biofísica Molecular - IBILCE

Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

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.

Kinetics and product distribution of p-nitrophenyl phosphate dianion solvolysis in ternary DMSO/alcohol/water mixtures are compatible with metaphosphate formation

The rate of solvolysis of p-nitrophenyl phosphate (PNPP) dianion in DMSO/water strongly decreases by increasing water concentration. Addition of linear alcohols (methanol, propanol, butanol, pentanol, and hexanol) at constant DMSO/water molar ratio produced an even sharper rate decrease. Alkyl phosphate formation, resulting from PNPP solvolysis in ternary DMSO/water/alcohol mixtures, increased with alcohol concentration and was essentially temperature independent. Methanol and hexanol were the poorest nucleophiles under all conditions. Activation energies and enthalpies for solvolysis in ternary mixtures were similar and entropies varied with alcohol concentration. Taken together these results can be best interpreted in terms of a dissociative mechanism with the intervention of metaphosphate. Copyright (C) 2011 John Wiley & Sons, Ltd.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

X-ray diffraction and theoretical investigation of the Gedunin crystal structure

The Gedunin compound (C28H34O6) is a natural product extracted from Trichilia pallida that has shown a wide activity. The crystallographic structure shows two conformers in the asymmetric unit, which differ in a rotation of the furan group. To understand this molecular arrangement, the density functional calculations. Molecular Electrostatic Potential (MEP) and thermodynamic function calculation have been performed at the B3LYP/6-311++g(d,p) level. Both conformers were optimized and the agreement with the experimental structure was very good, making possible further theoretical analysis of the structure. The inter-conversion between two conformers depends on the energy barrier. This process is studied in the vacuum and shows two transition states with a low energetic barrier for a potential energy curve scanning rigid around furan group: 4.37 kcal/mol and 16.52 kcal/mol. As the first transition state has a notably lower energetic barrier, the preferred inter-conversion pathway between the conformers involves the first rather than the second transition state. Understanding this transition state in detail led us to perform its optimization, showing an energetic barrier around 3.66 kcal/mol. The negative free energy and low enthalpy confirm that the process is spontaneous and exothermic. The results show that this requirement makes the existence of the two conformers in the asymmetric unit possible. The structure of molecules in the asymmetric unit is better understood when the MEP is used on the interaction between molecules. For Gedunin, both molecules have shown MEP with well-defined regions, and this behavior contributes to the observed link between molecules and for the negative regions complementing positive regions of another molecule. (C) 2011 Elsevier B.V. All rights reserved.

Design and synthesis of novel non peptidomimtic beta-secretase inhibitors in the treatment of Alzheimer's disease

The aspartic protease BACE1 (β-amyloid precursor protein cleaving enzyme, β-secretase) is recognized as one of the most promising targets in the treatment of Alzheimer's disease (AD). The accumulation of β-amyloid peptide (Aβ) in the brain is a major factor in the pathogenesis of AD. Aβ is formed by initial cleavage of β-amyloid precursor protein (APP) by β-secretase, therefore BACE1 inhibition represents one of the therapeutic approaches to control progression of AD, by preventing the abnormal generation of Aβ. For this reason, in the last decade, many research efforts have focused at the identification of new BACE1 inhibitors as drug candidates. Generally, BACE1 inhibitors are grouped into two families: substrate-based inhibitors, designed as peptidomimetic inhibitors, and non-peptidomimetic ones. The research on non-peptidomimetic small molecules BACE1 inhibitors remains the most interesting approach, since these compounds hold an improved bioavailability after systemic administration, due to a good blood-brain barrier permeability in comparison to peptidomimetic inhibitors. Very recently, our research group discovered a new promising lead compound for the treatment of AD, named lipocrine, a hybrid derivative between lipoic acid and the AChE inhibitor (AChEI) tacrine, characterized by a tetrahydroacridinic moiety. Lipocrine is one of the first compounds able to inhibit the catalytic activity of AChE and AChE-induced amyloid-β aggregation and to protect against reactive oxygen species. Due to this interesting profile, lipocrine was also evaluated for BACE1 inhibitory activity, resulting in a potent lead compound for BACE1 inhibition. Starting from this interesting profile, a series of tetrahydroacridine analogues were synthesised varying the chain length between the two fragments. Moreover, following the approach of combining in a single molecule two different pharmacophores, we designed and synthesised different compounds bearing the moieties of known AChEIs (rivastigmine and caproctamine) coupled with lipoic acid, since it was shown that dithiolane group is an important structural feature of lipocrine for the optimal inhibition of BACE1. All the tetrahydroacridines, rivastigmine and caproctamine-based compounds, were evaluated for BACE1 inhibitory activity in a FRET (fluorescence resonance energy transfer) enzymatic assay (test A). With the aim to enhancing the biological activity of the lead compound, we applied the molecular simplification approach to design and synthesize novel heterocyclic compounds related to lipocrine, in which the tetrahydroacridine moiety was replaced by 4-amino-quinoline or 4-amino-quinazoline rings. All the synthesized compounds were also evaluated in a modified FRET enzymatic assay (test B), changing the fluorescent substrate for enzymatic BACE1 cleavage. This test method guided deep structure-activity relationships for BACE1 inhibition on the most promising quinazoline-based derivatives. By varying the substituent on the 2-position of the quinazoline ring and by replacing the lipoic acid residue in lateral chain with different moieties (i.e. trans-ferulic acid, a known antioxidant molecule), a series of quinazoline derivatives were obtained. In order to confirm inhibitory activity of the most active compounds, they were evaluated with a third FRET assay (test C) which, surprisingly, did not confirm the previous good activity profiles. An evaluation study of kinetic parameters of the three assays revealed that method C is endowed with the best specificity and enzymatic efficiency. Biological evaluation of the modified 2,4-diamino-quinazoline derivatives measured through the method C, allow to obtain a new lead compound bearing the trans-ferulic acid residue coupled to 2,4-diamino-quinazoline core endowed with a good BACE1 inhibitory activity (IC50 = 0.8 mM). We reported on the variability of the results in the three different FRET assays that are known to have some disadvantages in term of interference rates that are strongly dependent on compound properties. The observed results variability could be also ascribed to different enzyme origin, varied substrate and different fluorescent groups. The inhibitors should be tested on a parallel screening in order to have a more reliable data prior to be tested into cellular assay. With this aim, preliminary cellular BACE1 inhibition assay carried out on lipocrine confirmed a good cellular activity profile (EC50 = 3.7 mM) strengthening the idea to find a small molecule non-peptidomimetic compound as BACE1 inhibitor. In conclusion, the present study allowed to identify a new lead compound endowed with BACE1 inhibitory activity in submicromolar range. Further lead optimization to the obtained derivative is needed in order to obtain a more potent and a selective BACE1 inhibitor based on 2,4-diamino-quinazoline scaffold. A side project related to the synthesis of novel enzymatic inhibitors of BACE1 in order to explore the pseudopeptidic transition-state isosteres chemistry was carried out during research stage at Università de Montrèal (Canada) in Hanessian's group. The aim of this work has been the synthesis of the δ-aminocyclohexane carboxylic acid motif with stereochemically defined substitution to incorporating such a constrained core in potential BACE1 inhibitors. This fragment, endowed with reduced peptidic character, is not known in the context of peptidomimetic design. In particular, we envisioned an alternative route based on an organocatalytic asymmetric conjugate addition of nitroalkanes to cyclohexenone in presence of D-proline and trans-2,5-dimethylpiperazine. The enantioenriched obtained 3-(α-nitroalkyl)-cyclohexanones were further functionalized to give the corresponding δ-nitroalkyl cyclohexane carboxylic acids. These intermediates were elaborated to the target structures 3-(α-aminoalkyl)-1-cyclohexane carboxylic acids in a new readily accessible way.