A novel assay for drug-DNA binding mode, affinity, and exclusion number: scanning force microscopy.

Determining the mode-of-binding of a DNA ligand is not always straightforward. Here, we establish a scanning force microscopic assay for mode-of-binding that is (i) direct: lengths of individual DNA-ligand complexes are directly measured; (ii) rapid: there are no requirements for staining or elaborate sample preparation; and (iii) unambiguous: an observed increase in DNA length upon addition of a ligand is definitive evidence for an intercalative mode-of-binding. Mode-of-binding, binding affinity, and site-exclusion number are readily determined from scanning force microscopy measurements of the changes in length of individual drug-DNA complexes as a function of drug concentration. With this assay, we resolve the ambiguity surrounding the mode of binding of 2,5-bis(4-amidinophenyl) furan (APF) to DNA and show that it binds to DNA by nonintercalative modes. APF is a member of an important class of aromatic dicationic drugs that show significant activity in the treatment of Pneumocystis carinii pneumonia, an opportunistic infection that is the leading cause of death in AIDS patients.

Crystal structure of D-amino acid oxidase: a case of active site mirror-image convergent evolution with flavocytochrome b2.

D-amino acid oxidase is the prototype of the FAD-dependent oxidases. It catalyses the oxidation of D-amino acids to the corresponding alpha-ketoacids. The reducing equivalents are transferred to molecular oxygen with production of hydrogen peroxide. We have solved the crystal structure of the complex of D-amino acid oxidase with benzoate, a competitive inhibitor of the substrate, by single isomorphous replacement and eightfold averaging. Each monomer is formed by two domains with an overall topology similar to that of p-hydroxybenzoate hydroxylase. The benzoate molecule lays parallel to the flavin ring and is held in position by a salt bridge with Arg-283. Analysis of the active site shows that no side chains are properly positioned to act as the postulated base required for the catalytic carboanion mechanism. On the contrary, the benzoate binding mode suggests a direct transfer of the substrate alpha-hydrogen to the flavin during the enzyme reductive half-reaction.The active site Of D-amino acid oxidase exhibits a striking similarity with that of flavocytochrome b2, a structurally unrelated FMN-dependent flavoenzyme. The active site groups (if these two enzymes are in fact superimposable once the mirror-image of the flavocytochrome b2 active site is generated with respect to the flavin plane. Therefore, the catalytic sites of D-amino acid oxidase and flavocytochrome b2 appear to have converged to a highly similar but enantiomeric architecture in order to catalvze similar reactions (oxidation of alpha-amino acids or alpha-hydroxy acids), although with opposite stereochemistry.

Recognition principle of the TAP transporter disclosed by combinatorial peptide libraries

Transport of peptides across the membrane of the endoplasmic reticulum for assembly with MHC class I molecules is an essential step in antigen presentation to cytotoxic T cells. This task is performed by the major histocompatibility complex-encoded transporter associated with antigen processing (TAP). Using a combinatorial approach we have analyzed the substrate specificity of human TAP at high resolution and in the absence of any given sequence context, revealing the contribution of each peptide residue in stabilizing binding to TAP. Human TAP was found to be highly selective with peptide affinities covering at least three orders of magnitude. Interestingly, the selectivity is not equally distributed over the substrate. Only the N-terminal three positions and the C-terminal residue are critical, whereas effects from other peptide positions are negligible. A major influence from the peptide backbone was uncovered by peptide scans and libraries containing d amino acids. Again, independent of peptide length, critical positions were clustered near the peptide termini. These approaches demonstrate that human TAP is selective, with residues determining the affinity located in distinct regions, and point to the role of the peptide backbone in binding to TAP. This binding mode of TAP has implications in an optimized repertoire selection and in a coevolution with the major histocompatibility complex/T cell receptor complex.

The Aer protein and the serine chemoreceptor Tsr independently sense intracellular energy levels and transduce oxygen, redox, and energy signals for Escherichia coli behavior

We identified a protein, Aer, as a signal transducer that senses intracellular energy levels rather than the external environment and that transduces signals for aerotaxis (taxis to oxygen) and other energy-dependent behavioral responses in Escherichia coli. Domains in Aer are similar to the signaling domain in chemotaxis receptors and the putative oxygen-sensing domain of some transcriptional activators. A putative FAD-binding site in the N-terminal domain of Aer shares a consensus sequence with the NifL, Bat, and Wc-1 signal-transducing proteins that regulate gene expression in response to redox changes, oxygen, and blue light, respectively. A double mutant deficient in aer and tsr, which codes for the serine chemoreceptor, was negative for aerotaxis, redox taxis, and glycerol taxis, each of which requires the proton motive force and/or electron transport system for signaling. We propose that Aer and Tsr sense the proton motive force or cellular redox state and thereby integrate diverse signals that guide E. coli to environments where maximal energy is available for growth.

(4-Aminomethyl)phenylguanidine derivatives as nonpeptidic highly selective inhibitors of human urokinase

Increased expression of the serine protease urokinase-type plasminogen activator (uPA) in tumor tissues is highly correlated with tumor cell migration, invasion, proliferation, progression, and metastasis. Thus inhibition of uPA activity represents a promising target for antimetastatic therapy. So far, only the x-ray crystal structure of uPA inactivated by H-Glu-Gly-Arg-chloromethylketone has been reported, thus limited data are available for a rational structure-based design of uPA inhibitors. Taking into account the trypsin-like arginine specificity of uPA, (4-aminomethyl)phenylguanidine was selected as a potential P1 residue and iterative derivatization of its amino group with various hydrophobic residues, and structure–activity relationship-based optimization of the spacer in terms of hydrogen bond acceptor/donor properties led to N-(1-adamantyl)-N′-(4-guanidinobenzyl)urea as a highly selective nonpeptidic uPA inhibitor. The x-ray crystal structure of the uPA B-chain complexed with this inhibitor revealed a surprising binding mode consisting of the expected insertion of the phenylguanidine moiety into the S1 pocket, but with the adamantyl residue protruding toward the hydrophobic S1′ enzyme subsite, thus exposing the ureido group to hydrogen-bonding interactions. Although in this enzyme-bound state the inhibitor is crossing the active site, interactions with the catalytic residues Ser-195 and His-57 are not observed, but their side chains are spatially displaced for steric reasons. Compared with other trypsin-like serine proteases, the S2 and S3/S4 pockets of uPA are reduced in size because of the 99-insertion loop. Therefore, the peculiar binding mode of the new type of uPA inhibitors offers the possibility of exploiting optimized interactions at the S1′/S2′ subsites to further enhance selectivity and potency. Because crystals of the uPA/benzamidine complex allow inhibitor exchange by soaking procedures, the structure-based design of new generations of uPA inhibitors can rely on the assistance of x-ray analysis.

Three-dimensional structure of human electron transfer flavoprotein to 2.1-Å resolution

Mammalian electron transfer flavoproteins (ETF) are heterodimers containing a single equivalent of flavin adenine dinucleotide (FAD). They function as electron shuttles between primary flavoprotein dehydrogenases involved in mitochondrial fatty acid and amino acid catabolism and the membrane-bound electron transfer flavoprotein ubiquinone oxidoreductase. The structure of human ETF solved to 2.1-Å resolution reveals that the ETF molecule is comprised of three distinct domains: two domains are contributed by the α subunit and the third domain is made up entirely by the β subunit. The N-terminal portion of the α subunit and the majority of the β subunit have identical polypeptide folds, in the absence of any sequence homology. FAD lies in a cleft between the two subunits, with most of the FAD molecule residing in the C-terminal portion of the α subunit. Alignment of all the known sequences for the ETF α subunits together with the putative FixB gene product shows that the residues directly involved in FAD binding are conserved. A hydrogen bond is formed between the N5 of the FAD isoalloxazine ring and the hydroxyl side chain of αT266, suggesting why the pathogenic mutation, αT266M, affects ETF activity in patients with glutaric acidemia type II. Hydrogen bonds between the 4′-hydroxyl of the ribityl chain of FAD and N1 of the isoalloxazine ring, and between αH286 and the C2-carbonyl oxygen of the isoalloxazine ring, may play a role in the stabilization of the anionic semiquinone. With the known structure of medium chain acyl-CoA dehydrogenase, we hypothesize a possible structure for docking the two proteins.

A β-hairpin structure in a 13-mer peptide that binds α-bungarotoxin with high affinity and neutralizes its toxicity

Snake-venom α-bungarotoxin is a member of the α-neurotoxin family that binds with very high affinity to the nicotinic acetylcholine receptor (AChR) at the neuromuscular junction. The structure of the complex between α-bungarotoxin and a 13-mer peptide (WRYYESSLEPYPD) that binds the toxin with high affinity, thus inhibiting its interactions with AChR with an IC50 of 2 nM, has been solved by 1H-NMR spectroscopy. The bound peptide folds into a β-hairpin structure created by two antiparallel β-strands, which combine with the already existing triple-stranded β-sheet of the toxin to form a five-stranded intermolecular, antiparallel β-sheet. Peptide residues Y3P, E5P, and L8P have the highest intermolecular contact area, indicating their importance in the binding of α-bungarotoxin; W1P, R2P, and Y4P also contribute significantly to the binding. A large number of characteristic hydrogen bonds and electrostatic and hydrophobic interactions are observed in the complex. The high-affinity peptide exhibits inhibitory potency that is better than any known peptide derived from AChR, and is equal to that of the whole α-subunit of AChR. The high degree of sequence similarity between the peptide and various types of AChRs implies that the binding mode found within the complex might possibly mimic the receptor binding to the toxin. The design of the high-affinity peptide was based on our previous findings: (i) the detection of a lead peptide (MRYYESSLKSYPD) that binds α-bungarotoxin, using a phage-display peptide library, (ii) the information about the three-dimensional structure of α-bungarotoxin/lead-peptide complex, and (iii) the amino acid sequence analysis of different AChRs.

Unraveling principles of lead discovery: from unfrustrated energy landscapes to novel molecular anchors.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as "recognition nuclei" is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Selenocysteine, identified as the penultimate C-terminal residue in human T-cell thioredoxin reductase, corresponds to TGA in the human placental gene.

The possible relationship of selenium to immunological function which has been suggested for decades was investigated in studies on selenium metabolism in human T cells. One of the major 75Se-labeled selenoproteins detected was purified to homogeneity and shown to be a homodimer of 55-kDa subunits. Each subunit contained about 1 FAD and at least 0.74 Se. This protein proved to be thioredoxin reductase (TR) on the basis of its catalytic activities, cross-reactivity with anti-rat liver TR antibodies, and sequence identities of several tryptic peptides with the published deduced sequence of human placental TR. Physicochemical characteristics of T-cell TR were similar to those of a selenocysteine (Secys)-containing TR recently isolated from human lung adenocarcinoma cells. The sequence of a 12-residue 75Se-labeled tryptic peptide from T-cell TR was identical with a C-terminal-deduced sequence of human placental TR except that Secys was present in the position corresponding to TGA, previously thought to be the termination codon, and this was followed by Gly-499, the actual C-terminal amino acid. The presence of the unusual conserved Cys-Secys-Gly sequence at the C terminus of TR in addition to the redox active cysteines of the Cys-Val-Asn-Val-Gly-Cys motif in the FAD-binding region may account for the peroxidase activity and the relatively low substrate specificity of mammalian TRs. The finding that T-cell TR is a selenoenzyme that contains Se in a conserved C-terminal region provides another example of the role of selenium in a major antioxidant enzyme system (i.e., thioredoxin-thioredoxin reductase), in addition to the well-known glutathione peroxidase enzyme system.

Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.

Multiple loop structures critical for ligand binding of the integrin α4 subunit in the upper face of the β-propeller mode 1

A non-I-domain integrin, α4β1, recognizes vascular cell adhesion molecule 1 (VCAM-1) and the IIICS portion of fibronectin. To localize regions of α4 critical for ligand binding, we swapped several predicted loops within or near the putative ligand-binding site of α4 (which spans repeats 2–5 of the seven N-terminal repeats) with the corresponding regions of α5. Swapping residues 112–131 in repeat 2, or residues 237–247 in repeat 4, completely blocked adhesion to immobilized VCAM-1 and connecting segment 1 (CS-1) peptide. However, swapping residues 40–52 in repeat 1, residues 151–164 in repeat 3, or residues 282–288 (which contain a putative cation binding motif) in repeat 5 did not affect or only slightly reduced adhesion to these ligands. The binding of several function-blocking antibodies is blocked by swapping residues 112–131, 151–164, and 186–191 (which contain previously identified residues critical for ligand binding, Tyr-187 and Gly-190). These results are consistent with the recently published β-propeller folding model of the integrin α4 subunit [Springer, T. A. (1997) Proc. Natl. Acad. Sci. USA 94, 65–72], in which seven four-stranded β-sheets are arranged in a torus around a pseudosymmetric axis. The regions of α4 critical for ligand binding are adjacent to each other and are located in the upper face, the predicted ligand-binding site, of the β-propeller model, although they are not adjacent in the primary structure.

The mode of action and the structure of a herbicide in complex with its target: binding of activated hydantocidin to the feedback regulation site of adenylosuccinate synthetase.

(+)-Hydantocidin, a recently discovered natural spironucleoside with potent herbicidal activity, is shown to be a proherbicide that, after phosphorylation at the 5' position, inhibits adenylosuccinate synthetase, an enzyme involved in de novo purine synthesis. The mode of binding of hydantocidin 5'-monophosphate to the target enzyme was analyzed by determining the crystal structure of the enzyme-inhibitor complex at 2.6-A resolution. It was found that adenylosuccinate synthetase binds the phosphorylated compound in the same fashion as it does adenosine 5'-monophosphate, the natural feedback regulator of this enzyme. This work provides the first crystal structure of a herbicide-target complex reported to date.

Batimastat, a potent matrix mealloproteinase inhibitor, exhibits an unexpected mode of binding.

Matrix metalloproteinase enzymes have been implicated in degenerative processes like tumor cell invasion, metastasis, and arthritis. Specific metalloproteinase inhibitors have been used to block tumor cell proliferation. We have examined the interaction of batimastat (BB-94) with a metalloproteinase [atrolysin C (Ht-d), EC 3.4.24.42] active site at 2.0-angstroms resolution (R = 16.8%). The title structure exhibits an unexpected binding geometry, with the thiophene ring deeply inserted into the primary specificity site. This unprecedented binding geometry dramatizes the significance of the cavernous primary specificity site, pointing the way for the design of a new generation of potential antitumor drugs.