Determining binding sites in protein-nucleic acid complexes by cross-saturation.

Cross-saturation experiments have been shown to give accurate information regarding the interacting surfaces in protein-protein and protein-RNA complexes. The rate of magnetization transfer depends on a number of factors including geometry, spin topology, and effective correlation times. To assess the influence of these variables on such experiments, and to determine the range of applicability of the technique, we have simulated the time-course of magnetization transfer across the interface in a variety of protein-nucleic acid complexes (434 Cro, SRY, MetJ and U1A), each having a different binding geometry. The simulations have been carried out primarily to provide information about the experimentally accessible targets for selective saturation, such as the anomeric protons and the imino protons of the nucleic acid. Saturation of either of these groups of signals leads to partial excitation throughout the nucleic acid molecule, and the resulting transfer of saturation to the labelled protein moiety can be readily detected by the reduction in intensity of particular peaks in the HSQC spectrum of the protein. The simulations show that information can be obtained about the residues in contact with the nucleic acid without any need for deuteration. Experimental cross-saturation data have been obtained from the Mbp1-DNA complex and interpreted in conjunction with the simulations to map out the binding surface in detail.

Solution structure, hydrodynamics and thermodynamics of the UvrB C-terminal domain.

The solution structure, thermodynamic stability and hydrodynamic properties of the 55-residue C-terminal domain of UvrB that interacts with UvrC during excision repair in E. coli have been determined using a combination of high resolution NMR, ultracentrifugation, 15N NMR relaxation, gel permeation, NMR diffusion, circular dichroism and differential scanning calorimetry. The subunit molecular weight is 7,438 kDa., compared with 14.5+/-1.0 kDa. determined by equilibrium sedimentation, indicating a dimeric structure. The structure determined from NMR showed a stable dimer of anti-parallel helical hairpins that associate in an unusual manner, with a small and hydrophobic interface. The Stokes radius of the protein decreases from a high plateau value (ca. 22 A) at protein concentrations greater than 4 microM to about 18 A at concentrations less than 0.1 microM. The concentration and temperature-dependence of the far UV circular dichroism show that the protein is thermally stable (Tm ca. 71.5 degrees C at 36 microM). The simplest model consistent with these data was a dimer dissociating into folded monomers that then unfolds co-operatively. The van't Hoff enthalpy and dissociation constant for both transition was derived by fitting, with deltaH1=23 kJ mol(-1). K1(298)=0.4 microM and deltaH2= 184 kJ mol(-1). This is in good agreement with direct calorimetric analysis of the thermal unfolding of the protein, which gave a calorimetric enthalpy change of 181 kJ mol(-1) and a van't Hoff enthalpy change of 354 kJ mol(-1), confirming the dimer to monomer unfolding. The thermodynamic data can be reconciled with the observed mode of dimerisation. 15N NMR relaxation measurements at 14.1 T and 11.75 T confirmed that the protein behaves as an asymmetric dimer at mM concentrations, with a flexible N-terminal linker for attachment to the remainder of the UvrB protein. The role of dimerisation of this domain in the excision repair mechanism is discussed.

Ca2+-independent binding of an EF-hand domain to a novel motif in the alpha-actinin-titin complex.

The interaction between alpha-actinin and titin, two modular muscle proteins, is essential for sarcomere assembly. We have solved the solution structure of a complex between the calcium-insensitive C-terminal EF-hand domain of alpha-actinin-2 and the seventh Z-repeat of titin. The structure of the complex is in a semi-open conformation and closely resembles that of myosin light chains in their complexes with heavy chain IQ motifs. However, no IQ motif is present in the Z-repeat, suggesting that the semi-open conformation is a general structural solution for calcium-independent recognition of EF-hand domains.

Towards a structural understanding of Friedreich's ataxia: the solution structure of frataxin.

