Pathways of ligand clearance in acetylcholinesterase by multiple copy sampling

The clearance of seven different ligands from the deeply buried active-site of Torpedo californica acetylcholinesterase is investigated by combining multiple copy sampling molecular dynamics simulations, with the analysis of protein-ligand interactions, protein motion and the electrostatic potential sampled by the ligand copies along their journey outwards. The considered ligands are the cations ammonium, methylammonium, and tetramethylammonium, the hydrophobic methane and neopentane, and the anionic product acetate and its neutral form, acetic acid. We find that the pathways explored by the different ligands vary with ligand size and chemical properties. Very small ligands, such as ammonium and methane, exit through several routes. One involves the main exit through the mouth of the enzyme gorge, another is through the so-called back door near Trp84, and a third uses a side door at a direction of approximately 45 degrees to the main exit. The larger polar ligands, methylammonium and acetic acid, leave through the main exit, but the bulkiest, tetramethylammonium and neopentane, as well as the smaller acetate ion, remain trapped in the enzyme gorge during the time of the simulations. The pattern of protein-ligand contacts during the diffusion process is highly non-random and differs for different ligands. A majority is made with aromatic side-chains, but classical H-bonds are also formed. In the case of acetate, but not acetic acid, the anionic and neutral form, respectively, of one of the reaction products, specific electrostatic interactions with protein groups, seem to slow ligand motion and interfere with protein flexibility; protonation of the acetate ion is therefore suggested to facilitate clearance. The Poisson-Boltzmann formalism is used to compute the electrostatic potential of the thermally fluctuating acetylcholinesterase protein at positions actually visited by the diffusing ligand copies. Ligands of different charge and size are shown to sample somewhat different electrostatic potentials during their migration, because they explore different microscopic routes. The potential along the clearance route of a cation such as methylammonium displays two clear minima at the active and peripheral anionic site. We find moreover that the electrostatic energy barrier that the cation needs to overcome when moving between these two sites is small in both directions, being of the order of the ligand kinetic energy. The peripheral site thus appears to play a role in trapping inbound cationic ligands as well as in cation clearance, and hence in product release.     

Acetylcholinesterase: Role of the enzyme's charge distribution in steering charged ligands toward the active site

The electrostatic steering of charged ligands toward the active site of Torpedo californica acetylcholinesterase is investigated by Brownian dynamics simulations of wild type enzyme and several mutated forms, in which some normally charged residues are neutralized. The simulations reveal that the total ligand influx through a surface of 42 Å radius centered in the enzyme monomer and separated from the protein surface by I-14 Å is not significantly influenced by electrostatic interactions. Electrostatic effects are visible for encounters with a surface of 32 Å radius, which is partially hidden inside the protein, but mostly within the solvent. A clear accumulation of encounter events for that sphere is observed in the area directly above the entrance to the active site gorge. In this area, the encounter events are increased by 40% compared to the case of a neutral ligand. However, the differences among the encounter rates for the various mutants considered here are not pronounced, all rate constants being within ±10% of the average value. The enzyme charge distribution becomes more important as the charged ligand moves toward the bottom of the gorge, where the active site is located. We show that neither the enzyme's total charge, nor its dipole moment, fully account for the electrostatic steering of ligand to the active site. Higher moments of the enzyme's charge distribution are also important. However, for a series of mutations for which the direction of the enzyme dipole moment is constant within a few degrees, one observes a gradual decrease in the diffusional encounter rate constant with the number of neutralized residues. On the other hand, for other mutants that change the direction of the dipole moment from that of the wild type, the calculated encounter rate constants can be very close to that of the wild type. The present work yields two new insights to the kinetics of acetylcholinesterase. First, evolution appears to have built a redundant electrostatic steering capability into this important enzyme through the overall distribution of its thousands of partially charged atoms. And second, roughly half of the rate enhancement due to electrostatics arises from steering of the substrate outside the enzyme; the other half of the rate enhancement arises from improved trapping of the substrate after it has entered the gorge. The computational results reproduce qualitatively, and help to rationalize, many surprising experimental results obtained recently for human acetylcholinesterase. © 1996 John Wiley & Sons, Inc.     

Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant k_on for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction. © 1998 John Wiley & Sons, Inc. Biopoly 46: 465–474, 1998     

