Free energies at QM accuracy from force fields via multimap targeted estimation

Accurate predictions of ligand binding affinities would greatly accelerate the first stages of drug discovery campaigns. However, using highly accurate interatomic potentials based on quantum mechanics (QM) in free energy methods has been so far largely unfeasible due to their prohibitive computational cost. Here, we present an efficient method to compute QM free energies from simulations using cheap reference potentials, such as force fields (FFs). This task has traditionally been out of reach due to the slow convergence of computing the correction from the FF to the QM potential. To overcome this bottleneck, we generalize targeted free energy methods to employ multiple maps-implemented with normalizing flow neural networks (NNs)-that maximize the overlap between the distributions. Critically, the method requires neither a separate expensive training phase for the NNs nor samples from the QM potential. We further propose a one-epoch learning policy to efficiently avoid overfitting, and we combine our approach with enhanced sampling strategies to overcome the pervasive problem of poor convergence due to slow degrees of freedom. On the drug-like molecules in the HiPen dataset, the method accelerates the calculation of the free energy difference of switching from an FF to a DFTB3 potential by three orders of magnitude compared to standard free energy perturbation and by a factor of eight compared to previously published nonequilibrium calculations. Our results suggest that our method, in combination with efficient QM/MM calculations, may be used in lead optimization campaigns in drug discovery and to study protein-ligand molecular recognition processes.

Hexamethylene amiloride binds the SARS-CoV-2 envelope protein at the protein-lipid interface

The SARS-CoV-2 envelope (E) protein forms a five-helix bundle in lipid bilayers whose cation-conducting activity is associated with the inflammatory response and respiratory distress symptoms of COVID-19. E channel activity is inhibited by the drug 5-(N,N-hexamethylene) amiloride (HMA). However, the binding site of HMA in E has not been determined. Here we use solid-state NMR to measure distances between HMA and the E transmembrane domain (ETM) in lipid bilayers. 13 C, 15 N-labeled HMA is combined with fluorinated or 13 C-labeled ETM. Conversely, fluorinated HMA is combined with 13 C, 15 N-labeled ETM. These orthogonal isotopic labeling patterns allow us to conduct dipolar recoupling NMR experiments to determine the HMA binding stoichiometry to ETM as well as HMA-protein distances. We find that HMA binds ETM with a stoichiometry of one drug per pentamer. Unexpectedly, the bound HMA is not centrally located within the channel pore, but lies on the lipid-facing surface in the middle of the TM domain. This result suggests that HMA may inhibit the E channel activity by interfering with the gating function of an aromatic network. These distance data are obtained under much lower drug concentrations than in previous chemical shift perturbation data, which showed the largest perturbation for N-terminal residues. This difference suggests that HMA has higher affinity for the protein-lipid interface than the channel pore. These results give insight into the inhibition mechanism of HMA for SARS-CoV-2 E.

Three Decades of REDOR in Protein Science: A Solid-State NMR Technique for Distance Measurement and Spectral Editing

Solid-state NMR (ss-NMR) is a powerful tool to investigate noncrystallizable, poorly soluble molecular systems, such as membrane proteins, amyloids, and cell walls, in environments that closely resemble their physical sites of action. Rotational-echo double resonance (REDOR) is an ss-NMR methodology, which by reintroducing heteronuclear dipolar coupling under magic angle spinning conditions provides intramolecular and intermolecular distance restraints at the atomic level. In addition, REDOR can be exploited as a selection tool to filter spectra based on dipolar couplings. Used extensively as a spectroscopic ruler between isolated spins in site-specifically labeled systems and more recently as a building block in multidimensional ss-NMR pulse sequences allowing the simultaneous measurement of multiple distances, REDOR yields atomic-scale information on the structure and interaction of proteins. By extending REDOR to the determination of 1H-X dipolar couplings in recent years, the limit of measurable distances has reached ~15-20 Å, making it an attractive method of choice for the study of complex biomolecular assemblies. Following a methodological introduction including the most recent implementations, examples are discussed to illustrate the versatility of REDOR in the study of biological systems.

Live Cells versus Fixated Cells: Kinetic Measurements of Biomolecular Interactions with the LigandTracer Method and Surface Plasmon Resonance Microscopy

Cell-based kinetic studies of ligand or candidate drug binding to membrane proteins have produced affinity and kinetic values that are different from measurements using purified proteins. However, ligand binding to fixated cells whose membrane constituents (e.g., proteins and their glycosylated forms) are partially connected by a cross-linking reagent has not been compared to that to live cells. Under the same experimental conditions for the LigandTracer method, we measured the interactions of fluorophore-labeled lectins and antibody molecules with glycans at HFF cells and the human epithelial growth receptor 2 at SKBR3 cells, respectively. In conjunction with surface plasmon resonance microscopy, the effects of labels and cell/sub-cell heterogeneity on binding kinetics were investigated. Our results revealed that, for cell constituents whose structures and functions are not closely dependent on cell viability, the ligand binding kinetics at fixated cells is only slightly different from that at live cells. The altered kinetics is explained on the basis of a less mobile receptor confined in a local environment created by partially interconnected protein molecules. We show that cell/sub-cell heterogeneity and labels on the ligands can alter the binding reaction more significantly. Thus, fixating cells not only simplifies experimental procedures for drug screening and renders assays more robust but also provides reliable kinetic information about drug binding to cell constituents whose structures are not changed by chemical fixation.

Long-range communication between transmembrane- and nucleotide-binding domains does not depend on drug binding to mutant P-glycoprotein

In this study, the impact of four P-gp mutations (G185V, G830V, F978A and ΔF335) on drug-binding and efflux-related signal-transmission mechanism was comprehensively evaluated in the presence of ligands within the drug-binding pocket (DBP), experimentally related with changes in their drug efflux profiles. The severe repacking of the transmembrane helices (TMH), induced by mutations and exacerbated by the presence of ligands, indicates that P-gp is sensitive to perturbations in the transmembrane region. Alterations on drug-binding were also observed as a consequence of the TMH repacking, but were not always correlated with alterations on ligands binding mode and/or binding affinity. Finally, and although all P-gp variants holo systems showed considerable changes in the intracellular coupling helices/nucleotide-binding domain (ICH-NBD) interactions, they seem to be primarily induced by the mutation itself rather than by the presence of ligands within the DBP. The data further suggest that the changes in drug efflux experimentally reported are mostly related with changes on drug specificity rather than effects on signal-transmission mechanism. We also hypothesize that an increase in the drug-binding affinity may also be related with the decreased drug efflux, while minor changes in binding affinities are possibly related with the increased drug efflux observed in transfected cells.Communicated by Ramaswamy H. Sarma.

