Three and four dimensional-quantitative structure activity relationship (3D/4D-QSAR) analyses of CYP2D6 inhibitors.

Three- and four-dimensional quantitative structure activity relationship (3D/4D-QSAR) pharmacophore models of competitive inhibitors of CYP2D6 were constructed using data from our laboratory or the literature. The 3D-QSAR pharmacophore models of the common structural features of CYP2D6 inhibitors were built using the program Catalyst (Molecular Simulations, San Diego, CA, USA). These 3D-QSAR models were compared with 3D and 4D-QSAR partial least squares (PLS) models which were constructed using molecular surface-weighted holistic invariant molecular (MS-WHIM) descriptors of size and shape of inhibitors. The first Catalyst model was generated from multiple conformers of competitive inhibitors (n = 20) of CYP2D6 mediated bufurolol 1'-hydroxylation. This model demonstrated a correlation of observed and predicted Ki (apparent) values of r = 0.75. A second Catalyst model was constructed from literature derived Ki (apparent) values (n = 31) for the inhibition of CYP2D6. This model provided a correlation of observed and predicted inhibition for CYP2D6 of r = 0.91. Both Catalyst Ki pharmacophores were then validated by predicting the Ki (apparent) of a test set of known CYP2D6 inhibitors (n = 15). Ten out of 15 of these Ki (apparent) values were predicted to be within one log residual of the observed value using our CYP2D6 inhibitor model, while the literature model predicted nine out of 15 values. Similarly, 3D- and 4D-QSARs derived from PLS MS-WHIM for our dataset yielded predictable models as assessed using cross-validation. The corresponding cross-validated PLS MS-WHIM model for the literature dataset yielded a comparable 3D-QSAR and improved 4D-QSAR value. Such computational models will aid in future prediction of drug-drug interactions.

Development of dual-acting agents for thromboxane receptor antagonism and thromboxane synthase inhibition. 3. Synthesis and biological activities of oxazolecarboxamide-substituted omega-phenyl-omega-(3-pyridyl)alkenoic acid derivatives and related compounds.

A novel series of oxazolecarboxamide-substituted omega-phenyl-omega-(3-pyridyl)alkenoic acid derivatives was discovered as potent dual-acting agents to block the TXA2 receptor and to inhibit the thromboxane synthase (TRA/TSI). Synthesis, structure-activity relationship (SAR), and in vitro and in vivo pharmacology of this series of compounds are described. Modification of the series revolved around the oxazole moiety to increase the hydrophilicity of the compounds and to correlate the biological activity with lipophilicity of the compounds. The most potent in the series was (E)-7-[4-[4-[[(4-cyclohexylbutyl)amino]carbonyl]-2-oxazolyl] phenyl]-7 -(3-pyridyl)hept-6-enoic acid (14) with Kd = 9.9 +/- 0.4 nM for the thromboxane receptor antagonism and IC50 = 55.0 +/- 17.9 nM for thromboxane synthase inhibition. The compound 14 was a selective TRA/TSI which exhibited desirable characteristics for oral activity, "shunt" effect to elevate PGI2 level, and absence of agonist activity.

Synthesis of syn and anti isomers of 6-[[(hydroxyimino)phenyl]methyl]-1-[(1-methylethyl)sulfonyl]-1H-benzimidazol-2-amine. Inhibitors of rhinovirus multiplication.

The synthesis and antirhinovirus activity of syn and anti isomers of 6-[[(hydroxyimino)phenyl]methyl]-1-[(1-methylethyl)sulfonyl]-1H-benzimidazol-2-amine (4 and 5) are reported. The structural assignments of 4 and 5 are based upon 13C NMR spectra of both isomers and also X-ray analysis of 5. The anti-isomer 5 was more potent than the syn-isomer 4 when compared as an inhibitor of rhinovirus multiplication in vitro. Both isomers inhibited multiplication of 15 different serotypes of rhinovirus.

Use of conformationally restricted benzamidines as arginine surrogates in the design of platelet GPIIb-IIIa receptor antagonists.

The use of 5,6-bicyclic amidines as arginine surrogates in the design of a novel class of potent platelet glycoprotein IIb-IIIa receptor (GPIIb-IIIa) antagonists is described. The additional conformational restriction offered by the bicyclic nucleus results in 20-400-fold increases in potency compared to the freely flexible, acyclic benzamidine counterpart. The design, synthesis, structure-activity relationships (SAR), and in vitro activity of this novel class of GPIIb-IIIa antagonists are presented.