Long Route or Shortcut? A Molecular Dynamics Study of Traffic of Thiocholine within the Active-Site Gorge of Acetylcholinesterase

The principal role of acetylcholinesterase is termination of nerve impulse transmission at cholinergic synapses, by rapid hydrolysis of the neurotransmitter acetylcholine to acetate and choline. Its active site is buried at the bottom of a deep and narrow gorge, at the rim of which is found a second anionic site, the peripheral anionic site. The fact that the active site is so deeply buried has raised cogent questions as to how rapid traffic of substrate and products occurs in such a confined environment. Various theoretical and experimental approaches have been used to solve this problem. Here, multiple conventional molecular dynamics simulations have been performed to investigate the clearance of the product, thiocholine, from the active-site gorge of acetylcholinesterase. Our results indicate that thiocholine is released from the peripheral anionic site via random pathways, while three exit routes appear to be favored for its release from the active site, namely, along the axis of the active-site gorge, and through putative back- and side-doors. The back-door pathway is that via which thiocholine exits most frequently. Our results are in good agreement with kinetic and kinetic-crystallography studies. We propose the use of multiple molecular dynamics simulations as a fast yet accurate complementary tool in structural studies of enzymatic trafficking.     

Interaction kinetics of reversible inhibitors and substrates with acetylcholinesterase and its fasciculin 2 complex

Fasciculin 2 (Fas2), a three-fingered peptide of 61 amino acids, binds tightly to the peripheral site of acetylcholinesterases (AChE; EC ), occluding the entry portal into the active center gorge of the enzyme and inhibiting its catalytic activity. We investigated the mechanism of Fas2 inhibition by studying hydrolysis of cationic and neutral substrates and by determining the kinetics of interaction for fast equilibrating cationic and neutral reversible inhibitors with the AChE.Fas2 complex and free AChE. Catalytic parameters, derived by eliminating residual Fas2-resistant activity, reveal that Fas2 reduces k(cat)/K(m) up to 10(6)-fold for cationic substrates and less than 10(3)-fold for neutral substrates. Rate constants for association of reversible inhibitors with the active center of the AChE.Fas2 complex were reduced about 10(4)-fold for both cationic and neutral inhibitors, while dissociation rate constants were reduced 10(2)-to 10(3)-fold, compared with AChE alone. Rates of ligand association with both AChE and AChE.Fas2 complex were dependent on the protonation state of ionizable ligands but were also markedly reduced by protonation of enzyme residue(s) with pK(a) of 6.1-6.2. Linear free energy relationships between the equilibrium constant and the kinetic constants show that Fas2, presumably through an allosteric influence, markedly alters the position of the transition state in the reaction pathway. Since Fas2 complexation introduces an energetic barrier for hydrolysis of substrates that exceeds that found for association of reversible ligands, Fas2 influences catalytic parameters by a more complex mechanism than simple restriction of diffusional entry and exit from the active center. Conformational flexibility appears critical for facilitating ligand passage in the narrow active center gorge for both AChE and the AChE.Fas2 complex.     

Blind Prediction of Charged Ligand Binding Affinities in a Model Binding Site

