Molecular Recognition at the Active Site of Catechol-O-Methyltransferase: Energetically Favorable Replacement of a Water Molecule Imported by a Bisubstrate Inhibitor

The temperature-dependence of host-guest binding thermodynamics: experimental and simulation studies

The thermodynamic parameters of host-guest binding can be used to describe, understand, and predict molecular recognition events in aqueous systems. However, interpreting binding thermodynamics remains challenging, even for these relatively simple molecules, as they are determined by both direct and solvent-mediated host-guest interactions. In this contribution, we focus on the contributions of water to binding by studying binding thermodynamics, both experimentally and computationally, for a series of nearly rigid, electrically neutral host-guest systems and report the temperature-dependent thermodynamic binding contributions ΔGb(T), ΔHb(T), ΔSb(T), and ΔCp,b. Combining isothermal titration calorimetry (ITC) measurements with molecular dynamics (MD) simulations, we provide insight into the binding forces at play for the macrocyclic hosts cucurbit[n]uril (CBn, n = 7-8) and β-cyclodextrin (β-CD) with a range of guest molecules. We find consistently negative changes in heat capacity on binding (ΔCp,b) for all systems studied herein - as well as for literature host-guest systems - indicating increased enthalpic driving forces for binding at higher temperatures. We ascribe these trends to solvation effects, as the solvent properties of water deteriorate as temperature rises. Unlike the entropic and enthalpic contributions to binding, with their differing signs and magnitudes for the classical and non-classical hydrophobic effect, heat capacity changes appear to be a unifying and more general feature of host-guest complex formation in water. This work has implications for understanding protein-ligand interactions and other complex systems in aqueous environments.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

Unveiling the biopathway for the design of novel COMT inhibitors

Catechol-O-methyltransferase (COMT) is an enzyme responsible for the O-methylation of biologically active catechol-based molecules. It has been associated with several neurological disorders, especially Parkinson's disease (PD), because of its involvement in catecholamine metabolism, and has been considered an important therapeutic target for central nervous system disorders. In this review, we summarize the biophysical, structural, and therapeutical relevance of COMT; the medicinal chemistry behind the development of COMT inhibitors and the application of computer-aided design to support the design of novel molecules; current methodologies for the biosynthesis, isolation, and purification of COMT; and revise existing bioanalytical approaches for the assessment of enzymatic activity in several biological matrices.

Discovery of Small Molecules as Membrane-Bound Catechol-O-methyltransferase Inhibitors with Interest in Parkinson's Disease: Pharmacophore Modeling, Molecular Docking and In Vitro Experimental Validation Studies

