Fluoroaromatic-fluoroaromatic interactions between inhibitors bound in the crystal lattice of human carbonic anhydrase II

Development and application of decatungstate catalyzed C-H 18F- and 19F-fluorination, fluoroalkylation and beyond

Over the past few decades, photocatalytic C-H functionalization reactions have received increasing attention due to the often mild reaction conditions and complementary selectivities to conventional functionalization processes. Now, photocatalytic C-H functionalization is a widely employed tool, supporting activities ranging from complex molecule synthesis to late-stage structure-activity relationship studies. In this perspective, we will discuss our efforts in developing a photocatalytic decatungstate catalyzed C-H fluorination reaction as well as its practical application realized through collaborations with industry partners at Hoffmann-La Roche and Merck, and extension to radiofluorination with radiopharmaceutical chemists and imaging experts at TRIUMF and the BC Cancer Agency. Importantly, we feel that our efforts address a question of utility posed by Professor Tobias Ritter in "Late-Stage Fluorination: Fancy Novelty or Useful Tool?" (ACIE, 2015, 54, 3216). In addition, we will discuss decatungstate catalyzed C-H fluoroalkylation and the interesting electrostatic effects observed in decatungstate-catalyzed C-H functionalization. We hope this perspective will inspire other researchers to explore the use of decatungstate for the purposes of photocatalytic C-H functionalization and further advance the exploitation of electrostatic effects for both rate acceleration and directing effects in these reactions.

Transient Formation of Hemiketals from Pentafluoro-gem-diols in the Presence of Alcohols

Pentafluoro-gem-diols have emerged as a source of reactive intermediates for synthesizing fluorinated molecules. When pentafluoro-gem-diols were exposed to alcohols as solvents, the formation of transient hemiketals was detected by 19F NMR. The conversion rates to hemiketals were found to be higher with primary alcohols than with secondary or fluorinated alcohols. These findings provide valuable insight for developing novel techniques to construct intricate fluorinated structures using pentafluoro-gem-diols.

Aromatic Ring Fluorination Patterns Modulate Inhibitory Potency of Fluorophenylhydroxamates Complexed with Histone Deacetylase 6

Bavarostat (EKZ-001) is a selective inhibitor of histone deacetylase 6 (HDAC6) that contains a meta-fluorophenylhydroxamate Zn2+-binding group. The recently determined crystal structure of its complex with HDAC6 from Danio rerio (zebrafish) revealed that the meta-fluoro substituent binds exclusively in an aromatic crevice defined by F583 and F643 rather than being oriented out toward solvent. To explore the binding of inhibitor C-F groups in this fluorophilic crevice, we now report a series of 10 simple fluorophenylhydroxamates bearing one or more fluorine atoms with different substitution patterns. Inhibitory potencies against human and zebrafish HDAC6 range widely from 121 to >30,000 nM. The best inhibitory potency is measured for meta-difluorophenylhydroxamate (5) with IC50 = 121 nM against human HDAC6; the worst inhibitory potencies are measured for ortho-fluorophenylhydroxamate (1) as well as fluorophenylhydroxamates 4, 7, 9, and 10, although there are some variations in activity trends against human and zebrafish HDAC6. These studies show that aromatic ring fluorination at the meta position(s) does not improve inhibitory activity against human HDAC6 relative to the nonfluorinated parent compound phenylhydroxamate (IC50 = 120 nM), but meta-fluorination does not seriously compromise inhibitory activity either. Crystal structures of selected zebrafish HDAC6-fluorophenylhydroxamate complexes reveal that the fluoroaromatic ring is uniformly accommodated in the F583-F643 aromatic crevice, so ring fluorination does not perturb the inhibitor binding conformation. However, hydroxamate-Zn2+ coordination is bidentate for some inhibitors and monodentate for others. These studies will inform design strategies underlying the design of 18F-labeled HDAC6 inhibitors intended for positron emission tomography.

Computational Prediction of the Binding Pose of Metal-Binding Pharmacophores

