Docking ligands into flexible and solvated macromolecules. 6. Development and application to the docking of HDACs and other zinc metalloenzymes inhibitors

Elucidation of the binding mode of organic polysulfides on the human TRPA1 receptor

Introduction: Previous studies have established that endogenous inorganic polysulfides have significant biological actions activating the Transient Receptor Potential Ankyrin 1 (TRPA1) receptor. Organic polysulfides exert similar effects, but they are much more stable molecules, therefore these compounds are more suitable as drugs. In this study, we aimed to better understand the mechanism of action of organic polysulfides by identification of their binding site on the TRPA1 receptor. Methods: Polysulfides can readily interact with the thiol side chain of the cysteine residues of the protein. To investigate their role in the TRPA1 activation, we replaced several cysteine residues by alanine via site-directed mutagenesis. We searched for TRPA1 mutant variants with decreased or lost activating effect of the polysulfides, but with other functions remaining intact (such as the effects of non-electrophilic agonists and antagonists). The binding properties of the mutant receptors were analyzed by in silico molecular docking. Functional changes were tested by in vitro methods: calcium sensitive fluorescent flow cytometry, whole-cell patch-clamp and radioactive calcium-45 liquid scintillation counting. Results: The cysteines forming the conventional binding site of electrophilic agonists, namely C621, C641 and C665 also bind the organic polysulfides, with the key role of C621. However, only their combined mutation abolished completely the organic polysulfide-induced activation of the receptor. Discussion: Since previous papers provided evidence that organic polysulfides exert analgesic and anti-inflammatory actions in different in vivo animal models, we anticipate that the development of TRPA1-targeted, organic polysulfide-based drugs will be promoted by this identification of the binding site.

Docking Ligands into Flexible and Solvated Macromolecules. 8. Forming New Bonds─Challenges and Opportunities

Over the years, structure-based design programs and specifically docking small molecules to proteins have become prominent in drug discovery. However, many of these computational tools have been developed to primarily dock enzyme inhibitors (and ligands to other protein classes) relying heavily on hydrogen bonds and electrostatic and hydrophobic interactions. In reality, many drug targets either feature metal ions, can be targeted covalently, or are simply not even proteins (e.g., nucleic acids). Herein, we describe several new features that we have implemented into Fitted to broaden its applicability to a wide range of covalent enzyme inhibitors and to metalloenzymes, where metal coordination is essential for drug binding. This updated version of our docking program was tested for its ability to predict the correct binding mode of drug-sized molecules in a large variety of proteins. We also report new datasets that were essential to demonstrate areas of success and those where additional efforts are required. This resource could be used by other program developers to assess their own software.

Inhibition of Collagenase Q1 of Bacillus cereus as a Novel Antivirulence Strategy for the Treatment of Skin‐Wound Infections

Despite the progress in surgical techniques and antibiotic prophylaxis, opportunistic wound infections with Bacillus cereus remain a public health problem. Secreted toxins are one of the main factors contributing to B. cereus pathogenicity. A promising strategy to treat such infections is to target these toxins and not the bacteria. Although the exoenzymes produced by B. cereus are thoroughly investigated, little is known about the role of B. cereus collagenases in wound infections.

In this report, the collagenolytic activity of secreted collagenases (Col) is characterized in the B. cereus culture supernatant (csn) and its isolated recombinantly produced ColQ1 is characterized. The data reveals that ColQ1 causes damage on dermal collagen (COL). This results in gaps in the tissue, which might facilitate the spread of bacteria. The importance of B. cereus collagenases is also demonstrated in disease promotion using two inhibitors. Compound 2 shows high efficacy in peptidolytic, gelatinolytic, and COL degradation assays. It also preserves the fibrillar COLs in skin tissue challenged with ColQ1, as well as the viability of skin cells treated with B. cereus csn. A Galleria mellonella model highlights the significance of collagenase inhibition in vivo.

Prerequisite Binding Modes Determine the Dynamics of Action of Covalent Agonists of Ion Channel TRPA1