Predicting absolute protein–ligand binding affinities remains a frontier challenge in ligand discovery and design. This becomes more difficult when ionic interactions are involved because of the large opposing solvation and electrostatic attraction energies. In a blind test, we examined whether alchemical free-energy calculations could predict binding affinities of 14 charged and 5 neutral compounds previously untested as ligands for a cavity binding site in cytochrome c peroxidase. In this simplified site, polar and cationic ligands compete with solvent to interact with a buried aspartate. Predictions were tested by calorimetry, spectroscopy, and crystallography. Of the 15 compounds predicted to bind, 13 were experimentally confirmed, while 4 compounds were false negative predictions. Predictions had a root-mean-square error of 1.95 kcal/mol to the experimental affinities, and predicted poses had an average RMSD of 1.7 Å to the crystallographic poses. This test serves as a benchmark for these thermodynamically rigorous calculations at predicting binding affinities for charged compounds and gives insights into the existing sources of error, which are primarily electrostatic interactions inside proteins. Our experiments also provide a useful set of ionic binding affinities in a simplified system for testing new affinity prediction methods.     

The role of electrostatic interactions in the absorption of ligands to the active sites of cholinesterases,as indicated by molecular modeling data

The effect of electrostatic interactions on the absorption of the positively charged reversible inhibitor tetram-ethylammonium, its neutral structural analogue neopentane C(CH₃)₄, and the natural substrate acethylcholine to the active sites of acetylcholinesterase and butyrylcholinesterase has been studied by molecular modeling methods. It has been shown that the dominant absorption of positively charged ligands is due to the anchoring of the cationic group of the ligand in the anionic subsite of both enzymes through the interaction of the π-cation with the benzene ring of tryptophan. The correlation between the free energy of complex formation and the catalytic activity of charged tetramethylammonium has been revealed for both enzymes. It has been shown that the effective binding of the acetylcholine molecule requires the additional activation of the enzyme.     

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.     

Structure and dynamics of the active site gorge of acetylcholinesterase: synergistic use of molecular dynamics simulation and X-ray crystallography.

The active site of acetylcholinesterase (AChE) from Torpedo californica is located 20 A from the enzyme surface at the bottom of a narrow gorge. To understand the role of this gorge in the function of AChE, we have studied simulations of its molecular dynamics. When simulations were conducted with pure water filling the gorge, residues in the vicinity of the active site deviated quickly and markedly from the crystal structure. Further study of the original crystallographic data suggests that a bis-quaternary decamethonium (DECA) ion, acquired during enzyme purification, residues in the gorge. There is additional electron density within the gorge that may represent small bound cations. When DECA and 2 cations are placed within the gorge, the simulation and the crystal structure are dramatically reconciled. The small cations, more so than DECA, appear to stabilize part of the gorge wall through electrostatic interactions. This part of the gorge wall is relatively thin and may regulate substrate, product, and water movement through the active site.     

Mechanism of reactant and product dissociation from the anthrax edema factor: a locally enhanced sampling and steered molecular dynamics study

The anthrax edema factor is a toxin overproducing damaging levels of cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi) from ATP. Here, mechanisms of dissociation of ATP and products (cAMP, PPi) from the active site are studied using locally enhanced sampling (LES) and steered molecular dynamics simulations. Various substrate conformations and ionic binding modes found in crystallographic structures are considered. LES simulations show that PPi and cAMP dissociate through different solvent accessible channels, while ATP dissociation requires significant active site exposure to solvent. The ionic content of the active site directly affects the dissociation of ATP and products. Only one ion dissociates along with ATP in the two-Mg(2+) binding site, suggesting that the other ion binds EF prior to ATP association. Dissociation of reaction products cAMP and PPi is impaired by direct electrostatic interactions between products and Mg(2+) ions. This provides an explanation for the inhibitory effect of high Mg(2+) concentrations on EF enzymatic activity. Breaking of electrostatic interactions is dependent on a competitive binding of water molecules to the ions, and thus on the solvent accessibility of the active site. Consequently, product dissociation seems to be a two-step process. First, ligands are progressively solvated while preserving the most important electrostatic interactions, in a process that is dependent on the flexibility of the active site. Second, breakage of the electrostatic bonds follows, and ligands diffuse into solvent. In agreement with this mechanism, product protonation facilitates dissociation.     

Diffusion of Particles in the Extracellular Matrix: The Effect of Repulsive Electrostatic Interactions