Computational modeling of inhibitors for metalloenzymes in virtual drug development campaigns has proven challenging. To overcome this limitation, a technique for predicting the binding pose of metal-binding pharmacophores (MBPs) is presented. Using a combination of density functional theory (DFT) calculations and docking using a genetic algorithm, inhibitor binding was evaluated in silico and compared with inhibitor-enzyme cocrystal structures. The predicted binding poses were found to be consistent with the cocrystal structures. The computational strategy presented represents a useful tool for predicting metalloenzyme-MBP interactions.

Small-molecule MDM2/X inhibitors and PROTAC degraders for cancer therapy: advances and perspectives

Blocking the MDM2/X–P53 protein–protein interaction has been widely recognized as an attractive therapeutic strategy for the treatment of cancers. Numerous small-molecule MDM2 inhibitors have been reported since the release of the structure of the MDM2–P53 interaction in 1996, SAR405838, NVP-CGM097, MK-8242, RG7112, RG7388, DS-3032b, and AMG232 currently undergo clinical evaluation for cancer therapy. This review is intended to provide a comprehensive and updated overview of MDM2 inhibitors and proteolysis targeting chimera (PROTAC) degraders with a particular focus on how these inhibitors or degraders are identified from starting points, strategies employed, structure–activity relationship (SAR) studies, binding modes or co-crystal structures, biochemical data, mechanistic studies, and preclinical/clinical studies. Moreover, we briefly discuss the challenges of designing MDM2/X inhibitors for cancer therapy such as dual MDM2/X inhibition, acquired resistance and toxicity of P53 activation as well as future directions.

Exploring Structural Determinants of Inhibitor Affinity and Selectivity in Complexes with Histone Deacetylase 6

Inhibition of histone deacetylase 6 (HDAC6) has emerged as a promising therapeutic strategy for the treatment of cancer, chemotherapy-induced peripheral neuropathy, and neurodegenerative disease. The recent X-ray crystal structure determination of HDAC6 enables an understanding of structural features directing affinity and selectivity in the active site. Here, we present the X-ray crystal structures of five HDAC6-inhibitor complexes that illuminate key molecular features of the inhibitor linker and capping groups that facilitate and differentiate binding to HDAC6. In particular, aromatic and heteroaromatic linkers nestle within an aromatic cleft defined by F583 and F643, and different aromatic linkers direct the capping group toward shallow pockets defined by the L1 loop, the L2 loop, or somewhere in between these pockets. These results expand our understanding of factors contributing to the selective inhibition of HDAC6, particularly regarding interactions that can be targeted in the region of the L2 pocket.

Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

The shielding effect of metal complexes on the binding affinities of ligands to metalloproteins

The contributions of metal–ligand interactions to the ligand binding affinities are largely reduced by the shielding effects of metal complexes.

Effect of π-π interaction in Bergman cyclisation

The role of π-π interactions in controlling the reactivity and selectivity of a chemical reaction is only recently being explored, even though their ubiquitous role in the structural aspects is well known. We have studied Bergman cyclisation focusing on the effect of π-π interactions on the activation barrier and the variation of π-π interactions along the reaction coordinate. We used enediyne substrates that contain phenyl groups connected to the reaction centres (C1 and C6 atoms), separated by 0, 1 and 2 linker groups. The main difference between the substrates is that the Ph groups enjoy different flexibility to accommodate the changes occurring during the progress of the reaction. The path length of the minimum energy path is increased - shortest in the least flexible substrate (a) and longer in the more flexible ones (c, d and e). We calculated the interaction between the Ph groups, the π-π interaction, using BP86-D3BJ, B3LYP-D3BJ, M06-2X, B2PLYP-D3BJ, SCS-MP2, and SAPT. The BP86-D3BJ was found to be sufficiently accurate with a mean absolute deviation of 0.26 kcal mol(-1) with respect to the SAPT2+3 values. The variation in the π-π interaction shows different behaviour in a-e, and this can be correlated with the flexibility of the Ph groups to orient themselves to maintain the optimal relative orientation while conforming to the changes in the reaction coordinate. We analysed the relative orientation of the phenyl groups using certain geometric parameters that showed that when Ph groups can attain a relative orientation close to that of the free dimer, the interaction is maximum. Energy decomposition analysis using SAPT showed that the dispersive interaction is the major contributor (50-60%) to the attractive forces. The π-π interactions influenced the overall activation energy, either by destabilising the substrates or by stabilising the TS - resulting in a variation of about 3.5 kcal mol(-1) in activation energies in a-e. The effect of substituents of different electronic nature was assessed which showed that electron donating and electron withdrawing substituents increase the π-π interactions; however, the TS is more stabilised and hence activation energies are increased.