Pre-Training of Equivariant Graph Matching Networks with Conformation Flexibility for Drug Binding

The latest biological findings observe that the motionless "lock-and-key" theory is not generally applicable and that changes in atomic sites and binding pose can provide important information for understanding drug binding. However, the computational expenditure limits the growth of protein trajectory-related studies, thus hindering the possibility of supervised learning. A spatial-temporal pre-training method based on the modified equivariant graph matching networks, dubbed ProtMD which has two specially designed self-supervised learning tasks: atom-level prompt-based denoising generative task and conformation-level snapshot ordering task to seize the flexibility information inside molecular dynamics (MD) trajectories with very fine temporal resolutions is presented. The ProtMD can grant the encoder network the capacity to capture the time-dependent geometric mobility of conformations along MD trajectories. Two downstream tasks are chosen to verify the effectiveness of ProtMD through linear detection and task-specific fine-tuning. A huge improvement from current state-of-the-art methods, with a decrease of 4.3% in root mean square error for the binding affinity problem and an average increase of 13.8% in the area under receiver operating characteristic curve and the area under the precision-recall curve for the ligand efficacy problem is observed. The results demonstrate a strong correlation between the magnitude of conformation's motion in the 3D space and the strength with which the ligand binds with its receptor.

Structure of ABCB1/P-Glycoprotein in the Presence of the CFTR Potentiator Ivacaftor

ABCB1/P-glycoprotein is an ATP binding cassette transporter that is involved in the clearance of xenobiotics, and it affects the disposition of many drugs in the body. Conformational flexibility of the protein within the membrane is an intrinsic part of its mechanism of action, but this has made structural studies challenging. Here, we have studied different conformations of P-glycoprotein simultaneously in the presence of ivacaftor, a known competitive inhibitor. In order to conduct this, we used high contrast cryo-electron microscopy imaging with a Volta phase plate. We associate the presence of ivacaftor with the appearance of an additional density in one of the conformational states detected. The additional density is in the central aqueous cavity and is associated with a wider separation of the two halves of the transporter in the inward-facing state. Conformational changes to the nucleotide-binding domains are also observed and may help to explain the stimulation of ATPase activity that occurs when transported substrate is bound in many ATP binding cassette transporters.

Mutation in Abl kinase with altered drug-binding kinetics indicates a novel mechanism of imatinib resistance

Protein kinase inhibitors are potent anticancer therapeutics. For example, the Bcr-Abl kinase inhibitor imatinib decreases mortality for chronic myeloid leukemia by 80%, but 22 to 41% of patients acquire resistance to imatinib. About 70% of relapsed patients harbor mutations in the Bcr-Abl kinase domain, where more than a hundred different mutations have been identified. Some mutations are located near the imatinib-binding site and cause resistance through altered interactions with the drug. However, many resistance mutations are located far from the drug-binding site, and it remains unclear how these mutations confer resistance. Additionally, earlier studies on small sets of patient-derived imatinib resistance mutations indicated that some of these mutant proteins were in fact sensitive to imatinib in cellular and biochemical studies. Here, we surveyed the resistance of 94 patient-derived Abl kinase domain mutations annotated as disease relevant or resistance causing using an engagement assay in live cells. We found that only two-thirds of mutations weaken imatinib affinity by more than twofold compared to Abl wild type. Surprisingly, one-third of mutations in the Abl kinase domain still remain sensitive to imatinib and bind with similar or higher affinity than wild type. Intriguingly, we identified three clinical Abl mutations that bind imatinib with wild type-like affinity but dissociate from imatinib considerably faster. Given the relevance of residence time for drug efficacy, mutations that alter binding kinetics could cause resistance in the nonequilibrium environment of the body where drug export and clearance play critical roles.

Functionalization of Human Serum Albumin by Tyrosine Click

Human serum albumin (HSA) is a promising drug delivery carrier. Although covalent modification of Cys34 is a well-established method, it is desirable to develop a novel covalent modification method that targets residues other than cysteine to introduce multiple functions into a single HSA molecule. We developed a tyrosine-selective modification of HSA. Three tyrosine selective modification methods, hemin-catalyzed, horseradish peroxidase (HRP)-catalyzed, and laccase-catalyzed reactions were performed, and the modification efficiencies and modification sites of the modified HSAs obtained by these methods were evaluated and compared. We found that the laccase-catalyzed method could efficiently modify the tyrosine residue of HSA under mild reaction conditions without inducing oxidative side reactions. An average of 2.2 molecules of functional groups could be introduced to a single molecule of HSA by the laccase method. Binding site analysis using mass spectrometry suggested Y84, Y138, and Y401 as the main modification sites. Furthermore, we evaluated binding to ibuprofen and found that, unlike the conventional lysine residue modification, the inhibition of drug binding was minimal. These results suggest that tyrosine-residue selective chemical modification is a promising method for covalent drug attachment to HSA.

Production of a Human Histamine Receptor for NMR Spectroscopy in Aqueous Solutions

G protein-coupled receptors (GPCRs) bind a broad array of extracellular molecules and transmit intracellular signals that initiate physiological responses. The signal transduction functions of GPCRs are inherently related to their structural plasticity, which can be experimentally observed by spectroscopic techniques. Nuclear magnetic resonance (NMR) spectroscopy in particular is an especially advantageous method to study the dynamic behavior of GPCRs. The success of NMR studies critically relies on the production of functional GPCRs containing stable-isotope labeled probes, which remains a challenging endeavor for most human GPCRs. We report a protocol for the production of the human histamine H₁ receptor (H₁R) in the methylotrophic yeast Pichia pastoris for NMR experiments. Systematic evaluation of multiple expression parameters resulted in a ten-fold increase in the yield of expressed H₁R over initial efforts in defined media. The expressed receptor could be purified to homogeneity and was found to respond to the addition of known H₁R ligands. Two-dimensional transverse relaxation-optimized spectroscopy (TROSY) NMR spectra of stable-isotope labeled H₁R show well-dispersed and resolved signals consistent with a properly folded protein, and ¹⁹F-NMR data register a response of the protein to differences in efficacies of bound ligands.

Using cryo-EM to understand antimycobacterial resistance in the catalase-peroxidase (KatG) from Mycobacterium tuberculosis