Diffusive transport of macromolecules and nanoparticles in charged fibrous media is of interest in many biological applications, including drug delivery and separation processes. Experimental findings have shown that diffusion can be significantly hindered by electrostatic interactions between the diffusing particle and charged components of the extracellular matrix. The implications, however, have not been analyzed rigorously. Here, we present a mathematical framework to study the effect of charge on the diffusive transport of macromolecules and nanoparticles in the extracellular matrix of biological tissues. The model takes into account steric, hydrodynamic, and electrostatic interactions. We show that when the fiber size is comparable to the Debye length, electrostatic forces between the fibers and the particles result in slowed diffusion. However, as the fiber diameter increases the repulsive forces become less important. Our results explain the experimental observations that neutral particles diffuse faster than charged particles. Taken together, we conclude that optimal particles for delivery to tumors should be initially cationic to target the tumor vessels and then change to neutral charge after exiting the blood vessels.     

The fasciculin-acetylcholinesterase interaction

Acetylcholinesterase (AChE) terminates the action of the neurotransmitter acetylcholine at cholinergic synapses in the central and peripheral nervous systems. Fasciculins, which belong to the family of "three-fingered" snake toxins, selectively inhibit mammalian AChEs with Ki values in the picomolar range. In solution, the cationic fasciculin appears to bind to the enzyme's peripheral anionic site, located near the mouth of the gorge leading to the active center, to inhibit catalysis either allosterically or by creating an electrostatic barrier at the gorge entry (or both). Yet the crystal structure of the fasciculin-mouse AChE complex, which shows that the central loop of fasciculin fits snugly at the entrance of the gorge, suggests that the mode of action of fasciculin is steric occlusion of substrate access to the active center. Mutagenesis of the fasciculin molecule, undertaken to establish a functional map of the binding surfaces, identified determinants common to those identified by the structural approach and revealed that only a few of the many fasciculin residues residing at the complex interface provide the strong contacts required for high affinity binding and enzyme inhibition. However, it did not reconcile the disparity between the kinetic and structural data. Finally, the crystal structure of mouse AChE without bound fasciculin shows a tetrameric assembly of subunits; within the tetramer, a short loop at the surface of a subunit associates with the peripheral site of a facing subunit and sterically occludes the entrance of the active center gorge. The position and complementarity of the peripheral site-occluding loop mimic the characteristics of the central loop of fasciculin bound to AChE. This suggests not only that the peripheral site of AChE is a site for association of heterologous proteins with interactive surface loops, but also that endogenous peptidic ligands of AChE sharing structural features with the fasciculin molecule might exist.     

Electrostatic attraction by surface charge does not contribute to the catalytic efficiency of acetylcholinesterase.

Acetylcholinesterases (AChEs) are characterized by a high net negative charge and by an uneven surface charge distribution, giving rise to a negative electrostatic potential extending over most of the molecular surface. To evaluate the contribution of these electrostatic properties to the catalytic efficiency, 20 single- and multiple-site mutants of human AChE were generated by replacing up to seven acidic residues, vicinal to the rim of the active-center gorge (Glu84, Glu285, Glu292, Asp349, Glu358, Glu389 and Asp390), by neutral amino acids. Progressive simulated replacement of these charged residues results in a gradual decrease of the negative electrostatic potential which is essentially eliminated by neutralizing six or seven charges. In marked contrast to the shrinking of the electrostatic potential, the corresponding mutations had no significant effect on the apparent bimolecular rate constants of hydrolysis for charged and non-charged substrates, or on the Ki value for a charged active center inhibitor. Moreover, the kcat values for all 20 mutants are essentially identical to that of the wild type enzyme, and the apparent bimolecular rate constants show a moderate dependence on the ionic strength, which is invariant for all the enzymes examined. These findings suggest that the surface electrostatic properties of AChE do not contribute to the catalytic rate, that this rate is probably not diffusion-controlled and that long-range electrostatic interactions play no role in stabilization of the transition states of the catalytic process.     

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.     

