A comparative docking study and the design of potentially selective MMP inhibitors

In-silico strategies for identification of potent inhibitor for MMP-1 to prevent metastasis of breast cancer

Matrix Metalloproteinase-1 (MMP-1) has been often upregulated in advanced breast cancers, known to participate in ECM degradation, migration, invasion, thus leading to metastasis. Due to these effects, the condition is often reported to inversely correlate with survival in advanced breast cancers. In the present study, in-silico method was adopted based on selective non zinc binding inhibitors of MMP-1. ADME properties were predicted for PASS filtered compounds and docking calculations were performed using Glide XP and IFD protocols of Schrodinger program. We identified six ligands as potent inhibitors and validated by observing structures and the interactions of MMP-1. The identified hits were validated using molecular dynamics simulation studies. Electronic structure analysis was performed for two top hit compounds myricetin and quercetin using density function theory (DFT) at B3LYP/6-31**G level to understand their molecular reactivity. Finally, one compound myricetin has emerged as the structurally stable compound with -7.801 kcal/mol and reasonable pose inside the binding site. Molecular dynamics results indicated that myricetin forms a stable interaction with the key amino acid residues such as Glu209, Glu219, Tyr240 and Pro238. In addition, it did not form any binding with the catalytic zinc at its active site. The interaction pattern of myricetin at its substrate binding site exhibited to be potent MMP-1 inhibitor. DFT study also showed that it has more potent inhibitory effect and solubility. These factors altogether show that myricetin could be considered as the best among the compounds evaluated in inhibiting MMP-1 thereby preventing metastasis of breast cancer. Communicated by Ramaswamy H. Sarma.

Design of Potent Panobinostat Histone Deacetylase Inhibitor Derivatives: Molecular Considerations for Enhanced Isozyme Selectivity between HDAC2 and HDAC8

Histone Deacetylases (HDACs) are an important family of 18 isozymes, which are being pursued as drug targets for many types of disorders. HDAC2 and HDAC8 are two of the isozymes, which have been identified as drug targets for the design of anti-cancer, neurodegenerative, immunological, and anti-parasitic agents. Design of potent HDAC2 and HDAC8 inhibitors will be useful for the therapeutic advances in many disorders. This work was undertaken to develop potent HDAC2 and HDAC8 inhibitors. A docking study was performed comparing panobinostat derivatives in both HDAC2 and HDAC8. Six of our derivatives showed stronger binding to HDAC2 than panobinostat, and two of our derivatives showed stronger binding to HDAC8 than panobinostat. We evaluated the molecular features, which improved potency of our inhibitors over panobinostat and also identified another molecular consideration, which could be used to enhance histone deacetylase inhibitor (HDACi) selectivity towards either the HDAC2 or HDAC8 isozymes. The results of this work can be used to assist future design of more potent and selective HDACi for HDAC2 and HDAC8.

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.

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

Metalloenzymes are ubiquitous proteins which feature one or more metal ions either directly involved in the enzymatic activity and/or structural properties (i.e., zinc fingers). Several members of this class take advantage of the Lewis acidic properties of zinc ions to carry out their various catalytic transformations including isomerization or amide cleavage. These enzymes have been validated as drug targets for a number of diseases including cancer; however, despite their pharmaceutical relevance and the availability of crystal structures, structure-based drug design methods have been poorly and indirectly parametrized for these classes of enzymes. More specifically, the metal coordination component and proton transfers of the process of drugs binding to metalloenzymes have been inadequately modeled by current docking programs, if at all. In addition, several known issues, such as coordination geometry, atomic charge variability, and a potential proton transfer from small molecules to a neighboring basic residue, have often been ignored. We report herein the development of specific functions and parameters to account for zinc-drug coordination focusing on the above-listed phenomena and their impact on docking to zinc metalloenzymes. These atom-type-dependent but atomic charge-independent functions implemented into Fitted 3.1 enable the simulation of drug binding to metalloenzymes, considering an acid-base reaction with a neighboring residue when necessary with good accuracy.