Decorating tetrathiafulvalene (TTF) with fluorinated phenyls through sulfur bridges: facile synthesis, properties, and aggregation through fluorine interactions

Tetrathiafulvalene derivatives (TTF1-TTF9) bearing fluorinated phenyl groups attached through the sulfur bridges have been synthesized by employing a copper-mediated C-S coupling reaction of C6 H5-x Fx I (x=1, 2, 5) and a zinc-thiolate complex, (TBA)2 [Zn(DMIT)2 ] (TBA=tetrabutyl ammonium, DMIT=1,3-dithiole-2-thione-4,5-dithiolate), as the key step. Particularly, the selective synthesis of C6 F5 -substituted (TTF8) and C6 F4 -fused (TTF9) TTFs from C6 F5 I is disclosed. The physicochemical properties and crystal structures of these TTFs are fully investigated by UV/Vis absorption spectra, cyclic voltammetry, molecular orbital calculation, and single-crystal X-ray diffraction. The exchange of hydrogen versus fluorine on the peripheral phenyl groups show a notable influence on both the electronic and crystallographic natures of the resulting TTFs: 1) lowering both the HOMO and the LUMO energy levels, 2) modulating the electrochemical properties by regioselective and/or the degree of fluorination, 3) enhancing the driving forces of stacking by multiple fluorine interactions (F⋅⋅⋅S, CF⋅⋅⋅π/πF , CF⋅⋅⋅FC, and CF⋅⋅⋅H). This work indicates that the decoration with fluorinated phenyls holds promise to produce functional TTFs with novel electronic and aggregation features.

Anti-parallel sheet structures of side-chain-free γ-, δ-, and ε-dipeptides stabilized by benzene–pentafluorobenzene stacking

Benzene–pentafluorobenzene stacking can guide ω-amino acid dipeptides to arrange in an anti-parallel manner.

Introduction of fluorine and fluorine-containing functional groups

Over the past decade, the most significant, conceptual advances in the field of fluorination were enabled most prominently by organo- and transition-metal catalysis. The most challenging transformation remains the formation of the parent C-F bond, primarily as a consequence of the high hydration energy of fluoride, strong metal-fluorine bonds, and highly polarized bonds to fluorine. Most fluorination reactions still lack generality, predictability, and cost-efficiency. Despite all current limitations, modern fluorination methods have made fluorinated molecules more readily available than ever before and have begun to have an impact on research areas that do not require large amounts of material, such as drug discovery and positron emission tomography. This Review gives a brief summary of conventional fluorination reactions, including those reactions that introduce fluorinated functional groups, and focuses on modern developments in the field.

Crystal structure of arginase from Leishmania mexicana and implications for the inhibition of polyamine biosynthesis in parasitic infections

Arginase from parasitic protozoa belonging to the genus Leishmania is a potential drug target for the treatment of leishmaniasis because this binuclear manganese metalloenzyme catalyzes the first committed step in the biosynthesis of polyamines that enable cell growth and survival. The high resolution X-ray crystal structures of the unliganded form of Leishmania mexicana arginase (LmARG) and four inhibitor complexes are now reported. These complexes include the reactive substrate analogue 2(S)-amino-6-boronohexanoic acid (ABH) and the hydroxylated substrate analogue nor-N(ω)-hydroxy-l-arginine (nor-NOHA), which are the most potent arginase inhibitors known to date. Comparisons of the LmARG structure with that of the archetypal arginase, human arginase I, reveal that all residues important for substrate binding and catalysis are strictly conserved. However, three regions of tertiary structure differ between the parasitic enzyme and the human enzyme corresponding to the G62 - S71, L161 - C172, and I219 - V230 segments of LmARG. Additionally, variations are observed in salt link interactions that stabilize trimer assembly in LmARG. We also report biological studies in which we demonstrate that localization of LmARG to the glycosome, a unique subcellular organelle peculiar to Leishmania and related parasites, is essential for robust pathogenesis.

