Target-Related Applications of First Principles Quantum Chemical Methods in Drug Design

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

The N-terminal nucleophile (Ntn) hydrolases are a superfamily of enzymes specialized in the hydrolytic cleavage of amide bonds. Even though several members of this family are emerging as innovative drug targets for cancer, inflammation, and pain, the processes through which they catalyze amide hydrolysis remains poorly understood. In particular, the catalytic reactions of cysteine Ntn-hydrolases have never been investigated from a mechanistic point of view. In the present study, we used free energy simulations in the quantum mechanics/molecular mechanics framework to determine the reaction mechanism of amide hydrolysis catalyzed by the prototypical cysteine Ntn-hydrolase, conjugated bile acid hydrolase (CBAH). The computational analyses, which were confirmed in water and using different CBAH mutants, revealed the existence of a chair-like transition state, which might be one of the specific features of the catalytic cycle of Ntn-hydrolases. Our results offer new insights on Ntn-mediated hydrolysis and suggest possible strategies for the creation of therapeutically useful inhibitors.

The increasing role of QM/MM in drug discovery

Since its first appearance in 1976, the quantum mechanics/molecular mechanics (QM/MM) approach has mostly been used to study the chemical reactions of enzymes, which are frequently the target of drug discovery programs. In principle, a detailed understanding of the enzymatic mechanism should help researchers to design a potent enzyme inhibitor or new drug. However, QM/MM has not yet had a widespread impact on structure-based drug design. This is mostly due to its high computational cost. We expect this to change with the recent and extraordinary increases in computational power, and with the availability of more efficient algorithms for QM/MM calculations. Here, we report on some representative examples of QM/MM studies, including our own research, of pharmaceutically relevant enzymes, such as ribonuclease H and fatty acid amide hydrolase (FAAH). We aim to show how QM/MM has traditionally been used to study enzymatic catalysis. In this regard, we discuss its potential to become a routinely used drug design tool. To support this, we also discuss selected computational studies where QM/MM insights have been helpful in improving the potency of covalent inhibitors of FAAH.

Software and resources for computational medicinal chemistry

Computer-aided drug design plays a vital role in drug discovery and development and has become an indispensable tool in the pharmaceutical industry. Computational medicinal chemists can take advantage of all kinds of software and resources in the computer-aided drug design field for the purposes of discovering and optimizing biologically active compounds. This article reviews software and other resources related to computer-aided drug design approaches, putting particular emphasis on structure-based drug design, ligand-based drug design, chemical databases and chemoinformatics tools.

Application of a SCC-DFTB QM/MM approach to the investigation of the catalytic mechanism of fatty acid amide hydrolase

Self-consistent charge density functional tight binding (SCC-DFTB) is a promising method for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of enzyme-catalyzed reactions. The acylation reaction of fatty acid amide hydrolase (FAAH), a promising drug target, was investigated by applying a SCC-DFTB/CHARMM27 scheme. Calculated potential energy barriers resulted in reasonable agreement with experiments for oleamide (OA) and oleoylmethyl ester (OME) substrates, outperforming previous calculations performed at the PM3/CHARMM22 level. Furthermore, the experimental preference of FAAH in hydrolyzing OA faster than OME was adequately reproduced by calculations. All these findings indicate that the SCC-DFTB/CHARMM27 approach can be successfully applied to mechanistic investigations of FAAH-catalyzed reactions.

A reliable docking/scoring scheme based on the semiempirical quantum mechanical PM6-DH2 method accurately covering dispersion and H-bonding: HIV-1 protease with 22 ligands

In this study, we introduce a fast and reliable rescoring scheme for docked complexes based on a semiempirical quantum mechanical PM6-DH2 method. The method utilizes a PM6-based Hamiltonian with corrections for dispersion energy and hydrogen bonds. The total score is constructed as the sum of the PM6-DH2 interaction enthalpy, the empirical force field (AMBER) interaction entropy, and the sum of the deformation (PM6-DH2, SMD) and the desolvation (SMD) energies of the ligand. The main advantage of the procedure is the fact that we do not add any empirical parameter for either an individual component of the total score or an individual protein-ligand complex. This rescoring method is applied to a very challenging system, namely, the HIV-1 protease with a set of ligands. As opposed to the conventional DOCK procedure, the PM6-DH2 rescoring based on all of the terms distinguishes between binders and nonbinders and provides a reliable correlation of the theoretical and experimental binding free energies. Such a dramatic improvement, resulting from the PM6-DH2 rescoring of all the complexes, provides a valuable yet inexpensive tool for rational drug discovery and de novo ligand design.

