Serdar Durdagi

Kinetic and docking studies of phenol-based inhibitors of carbonic anhydrase isoforms I, II, IX and XII evidence a new binding mode within the enzyme active site

Carbonic anhydrases (CAs, EC 4.2.1.1) are inhibited by sulfonamides, inorganic anions, phenols, coumarins (acting as prodrugs) and polyamines. A novel class of CA inhibitors (CAIs), interacting with the CA isozymes I, II (cytosolic) and IX, XII (transmembrane, tumor-associated) in a different manner, is reported here. Kinetic measurements allowed us to identify hydroxy-/methoxy-substituted benzoic acids as well as di-/tri-methoxy benzenes as submicromolar-low micromolar inhibitors of the four CA isozymes. Molecular docking studies of a set of such inhibitors within CA I and II allowed us to understand the inhibition mechanism. This new class of inhibitors binds differently compared to all other classes of inhibitors known to date: they were found between the phenol-binding site and the coumarin-binding site, filling thus the middle of the enzyme cavity. They exploit different interactions with amino acid residues and water molecules from the CA active site compared to other classes of inhibitors, offering the possibility to design CAIs with an interesting inhibition profile compared to the clinically used sulfonamides/sulfamates.

Discovery of potent carbonic anhydrase and acetylcholine esterase inhibitors: Novel sulfamoylcarbamates and sulfamides derived from acetophenones

In this study, several novel sulfamides were synthesized and evaluated for their acetylcholine esterase (AChE) and human carbonic anhydrase I, and II isoenzymes (hCA I and II) inhibition profiles. Reductive amination of methoxyacetophenones was used for the synthesis of amines. Amines were converted to sulfamoylcarbamates with chlorosulfonyl isocyanate (CSI) in the presence of BnOH. Pd-C catalyzed hydrogenolysis of sulfamoylcarbamates afforded sulfamides. These novel compounds were good inhibitors of the cytosolic hCA I, and hCA II with Ki values in the range of 45.9±8.9-687.5±84.3 pM for hCA I, and 48.80±8.2-672.2±71.9pM for hCA II. The inhibitory effects of the synthesized novel compounds on AChE were also investigated. The Ki values of these compounds were in the range of 4.52±0.61-38.28±6.84pM for AChE. These results show that hCA I, II, and AChE were effectively inhibited by the novel sulfamoylcarbamates 17-21 and sulfamide derivatives 22-26. All investigated compounds were docked within the active sites of the corresponding enzymes revealing the reasons of the effective inhibitory activity.

Carbonic anhydrase inhibitory properties of novel benzylsulfamides using molecular modeling and experimental studies

In this study, a series of sulfamoyl carbamates and sulfamide derivatives were synthesized. Six commercially available benzyl amines and BnOH were reacted with chlorosulfonyl isocyanate (CSI) to give sulfamoyl carbamates. Pd-C catalyzed hydrogenolysis reactions of carbamates afforded sulfamides. The inhibition effects of novel benzylsulfamides on the carbonic anhydrase I, and II isoenzymes (CA I, and CA II) purified from fresh human blood red cells were determined by Sepharose-4B-L-Tyrosine-sulfanilamide affinity chromatography. In vitro studies were shown that all of novel synthesized benzylsulfamide analogs inhibited, concentration dependently, both hCA isoenzyme activities. The novel benzylsulfamide compounds investigated here exhibited nanomolar inhibition constants against the two isoenzymes. Ki values were in the range of 28.48±0.01-837.09±0.19nM and 112.01±0.01-268.01±0.22nM for hCAI and hCA II isoenzymes, respectively. Molecular modeling approaches were also applied for studied compounds.

Combined Receptor and Ligand-Based Approach to the Universal Pharmacophore Model Development for Studies of Drug Blockade to the hERG1 Pore Domain

