We present a density functional theory (DFT) hybrid quantum mechanical/molecular mechanical (QM-MM) implementation developed for simulations of reactions in complex environments. It is particularly suited to study enzyme active sites or solutes in condensed phases. The method combines a QM description of the solute with a MM treatment of the environment. The QM fragment is treated using DFT as implemented in the computationally efficient program SIESTA, while the environment is treated using the Wang et al. Amber force field parametrization. We applied our new QM-MM scheme to study the conversion of chorismate to prephenate by computing the reaction energy profile in vacuo, aqueous solution and in the active site of the B. subtilis chorismate mutase enzyme. We have performed calculations for two different choices of the QM subsystem in the enzyme simulations: including only the substrate moiety and the substrate plus the charged side chains glu78 and arg90, respectively. In both cases, our results are in good agreement with experiment. The catalytic activity achieved by chorismate mutase relative to the uncatalyzed reaction in solution is due to both a minor destabilization of the substrate molecule by compression and a major electrostatic stabilization of the transition state, which reduce the activation energy of the reaction.

Mycobacterium tuberculosis is the causative agent of human tuberculosis. The nitric oxide reaction with oxy-truncated hemoglobin N (trHbN) has been proposed to be responsible for the resistance mechanism by which this microorganism can evade the toxic effects of NO. In this work, we explore the molecular basis of the NO detoxification mechanism using a combination of classical and hybrid quantum-classical (QM-MM) simulation techniques. We have investigated the structural flexibility of the protein, the ligand affinity properties, and the nitric oxide reaction with coordinated O2. The analysis of the classical MD trajectory allowed us to identify Phe62 as the gate of the main channel for ligand diffusion to the active site. Moreover, the opening of the channel stems from the interplay between collective backbone motions and local rearrangements in the side chains of the residues that form the bottleneck of the tunnel. Even though the protein environment is not found to make a significant contribution to the heme moiety catalyzed reaction, the binding site influences the physiological function of the enzyme at three different levels. First, by isolating the intermediates formed in the reaction, it prevents nondesired reactions from proceeding. Second, it modulates the ligand (O2, NO) affinity of the protein, which can be ascribed to both distal and proximal effects. Finally, the stabilization of the Tyr33-Gln58 pair upon O2 binding might alter the essential dynamics of the protein, leading in turn to a mechanism for ligand-induced regulation.

The molecular basis of the hydroxylation reaction of the Calpha of a C-terminal glycine catalyzed by peptidylglycine alpha-hydroxylating monooxygenase (PHM) was investigated using hybrid quantum-classical (QM-MM) computational techniques. We have identified the most reactive oxygenated species and presented new insights into the hydrogen abstraction (H-abstraction) mechanism operative in PHM. Our results suggest that O(2) binds to Cu(B) to generate Cu(B)(II)-O(2)(.-) followed by electron transfer (ET) from Cu(A) to form Cu(B)(I)-O(2)(.-). The computed potential energy profiles for the H-abstraction reaction for Cu(B)(II)-O(2)(.-), Cu(B)(I)-O(2)(.-), and [Cu(B)(II)-OOH](+) species indicate that none of these species can be responsible for abstraction. However, the latter species can spontaneously form [Cu(B)O](+2) (which consists of a two-unpaired-electrons [Cu(B)O](+) moiety ferromagnetically coupled with a radical cation located over the three Cu(B) ligands, in the quartet spin ground state) by abstracting a proton from the surrounding solvent. Both this monooxygenated species and the one obtained by reduction with ascorbate, [Cu(B)O](+), were found to be capable of carrying out the H-abstraction; however, whereas the former abstracts the hydrogen atom concertedly with almost no activation energy, the later forms an intermediate that continues the reaction by a rebinding step. We propose that the active species in H-abstraction in PHM is probably [Cu(B)O](+2) because it is formed exothermically and can concertedly abstract the substrate HA atom with the lower overall activation energy. Interestingly, this species resembles the active oxidant in cytochrome P450 enzymes, Compound I, suggesting that both PHM and cytochrome P450 enzymes may carry out substrate hydroxylation by using a similar mechanism.

