Free energy perturbation approach to the critical assessment of selective cyclooxygenase-2 inhibitors

Structural and Chemical Biology of the Interaction of Cyclooxygenase with Substrates and Non-Steroidal Anti-Inflammatory Drugs

Cyclooxgenases are key enzymes of lipid signaling. They carry out the first step in the production of prostaglandins, important mediators of inflammation, pain, cardiovascular disease, and cancer, and they are the molecular targets for nonsteroidal anti-inflammatory drugs, which are among the oldest and most chemically diverse set of drugs known. Homodimeric proteins that behave as allosterically modulated, functional heterodimers, the cyclooxygenases exhibit complex kinetic behavior, requiring peroxide-dependent activation and undergoing suicide inactivation. Due to their important physiological and pathophysiological roles and keen interest on the part of the pharmaceutical industry, the cyclooxygenases have been the focus of a vast array of structural studies, leading to the publication of over 80 crystal structures of the enzymes in complex with substrates or inhibitors supported by a wealth of functional data generated by site-directed mutation experiments. In this review, we explore the chemical biology of the cyclooxygenases through the lens of this wealth of structural and functional information. We identify key structural features of the cyclooxygenases, break down their active site into regional binding pockets to facilitate comparisons between structures, and explore similarities and differences in the binding modes of the wide variety of ligands (both substrates and inhibitors) that have been characterized in complex with the enzymes. Throughout, we correlate structure with function whenever possible. Finally, we summarize what can and cannot be learned from the currently available structural data and discuss the critical intriguing questions that remain despite the wealth of information that has been amassed in this field.

Homologation: A Versatile Molecular Modification Strategy to Drug Discovery

Homologation is a concept introduced by Gerhard in 1853 to describe a homologous series in organic chemistry. Since then, the concept has been adapted and used in medicinal chemistry as one of the most important strategies for molecular modification. The homologation types, their influence on physico-chemical properties and molecular conformation are presented and discussed. Its application in lead-identification and lead optimization steps, as well as its impact on pharmacodynamics/pharmacokinetic properties and on protein structure is highlighted from selected examples. • Homologation: definition and types • Homologous series in nature • Comparative physico-chemical and conformational properties • Application in lead-identification and lead-optimization • Impact on pharmacodynamic property • Impact on pharmacokinetic property • Impact on protein structure • Concluding remarks • Acknowledgment • References.

Oxicams, a class of nonsteroidal anti-inflammatory drugs and beyond

Oxicams are a class of nonsteroidal anti-inflammatory drugs (NSAIDs) structurally related to the enolic acid class of 4-hydroxy-1,2-benzothiazine carboxamides. They are used clinically to treat both acute and chronic inflammation by inhibiting the activity of the two cyclooxygenase (COX) isoforms, COX-1 and COX-2. Oxicams are structurally distinct from all other NSAIDs, exhibiting a novel binding pose in the COX active site. The 4-hydroxyl group on the thiazine ring partners with Ser-530 via hydrogen bonding while two coordinated water molecules mediate a polar interaction between the oxicam and COX. The rotation of Leu-531 in the complex opens a new pocket, which is not used for binding other NSAIDs to the enzyme. This structure provides the basis for understanding documented structure-activity relationships within the oxicam class. In addition, from the oxicam template, a series of potent microsomal prostaglandin E synthase-1 (mPGES-1) inhibitors represents a new direction for drug development. Here, we review the major route of oxicam synthesis and structure-activity for COX inhibition, as well as recent advances in oxicam-mediated mPGES-1 inhibition.

Oxicams bind in a novel mode to the cyclooxygenase active site via a two-water-mediated H-bonding Network

Oxicams are widely used nonsteroidal anti-inflammatory drugs (NSAIDs), but little is known about the molecular basis of the interaction with their target enzymes, the cyclooxygenases (COX). Isoxicam is a nonselective inhibitor of COX-1 and COX-2 whereas meloxicam displays some selectivity for COX-2. Here we report crystal complexes of COX-2 with isoxicam and meloxicam at 2.0 and 2.45 angstroms, respectively, and a crystal complex of COX-1 with meloxicam at 2.4 angstroms. These structures reveal that the oxicams bind to the active site of COX-2 using a binding pose not seen with other NSAIDs through two highly coordinated water molecules. The 4-hydroxyl group on the thiazine ring partners with Ser-530 via hydrogen bonding, and the heteroatom of the carboxamide ring of the oxicam scaffold interacts with Tyr-385 and Ser-530 through a highly coordinated water molecule. The nitrogen atom of the thiazine and the oxygen atom of the carboxamide bind to Arg-120 and Tyr-355 via another highly ordered water molecule. The rotation of Leu-531 in the structure opens a novel binding pocket, which is not utilized for the binding of other NSAIDs. In addition, a detailed study of meloxicam·COX-2 interactions revealed that mutation of Val-434 to Ile significantly reduces inhibition by meloxicam due to subtle changes around Phe-518, giving rise to the preferential inhibition of COX-2 over COX-1.

Accurate calculation of mutational effects on the thermodynamics of inhibitor binding to p38α MAP kinase: a combined computational and experimental study