BACKGROUND: Lesions in the gene for frataxin, a nuclear-encoded mitochondrial protein, cause the recessively inherited condition Friedreich's ataxia. It is thought that the condition arises from disregulation of mitochondrial iron homeostasis, with concomitant oxidative damage leading to neuronal death. Very little is, as yet, known about the biochemical function of frataxin. RESULTS: Here, we show that the mature form of recombinant frataxin behaves in solution as a monodisperse species that is composed of a 15-residue-long unstructured N terminus and an evolutionarily conserved C-terminal region that is able to fold independently. The structure of the C-terminal domain consists of a stable seven-stranded antiparallel beta sheet packing against a pair of parallel helices. The structure is compact with neither grooves nor cavities, features that are typical of iron-binding modules. Exposed evolutionarily conserved residues cover a broad area and all cluster on the beta-sheet face of the structure, suggesting that this is a functionally important surface. The effect of two clinically occurring mutations on the fold was checked experimentally. When the mature protein was titrated with iron, no tendency to iron-binding or to aggregation was observed. CONCLUSIONS: Knowledge of the frataxin structure provides important guidelines as to the nature of the frataxin binding partner. The absence of all the features expected for an iron-binding activity, the large conserved area on its surface and lack of evidence for iron-binding activity strongly support an indirect involvement of frataxin in iron metabolism. The effects of point mutations associated with Friedreich's ataxia can be rationalised by knowledge of the structure and suggest possible models for the occurrence of the disease in compound heterozygous patients.

The influence of DNA binding on the backbone dynamics of the yeast cell-cycle protein Mbp1.

Mbp1 is a transcription factor involved in the regulation of the cell cycle in yeast. The N-terminus of this protein contains a DNA binding domain that includes a winged helix-turn-helix motif. The C-terminal 24 residues of this domain (the 'tail') are disordered in the crystal state, but are important for DNA binding. We have measured 15N NMR relaxation rates at 11.75 and 14.1 T to determine the dynamics of the free protein and in its complex with a specific DNA duplex. The dynamics data were quantitatively analysed using both spectral density mapping and the Lipari-Szabo formalism including the effects of chemical exchange and rotational anisotropy. A detailed analysis has been made of the effect of anisotropy, exchange and experimental precision on the recovered motional parameters. The backbone NH relaxation is affected by motions on a variety of time scales from millisecond to tens of picoseconds. The relaxation data show a structured core of 100 residues corresponding to that observed in the crystal state. Within the core of the protein, two regions on either side of the putative recognition helix (helix B) show slow (ca. 0.2 ms) conformational exchange dynamics that are quenched upon DNA binding. The C-terminal 24 residues are generally more dynamic than in the core. However, in the free protein, a stretch of approximately 8 residues in the middle of the tail show relaxation behaviour similar to that in the core, indicating a structured region. NOEs between Ala 114 in this structured part of the tail and residues in the N-terminal beta strand of the core of the protein demonstrate that the tail folds back onto the core of the protein. In the complex with DNA, the structured part of the tail extends by ca. 3 residues. These data provide a framework for understanding the biochemical data on the mechanism and specificity of DNA binding.

Identification of platination sites on human serum transferrin using (13)C and (15)N NMR spectroscopy.

Reactions between various apo and metal-bound forms of human serum transferrin (80 kDa) and the recombinant N-lobe (40 kDa) with [Pt(en)Cl(2)] or cis-[PtCl(2)(NH(3))(2)] have been investigated in solution via observation of [(1)H,(15)N] NMR resonances of the Pt complexes, [(1)H,(13)C] resonances of the eCH(3) groups of the protein methionine residues, and by chromatographic analysis of single-site methionine mutants. For the whole protein, the preferred Pt binding site appears to be Met256. Additional binding occurs at the other surface-exposed methionine (Met499), which is platinated at a slower rate than Met256. In contrast, binding of similar Pt compounds to the N-lobe of the protein occurs at Met313, rather than Met256. Met313 is buried in the interlobe contact region of intact transferrin. After loss of one chloride ligand from Pt and binding to methionine sulfur of the N-lobe, chelate-ring closure appears to occur with binding to a deprotonated backbone amide nitrogen, and the loss of the other chloride ligand. Such chelate-ring closure was not observed during reactions of the whole protein, even after several days.

Solution structure of an EGF module pair from the Plasmodium falciparum merozoite surface protein 1.

The solution structure of the 96-residue C-terminal fragment of the merozoite surface protein 1 (MSP-1) from Plasmodium falciparum has been determined using nuclear magnetic resonance (NMR) spectroscopic measurements on uniformly13C/15N-labelled protein, efficiently expressed in the methylotrophic yeast Komagataella (Pichia) pastoris. The structure has two domains with epidermal growth factor (EGF)-like folds with a novel domain interface for the EGF domain pair interactions, formed from a cluster of hydrophobic residues. This gives the protein a U-shaped overall structure with the N-terminal proteolytic processing site close to the C-terminal glycosyl phosphatidyl inositol (GPI) membrane anchor site, which is consistent with the involvement of a membrane-bound proteinase in the processing of MSP-1 during erythrocyte invasion. This structure, which is the first protozoan EGF example to be determined, contrasts with the elongated structures seen for EGF-module pairs having shared Ca2+-ligation sites at their interface, as found, for example, in fibrillin-1. Recognition surfaces for antibodies that inhibit processing and invasion, and antibodies that block the binding of these inhibitory antibodies, have been mapped on the three-dimensional structure by considering specific MSP-1 mutants.

