New insights about the monomer and homodimer structures of the human AOX1.

Human aldehyde oxidase (hAOX1) is a molybdenum dependent enzyme that plays an important role in the metabolism of various compounds either endogenous or xenobiotics. Due to its promiscuity, hAOX1 plays a major role in the pharmacokinetics of many drugs and therefore has gathered a lot of attention from the scientific community and, particularly, from the pharmaceutical industry. In this work, homology modelling, molecular docking and molecular dynamics simulations were used to study the structure of the monomer and dimer of human AOX. The results with the monomer of hAOX1 allowed to shed some light on the role played by thioridazine and two malonate ions that are co-crystalized in the recent X-ray structure of hAOX1. The results show that these molecules endorse several conformational rearrangements in the binding pocket of the enzyme and these changes have an impact in the active site topology as well as in the stability of the substrate (phthalazine). The results show that the presence of both molecules open two gates located at the entrance of the binding pocket, from which results the flooding of the active site. They also endorse several modifications in the shape of the binding pocket (namely the position of Lys893) that, together with the presence of the solvent molecules, favour the release of the substrate to the solvent. Further insights were also obtained with the assembled homodimer of hAOX1. The allosteric inhibitor (THI) binds closely to the region where the dimerization of both monomers occur. These findings suggest that THI can interfere with protein dimerization.

Protocol for Computational Enzymatic Reactivity Based on Geometry Optimisation.

Enzymes play a biologically essential role in performing and controlling an important share of the chemical processes occurring in life. However, despite their critical role in nature, attaining a clear understanding of the way an enzyme acts is still cumbersome. Computational enzymology is playing an increasingly important role in this field of research. It allows the elucidation of a complete and detailed mechanism of an enzymatic reaction, including the characterization of reaction intermediates and transition states from both structural and energetic points of view, which is something that no other single experiment can produce alone. In this review, we present a general computational strategy to study enzymatic mechanisms based on adiabatic mapping and free geometry optimization. These methods allow chemical reactions to be studied with high theoretical levels, and allow a more exhaustive exploration of the chemical reactional space than other available methods, albeit being limited to the extent that they explore the enzyme conformational space. Special attention is given to the choice of the theoretical levels, as well as describing the model systems that are currently used to study enzymatic reactions. With this, we aim to provide a good introduction for non-specialised users in this field of research. We also provide a selection of hand-picked examples from our own work that illustrate the power of computational enzymology to study catalytic mechanisms. Some of these studies constitute pioneering work in the field that were later validated by experimental means.

The mechanism of the Ser-(cis)Ser-Lys catalytic triad of peptide amidases.

In this paper, we report a theoretical investigation of the catalytic mechanism of peptide amidases that involve a Ser-(cis)Ser-Lys catalytic triad. Previous suggestions propose that these enzymes should follow a distinct catalytic mechanism from the one that is present in the classic Ser-His-Asp catalytic triad. The theoretical and computational results obtained in this work indicate the opposite idea, showing that both mechanisms are very similar and only few differences are observed between both reactions. The results reveal that the different alignment of the Ser-(cis)Ser-Lys catalytic triad in relation to the classical Ser-His-Asp triad may provide a better stabilisation of the reaction intermediates, and therefore make these enzymes catalytically more efficient. The catalytic mechanism has been determined at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory and requires five sequential steps instead of the two that are generally proposed: (i) nucleophilic attack of serine on the carbonyl group of the substrate, forming the first tetrahedral intermediate, (ii) formation of an acyl-enzyme complex, (ii) release of an ammonia product, (iv) nucleophilic attack of a water molecule forming the second tetrahedral intermediate, and (iv) the release of the product of the reaction, the carboxylic acid. The computational results suggest that the rate-limiting step is the first one that requires an activation free energy of 15.93 kcal mol-1. This result agrees very well with the available experimental data that predict a reaction rate of 2200 s-1, which corresponds to a free energy barrier of 14 kcal mol-1.

Amino acid deprivation using enzymes as a targeted therapy for cancer and viral infections.

INTRODUCTION: Amino acid depletion in the blood serum is currently being exploited and explored for therapies in tumors or viral infections that are auxotrophic for a certain amino acid or have a metabolic defect and cannot produce it. The success of these treatments is because normal cells remain unaltered since they are less demanding and/or can synthesize these compounds in sufficient amounts for their needs by other mechanisms. Areas covered: This review is focused on amino acid depriving enzymes and their formulations that have been successfully used in the treatment of several types of cancer and viral infections. Particular attention will be given to the enzymes L-asparaginase, L-arginase, L-arginine deiminase, and L-methionine-γ-lyase. Expert opinion: The immunogenicity and other toxic effects are perhaps the major limitations of these therapies, but they have been successfully decreased either through the expression of these enzymes from other organisms, recombination processes, pegylation of the selected enzymes or by specific mutations in the proteins. In 2006, FDA has already approved the use of L-asparaginase in the treatment of acute lymphoblastic leukemia. Other enzymes and in particular L-arginase, L-arginine deiminase, and L-methioninase have been showing promising results in vitro and in vivo studies.