Anacardic acid inhibits the catalytic activity of matrix metalloproteinase-2 and matrix metalloproteinase-9

Cashew nut shell liquid (CNSL) has been used in traditional medicine for the treatment of a wide variety of pathophysiological conditions. To further define the mechanism of CNSL action, we investigated the effect of cashew nut shell extract (CNSE) on two matrix metalloproteinases, MMP-2/gelatinase A and MMP-9/gelatinase B, which are known to have critical roles in several disease states. We observed that the major constituent of CNSE, anacardic acid, markedly inhibited the gelatinase activity of 3T3-L1 cells. Our gelatin zymography studies on these two secreted gelatinases, present in the conditioned media from 3T3-L1 cells, established that anacardic acid directly inhibited the catalytic activities of both MMP-2 and MMP-9. Our docking studies suggested that anacardic acid binds into the MMP-2/9 active site, with the carboxylate group of anacardic acid chelating the catalytic zinc ion and forming a hydrogen bond to a key catalytic glutamate side chain and the C15 aliphatic group being accommodated within the relatively large S1' pocket of these gelatinases. In agreement with the docking results, our fluorescence-based studies on the recombinant MMP-2 catalytic core domain demonstrated that anacardic acid directly inhibits substrate peptide cleavage in a dose-dependent manner, with an IC₅₀ of 11.11 μM. In addition, our gelatinase zymography and fluorescence data confirmed that the cardol-cardanol mixture, salicylic acid, and aspirin, all of which lack key functional groups present in anacardic acid, are much weaker MMP-2/MMP-9 inhibitors. Our results provide the first evidence for inhibition of gelatinase catalytic activity by anacardic acid, providing a novel template for drug discovery and a molecular mechanism potentially involved in CNSL therapeutic action.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

Docking and quantitative structure-activity relationship studies for 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline, 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline, and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine derivatives as c-Met kinase inhibitors

We have performed docking of 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline (FPTA), 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline (FPPA), and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine (AFPP) derivatives complexed with c-Met kinase to study the orientations and preferred active conformations of these inhibitors. The study was conducted on a selected set of 103 compounds with variations both in structure and activity. Docking helped to analyze the molecular features which contribute to a high inhibitory activity for the studied compounds. In addition, the predicted biological activities of the c-Met kinase inhibitors, measured as IC(50) values were obtained by using quantitative structure-activity relationship (QSAR) methods: Comparative molecular similarity analysis (CoMSIA) and multiple linear regression (MLR) with topological vectors. The best CoMSIA model included steric, electrostatic, hydrophobic, and hydrogen bond-donor fields; furthermore, we found a predictive model containing 2D-autocorrelation descriptors, GETAWAY descriptors (GETAWAY: Geometry, Topology and Atom-Weight AssemblY), fragment-based polar surface area (PSA), and MlogP. The statistical parameters: cross-validate correlation coefficient and the fitted correlation coefficient, validated the quality of the obtained predictive models for 76 compounds. Additionally, these models predicted adequately 25 compounds that were not included in the training set.

Docking of sialic acid analogues against influenza A hemagglutinin: a correlational study between experimentally measured and computationally estimated affinities

A molecular docking tool of AutoDock3.05 was evaluated for its ability to reproduce experimentally determined affinities of various sialic acid analogues toward hemagglutinin of influenza A virus. With the exception of those with a C6-modified glycerol side chain, the experimental binding affinities of most sialic acid analogues (C2, C4 and C5-substituted) determined by viral hemadsorption inhibition assay, hemagglutination inhibition assay and nuclear magnetic resonance correlated well with the computationally estimated free energy of binding. Sialic acid analogues with modified glycerol side chains showed only poor correlation between the experimentally determined hemagglutinin inhibitor affinities and AutoDock3.05 scores, suggesting high mobility of the glutamic acid side chain at the glycerol binding pocket, which is difficult to simulate using a flexi-rigid molecular docking approach. In conclusion, except for some glycerol-substituted sialic acid analogues, the results showed the effectiveness of AutoDock3.05 searching and scoring functions in estimating affinities of sialic acid analogues toward influenza A hemagglutinin, making it a reliable tool for screening a database of virtually designed sialic acid analogues for hemagglutinin inhibitors.

Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go

