PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71Linked to original sources

Chlorella virus-encoded deoxyuridine triphosphatases exhibit different temperature optima.

A putative deoxyuridine triphosphatase (dUTPase) gene from chlorella virus PBCV-1 was cloned, and the recombinant protein was expressed in Escherichia coli. The recombinant protein has dUTPase activity and requires Mg(2+) for optimal activity, while it retains some activity in the presence of other divalent cations. Kinetic studies of the enzyme revealed a K(m) of 11.7 microM, a turnover k(cat) of 6.8 s(-1), and a catalytic efficiency of k(cat)/K(m) = 5.8 x 10(5) M(-1) s(-1). dUTPase genes were cloned and expressed from two other chlorella viruses IL-3A and SH-6A. The two dUTPases have similar properties to PBCV-1 dUTPase except that IL-3A dUTPase has a lower temperature optimum (37 degrees C) than PBCV-1 dUTPase (50 degrees C). The IL-3A dUTPase differs from the PBCV-1 enzyme by nine amino acids, including two amino acid substitutions, Glu81-->Ser81 and Thr84-->Arg84, in the highly conserved motif III of the proteins. To investigate the difference in temperature optima between the two enzymes, homology modeling and docking simulations were conducted. The results of the simulation and comparisons of amino acid sequence suggest that adjacent amino acids are important in the temperature optima. To confirm this suggestion, three site-directed amino acid substitutions were made in the IL-3A enzyme: Thr84-->Arg84, Glu81-->Ser81, and Glu81-->Ser81 plus Thr84-->Arg84. The single substitutions affected the optimal temperature for enzyme activity. The temperature optimum increased from 37 to 55 degrees C for the enzyme containing the two amino acid substitutions. We postulate that the change in temperature optimum is due to reduction in charge and balkiness in the active cavity that allows more movement of the ligand and protein before the enzyme and substrate complex is formed.

Amino Acid Motifs↗

Identification of novel inhibitors of BCR-ABL tyrosine kinase via virtual screening.

Inhibition of BCR-ABL tyrosine kinase activity has shown to be essential for the treatment of chronic myelogenous leukemia (CML). However, drug resistance has quickly arisen in recent clinical trials for STI571 (Gleevec), which is the first approved drug of CML by inhibiting ABL tyrosine kinase. It is desirable to develop new types of ABL tyrosine kinase inhibitors that may overcome this drug resistance problem. Here we present the discovery of novel inhibitors targeted at the catalytic domain of ABL tyrosine kinase by using three-dimensional database searching techniques. From a database containing 200,000 commercially available compounds, the top 1000 compounds with the best DOCK energy score were selected and subjected to structural diversity and drug likeness analysis, 15 compounds were submitted for biological assay. Eight out of the 15 showed inhibitory activity against K562 cells with IC(50) value ranging from 10 to 200 microM. Two promising compounds showed inhibition in further ABL tyrosine phosphorylation assay. It is anticipated that those two compounds can serve as lead compounds for further drug design and optimization.

Binding Sites↗

Optimised amino acid specific weighting factors for unbound protein docking.

BACKGROUND: One of the most challenging aspects of protein-protein docking is the inclusion of flexibility into the docking procedure. We developed a postfilter where the grid-representation of proteins for docking is extended by an optimised weighting factor for each amino acid. RESULTS: For up to 86% of the evaluated complexes a near-native structure was within the top 5% of the ranked prediction output. The weighting factors obtained by the optimisation procedure correlate to a certain extent with the flexibility of the amino acids, their hydrophobicity and with their propensity to be in the interface. CONCLUSION: Use of the optimised amino acid specific parameters yields a strong increase of near-native structures on the first ranks of the prediction.

Algorithms↗

Structure of a conjugating enzyme-ubiquitin thiolester intermediate reveals a novel role for the ubiquitin tail.