The catalytic mechanism of carboxylesterases: a computational study.

The catalytic mechanism of carboxylesterases (CEs, EC 3.1.1.1) is explored by computational means. CEs hydrolyze ester, amide, and carbamate bonds found in xenobiotics and endobiotics. They can also perform transesterification, a reaction important, for instance, in cholesterol homeostasis. The catalytic mechanisms with three different substrates (ester, thioester, and amide) have been established at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory. It was found that the reactions proceed through a mechanism involving four steps instead of two as is generally proposed: (i) nucleophilic attack of serine to the substrate, forming the first tetrahedral intermediate, (ii) formation of the acyl-enzyme complex concomitant with the release of the alcohol product, (iii) nucleophilic attack of a water or alcohol molecule forming the second tetrahedral intermediate, and (iv) the release of the second product of the reaction. The results agree very well with the available experimental data and show that the hydrolytic and the transesterification reactions are competitive processes when the substrate is an ester. In all the other studied substrates (thioester or amide), the hydrolytic and transesterification process are less favorable and some of them might not even take place under in vivo conditions.

CHEM-PATH-TRACKER: An automated tool to analyze chemical motifs in molecular structures.

In this article, we propose a method for locating functionally relevant chemical motifs in protein structures. The chemical motifs can be a small group of residues or structure protein fragments with highly conserved properties that have important biological functions. However, the detection of chemical motifs is rather difficult because they often consist of a set of amino acid residues separated by long, variable regions, and they only come together to form a functional group when the protein is folded into its three-dimensional structure. Furthermore, the assemblage of these residues is often dependent on non-covalent interactions among the constituent amino acids that are difficult to detect or visualize. To simplify the analysis of these chemical motifs and give access to a generalized use for all users, we developed chem-path-tracker. This software is a VMD plug-in that allows the user to highlight and reveal potential chemical motifs requiring only a few selections. The analysis is based on atoms/residues pair distances applying a modified version of Dijkstra's algorithm, and it makes possible to monitor the distances of a large pathway, even during a molecular dynamics simulation. This tool turned out to be very useful, fast, and user-friendly in the performed tests. The chem-path-tracker package is distributed as an independent platform and can be found at http://www.fc.up.pt/PortoBioComp/database/doku.php?id=chem-path-tracker.

Discovery of new druggable sites in the anti-cholesterol target HMG-CoA reductase by computational alanine scanning mutagenesis.

The enzyme 3-hydroxy-3-methyl-glutaryl-CoA reductase (HMG-CoA-R) is the fundamental target for the treatment of hypercholesterolemia nowadays. The HMG-CoA-R clinical active site inhibitors (statins) are among the most widespread and profitable drugs ever sold but their side effects (myopathies, sometimes severe) still limit their use, which makes the finding of alternatives to statins a field of intense research. In this line, we address here a new strategy for inhibiting the homotetrameric HMG-CoA-R. The enzyme consists of a "dimer of dimers", each dimer having two active sites. We pursue here the inhibition of enzyme oligomerization, through drug binding to the dimer interface. We have computationally mutated 232 interfacial residues by alanine and calculated the loss in binding free energy among the monomers that build up each dimer of the homotetramer. This led to the identification of the (ten) key residues for the formation of the active dimer (Glu528, Ile531, Met534, Tyr644, Glu665, Asn686, Lys692, Lys735, Met742, and Val863). The results show that these residues are located in two specific spots of the protein with a cleft shape, whose shape and size is favorable for small drug binding. It is expectable that small molecules specifically bound to these druggable pockets will have a major effect on the oligomerization of the protein or/and in active site formation. This paves the way for the discovery of new families of inhibitors of HMG-CoA-R.

Gemcitabine: a critical nucleoside for cancer therapy.

Gemcitabine (dFdC, 2',2'-difluorodeoxycytidine) is a deoxycytidine nucleoside analogue of deoxycytidine in which two fluorine atoms have been inserted into the deoxyribose ring. Like other nucleoside analogues, gemcitabine is a prodrug. It is inactive in its original form, and depends on the intracellular machinery to gain pharmacological activity. What makes gemcitabine different from other nucleoside analogues is that it is actively transported across the cell membrane, it is phosphorylated more efficiently and it is eliminated at a slower rate. These differences, together with self-potentiation mechanisms, masked DNA chain termination and extensive inhibitory efficiency against several enzymes, are the source of gemcitabine's cytotoxic activity against a wide variety of tumors. This unique combination of metabolic properties and mechanistic characteristics is only found in very few other anticancer drugs, and both the FDA and the EMEA have already approved its use for clinical purposes, for the treatment of several types of tumors. In spite of the promising results associated with gemcitabine, the knowledge of its mode of action and of the enzymes it interacts with is still not fully documented. In this article we propose to review all these aspects and summarize the path of gemcitabine inside the cell.