Accelerating the drug discovery process requires predictive computational protocols capable of reducing or simplifying the synthetic and/or combinatorial challenge. Docking-based virtual screening methods have been developed and successfully applied to a number of pharmaceutical targets. In this review, we first present the current status of docking and scoring methods, with exhaustive lists of these. We next discuss reported comparative studies, outlining criteria for their interpretation. In the final section, we describe some of the remaining developments that would potentially lead to a universally applicable docking/scoring method.

Crystal Structures of MMP-9 Complexes with Five Inhibitors: Contribution of the Flexible Arg424 Side-chain to Selectivity

Human matrix metalloproteinase 9 (MMP-9), also called gelatinase B, is particularly involved in inflammatory processes, bone remodelling and wound healing, but is also implicated in pathological processes such as rheumatoid arthritis, atherosclerosis, tumour growth, and metastasis. We have prepared the inactive E402Q mutant of the truncated catalytic domain of human MMP-9 and co-crystallized it with active site-directed synthetic inhibitors of different binding types. Here, we present the X-ray structures of five MMP-9 complexes with gelatinase-specific, tight binding inhibitors: a phosphinic acid (AM-409), a pyrimidine-2,4,6-trione (RO-206-0222), two carboxylate (An-1 and MJ-24), and a trifluoromethyl hydroxamic acid inhibitor (MS-560). These compounds bind by making a compromise between optimal coordination of the catalytic zinc, favourable hydrogen bond formation in the active-site cleft, and accommodation of their large hydrophobic P1' groups in the slightly flexible S1' cavity, which exhibits distinct rotational conformations of the Pro421 carbonyl group in each complex. In all these structures, the side-chain of Arg424 located at the bottom of the S1' cavity is not defined in the electron density beyond C(gamma), indicating its mobility. However, we suggest that the mobile Arg424 side-chain partially blocks the S1' cavity, which might explain the weaker binding of most inhibitors with a long P1' side-chain for MMP-9 compared with the closely related MMP-2 (gelatinase A), which exhibits a short threonine side-chain at the equivalent position. These novel structural details should facilitate the design of more selective MMP-9 inhibitors.

QSAR modeling of matrix metalloproteinase inhibition by N-hydroxy-alpha-phenylsulfonylacetamide derivatives

The main molecular features which determine the selectivity of a set of 80 N-hydroxy-alpha-phenylsulfonylacetamide derivatives (HPSAs) in the inhibition of three matrix metalloproteinases (MMP-1, MMP-9, and MMP-13) have been identified by using linear and nonlinear predictive models. The molecular information has been encoded in 2D autocorrelation descriptors, obtained from different weighting schemes. The linear models were built by multiple linear regression (MLR) combined with genetic algorithm (GA), and a robust QSAR mapping paradigm. The Bayesian-regularized genetic neural network (BRGNN) was employed for nonlinear modeling. In such approaches each model could have its own set of input variables. All models were predictive according to internal and external validation experiments; but the best results correspond to nonlinear ones. The 2D autocorrelation space brings different descriptors for each MMP inhibition, and suggests the atomic properties relevant for the inhibitors to interact with each MMP active site. On the basis of the current results, the reported models have the potential to discover new potent and selective inhibitors and bring useful molecular information about the ligand specificity for MMP S(1)(') and S(2)(') subsites.

Evaluation of docking programs for predicting binding of Golgi α-mannosidase II inhibitors: A comparison with crystallography