Carbonic Anhydrase Inhibitors: X-ray Crystallographic Structure of the Adduct of Human Isozyme II with the Perfluorobenzoyl Analogue of Methazolamide. Implications for the Drug Design of Fluorinated Inhibitors

The X-ray crystal structure for the adduct of human carbonic anhydrase (hCA) II with 4-methyl-5-perfluorophenylcarboximido-delta2-1,3,4-thiadiazoline-2-sulfonamide (PFMZ), a topically acting antiglaucoma sulfonamide, has been resolved at a resolution of 1.8 A. This compound is almost 10 times more effective as a hCA II inhibitor (KI of 1.5 nM) compared to the lead molecule, methazolamide, a clinically used drug (KI of 14 nM). Its binding to the enzyme active site is similar to that of other sulfonamide inhibitors, considering the interactions of the sulfonamide zinc anchoring group and thiadiazoline ring contacts, but differs considerably when the perfluorobenzoylimino fragment of the molecule is analyzed. Indeed, several unprecedented strong hydrogen bonds involving the imino nitrogen, carbonyl oxygen, a fluorine atom in the ortho position of the inhibitor, and two water molecules, as well as Gln 92 of the enzyme active site were seen. A stacking interaction of the perfluorophenyl ring of the inhibitor and the aromatic ring of Phe 131 was also observed for the first time in a CA-sulfonamide adduct. All these findings prove that more potent CA inhibitors incorporating perfluoroaryl/alkyl tails may be designed, with potentially improved antiglaucoma properties, in view of the new types of interactions seen here between the enzyme and the perfluorobenzoylated analogue of methazolamide.

A general approach for developing system-specific functions to score protein-ligand docked complexes using support vector inductive logic programming

Despite the increased recent use of protein-ligand and protein-protein docking in the drug discovery process due to the increases in computational power, the difficulty of accurately ranking the binding affinities of a series of ligands or a series of proteins docked to a protein receptor remains largely unsolved. This problem is of major concern in lead optimization procedures and has lead to the development of scoring functions tailored to rank the binding affinities of a series of ligands to a specific system. However, such methods can take a long time to develop and their transferability to other systems remains open to question. Here we demonstrate that given a suitable amount of background information a new approach using support vector inductive logic programming (SVILP) can be used to produce system-specific scoring functions. Inductive logic programming (ILP) learns logic-based rules for a given dataset that can be used to describe properties of each member of the set in a qualitative manner. By combining ILP with support vector machine regression, a quantitative set of rules can be obtained. SVILP has previously been used in a biological context to examine datasets containing a series of singular molecular structures and properties. Here we describe the use of SVILP to produce binding affinity predictions of a series of ligands to a particular protein. We also for the first time examine the applicability of SVILP techniques to datasets consisting of protein-ligand complexes. Our results show that SVILP performs comparably with other state-of-the-art methods on five protein-ligand systems as judged by similar cross-validated squares of their correlation coefficients. A McNemar test comparing SVILP to CoMFA and CoMSIA across the five systems indicates our method to be significantly better on one occasion. The ability to graphically display and understand the SVILP-produced rules is demonstrated and this feature of ILP can be used to derive hypothesis for future ligand design in lead optimization procedures. The approach can readily be extended to evaluate the binding affinities of a series of protein-protein complexes.

Fluorine in pharmaceuticals: looking beyond intuition

Fluorine substituents have become a widespread and important drug component, their introduction facilitated by the development of safe and selective fluorinating agents. Organofluorine affects nearly all physical and adsorption, distribution, metabolism, and excretion properties of a lead compound. Its inductive effects are relatively well understood, enhancing bioavailability, for example, by reducing the basicity of neighboring amines. In contrast, exploration of the specific influence of carbon-fluorine single bonds on docking interactions, whether through direct contact with the protein or through stereoelectronic effects on molecular conformation of the drug, has only recently begun. Here, we review experimental progress in this vein and add complementary analysis based on comprehensive searches in the Cambridge Structural Database and the Protein Data Bank.