Resolution advances in cryoelectron microscopy (cryo-EM) now offer the possibility to visualize structural effects of naturally occurring resistance mutations in proteins and also of understanding the binding mechanisms of small drug molecules. In Mycobacterium tuberculosis the multifunctional heme enzyme KatG is indispensable for activation of isoniazid (INH), a first-line pro-drug for treatment of tuberculosis. We present a cryo-EM methodology for structural and functional characterization of KatG and INH resistance variants. The cryo-EM structure of the 161 kDa KatG dimer in the presence of INH is reported to 2.7 Å resolution allowing the observation of potential INH binding sites. In addition, cryo-EM structures of two INH resistance variants, identified from clinical isolates, W107R and T275P, are reported. In combination with electronic absorbance spectroscopy our cryo-EM approach reveals how these resistance variants cause disorder in the heme environment preventing heme uptake and retention, providing insight into INH resistance.

•

A cryo-EM structure to 2.7 Å resolution of M. tuberculosis KatG with isoniazid

•

Cryo-EM is able to visualize multiple dynamic binding modes of isoniazid to KatG

•

Structural disorder in isoniazid resistance mutations is observed

•

Structural disorder of the resistance mutations results in the lack of heme retention

Solution structures of the Shewanella woodyi H-NOX protein in the presence and absence of soluble guanylyl cyclase stimulator IWP-051

Heme-nitric oxide/oxygen binding (H-NOX) domains bind gaseous ligands for signal transduction in organisms spanning prokaryotic and eukaryotic kingdoms. In the bioluminescent marine bacterium Shewanella woodyi (Sw), H-NOX proteins regulate quorum sensing and biofilm formation. In higher animals, soluble guanylyl cyclase (sGC) binds nitric oxide with an H-NOX domain to induce cyclase activity and regulate vascular tone, wound healing and memory formation. sGC also binds stimulator compounds targeting cardiovascular disease. The molecular details of stimulator binding to sGC remain obscure but involve a binding pocket near an interface between H-NOX and coiled-coil domains. Here, we report the full NMR structure for CO-ligated Sw H-NOX in the presence and absence of stimulator compound IWP-051, and its backbone dynamics. Nonplanar heme geometry was retained using a semi-empirical quantum potential energy approach. Although IWP-051 binding is weak, a single binding conformation was found at the interface of the two H-NOX subdomains, near but not overlapping with sites identified in sGC. Binding leads to rotation of the subdomains and closure of the binding pocket. Backbone dynamics are similar across both domains except for two helix-connecting loops, which display increased dynamics that are further enhanced by compound binding. Structure-based sequence analyses indicate high sequence diversity in the binding pocket, but the pocket itself appears conserved among H-NOX proteins. The largest dynamical loop lies at the interface between Sw H-NOX and its binding partner as well as in the interface with the coiled coil in sGC, suggesting a critical role for the loop in signal transduction.

Large-scale binding affinity calculations on commodity compute clouds

In recent years, it has become possible to calculate binding affinities of compounds bound to proteins via rapid, accurate, precise and reproducible free energy calculations. This is imperative in drug discovery as well as personalized medicine. This approach is based on molecular dynamics (MD) simulations and draws on sequence and structural information of the protein and compound concerned. Free energies are determined by ensemble averages of many MD replicas, each of which requires hundreds of cores and/or GPU accelerators, which are now available on commodity cloud computing platforms; there are also requirements for initial model building and subsequent data analysis stages. To automate the process, we have developed a workflow known as the binding affinity calculator. In this paper, we focus on the software infrastructure and interfaces that we have developed to automate the overall workflow and execute it on commodity cloud platforms, in order to reliably predict their binding affinities on time scales relevant to the domains of application, and illustrate its application to two free energy methods.

The 3D Brain Unit Network Model to Study Spatial Brain Drug Exposure under Healthy and Pathological Conditions

PURPOSE

We have developed a 3D brain unit network model to understand the spatial-temporal distribution of a drug within the brain under different (normal and disease) conditions. Our main aim is to study the impact of disease-induced changes in drug transport processes on spatial drug distribution within the brain extracellular fluid (ECF).

METHODS

The 3D brain unit network consists of multiple connected single 3D brain units in which the brain capillaries surround the brain ECF. The model includes the distribution of unbound drug within blood plasma, coupled with the distribution of drug within brain ECF and incorporates brain capillaryblood flow, passive paracellular and transcellular BBB transport, active BBB transport, brain ECF diffusion, brain ECF bulk flow, and specific and nonspecific brain tissue binding. All of these processes may change under disease conditions.

RESULTS

We show that the simulated disease-induced changes in brain tissue characteristics significantly affect drug concentrations within the brain ECF.

CONCLUSIONS

We demonstrate that the 3D brain unit network model is an excellent tool to gain understanding in the interdependencies of the factors governing spatial-temporal drug concentrations within the brain ECF. Additionally, the model helps in predicting the spatial-temporal brain ECF concentrations of existing drugs, under both normal and disease conditions.

Pathway and mechanism of drug binding to chemokine receptors revealed by accelerated molecular simulations

Background: Chemokine GPCRs play key roles in biology and medicine. Particularly, CXCR4 promotes cancer metastasis and facilitate HIV entry into host cells. Plerixafor (PLX) is a CXCR4 drug, but the pathway and binding site of PLX in CXCR4 remain unknown. Results & methodology: We have performed molecular docking and all-atom simulations using Gaussian accelerated molecular dynamics (GaMD), which are consistent with previous mutation experiments, suggesting that PLX binds to the orthosteric site of CXCR4 as an antagonist. The GaMD simulations further revealed an intermediate allosteric binding site at the extracellular mouth of CXCR4. Conclusion: The newly identified allosteric site can be targeted for novel drug design targeting CXCR4 and other chemokine receptors.

Structural Basis for How Biologic Medicines Bind their Targets in Psoriasis Therapy

As biologic therapies become first line treatments for many inflammatory disorders, it becomes increasingly important for the practicing physician to be familiar with how these drugs function at the molecular level. This information is useful in making therapeutic decisions and helping patients understand their treatment options. It is critical to patient safety and clinical response that the molecular differences between these drugs inform prescribing practices. To this end, we present and analyze the available structural biology information about the biologics used in the treatment of psoriasis including inhibitors of tumor necrosis factor alpha (TNFα), interleukin-17 (IL-17), and interleukin-23 (IL-23). We describe and analyze the molecular surface character of known binding epitopes for medications in these classes, showing that significant differences exist in epitope location, hydrophobicity, and charge. Some of these differences can be correlated with clinical data, but our analysis ultimately points to the need for more structural information to allow for a better understanding of the structure-function relationship of biologic therapies.