Transient receptor potential ankyrin 1 (TRPA1) is a transmembrane protein channeling the influx of calcium ions. As a polymodal nocisensor, TRPA1 can be activated by thermal, mechanical stimuli and a wide range of chemically damaging molecules including small volatile environmental toxicants and endogenous algogenic lipids. After activation by such compounds, the ion channel opens up, its central pore widens allowing calcium influx into the cytosol inducing signal transduction pathways. Afterwards, the calcium influx desensitizes irritant evoked responses and results in an inactive state of the ion channel. Recent experimental determination of structures of apo and holo forms of TRPA1 opened the way towards the design of new agonists, which can activate the ion channel. The present study is aimed at the elucidation of binding dynamics of agonists using experimental structures of TRPA1-agonist complexes at the atomic level applying molecular docking and dynamics methods accounting for covalent and non-covalent interactions. Following a test of docking methods focused on the final, holo structures, prerequisite binding modes were detected involving the apo forms. It was shown how reversible interactions with prerequisite binding sites contribute to structural changes of TRPA1 leading to covalent bonding of agonists. The proposed dynamics of action allowed a mechanism-based forecast of new, druggable binding sites of potent agonists.

Harnessing the Role of HDAC6 in Idiopathic Pulmonary Fibrosis: Design, Synthesis, Structural Analysis, and Biological Evaluation of Potent Inhibitors

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease characterized by a progressive-fibrosing phenotype. IPF has been associated with aberrant HDAC activities confirmed by our immunohistochemistry studies on HDAC6 overexpression in IPF lung tissues. We herein developed a series of novel hHDAC6 inhibitors, having low inhibitory potency over hHDAC1 and hHDAC8, as potential pharmacological tools for IPF treatment. Their inhibitory potency was combined with low in vitro and in vivo toxicity. Structural analysis of 6h and structure-activity relationship studies contributed to the optimization of the binding mode of the new molecules. The best-performing analogues were tested for their efficacy in inhibiting fibrotic sphere formation and cell viability, proving their capability in reverting the IPF phenotype. The efficacy of analogue 6h was also determined in a validated human lung model of TGF-β1-dependent fibrogenesis. The results highlighted in this manuscript may pave the way for the identification of first-in-class molecules for the treatment of IPF.

The latest automated docking technologies for novel drug discovery

INTRODUCTION

Molecular docking has been consolidated as one of the most important methods in the molecular modeling field. It has been recognized as a prominent tool in the study of protein-ligand complexes, to describe intermolecular interactions, to accurately predict poses of multiple ligands, to discover novel promising bioactive compounds. Molecular docking methods have evolved in terms of their accuracy and reliability; but there are pending issues to solve for improving the connection between the docking results and the experimental evidence.

AREAS COVERED

In this article, the author reviews very recent innovative molecular docking applications with special emphasis on reverse docking, treatment of protein flexibility, the use of experimental data to guide the selection of docking poses, the application of Quantum mechanics(QM) in docking, and covalent docking.

EXPERT OPINION

There are several issues being worked on in recent years that will lead to important breakthroughs in molecular docking methods in the near future These developments are related to more efficient exploration of large datasets and receptor conformations, advances in electronic description, and the use of structural information for guiding the selection of results.

Discovery of 1,2,4-oxadiazole-Containing hydroxamic acid derivatives as histone deacetylase inhibitors potential application in cancer therapy

In this study, a series of novel HDAC inhibitors, using 1,2,4-oxadiazole-containing as the cap group, were synthesized and evaluated in vitro. Compound 14b, N-hydroxy-2-(methyl((3-(1-(4-methylbenzyl)piperidin-4-yl)-1,2,4-oxadiazol-5-yl)methyl)amino)pyrimidine-5-carboxamide, displayed the most potent histone deacetylase (HDAC) inhibition, especially against HDAC1, 2, and 3 with IC₅₀ values of 1.8, 3.6 and 3.0 nM, respectively. In vitro antiproliferative studies confirmed that 14b was more potent than SAHA, with IC₅₀ values against 12 types of cancer cell lines ranging from 9.8 to 44.9 nM. The results of Western blot assays showed that compound 14b can significantly up-regulate the acetylation of the biomarker his-H₃ and molecular docking analyses revealed the mode of action of compound 14b against HDAC1. The results of flow-cytometry analysis suggested that compound 14b induces cell cycle arrest at the G1 phase and has apoptotic effects. Further investigation of the activity of 14b on the primary cells of three patients, showed IC₅₀ values of 21.3, 61.1, and 77.4 nM. More importantly, an oral bioavailability of up to 53.52% was observed for 14b. An in vivo pharmacodynamic evaluation demonstrated that compound 14b can significantly inhibit tumor growth in a Daudi Burkitt's lymphoma xenograft model, with tumor inhibition rates of 53.8 and 46.1% observed at 20 and 10 mg/kg when administered p.o. and i.v., respectively. These results indicate that compound 14b may be a suitable lead for further evaluation and development as an HDAC inhibitor and a potent anticancer agent.