1,2-Dibenzamidobenzene inhibitors of human factor Xa.

High-throughput screening of a combinatorial library of diamidophenols yielded lead compounds with the ability to inhibit human factor Xa (fXa) at micromolar concentrations (e.g. compound 4, fXa apparent K(ass) = 0.64 x 10(6) L/mol). SAR studies in this novel structural series of fXa inhibitors showed that the phenolic hydroxyl group was not essential for activity. The best activity was found in substituted 1,2-dibenzamidobenzenes in which the phenyl group of one benzoyl group (A-ring) was substituted in the 4-position with relatively small lipophilic or polarizable groups such as methoxy, vinyl, or chloro and the phenyl group of the other benzoyl group (B-ring) was substituted in the 4-position with larger lipophilic groups such as tert-butyl or dimethylamino. The central phenyl ring (C-ring) tolerated a wide variety of substituents, but methoxy, methanesulfonamido, hydroxyl, and carboxyl substitution produced slightly higher levels of activity than other substituents when present in combination with favorable B-ring substitution. Methylation of the amide nitrogen atoms was found to greatly decrease activity. Compound 12 is the highest affinity fXa inhibitor in this group of compounds, having fXa apparent K(ass) = 25.5 x 10(6) L/mol, about 40x more active than the original lead. This lead series does not show potent inhibition of human thrombin. A model for the binding of these ligands to the fXa active site is proposed. The model is consistent with the observed SAR and can serve to guide future SAR studies.

Depletion of cardiac norepinephrine in rats and mice by 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP).

1-Methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) is a commercially available chemical reagent. Although little has been known about its biological effects, recently MPTP has been reported to cause irreversible Parkinson's disease-like symptoms in humans and in monkeys. We describe here another pharmacologic effect of MPTP, the ability to deplete cardiac norepinephrine in rats and mice. In mice, cardiac norepinephrine concentration decreased within 1 hr, was maximally depleted at 24 hr, and recovered by 4-7 days after i.p. injection of a 32 mg/kg dose of MPTP. The depletion was antagonized by desipramine pretreatment, as was norepinephrine depletion by tyramine. In rats, cardiac norepinephrine depletion by 10-30 mg/kg, i.p., doses of MPTP was accompanied by depletion of cardiac dopamine and of norepinephrine in the mesenteric artery. In rats and in mice, norepinephrine in brain was affected to a smaller degree than was norepinephrine in heart, and dopamine in brain was depleted very little if at all. In spontaneously hypertensive rats, the depletion of cardiac norepinephrine was associated with a marked antihypertensive effect. The p-hydroxy analog of MPTP did not deplete cardiac norepinephrine in rats, indicating that its possible formation as a metabolite of MPTP was not involved in the depletion of cardiac norepinephrine. These findings extend the spectrum of known pharmacologic effects of MPTP.

Progress in predicting human ADME parameters in silico.

Understanding the development of a scientific approach is a valuable exercise in gauging the potential directions the process could take in the future. The relatively short history of applying computational methods to absorption, distribution, metabolism and excretion (ADME) can be split into defined periods. The first began in the 1960s and continued through the 1970s with the work of Corwin Hansch et al. Their models utilized small sets of in vivo ADME data. The second era from the 1980s through 1990s witnessed the widespread incorporation of in vitro approaches as surrogates of in vivo ADME studies. These approaches fostered the initiation and increase in interpretable computational ADME models available in the literature. The third era is the present were there are many literature data sets derived from in vitro data for absorption, drug-drug interactions (DDI), drug transporters and efflux pumps [P-glycoprotein (P-gp), MRP], intrinsic clearance and brain penetration, which can theoretically be used to predict the situation in vivo in humans. Combinatorial synthesis, high throughput screening and computational approaches have emerged as a result of continual pressure on pharmaceutical companies to accelerate drug discovery while decreasing drug development costs. The goal has become to reduce the drop-out rate of drug candidates in the latter, most expensive stages of drug development. This is accomplished by increasing the failure rate of candidate compounds in the preclinical stages and increasing the speed of nomination of likely clinical candidates. The industry now understands the reasons for clinical failure other than efficacy are mainly related to pharmacokinetics and toxicity. The late 1990s saw significant company investment in ADME and drug safety departments to assess properties such as metabolic stability, cytochrome P-450 inhibition, absorption and genotoxicity earlier in the drug discovery paradigm. The next logical step in this process is the evaluation of higher throughput data to determine if computational (in silico) models can be constructed and validated from it. Such models would allow an exponential increase in the number of compounds screened virtually for ADME parameters. A number of researchers have started to utilize in silico, in vitro and in vivo approaches in parallel to address intestinal permeability and cytochrome P-450-mediated DDI. This review will assess how computational approaches for ADME parameters have evolved and how they are likely to progress.