BACKGROUND: Ubiquitin-conjugating enzymes (E2s) are central enzymes involved in ubiquitin-mediated protein degradation. During this process, ubiquitin (Ub) and the E2 protein form an unstable E2-Ub thiolester intermediate prior to the transfer of ubiquitin to an E3-ligase protein and the labeling of a substrate for degradation. A series of complex interactions occur among the target substrate, ubiquitin, E2, and E3 in order to efficiently facilitate the transfer of the ubiquitin molecule. However, due to the inherent instability of the E2-Ub thiolester, the structural details of this complex intermediate are not known. RESULTS: A three-dimensional model of the E2-Ub thiolester intermediate has been determined for the catalytic domain of the E2 protein Ubc1 (Ubc1(Delta450)) and ubiquitin from S. cerevisiae. The interface of the E2-Ub intermediate was determined by kinetically monitoring thiolester formation by 1H-(15)N HSQC spectra by using combinations of 15N-labeled and unlabeled Ubc1(Delta450) and Ub proteins. By using the surface interface as a guide and the X-ray structures of Ub and the 1.9 A structure of Ubc1(Delta450) determined here, docking simulations followed by energy minimization were used to produce the first model of a E2-Ub thiolester intermediate. CONCLUSIONS: Complementary surfaces were found on the E2 and Ub proteins whereby the C terminus of Ub wraps around the E2 protein terminating in the thiolester between C88 (Ubc1(Delta450)) and G76 (Ub). The model supports in vivo and in vitro experiments of E2 derivatives carrying surface residue substitutions. Furthermore, the model provides insights into the arrangement of Ub, E2, and E3 within a ternary targeting complex.

Computer Simulation↗

Bioactive conformation of a potent stromelysin inhibitor determined by X-nucleus filtered and multidimensional NMR spectroscopy.

The biologically active conformation of a novel, very potent, nonpeptidic stromelysin inhibitor was determined by X-nucleus filtered and multidimensional NMR spectroscopy. This bound conformer was subsequently docked into the stromelysin catalytic domain (SCD) using intermolecular distance constraints derived from NOE data. The complex showed the S1' pocket of stromelysin to be the major site of enzyme-inhibitor interaction with other portions of the inhibitor spanning the S2' and S1 binding sites. Theoretical predictions of SCD-inhibitor binding from molecular modeling studies were consistent with the NMR data. Comparison of modeled enzyme-inhibitor complexes for stromelysin and collagenase revealed an alternate binding mode for the inhibitor in collagenase, suggesting a similar binding interaction might also be possible for stromelysin. The NMR results, however, revealed a single SCD-inhibitor binding mode and provided a structural template for the design of more potent stromelysin inhibitors.

Binding Sites↗

The molecular structure of wild-type and a mutant Fis protein: relationship between mutational changes and recombinational enhancer function or DNA binding.

The 98-amino acid Fis protein from Escherichia coli functions in a variety of reactions, including promotion of Hin-mediated site-specific DNA inversion when bound to an enhancer sequence. It is unique among site-specific DNA-binding proteins in that it binds to a large number of different DNA sequences, for which a consensus sequence is difficult to establish. X-ray crystal structure analyses have been carried out at 2.3 A resolution for wild-type Fis and for an Arg-89----Cys mutant that does not stimulate DNA inversion. Each monomer of the Fis dimer has four alpha-helices, A-D; the first 19 residues are disordered in the crystal. The end of each C helix is hydrogen bonded to the beginning of helix B' from the opposite subunit in what effectively is one long continuous, although bent, helix. The four helices, C, B', C', and B, together define a platform through the center of the Fis molecule: helices A and A' are believed to be involved with Hin recombinase on one side, and helices D and D' interact with DNA lying on the other side of the platform. Helices C and D of each subunit comprise a helix-turn-helix (HTH) DNA-binding element. The spacing of these two HTH elements in the dimer, 25 A, is too short to allow insertion into adjacent major grooves of a straight B-DNA helix. However, bending the DNA at discrete points, to an overall radius of curvature of 62 A, allows efficient docking of a B-DNA helix with the Fis molecule. The proposed complex explains the experimentally observed patterns of methylation protection and DNase I cleavage hypersensitivity. The x-ray structure accounts for the effects of mutations in the Fis sequence. Those that affect DNA inversion but not DNA binding are located within the N-terminal disordered region and helix A. This inversion activation domain is physically separated in the Fis molecule from the HTH elements and may specify a region of contact with the Hin recombinase. In contrast, mutations that affect HTH helices C and D, or interactions of these with helix B, have the additional effect of decreasing or eliminating binding to DNA.

Amino Acid Sequence↗

Application of the in combo screening approach for the discovery of non-alkaloid acetylcholinesterase inhibitors from Cichorium intybus.