Predicting Positions of Bridging Water Molecules in Nucleic Acid-Ligand Complexes

Over the past two decades, interests in DNA and RNA as drug targets have been growing rapidly. Following the trends observed with protein drug targets, computational approaches for drug design have been developed for this new class of molecules. Our efforts toward the development of a universal docking program, Fitted, led us to focus on nucleic acids. Throughout the development of this docking program, efforts were directed toward displaceable water molecules which must be accurately located for optimal docking-based drug discovery. However, although there is a plethora of methods to place water molecules in and around protein structures, there is, to the best of our knowledge, no such fully automated method for nucleic acids, which are significantly more polar and solvated than proteins. We report herein a new method, Splash'Em (Solvation Potential Laid around Statistical Hydration on Entire Macromolecules) developed to place water molecules within the binding cavity of nucleic acids. This fast method was shown to have high agreement with water positions in crystal structures and will therefore provide essential information to medicinal chemists.

In search of constrained FTY720 and phytosphingosine analogs as dual acting anticancer agents targeting metabolic and epigenetic pathways

A series of compounds containing pyrrolidine and pyrrolizidine cores with appended hydrophobic substituents were prepared as constrained analogs of FTY720 and phytosphingosine. The effect of these compounds on the viability of cancer cells, on downregulation of the nutrient transport systems, and on their ability to cause vacuolation was studied. An attempt to inhibit HDACs with some phosphate esters of our analogs was thwarted by our failure to reproduce the reported inhibitory action of FTY720-phosphate.

Atom Types Independent Molecular Mechanics Method for Predicting the Conformational Energy of Small Molecules

We previously implemented a well-known qualitative chemical principle into an accurate quantitative model computing relative potential energies of conformers. According to this principle, hyperconjugation strength correlates with electronegativity of donors and acceptors. While this earlier version of our model applies to σ bonds, lone pairs, disregarded in this earlier version, also have a major impact on the conformational preferences of molecules. Among the well-established principles used by organic chemists to rationalize some organic chemical behaviors are the anomeric effect, the alpha effect, basicity, and nucleophilicity. These effects are directly related to the presence of lone pairs. We report herein our effort to incorporate lone pairs into our model to extend its applicability domain to any saturated small molecules. The developed model H-TEQ 2 has been validated on a wide variety of molecules from polyaromatic molecules to carbohydrates and molecules with high heteroatoms/carbon ratios. Interestingly, this method, in contrast to common force field-based methods, does not rely on atom types and is virtually applicable to any organic molecules.

Roles of the Active Site Zn(II) and Residues in Substrate Discrimination by Threonyl-tRNA Synthetase: An MD and QM/MM Investigation

