Poly(triazolyl methacrylate) glycopolymers as potential targeted unimolecular nanocarriers

Synthetic glycopolymers are increasingly investigated as multivalent ligands for a range of biological and biomedical applications. This study indicates that glycopolymers with a fine-tuned balance between hydrophilic sugar pendant units and relatively hydrophobic polymer backbones can act as single-chain targeted nanocarriers for low molecular weight hydrophobic molecules. Non-covalent complexes formed from poly(triazolyl methacrylate) glycopolymers and low molecular weight hydrophobic guest molecules were characterised through a range of analytical techniques - DLS, SLS, TDA, fluorescence spectroscopy, surface tension analysis - and molecular dynamics (MD) modelling simulations provided further information on the macromolecular characteristics of these single chain complexes. Finally, we show that these nanocarriers can be utilised to deliver a hydrophobic guest molecule, Nile red, to both soluble and surface-immobilised concanavalin A (Con A) and peanut agglutinin (PNA) model lectins with high specificity, showing the potential of non-covalent complexation with specific glycopolymers in targeted guest-molecule delivery.

Development and evaluation of human AP endonuclease inhibitors in melanoma and glioma cell lines

Aims:

Modulation of DNA base excision repair (BER) has the potential to enhance response to chemotherapy and improve outcomes in tumours such as melanoma and glioma. APE1, a critical protein in BER that processes potentially cytotoxic abasic sites (AP sites), is a promising new target in cancer. In the current study, we aimed to develop small molecule inhibitors of APE1 for cancer therapy.

Methods:

An industry-standard high throughput virtual screening strategy was adopted. The Sybyl8.0 (Tripos, St Louis, MO, USA) molecular modelling software suite was used to build inhibitor templates. Similarity searching strategies were then applied using ROCS 2.3 (Open Eye Scientific, Santa Fe, NM, USA) to extract pharmacophorically related subsets of compounds from a chemically diverse database of 2.6 million compounds. The compounds in these subsets were subjected to docking against the active site of the APE1 model, using the genetic algorithm-based programme GOLD2.7 (CCDC, Cambridge, UK). Predicted ligand poses were ranked on the basis of several scoring functions. The top virtual hits with promising pharmaceutical properties underwent detailed in vitro analyses using fluorescence-based APE1 cleavage assays and counter screened using endonuclease IV cleavage assays, fluorescence quenching assays and radiolabelled oligonucleotide assays. Biochemical APE1 inhibitors were then subjected to detailed cytotoxicity analyses.

Results:

Several specific APE1 inhibitors were isolated by this approach. The IC₅₀ for APE1 inhibition ranged between 30 nM and 50 μM. We demonstrated that APE1 inhibitors lead to accumulation of AP sites in genomic DNA and potentiated the cytotoxicity of alkylating agents in melanoma and glioma cell lines.

Conclusions:

Our study provides evidence that APE1 is an emerging drug target and could have therapeutic application in patients with melanoma and glioma.

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

The mysteries of telomere structure and recognition: could radioprobing help?

PURPOSE

Telomeres are specialized DNA-protein complexes found at the ends of eukaryotic chromosomes. In normal somatic cells these become shorter with each cell division and appear to control their replicative lifespan. However almost all tumours show activation of the enzyme telomerase, a specialised reverse transcriptase/DNA polymerase, that can add new telomeric repeats to the ends of chromosomes and this appears to be a key factor in the cell immortalization process. Consequently there is much current interest in the potential for inhibitors of telomere extension in the treatment of cancer. Several groups have found that it is possible to produce inhibitory molecules that target the telomeric repeat (substrate) DNA rather than the telomerase enzyme itself. This is thought to work because it has been found that in vitro, these DNA sequences can fold up into a four-stranded (quadruplex) structure that the drugs recognise and stabilize, but which is not recognised by the enzyme. However, while medicinal chemists continue to base rational design programs on this hypothesis, there is currently very little evidence that these structures form in vivo, and that in vivo the drugs work by binding to them. To have incontrovertible evidence of where and how these telomerase inhibitors and DNA interact is therefore a pressing concern for a basic understanding of their mechanism of action and effective drug development.

MATERIALS AND METHODS

Radioprobing represents a valuable new approach to the study of DNA structures. Recently we have shown through computer simulations of radioprobing that the technique is a remarkably sensitive probe of quite fine details of DNA conformation. Here we report on our simulations of the binding of a radiolabelled telomerase inhibitor, related to a class of novel inhibitors under development at Nottingham, to a variety of possible structures for telomeric DNA.

RESULTS AND CONCLUSIONS

The predicted cleavage patterns prove to be very sensitive to the DNA structure, and the mode of binding of the drug. These results suggest that radioprobing experiments should be able to provide unambiguous evidence as to the 'true' nature of the telomere-drug complexes, and so aid the rational design programme.

Supercoiling and denaturation of DNA loops

We study the statistical mechanics of small DNA loops emphasizing the competition between elasticity, supercoiling, and denaturation. Motivated by recent experiments and atomistic molecular dynamics simulation, we propose a new coarse-grained phenomenological model of DNA. We extend the classical elastic rod models to include the possibility of denaturation and nonlinear twist elasticity. Using this coarse-grained model, we obtain a phase diagram in terms of fractional overtwist and loop size that can be used to rationalize a number of experimental results which have also been confirmed by atomistic simulations.

A simple physical description of DNA dynamics: quasi-harmonic analysis as a route to the configurational entropy

It has become increasingly apparent that the dynamic as well as the structural properties of biological macromolecules are important to their function. However, information concerning molecular flexibility can be difficult to obtain experimentally at the atomic level. Computer modelling techniques such as molecular dynamics (MD) have therefore proved invaluable in advancing our understanding of biomolecular flexibility. This paper describes how a combination of atomistic MD simulations and quasi-harmonic analysis can be used to describe the dynamics of duplex DNA, with a particular emphasis on methods for calculating differences in configurational entropies. We demonstrate that DNA possesses remarkably simple mechanical properties relative to globular proteins, making it an ideal system for exploring biomolecular flexibility in general. Our results also highlight the importance of solvent viscosity in determining the dynamic behaviour of DNA in aqueous solution.

Auger electrons--a nanoprobe for structural, molecular and cellular processes

This paper provides a brief review of recently published work on biophysical and biological aspects of Auger processes. Three specific questions have been considered. (1) Does charge neutralisation contribute to molecular damage such as DNA strand breaks? (2) How many DNA double strand breaks are produced by a single decay of DNA bound (125)I? (3) What is the correlation between number of gammaH2AX foci and number of double strand breaks (DSB)? The paper also gives preliminary reports on two new calculations: (a) calculation of the spectrum of Auger electrons released during decay of (124)I and (b) the use of Auger electrons in the decay of (125)I as a probing agent of novel DNA structures.

