Protein crystallography and examples of its applications in medicinal chemistry

The field of protein crystallography inspires and enthrals, whether it be for the beauty and symmetry of a perfectly formed protein crystal, the unlocked secrets of a novel protein fold, or the precise atomic-level detail yielded from a protein-ligand complex. Since 1958, when the first protein structure was solved, there have been tremendous advances in all aspects of protein crystallography, from protein preparation and crystallisation through to diffraction data measurement and structure refinement. These advances have significantly reduced the time required to solve protein crystal structures, while at the same time substantially improving the quality and resolution of the resulting structures. Moreover, the technological developments have induced researchers to tackle ever more complex systems, including ribosomes and intact membrane-bound proteins, with a reasonable expectation of success. In this review, the steps involved in determining a protein crystal structure are described and the impact of recent methodological advances identified. Protein crystal structures have proved to be extraordinarily useful in medicinal chemistry research, particularly with respect to inhibitor design. The precise interaction between a drug and its receptor can be visualised at the molecular level using protein crystal structures, and this information then used to improve the complementarity and thus increase the potency and selectivity of an inhibitor. The use of protein crystal structures in receptor-based drug design is highlighted by (i) HIV protease, (ii) influenza virus neuraminidase and (iii) prostaglandin H-2-synthetase. These represent, respectively, examples of protein crystal structures that (i) influenced the design of drugs currently approved for use in the treatment of HIV infection, (ii) led to the design of compounds currently in clinical trials for the treatment of influenza infection and (iii) could enable the design of highly specific non-steroidal anti-inflammatory drugs that lack the common side-effects of this drug class.

Hepatic disposition and metabolite kinetics of a homologous series of diflunisal esters

The hepatic disposition and metabolite kinetics of a homologous series of diflunisal O-acyl esters (acetyl, butanoyl, pentanoyl, anti hexanoyl) were determined using a single-pass perfused in situ rat liver preparation. The experiments were conducted using 2% BSA Krebs-Henseleit buffer (pH 7.4), and perfusions were performed at 30 mL/min in each liver. O-Acyl esters of diflunisal and pregenerated diflunisal were injected separately into the portal vein. The venous outflow samples containing the esters and metabolite diflunisal were analyzed by high performance liquid chromatography (HPLC). The normalized outflow concentration-time profiles for each parent ester and the formed metabolite, diflunisal, were analyzed using statistical moments analysis and the two-compartment dispersion model. Data (presented as mean +/- standard error for triplicate experiments) was compared using ANOVA repeated measures, significance level P < 0.05. The hepatic availability (AUC'), the fraction of the injected dose recovered in the outflowing perfusate, for O-acetyldiflunisal (C2D = 0.21 +/- 0.03) was significantly lower than the other esters (0.34-0.38). However, R-N/f(u), the removal efficiency number R-N divided by the unbound fraction in perfusate f(u), which represents the removal efficiency of unbound ester by the liver, was significantly higher for the most lipophilic ester (O-hexanoyldiflunisal, C6D = 16.50 +/- 0.22) compared to the other members of the series (9.57 to 11.17). The most lipophilic ester, C6D, had the largest permeability surface area (PS) product (94.52 +/- 38.20 mt min-l g-l liver) and tissue distribution value VT (35.62 +/- 11.33 mL g(-1) liver) in this series. The MTT of these O-acyl esters of diflunisal were not significantly different from one another. However, the metabolite diflunisal MTTs tended to increase with the increase in the parent ester lipophilicity (11.41 +/- 2.19 s for C2D to 38.63 +/- 9.81 s for C6D). The two-compartment dispersion model equations adequately described the outflow profiles for the parent esters and the metabolite diflunisal formed from the O-acyl esters of diflunisal in the liver.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

Oxidation of oxa and thia fatty acids and related compounds catalysed by 5-and 15-lipoxygenase

The modified fatty acids, (Z,Z,Z)-(octadeca-6,9,12-trienyloxy)acetic acid, (Z,Z,Z)-(octadeca-9,12,15-trienyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)acetic acid, 3-[(all-Z)-(eicosa-5,8,11,14-tetraenylthio)]propionic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)succinic acid, N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]glycine and N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]aspartic acid, all react with soybean 15-lipoxygenase. The products were treated with triphenylphosphine to give alcohols, which were isolated using HPLC. Analysis of the alcohols using negative ion tandem electrospray mass spectrometry, and by comparison with compounds obtained by autoxidation of arachidonic acid, shows that each enzyme catalysed oxidation occurs at the omega -6 position of the substrate. In a similar fashion, it has been found that (Z,Z,Z)-(octadeca-6,9,12-trienyloxy)acetic acid, (Z,Z,Z)-(octadeca-9,12,15-trienyloxy)acetic acid, (all-Z)-(eicosa-5,8,11,14-tetraenylthio)acetic acid and N-[(all-Z)-(eicosa-5,8, 11.14-tetraenylthio)]propionic acid each undergoes regioselective oxidation at the carboxyl end of the polyene moiety on treatment with potato 5-lipoxygenase. Neither (all-Z)-(eicosa-5,8,11,14-tetraenylthio)succinic acid nor N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]aspartic acid reacts in the presence of this enzyme, while N-[(all-Z)-(eicosa-5,8,11,14-tetraenoyl)]glycine affords the C11' oxidation product. The alcohol derived from (Z,Z,Z)-(octadeca-6,9, 12-trienyloxy)acetic acid using the 15-lipoxygenase reacts at the C6' position with the 5-lipoxygenase. (C) 2001 Elsevier Science Ltd. All rights reserved.