Threonyl-tRNA synthetase (ThrRS) is a Zn(II) containing enzyme that catalyzes the activation of threonine and its subsequent transfer to the cognate tRNA. This process is accomplished with remarkable fidelity, with ThrRS being able to discriminate its cognate substrate from similar analogues such as serine and valine. Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) methods have been used to elucidate the role of Zn(II) in the aminoacylation mechanism of ThrRS. More specifically, the role of Zn(II) and active site residues in ThrRS's ability to discriminate between its cognate substrate l-threonine and the noncognate l-serine, l-valine, and d-threonine has been examined. The present results suggest that a role of the Zn(II) ion, with its Lewis acidity, is to facilitate deprotonation of the side chain hydroxyl groups of the aminoacyl moieties of cognate Thr-AMP and noncognate Ser-AMP substrates. In their deprotonated forms, these substrates are able to adopt a conformation preferable for aminoacyl transfer from aa-AMP onto the Ado-3'OH of the tRNA^Thr cosubstrate. Relative to the neutral substrates, when the substrates are first deprotonated with the assistance of the Zn(II) ion, the barrier for the rate-limiting step is decreased significantly by 42.0 and 39.2 kJ mol^-1 for l-Thr-AMP and l-Ser-AMP, respectively. An active site arginyl also plays a key role in stabilizing the buildup of negative charge on the substrate's bridging phosphate oxygen during the mechanism. For the enantiomeric substrate analogue d-Thr-AMP, product formation is highly disfavored, and as a result, the reverse reaction has a very low barrier of 16.0 kJ mol^-1.

Inhibitors of Influenza Virus Polymerase Acidic (PA) Endonuclease: Contemporary Developments and Perspectives

Influenza virus (IFV) causes periodic global influenza pandemics, resulting in substantial socioeconomic loss and burden on medical facilities. Yearly variation in the effectiveness of vaccines, slow responsiveness to vaccination in cases of pandemic IFV, and emerging resistance to available drugs highlight the need to develop additional small-molecular inhibitors that act on IFV proteins. One promising target is polymerase acidic (PA) endonuclease, which is a bridged dinuclear metalloenzyme that plays a crucial role in initiating IFV replication. During the past decade, intensive efforts have been made to develop small-molecular inhibitors of this endonuclease as candidate agents for treatment of IFV infection. Here, we review the current status of development of PA endonuclease inhibitors and we discuss the applicability of newer medicinal-chemistry strategies for the discovery more potent, selective, and safer inhibitors.

Captopril/enalapril inhibit promiscuous esterase activity of carbonic anhydrase at micromolar concentrations: An in vitro study

The inhibitory activity of captopril, a thiol-containing competitive inhibitor of the angiotensin-converting enzyme, ACE, against esterase activity of carbonic anhydrase, CA was investigated. This small molecule, as well as enalapril, was selected in order to represents both thiol and carboxylate, as two well-known metal binding functional groups of metalloprotein inhibitors. Since captopril, has also been observed to inhibit other metalloenzymes such as tyrosinase and metallo-beta lactamase through binding to the catalytic metal ions and regarding CA as a zinc-containing metallo-enzyme, in the current study, we set out to determine whether captopril/enalapril inhibit CA esterase activity of the purified human CA II or not? Then, we revealed the inhibitors' potencies (IC₅₀, K_i and K_diss values) and also mode of inhibition. Our results also showed that enalapril is more potent CA inhibitor than captopril. Since enalapril represents no sulfhydryl moiety, thus carboxylate group may have a determinant role in inhibiting of CA esterase activity, the conclusion confirmed by molecular docking studies. Additionally, since CA inhibitory potencies of captopril/enalapril were much lower than those of classic sulfonamide drugs, the findings of the current study may explain why these drugs exhibit no effective CA inhibition at the concentrations reached in vivo and also may shed light on the way of generating new class of inhibitors that will discriminately inhibit various CA isoforms.

Medicinal Chemistry Projects Requiring Imaginative Structure-Based Drug Design Methods