Because of the direct correlation of cholinergic deficit and the severity of dementia, Alzheimer's disease is preferentially treated with acetylcholinesterase (AChE) inhibitors to supplement the acetylcholine level. In this study we focused on non-alkaloid AChE inhibitors from natural sources in order to discover new lead structures. In the course of in vitro extract screening of Tyrolean plants using an enzyme assay with Ellman's reagent, the dichloromethane extract of chicory roots (Cichorium intybus L.) showed a pronounced inhibitory effect on AChE. At a concentration of 1 mg extract/ml an inhibition of 70% was measured. Based on a 3D multi-conformational molecular-structure database consisting of secondary metabolites from C. intybus known from the relevant literature, virtual screening filtering experiments were conducted using both a feature-based pharmacophore model and a docking procedure. Some low molecular weight sesquiterpenoids exhibited distinct interactions with the pharmacophore model. In order to verify the applicability of this computer-aided strategy, an activity-guided fractionation of the chicory root extract was performed, which resulted in the isolation of two sesquiterpene lactones, 8-deoxylactucin and lactucopicrin, showing significant and dose-dependent inhibitory activity on AChE (IC(50) of 308.1 microM [CI(95) 243.9 - 405.3 microM] and 150.3 microM [CI(95) 100.8 - 188.1 microM], respectively). The two isolates were correctly predicted within the virtual screening process which corroborates the potential of the computer-assisted in combo screening approach for the discovery of the anti-cholinesterase compounds from C. intybus.

Acetylcholinesterase↗

Atypical muscarinic allosteric modulation: cooperativity between modulators and their atypical binding topology in muscarinic M2 and M2/M5 chimeric receptors.

The binding and function of muscarinic acetylcholine receptors can be modulated allosterically. Some allosteric muscarinic ligands are "atypical", having steep concentration-effect curves and not interacting competitively with "typical" allosteric modulators. For atypical agents, a second allosteric site has been proposed. Different approaches have been used to gain further insight into the interaction with M2 receptors of two atypical agents, tacrine and the bispyridinium compound 4,4'-bis-[(2,6-dichloro-benzyloxy-imino)-methyl]-1,1'-propane-1,3-diyl-bispyridinium dibromide (Duo3). Interaction studies, using radioligand binding assays and the allosteric ligands obidoxime, Mg2+, and the new tool hexamethonium to antagonize the allosteric actions of the atypical ligands, showed different modes of interaction for tacrine and Duo3 at M2 receptors. A negatively cooperative interaction was observed between hexamethonium and tacrine (but not Duo3). A tacrine dimer that exhibited increased allosteric potency relative to tacrine but behaved like a typical allosteric modulator was competitively inhibited by hexamethonium. M2/M5-receptor mutants revealed a dependence of tacrine and Duo3 affinity on different receptor epitopes. This was confirmed by docking simulations using a three-dimensional model of the M2 receptor. These showed that the allosteric site could accommodate two molecules of tacrine simultaneously but only one molecule of Duo3, which binds in different mode from typical allosteric agents. Therefore, the atypical actions of tacrine and Duo3 involve different modes of receptor interaction, but their sites of attachment seem to be the "common" allosteric binding domain at the entrance to the orthosteric ligand binding pocket of the M2-receptor. Additional complex behavior may be rationalized by allosteric interactions transmitted within a receptor dimer.

Allosteric Regulation↗

A structural model of EmrE, a multi-drug transporter from Escherichia coli.

Using a recently reported computational method, we describe an approach to model the structure of EmrE, a proton coupled multi-drug transporter of Escherichia coli. EmrE is the smallest ion-coupled transporter known; it functions as an oligomer and each monomer comprises four transmembrane segments. Because of its size, EmrE provides a unique experimental paradigm. The computational method does not afford a unique solution for the monomer. The experimental constraints available were used to select the most likely structure and to dock two monomers together to yield a dimer. The model is further validated by modeling of Hsmr, an EmrE homolog with a remarkable amino acid composition with over 40% of Ala and Val. The Hsmr model is similar to that of EmrE, with the majority of the Ala or Val residues facing the lipid. In addition, the model of EmrE features a putative substrate-binding site very similar to that observed in BmrR, a transcription activator of multi-drug transporters, with a similar substrate profile. The two crucial residues that couple proton fluxes with substrate binding in the homo-dimer of EmrE, Glu-14, have a spatial arrangement that agrees with proposed molecular mechanisms of transport.

Alanine↗

Refined solution structure of omega-conotoxin GVIA: implications for calcium channel binding.