Cooperativity in drug-DNA recognition: a molecular dynamics study

NMR studies have shown that the minor groove-binding ligand Hoechst 33258 binds to the two T4/A4 tracts within the duplex d(CTTTTCGAAAAG)2 in a highly cooperative manner, such that in titration experiments no intermediate 1:1 complex can be detected. The NMR-derived structures of the free DNA and the 2:1 complex have been obtained, but can shed little light on what the origins of this cooperativity may be. Here we present the results of a series of molecular dynamics simulations on the free DNA, the 1:1 complex, and the 2:1 complex, which have been designed to enable us to calculate thermodynamic parameters associated with the molecular recognition events. The results of the molecular dynamics studies confirm that structural factors alone cannot explain the cooperativity observed, indeed when enthalpic and hydration factors are looked at in isolation, the recognition process is predicted to be slightly anticooperative. However, when changes in configurational entropy are taken into account as well, the overall free energy differences are such that the calculated cooperativity is in good agreement with that observed experimentally. The results indicate the power of molecular dynamics methods to provide reasonable explanations for phenomena that are difficult to explain on the basis of static models alone, and provide a nice example of the concept of "allostery without conformational change".

Enantioselectivity of some 1-[(benzofuran-2-yl) phenylmethyl] imidazoles as aromatase (P450AROM) inhibitors

The enantioselectivity ratio ((+)-:(-)-forms) of three substituted 1-[(benzofuran-2-yl) phenylmethyl] imidazoles as inhibitors of aromatase (P450AROM) was 2.16, 12.3 and 1.0 for the 4-methyl-, 4-fluoro- and 4-chloro-substituted compounds, respectively. The (+/-)-compounds were all >1000 times more potent than (+/-)-aminoglutethimide (IC50 = 12 x 10(3) nM). High potency (5.3-65.0 nM) for all the enantiomers studied is unusual since activity usually resides in one form for chiral inhibitors of P450AROM. The 4-methyl derivative was fitted into the model [Furet, P., Batzl, C., Bhatnager, A.S., Francotte, E., Rihs, G. and Lang, M. (1993) J. Med. Chem. 36, pp. 1393-1400] for binding of S-(-)-fadrazole to the active site and the (R)- and (S)- forms both gave a good fitting pattern with (S)-(-)-fadrazole so accounting for their close activity. Docking of both forms into the active site model for P450AROM [Laughton, C.A., Zvelebil, M.J.J.M. and Neidel, S. (1993) J. Steroid Biochem. Mol. Biol. 44, pp. 399-407], using the orientation of (S)-(-)-fadrazole, gave similar strong binding along the position of the C and D rings of the steroid substrate and in the hydrophobic cavity below the A/B rings. The site was probed for group size accommodation using the less potent 4-phenyl analogue (IC50(+/-) = 242 nM): the (S)-form showed restricted access to the region under the A ring due to the extended bulk of the biphenyl group.

Force-induced melting of a short DNA double helix

The dynamic behaviour of DNA is of fundamental importance to many cellular processes. One principal characteristic, central to transcription and replication, is the ability of the duplex to "melt". It has recently been shown that dynamic force spectroscopy provides information about the energetics of biomolecular dissociation. We have employed this technique to investigate the unbinding of single dodecanucleotide molecules. To separate the duplex to single-stranded DNA, forces ranging from 17 to 40 pN were required over a range of loading rates. Interpretation of the dependence of melting force on loading rate revealed that the energy barrier to rupture is between 9 and 13 kcal mol-1 in height and situated 0.58 nm from an intermediate structural state. Thermal melting studies show that, prior to dissociation, the oligonucleotide underwent a transition which required between 7 and 11 kcal mol-1 in energy. Through combined dynamic force spectroscopy and thermal melting studies we show the derivation of an energy landscape to dissociate a 12-mer duplex. Until very recently, this type of information was only accessible by computational analysis. Additionally, the force spectroscopy data allow an estimation of the kinetics of duplex formation and melting.

Structural studies on bioactive compounds. 34. Design, synthesis, and biological evaluation of triazenyl-substituted pyrimethamine inhibitors of Pneumocystis carinii dihydrofolate reductase

The triazenyl-pyrimethamine derivative 3a (TAB), a potent and selective inhibitor of Pneumocystis carinii DHFR, was selected as the starting point for a lead optimization study. Molecular modeling studies, corroborated by a recent crystal structure determination of the ternary complex of P. carinii DHFR--NADPH bound to TAB, predicted that modifications to the acetoxy residue of the lead inhibitor could exploit binding opportunities in the vicinity of an active site pocket bounded by residues Ile33, Lys37, and Leu72. Substitutions in the benzyl moiety with electron-donating and electron-withdrawing groups were predicted to probe face-edge interactions with amino acid Phe69 unique to the P. carinii enzyme. New triazenes 10a--v and 12a--f were prepared by coupling the diazonium tetrafluoroborate salt 6b of aminopyrimethamine with substituted benzylamines or phenethylamines. The most potent of the new inhibitors against P. carinii DHFR was the naphthylmethyl-substituted triazene 10t (IC(50): 0.053 microM), but a more substantial increase in potency against the rat liver DHFR led to a reduction in selectivity (ratio rat liver DHFR IC(50)/P. carinii DHFR IC(50): 5.36) compared to the original lead structure 3a (ratio rat liver DHFR IC(50)/P. carinii DHFR IC(50): 114).

Dynamic model of base pair breathing in a DNA chain with a defect

We model and analyze a short section of a DNA chain with a defect, with the aim of understanding how the frequency, amplitude, and localization of breathing events depend on the strength of the bonds between base pairs, both along the chain and between the chains. Our results show that the presence of a defect in the chain permits the existence of a localized breather mode. The models we analyze are linear and hence solvable, with solvability extending to the statistical mechanics formulation of the problem. Parameter values for the interaction energy of a base with its nearest neighbors are obtained from AMBER. The results indicate good agreement with both the amplitude and the number of base pairs affected by defect-induced breathing motion.

Distribution of strand breaks produced by Auger electrons in decay of 125I in triplex DNA