Structural Fluctuations in Enzyme-Catalyzed Reactions: Determinants of Reactivity in Fatty Acid Amide Hydrolase from Multivariate Statistical Analysis of Quantum Mechanics/Molecular Mechanics Paths

The effects of structural fluctuations, due to protein dynamics, on enzyme activity are at the heart of current debates on enzyme catalysis. There is evidence that fatty acid amide hydrolase (FAAH) is an enzyme for which reaction proceeds via a high-energy, reactive conformation, distinct from the predominant enzyme-substrate complex (Lodola et al. Biophys. J. 2007, 92, L20-22). Identifying the structural causes of differences in reactivity between conformations in such complex systems is not trivial. Here, we show that multivariate analysis of key structural parameters can identify structural determinants of barrier height by analysis of multiple reaction paths. We apply a well-tested quantum mechanics/molecular mechanics (QM/MM) method to the first step of the acylation reaction between FAAH and oleamide substrate for 36 different starting structures. Geometrical parameters (consisting of the key bond distances that change during the reaction) were collected and used for principal component analysis (PCA), partial least-squares (PLS) regression analysis, and multiple linear regression (MLR) analysis. PCA indicates that different "families" of enzyme-substrate conformations arise from QM/MM molecular dynamics simulation and that rarely sampled, catalytically significant conformational states can be identified. PLS and MLR analyses allowed the construction of linear regression models, correlating the calculated activation barriers with simple geometrical descriptors. These analyses reveal the presence of two fully independent geometrical effects, explaining 78% of the variation in the activation barrier, which are directly correlated with transition-state stabilization (playing a major role in catalysis) and substrate binding. These results highlight the power of statistical approaches of this type in identifying crucial structural features that contribute to enzyme reactivity.

Computational characterization for catalytic activities of human CD38's wild type, E226 and E146 mutants

A series of the complexes of human CD38's wild type, E226 and E146 mutants as well have been simulated. The biosoftwares well simulate the penetration of nicotinamide-adenine-dinucleotide (NAD) into the active site. The nicotinamide end of NAD penetrates deep into the active site consistent with cleavage of the nicotinamide-glycosidic bond which is the first step of catalysis creating a Michaelis complex regarded as the intermediate product of NAD cyclase and hydrolysis reaction. The breaking down hydrogen bond between 2'-3' OH ribosyl and the residues replaced Glu(226) makes NAD to be less constrained in active site and nicotinamide (NA) becomes more difficult to be cleaved and eliminates the mutant catalytic activities. The large majority of the substrate NAD is hydrolyzed to ADPR while the conversion of NAD to cADPR is not the dominant reaction catalyzed by wild-type human CD38. The more strongly kept ribosyl group by hydrogen bonds the more NADase and the less cyclase activity. Breaking hydrogen bonds of ribosyl 2'- and 3'-OH by mutation will loosen it to promote the cyclase. The cyclic adenosine diphosphate-ribose (cADPR) could also penetrate deeply into active site to make some hydrogen bonds with Glu(146) and Glu(226); however, its docking poses are affected by a residue located at the entrance of the catalytic pocket (Lys(129)). These results are in good agreement with the previous crystallographic analysis and the experiments quantified the catalytic activities of human CD38 and its mutants.

Insights into docking and scoring neuronal alpha4beta2 nicotinic receptor agonists using molecular dynamics simulations and QM/MM calculations

A combined quantum mechanical (QM)-polarized docking and molecular dynamics approach to study the binding mode and to predict the binding affinity of ligands acting at the alpha4beta2-nAChR is presented. The results obtained in this study indicate that the quantum mechanical/molecular mechanics docking protocol well describes the charge-driven interactions occurring in the binding of nicotinic agonists, and it is able to represent the polarization effects on the ligand exerted by the surrounding atoms of the receptor at the binding site. This makes it possible to properly score agonists of alpha4beta2-nAChR and to reproduce the experimental binding affinity data with good accuracy, within a mean error of 2.2 kcal/mol. Moreover, applying the QM-polarized docking to an ensemble of nAChR conformations obtained from MD simulations enabled us to accurately capture nAChR-ligand induced-fit effects.

Water-chloride and water-bromide hydrogen-bonded networks: influence of the nature of the halide ions on the stability of the supramolecular assemblies