The polypeptide omega-conotoxin GVIA (GVIA) is an N-type calcium channel blocker from the venom of Conus geographus, a fish-hunting cone shell. Here we describe a high-resolution solution structure of this member of the 'inhibitor cystine knot' protein family. The structure, based on NMR data acquired at 600 MHz, has mean pairwise RMS differences of 0.25 +/- 0.06 and 1.07 +/- 0.14 A over the backbone heavy atoms and all heavy atoms, respectively. The solvent-accessible side chains are better defined than in previously published structures and provide an improved basis for docking GVIA with models of the calcium channel. Moreover, some side chain interactions important in GVIA folding in vitro and in stabilizing the native structure are defined clearly in the refined structure. Two qualitatively different backbone conformations in the segment from Thr11 to Asn14 persisted in the restrained simulated annealing calculations until a small number of lower bound constraints was included to prevent close contacts from occurring that did not correspond with peaks in the NOESY spectrum. It is possible that GVIA is genuinely flexible at this segment, spending a finite time in the alternative conformation, and this may influence its interaction with the calcium channel.

Amino Acid Sequence↗

Modelling of the interaction of verotoxin-1 (VT1) with its glycolipid receptor, globotriaosylceramide (Gb3).

Possible binding sites for the glycolipid globotriaosylceramide (Gal alpha 1-->4Gal beta 1-->4Glc beta 1-->1 Cer; Gb3) on the B-subunits of verotoxin-1 (VT1) were explored using binding data for specifically mutated verotoxins and by computational docking of favoured conformers of Gb3 with the crystal structure of VT1. Calculations using the GRID program suggested a site with favourable hydrophobic interactions at the exposed side chain of Phe30. One of the favoured conformers of Gb3 was docked into this site, with the hydrophobic face of the internal Gal beta residue in contact with the side chain of Phe30. After energy minimization, the two terminal saccharide residues of Gb3 (Gal alpha and Gal beta) showed favourable interactions with the toxin. In the proposed model of the complex, the terminal Gal alpha of Gb3 is located in proximity to aspartates 16-18 of VT1. The model is in agreement with available experimental binding data for the interaction of globoglycolipids with different naturally occurring and mutated verotoxins.

Animals↗

Heteroligation of a mouse monoclonal IgE antibody (La2) with small molecules, analysed by computer-aided automated docking.

A mouse monoclonal anti-TNP IgE antibody (IgE-La2) was screened by a competitive-binding ELISA with a random pool of over 2000 small molecules, mostly drugs, drug derivatives and metabolites. Thirteen of these (naproxene, beta-carboxy-chi-naphthol, oxolinic acid, hymecromone, 8-aminoquinoline, beta-naphthylamine, chi-nitrilo-cinnamic acid, 1,5-diaminonaphthaline, prolonium iodide, diaspirin, 3,4,5-trimethoxy-cinnamic acid, cycrimine, hemimellitic acid) were found to bind as strongly, or stronger, to the antibody as the immunizing hapten. We have used a Monte Carlo search technique for simulated docking of the DNP and non-DNP ligands to a model of the Fv region of IgE(La2). The validity of structural predictions made by the AutoDock program were tested on IgG(ANO2), the three-dimensional structure of which had been obtained previously by X-ray crystallography and 2D-NMR. The rms differences between the experimentally determined and auto-docked complexes in the energetically most favored binding modes were 0.31-0.44 A. Evaluation of structures of IgE(La2)-ligand complexes [including 2,4-dinitrophenol (DNP), 16 DNP amino acids, and the 13 non-DNP ligands listed above] obtained by computer-aided automated docking, suggested the existence of two subsites within an approximately 12 x 18 A2 groove extending between the H and L CDRs. Some of the ligands (DNP-Glu, 8-aminoquinoline, prolonium-I, beta-naphthylamine) were found to bind exclusively to subsite 1, others (DNP-Ala, chi-nitrilo-cinnamic acid, hemimellitic acid, beta-carboxy-chi-naphthol) to subsite 2. The majority of DNP amino acids and other ligands (oxolinic acid, 3,4,5-trimethoxy-cinnamic acid, diaspirin, [R]-cycrimine) were found to occupy an overlapping area including subsites 1 and 2, while some of the compounds (DNP-Asn, DNP-Pro, hymecromone, 1,5-naphthylenediamine) were predicted to interact with either of these subsites with comparable probabilities. When all of the docked La2-ligand complexes were taken into account, five tyrosine residues (H33, L32, L91, L92, L96) were found to provide the majority (53.4%) of all observed contact points. Thus, a multitude of interactions with aromatic residues, and a combinatorial type of interaction within the binding region, seem to be the major factors to explain the mechanism of heteroligation by IgE(La2).