Structure and dynamics in solution of the complex of Lactobacillus casei dihydrofolate reductase with the new lipophilic antifolate drug trimetrexate.

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase and the anticancer drug trimetrexate. Two thousand seventy distance, 345 dihedral angle, and 144 hydrogen bond restraints were obtained from analysis of multidimensional NMR spectra recorded for complexes containing 15N-labeled protein. Simulated annealing calculations produced a family of 22 structures fully consistent with the constraints. Several intermolecular protein-ligand NOEs were obtained by using a novel approach monitoring temperature effects of NOE signals resulting from dynamic processes in the bound ligand. At low temperature (5 degrees C) the trimethoxy ring of bound trimetrexate is flipping sufficiently slowly to give narrow signals in slow exchange, which give good NOE cross peaks. At higher temperature these broaden and their NOE cross peaks disappear thus allowing the signals in the lower-temperature spectrum to be identified as NOEs involving ligand protons. The binding site for trimetrexate is well defined and this was compared with the binding sites in related complexes formed with methotrexate and trimethoprim. No major conformational differences were detected between the different complexes. The 2,4-diaminopyrimidine-containing moieties in the three drugs bind essentially in the same binding pocket and the remaining parts of their molecules adapt their conformations such that they can make effective van der Waals interactions with essentially the same set of hydrophobic amino acids, the side-chain orientations and local conformations of which are not greatly changed in the different complexes (similar chi1 and chi2 values).

High resolution structure of the N-terminal domain of tissue inhibitor of metalloproteinases-2 and characterization of its interaction site with matrix metalloproteinase-3.

The high resolution structure of the N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP-2) in solution has been determined using multidimensional heteronuclear NMR spectroscopy, with the structural calculations based on an extensive set of constraints, including 3132 nuclear Overhauser effect-based distance constraints, 56 hydrogen bond constraints, and 220 torsion angle constraints (an average of 26.9 constraints/residue). The core of the protein consists of a five-stranded beta-barrel that is homologous to the beta-barrel found in the oligosaccharide/oligonucleotide binding protein fold. The binding site for the catalytic domain of matrix metalloproteinases-3 (N-MMP-3) on N-TIMP-2 has been mapped by determining the changes in chemical shifts on complex formation for signals from the protein backbone (15N, 13C, and 1H). This approach identified a discrete N-MMP-3 binding site on N-TIMP-2 composed of the N terminus of the protein and the loops between beta-strands AB, CD, and EF. The beta-hairpin formed from strands A and B in N-TIMP-2 is significantly longer than the equivalent structure in TIMP-1, allowing it to make more extensive binding interactions with the MMP catalytic domain. A detailed comparison of the N-TIMP-2 structure with that of TIMP-1 bound to N-MMP-3 (Gomis-Ruth, F.-X., Maskos, K., Betz, M., Bergner, A., Huber, R., Suzuki, K., Yoshida, N., Nagase, H. , Brew, K., Bourne, G. P., Bartunik, H. & Bode, W. (1997) Nature 389, 77-80) revealed that the core beta-barrels are very similar in topology but that the loop connecting beta-strands CD (P67-C72) would need to undergo a large conformational change for TIMP-2 to bind in a similar manner to TIMP-1.

Solution structure of the B-Myb DNA-binding domain: a possible role for conformational instability of the protein in DNA binding and control of gene expression.

Double- and triple-resonance heteronuclear NMR spectroscopy have been used to determine the high-resolution solution structure of the minimal B-Myb DNA-binding domain (B-MybR2R3) and to characterize the specific complex formed with a synthetic DNA fragment corresponding to the Myb target site on the Myb-regulated gene tom-1. B-MybR2R3 is shown to consist of two independent protein domains (R2 and R3) joined by a short linker, which have strikingly different tertiary structures despite significant sequence similarities. In addition, the C-terminal region of B-Myb R2 is confirmed to have a poorly defined structure, reflecting the existence of multiple conformations in slow to intermediate exchange. This contrasts with the tertiary structure reported for c-MybR2R3, in which both R2 and R3 have the same fold and the C-terminal region of R2 forms a stable, well-defined helix [Ogata, K., et al. (1995) Nat. Struct. Biol. 2, 309-320]. The NMR data suggest there are extensive contacts between B-MybR2R3 and its DNA target site in the complex and are consistent with a significant conformational change in the protein on binding to DNA, with one possibility being the formation of a stable helix in the C-terminal region of R2. In addition, conformational heterogeneity identified in R2 of B-MybR2R3 bound to the tom-1-A target site may play an important role in the control of gene expression by Myb proteins.