Computational methods for docking small molecules to proteins are prominent in drug discovery. There are hundreds, if not thousands, of documented examples-and several pertinent cases within our research program. Fifteen years ago, our first docking-guided drug design project yielded nanomolar metalloproteinase inhibitors and illustrated the potential of structure-based drug design. Subsequent applications of docking programs to the design of integrin antagonists, BACE-1 inhibitors, and aminoglycosides binding to bacterial RNA demonstrated that available docking programs needed significant improvement. At that time, docking programs primarily considered flexible ligands and rigid proteins. We demonstrated that accounting for protein flexibility, employing displaceable water molecules, and using ligand-based pharmacophores improved the docking accuracy of existing methods-enabling the design of bioactive molecules. The success prompted the development of our own program, Fitted, implementing all of these aspects. The primary motivation has always been to respond to the needs of drug design studies; the majority of the concepts behind the evolution of Fitted are rooted in medicinal chemistry projects and collaborations. Several examples follow: (1) Searching for HDAC inhibitors led us to develop methods considering drug-zinc coordination and its effect on the pKa of surrounding residues. (2) Targeting covalent prolyl oligopeptidase (POP) inhibitors prompted an update to Fitted to identify reactive groups and form bonds with a given residue (e.g., a catalytic residue) when the geometry allows it. Fitted-the first fully automated covalent docking program-was successfully applied to the discovery of four new classes of covalent POP inhibitors. As a result, efficient stereoselective syntheses of a few screening hits were prioritized rather than synthesizing large chemical libraries-yielding nanomolar inhibitors. (3) In order to study the metabolism of POP inhibitors by cytochrome P450 enzymes (CYPs)-for toxicology studies-the program Impacts was derived from Fitted and helped us to reveal a complex metabolism with unforeseen stereocenter isomerizations. These efforts, combined with those of other docking software developers, have strengthened our understanding of the complex drug-protein binding process while providing the medicinal chemistry community with useful tools that have led to drug discoveries. In this Account, we describe our contributions over the past 15 years-within their historical context-to the design of drug candidates, including BACE-1 inhibitors, POP covalent inhibitors, G-quadruplex binders, and aminoglycosides binding to nucleic acids. We also remark the necessary developments of docking programs, specifically Fitted, that enabled structure-based design to flourish and yielded multiple fruitful, rational medicinal chemistry campaigns.

Structure of 'linkerless' hydroxamic acid inhibitor-HDAC8 complex confirms the formation of an isoform-specific subpocket

Histone deacetylases (HDACs) catalyze the hydrolysis of acetylated lysine side chains in histone and non-histone proteins, and play a critical role in the regulation of many biological processes, including cell differentiation, proliferation, senescence, and apoptosis. Aberrant HDAC activity is associated with cancer, making these enzymes important targets for drug design. In general, HDAC inhibitors (HDACi) block the proliferation of tumor cells by inducing cell differentiation, cell cycle arrest, and/or apoptosis, and comprise some of the leading therapies in cancer treatments. To date, four HDACi have been FDA approved for the treatment of cancers: suberoylanilide hydroxamic acid (SAHA, Vorinostat, Zolinza®), romidepsin (FK228, Istodax®), belinostat (Beleodaq®), and panobinostat (Farydak®). Most current inhibitors are pan-HDACi, and non-selectively target a number of HDAC isoforms. Six previously reported HDACi were rationally designed, however, to target a unique sub-pocket found only in HDAC8. While these inhibitors were indeed potent against HDAC8, and even demonstrated specificity for HDAC8 over HDACs 1 and 6, there were no structural data to confirm the mode of binding. Here we report the X-ray crystal structure of Compound 6 complexed with HDAC8 to 1.98Å resolution. We also describe the use of molecular docking studies to explore the binding interactions of the other 5 related HDACi. Our studies confirm that the HDACi induce the formation of and bind in the HDAC8-specific subpocket, offering insights into isoform-specific inhibition.

Inside HDACs with more selective HDAC inhibitors

Inhibitors of histone deacetylases (HDACs) are nowadays part of the therapeutic arsenal mainly against cancers, with four compounds approved by the Food and Drug Administration. During the last five years, several groups have made continuous efforts to improve this class of compounds, designing more selective compounds or compounds with multiple capacities. After a survey of the HDAC biology and structures, this review summarizes the results of the chemists working in this field, and highlights when possible the behavior of the molecules inside their targets.

Development of Purine-Based Hydroxamic Acid Derivatives: Potent Histone Deacetylase Inhibitors with Marked in Vitro and in Vivo Antitumor Activities