Valproic acid interactions with the NavMs voltage-gated sodium channel

Valproic acid (VPA) is an anticonvulsant drug that is also used to treat migraines and bipolar disorder. Its proposed biological targets include human voltage-gated sodium channels, among other membrane proteins. We used the prokaryotic NavMs sodium channel, which has been shown to be a good exemplar for drug binding to human sodium channels, to examine the structural and functional interactions of VPA. Thermal melt synchrotron radiation circular dichroism spectroscopic binding studies of the full-length NavMs channel (which includes both pore and voltage sensor domains), and a pore-only construct, undertaken in the presence and absence of VPA, indicated that the drug binds to and destabilizes the channel, but not the pore-only construct. This is in contrast to other antiepileptic compounds that have previously been shown to bind in the central hydrophobic core of the pore region of the channel, and that tend to increase the thermal stability of both pore-only constructs and full-length channels. Molecular docking studies also indicated that the VPA binding site is associated with the voltage sensor, rather than the hydrophobic cavity of the pore domain. Electrophysiological studies show that VPA influences the block and inactivation rates of the NavMs channel, although with lower efficacy than classical channel-blocking compounds. It thus appears that, while VPA is capable of binding to these voltage-gated sodium channels, it has a very different mode and site of action than other anticonvulsant compounds.

Interaction of Aripiprazole With Human α1-Acid Glycoprotein

We recently reported that aripiprazole binds strongly to human albumin. In continuing our investigations, we investigated the mechanism responsible for the binding and the related interactions of aripiprazole with α1-acid glycoprotein (AGP). The extrinsic Cotton effects for the binding of aripiprazole and its derivatives to AGP were generated, but the magnitudes of the induced circular dichroism intensities did not correlate with those for the binding affinities. It therefore appears that the binding mode of aripiprazole with AGP is somewhat complicated, compared with that of albumin. Isothermal titration calorimetry data obtained for the binding of aripiprazole with AGP were different from that for albumin systems in that the 3 driving reactions, entropy-driven, enthalpy-driven, and the entropy-enthalpy mixed type, were all found for the AGP system, but not albumin. Moreover, the weak binding mode of aripiprazole with the 2 proteins were supported by a molecular docking model analysis. The concentration of albumin in plasma is about 50 times higher than those of AGP, but AGP levels in plasma are increased by about 10 times under inflammatory disease. Therefore, the involvement of these 2 plasma proteins should be considered in more depth for understanding the pharmacokinetics of aripiprazole.

Further Evidence Regarding the Important Role of Chlorine Atoms of Aripiprazole on Binding to the Site II Area of Human Albumin

Previously, we reported on the high-affinity binding of aripiprazole (ARP), an antipsychotic drug, to human albumin and the role of the chlorine atom of ARP on this binding. In this study, we investigated the binding mode of ARP to human albumin in detail using ARP derivatives and several animal-derived albumins. ARP bound strongly to human and dog albumin. The circular dichroism (CD) spectra of ARP bound to human and dog albumin were also similar. Deschloro-ARP bound less strongly to all of the albumin species compared to ARP, and the shapes of CD spectra were similar for all albumin species. CD spectra of dimethyl-ARP, for which chlorine atoms were substituted methyl groups, were quite similar to that of deschloro-ARP. In displacement experiments, competitive binding was observed between ARP and deschloro-ARP. These results suggest that the chlorine atoms in ARP are involved in the binding modes of ARP for human and dog albumins, whereas ARP and deschloro-ARP appear to share the same binding region in site II. The aforementioned results imply that compounds having a chlorine atom bind more strongly to plasma proteins, resulting in a long blood retention time. Therefore, findings reported here may provide the basically useful data for drug design.

Dynamics of human protein kinase Aurora A linked to drug selectivity

Protein kinases are major drug targets, but the development of highly-selective inhibitors has been challenging due to the similarity of their active sites. The observation of distinct structural states of the fully-conserved Asp-Phe-Gly (DFG) loop has put the concept of conformational selection for the DFG-state at the center of kinase drug discovery. Recently, it was shown that Gleevec selectivity for the Tyr-kinase Abl was instead rooted in conformational changes after drug binding. Here, we investigate whether protein dynamics after binding is a more general paradigm for drug selectivity by characterizing the binding of several approved drugs to the Ser/Thr-kinase Aurora A. Using a combination of biophysical techniques, we propose a universal drug-binding mechanism, that rationalizes selectivity, affinity and long on-target residence time for kinase inhibitors. These new concepts, where protein dynamics in the drug-bound state plays the crucial role, can be applied to inhibitor design of targets outside the kinome.

Protein kinases are a family of enzymes found in all living organisms. These enzymes help to control many biological processes, including cell division. When particular protein kinases do not work correctly, cells may start to divide uncontrollably, which can lead to cancer. One example is the kinase Aurora A, which is over-active in many common human cancers. As a result, researchers are currently trying to design drugs that reduce the activity of Aurora A in the hope that these could form new anticancer treatments.

In general, drugs are designed to be as specific in their action as possible to reduce the risk of harmful side effects to the patient. Designing a drug that affects a single protein kinase, however, is difficult because there are hundreds of different kinases in the body, all with similar structures. Because drugs often work by binding to specific structural features, a drug that targets one protein kinase can often alter the activity of a large number of others too.

Gleevec is a successful anti-leukemia drug that specifically works on one target kinase, producing minimal side effects. It was recently discovered that the drug works through a phenomenon called ‘induced fit’. This means that after the drug binds it causes a change in the enzyme’s overall shape that alters the activity of the enzyme. The shape change is complex, and so even small structural differences can change the effect of a particular drug.

Do other drugs that target other protein kinases also produce induced fit effects? To find out, Pitsawong, Buosi, Otten, Agafonov et al. studied how three anti-cancer drugs interact with Aurora A: two drugs specifically designed to switch off Aurora A, and Gleevec (which does not target Aurora A).

The two drugs that specifically target Aurora A were thought to work by targeting one structural feature of the enzyme. However, the biochemical and biophysical experiments performed by Pitsawong et al. revealed that these drugs instead work through an induced fit effect. By contrast, Gleevec did not trigger an induced fit on Aurora A and so bound less tightly to it.

In light of these results, Pitsawong et al. suggest that future efforts to design drugs that target protein kinases should focus on exploiting the induced fit process. This will require more research into the structure of particular kinases.

Binding of Small Molecule Drugs to Porcine Vitreous Humor

