X-ray analysis of azido-thymidine diphosphate binding to nucleoside diphosphate kinase

To be effective as antiviral agent, AZT (3′-azido-3′-deoxythymidine) must be converted to a triphosphate derivative by cellular kinases. The conversion is inefficient and, to understand why AZT diphosphate is a poor substrate of nucleoside diphosphate (NDP) kinase, we determined a 2.3-Å x-ray structure of a complex with the N119A point mutant of Dictyostelium NDP kinase. It shows that the analog binds at the same site and, except for the sugar ring pucker which is different, binds in the same way as the natural substrate thymidine diphosphate. However, the azido group that replaces the 3′OH of the deoxyribose in AZT displaces a lysine side chain involved in catalysis. Moreover, it is unable to make an internal hydrogen bond to the oxygen bridging the β- and γ-phosphate, which plays an important part in phosphate transfer.

A zinc finger directory for high-affinity DNA recognition

We have used two monovalent phage display libraries containing variants of the Zif268 DNA-binding domain to obtain families of zinc fingers that bind to alterations in the last 4 bp of the DNA sequence of the Zif268 consensus operator, GCG TGGGCG. Affinity selection was performed by altering the Zif268 operator three base pairs at a time, and simultaneously selecting for sets of 16 related DNA sequences. In this way, only four experiments were required to select for all possible 64 combinations of DNA triplet sequences. The results show that (i) for high-affinity DNA binding in the range observed for the Zif268 wild-type complex (Kd = 0.5–5 nM), finger 1 specifically requires the arginine at the carboxy terminus of its recognition helix that forms a bidentate hydrogen-bond with the guanine base (G) in the crystal structure of Zif268 complexed to its DNA operator sequence GCG TGG GCG; (ii) when the guanine base (G) is replaced by A, C, or T, a lower-affinity family (Kd ⩾ 50 nM) can be detected that shows an overall tendency to bind G-rich DNA; (iii) the residues at position 2 on the finger 2 recognition helix do not appear to interact strongly with the complementary 5′ base in the finger 1 binding site; and (iv) unexpected substitutions at the amino terminus of finger 1 can occasionally result in specificity for the 3′ base in the finger 1 binding site. A DNA recognition directory was constructed for high-affinity zinc fingers that recognize all three bases in a DNA triplet for seven sequences of the type GNN. Similar approaches may be applied to other zinc fingers to broaden the scope of the directory.

Specific mutations in a viral RNA pseudoknot drastically change ribosomal frameshifting efficiency

Many viruses regulate protein synthesis by −1 ribosomal frameshifting using an RNA pseudoknot. Frameshifting is vital for viral reproduction. Using the information gained from the recent high-resolution crystal structure of the beet western yellow virus pseudoknot, a systematic mutational analysis has been carried out in vitro and in vivo. We find that specific nucleotide tertiary interactions at the junction between the two stems of the pseudoknot are crucial. A triplex is found between stem 1 and loop 2, and triplex interactions are required for frameshifting function. For some mutations, loss of one hydrogen bond is sufficient to abolish frameshifting. Furthermore, mutations near the 5′ end of the pseudoknot can increase frameshifting by nearly 300%, possibly by modifying ribosomal contacts. It is likely that the selection of suitable mutations can thus allow viruses to adjust frameshifting efficiencies and thereby regulate protein synthesis in response to environmental change.

A physical basis for protein secondary structure