The solution structure of the complex of Lactobacillus casei dihydrofolate reductase with methotrexate.

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase (18.3 kDa, 162 amino acid residues) formed with the anticancer drug methotrexate using 2531 distance, 361 dihedral angle and 48 hydrogen bond restraints obtained from analysis of multidimensional NMR spectra. Simulated annealing calculations produced a family of 21 structures fully consistent with the constraints. The structure has four alpha-helices and eight beta-strands with two other regions, comprising residues 11 to 14 and 126 to 127, also interacting with each other in a beta-sheet manner. The methotrexate binding site is very well defined and the structure around its glutamate moiety was improved by including restraints reflecting the previously determined specific interactions between the glutamate alpha-carboxylate group with Arg57 and the gamma-carboxylate group with His28. The overall fold of the binary complex in solution is very similar to that observed in the X-ray studies of the ternary complex of L. casei dihydrofolate reductase formed with methotrexate and NADPH (the structures of the binary and ternary complexes have a root-mean-square difference over the backbone atoms of 0.97 A). Thus no major conformational change takes place when NADPH binds to the binary complex. In the binary complex, the loop comprising residues 9 to 23 which forms part of the active site has been shown to be in the "closed" conformation as defined by M. R. Sawaya & J. Kraut, who considered the corresponding loops in crystal structures of complexes of dihydrofolate reductases from several organisms. Thus the absence of the NADPH does not result in the "occluded" form of the loop as seen in crystal studies of some other dihydrofolate reductases in the absence of coenzyme. Some regions of the structure in the binary complex which form interaction sites for NADPH are less well defined than other regions. However, in general terms, the NADPH binding site appears to be essentially pre-formed in the binary complex. This may contribute to the tighter binding of coenzyme in the presence of methotrexate.

Mapping the binding site for matrix metalloproteinase on the N-terminal domain of the tissue inhibitor of metalloproteinases-2 by NMR chemical shift perturbation.

Changes in the NMR chemical shift of backbone amide nuclei (1H and 15N) have been used to map the matrix metalloproteinase (MMP) binding site on the N-terminal domain of the tissue inhibitor of metalloproteinase-2 (N-TIMP-2). Amide chemical shift changes were measured on formation of a stable complex with the catalytic domain of stromelysin-1 (N-MMP-3). Residues with significantly shifted amide signals mapped specifically to a broad site covering one face of the molecule. This site (the MMP binding site) consists primarily of residues 1-11, 27-41, 68-73, 87-90, and 97-104. The site overlaps with the OB-fold binding site seen in other proteins that share the same five-stranded beta-barrel topology. Sequence conservation data and recent site-directed mutagenesis studies are discussed in relation to the MMP binding site identified in this work.

Trp128Tyr mutation in the N-lobe of recombinant human serum transferrin: 1H- and 15N-NMR and metal binding studies.

The conserved Trp residue within helix 5 of the N-lobe of human serum transferrin (hTF/2N, 40 kDa) has been mutated to Tyr. NMR and CD spectra and energy calculations show that the mutation causes little perturbation of the overall structure of hTF/2N although the chelating agent Tiron removed Fe3+ from the mutant protein about three times faster than from wild-type hTF/2N. 1H-NMR resonances of residues in the Leu122-Trp128-Ile132 hydrophobic patch are assigned both by ring current calculations and with the aid of the mutation. [1H, 15N]-NMR resonances for 11 of the 14 Tyr residues were observed in the spectra of 15N-Tyr-hTF/2N and a resonance for Tyr128 was assignable in spectra of the mutant. The 15N resonance of Y128 was sensitive to oxalate and Ga3+ binding, and Ga3+ binding perturbed 15N resonances for most of the Tyr residues. Since these are well distributed over the N-lobe, it can be concluded that metal-induced structural changes are not merely local to the binding site.

Correlated bond rotations in interactions of arginine residues with ligand carboxylate groups in protein ligand complexes.