Mycobacterium tuberculosis, the causative agent of human tuberculosis, is forced into latency by nitric oxide produced by macrophages during infection. In response to nitrosative stress M. tuberculosis has evolved a defense mechanism that relies on the oxygenated form of "truncated hemoglobin" N (trHbN), formally acting as NO-dioxygenase, yielding the harmless nitrate ion. X-ray crystal structures have shown that trHbN hosts a two-branched protein matrix tunnel system, proposed to control diatomic ligand migration to the heme, as the rate-limiting step in NO conversion to nitrate. Extended molecular dynamics simulations (0.1 micros), employed here to characterize the factors controlling diatomic ligand diffusion through the apolar tunnel system, suggest that O2 migration in deoxy-trHbN is restricted to a short branch of the tunnel, and that O2 binding to the heme drives conformational and dynamical fluctuations promoting NO migration through the long tunnel branch. The simulation results suggest that trHbN has evolved a dual-path mechanism for migration of O2 and NO to the heme, to achieve the most efficient NO detoxification.

A novel technique for computing free energy profiles in enzymatic reactions using the multiple steering molecular dynamics approach in the context of an efficient QM-MM density functional scheme is presented. The conversion reaction of chorismate to prephenate catalyzed by the Bacillus subtilis enzyme chorismate mutase has been chosen as an illustrative example.

Targeting kinases is central to drug-based cancer therapy but remains challenging because the drugs often lack specificity, which may cause toxic side effects. Modulating side effects is difficult because kinases are evolutionarily and hence structurally related. The lack of specificity of the anticancer drug imatinib enables it to be used to treat chronic myeloid leukemia, where its target is the Bcr-Abl kinase, as well as a proportion of gastrointestinal stromal tumors (GISTs), where its target is the C-Kit kinase. However, imatinib also has cardiotoxic effects traceable to its impact on the C-Abl kinase. Motivated by this finding, we made a modification to imatinib that hampers Bcr-Abl inhibition; refocuses the impact on the C-Kit kinase; and promotes inhibition of an additional target, JNK, a change that is required to reinforce prevention of cardiotoxicity. We established the molecular blueprint for target discrimination in vitro using spectrophotometric and colorimetric assays and through a phage-displayed kinase screening library. We demonstrated controlled inhibitory impact on C-Kit kinase in human cell lines and established the therapeutic impact of the engineered compound in a novel GIST mouse model, revealing a marked reduction of cardiotoxicity. These findings identify the reengineered imatinib as an agent to treat GISTs with curbed side effects and reveal a bottom-up approach to control drug specificity.

No abstract

Heme proteins are found in all living organisms and are capable of performing a wide variety of tasks, requiring in many cases the binding of diatomic ligands, namely, O(2), CO, and/or NO. Therefore, subtle regulation of these diatomic ligands' affinity is one of the key issues for determining a heme protein's function. This regulation is achieved through direct H-bond interactions between the bound ligand and the protein, and by subtle tuning of the intrinsic heme group reactivity. In this work, we present an investigation of the proximal regulation of oxygen affinity in Fe(II) histidine coordinated heme proteins by means of computer simulation. Density functional theory calculations on heme model systems are used to analyze three proximal effects: charge donation, rotational position, and distance to the heme porphyrin plane of the proximal histidine. In addition, hybrid quantum-classical (QM-MM) calculations were performed in two representative proteins: myoglobin and leghemoglobin. Our results show that all three effects are capable of tuning the Fe-O(2) bond strength in a cooperative way, consistently with the experimental data on oxygen affinity. The proximal effects described herein could operate in a large variety of O(2)-binding heme proteins-in combination with distal effects-and are essential to understand the factors determining a heme protein's O(2) affinity.