Percutaneous absorption of sunscreen agents from liquid paraffin: Self-association of octyl salicylate and effects on skin flux

This study provides an investigation of the availability of octyl salicylate (OS), a common sunscreen agent, from liquid paraffin and the effect of OS on skin permeability. A model membrane system to isolate the vehicle effect from membrane permeability has been developed. Partitioning of OS between liquid paraffin and aqueous receptor phases was conducted. Partition coefficients increased with increase in OS concentration. A range of OS concentrations in liquid paraffin was diffused across human epidermis and synthetic membranes into 4% bovine serum albumin in phosphate-buffered saline and 50% ethanol. Absorption profiles of OS obtained from silicone and low-density polyethylene (LDPE) membranes were similar to each other but higher than for the high-density polyethylene [HDPE (3 times)] membrane and human epidermis (15 times). The steady state fluxes and apparent permeability coefficients (K-p') obtained from the diffusion studies showed the same trends with all membranes, except for the HDPE membrane which showed greater increase in flux and K-p' at concentrations above 30%. IR spectra showed that several bands of OS were shifted with concentrations, and the molecular models further suggested that the main contribution to the self-association is from non-1,4 van der Waals interactions.

Bioactivation of Phenytoin by Human Cytochrome P450: Characterization of the Mechanism and Targets of Covalent Adduct Formation

The cytochrome P450-dependent covalent binding of radiolabel derived fi om phenytoin (DPH) and its phenol and catechol metabolites, 5-(4'-hydroxyphenyl)-5-phenylhydantoin (HPPH) and 5-(3',4'-dihydroxyphenyl)-5-phenylhydantoin (CAT), was examined in liver microsomes. Radiolabeled HPPH and CAT and unlabeled CAT were obtained from microsomal incubations and isolated by preparative HPLC. NADPH-dependent covalent binding was demonstrated in incubations of human liver microsomes with HPPH. When CAT was used as substrate, covalent adduct formation was independent of NADPH, was enhanced in the presence of systems generating reactive oxygen species, and was diminished under anaerobic conditions or in the presence of cytoprotective reducing agents. Fluorographic analysis showed that radiolabel derived from DPH and HPPH was selectively associated with proteins migrating with approximate relative molecular weights of 57-59 kDa and at the dye front (molecular weights < 23 kDa) on denaturing gels. Lower levels of radiolabel were distributed throughout the molecular weight range. In contrast, little selectivity was seen in covalent adducts formed from CAT. HPPH was shown to be a mechanism-based inactivator of P450, supporting the contention that a cytochrome P450 is one target of covalent binding. These results suggest that covalent binding of radiolabel derived from DPH in rat and human Liver microsomes occurs via initial P450-dependent catechol formation followed by spontaneous oxidation to quinone and semiquinone derivatives that ultimately react with microsomal protein. Targets for covalent binding may include P450s, though the catechol appears to be sufficiently stable to migrate out of the P450 active site to form adducts with other proteins. In conclusion, we have demonstrated that DPH can be bioactivated in human liver to metabolites capable of covalently binding to proteins. The relationship of adduct formation to DPH-induced hypersensitivity reactions remains to be clarified.

Photodegradation of irinotecan (CPT-11) in aqueous solutions: Identification of fluorescent products and influence of solution composition

The photodegradation of irinotecan (CPT-11), the semisynthetic derivative of the antitumor alkaloid 20(S)-camptothecin, has been investigated. The drug was exposed to laboratory light for up to 5 days in 0.9% saline solution (pH 8.5). Five significant photodegradation products were observed and a high-performance liquid chromatography (HPLC) assay was employed to isolate them from CPT-11 using gradient conditions. The structures were elucidated by nuclear magnetic resonance spectroscopy and tandem mass spectrometry and shown to be the result of extensive modifications of the lactone ring of CPT-11. Three of the compounds were found to belong to the mappicine group of alkaloids. In addition, the effect of light on the stability of CPT-11 in aqueous solutions and biological fluids was also assessed, Potassium phosphate buffers (0.05 M, pH 5.0-8.2) and saline, plasma, urine, and bile solutions containing 20 mu M CPT-11 were equilibrated in the dark for 24 h before being exposed to laboratory light for up to 171 h at ambient temperature. Four of the five identified photodegradation products were observed and quantitated by isocratic HPLC, using a different detection mode (fluorescence) than the one used for gradient elution, In general, CPT-11 was found to be unstable under neutral and alkaline conditions for all solutions investigated, with the exception of bile. We conclude that CPT-11 is photolabile and that care should be taken to protect samples, particularly those intended for the isolation and identification of novel metabolites of CPT-11.

Diffusion modeling of percutaneous absorption kinetics: 2. Finite vehicle volume and solvent deposited solids

The diffusion model for percutaneous absorption is developed for the specific case of delivery to the skin being limited by the application of a finite amount of solute. Two cases are considered; in the first, there is an application of a finite donor (vehicle) volume, and in the second, there are solvent-deposited solids and a thin vehicle with a high partition coefficient. In both cases, the potential effect of an interfacial resistance at the stratum corneum surface is also considered. As in the previous paper, which was concerned with the application of a constant donor concentration, clearance limitations due to the viable eqidermis, the in vitro sampling rate, or perfusion rate in vivo are included. Numerical inversion of the Laplace domain solutions was used for simulations of solute flux and cumulative amount absorbed and to model specific examples of percutaneous absorption of solvent-deposited solids. It was concluded that numerical inversions of the Laplace domain solutions for a diffusion model of the percutaneous absorption, using standard scientific software (such as SCIENTIST, MicroMath Scientific software) on modern personal computers, is a practical alternative to computation of infinite series solutions. Limits of the Laplace domain solutions were used to define the moments of the flux-time profiles for finite donor volumes and the slope of the terminal log flux-time profile. The mean transit time could be related to the diffusion time through stratum corneum, viable epidermal, and donor diffusion layer resistances and clearance from the receptor phase. Approximate expressions for the time to reach maximum flux (peak time) and maximum flux were also derived. The model was then validated using reported amount-time and flux-time profiles for finite doses applied to the skin. It was concluded that for very small donor phase volume or for very large stratum corneum-vehicle partitioning coefficients (e.g., for solvent deposited solids), the flux and amount of solute absorbed are affected by receptor conditions to a lesser extent than is obvious for a constant donor constant donor concentrations. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:504-520, 2001.

Elucidation of reaction scheme describing malondialdehyde-acetaldehyde-protein adduct formation

Malondialdehyde and acetaldehyde react together with proteins and form hybrid protein conjugates designated as MAA adducts, which have been detected in livers of ethanol-fed animals. Our previous studies have shown that MAA adducts are comprised of two distinct products. One adduct is composed of two molecules of malondialdehyde and one molecule of acetaldehyde and was identified as the 4-methpl-1,4-dihydropyridine-3,5-dicarbaldehyde derivative of an amino group (MHHDC adduct). The other adduct is a 1:1 adduct of malondialdehyde and acetaldehyde and was identified as the 2-formyl-3-(alkylamino)butanal derivative of an amino group (FAAB adduct). In this study, information on the mechanism of MAA adduct formation was obtained, focusing on whether the FAAB adduct serves as a precursor for the MDHDC adduct. Upon the basis of chemical analysis and NMR spectroscopy, two initial reaction steps appear to be a prerequisite for MDHDC formation. One step involves the reaction of one molecule of malondialdehyde and one of acetaldehyde with an amino group of a protein to form the FAAB product, while the other step involves the generation of a malondialdehyde-enamine. It appears that generation of the MDHDC adduct requires the FAAB moiety to be transferred to the nitrogen of the MDA-enamine. For efficient reaction of FAAB with the enamine to take place, additional experiments indicated that these two intermediates likely must be in positions on the protein of close proximity to each other. Further studies showed that the incubation of liver proteins from ethanol-fed rats with MDA resulted in a marked generation of MDHDC adducts, indicating the presence of a pool of FAAB adducts in the liver of ethanol-fed animals. Overall, these findings show that MDHDC-protein adduct formation occurs via the reaction of the FAAB moiety with a malondialdehyde-enamine, and further suggest that a similar mechanism may be operative in vivo in the liver during prolonged ethanol consumption.

Ring flexibility within tricyclic antidepressant drugs

The internal flexibility of the central seven-membered ring of a series of tricyclic antidepressant drugs (TCAs), imipramine {l}, amitriptyline {2}, doxepin {3}, and dothiepin {4}, has been investigated by H-1 and C-13 nuclear magnetic (NMR) techniques. Two dynamic processes were examined: ring inversion and bridge flexing. H-1 NMR lineshape analysis was used to obtain ring inversion barriers for 2-4. These studies yielded energy barriers of 14.3, 16.7, and 15.7 +/- 0.6 kcal/mol for the hydrochloride salts of doxepin, dothiepin, and amitriptyline, respectively. The barriers for the corresponding free bases were lower by 0.6 kcal/mol on average. (CT1)-C-13 relaxation measurements were used to determine the degree of bridge flexing associated with the central seven-membered ring for all four compounds. By fitting the T-1 data to a two-state jump model, lifetimes and amplitudes of rapid bridge flexing motions were determined. The results show that imipramine has the fastest rate of bridge flexing, followed by amitriptyline, doxepin, and dothiepin. The pharmacological profiles of the TCAs are complex and they interact with many receptor sites, resulting in numerous side effects and a general lack of understanding of their precise mode of action in different anxiety-related disorders. They all have similar three-dimensional structures, which makes it difficult to rationalize their differing relative potency in different assays/clinical settings. However, the clear finding here that there are significantly different degrees of internal mobility suggests that molecular dynamics should be an additional factor considered when trying to understand the mode of action of this clinically important family of molecules. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:713-721, 2001.

Oxidizing species in the mechanism of cytochrome P450

This review discusses the mechanisms of oxygen activation by cytochrome P450 enzymes, the possible catalytic roles of the various iron-oxygen species formed in the catalytic cycle, and progress in understanding the mechanisms of hydrocarbon hydroxylation, heteroatom oxidation, and olefin epoxidation. The focus of the review is on recent results, but earlier work is discussed as appropriate. The literature through to February 2002 is surveyed, and 175 referenced are cited.

Syntheses of polycationic dendrimers on lipophilic peptide core for complexation and transport of oligonucleotides

Synthesis of novel polycationic lipophilic peptide core(s) was accomplished and these agents successfully transfected human retinal pigment epithelium cells with ODN1 upon complexation with the oligonucleotide. The level of transfection was indirectly measured by the decreased production of the protein hVEGF (human vascular endothelial growth factor) in comparison to the transfection agent cytofectin GSV(TM). (C) 2002 Elsevier Science Ltd. All rights reserved.

Bioactivation of tamoxifen by recombinant human cytochrome P450 enzymes

Tamoxifen is a major drug used for adjuvant chemotherapy of breast cancer; however, its use has been associated with a small but significant increase in risk of endometrial cancer. In rats, tamoxifen is a hepatocarcinogen, and DNA adducts have been observed in both rat and human tissues. Tamoxifen has been shown previously to be metabolized to reactive products that have the potential to form protein and DNA adducts. Previous studies have suggested a role for P450 3A4 in protein adduct formation in human liver microsomes, via a catechol intermediate; however, no clear correlation was seen between P450 3A4 content of human liver microsomes and adduct formation. In the present study, we investigated the P450 forms responsible for covalent drug-protein adduct formation and the possibility that covalent adduct formation might occur via alternative pathways to catechol formation. Recombinant P450 3A4 catalyzed adduct formation, and this correlated with the level of uncoupling in the P450 incubation, consistent with a role of reactive oxygen species in potentiating adduct formation after enzymatic formation of the catechol metabolite. Whereas P450s 1AI, 2D6, and 3A5 generated catechol metabolite, no covalent adduct formation was observed with these forms. By contrast, P450 2136, 2C19, and rat liver microsomes catalyzed drug-protein adduct formation but not catechol formation. Drug protein adducts formed specifically with P450 3A4 in incubations using membranes isolated from bacteria expressing P450 3A4 and reductase, as well as in reconstitutions of purified 3A4, suggesting that the electrophilic species reacted preferentially with the P450 enzymes concerned.

D-Tyrosine as a chiral precusor to potent inhibitors of human nonpancreatic secretory phospholipase A(2) (IIa) with antiinflammatory activity

Few reported inhibitors of secretory phospholipase A(2) enzymes inhibit the IIa human isoform (hnpsPLA(2)-IIa) noncovalently at submicromolar concentrations. Herein, the simple chiral precursor D-tyrosine was derivastised to give a series of potent new inhibitors of hnpsPLA(2)-IIa. A 2.2-Angstrom crystal structure shows an inhibitor bound in the active site of the enzyme, chelated to a Ca2+ ion through carboxylate and amide oxygen atoms, H bonded through an amide NH group to His48, with multiple hydrophobic contacts and a T-shaped aromatic-group-His6 interaction. Antiinflammatory activity is also demonstrated for two compounds administered orally to rats.