Role of charged residues in the catalytic mechanism of hepatitis C virus NS3 protease: electrostatic precollision guidance and transition-state stabilization

Maturational cleavage of the hepatitis C virus polyprotein involves the viral chymotrypsin-like serine protease NS3. The substrate binding site of this enzyme is unusually flat and featureless. We here show that NS3 has a highly asymmetric charge distribution that is characterized by strong positive potentials in the vicinity of its active site and in the S5/S6 region. Using electrostatic potential calculations, we identified determinants of this positive potential, and the role of six different residues was explored by site-directed mutagenesis. Mutation of residues in the vicinity of the active site led to changes in k(cat) values of a peptide substrate indicating that basic amino acids play a role in the stabilization of the transition state. Charge neutralization in the S5/S6 region increased the K(m) values of peptide substrates in a manner that depended on the presence of negatively charged residues in the P5 and P6 positions. K(i) values of hexapeptide acids spanning P6-P1 (product inhibitors) were affected by charge neutralization in both the active site region and the S5/S6 region. Pre-steady-state kinetic data showed that the electrostatic surface potential is used by this enzyme to enhance collision rates between peptidic ligands and the active site. Calculations of the interaction energies of protease-substrate or protease-inhibitor complexes showed that electrostatic interaction energies oppose the formation of a tightly bound complex due to an unfavorable change in the desolvation energy. We propose that desolvation costs are minimized by avoiding the formation of individual ion pair interactions through the use of clusters of positively charged residues in the generation of local electrostatic potentials.     

A neutral molecule in a cation-binding site: specific binding of a PEG-SH to acetylcholinesterase from Torpedo californica

The crystal structure of acetylcholinesterase from Torpedo californica complexed with the uncharged inhibitor, PEG-SH-350 (containing mainly heptameric polyethylene glycol with a terminal thiol group) is determined at 2.3 A resolution. This is an untypical acetylcholinesterase inhibitor, since it lacks the cationic moiety typical of the substrate (acetylcholine). In the crystal structure, the elongated ligand extends along the whole of the deep and narrow active-site gorge, with the terminal thiol group bound near the bottom, close to the catalytic site. Unexpectedly, the cation-binding site (formed by the faces of aromatic side-chains) is occupied by CH(2) groups of the inhibitor, which are engaged in C-H...pi interactions that structurally mimic the cation-pi interactions made by the choline moiety of acetylcholine. In addition, the PEG-SH molecule makes numerous other weak but specific interactions of the C-H...O and C-H...pi types.     

Homology Built Model of Acetylcholinesterase from Drosophila Melanogaster

Abstract

Acetylcholinesterases from Drosophila melanogaster and Torpedo marmorata possess 35% identical residues. We built a homology model of the Drosophila enzyme on the basis of the known three-dimensional structure of Torpedo acetylcholinesterase, which revealed an oval rim of the active site gorge with an additional hollow which could accept small charged ligands more firmly than the corresponding surface in the Torpedo enzyme. This difference at the peripheral site, together with the kinetics of W121A and W359L mutants, suggests coordinate action of important hydrophobic residues that form the active site gorge during the catalytic process. It may also account for the activation-inhibition kinetic pattern which is characteristic for the insect enzyme.     

Molecular docking study on the “back door” hypothesis for product clearance in acetylcholinesterase

Acetylcholinesterase (AChE) is one of the fastest enzymes known, even though the active site is buried inside the protein at the end of a 20-A deep narrow gorge. Among the great variety of crystal structures of this enzyme, both in the absence and presence of various ligands and proteins, the structure of a complex of AChE with the pseudo-irreversible inhibitor Mf268 is of particular interest, as it assists in the proposal of a back door for product clearance from the active site. Binding of Mf268 to AChE results in the carbamoylation of Ser200 and liberation of an eseroline-fragment as the leaving group. The crystal structure of the AChE-Mf268 complex, however, proves that eseroline has escaped from the enzyme, despite the fact that the Ser-bound inhibitor fragment blocks the gorge entrance. The existence of alternative routes other than through the gorge for product clearance has been postulated but is still controversially discussed in the literature, as an experimental proof for such a back door is still missing. We have used Monte Carlo-based molecular docking methods in order to examine possible alternative pathways that could allow eseroline to be released from the protein after being cleaved from the substrate by Ser200. Based on our results, a short channel at the bottom of the gorge seems to be the most probable back-door site, which begins at amino acid Trp84 and ends at the enzyme surface in a cavity close to amino acid Glu445. [Figure: see text].     