Golgi alpha-mannosidase II (GMII), a zinc-dependent glycosyl hydrolase, is a promising target for drug development in anti-tumor therapies. Using X-ray crystallography, we have determined the structure of Drosophila melanogaster GMII (dGMII) complexed with three different inhibitors exhibiting IC50's ranging from 80 to 1000 microM. These structures, along with those of seven other available dGMII/inhibitor complexes, were then used as a basis for the evaluation of seven docking programs (GOLD, Glide, FlexX, AutoDock, eHiTS, LigandFit, and FITTED). We found that small inhibitors could be accurately docked by most of the software, while docking of larger compounds (i.e., those with extended aromatic cycles or long aliphatic chains) was more problematic. Overall, Glide provided the best docking results, with the most accurately predicted binding around the active site zinc atom. Further evaluation of Glide's performance revealed its ability to extract active compounds from a benchmark library of decoys.

Linear and nonlinear QSAR study of N-hydroxy-2-[(phenylsulfonyl)amino]acetamide derivatives as matrix metalloproteinase inhibitors

The inhibitory activity (IC50) toward matrix metalloproteinases (MMP-1, MMP-2, MMP-3, MMP-9, and MMP-13) of N-hydroxy-2-[(phenylsulfonyl)amino]acetamide derivatives (HPSAAs) has been successfully modeled using 2D autocorrelation descriptors. The relevant molecular descriptors were selected by linear and nonlinear genetic algorithm (GA) feature selection using multiple linear regression (MLR) and Bayesian-regularized neural network (BRANN) approaches, respectively. The quality of the models was evaluated by means of cross-validation experiments and the best results correspond to nonlinear ones (Q2>0.7 for all models). Despite the high correlation between the studied compound IC50 values, the 2D autocorrelation space brings different descriptors for each MMP inhibition. On the basis of these results, these models contain useful molecular information about the ligand specificity for MMP S'1, S1, and S'2 pockets.

Amber force field implementation, molecular modelling study, synthesis and MMP-1/MMP-2 inhibition profile of (R)- and (S)-N-hydroxy-2-(N-isopropoxybiphenyl-4-ylsulfonamido)-3-methylbutanamides

Ab initio calculations (B3LYP/Lanl2DZ level of theory) were performed in this study to determine all the structural and catalytic zinc parameters required in order to study MMPs and their complexes with hydroxamate inhibitors by means of the AMBER force field. The parameters thus obtained were used in order to study the docking of some known MMPi (Batimastat, CGS 27023A and Prinomastat) and our previously described inhibitor a which had shown an inhibitory activity for MMP-1, and -2, with the aim of explaining the different selectivity. On this basis the two enantiomers (R)-b and (S)-b were designed and synthesized, as more potent MMP-2 inhibitors than our previously described inhibitor a. Between these two enantiomers the eutomer (R)-b proved to be 24.7 times and 15.3 times more potent than CGS 27023A and the parent compound a on MMP-2, maintaining a higher index of MMP-2/MMP-1 selectivity compared with CGS 27023A and the more potent inhibitor Prinomastat. The hydroxamate (R)-b can be considered as a progenitor of a new class of biphenylsulfonamido-based inhibitors that differ from compound a in the presence of an alkyl side chain on the C alpha atom, and show different potency and selectivity profiles on the two MMPs considered.

Combining pharmacophore search, automated docking, and molecular dynamics simulations as a novel strategy for flexible docking. Proof of concept: docking of arginine-glycine-aspartic acid-like compounds into the alphavbeta3 binding site

A novel and highly efficient flexible docking approach is presented where the conformations (internal degrees of freedom) and orientations (external degrees of freedom) of the ligands are successively considered. This hybrid method takes advantage of the synergistic effects of structure-based and ligand-based drug design techniques. Preliminary antagonist-derived pharmacophore determination provides the postulated bioactive conformation. Subsequent docking of this pharmacophore to the receptor crystal structure results in a postulated pharmacophore/receptor binding mode. Pharmacophore-oriented docking of antagonists is subsequently achieved by matching ligand interacting groups with pharmacophore points. Molecular dynamics in water refines the proposed complexes. To validate the method, arginine-glycine-aspartic acid (RGD) containing peptides, pseudopeptides, and RGD-like antagonists were docked to the crystal structure of alphavbeta3 holoprotein and apoprotein. The proposed directed docking was found to be more accurate, faster, and less biased with respect to the protein structure (holo and apoprotein) than DOCK, Autodock, and FlexX docking methods. The successful docking of an antagonist recently cocrystallized with the receptor to both apo and holoprotein is particularly appealing. The results summarized in this report illustrated the efficiency of our light CoMFA/rigid body docking hybrid method.