Long QT syndrome, LQTS, results in serious cardiovascular disorders, such as tachyarrhythmia and sudden cardiac death. A promiscuous binding of different drugs to the intracavitary binding site in the pore domain (PD) of human ether-a-go-go related gene (hERG) channels leads to a similar dysfunction, known as a drug-induced LQTS. Therefore, an assessment of the blocking ability for potent drugs is of great pragmatic value for molecular pharmacology and medicinal chemistry of hERGs. Thus, we attempted to create an in silico model aimed at blinded drug screening for their blocking ability to the hERG1 PD. Two distinct approaches to the drug blockage, ligand-based QSAR and receptor-based molecular docking methods, are combined for development of a universal pharmacophore model, which provides rapid assessment of drug blocking ability to the hERG1 channel. The best 3D-QSAR model (AAADR.7) from PHASE modeling was selected from a pool consisting of 44 initial candidates. The constructed model using 31 hERG blockers was validated with 9 test set compounds. The resulting model correctly predicted the pIC(50) values of test set compounds as true unknowns. To further evaluate the pharmacophore model, 14 hERG blockers with diverse hERG blocking potencies were selected from literature and they were used as additional external blind test sets. The resulting average deviation between in vitro and predicted pIC(50) values of external test set blockers is found as 0.29 suggesting that the model is able to accuretely predict the pIC(50) values as true unknowns. These pharmacophore models were merged with a previously developed atomistic receptor model for the hERG1 PD and exhibited a high consistency between ligand-based and receptor-based models. Therefore, the developed 3D-QSAR model provides a predictive tool for profiling candidate compounds before their synthesis. This model also indicated the key functional groups determining a high-affinity blockade of the hERG1 channel. To cross-validate consistency between the constructed hERG1 pore domain and the pharmacophore models, we performed docking studies using the homology model of hERG1. To understand how polar or nonpolar moieties of inhibitors stimulate channel inhibition, critical amino acid replacement (i.e., T623, S624, S649, Y652 and F656) at the hERG cavity was examined by in silico mutagenesis. The average docking score differences between wild type and mutated hERG channels was found to have the following order: F656A > Y652A > S624A > T623A > S649A. These results are in agreement with experimental data.

Computational design of novel fullerene analogues as potential HIV-1 PR inhibitors: Analysis of the binding interactions between fullerene inhibitors and HIV-1 PR residues using 3D QSAR, molecular docking and molecular dynamics simulations.

A series of experimentally reported as well as computationally designed monoadducts and bisadducts of [60]fullerene analogues have been used in order to analyze the binding interactions between fullerene based inhibitors and HIV-1 PR employing docking studies. MD simulations of ligand-free and the inhibitor bound HIV-1 PR systems complemented the above studies and provided proper input structure of HIV-1 PR in docking simulations. The obtained results revealed a different orientation of the beta-hairpin flaps at these two systems. In inhibitor bound system, the flaps of the enzyme are pulled in toward the bottom of the active site (the closed form) while, in ligand-free system flaps shifted away from the dual Asp25 catalytic site and this system adopts a semi-open form. The structural analysis of these systems at catalytic and flexible flap regions of the HIV-1 PR through the simulation, assisted in understanding the structural preferences of these regions, as well as, the adopted orientations of fullerene derivatives within the active site of the enzyme. Five different combinations of steroelectronic fields of 3D QSAR/CoMSIA models were obtained from the set of biologically evaluated and computationally designed fullerene derivatives (training set=43, test set=6) in order to predict novel compounds with improved inhibition effect. The best 3D QSAR/CoMSIA model yielded a cross validated r(2) value of 0.739 and a non-cross validated r(2) value of 0.993. The derived model indicated the importance of steric (42.6%), electrostatic (12.7%), H-bond donor (16.7%) and H-bond acceptor (28.0%) contributions. The derived contour plots together with de novo drug design were then used as pilot models for proposing the novel analogues with enhanced binding affinities. Such structures may trigger the interest of medicinal chemists for novel HIV-1 PR inhibitors possessing higher bioactivity.

Modeling of open, closed, and open-inactivated states of the hERG1 channel: structural mechanisms of the state-dependent drug binding.

The human ether-a-go-go related gene 1 (hERG1) K ion channel is a key element for the rapid component of the delayed rectified potassium current in cardiac myocytes. Since there are no crystal structures for hERG channels, creation and validation of its reliable atomistic models have been key targets in molecular cardiology for the past decade. In this study, we developed and vigorously validated models for open, closed, and open-inactivated states of hERG1 using a multistep protocol. The conserved elements were derived using multiple-template homology modeling utilizing available structures for Kv1.2, Kv1.2/2.1 chimera, and KcsA channels. Then missing elements were modeled with the ROSETTA De Novo protein-designing suite and further refined with all-atom molecular dynamics simulations. The final ensemble of models was evaluated for consistency to the reported experimental data from biochemical, biophysical, and electrophysiological studies. The closed state models were cross-validated against available experimental data on toxin footprinting with protein-protein docking using hERG state-selective toxin BeKm-1. Poisson-Boltzmann calculations were performed to determine gating charge and compare it to electrophysiological measurements. The validated structures offered us a unique chance to assess molecular mechanisms of state-dependent drug binding in three different states of the channel.