A physical theory of protein secondary structure is proposed and tested by performing exceedingly simple Monte Carlo simulations. In essence, secondary structure propensities are predominantly a consequence of two competing local effects, one favoring hydrogen bond formation in helices and turns, the other opposing the attendant reduction in sidechain conformational entropy on helix and turn formation. These sequence specific biases are densely dispersed throughout the unfolded polypeptide chain, where they serve to preorganize the folding process and largely, but imperfectly, anticipate the native secondary structure.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 × 106 M−1, two orders of magnitude higher than that for biotin (1 × 104 M−1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 × 104 M−1) was lower than predicted (2.9 × 105 M−1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

(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.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Tight Asp-85–Thr-89 association during the pump switch of bacteriorhodopsin

Unidirectional proton transport in bacteriorhodopsin is enforced by the switching machinery of the active site. Threonine 89 is located in this region, with its O—H group forming a hydrogen bond with Asp-85, the acceptor for proton transfer from the Schiff base of the retinal chromophore. Previous IR spectroscopy of [3-18O]threonine-labeled bacteriorhodopsin showed that the hydrogen bond of the O—D group of Thr-89 in D2O is strengthened in the K photocycle intermediate. Here, we show that the strength and orientation of this hydrogen bond remains unchanged in the L intermediate and through the M intermediate. Furthermore, a strong interaction between Asp-85 and the O—H (O—D) group of Thr-89 in M is indicated by a shift in the C⩵O stretching vibration of the former because of 18O substitution in the latter. Thus, the strong hydrogen bond between Asp-85 and Thr-89 in K persists through M, contrary to structural models based on x-ray crystallography of the photocycle intermediates. We propose that, upon photoisomerization of the chromophore, Thr-89 forms a tight, persistent complex with one of the side-chain oxygens of Asp-85 and is thereby precluded from participating in the switching process. On the other hand, the loss of hydrogen bonding at the other oxygen of Asp-85 in M may be related to the switching event.

Structural analysis of the binding modes of minor groove ligands comprised of disubstituted benzenes

Two-dimensional homonuclear NMR was used to characterize synthetic DNA minor groove-binding ligands in complexes with oligonucleotides containing three different A-T binding sites. The three ligands studied have a C2 axis of symmetry and have the same general structural motif of a central para-substituted benzene ring flanked by two meta-substituted rings, giving the molecules a crescent shape. As with other ligands of this shape, specificity seems to arise from a tight fit in the narrow minor groove of the preferred A-T-rich sequences. We found that these ligands slide between binding subsites, behavior attributed to the fact that all of the amide protons in the ligand backbone cannot hydrogen bond to the minor groove simultaneously.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Polar residues drive association of polyleucine transmembrane helices

Although many polar residues are directly involved in transmembrane protein functions, the extent to which they contribute to more general structural features is still unclear. Previous studies have demonstrated that asparagine residues can drive transmembrane helix association through interhelical hydrogen bonding [Choma, C., Gratkowski, H., Lear, J. D. & DeGrado, W. F. (2000) Nat. Struct. Biol. 7, 161–166; and Zhou, F. X., Cocco, M. J., Russ, W. P., Brunger, A. T. & Engelman, D. M. (2000) Nat. Struct. Biol. 7, 154–160]. We have studied the ability of other polar residues to promote helix association in detergent micelles and in biological membranes. Our results show that polyleucine sequences with Asn, Asp, Gln, Glu, and His, residues capable of being simultaneously hydrogen bond donors and acceptors, form homo- or heterooligomers. In contrast, polyleucine sequences with Ser, Thr, and Tyr do not associate more than the polyleucine sequence alone. The results therefore provide experimental evidence that interactions between polar residues in the helices of transmembrane proteins may serve to provide structural stability and oligomerization specificity. Furthermore, such interactions can allow structural flexibility required for the function of some membrane proteins.

An amyloid-forming peptide from the yeast prion Sup35 reveals a dehydrated β-sheet structure for amyloid

X-ray diffraction and other biophysical tools reveal features of the atomic structure of an amyloid-like crystal. Sup35, a prion-like protein in yeast, forms fibrillar amyloid assemblies intrinsic to its prion function. We have identified a polar peptide from the N-terminal prion-determining domain of Sup35 that exhibits the amyloid properties of full-length Sup35, including cooperative kinetics of aggregation, fibril formation, binding of the dye Congo red, and the characteristic cross-β x-ray diffraction pattern. Microcrystals of this peptide also share the principal properties of the fibrillar amyloid, including a highly stable, β-sheet-rich structure and the binding of Congo red. The x-ray powder pattern of the microcrystals, extending to 0.9-Å resolution, yields the unit cell dimensions of the well-ordered structure. These dimensions restrict possible atomic models of this amyloid-like structure and demonstrate that it forms packed, parallel-stranded β-sheets. The unusually high density of the crystals shows that the packed β-sheets are dehydrated, despite the polar character of the side chains. These results suggest that amyloid is a highly intermolecularly bonded, dehydrated array of densely packed β-sheets. This dry β-sheet could form as Sup35 partially unfolds to expose the peptide, permitting it to hydrogen-bond to the same peptide of other Sup35 molecules. The implication is that amyloid-forming units may be short segments of proteins, exposed for interactions by partial unfolding.

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

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

Toward a quantum-mechanical description of metal-assisted phosphoryl transfer in pyrophosphatase

The wealth of kinetic and structural information makes inorganic pyrophosphatases (PPases) a good model system to study the details of enzymatic phosphoryl transfer. The enzyme accelerates metal-complexed phosphoryl transfer 1010-fold: but how? Our structures of the yeast PPase product complex at 1.15 Å and fluoride-inhibited complex at 1.9 Å visualize the active site in three different states: substrate-bound, immediate product bound, and relaxed product bound. These span the steps around chemical catalysis and provide strong evidence that a water molecule (Onu) directly attacks PPi with a pKa vastly lowered by coordination to two metal ions and D117. They also suggest that a low-barrier hydrogen bond (LBHB) forms between D117 and Onu, in part because of steric crowding by W100 and N116. Direct visualization of the double bonds on the phosphates appears possible. The flexible side chains at the top of the active site absorb the motion involved in the reaction, which may help accelerate catalysis. Relaxation of the product allows a new nucleophile to be generated and creates symmetry in the elementary catalytic steps on the enzyme. We are thus moving closer to understanding phosphoryl transfer in PPases at the quantum mechanical level. Ultra-high resolution structures can thus tease out overlapping complexes and so are as relevant to discussion of enzyme mechanism as structures produced by time-resolved crystallography.