Carbonic anhydrase inhibitors: Inhibition of human, bacterial, and archaeal isozymes with benzene-1,3-disulfonamides—Solution and crystallographic studies

Three benzene-1,3-disulfonamide derivatives were investigated for their interaction with 12 mammalian alpha-carbonic anhydrases (CAs, EC 4.2.1.1), and three bacterial/archaeal CAs belonging to the alpha-, beta-, and gamma-CA class, respectively. X-ray crystal structure of the three inhibitors in complex with the dominant human isozyme CA II revealed a particular binding mode within the cavity. The sulfonamide group in meta-position to the Zn(2+)-coordinated SO(2)NH(2) moiety was oriented toward the hydrophilic side of the active site cleft, establishing hydrogen bonds with His64, Asn67, Gln92, and Thr200. The plane of the phenyl moiety of the inhibitors was rotated by 45 degrees and tilted by 10 degrees with respect to its most recurrent orientation in other CA II-sulfonamide complexes.

Understanding of Assembly Phenomena by Aromatic−Aromatic Interactions: Benzene Dimer and the Substituted Systems

Interactions involving aromatic rings are important in molecular/biomolecular assembly and engineering. As a consequence, there have been a number of investigations on dimers involving benzene or other substituted pi systems. In this Feature Article, we examine the relevance of the magnitudes of their attractive and repulsive interaction energy components in governing the geometries of several pi-pi systems. The geometries and the associated binding energies were evaluated at the complete basis set (CBS) limit of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using a least biased scheme for the given data set. The results for the benzene dimer indicate that the floppy T-shaped structure (center-to-center distance: 4.96 A, with an axial benzene off-centered above the facial benzene) is isoenergetic in zero-point-energy (ZPE) corrected binding energy (D0) to the displaced-stacked structure (vertical interplanar distance: 3.54 A). However, the T-shaped structure is likely to be slightly more stable (D0 approximately equal to 2.4-2.5 kcal/mol) if quadruple excitations are included in the coupled cluster calculations. The presence of substituents on the aromatic ring, irrespective of their electron withdrawing or donating nature, leads to an increase in the binding energy, and the displaced-stacked conformations are more stabilized than the T-shaped conformers. This explains the wide prevalence of displaced stacked structures in organic crystals. Despite that the dispersion energy is dominating, the substituent as well as the conformational effects are correlated to the electrostatic interaction. This electrostatic origin implies that the substituent effect would be reduced in polar solution, but important in apolar media, in particular, for assembling processes.

Outliers in SAR and QSAR: is unusual binding mode a possible source of outliers?

A lead optimization is usually carried out by structure-activity relationship (SAR) and/or quantitative structure-activity relationship (QSAR) studies. One of the assumptions in SAR and QSAR studies is that similar analogs bind to the same binding site in a similar binding mode. One often observes that there are outliers, especially in QSAR. However, most QSAR studies are carried out focusing their attention to the development of QSAR and leave the outliers without much attention. We searched a number of ligand-bound X-ray crystal structures from the protein structure database to find evidences that could indicate a possible source of outliers in SAR or QSAR. Our results show that unusual binding mode could be a source of outliers.

Ultrahigh resolution crystal structures of human carbonic anhydrases I and II complexed with "two-prong" inhibitors reveal the molecular basis of high affinity

The atomic-resolution crystal structures of human carbonic anhydrases I and II complexed with "two-prong" inhibitors are reported. Each inhibitor contains a benzenesulfonamide prong and a cupric iminodiacetate (IDA-Cu(2+)) prong separated by linkers of different lengths and compositions. The ionized NH(-) group of each benzenesulfonamide coordinates to the active site Zn(2+) ion; the IDA-Cu(2+) prong of the tightest-binding inhibitor, BR30, binds to H64 of CAII and H200 of CAI. This work provides the first evidence verifying the structural basis of nanomolar affinity measured for two-prong inhibitors targeting the carbonic anhydrases.

Substituent Effects on the Edge-to-Face Aromatic Interactions