Three-dimensional-quantitative structure activity relationship analysis of cytochrome P-450 3A4 substrates.

To gain a better understanding of the active site of cytochrome P-450 (CYP) 3A4, a three-dimensional-quantitative structure activity relationship model was constructed using the structures and K(m (apparent)) values of 38 substrates of human liver microsomal CYP3A4. This pharmacophore was built using the program Catalyst and consisted of four features: two hydrogen bond acceptors, one hydrogen bond donor, and one hydrophobic region. The pharmacophore demonstrated a fit value (r) of observed and expected K(m(apparent)) value of 0.67. The validity of the CYP3A4 substrate model was tested by twice permuting (randomizing) the activity values and substrate structures. The results of this validation procedure indicated that the original model was a significant representation of the features required of CYP3A4 substrates. The second validation method used the Catalyst model to predict the K(m(apparent)) values of a test set of structurally diverse substrates for CYP3A4 not included in the 38 molecules used to build the model. Two fitting algorithms included in this software were examined: fast fit and best fit. The fast fitting method resulted in predictions for all 12 substrates that were within 1 log unit for the residual [i.e., the difference between predicted and observed K(m(apparent))]. In contrast, the best fit algorithm poorly predicted the K(m (apparent)) values (i.e., residual >1 log unit) of 4 of 12 substrates. These poor fits with the best fit function suggest that the fast fit method within Catalyst is more representative of the observed K(m(apparent)) values for CYP3A4 substrates and enables good in silico prediction of this activity. A Catalyst common features pharmacophore was also constructed from three molecules known to activate their own metabolism included in the 38 molecules of the initial CYP3A4 model. This demonstrated that activators of CYP3A4 possess multiple hydrophobic regions that might correspond with a region in the active site away from the metabolic site.

Three- and four-dimensional quantitative structure activity relationship analyses of cytochrome P-450 3A4 inhibitors.

The program Catalyst was used to build three-dimensional quantitative structure activity relationship (3D-QSAR) pharmacophore models of the structural features common to competitive-type inhibitors of cytochrome P-450 (CYP) 3A4. These were compared with 3D- and four-dimensional (4D)-QSAR partial least-squares (PLS) models built using molecular surface-weighted holistic invariant molecular (MS-WHIM) descriptors for size and shape of the inhibitor. The Catalyst pharmacophore model generated from multiple conformers of competitive inhibitors of CYP3A4-mediated midazolam 1'-hydroxylation (n = 14) yielded a high correlation of observed and predicted Ki values of r = 0.91. Similarly, PLS MS-WHIM was used to produce 3D- and 4D-QSARs for this data set and produced models that were statistically predictable after cross-validation. Two additional Catalyst pharmacophores were constructed from literature Ki values (n = 32) derived from the inhibition of CYP3A-mediated cyclosporin A metabolism and IC50 data (n = 22) from the inhibition of CYP3A4-mediated quinine 3-hydroxylation. These Catalyst pharmacophores illustrated correlations of observed and predicted inhibition for CYP3A4 of r = 0.77 and 0.92, respectively. The corresponding 4D-QSARs generated by PLS MS-WHIM for these data sets were of comparable quality as judged by cross-validation. Both Ki pharmacophores generated with Catalyst were also validated by predicting the Ki(apparent) values of a test set of eight CYP3A4 inhibitors not included in either model. In seven of eight cases, the residuals of the predicted Ki(apparent) values were within 1 log unit of the observed values. The 3D- and 4D-QSAR models produced in this study suggest the utility of future in silico prediction of CYP3A4-mediated drug-drug interactions.

Three- and four-dimensional-quantitative structure activity relationship (3D/4D-QSAR) analyses of CYP2C9 inhibitors.