New N-arylsulfonyl-N-alkoxyaminoacetohydroxamic acids as selective inhibitors of gelatinase A (MMP-2)

New N-arylsulfonyl-substituted alkoxyaminoaceto hydroxamic acid derivatives of types 8 and 10 designed as oxa-analogues of known sulfonamide-based MMPi of types 2 and 7 were synthesized and tested for their inhibitory activities on some matrix metalloproteinases. The combination of a biphenylsulfonamide group with oxyamino oxygen in the pharmacophoric central skeleton of sulfonamide-based MMPi obtained in the new sulfonamides 10 seems to be able to give selectivity for MMP-2 over MMP-1. The most potent derivative of this type, 10a, shows similar anti-invasive properties to the analogue reference drug CGS27023A, 2, in an in vitro model of invasion on matrigel, carried out on cellular lines of fibrosarcoma HT1080 (tumoural cells over-expressing MMP-2 and MMP-9).

A practical approach to docking of zinc metalloproteinase inhibitors

Forty zinc-dependent metalloproteinase/ligand complexes with known crystal structures were re-docked using five docking/scoring approaches (DOCK, FlexX, DrugScore, GOLD, and AutoDock). Correct geometry of the coordination bonds between the ligand's zinc binding group (ZBG) and the catalytic zinc is important for docking accuracy and scoring reliability. More than 75% of docked poses with RMSD less than 2A were found to have appropriate ZBG binding, but for poor ZBG binding, about 95% of poses failed to dock correctly. Elimination of poses with inappropriate zinc binding resulted in better binding energy predictions that were further improved by dividing the ligands into subsets according to the ZBG (carboxylates, hydroxamates, and phosphorus containing groups). After a subset re-scoring using the regression functions obtained for individual subsets, DrugScore was able to explain 77% and the consensus scoring scheme X-CSCORE even 88% of variance in binding energies. The approach combining ZBG-based pose selection and subset re-scoring improved the hit rate in virtual screening for metalloproteinase inhibitors for all tested methods by 4-16%.

The inhibitor thiomandelic acid binds to both metal ions in metallo-beta-lactamase and induces positive cooperativity in metal binding

Thiomandelic acid is a simple, broad spectrum, and reasonably potent inhibitor of metallo-beta-lactamases, enzymes that mediate resistance to beta-lactam antibiotics. We report studies by NMR and perturbed angular correlation (PAC) spectroscopy of the mode of binding of the R and S enantiomers of thiomandelic acid, focusing on their interaction with the two metal ions in cadmium-substituted Bacillus cereus metallo-beta-lactamase. The 113Cd resonances are specifically assigned to the metals in the two individual sites on the protein by using 113Cd-edited 1H NMR spectra. Each enantiomer of thiomandelate produces large downfield shifts of both 113Cd resonances and changes in the PAC spectra, which indicate that they bind such that the thiol of the inhibitor bridges between the two metals. For R-thiomandelate, this is unambiguously confirmed by the observation of scalar coupling between Halpha of the inhibitor and both cadmium ions. The NMR and PAC spectra reveal that the two chiral forms of the inhibitor differ in the details of their coordination geometry. The complex with R-thiomandelate, but not that with the S-enantiomer, shows evidence in the PAC spectra of a dynamic process in the nanosecond time regime, the possible nature of which is discussed. The thiomandelate complex of the mononuclear enzyme can be detected only at low metal to enzyme stoichiometry; the relative populations of mononuclear and binuclear enzyme as a function of cadmium concentration provide clear evidence for positive cooperativity in metal ion binding in the presence of the inhibitor, in contrast to the negative cooperativity observed in the free enzyme.