The edge-to-face interactions for either axially or facially substituted benzenes are investigated by using ab initio calculations. The predicted maximum energy difference between substituted and unsubstituted systems is approximately 0.7 kcal/mol (approximately 1.2 kcal/mol if substituents are on both axially and facially substituted positions). In the case of axially substituted aromatic systems, the electron density at the para position is an important stabilizing factor, and thus the stabilization/destabilization by substitution is highly correlated to the electrostatic energy. This results in its subsequent correlation with the polarization and charge transfer. Thus, the stabilization/destabilization by substitution is represented by the sum of electrostatic energy and induction energy. On the other hand, the facially substituted aromatic system depends on not only the electron-donating ability responsible for the electrostatic energy but also the dispersion interaction and exchange repulsion. Although the dispersion energy is the most dominating interaction in both axial and facial substitutions, it is almost canceled by the exchange repulsion in the axial substitution, whereas in the facial substitution, together with the exchange repulsion it augments the electrostatic energy. The systems with electron-accepting substituents (NO2, CN, Br, Cl, F) favor the axial substituent conformation, while those with electron-donating substituents (NH2, CH3, OH) favor the facial substituent conformation. The interactions for the T-shape complex systems of an aromatic ring with other counterpart such as H2, H2O, HCl, and HF are also studied.

Evaluation of library ranking efficacy in virtual screening

We present the results of a comprehensive study in which we explored how the docking procedure affects the performance of a virtual screening approach. We used four docking engines and applied 10 scoring functions to the top-ranked docking solutions of seeded databases against six target proteins. The scores of the experimental poses were placed within the total set to assess whether the scoring function required an accurate pose to provide the appropriate rank for the seeded compounds. This method allows a direct comparison of library ranking efficacy. Our results indicate that the LigandFit/Ligscore1 and LigandFit/GOLD docking/scoring combinations, and to a lesser degree FlexX/FlexX, Glide/Ligscore1, DOCK/PMF (Tripos implementation), LigandFit1/Ligscore2 and LigandFit/PMF (Tripos implementation) were able to retrieve the highest number of actives at a 10% fraction of the database when all targets were looked upon collectively. We also show that the scoring functions rank the observed binding modes higher than the inaccurate poses provided that the experimental poses are available. This finding stresses the discriminatory ability of the scoring algorithms, when better poses are available, and suggests that the number of false positives can be lowered with conformers closer to bioactive ones.

Carbonic anhydrase inhibitors: X-ray crystallographic structure of the adduct of human isozyme II with a topically acting antiglaucoma sulfonamide

The X-ray crystal structure for the adduct of human carbonic anhydrase (hCA) II with a topically acting antiglaucoma sulfonamide (the 2-N,N-diethylaminoethylamide of 5-(4-carboxybenzenesulfonamido-1,3,4-thiadiazole-2-sulfonamide), has been resolved at a resolution of 1.6A. This compound is a very potent inhibitor of the physiologically most relevant isozyme hCA II for the secretion of aqueous humor within the eye K(I) of 1.4 nM), and in animal models of glaucoma showed very effective intraocular pressure (IOP) lowering after topical administration. Surprisingly, the inhibitor bound within the enzyme active site is in the sulfonylimido-4H- delta(2)-1,3,4-thiadiazoline tautomeric form. The inhibitor is directly bound to the Zn(II) ion of the enzyme through the deprotonated primary sulfonamide moiety, participating to the classical hydrogen bond network involving residues of the zinc-binding function and Thr 199 and Glu 106. The 1,3,4-thiadiazoline fragment of the inhibitor makes two hydrogen bonds with the active site residue Thr 200, the secondary sulfonamide moiety makes two hydrogen bonds involving a water molecule and the residue Gln 92, whereas the phenyl ring of the inhibitor participates to an edge-to-face interaction with the phenyl ring of Phe 131, the two cycles being almost perfectly perpendicular to each other. The tertiary amine fragment of the carboxamido tail and the carboxamido moiety itself make hydrogen bonds with water molecules present at the rim of the active site entrance and van der Waals contacts with His 4, Trp 5, and Phe 20. All these multiple interactions never evidenced previously in CA-sulfonamide complexes, explain the very high affinity of this inhibitor for the hCA II active site and may allow further optimization of this class of inhibitors.

Fluorine in medicinal chemistry

