Logo: to the web site of Uppsala University

Cyclooxygenase-1 (COX-1) is one of the main targets of most pain-relieving pharmaceuticals. Although the enzyme is well characterized, it is known to be a difficult target for automated molecular docking and scoring. We collected from the literature a structurally diverse set of 45 nonsteroidal anti-inflammatory drugs (NSAIDs) and COX-2-selective inhibitors (coxibs) with a wide range of binding affinities for COX-1. The binding of this data set to a homology model of human COX-1 was analyzed with different combinations of molecular docking algorithms, scoring functions, and the linear interaction energy (LIE) method for estimating binding affinities. It is found that the computational protocols for estimation of binding affinities are extremely sensitive to the initial orientations of the ligands in the binding pocket. To overcome this limitation, we propose a systematic exploration of docking poses using the LIE calculations as a postscoring function. This scheme yields predictions in excellent agreement with experiment, with a mean unsigned error of 0.9 kcal/mol for binding free energies and structures of high quality. A significant improvement of the results is also seen when averaging over experimental data from several independent measurements.

Nonsteroidal anti-inflammatory drugs (NSAIDs) are widely used for the treatment of pain, fever, inflammation, and some types of cancers. Their mechanism of action is the inhibition of isoforms 1 and 2 of the enzyme cyclooxygenase (COX-1 and COX-2, respectively). However, both nonselective and selective NSAIDs may have side effects that include gastric intestinal bleeding, peptic ulcer formation, kidney problems, and occurrences of myocardial infarction. The search for selective high-affinity COX inhibitors resulted in a number of compounds characterized by a slow, tight-binding inhibition that occurs in a two-step manner. It has been suggested that the final, only very slowly reversible, tight-binding event is the result of conformational changes in the enzyme. However, the nature of these conformational changes has remained elusive. Here we explore the structural determinants of the tight-binding phenomenon in COX-1 with molecular dynamics and free energy simulations. The calculations reveal how different classes of inhibitors affect the equilibrium between two conformational substates of the enzyme in distinctly different ways. The class of tight-binding inhibitors is found to exclusively stabilize an otherwise unfavorable enzyme conformation and bind significantly stronger to this state than to that normally observed in crystal structures. By also computing free energies of binding to the two enzyme conformations for 16 different NSAIDs, we identify an induced-fit mechanism and the key structural features associated with high-affinity tight binding. These results may facilitate the rational development of new COX inhibitors with improved selectivity profiles.

Time-dependent inhibition: of the cyclooxygenases (COX) by a range of nonsteroidal-anti-inflammatory drugs has been described since the first experimental assays of COX were performed. Slow tight-binding inhibitors of COX-1 bind in a two-step mechanism in which the EI -> EI* transition is slow and practically irreversible. Since then, various properties of the inhibitors have :been proposed to cause or affect the time dependency. Conformational changes :1-7, in the enzyme have also been proposed to cause the time, dependency, but no particular structural feature has been identified. Here, we investigated a series of inhibitors of COX 1 that are either time-independent or time-dependent using a combination of molecular dynamics simulations, binding free energy calculations, and potential of mean force calculations. We find that the time-dependent inhibitors stabilize a conformational change in the enzyme mainly identified by the rotation of a leucine Side chain adjacent to the binding pocket. The induced conformation has been previously Shown to be essential for the high binding affinities of tight-binding inhibitors in COX-1. The results of this work show that the structural features of the enzyme involved in both time-dependent and tight binding inhibition are identical and further;identify a structural mechanism responsible for the transition between the two enzyme inhibitor complexes characteristic of slow tight-binding COX-1 inhibitors.

Nonsteroidal anti-inflammatory drugs (NSAIDs) inhibit cyclooxygenase (COX) 1 and 2 with varying degrees of selectivity. A group of COX-2 selective inhibitors —coxibs— bind in a time-dependent manner through a three-step mechanism, utilizing a side-pocket in the binding site. Coxibs have been extensively probed to identify the structural features regulating the slow tight-binding mechanism responsible for COX-2 selectivity. In this study, we further probe a structurally and kinetically diverse data set of COX inhibitors by molecular dynamics and free energy simulations. We found that the features regulating the high affinities associated with time-dependency in COX depend on the inhibitor kinetics. In particular, most time-dependent inhibitors share a common structural mechanism, consisting in an induced-fit rotation of the side-chain of Leu531 in the main binding pocket. The high affinities of two-step slow tight-binding inhibitors and some slow reversible inhibitors can be thus explained by the increased space in the main binding pocket after this rotation. Coxibs that belong to a separate class of slow tight-binding inhibitors benefit more from the displacement of the neighboring side-chain of Arg513, exclusive of the COX-2 side-pocket. This displacement further stabilizes the aforementioned rotation of Leu531, and can explain the selectivity of coxibs for COX-2.

Angiotensin IV (Ang IV) and related peptide analogues, as well as non-peptide inhibitors of insulin-regulated aminopeptidase (IRAP), have previously been shown to enhance memory and cognition in animal models. Furthermore, the endogenous IRAP substrates oxytocin and vasopressin are known to facilitate learning and memory. In this study, the two recently synthesized 13-membered macrocylic competitive IRAP inhibitors HA08 and HA09, which were designed to mimic the N-terminal of oxytocin and vasopressin, were assessed and compared based on their ability to bind to the IRAP active site, and alter dendritic spine density in rat hippocampal primary cultures. The binding modes of the IRAP inhibitors HA08, HA09 and of Ang IV in either the extended or γ-turn conformation at the C-terminal to human IRAP were predicted by docking and molecular dynamics (MD) simulations. The binding free energies calculated with the linear interaction energy (LIE) method, which are in excellent agreement with experimental data and simulations, have been used to explain the differences in activities of the IRAP inhibitors, both of which are structurally very similar, but differ only with regard to one stereogenic center. In addition, we show that HA08, which is 100-fold more potent than the epimer HA09, can enhance dendritic spine number and alter morphology, a process associated with memory facilitation. Therefore, HA08, one of the most potent IRAP inhibitors known today, may serve as a suitable starting point for medicinal chemistry programs aided by MD simulations aimed at discovering more drug-like cognitive enhancers acting via augmenting synaptic plasticity.

The zinc metallopeptidase insulin regulated aminopeptidase (IRAP), which is highly expressed in the hippocampus and other brain regions associated with cognitive function, has been identified as a high-affinity binding site of the hexapeptide angiotensin IV (Ang IV). This hexapeptide is thought to facilitate learning and memory by binding to the catalytic site of IRAP to inhibit its enzymatic activity. In support of this hypothesis, low molecular weight, nonpeptide specific inhibitors of TRAP have been shown to enhance memory in rodent models. Recently, it was demonstrated that linear and macrocyclic Ang IV-derived peptides can alter the shape and increase the number of dendritic spines in hippocampal cultures, properties associated with enhanced cognitive performance. After screening a library of 10 500 drug like substances for their ability to inhibit IRAP, we identified a series of low molecular weight aryl sulfonamides, which exhibit no structural similarity to Ang IV, as moderately potent IRAP inhibitors:A structural and biological characterization of three of these aryl sulfonamides was performed. Their binding modes to human IRAP were explored by docking calculations combined with molecular dynamics simulations and binding affinity estimations using the linear interaction energy method. Two alternative binding modes emerged from this analysis, both of which correctly rank the ligands according to their experimental binding affinities for this series of compounds. Finally, we show that two of these drug-like IRAP inhibitors can alter dendritic spine morphology and increase spine density in primary cultures of hippocampal neurons.