Pharmacokinetics in the posterior eye segment has therapeutic implications due to the importance of retinal diseases in ophthalmology. In principle, drug binding to the components of the vitreous, such as proteins, collagen, or glycosaminoglycans, could prolong ocular drug retention and modify levels of pharmacologically active free drug in the posterior eye segment. Since drug binding in the vitreous has been investigated only sparsely, we studied vitreal drug binding of 35 clinical small molecule drugs. Isolated homogenized porcine vitreous and the drugs were placed in a two-compartment dialysis system that was used to separate the bound and unbound drug. Free drug concentrations and binding percentages were quantitated using LC-MS/MS. Drug binding levels varied between 21 and 74% in the fresh vitreous and 0 and 64% in the frozen vitreous. The vitreal binding percentages did not correlate with those in plasma. Our data-based pharmacokinetic simulations suggest that vitreal binding of small molecule drugs has only a modest influence on the AUC of free drug or drug half-life in the vitreous. Therefore, it is likely that vitreal binding is not a major reason for interindividual variability in ocular drug responses or drug-drug interactions.

"To Be or Not to Be" Protonated: Atomic Details of Human Carbonic Anhydrase-Clinical Drug Complexes by Neutron Crystallography and Simulation

Human carbonic anhydrases (hCAs) play various roles in cells, and have been drug targets for decades. Sequence similarities of hCA isoforms necessitate designing specific inhibitors, which requires detailed structural information for hCA-inhibitor complexes. We present room temperature neutron structures of hCA II in complex with three clinical drugs that provide in-depth analysis of drug binding, including protonation states of the inhibitors, hydration water structure, and direct visualization of hydrogen-bonding networks in the enzyme's active site. All sulfonamide inhibitors studied bind to the Zn metal center in the deprotonated, anionic, form. Other chemical groups of the drugs can remain neutral or be protonated when bound to hCA II. MD simulations have shown that flexible functional groups of the inhibitors may alter their conformations at room temperature and occupy different sub-sites. This study offers insights into the design of specific drugs to target cancer-related hCA isoform IX.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

2-(m-Azidobenzoyl)taxol binds differentially to distinct β-tubulin isotypes

There are seven β-tubulin isotypes present in distinct quantities in mammalian cells of different origin. Altered expression of β-tubulin isotypes has been reported in cancer cell lines resistant to microtubule stabilizing agents (MSAs) and in human tumors resistant to Taxol. To study the relative binding affinities of MSAs, tubulin from different sources, with distinct β-tubulin isotype content, were specifically photolabeled with a tritium-labeled Taxol analog, 2-(m-azidobenzoyl)taxol, alone or in the presence of MSAs. The inhibitory effects elicited by these MSAs on photolabeling were distinct for β-tubulin from different sources. To determine the exact amount of drug that binds to different β-tubulin isotypes, bovine brain tubulin was photolabeled and the isotypes resolved by high-resolution isoelectrofocusing. All bands were analyzed by mass spectrometry following cyanogen bromide digestion, and the identity and relative quantity of each β-tubulin isotype determined. It was found that compared with other β-tubulin isotypes, βIII-tubulin bound the least amount of 2-(m-azidobenzoyl)taxol. Analysis of the sequences of β-tubulin near the Taxol binding site indicated that, in addition to the M-loop that is known to be involved in drug binding, the leucine cluster region of βIII-tubulin contains a unique residue, alanine, at 218, compared with other isotypes that contain threonine. Molecular dynamic simulations indicated that the frequency of Taxol-accommodating conformations decreased dramatically in the T218A variant, compared with other β-tubulins. Our results indicate that the difference in residue 218 in βIII-tubulin may be responsible for inhibition of drug binding to this isotype, which could influence downstream cellular events.

Neutron structure of human carbonic anhydrase II in complex with methazolamide: mapping the solvent and hydrogen-bonding patterns of an effective clinical drug

Carbonic anhydrases (CAs; EC 4.2.1.1) catalyze the interconversion of CO2 and HCO3-, and their inhibitors have long been used as diuretics and as a therapeutic treatment for many disorders such as glaucoma and epilepsy. Acetazolamide (AZM) and methazolamide (MZM, a methyl derivative of AZM) are two of the classical CA inhibitory drugs that have been used clinically for decades. The jointly refined X-ray/neutron structure of MZM in complex with human CA isoform II (hCA II) has been determined to a resolution of 2.2 Å with an Rcryst of ∼16.0%. Presented in this article, along with only the second neutron structure of a clinical drug-bound hCA, is an in-depth structural comparison and analyses of differences in hydrogen-bonding network, water-molecule orientation and solvent displacement that take place upon the binding of AZM and MZM in the active site of hCA II. Even though MZM is slightly more hydrophobic and displaces more waters than AZM, the overall binding affinity (Ki) for both of the drugs against hCA II is similar (∼10 nM). The plausible reasons behind this finding have also been discussed using molecular dynamics and X-ray crystal structures of hCA II-MZM determined at cryotemperature and room temperature. This study not only allows a direct comparison of the hydrogen bonding, protonation states and solvent orientation/displacement of AZM and MZM, but also shows the significant effect that the methyl derivative has on the solvent organization in the hCA II active site.

Differential efficacy of GoSlo-SR compounds on BKα and BKαγ1-4 channels

Large conductance, voltage and Ca²⁺ activated K⁺ channels (BK channels) are abundantly expressed throughout the body and are important regulators of smooth muscle tone and neuronal excitability. Their dysfunction is implicated in various diseases including overactive bladder, hypertension and erectile dysfunction. Therefore, BK channel openers bear significant therapeutic potential to treat the above diseases. GoSlo-SR compounds were designed to be potent and efficacious BK channel openers. Although their structural activity relationships, activation in both BKα and BKαβ channels and the hypothetical mode of action of these compounds has been studied in detail in recent years, their effectiveness to open the BKαγ channels still remains unexplored. In this study, we have examined the efficacy of 3 closely related GoSlo-SR openers, GoSlo-SR-5-6 (SR-5-6), GoSlo-SR-5-44 (SR-5-44) and GoSlo-SR-5-130 (SR-5-130) using inside out patches on BKα channels coexpressed with 4 different LRRC (γ_1-4) subunits in HEK293 cells. Our data suggests that the activation effects due to SR-5-6 were not significantly affected in the presence of γ_1-4 subunits. Interestingly, the effects of more efficacious BK channel opener SR-5-44 were altered by different γ subunits. In cells expressing BKα channels, the shift in V_1/2 (ΔV_1/2) induced by SR-5-44 (3 μM) was -76 ± 3 mV, whereas it was significantly reduced by ∼70 % in BKαγ₁ channels (ΔV_1/2= -23 ± 3, p < 0.001, ANOVA). In BKαγ₂ channels the ΔV_1/2 was -36 ± 1 mV, which was less than that observed in BKαγ₃ and BKαγ₄ channels where the ΔV_1/2 was -47 ± 5 mV, and -82 ± 5 mV, respectively. Additionally, the excitatory effects of a 'β specific' BK channel opener, SR-5-130 were only partially restored in the patches containing BKαγ_1-4 channels. Together this data highlights that subtle modifications in GoSlo-SR structures alter their effectiveness on BK channels with accessory γ subunits and this study might provide a scaffold for the development of more tissue specific BK channel openers.