Fluorinated compounds are synthesized in pharmaceutical research on a routine basis and many marketed compounds contain fluorine. The present review summarizes some of the most frequently employed strategies for using fluorine substituents in medicinal chemistry. Quite often, fluorine is introduced to improve the metabolic stability by blocking metabolically labile sites. However, fluorine can also be used to modulate the physicochemical properties, such as lipophilicity or basicity. It may exert a substantial effect on the conformation of a molecule. Increasingly, fluorine is used to enhance the binding affinity to the target protein. Recent 3D-structure determinations of protein complexes with bound fluorinated ligands have led to an improved understanding of the nonbonding protein-ligand interactions that involve fluorine.

Carbonic anhydrase inhibitors: the first selective, membrane-impermeant inhibitors targeting the tumor-associated isozyme IX

The inhibition of the tumor-associated transmembrane carbonic anhydrase IX (CA IX) isozyme possessing an extracellular active site has been investigated with a series of positively-charged, pyridinium derivatives of sulfanilamide, homosulfanilamide and 4-aminoethylbenzenesulfonamide. Inhibition data for the physiologically relevant isozymes I and II (cytosolic forms) and IV (membrane-bound) were also provided for comparison. A very interesting inhibition profile against CA IX with these sulfonamides has been observed. Several nanomolar (K(i)'s in the range of 6-54 nM) CA IX inhibitors have also been detected. Because CA IX is a highly active isozyme predominantly expressed in tumor tissues with bad prognosis of disease progression, this finding is very promising for the potential design of CA IX-specific inhibitors with applications as anti-tumor agents. This is the first report of inhibitors that may selectively target CA IX, due to their membrane-impermeability and high affinity for this clinically relevant isozyme.

Carbonic anhydrase inhibitors: X-ray crystallographic structure of the adduct of human isozyme II with EMATE, a dual inhibitor of carbonic anhydrases and steroid sulfatase

The X-ray crystal structure for the adduct of human carbonic anhydrase II (hCA II) with estrone-3-O-sulfamate (EMATE), an antiendocrine agent showing both CA and estrone sulfatase inhibitory properties, has been resolved at a resolution of 1.5A. Its binding to the enzyme is similar to that of other sulfamates/sulfonamides, considering the interactions of the zinc anchoring group, but differs considerably when the steroidal scaffold of the inhibitor is analyzed. This part of the inhibitor interacts only within the hydrophobic half of the CA active site, interacting with residues Val 121, Phe 131, Val 135 and Pro 202, and leaving the hydrophilic half able to accommodate several water molecules not present in the uncomplexed enzyme. In addition, a very short bond of 1.78A between the zinc ion and the coordinated nitrogen atom of the sulfamate moiety is observed, which may explain the high affinity of this inhibitor for hCA II (K(i) of 10nM).

A tree-based algorithm for determining the effects of solvation on the structure of salivary gland tripeptide NH3+-D-PHE-D-GLU-GLY-COO-

A D-enantiomeric analog of the submandibular gland rat-1 tripeptide FEG (Seq: NH(3)(+)-Phe-Glu-Gly-COO(-)) called feG (Seq: NH(3)(+)-D-Phe-D-Glu-Gly-COO(-)) was examined by molecular dynamics simulations in water. Previous in vacuo simulations suggested a conformation consisting predominantly of interactions between the Phe side chain and glutamyl-carboxyl group and a carboxyl/amino termini interaction. The solvated peptide was simulated using two approaches which were compared-a single 400-ns simulation and a "simulation tree." The "tree" approach utilized 45 10-ns simulations with different conformations used as initial structures for given trajectories. We demonstrate that multiple short duration simulations are able to describe the same conformational space as that described by longer simulations. Furthermore, previously described in vacuo interactions were confirmed with amendments: the previously described head-to-tail arrangement of the amino and carboxyl termini, was not observed; the interaction between the glutamyl carboxyl and Phe side chain describes only one of a continuum of conformations present wherein the aromatic residue remains in close proximity to the glutamyl carbonyl group, and also interacts with either of the two available carboxyl groups. Finally, utilizing only two separate 10-ns trajectories, we were able to better describe the conformational space than a single 60-ns trajectory, realizing a threefold decrease in the computational complexity of the problem.