Interaction of nonsteroidal anti-inflammatory drugs with membranes: in vitro assessment and relevance for their biological actions.

Nonsteroidal anti-inflammatory drugs (NSAIDs) are among the most commonly used drugs in the world due to their anti-inflammatory, analgesic and antipyretic properties. Nevertheless, the consumption of these drugs is still associated with the occurrence of a wide spectrum of adverse effects. Regarding the major role of membranes in cellular events, the hypothesis that the biological actions of NSAIDs may be related to their effect at the membrane level has triggered the in vitro assessment of NSAIDs-membrane interactions. The use of membrane mimetic models, cell cultures, a wide range of experimental techniques and molecular dynamics simulations has been providing significant information about drugs partition and location within membranes and also about their effect on diverse membrane properties. These studies have indeed been providing evidences that the effect of NSAIDs at membrane level may be an additional mechanism of action and toxicity of NSAIDs. In fact, the pharmacokinetic properties of NSAIDs are closely related to the ability of these drugs to interact and overcome biological membranes. Moreover, the therapeutic actions of NSAIDs may also result from the indirect inhibition of cyclooxygenase due to the disturbing effect of NSAIDs on membrane properties. Furthermore, increasing evidences suggest that the disordering effects of these drugs on membranes may be in the basis of the NSAIDs-induced toxicity in diverse organ systems. Overall, the study of NSAIDs-membrane interactions has proved to be not only important for the better understanding of their pharmacological actions, but also for the rational development of new approaches to overcome NSAIDs adverse effects.

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

The discovery of selective cyclooxygenase-2 (COX-2) inhibitors represents a major achievement of the efforts over the past few decades to develop therapeutic treatments for inflammation. To gain insights into designing new COX-2-selective inhibitors, we address the energetic and structural basis for the selective inhibition of COX isozymes by means of a combined computational protocol involving docking experiment, force field design for the heme prothetic group, and free energy perturbation (FEP) simulation. We consider both COX-2- and COX-1-selective inhibitors taking the V523I mutant of COX-2 to be a relevant structural model for COX-1 as confirmed by a variety of experimental and theoretical evidences. For all COX-2-selective inhibitors under consideration, we find that free energies of binding become less favorable as the receptor changes from COX-2 to COX-1, due to the weakening and/or loss of hydrogen bond and hydrophobic interactions that stabilize the inhibitors in the COX-2 active site. On the other hand, COX-1-selective oxicam inhibitors gain extra stabilization energy with the change of residue 523 from valine to isoleucine because of the formations of new hydrogen bonds in the enzyme-inhibitor complexes. The utility of the combined computational approach, as a valuable tool for in silico screening of COX-2-selective inhibitors, is further exemplified by identifying the physicochemical origins of the enantiospecific selective inhibition of COX-2 by alpha-substituted indomethacin ethanolamide inhibitors.

A free energy based computational pathway from chemical templates to lead compounds: a case study of COX-2 inhibitors.

Automation of lead compound design in silico given the structure of the protein target and a definition of its active site vies for the top of the wish list in any drug discovery programme. We present here an enumeration of steps starting from chemical templates and propose a solution at the state of the art, in the form of a system independent comprehensive computational pathway. This methodology is illustrated with cyclooxygenase-2 (COX-2) as a target. We built candidate molecules including a few Non Steroidal Anti-inflammatory Drugs (NSAIDs) from chemical templates, passed them through empirical filters to assess drug-like properties, optimized their geometries, derived partial atomic charges via quantum calculations, performed Monte Carlo docking, carried out molecular mechanics and developed free energy estimates with Molecular Mechanics Generalized Born Solvent Accessibility (MMGBSA) methodology for each of the candidate molecules. For the case of aspirin, we also conducted molecular dynamics on the enzyme, the drug and the complex with explicit solvent followed by binding free energy analysis. Collectively, the results obtained from the above studies viz. sorting of drugs from non-drugs, semi-quantitative estimates of binding free energies, amply demonstrate the viability of the strategy proposed for lead selection/design for biomolecular targets.

Recognition of cyclooxygenase-2 (COX-2) active site by NSAIDs: a computer modelling study.

The energetics and models of COX-2 complexed with nonsteroidal anti-inflammatory drugs (NSAIDs) having different degrees of selectivity for two isoforms of COX (COX-2 and COX-1) have been studied using computer modelling approach. The models are obtained for complexes of NS398 (NS), a selective COX-2 inhibitor; indoprofen (Ind), a non-selective inhibitor; di-tert-butylbenzofurans (DHDMBFs) with substituents at the 5th position: CONH(CH2)2OMe (BF1), CONH-c-Pr (BF2), 3-methylene-gamma-butyrolactonyl (BF3) and oxicams namely, meloxicam (Mel), piroxicam (Pir) and tenoxicam (Ten). These were optimized using molecular mechanics (MM) and molecular dynamics (MD) techniques. The binding energies and structures were compared with pharmacological parameters and available results with COX-1. In case of NS a larger difference in the binding energies between COX-2 and COX-1 was noticed as compared to that of Ind. It also had stronger interaction with His90 and Tyr355 which is considered important for COX-2 selectivity. There was a difference in the compactness at the channel entrance between COX-2 selective and non-selective ligands. Models with DHDMBFs and oxicams showed a similar correlation. The results were used to design a peptide inhibitor, Tyr-Arg-Cys-Ala-delta Phe-Cys (Pept) which could fit better in the COX-2 cavity. As per our MD simulation results this peptide inhibitor showed both higher activity and COX-2 selectivity.

Rationale for the observed COX-2/COX-1 selectivity of celecoxib from Monte Carlo simulations.

Computational studies have yielded an analysis of the contributions to the free energy difference between the binding of celecoxib to COX-1 and to COX-2. The energetic and structural results point to the Ile to Val mutation at residue 523 as the key contributor to COX-2 selectivity; unfavorable steric contact between a sulfonamide oxygen and the delta methyl group of Ile523 destabilizes the complex with COX-1. The His to Arg change at residue 513 is less significant.