Elucidation of structural and functional properties of albumin bound to gold nanoparticles

Nanoparticle-albumin complexes are being designed for targeted drug delivery and imaging. However, the changes in the functional properties of albumin due to adsorption on nanoparticles remain elusive. Thus, the objective of this work was to elucidate the structural and functional properties of human and bovine serum albumin bound to negatively charged gold nanoparticles (GNPs). Fluorescence data demonstrated static quenching of albumin by GNP with the quenching of buried as well as surface tryptophan in BSA. The binding process was enthalpy and entropy-driven in HSA and BSA, respectively. At lower concentrations of GNP there was a higher affinity for tryptophan, whereas at higher concentrations both tryptophan and tyrosine participated in the interaction. Synchronous fluorescence spectra revealed that the microenvironment of tryptophan in HSA turned more hydrophilic upon exposure to GNP. The α-helical content of albumin was unaltered by GNP. Approximately 37 and 23% reduction in specific activity of HSA and BSA was observed due to GNP binding. In presence of warfarin and ibuprofen the binding constants of albumin-GNP complexes were altered. A very interesting observation not reported so far is the retained antioxidant activity of albumin in presence of GNP i.e. we believe that GNPs did not bind to the free sulfhydryl groups of albumin. However enhanced levels of copper binding were observed. We have also highlighted the differential response in albumin due to gold and silver nanoparticles which could be attributed to differences in the charge of the nanoparticle.

A functional NMR for membrane proteins: dynamics, ligand binding, and allosteric modulation

By nature of conducting ions, transporting substrates and transducing signals, membrane channels, transporters and receptors are expected to exhibit intrinsic conformational dynamics. It is therefore of great interest and importance to understand the various properties of conformational dynamics acquired by these proteins, for example, the relative population of states, exchange rate, conformations of multiple states, and how small molecule ligands modulate the conformational exchange. Because small molecule binding to membrane proteins can be weak and/or dynamic, structural characterization of these effects is very challenging. This review describes several NMR studies of membrane protein dynamics, ligand-induced conformational rearrangements, and the effect of ligand binding on the equilibrium of conformational exchange. The functional significance of the observed phenomena is discussed.

Mutant bacterial sodium channels as models for local anesthetic block of eukaryotic proteins

Voltage gated sodium channels are the target of a range of local anesthetic, anti-epileptic and anti-arrhythmic compounds. But, gaining a molecular level understanding of their mode of action is difficult as we only have atomic resolution structures of bacterial sodium channels not their eukaryotic counterparts. In this study we used molecular dynamics simulations to demonstrate that the binding sites of both the local anesthetic benzocaine and the anti-epileptic phenytoin to the bacterial sodium channel NavAb can be altered significantly by the introduction of point mutations. Free energy techniques were applied to show that increased aromaticity in the pore of the channel, used to emulate the aromatic residues observed in eukaryotic Nav1.2, led to changes in the location of binding and dissociation constants of each drug relative to wild type NavAb. Further, binding locations and dissociation constants obtained for both benzocaine (660 μM) and phenytoin (1 μM) in the mutant channels were within the range expected from experimental values obtained from drug binding to eukaryotic sodium channels, indicating that these mutant NavAb may be a better model for drug binding to eukaryotic channels than the wild type.

Identification of residues in ABCG2 affecting protein trafficking and drug transport, using co-evolutionary analysis of ABCG sequences

ABCG2 is an ABC (ATP-binding cassette) transporter with a physiological role in urate transport in the kidney and is also implicated in multi-drug efflux from a number of organs in the body. The trafficking of the protein and the mechanism by which it recognizes and transports diverse drugs are important areas of research. In the current study, we have made a series of single amino acid mutations in ABCG2 on the basis of sequence analysis. Mutant isoforms were characterized for cell surface expression and function. One mutant (I573A) showed disrupted glycosylation and reduced trafficking kinetics. In contrast with many ABC transporter folding mutations which appear to be 'rescued' by chemical chaperones or low temperature incubation, the I573A mutation was not enriched at the cell surface by either treatment, with the majority of the protein being retained in the endoplasmic reticulum (ER). Two other mutations (P485A and M549A) showed distinct effects on transport of ABCG2 substrates reinforcing the role of TM helix 3 in drug recognition and transport and indicating the presence of intracellular coupling regions in ABCG2.

Structure of the mycobacterial ATP synthase Fo rotor ring in complex with the anti-TB drug bedaquiline

Multidrug-resistant tuberculosis (MDR-TB) is more prevalent today than at any other time in human history. Bedaquiline (BDQ), a novel Mycobacterium-specific adenosine triphosphate (ATP) synthase inhibitor, is the first drug in the last 40 years to be approved for the treatment of MDR-TB. This bactericidal compound targets the membrane-embedded rotor (c-ring) of the mycobacterial ATP synthase, a key metabolic enzyme required for ATP generation. We report the x-ray crystal structures of a mycobacterial c9 ring without and with BDQ bound at 1.55- and 1.7-Å resolution, respectively. The structures and supporting functional assays reveal how BDQ specifically interacts with the rotor ring via numerous interactions and thereby completely covers the c-ring's ion-binding sites. This prevents the rotor ring from acting as an ion shuttle and stalls ATP synthase operation. The structures explain how diarylquinoline chemicals specifically inhibit the mycobacterial ATP synthase and thus enable structure-based drug design of next-generation ATP synthase inhibitors against Mycobacterium tuberculosis and other bacterial pathogens.

Local anesthetic and antiepileptic drug access and binding to a bacterial voltage-gated sodium channel