Inhibitors tethered near the acetylcholinesterase active site serve as molecular rulers of the peripheral and acylation sites

The acetylcholinesterase (AChE) active site consists of a narrow gorge with two separate ligand binding sites: an acylation site (or A-site) at the bottom of the gorge where substrate hydrolysis occurs and a peripheral site (or P-site) at the gorge mouth. AChE is inactivated by organophosphates as they pass through the P-site and phosphorylate the catalytic serine in the A-site. One strategy to protect against organophosphate inactivation is to design cyclic ligands that will bind specifically to the P-site and block the passage of organophosphates but not acetylcholine. To accelerate the process of identifying cyclic compounds with high affinity for the AChE P-site, we introduced a cysteine residue near the rim of the P-site by site-specific mutagenesis to generate recombinant human H287C AChE. Compounds were synthesized with a highly reactive methanethiosulfonyl substituent and linked to this cysteine through a disulfide bond. The advantages of this tethering were demonstrated with H287C AChE modified with six compounds, consisting of cationic trialkylammonium, acridinium, and tacrine ligands with tethers of varying length. Modification by ligands with short tethers had little effect on catalytic properties, but longer tethering resulted in shifts in substrate hydrolysis profiles and reduced affinity for acridinium affinity resin. Molecular modeling calculations indicated that cationic ligands with tethers of intermediate length bound to the P-site, whereas those with long tethers reached the A-site. These binding locations were confirmed experimentally by measuring competitive inhibition constants KI2 for propidium and tacrine, inhibitors specific for the P- and A-sites, respectively. Values of KI2 for propidium increased 30- to 100-fold when ligands had either intermediate or long tethers. In contrast, the value of KI2 for tacrine increased substantially only when ligands had long tethers. These relative changes in propidium and tacrine affinities thus provided a sensitive molecular ruler for assigning the binding locations of the tethered cations.     

Determinants of cytochrome P450 2C8 substrate binding: structures of complexes with montelukast, troglitazone, felodipine, and 9-cis-retinoic acid

Although a crystal structure and a pharmacophore model are available for cytochrome P450 2C8, the role of protein flexibility and specific ligand-protein interactions that govern substrate binding are poorly understood. X-ray crystal structures of P450 2C8 complexed with montelukast (2.8 A), troglitazone (2.7 A), felodipine (2.3 A), and 9-cis-retinoic acid (2.6 A) were determined to examine ligand-protein interactions for these chemically diverse compounds. Montelukast is a relatively large anionic inhibitor that exhibits a tripartite structure and complements the size and shape of the active-site cavity. The inhibitor troglitazone occupies the upper portion of the active-site cavity, leaving a substantial part of the cavity unoccupied. The smaller neutral felodipine molecule is sequestered with its dichlorophenyl group positioned close to the heme iron, and water molecules fill the distal portion of the cavity. The structure of the 9-cis-retinoic acid complex reveals that two substrate molecules bind simultaneously in the active site of P450 2C8. A second molecule of 9-cis-retinoic acid is located above the proximal molecule and can restrain the position of the latter for more efficient oxygenation. Solution binding studies do not discriminate between cooperative and noncooperative models for multiple substrate binding. The complexes with structurally distinct ligands further demonstrate the conformational adaptability of active site-constituting residues, especially Arg-241, that can reorient in the active-site cavity to stabilize a negatively charged functional group and define two spatially distinct binding sites for anionic moieties of substrates.