Nanoscale enzyme inhibitors: fullerenes inhibit carbonic anhydrase by occluding the active site entrance.

We investigated a series of derivatized fullerenes possessing alcohol, amine, and amino acid pendant groups as inhibitors of the zinc enzymes carbonic anhydrases (CAs, EC 4.2.1.1). We discovered that fullerenes bind CAs with submicromolar-low micromolar affinity, despite the fact that these compounds do not possess moieties normally associated with CA inhibitors such as the sulfonamides and their isosteres, or the coumarins. The 13 different mammalian CA isoforms showed a diverse inhibition profile with these compounds. By means of computational methods we assessed the inhibition mechanism as being due to occlusion of the active site entrance by means of the fullerene cage (possessing dimension of the same order of magnitude as the opening of the enzyme cavity, of 1nm). The pendant moieties to the fullerene cage make interactions with amino acid residues from the active site, among which His64, His94, His96, Val121, and Thr200. Fullerenes thus represent a totally new class of nanoscale CA inhibitors which may show applications for targeting physiologically relevant isoforms, such as the dominant CA II and the tumor-associated CA IX.

Acetylcholinesterase and carbonic anhydrase inhibitory properties of novel urea and sulfamide derivatives incorporating dopaminergic 2-aminotetralin scaffolds.

In the present study a series of urea and sulfamide compounds incorporating the tetralin scaffolds were synthesized and evaluated for their acetylcholinesterase (AChE), human carbonic anhydrase (CA, EC 4.2.1.1) isoenzyme I, and II (hCA I and hCA II) inhibitory properties. The urea and their sulfamide analogs were synthesized from the reactions of 2-aminotetralins with N,N-dimethylcarbamoyl chloride and N,N-dimethylsulfamoyl chloride, followed by conversion to the corresponding phenols via O-demethylation with BBr3. The novel urea and sulfamide derivatives were tested for inhibition of hCA I, II and AChE enzymes. These derivatives exhibited excellent inhibitory effects, in the low nanomolar range, with Ki values of 2.61-3.69nM against hCA I, 1.64-2.80nM against hCA II, and in the range of 0.45-1.74nM against AChE. In silico techniques such as, atomistic molecular dynamics (MD) and molecular docking simulations, were used to understand the scenario of the inhibition mechanism upon approaching of the ligands into the active site of the target enzymes. In light of the experimental and computational results, crucial amino acids playing a role in the stabilization of the enzyme-inhibitor adducts were identified.

3D QSAR CoMFA/CoMSIA, molecular docking and molecular dynamics studies of fullerene-based HIV-1 PR inhibitors

For the first time, a set of experimentally reported [60] fullerene derivatives were subjected to the 3D-QSAR/CoMFA and CoMSIA studies. The aim of this study is to propose a series of novel [60] fullerene-based inhibitors with optimal binding affinity for the HIV-1 PR enzyme. The position of the template molecule at the cavity of HIV-1 PR was optimized and 3D QSAR models were developed. Relative contributions of steric/electrostatic fields of the 3D-QSAR/CoMFA and CoMSIA models have shown that steric effects govern the bioactivity of the compounds, but electrostatic interactions play also an important role.The de novo drug design Leapfrog simulations provided a series of novel compounds with predicted improved inhibition effect.

Strategies in the rational drug design.

Rational design is applied in the discovery of novel lead drugs. Its rapid development is mainly attributed to the tremendous advancements in the computer science, statistics, molecular biology, biophysics, biochemistry, medicinal chemistry, pharmacokinetics and pharmacodynamics experienced in the last few decades. The promising feature that characterizes the application of rational drug design is that it uses for developing potential leads in drug discovery all known theoretical and experimental knowledge of the system under study. The utilization of the knowledge of the molecular basis of the system ultimately aims to reduce human power cost, time saving and laboratory expenses in the drug discovery. In this review paper various strategies applied for systems which include: (i) absence of knowledge of the receptor active site; (ii) the knowledge of a homology model of a receptor, (iii) the knowledge of the experimentally determined (i.e. X-ray crystallography, NMR spectroscopy) coordinates of the active site of the protein in absence and (iv) the presence of the ligand will be analyzed.