The 1H/15N HSQC NMR spectra of complexes of Lactobacillus casei dihydrofolate reductase containing methotrexate recorded at 1 degree C show four resolved signals for the four NH(eta) protons of the Arg57 residue. This is consistent with hindered rotation in the guanidino group resulting from interactions with the alpha-carboxylate of methotrexate. Increasing the temperature causes exchange line-broadening and coalescence of signals. Rotation rates for the N(epsilon)C(zeta) and C(zeta)N(eta) bonds have been calculated from lineshape analysis and from zz-HSQC exchange experiments. The interactions between the methotrexate alpha-carboxylate group and the Arg57 guanidino group decrease the rotation rates for the N(epsilon)C(zeta) bond by about a factor of 10 and those for the C(zeta)N(eta) bonds by more than a factor of 100 with respect to their values in free arginine. Furthermore, the relative rates of rotation about these two bonds are reversed in the protein complexes compared with their values in free arginine indicating that there are concerted rotations about the N(epsilon)C(zeta) bond of the Arg57 guanidino group and the C'C(alpha) bond of the glutamate alpha-carboxylate group of methotrexate.

High-resolution solution structure of human pNR-2/pS2: a single trefoil motif protein.

pNR-2/pS2 is a 60 residue extracellular protein, which was originally discovered in human breast cancer cells, and subsequently found in other tumours and normal gastric epithelial cells. We have determined the three-dimensional solution structure of a C58S mutant of human pNR-2/pS2 using 639 distance and 137 torsion angle constraints obtained from analysis of multidimensional NMR spectra. A series of simulated annealing calculations resulted in the unambiguous determination of the protein's disulphide bonding pattern and produced a family of 19 structures consistent with the constraints. The peptide contains a single "trefoil" sequence motif, a region of about 40 residues with a characteristic sequence pattern, which has been found, either singly or as a repeat, in about a dozen extracellular proteins. The trefoil domain contains three disulphide bonds, whose 1-5, 2-4 and 3-6 cysteine pairings form the structure into three closely packed loops with only a small amount of secondary structure, which consists of a short alpha-helix packed against a two-stranded antiparallel beta-sheet. The structure of the domain is very similar to those of the two trefoil domains that occur in porcine spasmolytic polypeptide (PSP), the only member of the trefoil family whose three-dimensional structure has been previously determined. Outside the trefoil domain, which forms the compact "head" of the molecule, the N and C-terminal strands are closely associated, forming an extended "tail", which has some beta-sheet character for part of its length and which becomes more disordered towards the termini as indicated by (15)N{(1)H} NOEs. We have considered the structural implications of the possible formation of a native C58-C58 disulphide-bonded homodimer. Comparison of the surface features of pNR-2/pS2 and PSP, and consideration of the sequences of the other human trefoil domains in the light of these structures, illuminates the possible role of specific residues in ligand/receptor binding.

Hydration and solution structure of d(CGCAAATTTGCG)2 and its complex with propamidine from NMR and molecular modelling.

The hydration of the d(CGCAAATTTGCG)2 duplex and its complex with a propamidine reporter ligand has been examined in aqueous solution by two-dimensional NMR at two spectrometer frequencies and three temperatures. Quantitative analysis of ROESY and NOESY cross-peaks showed effective correlation times of approximately 0.5 ns at 283 K for DNA-water interactions in the major groove. In some cases the sign of the NOE inverts on changing either the temperature or spectrometer frequency. Larger effective correlation times of approximately 1 ns were observed for water interactions with A5(H2) and A6(H2) atoms located in the minor groove. Interproton NOEs and changes in chemical shifts showed that propamidine binds in the minor groove 5'-AATTT region of the host duplex, but does not displace waters adjacent to either A5(H2) or A6(H2). In the complex, the effective correlation times of these waters increase more than two-fold, possibly as a result of stabilisation due to H-bonded interaction with the amidine groups of the ligand. Hydration of the bound molecule was also found, suggesting that water may contribute to the DNA binding process for bis(amidine) drugs. Structure refinement by a NOE-restrained dynamic annealing procedure revealed that ligand binding is non-centrosymmetric with respect to the duplex, in accordance with the energetically favoured 5'-ATT (= 5'-AAT) sites predicted by analytical molecular modelling. In particular, the bound propamidine spans 3-4 base pairs in the A6-T7-T8 tract and makes close H-bonded contacts with A(N3/O4) acceptors positioned close to the minor groove floor. The refined NMR structure for the DNA-propamidine complex is compared with that determined recently using X-ray crystallographic methods.