Amino Acid Sequence↗

Artificial transfer hydrogenases based on the biotin-(strept)avidin technology: fine tuning the selectivity by saturation mutagenesis of the host protein.

Incorporation of biotinylated racemic three-legged d6-piano stool complexes in streptavidin yields enantioselective transfer hydrogenation artificial metalloenzymes for the reduction of ketones. Having identified the most promising organometallic catalyst precursors in the presence of wild-type streptavidin, fine-tuning of the selectivity is achieved by saturation mutagenesis at position S112. This choice for the genetic optimization site is suggested by docking studies which reveal that this position lies closest to the biotinylated metal upon incorporation into streptavidin. For aromatic ketones, the reaction proceeds smoothly to afford the corresponding enantioenriched alcohols in up to 97% ee (R) or 70% (S). On the basis of these results, we suggest that the enantioselection is mostly dictated by CH/pi interactions between the substrate and the eta6-bound arene. However, these enantiodiscriminating interactions can be outweighed in the presence of cationic residues at position S112 to afford the opposite enantiomers of the product.

Biotin↗

Site-saturation mutagenesis and three-dimensional modelling of ROB-1 define a substrate binding role of Ser130 in class A beta-lactamases.

Site-saturation mutagenesis was performed on the class A ROB-1 beta-lactamase at conserved Ser130, which is centrally located in the antibiotic binding site where it can participate in both protein-protein and protein-substrate hydrogen bonding. Mutation Thr130 gave a beta-lactamase hydrolysing penicillins and cephalosporins but which showed a 3-fold lower affinity (Km) for ampicillin and cephalexin, and a 30-fold lower hydrolytic (Vmax) activity for ampicillin. In contrast, the hydrolytic activity for cephalexin was similar to the wild-type for the Thr130 mutation. Mutation Gly130 gave a beta-lactamase hydrolysing only penicillins with an affinity and hydrolysis activity for these compounds approximately 15-fold lower than the wild-type, but no detectable activity against cephalosporins. Mutation Ala130 produced an enzyme capable of hydrolysing penicillins only at a low rate. Modelling the ROB-1 active site was done from the refined 2 A X-ray structure of the homologous Bacillus licheniformis beta-lactamase. Ampicillin and cephalexin were docked into the active site and were energy minimized with the CVFF empirical force field. Dockings were stable only when Ser70 was made anionic and Glu166 was made neutral. Interaction energies and distances were calculated for fully hydrated pre-acylation complexes with the Ser, Thr, Gly and Ala130 enzymes. The catalytic data from all mutations and the computed interactions from modelling confirmed that the Ser130 has a structural as well as a functional role in binding and hydrolysis of penicillins. This highly conserved residue also plays a substrate specificity role by hydrogen binding the carboxylic acid group of cephalosporins more tightly than penicillins.

Ampicillin↗

Receptor-based 3D QSAR analysis of estrogen receptor ligands--merging the accuracy of receptor-based alignments with the computational efficiency of ligand-based methods.

One of the major challenges in computational approaches to drug design is the accurate prediction of binding affinity of biomolecules. In the present study several prediction methods for a published set of estrogen receptor ligands are investigated and compared. The binding modes of 30 ligands were determined using the docking program AutoDock and were compared with available X-ray structures of estrogen receptor-ligand complexes. On the basis of the docking results an interaction energy-based model, which uses the information of the whole ligand-receptor complex, was generated. Several parameters were modified in order to analyze their influence onto the correlation between binding affinities and calculated ligand-receptor interaction energies. The highest correlation coefficient (r2 = 0.617, q2Loo = 0.570) was obtained considering protein flexibility during the interaction energy evaluation. The second prediction method uses a combination of receptor-based and 3D quantitative structure-activity relationships (3D QSAR) methods. The ligand alignment obtained from the docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection, a significant and robust model was obtained (r2 = 0.991, q2LOO = 0.921). The predictive ability of the established model was further evaluated by using a test set of six additional compounds. The comparison with the generated interaction energy-based model and with a traditional CoMFA model obtained using a ligand-based alignment (r2 = 0.951, q2L00 = 0.796) indicates that the combination of receptor-based and 3D QSAR methods is able to improve the quality of the underlying model.

Binding Sites↗

Modelling study of protein kinase inhibitors: binding mode of staurosporine and origin of the selectivity of CGP 52411.