Voltage-gated sodium (Nav) channels are important targets in the treatment of a range of pathologies. Bacterial channels, for which crystal structures have been solved, exhibit modulation by local anesthetic and anti-epileptic agents, allowing molecular-level investigations into sodium channel-drug interactions. These structures reveal no basis for the "hinged lid"-based fast inactivation, seen in eukaryotic Nav channels. Thus, they enable examination of potential mechanisms of use- or state-dependent drug action based on activation gating, or slower pore-based inactivation processes. Multimicrosecond simulations of NavAb reveal high-affinity binding of benzocaine to F203 that is a surrogate for FS6, conserved in helix S6 of Domain IV of mammalian sodium channels, as well as low-affinity sites suggested to stabilize different states of the channel. Phenytoin exhibits a different binding distribution owing to preferential interactions at the membrane and water-protein interfaces. Two drug-access pathways into the pore are observed: via lateral fenestrations connecting to the membrane lipid phase, as well as via an aqueous pathway through the intracellular activation gate, despite being closed. These observations provide insight into drug modulation that will guide further developments of Nav inhibitors.

Mathematical and computational models of drug transport in tumours

The ability to predict how far a drug will penetrate into the tumour microenvironment within its pharmacokinetic (PK) lifespan would provide valuable information about therapeutic response. As the PK profile is directly related to the route and schedule of drug administration, an in silico tool that can predict the drug administration schedule that results in optimal drug delivery to tumours would streamline clinical trial design. This paper investigates the application of mathematical and computational modelling techniques to help improve our understanding of the fundamental mechanisms underlying drug delivery, and compares the performance of a simple model with more complex approaches. Three models of drug transport are developed, all based on the same drug binding model and parametrized by bespoke in vitro experiments. Their predictions, compared for a 'tumour cord' geometry, are qualitatively and quantitatively similar. We assess the effect of varying the PK profile of the supplied drug, and the binding affinity of the drug to tumour cells, on the concentration of drug reaching cells and the accumulated exposure of cells to drug at arbitrary distances from a supplying blood vessel. This is a contribution towards developing a useful drug transport modelling tool for informing strategies for the treatment of tumour cells which are 'pharmacokinetically resistant' to chemotherapeutic strategies.

Complex Interplay between the P-Glycoprotein Multidrug Efflux Pump and the Membrane: Its Role in Modulating Protein Function

Multidrug resistance in cancer is linked to expression of the P-glycoprotein multidrug transporter (Pgp, ABCB1), which exports many structurally diverse compounds from cells. Substrates first partition into the bilayer and then interact with a large flexible binding pocket within the transporter's transmembrane regions. Pgp has been described as a hydrophobic vacuum cleaner or an outwardly directed drug/lipid flippase. Recent X-ray crystal structures have shed some light on the nature of the drug-binding pocket and suggested routes by which substrates can enter it from the membrane. Detergents have profound effects on Pgp function, and several appear to be substrates. Biochemical and biophysical studies in vitro, some using purified reconstituted protein, have explored the effects of the membrane environment. They have demonstrated that Pgp is involved in a complex relationship with its lipid environment, which modulates the behavior of its substrates, as well as various functions of the protein, including ATP hydrolysis, drug binding, and drug transport. Membrane lipid composition and fluidity, phospholipid headgroup and acyl chain length all influence Pgp function. Recent studies focusing on thermodynamics and kinetics have revealed some important principles governing Pgp-lipid and substrate-lipid interactions, and how these affect drug-binding and transport. In some cells, Pgp is associated with cholesterol-rich microdomains, which may modulate its functions. The relationship between Pgp and cholesterol remains an open question; however, it clearly affects several aspects of its function in addition to substrate-membrane partitioning. The action of Pgp modulators appears to depend on their membrane permeability, and membrane fluidizers and surfactants reverse drug resistance, likely via an indirect mechanism. A detailed understanding of how the membrane affects Pgp substrates and Pgp's catalytic cycle may lead to new strategies to combat clinical drug resistance.

State-Dependent Inhibition of Sodium Channels by Local Anesthetics: A 40-Year Evolution

Knowledge about the mechanism of impulse blockade by local anesthetics has evolved over the past four decades, from the realization that Na⁺ channels were inhibited to affect the impulse blockade to an identification of the amino acid residues within the Na⁺ channel that bind the local anesthetic molecule. Within this period appreciation has grown of the state-dependent nature of channel inhibition, with rapid binding and unbinding at relatively high affinity to the open state, and weaker binding to the closed resting state. Slow binding of high affinity for the inactivated state accounts for the salutary therapeutic as well as the toxic actions of diverse class I anti-arrhythmic agents, but may have little importance for impulse blockade, which requires concentrations high enough to block the resting state. At the molecular level, residues on the S6 transmembrane segments in three of the homologous domains of the channel appear to contribute to the binding of local anesthetics, with some contribution also from parts of the selectivity filter. Binding to the inactivated state, and perhaps the open state, involves some residues that are not identical to those that bind these drugs in the resting state, suggesting spatial flexibility in the "binding site". Questions remaining include the mechanism that links local anesthetic binding with the inhibition of gating charge movements, and the molecular nature of the theoretical "hydrophobic pathway" that may be critical for determining the recovery rates from blockade of closed channels, and thus account for both therapeutic and cardiotoxic actions.

Missing fragments: detecting cooperative binding in fragment-based drug design

The aim of fragment-based drug design (FBDD) is to identify molecular fragments that bind to alternate subsites within a given binding pocket leading to cooperative binding when linked. In this study, the binding of fragments to human phenylethanolamine N-methyltransferase is used to illustrate how (a) current protocols may fail to detect fragments that bind cooperatively, (b) theoretical approaches can be used to validate potential hits, and (c) apparent false positives obtained when screening against cocktails of fragments may in fact indicate promising leads.

Photoaffinity labeling of the multidrug resistance protein 2 (ABCC2/cMOAT) with a photoreactive analog of LTC(4)

Several studies have shown that the multidrug resistant protein MRP2 mediates the transport of chemotherapeutic drugs and normal cell metabolites, including Leukotriene C (LTC(4)); however direct binding of the LTC(4) to MRP2 has not been demonstrated. In this study, a photoreactive analog of LTC(4) (IAALTC(4)) was used to demonstrate its direct binding to MRP2. Our results show specific photoaffinity labeling of MRP2 with IAALTC(4) in plasma membranes from MDCKII(MRP2) cells. The photoaffinity labeling signal of MRP2 with IAALTC(4) was much lower than that of MRP1, consistent with previous studies whereby the measured K(m) values of MRP1 and MRP2 for LTC(4) were 1 μM and 0.1 μM LTC(4), respectively. Competition of IAALTC(4) photoaffinity labeling to MRP2 with MK571, a well characterized inhibitor of MRP2 function, showed ~75% reduction in binding in the presence of 50 μM excess MK571. Interestingly, unmodified LTC(4) enhanced the photoaffinity labeling of IAALTC(4) to MRP2, whereas excess GSH and Quercetin had no significant effect. Mild tryptic digestion of photoaffinity labeled MRP2 revealed several photoaffinity labeled peptides that localized the IAALTC(4) binding to a 15 kDa amino acid sequence that contains transmembrane 16 and 17. Together these results provide the first demonstration of direct LTC(4) binding to MRP2.