A pharmacophore-based virtual screening methodology was used to discover new catechol-O-methyltransferase (COMT) inhibitors with interest in Parkinson’s disease therapy. To do so, pharmacophore models were constructed using the structure of known inhibitors and then they were used in a screening in the ZINCPharmer database to discover hit molecules with the desired structural moieties and drug-likeness properties. Following this, the 50 best ranked molecules were submitted to molecular docking to better understand their atomic interactions and binding poses with the COMT (PDB#6I3C) active site. Additionally, the hits’ ADMET properties were also studied to improve the obtained results and to select the most promising compounds to advance for in-vitro studies. Then, the 10 compounds selected were purchased and studied regarding their in-vitro inhibitory potency on human recombinant membrane-bound COMT (MBCOMT), as well as their cytotoxicity in rat dopaminergic cells (N27) and human dermal fibroblasts (NHDF). Of these, the compound ZIN27985035 displayed the best results: For MBCOMT inhibition an IC₅₀ of 17.6 nM was determined, and low cytotoxicity was observed in both cell lines (61.26 and 40.32 μM, respectively). Therefore, the promising results obtained, combined with the structure similarity with commercial COMT inhibitors, can allow for the future development of a potential new Parkinson’s disease drug candidate with improved properties.

Some important inhibitors and mechanisms of rheumatoid arthritis

Rheumatoid arthritis is a chronic disease that seriously affects human health and quality of life, and it is one of the main causes of labor loss and disability. Many countries have listed rheumatoid arthritis as one of the national a key diseases to tackle. The pathogenesis of RA in humans is still unknown, and medical researchers believe that the pathogenesis of RA may be the result of a combination of genetic and environmental factors. RA is an incurable condition that can only be controlled and treated with conventional drugs. In this paper, the pathologic features and pathogenesis of RA were introduced, and the research progress of new anti-rheumatoid arthritis chemical drugs in recent years was reviewed.

Nucleoside-modified AdoMet analogues for differential methyltransferase targeting

Methyltransferases (MTases) modify a wide range of biomolecules using S-adenosyl-l-methionine (AdoMet) as the cosubstrate. Synthetic AdoMet analogues are powerful tools to site-specifically introduce a variety of functional groups and exhibit potential to be converted only by distinct MTases. Extending the size of the substituent at the sulfur/selenium atom provides selectivity among MTases but is insufficient to discriminate between promiscuous MTases. We present a panel of AdoMet analogues differing in the nucleoside moiety (NM-AdoMets). These NM-AdoMets were efficiently produced by a previously uncharacterized methionine adenosyltransferase (MAT) from methionine and ATP analogues, such as ITP and N6-propargyl-ATP. The N6-modification changed the relative activity of three representative MTases up to 13-fold resulting in discrimination of substrates for the methyl transfer and could also be combined with transfer of allyl and propargyl groups.

Targeting Metalloenzymes for Therapeutic Intervention

Metalloenzymes are central to a wide range of essential biological activities, including nucleic acid modification, protein degradation, and many others. The role of metalloenzymes in these processes also makes them central for the progression of many diseases and, as such, makes metalloenzymes attractive targets for therapeutic intervention. Increasing awareness of the role metalloenzymes play in disease and their importance as a class of targets has amplified interest in the development of new strategies to develop inhibitors and ultimately useful drugs. In this Review, we provide a broad overview of several drug discovery efforts focused on metalloenzymes and attempt to map out the current landscape of high-value metalloenzyme targets.

Oxindole-based SYK and JAK3 dual inhibitors for rheumatoid arthritis: designing, synthesis and biological evaluation

Aim: Autoimmune disorders have complex pathophysiology and focus is laid on the development of multitargeted agents. Two well-established kinases: SYK and JAK3, were considered to design dual inhibitors as potential therapeutics using various molecular-modeling approaches. Mehodology: Pharmacophore models for SYK and JAK3 were generated using oxindole-based inhibitors. Furthermore, an in-house database was designed that was screened against the best selected models. The obtained hits were employed for docking analysis and subjected to MM-GBSA analysis and molecular dynamic simulation. Results: Top five oxindole derivatives were synthesized and evaluated for in vitro SYK and JAK3 activity. The most active compound 4a was evaluated for in vivo antiarthritic activity. It showed significant anti-arthritic activity. Conclusion: Thus, the designed inhibitors resulted in potential therapeutic agents for rheumatoid arthritis.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

SAM/SAH Analogs as Versatile Tools for SAM-Dependent Methyltransferases

S-Adenosyl-L-methionine (SAM) is a sulfonium molecule with a structural hybrid of methionine and adenosine. As the second largest cofactor in the human body, its major function is to serve as methyl donor for SAM-dependent methyltransferases (MTases). The resultant transmethylation of biomolecules constitutes a significant biochemical mechanism in epigenetic regulation, cellular signaling, and metabolite degradation. Recently, numerous SAM analogs have been developed as synthetic cofactors to transfer the activated groups on MTase substrates for downstream ligation and identification. Meanwhile, new compounds built upon or derived from the SAM scaffold have been designed and tested as selective inhibitors for important MTase targets. Here, we summarized the recent development and application of SAM analogs as chemical biology tools for MTases.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

Mapping the conformational space accessible to catechol-O-methyltransferase

Methylation catalysed by catechol-O-methyltransferase (COMT) is the main pathway of catechol neurotransmitter deactivation in the prefrontal cortex. Low levels of this class of neurotransmitters are held to be causative of diseases such as schizophrenia, depression and Parkinson's disease. Inhibition of COMT may increase neurotransmitter levels, thus offering a route for treatment. Structure-based drug design hitherto seems to be based on the closed enzyme conformation. Here, a set of apo, semi-holo, holo and Michaelis form crystal structures are described that define the conformational space available to COMT and that include likely intermediates along the catalytic pathway. Domain swaps and sizeable loop movements around the active site testify to the flexibility of this enzyme, rendering COMT a difficult drug target. The low affinity of the co-substrate S-adenosylmethionine and the large conformational changes involved during catalysis highlight significant energetic investment to achieve the closed conformation. Since each conformation of COMT is a bona fide target for inhibitors, other states than the closed conformation may be promising to address. Crystallographic data for an alternative avenue of COMT inhibition, i.e. locking of the apo state by an inhibitor, are presented. The set of COMT structures may prove to be useful for the development of novel classes of inhibitors.

Structure-based drug design of catechol-O-methyltransferase inhibitors for CNS disorders

Catechol-O-methyltransferase (COMT) is of great importance in pharmacology because it catalyzes the metabolism (methylation) of endogenous and xenobiotic catechols. Moreover, inhibition of COMT is the drug target in the management of central nervous system (CNS) disorders such as Parkinson's disease due to its role in regulation of the dopamine level in the brain. The X-ray crystal structures for COMT have been available since 1994. The active sites for cofactor and substrate/inhibitor binding are well resolved to an atomic level, providing valuable insights into the catalytic mechanisms as well as the role of magnesium ions in catalysis. Determination of how the substrates/inhibitors bind to the protein leads to a structure-based approach that has resulted in potent and selective inhibitors. This review focuses on the design of two types of inhibitors (nitrocatechol-type and bisubstrate inhibitors) for COMT using the protein structures.

Flooding enzymes: quantifying the contributions of interstitial water and cavity shape to ligand binding using extended linear response free energy calculations

Glutamate racemase (GR) is a cofactor independent amino acid racemase that has recently garnered increasing attention as an antimicrobial drug target. There are numerous high resolution crystal structures of GR, yet these are invariably bound to either d-glutamate or very weakly bound oxygen-based salts. Recent in silico screens have identified a number of new competitive inhibitor scaffolds, which are not based on d-Glu, but exploit many of the same hydrogen bond donor positions. In silico studies on 1-H-benzimidazole-2-sulfonic acid (BISA) show that the sulfonic acid points to the back of the GR active site, in the most buried region, analogous to the C2-carboxylate binding position in the GR-d-glutamate complex. Furthermore, BISA has been shown to be the strongest nonamino acid competitive inhibitor. Previously published computational studies have suggested that a portion of this binding strength is derived from complexation with a more closed active site, relative to weaker ligands, and in which the internal water network is more isolated from the bulk solvent. In order to validate key contacts between the buried sulfonate moiety of BISA and moieties in the back of the enzyme active site, as well as to probe the energetic importance of the potentially large number of interstitial waters contacted by the BISA scaffold, we have designed several mutants of Asn75. GR-N75A removes a key hydrogen bond donor to the sulfonate of BISA, but also serves to introduce an additional interstitial water, due to the newly created space of the mutation. GR- N75L should also show the loss of a hydrogen bond donor to the sulfonate of BISA, but does not (a priori) seem to permit an additional interstitial water contact. In order to investigate the dynamics, structure, and energies of this water-mediated complexation, we have employed the extended linear response (ELR) approach for the calculation of binding free energies to GR, using the YASARA2 knowledge based force field on a set of ten GR complexes, and yielding an R-squared value of 0.85 and a RMSE of 2.0 kJ/mol. Surprisingly, the inhibitor set produces a uniformly large interstitial water contribution to the electrostatic interaction energy (⟨V^el⟩), ranging from 30 to >50%, except for the natural substrate (d-glutamate), which has only a 7% contribution of ⟨V^el⟩ from water. The broader implications for predicting and exploiting significant interstitial water contacts in ligand–enzyme complexation are discussed.

Optimization of 3,5-dimethylisoxazole derivatives as potent bromodomain ligands

The bromodomain protein module, which binds to acetylated lysine, is emerging as an important epigenetic therapeutic target. We report the structure-guided optimization of 3,5-dimethylisoxazole derivatives to develop potent inhibitors of the BET (bromodomain and extra terminal domain) bromodomain family with good ligand efficiency. X-ray crystal structures of the most potent compounds reveal key interactions required for high affinity at BRD4(1). Cellular studies demonstrate that the phenol and acetate derivatives of the lead compounds showed strong antiproliferative effects on MV4;11 acute myeloid leukemia cells, as shown for other BET bromodomain inhibitors and genetic BRD4 knockdown, whereas the reported compounds showed no general cytotoxicity in other cancer cell lines tested.

From lin-benzoguanines to lin-benzohypoxanthines as ligands for Zymomonas mobilis tRNA-guanine transglycosylase: replacement of protein-ligand hydrogen bonding by importing water clusters

The foodborne illness shigellosis is caused by Shigella bacteria that secrete the highly cytotoxic Shiga toxin, which is also formed by the closely related enterohemorrhagic Escherichia coli (EHEC). It has been shown that tRNA-guanine transglycosylase (TGT) is essential for the pathogenicity of Shigella flexneri. Herein, the molecular recognition properties of a guanine binding pocket in Zymomonas mobilis TGT are investigated with a series of lin-benzohypoxanthine- and lin-benzoguanine-based inhibitors that bear substituents to occupy either the ribose-33 or the ribose-34 pocket. The three inhibitor scaffolds differ by the substituent at C(6) being H, NH(2), or NH-alkyl. These differences lead to major changes in the inhibition constants, pK(a) values, and binding modes. Compared to the lin-benzoguanines, with an exocyclic NH(2) at C(6), the lin-benzohypoxanthines without an exocyclic NH(2) group have a weaker affinity as several ionic protein-ligand hydrogen bonds are lost. X-ray cocrystal structure analysis reveals that a new water cluster is imported into the space vacated by the lacking NH(2) group and by a conformational shift of the side chain of catalytic Asp102. In the presence of an N-alkyl group at C(6) in lin-benzoguanine ligands, this water cluster is largely maintained but replacement of one of the water molecules in the cluster leads to a substantial loss in binding affinity. This study provides new insight into the role of water clusters at enzyme active sites and their challenging substitution by ligand parts, a topic of general interest in contemporary structure-based drug design.

(+)-syn-Benzotriborneol an enantiopure C3-symmetric receptor for water

(+)-syn-Benzotriborneol forms stable complexes with one molecule of water. This is due to the ability of the host to form three hydrogen bonds with water, to act simultaneously as a hydrogen-bond acceptor and donor, and to a perfect geometrical match between the pair. We report experimental (X-ray and neutron diffraction, VT NMR, DSC, TGA) and stereochemical studies carried out to elucidate and quantify the molecular and thermodynamic aspects of this supramolecular complex.

Molecular recognition at the active site of factor Xa: cation-π interactions, stacking on planar peptide surfaces, and replacement of structural water

Factor Xa, a serine protease from the blood coagulation cascade, is an ideal enzyme for molecular recognition studies, as its active site is highly shape-persistent and features distinct, concave sub-pockets. We developed a family of non-peptidic, small-molecule inhibitors with a central tricyclic core orienting a neutral heterocyclic substituent into the S1 pocket and a quaternary ammonium ion into the aromatic box in the S4 pocket. The substituents were systematically varied to investigate cation-π interactions in the S4 pocket, optimal heterocyclic stacking on the flat peptide walls lining the S1 pocket, and potential water replacements in both the S1 and the S4 pockets. Structure-activity relationships were established to reveal and quantify contributions to the binding free enthalpy, resulting from single-atom replacements or positional changes in the ligands. A series of high-affinity ligands with inhibitory constants down to K(i)=2 nM were obtained and their proposed binding geometries confirmed by X-ray co-crystal structures of protein-ligand complexes.

Halogen bonding at the active sites of human cathepsin L and MEK1 kinase: efficient interactions in different environments

In two series of small-molecule ligands, one inhibiting human cathepsin L (hcatL) and the other MEK1 kinase, biological affinities were found to strongly increase when an aryl ring of the inhibitors is substituted with the larger halogens Cl, Br, and I, but to decrease upon F substitution. X-ray co-crystal structure analyses revealed that the higher halides engage in halogen bonding (XB) with a backbone C=O in the S3 pocket of hcatL and in a back pocket of MEK1. While the S3 pocket is located at the surface of the enzyme, which provides a polar environment, the back pocket in MEK1 is deeply buried in the protein and is of pronounced apolar character. This study analyzes environmental effects on XB in protein-ligand complexes. It is hypothesized that energetic gains by XB are predominantly not due to water replacements but originate from direct interactions between the XB donor (Caryl-X) and the XB acceptor (C=O) in the correct geometry. New X-ray co-crystal structures in the same crystal form (space group P2(1)2(1)2(1)) were obtained for aryl chloride, bromide, and iodide ligands bound to hcatL. These high-resolution structures reveal that the backbone C=O group of Gly61 in most hcatL co-crystal structures maintains water solvation while engaging in XB. An aryl-CF3-substituted ligand of hcatL with an unexpectedly high affinity was found to adopt the same binding geometry as the aryl halides, with the CF3 group pointing to the C=O group of Gly61 in the S3 pocket. In this case, a repulsive F2C-F⋅⋅⋅O=C contact apparently is energetically overcompensated by other favorable protein-ligand contacts established by the CF3 group.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.