In this study we investigate the possibility of using Auger electrons as a probing agent for the study of structures of nucleic acids. To this end, we present the distribution of breaks produced in strands of a DNA duplex and a triplex-forming oligonucleotide (TFO) carrying Auger emitting radionuclide 125I. The method of calculation includes use of a molecular model of plasmid DNA duplex with bound TFO carrying a labelled 125I at position C5 of a single deoxycytosine residue, a source of Auger spectra, Monte Carlo electron track structure and the ensuing chemistry codes, to simulate the distribution of breaks produced in both strands of a plasmid DNA. Frequencies of fragment length distributions were obtained for the TFO, the purine and the pyrimidine strands. The frequency of breaks in the purine strand showed good correlation with the published experimental results, while that for the pyrimidine strand is lower by a factor of 3. It is concluded that the true structure of triplex DNA may not be purely of B-form.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)(2) and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar-phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving 'lock and key' recognition.

Evaluation of the mechanism of aromatase cytochrome P450. A site-directed mutagenesis study

Aromatase (CYP19) catalyzes three consecutive hydroxylation reactions converting C19 androgens to aromatic C18 estrogenic steroids. In this study, five human aromatase mutants (E302D, S478A, S478T, H480K, and H480Q) were prepared using a mammalian cell expression system. These mutants were evaluated by enzyme kinetic analysis, inhibitory profile studies, and reaction intermediate measurements. Three steroidal inhibitors [4-hydroxyandrostenedione (4-OHA), 7alpha-(4'-amino)phenylthio-1,4-androstandiene-3,17-dione (7alpha-APTADD), and bridge (2,19-methyleneoxy) androstene-3,17-dione (MDL 101003)], and four nonsteroidal inhibitors [aminoglutethimide (AG), CGS 20267, ICI D1033, and vorozole (R83842)] were used in the inhibitory profile studies. Our computer model of aromatase suggests that Glu302 is situated in the conserved I-helix region and located near the C-19 position of the steroid substrate. The model was supported by significant changes in kinetic parameters and a sevenfold increase in the Ki value of MDL 101,003 for the mutant E302D. As S478A was found to have kinetic properties similar to the wild-type enzyme and a much higher activity than S478T, Ser478 is thought to be situated in a rather restricted environment. There was a 10-fold increase in the Ki value of 7alpha-APTADD for S478T over that for the wild-type enzyme, suggesting that Ser478 might be near the C-7 position of the substrate. The reaction intermediate analysis revealed that significantly more 19-ol intermediate was generated by both S478A and S478T than the wild-type enzyme. These results would support a hypothesis that Ser478 plays a role in the first and second hydroxylation reactions. A positive charged amino acid is preferred at position 480 as shown by the fact that H480K has a significantly higher activity than H480Q. The Ki value of 4-OHA for H480Q was found to be three times that of the wild-type enzyme. In addition, significantly more 19-ol and 19-al intermediates were detected for both mutants H480K and H480Q than for the wild-type enzyme. Evaluation of the two mutations at His480 allows us to propose that this residue may participate in the aromatization reaction (the third step) by acting as a hydrogen bond donor for the C-3 keto group of the substrate. Furthermore, new products were generated when the enzyme was mutated at Ser478 and His480. Thus, these two residues must play an important role in the catalysis and are likely closer to the substrate binding site than previously predicted.

A method for radioprobing DNA structures using Auger electrons

PURPOSE

To present a new method for radioprobing a DNA triple helix structure by Auger electrons emitted in the decay of 125I using theoretical/computational approaches.

MATERIALS AND METHODS

A Monte Carlo track structure method was used to simulate the damage to a triplex resulting from Auger electrons emitted in the decay of an incorporated 125I atom in plasmid DNA. Comparison of the theoretical frequency distributions of single-strand breaks induced on the Pu and Py strands with the experimental data and a knowledge of the distances from the strand breaks to the iodine provide information on the structures otherwise difficult to obtain with X-ray crystallography.

RESULTS

In comparing theoretical frequency distributions of single-strand breaks with the experimental data it is found that the results are very sensitive to the conformation of the triplex model used. It is found that the best fit to the experimental data results from using a hybrid triplex model, in which the base-step geometry is A-like, while the sugar puckers adopt the B-like C2'-endo conformation.

CONCLUSIONS

The approach and technique presented here represent a valuable new addition to the methods available for DNA structure determination since they provide information on medium-range structure otherwize difficult to obtain in the absence of X-ray crystallography. It is concluded that currently accepted models for triplex structure are not optimal, and a modified structure is proposed that fits the radioprobing results better, while maintaining agreement with the fibre diffraction and NMR data. Although the method has proved to be very useful for scoring alternative trial solutions, further studies combining experimental data from multiple iodine positions with track structure modelling are required for directing structural optimization.

Atomic force microscopy studies of intercalation-induced changes in plasmid DNA tertiary structure

Structural transitions in the tertiary structure of plasmid DNA have been investigated using atomic force microscopy. Changes in superhelical stress were induced by ethidium bromide intercalation, and conformational effects monitored by recording topographic images from DNA complexes of various ethidium bromide:base pair stoichiometry. Significant changes in the tertiary structure of individual DNA molecules were observed with increasing ethidium bromide concentration. The first distinct conformational transition was from a predominantly relaxed structure to one consisting solely of toroidal supercoils. A further increase in ethidium bromide concentration resulted in the formation of regions of plectonemic supercoiling. The ratio of plectonemic:toroidal supercoiling gradually increased until an extremely tightly interwound structure of solely plectonemic supercoiling was finally adopted. The toroidal form of supercoiling observed in this study is unusual as both atomic force microscopy and electron microscopy techniques have previously shown that plectonemic supercoiling is the predominant form adopted by plasmid DNA.

Structural studies on bioactive compounds. 30. Crystal structure and molecular modeling studies on the Pneumocystis carinii dihydrofolate reductase cofactor complex with TAB, a highly selective antifolate

The crystal structure of the ternary complex of NADPH, the potent antifolate [2, 4-diamino-5-¿3-[3-(2-acetyloxyethyl)-3-benzyltriazen-1-yl]-4 -chloroph enyl¿-6-ethylpyrimidine] (TAB, 1) and Pneumocystis carinii dihydrofolate reductase (pcDHFR), refined to 2.1 A resolution, reveals that TAB binds similar to the antifolates trimethoprim and methotrexate. These data also reveal multiple conformations for the binding geometry of TAB with two preferred orientations of the acetyloxy and benzyl groups that results from a 180 degrees rotation about the N2-N3 triazenyl bond. The methyl of the acetyloxy and benzyl ring of TAB probes large hydrophobic regions of the p-aminobenzoyl folate binding pocket of the active site, in particular the region near Phe69, which is unique to the pcDHFR sequence. These results confirm prior molecular modeling investigations of the binding of TAB to pcDHFR that identified four low-energy binding geometries, two involving rotations about the terminal N(2)-N(3) triazenyl linkage and two involving atropisomerism about the pivotal pyrimethamine-phenyl bond. The primary differences in the molecular dynamics (MD) models and those observed in this crystal complex result from small conformational changes in active-site residues on energy minimization. However, two MD models place the acetyloxy and benzyl ring groups in a region of the active site between the cofactor-binding region and the p-aminobenzoyl folate pocket; an orientation never observed in any DHFR crystal structure to date. These conformers interact with solvent near the enzyme surface and are probably not observed due to the loss of specific hydrogen bonds with the enzyme. The high species pcDHFR selectivity of TAB could be the result of ligand flexibility that enables multiple binding orientations at the enzyme active site. Further modification of the acetyloxy region of TAB could increase its potency and selectivity for pcDHFR.

Structural studies on bioactive compounds. Part 29: palladium catalysed arylations and alkynylations of sterically hindered immunomodulatory 2-amino-5-halo-4,6-(disubstituted)pyrimidines

The immunological agent bropirimine 5 is a tetra-substituted pyrimidine with anticancer and interferon-inducing properties. Synthetic routes to novel 5-aryl analogues of bropirimine have been developed and their potential molecular recognition properties analysed by molecular modelling methods. Sterically challenged 2-amino-5-halo-6-phenylpyrimidin-4-ones (halo = Br or I) are poor substrates for palladium catalysed Suzuki cross-coupling reactions with benzeneboronic acid because the basic conditions of the reaction converts the amphoteric pyrimidinones to their unreactive enolic forms. Palladium-mediated reductive dehalogenation of the pyrimidinone substrates effectively competes with cross-coupling. 2-Amino-5-halo-4-methoxy-6-phenylpyrimidines can be converted to a range of 5-aryl derivatives with the 5-iodopyrimidines being the most efficient substrates. Hydrolysis of the 2-amino-5-aryl-4-methoxy-6-phenylpyrimidines affords the required pyrimidin-4-ones in high yields. Semi-empirical quantum mechanical calculations show how the nature of the 5-substituent influences the equilibrium between the 1H- and 3H-tautomeric forms, and the rotational freedom about the bond connecting the 6-phenyl group and the pyrimidine ring. Both of these factors may influence the biological properties of these compounds.

Solution structure and dynamics of the A-T tract DNA decamer duplex d(GGTAATTACC)2: implications for recognition by minor groove binding drugs

The structure of the DNA decamer duplex d(GGTAATTACC)(2) has been determined using NMR distance restraints and molecular dynamics simulations of 500 ps to 1 ns in aqueous solution at 300 K. Using both canonical A and canonical B starting structures [root-mean-square deviation (RMSD) 4.6 A; 1 A=10(-10) m], with and without experimental restraints, we show that all four simulations converge to a similar envelope of final conformations with B-like helical parameters (pairwise RMSD 1.27-2.03 A between time-averaged structures). While the two restrained simulations reach a stable trajectory after 300-400 ps, the unrestrained trajectories take longer to equilibrate. We have analysed the dynamic aspects of these structures (sugar pucker, helical twist, roll, propeller twist and groove width) and show that the minor groove width in the AATT core of the duplex fluctuates significantly, sampling both wide and narrow conformations. The structure does not have the highly pre-organized narrow minor groove generally regarded as essential for recognition and binding by small molecules, suggesting that ligand binding carries with it a significant component of 'induced-fit'. Our simulations show that there are significant differences in structure between the TpA step (where p=phosphate) and the ApA and ApT steps, where a large roll into the major groove at the TpA step appears to be an important factor in widening the minor groove at this position.

Molecular recognition between a new pentacyclic acridinium salt and DNA sequences investigated by optical spectroscopic techniques, proton nuclear magnetic resonance spectroscopy, and molecular modeling

A pentacyclic acridine, 1H-2,3-dihydroindolizino[7,6,5-kl]acridinium chloride (1), related in structure to tetra- and pentacyclic marine natural products, has previously been shown to induce apoptosis in breast and non-small-cell lung tumor cell lines and shows significant differences in biological potency and antitumor profile from other intercalating agents based on the acridine framework. We report on the molecular recognition of the acridinium salt with DNA, quantified by optical spectroscopic methods, and have compared these results with the clinical agent amsacrine (m-AMSA). The results point to an intercalative association between 1 and G-C-rich sequences of DNA. We have synthesized a hexamer duplex d(ACGCGT)2, presenting two potential 5'-CpG recipient sites, and have investigated in detail by NMR and molecular modeling methods the orientational preferences of 1, particularly with regard to the pyrrolidine ring system. On the basis of the intermolecular nuclear Overhauser effect (NOE) data, four possible intercalation models were considered; no single model produced a significantly better fit than any of the others. The best fit to the experimental data was obtained by considering a dynamic equilibrium between the different intercalated orientations with the drug maximizing pi-overlap with the G-C base pairs at the intercalation site. We found little evidence for any degree of groove specificity imparted by the pyrrolidine ring. If these simulations have biological relevance they suggest that, at most, the agent induces only a transitory hot spot in the DNA which, evidently, is sufficient to be sensed by damage-recognition mechanisms of the cell.

Production and characterization of an anti-(MUC1 mucin) recombinant diabody

A recombinant diabody fragment based on the anti-MUC1 monoclonal antibody, C595 has been produced in a bacterial expression system. Substitution of a 7-amino-acid linker sequence (Gly6Ser) for the original single-chain (sc)Fv 15-amino-acid linker (Gly4-Ser)3, using polymerase-chain-reaction-based strategies, forces variable heavy (V(H)) and light (V(L)) domains to pair with complementary domains on neighbouring scFv molecules, forming a scFv dimer (diabody). This recombinant protein shows similar binding characteristics to the parental C595 monoclonal antibody. The ability to bind to MUC1 mucin on carcinoma cell surfaces will allow its potential as a diagnostic and therapeutic reagent of clinical utility to be investigated.

DNA minor groove recognition by a tetrahydropyrimidinium analogue of hoechst 33258: NMR and molecular dynamics studies of the complex with d(GGTAATTACC)2

Hoechst 43254 (H43254), a 2,3,4,5-tetrahydropyrimidin-1-ium analogue of the bis-benzimidazole minor groove binding agent Hoechst 33258 (H33258), has been studied by NMR and restrained molecular dynamics in its complex with d(GGTAATTACC)2. We investigate the origin of the enhanced complex stability afforded by the replacement of the N-methylpiperazine ring of H33258 with the tetrahydropyrimidinium ring of H43254, the latter presenting the opportunity for specific minor groove-directed recognition through a pyrimidinium NH. A set of 25 drug-DNA NOEs define the binding site with some precision and are used as part of the structural analysis using restrained molecular dynamics simulations considering explicit solvation and the treatment of electrostatic interactions using the particle mesh Ewald method within AMBER 4.1. Starting with three different initial structures with the drug located at different sites in the groove (pairwise RMSD 4.3-12.6 A) we arrive at three very similar structures (pairwise RMSD 0.80-1.34 A) representing one converged binding site at the centre of the AATT tract. Two of the three structures show the tetrahydropyrimidinium ring to be suitably positioned for an -NH to adenine N3 hydrogen bond suggesting that electrostatic interactions may play an important role in the enhanced affinity as well as imparting additional A-T specificity. The NMR data show that the pyrimidinium NH interaction is dynamic since signal averaging from the two sides of the ring indicate rapid rotations in the bound form.

Molecular basis of the inhibition of human aromatase (estrogen synthetase) by flavone and isoflavone phytoestrogens: A site-directed mutagenesis study

Flavone and isoflavone phytoestrogens are plant chemicals and are known to be competitive inhibitors of cytochrome P450 aromatase with respect to the androgen substrate. Aromatase is the enzyme that converts androgen to estrogen; therefore, these plant chemicals are thought to be capable of modifying the estrogen level in women. In this study, the inhibition profiles of four flavones [chrysin (5, 7-dihydroxyflavone), 7,8-dihydroxyflavone, baicalein (5,6,7-trihydroxyflavone), and galangin (3,5,7-trihydroxyflavone)], two isoflavones [genistein (4,5,7-trihydroxyisoflavone) and biochanin A (5,7-dihydroxy-4-methoxyisoflavone)], one flavanone [naringenin (4, 5,7-trihydroxyflavanone)], and one naphthoflavone (alpha-naphthoflavone) on the wild-type and six human aromatase mutants (I133Y, P308F, D309A, T310S, I395F, and I474Y) were determined. In combination with computer modeling, the binding characteristics and the structure requirement for flavone and isoflavone phytoestrogens to inhibit human aromatase were obtained. These compounds were found to bind to the active site of aromatase in an orientation in which rings A and C mimic rings D and C of the androgen substrate, respectively. This study also provides a molecular basis as to why isoflavones are significantly poorer inhibitors of aromatase than flavones.

Binding characteristics of aromatase inhibitors and phytoestrogens to human aromatase

We have evaluated the binding characteristics of three steroidal inhibitors [4-hydroxyandrostenedione (4-OHA), 7alpha-(4'-amino)phenylthio-1,4-androstadiene-3,17-dione (7alpha-APTADD), and bridge (2,19-methyleneoxy) androstene-3,17-dione (MDL 101,003)], four nonsteroidal inhibitors [aminoglutethimide (AG), CGS 20267, ICI D1033, and vorozole (R83842)], and two flavone phytoestrogens (chrysin, and 7,8-dihydroxyflavone) to aromatase through a combination of computer modeling and inhibitory profile studies on the wild-type and six aromatase mutants (I133Y, P308F, D309A, T310S, I395F, and I474Y). We have generated two aromatase models based on the x-ray structures of cytochrome P450-cam and cytochrome P450bm3, respectively. A major difference between the cytochrome P450cam-based and cytochrome P450bm3-based models is in the predicted lengths of helices F and G. In the cytochrome P450cam-based model, helices F and G lie antiparallel and extend across the active-site face of the molecule from one edge to the center, so that the carboxyl-terminal residues of helix F and the N-terminal residues of helix G make a major contribution to the structure of the active site. In the cytochrome P450bm3-based model, both helices are longer and so extend almost all the way across the active-site face of the molecule. Considering the size of the androgen substrate, we evaluated our results mainly based on the cytochrome P450cam model. The mutations involved in this study are thought to be at or near the proposed active site pocket. The inhibitory profile analysis has produced very interesting results and provided a molecular basis as to how seven aromatase inhibitors with different structures bind to the active site of aromatase. Furthermore, the investigation reveals that phytoestrogens bind to the active site of aromatase in a different orientation from that in the estrogen receptor.

Homology modelling of the enzyme P450 17 alpha-hydroxylase/17,20-lyase--a target for prostate cancer chemotherapy--from the crystal structure of P450BM-3

Using homology modelling techniques the structure of cytochrome P450 17 alpha-hydroxylase/17-20 lyase (P450c17) has been predicted from the crystal structure of P450BM-3. The resulting structure has been compared to a previous model, built on the basis of homology to P450cam. Despite considerable structural differences between the two template structures, the two derived models for P450c17 show a high degree of similarity in the active-site region. As before, a bilobal active-site cavity is predicted, and we hypothesize that binding of the steroid in one lobe of the active site is associated with the hydroxylase activity of the enzyme, whilst binding in the other is associated with the lyase reaction.

Efficient triple helix formation by oligodeoxyribonucleotides containing alpha- or beta-2-amino-5-(2-deoxy-D-ribofuranosyl) pyridine residues

Triple helices containing C+xGxC triplets are destabilised at physiological pH due to the requirement for base protonation of 2'-deoxycytidine (dC), which has a pKa of 4.3. The C nucleoside 2-amino-5-(2'-deoxy-beta-D-ribofuranosyl)pyridine (beta-AP) is structurally analogous to dC but is considerably more basic, with a pKa of 5.93. We have synthesised 5'-psoralen linked oligodeoxyribonucleotides (ODNs) containing thymidine (dT) and either beta-AP or its alpha-anomer (alpha-AP) and have assessed their ability to form triplexes with a double-stranded target derived from standard deoxynucleotides (i.e. beta-anomers). Third strand ODNs derived from dT and beta-AP were found to have considerably higher binding affinities for the target than the corresponding ODNs derived from dT and either dC or 5-methyl-2'-deoxycytidine (5-Me-dC). ODNs containing dT and alpha-AP also showed enhanced triplex formation with the duplex target and, in addition are more stable in serum-containing medium than standard oligopyrimidine-derived ODNs or ODNs derived from dT and beta-AP. Molecular modelling studies showed that an alpha-anomeric AP nucleotide can be accommodated within an otherwise beta-anomeric triplex with only minor perturbation of the triplex structure. Molecular dynamics (MD) simulations on triplexes containing either the alpha- or beta-anomer of (N1-protonated) AP showed that in both cases the base retained two standard hydrogen bonds to its associated guanine when the 'A-type' model of the triplex was used as the start-point for the simulation, but that bifurcated hydrogen bonds resulted when the alternative 'B-type' triplex model was used. The lack of a differential stability between alpha-AP- and beta-AP-containing triplexes at pH >7, predicted from the behaviour of the B-type models, suggests that the A-type models are more appropriate.

Binding characteristics of seven inhibitors of human aromatase: a site-directed mutagenesis study

Aromatase, a cytochrome P450, catalyzes three consecutive hydroxylation reactions converting C19 androgens to aromatic C18 estrogenic steroids. In this study, eight human aromatase mutants (I133Y, I133W, F235L, I395F, I474Y, I474W, I474M, and I474N) were prepared to evaluate the active site and a proposed hydrophobic pocket of the enzyme that exists in an aromatase model based on the X-ray structure of cytochrome P450cam. In addition, the binding characteristics of three steroidal inhibitors [4-hydroxyandrostenedione, 7alpha-(4'-amino)phenylthio-1,4-androstandiene-3,17-dione, and bridge (2,19-methyleneoxy)androstene-3,17-dione (MDL 101,003)] and four nonsteroidal inhibitors [aminoglutethimide, CGS 20267, ICI D1033, and vorozole (R83842)] were investigated through inhibitory profile studies on the eight new and three previously generated mutants (P308F, D309A, and T310S). The latter analyses have provided a molecular basis regarding how seven aromatase inhibitors with different structures bind to the active site of aromatase.

A crystallographic and spectroscopic study of the complex between d(CGCGAATTCGCG)2 and 2,5-bis(4-guanylphenyl)furan, an analogue of berenil. Structural origins of enhanced DNA-binding affinity

2,5-Bis(4-guanylphenyl)furan ("furamidine") is a dicationic minor groove binding drug that has been shown to be more effective than pentamidine against the Pneumocystis carinii pathogen in an immunosuppressed rat model. It has a close structural similarity to the antitrypanosomal drug berenil, differing only on the replacement of the central triazene unit with a furan moiety. we have determined the crystal structure of the complex between furamidine and the DNA dodecamer d(CGCGAATTCGCG)2 and compared it with the same DNA sequence by UV-visible, fluorescence, and CD spectroscopy. Furamidine shows tighter binding to this sequence (Keq = 6.7 x 10(6)) than berenil (Keq = 6.6 x 10(5)). The crystal structure reveals that, unlike berenil, furamidine makes direct hydrogen bond interactions with this DNA sequence through both amidinium groups to O2 atoms of thymine bases and is more isohelical with the minor groove. Molecular mechanics calculations support the hypothesis that these differences result in the improved interaction energy between the ligand and the DNA.

Characteristics of triplex-directed photoadduct formation by psoralen-linked oligodeoxynucleotides

A triplex-forming oligopyrimidine has been attached at its 5'-end to a photoreactive psoralen derivative and used to target a sequence which forms part of the coding region of the human aromatase gene. The 20 base pair sequence is not a perfect triplex target since it contains three pyrimidine interruptions within the purine-rich strand. Despite this, we have detected triplex-directed photoadduct formation at pH 7.0 between the psoralen-linked oligonucleotide and a 30mer duplex representing the aromatase target. Photoadduct formation was found to be sensitive to pH, temperature, cation concentration and the base composition of the third strand. By varying the base sequence of the target duplex around the psoralen intercalation site, we have characterised the site and mode of psoralen intercalation. The attached psoralen has been found to intercalate at the triplex-duplex junction with a strong preference for one orientation. We have shown that the psoralen will bind at the junction even when there is a preferred TpA step at an adjacent site. We have also compared the binding affinity and photoreactivity of oligodeoxyribonucleotides linked to two different psoralen derivatives and found differences in the rate of crosslinking and the extent of crosslink formation. Finally, we have examined oligodeoxyribonucleotides which are attached to psoralen by polymethylene linkers of different lengths.

Detection and kinetic studies of triplex formation by oligodeoxynucleotides using real-time biomolecular interaction analysis (BIA)

Real-time biomolecular interaction analysis (BIA) has been applied to triplex formation between oligodeoxynucleotides. 5'-Biotinylated oligonucleotides were immobilised on the streptavidin-coated surface of a biosensor chip and subsequently hybridised to their complementary strand. Sequence-specific triplex formation was observed when a suitable third-strand oligopyrimidine was injected over the surface-bound duplex. In addition, a single-stranded oligonucleotide immobilised on the chip surface was able to capture a DNA duplex by triplex recognition. The presence of spermine increases the rate of association between the third strand and immobilised duplex, but at elevated spermine concentrations non-specific association is observed. A preliminary kinetic analysis of triplex formation at pH 5.2 by an 11mer third strand containing thymine, cytosine and uracil is reported. Values for the association and dissociation rate constants were determined to be (1.9 +/- 0.2) x 10(3) M-1 s-1 and (8.1 +/- 1.9) x 10(-5) s-1, respectively.

An approach to protein homology modelling based on an ensemble of NMR structures: application to the Sox-5 HMG-box protein

A new approach has been developed to reduce multiple protein structures obtained from NMR structure analysis to a smaller number of representative structures which still reflect the structural diversity of the data sets. The method, based on the clustering of similar structures, has been tested in the homology model building of the structure of Sox-5, a sequence-specific DNA-binding protein belonging to the high mobility group (HMG) nuclear proteins family. Sox (SRY box) genes are the autosomal genes related to the sex-determining SRY, Y chromosomal gene. The Sox-5 protein, encoded by one of the SRY-related genes, displays a 29% sequence identity with the HMG1 B-box domain whose structure, determined previously by NMR, has been used in our study to predict the structure of Sox-5. Two independent ensembles of HMG1 structures, each represented by closely related coordinate sets, were used. Nine representative structures for HMG1 were subsequently selected as starting points for the modelling of Sox-5. The model of the protein shows close similarity to the HMG1 fold, with differences at the secondary structure level located mainly in alpha-helices 1 and 3. A left-handed, three residue per turn polyproline II helix, forming a conserved polyproline II/alpha-helix supersecondary motif, was identified in the N-terminal region of Sox-5 and other HMG boxes.

Rationally designed analogues of tamoxifen with improved calmodulin antagonism

Computerized molecular modeling studies on the interactions of the antiestrogen tamoxifen (1) and its analogues bound to the calcium-binding protein calmodulin have guided the rational design of more potent antagonists. Compounds with either three or four methylene units in the basic side chain or slim lipophilic 4-substituents were expected to be more potent. All compounds were tested for antagonism of the calmodulin-dependent activity of cAMP phosphodiesterase and for binding affinity to the estrogen receptor from rat uteri. Some compounds were assayed for cytotoxicity against MCF-7 breast tumor cells in vitro. Introduction of lipophilic 4-substituents was accomplished by using palladium(0)-catalyzed coupling reactions with a 4-iodinated precursor. Both the 4-ethynyl (16 and 17) and 4-butyl (18 and 19) compounds were more potent calmodulin antagonists than tamoxifen. Extension of the basic aminoethoxy side chain of 4-iodotamoxifen (3) and idoxifene (2) ((E)-1-[4-[2-(N-pyrrolidino)ethoxy]phenyl]-1-(4-iodophenyl)-2-phen yl-1- butene) by one or two methylene units resulted in modest gains in calmodulin antagonism (10-13). All the compounds assayed retained estrogen receptor binding characteristics. The compound possessing the optimal combination of calmodulin antagonism and estrogen receptor binding was 12 ((E)-1-[4-[3-(N-pyrrolidino)propoxy]phenyl]-1-(4-iodophenyl)-2-phe nyl-1 - butene) (IC50 = 1.1 microM, RBA = 23). Correlation between calmodulin antagonism and cytotoxicity was demonstrated for selected compounds.

Mutagenesis study at a postulated hydrophobic region near the active site of aromatase cytochrome P450

Aromatase, a cytochrome P450, catalyzes three consecutive hydroxylation reactions converting C19 androgens to aromatic C18 estrogenic steroids. On the basis of a recent computer modeling of the active site of aromatase, a hydrophobic pocket, thought to be important for the binding of some aromatase inhibitors, was predicted to extend roughly in the plane of the steroid substrate, from the position that would be occupied by its C4 and C7 atoms. Four mutants, G121A, I125N, F235N, and I474F, were generated to test this model. Although the mutagenesis results have shown that the current model for the active site of aromatase almost certainly contains a number of errors, the results are in general very satisfactory in that they suggest how the model should be altered by local realignments of the aromatase sequence with that of cytochrome P450cam. Among the mutants, I474F is the most interesting one. Its Km value for androstenedione was found to be lower than the wild type enzyme, and the kinetic analysis exhibited a substrate inhibition-like kinetic profile through an "in-cell" assay. In addition, this mutation reduces the binding affinity of an aromatase inhibitor, 4-hydroxyandrostenedione, and increases the binding affinity of two aromatase inhibitors, aminoglutethimide and CGS 16949. This study demonstrates a useful approach, by a combination of computer modeling, site-directed mutagenesis, and inhibitor binding studies, to examine the structure of the active site of aromatase and the binding nature of various aromatase inhibitors.

A study of simulated annealing protocols for use with molecular dynamics in protein structure prediction

The method of simulated annealing can be of use in protein structure prediction by homology modelling where side chain conformations must be predicted. In this study an attempt has been made to optimize a molecular dynamics method for this purpose. Heating and cooling protocols to maximize the accuracy of the predictions have been developed. The optimized protocol involves cooling from 3000 to 0 K over 20 ps while simultaneously introducing the non-bonded energy term. The use of a 'soft' non-bonded interaction energy term in place of a standard 6-12 potential is found to be important. The reliability of the predictions has been analysed in terms of the environment of the residues (solvent accessibility) and the degree of uncertainty in the structure (number of unknown torsion angles). Depending on these factors the percentage of unknown side chain torsion angles that are correctly predicted within 30 degrees ranges from approximately 50 to 75%. Potential problems and limitations of the method are discussed.

Prediction of protein side-chain conformations from local three-dimensional homology relationships

A method for predicting the conformations of protein side-chains, starting from main-chain co-ordinates alone, is described. The method involves the comparison of the local environment of each residue whose side-chain conformation is to be predicted with a database of local environments for the same residue type constructed from an analysis of high-resolution protein structures. Local environments are described in terms of the residue type and location in space of residues that interact with the side-chain of interest. The best (most three-dimensionally homologous) few matches to each residue are then input to a Monte-Carlo procedure to give a final predicted structure. The method has been tested on a selection of eight proteins, ranging in size from 46 to 323 amino acid residues. The average side-chain atom root-mean-square deviation between the actual and predicted structures is 1.71 A taken over all residues, and 1.00 A if restricted to buried residues. Over all residues, an average of 59.8% of all side-chain dihedral angles are predicted within +/- 30 degrees of the crystal structure values. Considering buried residues only, this rises to 79.6%.

A detailed molecular model for human aromatase

Using a variety of techniques, including sequence alignment, secondary structure prediction, molecular mechanics and molecular dynamics, we have constructed a model for the three-dimensional structure of P-450arom (human aromatase) based on that of P-450cam, the only cytochrome P-450 enzyme for which the crystal structure is known. The predicted structure is found to be in good agreement with current experimental data; both direct, from site-directed mutagenesis studies, and indirect, from the consideration of the structures and activities of known substrates and inhibitors.

Prediction of the structure of the Y+.R-.R(+)-type DNA triple helix by molecular modelling

Molecular mechanics has been used to predict the structure of the Y+.R-.R(+)-type DNA triple helix, in which a second polypurine strand binds antiparallel to the homopurine strand of a homopurine/homopyrimidine stretch of duplex DNA. From calculations on the sequence d(C)10.d(G)10.d(G)10, two likely structures emerge. One has the glycosidic torsions of the third strand bases in the anti-conformation and Hoogsteen hydrogen-bonds to the purine strand of the duplex, the other has the third strand purines in the syn orientation and uses a reverse-Hoogsteen hydrogen-bonding pattern. Despite the large structural differences between these two types of triplex, calculations performed in vacuo with a distance-dependent dielectric constant to mimic the shielding effect of solvent show them to be energetically very similar, with the latter (syn) slightly preferred. However, if explicit solvent molecules are included in the calculation, the anti conformation is found to be much preferred. This difference in the results seems to stem from an underestimation of short-range electrostatic interactions in the in vacuo simulations. When TAA or TAT base triples are substituted for the sixth CGG triple in the sequence, it is found that, for the solvated model, the third strand base of the TAA triple prefers the syn orientation while that in the TAT triple retains a preference, though reduced, for the anti conformation.

A note on the conformational flexibility of the antiestrogenic drug tamoxifen: preferred conformations in the free state and bound to the protein calmodulin

The conformational properties of the antiestrogenic drug tamoxifen, a triphenylbut-1-ene derivative, have been studied using molecular mechanics. Four distinct conformers have been identified, and the energy barriers between them have been established. The orientation of the ethyl group substituent has been examined in particular, since the lowest-energy conformers have this group orientated 180 degrees away from its position in the crystal structures of tamoxifen and its derivatives. These differences have implications for the interactions of tamoxifen with the calcium-binding protein calmodulin; relevant results from a molecular-modelling study of this protein-drug complex are presented.

A molecular modeling study of the interactions between the antiestrogen drug tamoxifen and several derivatives, and the calcium-binding protein calmodulin

The interactions of the antiestrogenic drug tamoxifen with the calcium-binding protein calmodulin have been studied by computerized molecular modeling methods. Sites in both the N and C domains of the protein have been established, with one in the C domain having the highest calculated enthalpy of binding. The residues involved in the sites have been detailed. Modeling studies are reported for six tamoxifen derivatives, and their calculated enthalpies of binding are compared with the ability of the analogues to inhibit calmodulin-dependent cyclic AMP phosphodiesterase (PDE) (Rowlands et al. Biochem, Pharmacol. 1990, 40, 283-289). The poor binding properties of the piperazino and C-methyl derivatives are correctly predicted, whereas the superior affinity of 4-iodotamoxifen is not fully explained by the model.

Molecular dynamics simulation of the DNA triplex d(TC)5.d(GA)5.d(C+T)5

A molecular dynamics simulation of the DNA triple helix d(TC)5.d(GA)5.d(C+T)5 is described (C+ represents a protonated cytosine residue). The simulation has been performed using the program AMBER 3.1 and includes counterions and explicit solvent under periodic boundary conditions. Both the dynamic and time-averaged behaviour of the system has been analysed. Considerable deviations from the fibre-diffraction model for DNA triple helix structure are observed, including the repuckering of the purine strand sugars that has been identified in some nuclear magnetic resonance (n.m.r.) studies. The simulation suggests that this conformational change may be driven by the possibility of improved interactions between the phosphate groups of this strand and both the solvent and counterions. Several examples of a particular conformational transition are observed, involving correlated changes in the backbone angles alpha and gamma. These transitions provide a possible explanation for some unusual n.m.r. data that have been reported. The structure of the triple helix major groove also suggests an explanation for the observed stabilization of DNA triplexes by polyvalent cations, and their ability to interact with drugs that bind in the minor groove of DNA duplexes.

Inhibitors of the P450 enzymes aromatase and lyase. Crystallographic and molecular modeling studies suggest structural features of pyridylacetic acid derivatives responsible for differences in enzyme inhibitory activity

Derivatives of 4-pyridylacetic acid are known to be inhibitors of the cytochrome P450 enzymes aromatase and lyase (17 alpha-hydroxylase/C17-20lyase), and are therefore of interest in the treatment of hormone dependent breast and prostate cancers. We report the determination of the crystal structure of one such derivative, the 4-tert-butyl cyclohexyl ester, and molecular modeling studies on two related inhibitors, the cyclohexyl ester and its alpha-methyl derivative. These latter two compounds show a marked difference in their relative activities against aromatase and lyase. Two models are proposed for the interaction of these molecules with the target enzymes on the basis of their ability to adopt conformations that partially mimic steroid substrates. From these models an explanation can be advanced for the fact that, compared with the unmethylated analogue, the (racemic) alpha-methylated compound is seven times poorer as an inhibitor of aromatase but seven times better as an inhibitor of lyase. The model proposed for binding to aromatase places the alpha-carbon of the ester group in the position occupied by C(2) of steroid substrates. In contrast, that proposed for binding to lyase places this atom in the position occupied by C(17) of steroid substrates. The introduction of steric bulk at C(2) is known to be unfavorable for aromatase inhibition, while its introduction at C(17) may lead to a better mimicry of the steroid D-ring and so improve lyase inhibition.

Molecular-mechanics modelling of drug-DNA structures; the effects of differing dielectric treatment on helix parameters and comparison with a fully solvated structural model

This study analyses the influence that the nature of the dielectric constant has on the final structures obtained from in vacuo molecular mechanics calculations on a drug-DNA complex and compares these structures with the energy minimised complex including explicit solvent molecules. Minimisations have been performed on a proflavine-decanucleotide structure, where the drug was intercalated at the d(CpG) site of the d(GpApTpApCpGpApTpApC) decamer duplex, using two expressions for the dielectric constant: a distance-independent, epsilon ij = EPS, and a distance-dependent, epsilon ij = EPS*Rij, form and for values of EPS from 1 to 8. Significantly different structures are obtained for the distance-independent and the distance-dependent expressions of the dielectric constant. The use of a distance-independent dielectric constant leads to distorted structures, which are very sensitive to slight changes in the value of EPS. The use of a distance-dependent dielectric constant leads to less distorted and more stable structures. The effects on helical parameters are analysed in detail. The structures obtained for different values of EPS (within the distance-dependent formalism) seem to converge for values of EPS equal to 4 or greater. Based on these results a distance-dependent form of the dielectric with an EPS value of 4 is recommended in order to produce reliable refined nucleic acid structures by molecular mechanics. These conclusions have been supported by molecular-mechanics minimisation of the same structure with the inclusion of explicit water molecules and counter-ions.

A molecular model for the enzyme cytochrome P450(17 alpha), a major target for the chemotherapy of prostatic cancer

The enzyme cytochrome P450(17 alpha) catalyses two key steps in the biosynthesis of the androgens from pregnanes: the 17 alpha hydroxylation step and the subsequent 17-20 lyase reaction. Using a variety of techniques, including sequence alignment, secondary structure prediction, molecular mechanics and molecular dynamics, we have constructed a model for the three-dimensional structure of P450(17 alpha) based on that of P450cam, the only cytochrome P450 enzyme for which the crystal structure is known. The model suggests the possibility of two modes of binding of steroid substrates at the active site, perhaps reflecting the dual functionality of the enzyme.

Crystallographic and molecular modeling studies on 3-ethyl-3-(4-pyridyl)piperidine-2,6-dione and its butyl analogue, inhibitors of mammalian aromatase. Comparison with natural substrates: prediction of enantioselectivity for N-alkyl derivatives

Inhibitors of the cytochrome P450 enzyme aromatase, which is involved in the biosynthesis of estrogens from androgens, are of proven utility in the treatment of hormone-dependent breast cancer. The determination of the crystal structure of one such inhibitor, 3-ethyl-3-(4-pyridyl)piperidine-2,6-dione (2) and its 3-butyl analogue (3) is described. In the absence of three-dimensional structural information for the enzyme, conformational analysis and comparison with natural substrates has been performed in order to define possible "active" conformations. The enhanced inhibitory activity of 3 may be linked to hydrophobic interactions between the side chain and that portion of the enzyme that normally interacts with the B and C rings of a steroid substrate. Information gained from this study and previous studies by other workers has been combined in order to produce a hypothesis to explain the pattern of activity of N(1)-alkyl derivatives of 2. The successful application of this hypothesis to the prediction of the relative aromatase inhibitory activities of the two enantiomers of the N-octyl derivative (4) is described.