Protein structures by spallation neutron crystallography

The Protein Crystallography Station at Los Alamos Neutron Science Center is a high-performance beamline that forms the core of a capability for neutron macromolecular structure and function determination. This capability also includes the Macromolecular Neutron Crystallography (MNC) consortium between Los Alamos (LANL) and Lawrence Berkeley National Laboratories for developing computational tools for neutron protein crystallography, a biological deuteration laboratory, the National Stable Isotope Production Facility, and an MNC drug design consortium between LANL and Case Western Reserve University.

Detailed characterization of cysteine-less P-glycoprotein reveals subtle pharmacological differences in function from wild-type protein

1. Subtle alterations in the coupling of drug binding to nucleotide hydrolysis were observed following mutation of all seven endogenous cysteine residues to serines in the human multidrug resistance transporter, P-glycoprotein. Wild-type (wt) and the mutant (cys-less) forms of P-gp were expressed in Trichoplusia ni (High Five) cells and purified by metal affinity chromatography in order to undertake functional studies. 2. No significant differences were observed in substrate ([(3)H]-azidopine) binding to wt or cys-less P-gp. Furthermore, neither the transported substrate vinblastine, nor the modulator nicardipine, differed in their respective potencies to displace [(3)H]-azidopine from the wt or cys-less P-gp. These results suggest that respective binding sites for these drugs were unaffected by the introduced cysteine to serine substitutions. 3. The Michaelis-Menten characteristics of basal ATP hydrolysis of the two isoforms of P-gp were identical. The maximal ATPase activity in the presence of vinblastine was marginally reduced whilst the K(m) was unchanged in cys-less P-gp compared to control. However, cys-less P-gp displayed lower overall maximal ATPase activity (62%), a decreased K(m) and a lower degree of stimulation (76%) in the presence of the modulator nicardipine. 4. Therefore, the serine to cysteine mutations in P-gp may suggest that vinblastine and nicardipine transduce their effects on ATP hydrolysis through distinct conformational pathways. The wt and cys-less P-gp isoforms display similarity in their fundamental kinetic properties thereby validating the use of cys-less P-gp as a template for future cysteine-directed structure/function analysis.

The equilibrium and kinetic drug binding properties of the mouse P-gp1a and P-gp1b P-glycoproteins are similar

The gene encoding the multidrug resistance P-glycoprotein (P-gp) is duplicated in rodent species and the functional basis for this remains unresolved. Despite a high sequence similarity, the mouse P-gp1a and P-gp1b isoforms show distinct patterns of tissue distribution which suggest a specific role of the P-gp1b isoform in steroid transport. In the present study possible biochemical differences between the isoforms were directly investigated at the level of drug interaction. There was no detectable difference in the affinity or binding capacity of the two isoforms towards [3H]vinblastine at equilibrium. Similarly, the rate at which [3H]vinblastine associates with P-gp was indistinguishable between the two isoforms. Some modest differences were observed in the relative abilities of the multidrug-resistant (MDR) reversing agents CP100-356, nicardipine and verapamil to displace equilibrium [3H]vinblastine binding to P-gp1a and P-gp1b. The steroid hormone progesterone displayed a low affinity (Ki = 1.2 +/- 0.2 microM for P-gp1a and 3.5 +/- 0.5 microM for P-gp1b), suggesting an unlikely role as a physiological substrate. Thus the mouse isoforms do not appear to exhibit functional differences at the level of initial substrate interaction with protein.

The multi-drug resistance reversal agent SR33557 and modulation of vinca alkaloid binding to P-glycoprotein by an allosteric interaction

1. The interaction of the indolizin sulfone SR33557 with the multidrug resistance P-glycoprotein (P-gp), was used to explore the nature of drug binding site(s) on this transporter. The steady-state accumulation of [3H]-vinblastine in P-gp expressing CHrB30 cells was increased by SR33557 with greater potency than verapamil. Furthermore, SR33557 potentiated the affinity of verapamil to modulate vinblastine transport when added simultaneously. 2. Verapamil elicited a 1.5 to 2.5 fold stimulation of basal ATPase activity in CHrB30 membranes, whereas SR33557 and vinblastine inhibited activity, but only at relatively high concentrations. However, SR33557 and vinblastine decreased the Vmax but not the Km for verapamil stimulation of ATPase activity. This is indicative of a non-competitive interaction, most likely at distinct sites. 3. The specific [3H]-vinblastine binding to P-gp in CHrB30 cell membranes was displaced by SR33557 with an IC50 of 8.3 +/- 4.5 nM. Moreover, SR33557 caused a 3 fold increase in the dissociation rate of vinblastine binding to P-gp indicating a negative allosteric effect on the vinca alkaloid acceptor site. 4. These results demonstrate that SR33557 interacts with a site on P-gp which is distinct from, but allosterically linked to the vinca alkaloid site. The apparent broad substrate specificity displayed by P-gp may be explained by a multiple drug binding site model.

Addition of colchicine--tubulin complex to microtubule ends: the mechanism of substoichiometric colchicine poisoning

Colchicine blocks microtubule polymerization by an unusual substoichiometric poisoning mechanism. We have investigated the mechanism by which this poisoning occurs with several experimental approaches, and have found that colchicine acts by addition to microtubule ends as a colchicine-tubulin complex.

Ethidium bromide as a cooperative effector of a DNA structure

A salt-induced cooperative conformational transition of a synthetic DNA, poly(dG-dC), is reversed by addition of ethidium bromide. Binding of the dye at high salt concentrations is highly cooperative. Circular dichroism spectra of the complex and the kinetic data support a model for this cooperative binding that is formally equivalent to the "allosteric" one proposed for oligomeric proteins by Monod et al. Thus, double-helical DNA of at least one defined sequence can undergo a cooperative conformational change in solution, with simple salts and drug molecules as antagonistic effectors. Such transitions may be involved in regulatory phenomena operating directly at the level of nucleic acid structure.