The interaction of competitive type inhibitors with the active site of cytochrome P450 (CYP) 2C9 has been predicted using three- and four-dimensional quantitative structure activity relationship (3D-/4D-QSAR) models constructed using previously unreported and literature-derived data. 3D-QSAR pharmacophore models of the common structural features of CYP2C9 inhibitors were built using the program Catalyst and compared with 3D- and 4D-QSAR partial least-squares models, which use molecular surface-weighted holistic invariant molecular descriptors of the size and shape of inhibitors. The Catalyst models generated from multiple conformers of competitive inhibitors of CYP2C9 activities contained at least one hydrophobic and two hydrogen bond acceptor/donor regions. Catalyst model 1 was constructed with Ki(apparent) values for inhibitors of tolbutamide and diclofenac 4'-hydroxylation (n = 9). Catalyst model 2 was generated from literature Ki(apparent) values for (S)-warfarin 7-hydroxylation (n = 29), and Catalyst model 3 from literature IC50 values for tolbutamide 4-hydroxylation (n = 13). These three models illustrated correlation values of observed and predicted inhibition for CYP2C9 of r = 0.91, 0.89, and 0.71, respectively. Catalyst pharmacophores generated with Ki(apparent) values were validated by predicting the Ki(apparent) value of a test set of CYP2C9 inhibitors also derived from the literature (n = 14). Twelve of fourteen of these Ki(apparent) values were predicted to be within 1 log residual of the observed value using Catalyst model 1, whereas Catalyst model 2 predicted 10 of 14 Ki(apparent) values. The corresponding partial least-squares molecular surface-weighted holistic invariant molecular 3D- and 4D-QSAR models for all CYP2C9 data sets yielded predictable models as assessed using cross-validation. These 3D- and 4D-QSAR models of CYP inhibition will aid in future prediction of drug-drug interactions.

LY249933: a cardioselective 1,4-dihydropyridine with positive inotropic activity.

Compound LY249933 and its component diastereomers, (RR) and (SR), were studied for their vascular and cardiac effects in vitro and in vivo. In guinea pig cardiac ventricular membranes, LY249933, (RR), and (SR) potently displaced bound [3H]nitrendipine (Kd values = 2-6 nM). In isolated guinea pig right ventricular strips, LY249933 produced a small but significant increase in contraction, whereas (RR) substantially increased (-log EC50 (M) = 4.6 +/- 0.8) and (SR) decreased contraction (-log EC50 (M) = 4.1 +/- 0.8). In isolated canine cephalic vein, contracted with 80 mM KCl, an increase in contraction was produced by (RR), whereas relaxation was produced by LY249933 (-log EC50 (M) = 5.9 +/- 0.9) and (SR) (-log EC50 (M) = 6.0 +/- 0.7). At 20 mM KCl, (RR) increased, (SR) decreased, but LY249933 did not alter contraction. In anesthetized dogs, LY249933 (200 micrograms/kg/min, i.v.) increased dP/dt60, decreased heart rate, but did not change vascular resistance or rate pressure product. At the same dose, (RR) and (SR) both tended to increase dP/dt60 nonsignificantly, whereas (RR) increased and (SR) decreased vascular resistance. Both (RR) and (SR) tended to decrease heart rate nonsignificantly, whereas (RR) did not change and (SR) decreased rate pressure product. Thus, LY249933 produced potentially beneficial cardiovascular changes resulting from the combined actions of its (RR) and (SR) diastereomers that are postulated to be calcium agonist and antagonist, respectively.

Comparative QSAR studies of two series of 1,4-dihydropyridines as slow calcium channel blockers.

Quantitative structure activity analysis was applied to two series of dihydropyridine (DHP) calcium channel blocking agents. One series of compounds was composed of DHPs substituted in the 4-position with an ortho or meta nitro substituted phenyl ring. The second group consisted of DHPs substituted at the 4-position with a novel thieno [3,2-c] pyridine ring. Both series consisted of compounds with unsymmetrical ester substitutions on the dihydropyridine ring. The antihypertensive activity of the compounds were determined in a spontaneously hypertensive rat model. Regression analysis indicated the antihypertensive activity of an i.v. dose correlated with the calculated octanol/water coefficent (clogP). Regression analysis did not find correlation with the in vitro potency and the clogP values.