A major current challenge for drug design efforts focused on protein kinases is the development of drug resistance caused by spontaneous mutations in the kinase catalytic domain. The ubiquity of this problem means that it would be advantageous to develop fast, effective computational methods that could be used to determine the effects of potential resistance-causing mutations before they arise in a clinical setting. With this long-term goal in mind, we have conducted a combined experimental and computational study of the thermodynamic effects of active-site mutations on a well-characterized and high-affinity interaction between a protein kinase and a small-molecule inhibitor. Specifically, we developed a fluorescence-based assay to measure the binding free energy of the small-molecule inhibitor, SB203580, to the p38α MAP kinase and used it measure the inhibitor's affinity for five different kinase mutants involving two residues (Val38 and Ala51) that contact the inhibitor in the crystal structure of the inhibitor-kinase complex. We then conducted long, explicit-solvent thermodynamic integration (TI) simulations in an attempt to reproduce the experimental relative binding affinities of the inhibitor for the five mutants; in total, a combined simulation time of 18.5 μs was obtained. Two widely used force fields - OPLS-AA/L and Amber ff99SB-ILDN - were tested in the TI simulations. Both force fields produced excellent agreement with experiment for three of the five mutants; simulations performed with the OPLS-AA/L force field, however, produced qualitatively incorrect results for the constructs that contained an A51V mutation. Interestingly, the discrepancies with the OPLS-AA/L force field could be rectified by the imposition of position restraints on the atoms of the protein backbone and the inhibitor without destroying the agreement for other mutations; the ability to reproduce experiment depended, however, upon the strength of the restraints' force constant. Imposition of position restraints in corresponding simulations that used the Amber ff99SB-ILDN force field had little effect on their ability to match experiment. Overall, the study shows that both force fields can work well for predicting the effects of active-site mutations on small molecule binding affinities and demonstrates how a direct combination of experiment and computation can be a powerful strategy for developing an understanding of protein-inhibitor interactions.

HIV-1 NNRTIs: structural diversity, pharmacophore similarity, and implications for drug design

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) nowadays represent very potent and most promising anti-AIDS agents that specifically target the HIV-1 reverse transcriptase (RT). However, the effectiveness of NNRTI drugs can be hampered by rapid emergence of drug-resistant viruses and severe side effects upon long-term use. Therefore, there is an urgent need to develop novel, highly potent NNRTIs with broad spectrum antiviral activity and improved pharmacokinetic properties, and more efficient strategies that facilitate and shorten the drug discovery process would be extremely beneficial. Fortunately, the structural diversity of NNRTIs provided a wide space for novel lead discovery, and the pharmacophore similarity of NNRTIs gave valuable hints for lead discovery and optimization. More importantly, with the continued efforts in the development of computational tools and increased crystallographic information on RT/NNRTI complexes, structure-based approaches using a combination of traditional medicinal chemistry, structural biology, and computational chemistry are being used increasingly in the design of NNRTIs. First, this review covers two decades of research and development for various NNRTI families based on their chemical scaffolds, and then describes the structural similarity of NNRTIs. We have attempted to assemble a comprehensive overview of the general approaches in NNRTI lead discovery and optimization reported in the literature during the last decade. The successful applications of medicinal chemistry strategies, crystallography, and computational tools for designing novel NNRTIs are highlighted. Future directions for research are also outlined.

Proteomic profiling identifies cyclooxygenase-2-independent global proteomic changes by celecoxib in colorectal cancer cells

Celecoxib, a selective inhibitor of the enzyme cyclooxygenase-2 (COX-2), has been shown to be a promising chemoprevention agent. The chemopreventive efficacy of celecoxib is believed to be a consequence of its COX-2-dependent and COX-2-independent effects on a variety of cellular processes including proliferation, apoptosis, angiogenesis, and immunosurveillance. In an attempt to identify proteomic markers modulated by celecoxib that are independent of its inhibitory effect on COX-2, the colorectal cancer cell line HCT-116, a nonexpresser of COX-2, was treated with celecoxib. We used the powerful, state-of-the-art two-dimensional difference gel electrophoresis technology coupled with mass spectrometric sequencing to compare global proteomic profiles of HCT-116 cells before and after treatment with celecoxib. Among the differentially expressed proteins identified following celecoxib treatment were proteins involved in diverse cellular functions including glycolysis, protein biosynthesis, DNA synthesis, mRNA processing, protein folding, phosphorylation, redox regulation, and molecular chaperon activities. Our study presents a comprehensive analysis of large-scale celecoxib-modulated proteomic alterations, at least some of which may be mechanistically related to the COX-2-independent chemopreventive effect of celecoxib.

Combining docking and molecular dynamic simulations in drug design

A rational approach is needed to maximize the chances of finding new drugs, and to exploit the opportunities of potential new drug targets emerging from genomic and proteomic initiatives, and from the large libraries of small compounds now readily available through combinatorial chemistry. Despite a shaky early history, computer-aided drug design techniques can now be effective in reducing costs and speeding up drug discovery. This happy outcome results from development of more accurate and reliable algorithms, use of more thoughtfully planned strategies to apply them, and greatly increased computer power to allow studies with the necessary reliability to be performed. Our review focuses on applications and protocols, with the main emphasis on critical analysis of recent studies where docking calculations and molecular dynamics (MD) simulations were combined to dock small molecules into protein receptors. We highlight successes to demonstrate what is possible now, but also point out drawbacks and future directions. The review is structured to lead the reader from the simpler to more compute-intensive methods. Thus, while inexpensive and fast docking algorithms can be used to scan large compound libraries and reduce their size, more accurate but expensive MD simulations can be applied when a few selected ligand candidates remain. MD simulations can be used: during the preparation of the protein receptor before docking, to optimize its structure and account for protein flexibility; for the refinement of docked complexes, to include solvent effects and account for induced fit; to calculate binding free energies, to provide an accurate ranking of the potential ligands; and in the latest developments, during the docking process itself to find the binding site and correctly dock the ligand a priori.