In the present study, a series of novel histone deacetylase (HDAC) inhibitors using the morpholinopurine as the capping group were designed and synthesized. Several compounds demonstrated significant HDAC inhibitory activities and antiproliferative effects against diverse human tumor cell lines. Among them, compound 10o was identified as a potent class I and class IIb HDAC inhibitor with good pharmaceutical profile and druglike properties. Western blot analysis further confirmed that 10o more effectively increased acetylated histone H3 than panobinostat (LBH-589) and vorinostat (SAHA) at the same concentration in vitro. In in vivo efficacy evaluations of HCT116, MV4-11, Ramos, and MM1S xenograft models, 10o showed higher efficacy than SAHA or LBH-589 without causing significant loss of body weight and toxicity. All the results indicated that 10o could be a suitable candidate for treatment of both solid and hematological cancer.

Docking Screens for Novel Ligands Conferring New Biology

It is now plausible to dock libraries of 10 million molecules against targets over several days or weeks. When the molecules screened are commercially available, they may be rapidly tested to find new leads. Although docking retains important liabilities (it cannot calculate affinities accurately nor even reliably rank order high-scoring molecules), it can often can distinguish likely from unlikely ligands, often with hit rates above 10%. Here we summarize the improvements in libraries, target quality, and methods that have supported these advances, and the open access resources that make docking accessible. Recent docking screens for new ligands are sketched, as are the binding, crystallographic, and in vivo assays that support them. Like any technique, controls are crucial, and key experimental ones are reviewed. With such controls, docking campaigns can find ligands with new chemotypes, often revealing the new biology that may be docking's greatest impact over the next few years.

Efficiently functionalized oxacalix[4]arenes: Synthesis, characterization and exploration of their biological profile as novel HDAC inhibitors

A series of novel substituted oxacalix[4]arene has been synthesized and explored for their biological profile by evaluating anticancer, antifungal and antibacterial properties. The derivatives have been characterized by various spectroscopic techniques such as IR, (1)H NMR, (13)C NMR and Mass spectrometry. Many compounds showed strong inhibition (MIC) in the range of ∼0-50 μM with interesting cytotoxic activities against Hela cells in particular. The compounds were theoretically evaluated by docking studies as potential histone deacetylase inhibitors (HDACi). The study indicates that compounds bound adequately with HDAC, and hence complemented the experimental findings.

Synthesis, biological characterization and molecular modeling insights of spirochromanes as potent HDAC inhibitors