Three-dimensional quantitative structure-permeability relationship analysis for a series of inhibitors of rhinovirus replication.

Multiple three-dimensional quantitative structure-activity relationship (3D-QSAR) approaches were applied to predicting passive Caco-2 permeability for a series of 28 inhibitors of rhinovirus replication. Catalyst, genetic function approximation (GFA) with MS-WHIM descriptors, CoMFA, and VolSurf were all used for generating 3D-quantitative structure permeability relationships utilizing a training set of 19 molecules. Each of these approaches was then compared using a test set of nine molecules not present in the training set. Statistical parameters for the test set predictions (r(2) and leave-one-out q(2)) were used to compare the models. It was found that the Catalyst pharmacophore model was the most predictive (test set of predicted versus observed permeability, r(2) = 0.94). This model consisted of a hydrogen bond acceptor, hydrogen bond donor, and ring aromatic feature with a training set correlation of r(2) = 0.83. The CoMFA model consisted of three components with an r(2) value of 0.96 and produced good predictions for the test set (r(2) = 0.84). VolSurf resulted in an r(2) value of 0.76 and good predictions for the test set (r(2) = 0.83). Test set predictions with GFA/WHIM descriptors (r(2) = 0.46) were inferior when compared with the Catalyst, CoMFA, and VolSurf model predictions in this evaluation. In summary it would appear that the 3D techniques have considerable value in predicting passive permeability for a congeneric series of molecules, representing a valuable asset for drug discovery.

Three-dimensional quantitative structure activity relationship analyses of substrates for CYP2B6.

To begin to build an understanding of the interactions of CYP2B6 with substrates, two different 3-dimensional quantitative structure activity relationship (3D-QSAR) models were constructed using 16 substrates of B-lymphoblastoid expressed CYP2B6. A pharmacophore model was built using the program Catalyst, which was compared with a partial least-squares (PLS) model using molecular surface-weighted holistic invariant molecular (MS-WHIM) descriptors. The Catalyst model yielded a 3-dimensional model of the common structural features of CYP2B6 substrates, whereas PLS MS-WHIM generated a model based on statistical analyses of molecular descriptors for size and shape of the substrate. The pharmacophore model obtained with Catalyst consisted of three hydrophobes and one hydrogen bond acceptor region. The cross-validated PLS MS-WHIM model gave a good q2 value of 0.607. Size, positive electrostatic potential, hydrogen bonding acceptor capacity, and hydrophobicity were found to be the most relevant descriptors for the model. These models were then used to predict the Km (apparent) values of a test set of structurally diverse substrates for CYP2B6 not included in the model building, specifically lidocaine, amitriptyline, bupropion, arteether, and verapamil. Overall, both 3D-QSAR methods yielded satisfactory Km (apparent) value predictions for the majority of the molecules in the test set. However, PLS MS-WHIM was unable to reliably predict the Km (apparent) value for verapamil, whereas Catalyst did not predict the Km (apparent) value for lidocaine. In both of these cases the residual of the Km (apparent) value was greater than one log unit. The strengths and limitations of both of these 3D-QSAR approaches are discussed.

Diamino benzo[b]thiophene derivatives as a novel class of active site directed thrombin inhibitors: 3. Enhancing activity by imposing conformational restriction in the C-4" side chain.

The preparation and biological evaluation of a series of benzo[b]thiophene diamine thrombin inhibitors possessing conformationally restricted C-4" linkers are reported. Compared to the parent compounds 1a/b, the unsaturated derivatives 3a/b exhibited a modest twofold increase in thrombin inhibitory activity, while the more lipophilic carbocyclic ring containing analogs 4a/b affected an eightfold enhancement in potency.

Dibasic benzo[b]thiophene derivatives as a novel class of active site directed thrombin inhibitors. 2. Exploring interactions at the proximal (S2) binding site.

In an effort to increase the thrombin inhibitory activity of a novel series of inhibitors (i.e., 1a), substituents were incorporated at the C-3" position of the C-3 aryl ring (2). Consistent with the X-ray crystallography studies, small hydrophobic groups at the C-3" site (Br and Me) enhanced thrombin inhibitory activity by 8-fold. However, a few more hydrophilic substituents (NO2 and OMe) also enhanced the potency of the series. The biological results are discussed in terms of molecular modeling studies.