A model for the binding mode of the potent protein kinase inhibitor staurosporine is proposed. Using the information provided by the crystal structure of the cyclic-AMP-dependent protein kinase, it is suggested that staurosporine, despite a seemingly unrelated chemical structure, exploits the same key hydrogen-bond interactions as ATP, the cofactor of the protein kinases, in its binding mode. The structure-activity relationship of the inhibitor and a docking analysis give strong support to this hypothesis. The selectivity of the dianilinophthalimide inhibitor CGP 52411 towards the EGF-receptor protein tyrosine kinase is rationalized on the basis of the model. It is proposed that this selectivity originates in the occupancy, by one of the anilino moieties of the inhibitor, of the region of the enzyme cleft that normally binds the ribose ring of ATP, which appears to possess a marked lipophilic character in this kinase.

Alkaloids↗

In silico structure-based design of a potent and selective small peptide inhibitor of protein tyrosine phosphatase 1B, a novel therapeutic target for obesity and type 2 diabetes mellitus: a computer modeling approach.

Protein Tyrosine Phosphatase 1B (PTP1B) has been shown to be a negative regulator of insulin signaling by dephosphorylating key tyrosine residues within the regulatory domain of the beta-subunit of the insulin receptor. Recent gene knockout studies in mice have shown the mice to have increased insulin sensitivity and improved glucose tolerance. Furthermore, these mice also exhibited a resistance to diet induced obesity. Inhibitors of PTP1B would have the potential of enhancing insulin action by prolonging the phosphorylated state of the insulin receptor. In addition, recent clinical studies have shown that the haplotype ACTTCAG0 of the PTPN1 gene, which encodes PTP1B, is a major risk contributor to type 2 diabetes mellitus (T2DM). Thus, there is compelling evidence that small molecule inhibitors of PTP1B may be effective in treating insulin resistance at an early stage, thereby leading to a prevention strategy for T2DM and obesity. Based on the crystal structure of the complex of PTP1B with a known inhibitor, we have identified a tetrapeptide inhibitor with the sequence WKPD. Docking calculations indicate that this peptide is as potent as the existing inhibitors. Moreover, the peptide is also found to be selective for PTP1B with a greatly reduced potency against other biologically important protein tyrosine phosphatases such as PTP-LAR, Calcineurin, and the highly homologous T-Cell Protein Tyrosine Phosphatase (TCPTP). Thus the designed tetrapeptide is a suitable lead compound for the development of new drugs against type 2 diabetes and obesity.

Amino Acid Sequence↗

Identification of an extracytoplasmic region of H+,K(+)-ATPase labeled by a K(+)-competitive photoaffinity inhibitor.

The photoaffinity reagent 8-[(4-azidophenyl)-methoxy]-1-tritiomethyl-2, 3-dimethylimidazo-[1,2-alpha]pyridinium iodide ([3H]mDAZIP) has been synthesized and used to photoinactivate and label purified hog gastric H+,K(+)-ATPase. The specific (K(+)-sensitive) components of both photoinactivation and labeling showed dependences on inhibitor concentration consistent with covalent modification at an extracytoplasmic site of reversible K(+)-competitive binding in the dark. The maximum amount of specific labeling (1.2 nmol/mg) was similar to the number of phosphorylation sites measured (1.0 +/- 0.14 nmol/mg). Specific labeling was distributed 76% on the alpha chain, 18% on the beta chain, and 6% on undefined peptides. Various digestions with trypsin, protease V8, and thermolysin were employed to fragment the labeled enzyme. Gasphase sequencing of the radioactive peptides identified the major site of specific labeling to be within a region where only two stretches of amino acids (Leu105 to Ile126 and Leu139 to Phe155, designated H1 and H2, respectively) are predicted to span the membrane. This in turn suggested that the labeling site was located within or close to the proposed loop between them (Gln127 to Asn138). A computer-driven energy minimization protocol yielded a loop structure to which SCH 28080 (the parent structure of [3H]mDAZIP) could be docked. Conversely, modeling of the corresponding region of Na+,K(+)-ATPase (a homologous enzyme with much lower affinity for SCH 28080) yielded no apparent binding site. Similarities in the inhibition of H+,K(+)-ATPase by SCH 28080 and of Na+,K(+)-ATPase by ouabain lead to the hypothesis that, in each case, inhibitor binding to E2-P is associated with an increase in the hydrophobicity of the environment of the loop between H1 and H2.

Adenosine Triphosphatases↗