Computational Mechanistic Studies Addressed to the Transimination Reaction Present in All Pyridoxal 5'-Phosphate-Requiring Enzymes.

The pyridoxal-5'-phosphate-dependent enzymes (PLP enzymes) catalyze a myriad of biochemical reactions, being actively involved in the biosynthesis of amino acids and amino acid-derived metabolites as well as in the biosynthetic pathways of amino sugars and in the synthesis or catabolism of neurotransmitters. Although the scope of PLP-catalyzed reactions initially appears to be bewilderingly diverse, there is a simple unifying principle: In the resting state, the cofactor (PLP) is covalently bonded to the amino group of an active site lysine, forming an internal aldimine. Once the amino substrate interacts with the active site, a new Schiff base is generated, commonly referred to as the external aldimine. Only after this step, the mechanistic pathway for each PLP-catalyzed reaction diverges. In this paper, density functional methods have been applied to investigate this common step present in all PLP-dependent enzymes-the transimination reaction. The results indicate that the reaction involves three sequential steps: (i) formation of a tetrahedral intermediate with the active site lysine and the amino substrate bonded to the PLP cofactor; (ii) nondirect proton transfer between the amino substrate and the lysine residue; and (iii) formation of the external aldimine after the dissociation of the lysine residue. The overall reaction is exothermic (-12.0 kcal/mol), and the rate-limiting step is the second one with 12.6 kcal/mol for the activation energy.

Ribonucleotide reductase: a mechanistic portrait of substrate analogues inhibitors.

Ribonucleotide reductase (RNR) is the key enzyme in the biosynthesis of deoxyribonucleotides. Several different strategies for inactivation of RNRs have been reported, including the use of substrate analogues as mechanism-based inhibitors. This article undergoes a critical analysis on the current status of ribonucleotide reductase inhibitory mechanisms by substrate analogues highlighting experimental and theoretical/computational approaches. We have summarized a general portrait of the inhibitory mechanisms and classified the nucleoside analogue inhibitors in three main classes. The critical analysis undertaken will contribute in finding new and more effective ways of inhibiting RNR.

The effect of the sixth sulfur ligand in the catalytic mechanism of periplasmic nitrate reductase.

The catalytic mechanism of nitrate reduction by periplasmic nitrate reductases has been investigated using theoretical and computational means. We have found that the nitrate molecule binds to the active site with the Mo ion in the +6 oxidation state. Electron transfer to the active site occurs only in the proton-electron transfer stage, where the Mo(V) species plays an important role in catalysis. The presence of the sulfur atom in the molybdenum coordination sphere creates a pseudo-dithiolene ligand that protects it from any direct attack from the solvent. Upon the nitrate binding there is a conformational rearrangement of this ring that allows the direct contact of the nitrate with Mo(VI) ion. This rearrangement is stabilized by the conserved methionines Met141 and Met308. The reduction of nitrate into nitrite occurs in the second step of the mechanism where the two dimethyl-dithiolene ligands have a key role in spreading the excess of negative charge near the Mo atom to make it available for the chemical reaction. The reaction involves the oxidation of the sulfur atoms and not of the molybdenum as previously suggested. The mechanism involves a molybdenum and sulfur-based redox chemistry instead of the currently accepted redox chemistry based only on the Mo ion. The second part of the mechanism involves two protonation steps that are promoted by the presence of Mo(V) species. Mo(VI) intermediates might also be present in this stage depending on the availability of protons and electrons. Once the water molecule is generated only the Mo(VI) species allow water molecule dissociation, and, the concomitant enzymatic turnover.

MADAMM: a multistaged docking with an automated molecular modeling protocol.

Dealing with receptor flexibility in docking methodology is still a problem. The main reason behind this difficulty is the large number of degrees of freedom that have to be considered in this kind of calculations. In this paper, we present an automated procedure, called MADAMM, that allows flexibilization of both the receptor and the ligand during a multistaged docking with an automated molecular modeling protocol. We show that the orientation of particular residues at the interface between the protein and the ligand have a crucial influence on the way they interact during the docking process, and the standard docking methodologies failed to predict their correct mode of binding. We present some examples that demonstrate the capabilities of this approach when compared with traditional docking methodologies.