Two compounds, namely, [TTPH(2)](Cl)(2) x 4 H(2)O (1) and [TTPH(2)](Br)(2) x 4 H(2)O (2), (TTP = 4'-p-tolyl-2,2':6',2''-terpyridine) were synthesized from purely aqueous media and characterized by physical techniques. In the solid-state structures of these compounds, interesting supramolecular assemblies are observed. In 1, an unusual staircase-like architecture of the tape of edge-shared planar water hexamer is of importance, where the chloride ions are at the two edges of the tape. In 2, the polymeric nature of the water-bromide assembly is of interest, where discrete open-cube water octamers are doubly bridged by bromide ions. Semiempirical and DFT calculations confirm that the nature of the anion indeed affects the topology of the water-halide assemblies. We conclude that the protonated [TTPH(2)](2+) species can act as appropriate receptors for halide ions, which in turn act as a matrix for the formation of polymeric 1D water-halide assemblies.

Insights into the mechanism and inhibition of fatty acid amide hydrolase from quantum mechanics/molecular mechanics (QM/MM) modelling

FAAH (fatty acid amide hydrolase) is a promising target for the treatment of several central nervous system and peripheral disorders. Combined QM/MM (quantum mechanics/molecular mechanics) calculations have elucidated the role of its unusual catalytic triad in the hydrolysis of oleamide and oleoylmethyl ester substrates, and have identified the productive inhibitor-binding orientation for the carbamoylating compound URB524. These are potentially crucial insights for designing new covalent inhibitors of this drug target.

Interpretation of protein/ligand crystal structure using QM/MM calculations: case of HIV-1 protease/metallacarborane complex

Deltahedral metallacarborane compounds have recently been discovered as potent, specific, stable, and nontoxic inhibitors of HIV-1 protease (PR), the major target for AIDS therapy. The 2.15 A-resolution X-ray structure has exhibited a nonsymmetrical binding of the parental compound [Co(3+)-(C2B9H11)2](-) (GB-18) into PR dimer and a symmetrical arrangement in the crystal of two PR dimer complexes into a tetramer. In order to explore structural and energetic details of the inhibitor binding, quantum mechanics coupled with molecular mechanics approach was utilized. Realizing the close positioning of anionic inhibitors in the active site cavity, the possibility of an exchange of structural water molecules Wat50 and Wat128 by Na+ counterions was studied. The energy profiles for the rotation of the GB-18 molecules along their longitudinal axes in complex with PR were calculated. The results show that two Na+ counterions are present in the active site cavity and provide energetically favorable and unfavorable positions for carbon atoms within the carborane cages. Eighty-one rotamer combinations of four molecules of GB-18 bound to PR out of 4 x 10(5) are predicted to be highly populated. These results lay ground for further calculations of interaction energies between GB-18 and amino acids of PR active site and will make it possible to interpret computationally the binding of similar metallacarborane molecules to PR as well as to resistant PR variants. Moreover, this computational tool will allow the design of new, more potent metallacarborane-based HIV-1 protease inhibitors.

MurD ligase from E. coli: Tetrahedral intermediate formation study by hybrid quantum mechanical/molecular mechanical replica path method

MurD (UDP-N-acetylmuramoyl-L-alanine:D-glutamate ligase), a three-domain bacterial protein, catalyses a highly specific incorporation of D-glutamate to the cytoplasmic intermediate UDP-N-acetyl-muramoyl-L-alanine (UMA) utilizing ATP hydrolysis to ADP and P(i). This reaction is part of a biosynthetic path yielding bacterial peptidoglycan. On the basis of structural studies of MurD complexes, a stepwise catalytic mechanism was proposed that commences with a formation of the acyl-phosphate intermediate, followed by a nucleophilic attack of D-glutamate that, through the formation of a tetrahedral reaction intermediate and subsequent phosphate dissociation, affords the final product, UDP-N-acetyl-muramoyl-L-alanine-D-glutamate (UMAG). A hybrid quantum mechanical/molecular mechanical (QM/MM) molecular modeling approach was utilized, combining the B3LYP QM level of theory with empirical force field simulations to evaluate three possible reaction pathways leading to tetrahedral intermediate formation. Geometries of the starting structures based on crystallographic experimental data and tetrahedral intermediates were carefully examined together with a role of crucial amino acids and water molecules. The replica path method was used to generate the reaction pathways between the starting structures and the corresponding tetrahedral reaction intermediates, offering direct comparisons with a sequential kinetic mechanism and the available structural data for this enzyme. The acquired knowledge represents new and valuable information to assist in the ongoing efforts leading toward novel inhibitors of MurD as potential antibacterial drugs.

Predicting hydrogen-bond strengths from acid-base molecular properties. The pK(a) slide rule: toward the solution of a long-lasting problem

Unlike normal chemical bonds, hydrogen bonds (H-bonds) characteristically feature binding energies and contact distances that do not simply depend on the donor (D) and acceptor (:A) nature. Instead, their chemical context can lead to large variations even for a same donor-acceptor couple. As a striking example, the weak HO-H...OH(2) bond in neutral water changes, in acidic or basic medium, to the 6-fold stronger and 15% shorter [H(2)O...H...OH(2)](+) or [HO...H...OH](-) bonds. This surprising behavior, sometimes called the H-bond puzzle, practically prevents prediction of H-bond strengths from the properties of the interacting molecules. Explaining this puzzle has been the main research interest of our laboratory in the last 20 years. Our first contribution was the proposal of RAHB (resonance-assisted H-bond), a new type of strong H-bond where donor and acceptor are linked by a short pi-conjugated fragment. The RAHB discovery prompted new studies on strong H-bonds, finally leading to a general H-bond classification in six classes, called the six chemical leitmotifs, four of which include all known types of strong bonds. These studies attested to the covalent nature of the strong H-bond showing, by a formal valence-bond treatment, that weak H-bonds are basically electrostatic while stronger ones are mixtures of electrostatic and covalent contributions. The covalent component gradually increases as the difference of donor-acceptor proton affinities, DeltaPA, or acidic constants, DeltapK(a), approaches zero. At this limit, the strong and symmetrical D...H...A bonds formed can be viewed as true three-center-four-electron covalent bonds. These results emphasize the role PA/pK(a) equalization plays in strengthening the H-bond, a hypothesis often invoked in the past but never fully verified. In this Account, this hypothesis is reconsidered by using a new instrument, the pK(a) slide rule, a bar chart that reports in separate scales the pK(a)'s of the D-H proton donors and :A proton acceptors most frequently involved in D-H...:A bond formation. Allowing the two scales to shift so to bring selected donor and acceptor molecules into coincidence, the ruler permits graphical evaluation of DeltapK(a) and then empirical appreciation of the D-H...:A bond strength according to the pK(a) equalization principle. Reliability of pK(a) slide rule predictions has been verified by extensive comparison with two classical sources of H-bond strengths: (i) the gas-phase dissociation enthalpies of charged [X...H...X](-) and [X...H...X](+) bonds derived from the thermodynamic NIST Database and (ii) the geometries of more than 9500 H-bonds retrieved from the Cambridge Structural Database. The results attest that the pK(a) slide rule provides a reliable solution for the long-standing problem of H-bond-strength prediction and represents an efficient and practical tool for making such predictions directly accessible to all scientists.

Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism

Ribonuclease H (RNase H) belongs to the nucleotidyl-transferase (NT) superfamily and hydrolyzes the phosphodiester linkages that form the backbone of the RNA strand in RNA x DNA hybrids. This enzyme is implicated in replication initiation and DNA topology restoration and represents a very promising target for anti-HIV drug design. Structural information has been provided by high-resolution crystal structures of the complex RNase H/RNA x DNA from Bacillus halodurans (Bh), which reveals that two metal ions are required for formation of a catalytic active complex. Here, we use classical force field-based and quantum mechanics/molecular mechanics calculations for modeling the nucleotidyl transfer reaction in RNase H, clarifying the role of the metal ions and the nature of the nucleophile (water versus hydroxide ion). During the catalysis, the two metal ions act cooperatively, facilitating nucleophile formation and stabilizing both transition state and leaving group. Importantly, the two Mg(2+) metals also support the formation of a meta-stable phosphorane intermediate along the reaction, which resembles the phosphorane intermediate structure obtained only in the debated beta-phosphoglucomutase crystal (Lahiri, S. D.; et al. Science 2003, 299 (5615), 2067-2071). The nucleophile formation (i.e., water deprotonation) can be achieved in situ, after migration of one proton from the water to the scissile phosphate in the transition state. This proton transfer is actually mediated by solvation water molecules. Due to the highly conserved nature of the enzymatic bimetal motif, these results might also be relevant for structurally similar enzymes belonging to the NT superfamily.

Is quantum mechanics necessary for predicting binding free energy?

To take into account polarization effects, the linear interaction energy model with continuum electrostatic solvation (LIECE) is supplemented by the linear-scaling semiempirical quantum mechanical calculation of the intermolecular electrostatic energy (QMLIECE). QMLIECE and LIECE are compared on three enzymes belonging to different classes: the West Nile virus NS3 serine protease (WNV PR), the aspartic protease of the human immunodeficiency virus (HIV-1 PR), and the human cyclin-dependent kinase 2 (CDK2). QMLIECE is superior for 44 peptidic inhibitors of WNV PR because of the different amount of polarization due to the broad range of formal charges of the inhibitors (from 0 to 3). On the other hand, QMLIECE and LIECE show similar accuracy for 24 peptidic inhibitors of HIV-1 PR (20 neutral and 4 with one formal charge) and for 73 CDK2 inhibitors (all neutral). These results indicate that quantum mechanics is essential when the inhibitor/protein complexes have highly variable charge-charge interactions.

Atomistic insights into the inhibition of cysteine proteases: first QM/MM calculations clarifying the regiospecificity and the inhibition potency of epoxide- and aziridine-based inhibitors

Epoxides and aziridines are important building blocks for inhibitors of cysteine proteases which are promising drug targets for many diseases. In spite of the large amount of experimental data concerning inhibition potency, structure-activity relationships, and structural arrangements of enzyme-inhibitor complexes, little is known about the basic principles which connect the substitution pattern with the resulting activities. To shed some light on this issue which is essential for the rational design of improved compounds, we have studied the inhibition processes theoretically for various inhibitors using quantum mechanical/molecular mechanical hybrid approaches and classical molecular dynamics simulations. The careful analysis of the computational results allows insight into the interactions which govern the regio- and stereospecificity of the interactions. Known structure-activity relationships are rationalized in terms of the same interactions that determine the measured pH dependencies. Inconsistencies in existing X-ray structures are resolved through comparison with the computed structures, which leads to a reassessment of the factors that control the inhibition potency. Similarities and differences in the mode of action of epoxide- and aziridine-based inhibitors are elucidated. Finally the small reaction barriers computed for the irreversible step in E64 analogues call into question the commonly accepted two-step model of inhibition since the second, irreversible step is predicted to be so fast that suitably oriented enzyme-inhibitor complexes will react rather than dissociate and equilibrate.

Computational Study of the Phosphoryl Transfer Catalyzed by a Cyclin-Dependent Kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.

The role of quantum mechanics in structure-based drug design

Herein we will focus on the use of quantum mechanics (QM) in drug design (DD) to solve disparate problems from scoring protein-ligand poses to building QM QSAR models. Through the variational principle of QM we know that we can obtain a more accurate representation of molecular systems than classical models, and while this is not a matter of debate, it still has not been shown that the expense of QM approaches is offset by improved accuracy in DD applications. Objectively validating the improved applicability and performance of QM over classical-based models in DD will be the focus of research in the coming years along with research on the conformational sampling problem as it relates to protein-ligand complexes.

A Linear Scaling Study of Solvent−Solute Interaction Energy of Drug Molecules in Aqua Solution

Solvent-solute interaction energies for three well-known drug molecules in water solution are computed at the Hartree-Fock and B3LYP density functional theory levels using a linear scaling technique, which allows one to explicitly include in the model water molecules up to 14 A away from the solute molecule. The dependence of calculated interaction energies on the amount of included solvent has been examined. It is found that it is necessary to account for water molecules within an 8 A radius around the drug molecule to reach the saturated solvent interaction level. Effects of electron correlation and basis set on solvent-solute interaction energies are discussed.

Extending the Power of Quantum Chemistry to Large Systems with the Fragment Molecular Orbital Method

Following the brief review of the modern fragment-based methods and other approaches to perform quantum-mechanical calculations of large systems, the theoretical development of the fragment molecular orbital method (FMO) is covered in detail, with the emphasis on the physical properties, which can be computed with FMO. The FMO-based polarizable continuum model (PCM) for treating the solvent effects in large systems and the pair interaction energy decomposition analysis (PIEDA) are described in some detail, and a range of applications of FMO to biological studies is introduced. The factors determining the relative stability of polypeptide conformers (alpha-helix, beta-turn, and extended form) are elucidated using FMO/PCM and PIEDA, and the interactions in the Trp-cage miniprotein construct (PDB: 1L2Y) are analyzed using PIEDA.

Inhibitor Binding by Metallo-β-lactamase IMP-1 from Pseudomonas aeruginosa: Quantum Mechanical/Molecular Mechanical Simulations

The dynamics of the IMP-1 enzyme complexed with three prototypical inhibitors are investigated using a quantum mechanical/molecular mechanical (QM/MM) method based on the self-consistent-charge density-functional tight-binding model. The binding patterns of the inhibitors observed in X-ray diffraction experiments are well reproduced in 600 ps molecular dynamics simulations at room temperature. These inhibitors anchor themselves in the enzyme active site by direct coordination with the two zinc ions, displacing the hydroxide nucleophile that bridges the two zinc ions. In addition, they also interact with several active-site residues and those in two mobile loops. The excellent agreement with experimental structural data validates the QM/MM treatment used in our simulations.

Comparative assessment of theoretical methods for the determination of geometrical properties in biological zinc complexes

In the present study, we have compared the performance of the density functional theory (DFT) functionals B1B95, B3LYP, B97-2, BP86, and BPW91 with MP2 for geometry determination in biological mononuclear Zn complexes. A total of 15 different basis sets, of rather diverse complexity, were tested, several which included also three different types of common effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. In addition, the ability to describe mononuclear Zn biological systems using relatively simple models of the metal coordination sphere, comprising only the metal atom and a simplified representation of the ligands at the first coordination sphere, starting from a set of high-resolution X-ray crystallographic structures, is evaluated for 90 combinations of method/basis set. The results show that the use of such models allows for a relatively accurate description of the Zn-ligand bond lengths, although failing to correctly represent the topology of the metal coordination sphere (namely, the angles involving the metal atom) if constraints at the Calpha atoms are not considered. Globally, B3LYP had the best average performance in the test, closely followed by MP2, whereas B1B95 was the least accurate method. The study also points out B3LYP/CEP-121G and B3LYP/SDD, which use, respectively, the Steven-Basch-Krauss and the Stuttgart-Dresden effective-core potentials, as the best compromise between accuracy and CPU time for the geometrical characterization of metal-ligand bond lengths in Zn biological systems.

Impact of pyrophosphate and O-ethyl-substituted pyrophosphate groups on DNA structure

Design of novel DNA probes to inhibit specific repair pathways is important for basic science applications and for use as therapeutic agents. As shown previously, single pyrophosphate (PP) and O-ethyl-substituted pyrophosphate (SPP) modifications can inhibit the DNA glycosylase activities on damaged DNA. To understand the structural basis of this inhibition, the influence of the PP and SPP internucleotide groups on the helical parameters and geometry of a double-stranded DNA was studied by using molecular modeling tools including molecular dynamics and quantum mechanical-molecular mechanical (QM/MM) approaches. Native and locally modified PP- and SPP-containing DNA duplexes of dodecanucleotide d(C1G2C3G4A5A6T7T8C9G10C11G12) were simulated in aqueous solution. The energies and forces were computed by using the PBE0/6-31+G** approach in the QM part and the AMBER force-field parameters in the MM part. Analysis of the local base-pair helical parameters, internucleotide distances, and overall global structure at the located stationary points revealed a close similarity of the initial and modified duplexes, with only torsion angles of the main chain being altered in the vicinity of introduced chemical modification. Results show that the PP and SPP groups are built into a helix structure without elongation of the internucleotide distance due to flipping-out of phosphate group from the sugar-phosphate backbone. The mechanism of such embedding has only a minor impact on the base pairs stacking and Watson-Crick interactions. Biochemical studies revealed that the PP and SPP groups immediately 5', but not 3', to the 8-oxoguanosine (8oxodG) inhibit translesion synthesis by a DNA polymerase in vitro. These results suggest that subtle perturbations of the DNA backbone conformation influence processing of base lesions.

Antibiotic Binding to Dizinc β-Lactamase L1 from Stenotrophomonas maltophilia: SCC-DFTB/CHARMM and DFT Studies

A dizinc beta-lactamase (L1 from Stenotrophomonas maltophilia) complexed with an antibiotic compound (moxalactam) has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach. The QM region is described by the self-consistent charge-density functional tight binding (SCC-DFTB) model while the MM by CHARMM. The Michaelis complex, which is constructed from a recent X-ray structure of the L1 enzyme with the hydrolyzed moxalactam, is simulated by molecular dynamics. The simulation yields valuable insights into substrate-enzyme interaction, whose implications in the enzyme catalysis are discussed. Finally, the QM/MM results are compared with a high-level density functional theory study of a truncated active-site model and the agreement provides strong support for the SCC-DFTB treatment of the QM region.