In the last decades, inhibitors of histone deacetylases (HDAC) have become an important class of anti-cancer agents. In a previous study we described the synthesis of spiro[chromane-2,4'-piperidine]hydroxamic acid derivatives able to inhibit histone deacetylase enzymes. Herein, we present our exploration for new derivatives by replacing the piperidine moiety with various cycloamines. The goal was to obtain highly potent compounds with a good in vitro ADME profile. In addition, molecular modeling studies unravelled the binding mode of these inhibitors.

Activity prediction of substrates in NADH-dependent carbonyl reductase by docking requires catalytic constraints and charge parameterization of catalytic zinc environment

Molecular docking of substrates is more challenging compared to inhibitors as the reaction mechanism has to be considered. This becomes more pronounced for zinc-dependent enzymes since the coordination state of the catalytic zinc ion is of greater importance. In order to develop a predictive substrate docking protocol, we have performed molecular docking studies of diketone substrates using the catalytic state of carbonyl reductase 2 from Candida parapsilosis (CPCR2). Different docking protocols using two docking methods (AutoDock Vina and AutoDock4.2) with two different sets of atomic charges (AM1-BCC and HF-RESP) for catalytic zinc environment and substrates as well as two sets of vdW parameters for zinc ion were examined. We have selected the catalytic binding pose of each substrate by applying mechanism based distance criteria. To compare the performance of the docking protocols, the correlation plots for the binding energies of these catalytic poses were obtained against experimental Vmax values of the 11 diketone substrates for CPCR2. The best correlation of 0.73 was achieved with AutoDock4.2 while treating catalytic zinc ion in optimized non-bonded (NBopt) state with +1.01 charge on the zinc ion, compared to 0.36 in non-bonded (+2.00 charge on the zinc ion) state. These results indicate the importance of catalytic constraints and charge parameterization of catalytic zinc environment for the prediction of substrate activity in zinc-dependent enzymes by molecular docking. The developed predictive docking protocol described here is in principle generally applicable for the efficient in silico substrate spectra characterization of zinc-dependent ADH.

Role of computer-aided drug design in modern drug discovery

Drug discovery utilizes chemical biology and computational drug design approaches for the efficient identification and optimization of lead compounds. Chemical biology is mostly involved in the elucidation of the biological function of a target and the mechanism of action of a chemical modulator. On the other hand, computer-aided drug design makes use of the structural knowledge of either the target (structure-based) or known ligands with bioactivity (ligand-based) to facilitate the determination of promising candidate drugs. Various virtual screening techniques are now being used by both pharmaceutical companies and academic research groups to reduce the cost and time required for the discovery of a potent drug. Despite the rapid advances in these methods, continuous improvements are critical for future drug discovery tools. Advantages presented by structure-based and ligand-based drug design suggest that their complementary use, as well as their integration with experimental routines, has a powerful impact on rational drug design. In this article, we give an overview of the current computational drug design and their application in integrated rational drug development to aid in the progress of drug discovery research.

Computer-aided Molecular Design of Compounds Targeting Histone Modifying Enzymes

Growing evidences show that epigenetic mechanisms play crucial roles in the genesis and progression of many physiopathological processes. As a result, research in epigenetic grew at a fast pace in the last decade. In particular, the study of histone post-translational modifications encountered an extraordinary progression and many modifications have been characterized and associated to fundamental biological processes and pathological conditions. Histone modifications are the catalytic result of a large set of enzyme families that operate covalent modifications on specific residues at the histone tails. Taken together, these modifications elicit a complex and concerted processing that greatly contribute to the chromatin remodeling and may drive different pathological conditions, especially cancer. For this reason, several epigenetic targets are currently under validation for drug discovery purposes and different academic and industrial programs have been already launched to produce the first pre-clinical and clinical outcomes. In this scenario, computer-aided molecular design techniques are offering important tools, mainly as a consequence of the increasing structural information available for these targets. In this mini-review we will briefly discuss the most common types of known histone modifications and the corresponding operating enzymes by emphasizing the computer-aided molecular design approaches that can be of use to speed-up the efforts to generate new pharmaceutically relevant compounds.

AutoDock4(Zn): an improved AutoDock force field for small-molecule docking to zinc metalloproteins

Zinc is present in a wide variety of proteins and is important in the metabolism of most organisms. Zinc metalloenzymes are therapeutically relevant targets in diseases such as cancer, heart disease, bacterial infection, and Alzheimer's disease. In most cases a drug molecule targeting such enzymes establishes an interaction that coordinates with the zinc ion. Thus, accurate prediction of the interaction of ligands with zinc is an important aspect of computational docking and virtual screening against zinc containing proteins. We have extended the AutoDock force field to include a specialized potential describing the interactions of zinc-coordinating ligands. This potential describes both the energetic and geometric components of the interaction. The new force field, named AutoDock4Zn, was calibrated on a data set of 292 crystal complexes containing zinc. Redocking experiments show that the force field provides significant improvement in performance in both free energy of binding estimation as well as in root-mean-square deviation from the crystal structure pose. The new force field has been implemented in AutoDock without modification to the source code.

Exploring the influence of the protein environment on metal-binding pharmacophores

The binding of a series of metal-binding pharmacophores (MBPs) related to the ligand 1-hydroxypyridine-2-(1H)-thione (1,2-HOPTO) in the active site of human carbonic anhydrase II (hCAII) has been investigated. The presence and/or position of a single methyl substituent drastically alters inhibitor potency and can result in coordination modes not observed in small-molecule model complexes. It is shown that this unexpected binding mode is the result of a steric clash between the methyl group and a highly ordered water network in the active site that is further stabilized by the formation of a hydrogen bond and favorable hydrophobic contacts. The affinity of MBPs is dependent on a large number of factors including donor atom identity, orientation, electrostatics, and van der Waals interactions. These results suggest that metal coordination by metalloenzyme inhibitors is a malleable interaction and that it is thus more appropriate to consider the metal-binding motif of these inhibitors as a pharmacophore rather than a "chelator". The rational design of inhibitors targeting metalloenzymes will benefit greatly from a deeper understanding of the interplay between the variety of forces governing the binding of MBPs to active site metal ions.