Comparative Computational Analysis of Different Active Site Conformations and Substrates in a Chalcone Isomerase Catalyzed Reaction

AN AB INITIO STUDY ON THE STRUCTURE–CYTOTOXICITY RELATIONSHIP OF TERPENOID LACTONES BASED ON THE MICHAEL REACTION BETWEEN THEIR PHARMACOPHORES AND L-CYSTEINE-METHYLESTER-1

The cytotoxic effects of terpenoid lactones are attributed to the alkylation of biological nucleophiles, especially sulfhydryl groups in proteins, by the α,β-unsaturated carbonyl moiety of lactones through Michael reaction. Therefore, the cytotoxicity could be reflected by the reactivity of the pharmacophores. In this work, the Michael reaction between 12 α,β-unsaturated-carbonyl-containing small species, i.e. 10 analogues of the alpha methylene gamma butyrolactone moiety of andrographolide, one cyclopentenone, and one methylene–pentanolide, and L-cysteine-methylester-1 were investigated by ab initio methods to mimic the alkylation of proteins by terpenoid lactones. The trend in the calculated reaction free energies of the small species is qualitatively in accordance with the reported cytotoxicity of corresponding terpenoid lactones.

Biosynthesis of ethers: unusual or common natural events?

Ether bonds are found in a wide variety of natural products--mainly secondary metabolites--including lipids, oxiranes, terpenoids, flavonoids, polyketides, and carbohydrate derivatives, to name some representative examples. To furnish such a biodiversity of structures, a large number of different enzymes are involved in several different biosynthetic pathways. Depending on the compound and on the (micro) environment in which the reaction is performed, ethers are produced by very different (enzymatic) reactions, thus providing an impressive display of how Nature has combined evolution and thermodynamics to be able to produce a vast number of compounds. In addition, many of these compounds possess different biological activities of pharmacological interest. Moreover, some of these ethers (i.e., epoxides) have high chemical reactivity, and can be useful starting materials for further synthetic processes. This review aims to provide an overview of the different strategies that are found in Nature for the formation of these "bioethers". Both fundamental and practical insights of the biosynthetic processes will be discussed.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.