Determination of helix orientations in a flexible DNA by multi-frequency EPR spectroscopy.

Distance measurements are performed between a pair of spin labels attached to nucleic acids using Pulsed Electron-Electron Double Resonance (PELDOR, also called DEER) spectroscopy which is a complementary tool to other structure determination methods in structural biology. The rigid spin label Ç, when incorporated pairwise into two helical parts of a nucleic acid molecule, allows the determination of both the mutual orientation and the distance between those labels, since Ç moves rigidly with the helix to which it is attached. We have developed a two-step protocol to investigate the conformational flexibility of flexible nucleic acid molecules by multi-frequency PELDOR. In the first step, a library with a broad collection of conformers, which are in agreement with topological constraints, NMR restraints and distances derived from PELDOR, was created. In the second step, a weighted structural ensemble of these conformers was chosen, such that it fits the multi-frequency PELDOR time traces of all doubly Ç-labelled samples simultaneously. This ensemble reflects the global structure and the conformational flexibility of the two-way DNA junction. We demonstrate this approach on a flexible bent DNA molecule, consisting of two short helical parts with a five adenine bulge at the center. The kink and twist motions between both helical parts were quantitatively determined and showed high flexibility, in agreement with a Förster Resonance Energy Transfer (FRET) study on a similar bent DNA motif. The approach presented here should be useful to describe the relative orientation of helical motifs and the conformational flexibility of nucleic acid structures, both alone and in complexes with proteins and other molecules.

Comparative NMR analysis of an 80-residue G protein-coupled receptor fragment in two membrane mimetic environments.

Fragments of integral membrane proteins have been used to study the physical chemical properties of regions of transporters and receptors. Ste2p(G31-T110) is an 80-residue polypeptide which contains a portion of the N-terminal domain, transmembrane domain 1 (TM1), intracellular loop 1, TM2 and part of extracellular loop 1 of the α-factor receptor (Ste2p) from Saccharomyces cerevisiae. The structure of this peptide was previously determined to form a helical hairpin in lyso-palmitoylphosphatidyl-glycerol micelles (LPPG) [1]. Herein, we perform a systematic comparison of the structure of this protein fragment in micelles and trifluoroethanol (TFE):water in order to understand whether spectra recorded in organic:aqueous medium can facilitate the structure determination in a micellar environment. Using uniformly labeled peptide and peptide selectively protonated on Ile, Val and Leu methyl groups in a perdeuterated background and a broad set of 3D NMR experiments we assigned 89% of the observable atoms. NOEs and chemical shift analysis were used to define the helical regions of the fragment. Together with constraints from paramagnetic spin labeling, NOEs were used to calculate a transiently folded helical hairpin structure for this peptide in TFE:water. Correlation of chemical shifts was insufficient to transfer assignments from TFE:water to LPPG spectra in the absence of further information.

Torsion angle dynamics for NMR structure calculation with the new program DYANA.

The new program DYANA (DYnamics Algorithm for Nmr Applications) for efficient calculation of three-dimensional protein and nucleic acid structures from distance constraints and torsion angle constraints collected by nuclear magnetic resonance (NMR) experiments performs simulated annealing by molecular dynamics in torsion angle space and uses a fast recursive algorithm to integrate the equations of motions. Torsion angle dynamics can be more efficient than molecular dynamics in Cartesian coordinate space because of the reduced number of degrees of freedom and the concomitant absence of high-frequency bond and angle vibrations, which allows for the use of longer time-steps and/or higher temperatures in the structure calculation. It also represents a significant advance over the variable target function method in torsion angle space with the REDAC strategy used by the predecessor program DIANA. DYANA computation times per accepted conformer in the "bundle" used to represent the NMR structure compare favorably with those of other presently available structure calculation algorithms, and are of the order of 160 seconds for a protein of 165 amino acid residues when using a DEC Alpha 8400 5/300 computer. Test calculations starting from conformers with random torsion angle values further showed that DYANA is capable of efficient calculation of high-quality protein structures with up to 400 amino acid residues, and of nucleic acid structures.

Efficient computation of three-dimensional protein structures in solution from nuclear magnetic resonance data using the program DIANA and the supporting programs CALIBA, HABAS and GLOMSA.

A novel procedure for efficient computation of three-dimensional protein structures from nuclear magnetic resonance (n.m.r.) data in solution is described, which is based on using the program DIANA in combination with the supporting programs CALIBA, HABAS and GLOMSA. The first part of this paper describes the new programs DIANA. CALIBA and GLOMSA. DIANA is a new, fully vectorized implementation of the variable target function algorithm for the computation of protein structures from n.m.r. data. Its main advantages, when compared to previously available programs using the variable target function algorithm, are a significant reduction of the computation time, and a novel treatment of experimental distance constraints involving diastereotopic groups of hydrogen atoms that were not individually assigned. CALIBA converts the measured nuclear Overhauser effects into upper distance limits and thus prepares the input for the previously described program HABAS and for DIANA. GLOMSA is used for obtaining individual assignments for pairs of diastereotopic substituents by comparison of the experimental constraints with preliminary results of the structure calculations. With its general outlay, the presently used combination of the four programs is particularly user-friendly. In the second part of the paper, initial results are presented on the influence of the novel DIANA treatment of diastereotopic protons on the quality of the structures obtained, and a systematic study of the central processing unit times needed for the same protein structure calculation on a range of different, commonly available computers is described.

The NMR solution conformation of unligated human cyclophilin A.

The nuclear magnetic resonance (NMR) solution structure of free, unligated cyclophilin A (CypA), which is an 18 kDa protein from human T-lymphocytes that was expressed in Escherichia coli for the present study, was determined using multidimensional heteronuclear NMR techniques. Sequence-specific resonance assignments for 99.5% of all backbone amide protons and non-labile hydrogen atoms provided the basis for collection of an input of 4101 nuclear Overhauser enhancement (NOE) upper distance constraints and 371 dihedral angle constraints for distance geometry calculations and energy minimization with the programs DIANA and OPAL. The average RMSD values of the 20 best energy-refined NMR conformers relative to the mean coordinates are 0.49 A for the backbone atoms and 0.88 A for all heavy atoms of residues 2 to 165. The molecular architecture includes an eight-stranded antiparallel beta-barrel that is closed by two amphipathic alpha-helices. Detailed comparisons with the crystal structure of free CypA revealed subtle but significant conformational differences that can in most cases be related to lattice contacts in the crystal structure. 15N spin relaxation times and NMR lineshape analyses for CypA in the free form and complexed with cyclosporin A (CsA) revealed transitions of polypeptide loops surrounding the ligand-binding site from locally flexible conformations in the free protein, some of which include well-defined conformational equilibria, to well-defined spatial arrangements in the CypA-CsA complex. Compared to the crystal structure of free CypA, where the ligand-binding area is extensively involved in lattice contacts, the NMR structure presents a highly relevant reference for studies of changes in structure and internal mobility of the binding pocket upon ligand binding, and possible consequences of this conformational variability for calcineurin recognition by the CypA-CsA complex.

Nuclear magnetic resonance solution structure of dendrotoxin K from the venom of Dendroaspis polylepis polylepis.

The solution structure of dendrotoxin K (Toxin K), a protein consisting of one polypeptide chain with 57 residues purified from the venom of the black mamba, Dendroaspis polylepis polylepis, was determined by nuclear magnetic resonance (NMR) spectroscopy. On the basis of virtually complete sequence-specific 1H NMR assignments, including individual assignments for 38 pairs of diastereotopic substituents and side-chain amide protons, a total of 818 nuclear Overhauser effect distance constraints and 123 dihedral angle constraints were identified. Using this input, the solution structure of Toxin K was calculated with the program DIANA, and refined by restrained energy-minimization with a modified version of the program AMBER. The average root-mean-square deviation (r.m.s.d.) relative to the mean atomic co-ordinates of the 20 conformers selected to represent the solution structure is 0.31 A for all backbone atoms N, C alpha and C', and 0.90 A for all heavy-atoms of residues 2 to 56. The solution structure of Toxin K is very similar to the solution structure of the basic pancreatic trypsin inhibitor (BPTI) and the X-ray crystal structure of the alpha-dendrotoxin from Dendroaspis angusticeps (alpha-DTX), with r.m.s.d. values of 1.31 A and 0.92 A, respectively, for the backbone atoms of residues 2 to 56. Some local structural differences between Toxin K and BPTI are directly related to the fact that intermolecular interactions with two of the four internal molecules of hydration water in BPTI are replaced by intramolecular hydrogen bonds in Toxin K.

Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.

Nuclear magnetic resonance solution structure of hirudin(1-51) and comparison with corresponding three-dimensional structures determined using the complete 65-residue hirudin polypeptide chain.

The three-dimensional structure of the N-terminal 51-residue domain of recombinant hirudin in aqueous solution was determined by 1H nuclear magnetic resonance (NMR) spectroscopy, and the resulting high-quality solution structure was compared with corresponding structures obtained from studies with the intact, 65-residue polypeptide chain of hirudin. On the basis of 580 distance constraints derived from nuclear Overhauser effects and 109 dihedral angle constraints, a group of 20 conformers representing the solution structure of hirudin(1-51) was computed with the program DIANA and energy-minimized with a modified version of the program AMBER. Residues 3 to 30 and 37 to 48 form a well-defined molecular core with two antiparallel beta-sheets composed of residues 14 to 16 and 20 to 22, and 27 to 31 and 36 to 40, and three reverse turns at residues 8 to 11 (type II), 17 to 20 (type II') and 23 to 26 (type II). The average root-mean-square deviation of the individual NMR conformers relative to their mean co-ordinates is 0.38 A for the backbone atoms and 0.77 A for all heavy atoms of these residues. Increased structural disorder was found for the N-terminal dipeptide segment, the loop at residues 31 to 36, and the C-terminal tripeptide segment. The solution structure of hirudin(1-51) has the same molecular architecture as the corresponding polypeptide segment in natural hirudin and recombinant desulfatohirudin. It is also closely similar to the crystal structure of the N-terminal 51-residue segment of hirudin in a hirudin-thrombin complex, with root-mean-square deviations of the crystal structure relative to the mean solution structure of 0.61 A for the backbone atoms and 0.91 A for all heavy atoms of residues 3 to 30 and 37 to 48. Further coincidence is found for the loop formed by residues 31 to 36, which shows increased structural disorder in all available solution structures of hirudin, and of which residues 32 to 35 are not observable in the electron density map of the thrombin complex. Significant local structural differences between hirudin(1-51) in solution and hirudin in the crystalline thrombin complex were identified mainly for the N-terminal tripeptide segment and residues 17 to 21. These are further analyzed in an accompanying paper.

Structure determination of the Antp (C39----S) homeodomain from nuclear magnetic resonance data in solution using a novel strategy for the structure calculation with the programs DIANA, CALIBA, HABAS and GLOMSA.

The structure of a mutant Antennapedia homeodomain, Antp(C39----S), from Drosophila melanogaster was determined using a set of new programs introduced in the accompanying paper. An input dataset of 957 distance constraints and 171 dihedral angle constraints was collected using two-dimensional n.m.r. experiments with the 15N-labeled protein. The resulting high quality structure for Antp(C39----S), with an average root-mean-square deviation of 0.53 A between the backbone atoms of residues 7 to 59 in 20 energy-refined distance geometry structures and the mean structure, is nearly identical to the previously reported structure of the wild-type Antp homeodomain. The only significant difference is in the connection between helices III and IV, which was found to be less kinked than was indicated by the structure determination for Antp. The main emphasis of the presentation in this paper is on a detailed account of the practical use of a novel strategy for the computation of nuclear magnetic resonance structures of proteins with the combined use of the programs DIANA, CALIBA, HABAS and GLOMSA.

Impact of protein-protein contacts on the conformation of thrombin-bound hirudin studied by comparison with the nuclear magnetic resonance solution structure of hirudin(1-51).

The impact of protein-protein interactions on the conformation of the N-terminal hirudin domain consisting of residues 1 to 51 in the X-ray crystal structure of a hirudin-thrombin complex was investigated through comparisons with the nuclear magnetic resonance solution structure of hirudin(1-51). The close overall similarity observed between these two structures contrasts with the behavior of the C-terminal 17-residue polypeptide segment of hirudin, which is flexibly disordered in solution but exhibits a defined conformation in the complex with thrombin. Localized structural differences in the N-terminal domain include that residues 1 to 3 of hirudin in the crystalline complex form a hydrogen-bonding network with thrombin that is reminiscent of a parallel beta-sheet. Moreover, the backbone conformation of residues 17 to 20 in the complex does not contain the characteristic hydrogen bond observed for the type II' reverse turn in the solution structure, and the side-chains of Ser19 and Val21 have significantly different orientations in the two structures. Most of these structural changes can be related directly to thrombin-hirudin contacts, which may also be an important factor in the mechanism of hirudin action. In this context, it is of special interest that other residues that also make numerous contacts with thrombin, e.g. Thr4, Asp5 and Asn20, have identical conformations in free hirudin and in the complex.

NMR structure of the sea urchin (Strongylocentrotus purpuratus) metallothionein MTA.

The three-dimensional structure of [(113)Cd7]-metallothionein-A (MTA) of the sea urchin Strongylocentrotus purpuratus was determined by homonuclear(1)H NMR experiments and heteronuclear [(1)H, (113)Cd]-correlation spectroscopy. MTA is composed of two globular domains, an N-terminal four-metal domain of the amino acid residues 1 to 36 and a Cd4Cys11cluster, and a C-terminal three-metal domain including the amino acid residues 37 to 65 and a Cd3Cys9cluster. The structure resembles the known mammalian and crustacean metallothioneins, but has a significantly different connectivity pattern of the Cys-metal co-ordination bonds and concomitantly contains novel local folds of some polypeptide backbone segments. These differences can be related to variations of the Cys sequence positions and thus emphasize the special role of the cysteine residues in defining the structure of metallothioneins, both on the level of the domain architecture and the topology of the metal-thiolate clusters.

NMR solution structure of the recombinant tick anticoagulant protein (rTAP), a factor Xa inhibitor from the tick Ornithodoros moubata.

The solution structure of the recombinant tick anticoagulant protein (rTAP) was determined by 1H nuclear magnetic resonance (NMR) spectroscopy in aqueous solution at pH 3.6 and 36 degrees C. rTAP is a 60-residue protein functioning as a highly specific inhibitor of the coagulation protease factor Xa, which was originally isolated from the tick Ornithodoros moubata. Its regular secondary structure consists of a two-stranded antiparallel beta-sheet with residues 22-28 and 32-38, and an alpha-helix with residues 51-60. The relative orientation of these regular secondary structure elements has nearly identical counterparts in the bovine pancreatic trypsin inhibitor (BPTI). In contrast, the loop between the beta-sheet and the C-terminal alpha-helix as well as the N-terminal 20-residue segment preceding the beta-sheet adopt different three-dimensional folds in the two proteins. These observations are discussed with regard to the implication of different mechanisms of protease inhibition by rTAP and by Kunitz-type protein proteinase inhibitors.

Automated peak picking and peak integration in macromolecular NMR spectra using AUTOPSY.

A new approach for automated peak picking of multidimensional protein NMR spectra with strong overlap is introduced, which makes use of the program AUTOPSY (automated peak picking for NMR spectroscopy). The main elements of this program are a novel function for local noise level calculation, the use of symmetry considerations, and the use of lineshapes extracted from well-separated peaks for resolving groups of strongly overlapping peaks. The algorithm generates peak lists with precise chemical shift and integral intensities, and a reliability measure for the recognition of each peak. The results of automated peak picking of NOESY spectra with AUTOPSY were tested in combination with the combined automated NOESY cross peak assignment and structure calculation routine NOAH implemented in the program DYANA. The quality of the resulting structures was found to be comparable with those from corresponding data obtained with manual peak picking.

Conformational stability and activity of p73 require a second helix in the tetramerization domain.

p73 and p63, the two ancestral members of the p53 family, are involved in neurogenesis, epithelial stem cell maintenance and quality control of female germ cells. The highly conserved oligomerization domain (OD) of tumor suppressor p53 is essential for its biological functions, and its structure was believed to be the prototype for all three proteins. However, we report that the ODs of p73 and p63 differ from the OD of p53 by containing an additional alpha-helix that is not present in the structure of the p53 OD. Deletion of this helix causes a dissociation of the OD into dimers; it also causes conformational instability and reduces the transcriptional activity of p73. Moreover, we show that ODs of p73 and p63 strongly interact and that a large number of different heterotetramers are supported by the additional helix. Detailed analysis shows that the heterotetramer consisting of two homodimers is thermodynamically more stable than the two homotetramers. No heterooligomerization between p53 and the p73/p63 subfamily was observed, supporting the notion of functional orthogonality within the p53 family.

Improved efficiency of protein structure calculations from NMR data using the program DIANA with redundant dihedral angle constraints.

A new strategy for NMR structure calculations of proteins with the variable target function method (Braun, W. and Go, N. (1985) J. Mol. Biol., 186, 611) is described, which makes use of redundant dihedral angle constraints (REDAC) derived from preliminary calculations of the complete structure. The REDAC approach reduces the computation time for obtaining a group of acceptable conformers with the program DIANA 5-100-fold, depending on the complexity of the protein structure, and retains good sampling of conformation space.

Automated combined assignment of NOESY spectra and three-dimensional protein structure determination.

A procedure for automated protein structure determination is presented that is based on an iterative procedure during which the NOESY peak list assignment and the structure calculation are performed simultaneously. The input consists of a list of NOESY peak positions and a list of chemical shifts as obtained from sequence-specific resonance assignment. For the present applications of this approach the previously introduced NOAH routine was implemented in the distance geometry program DIANA. As an illustration, experimental 2D and 3D NOESY cross-peak lists of six proteins have been analyzed, for which complete sequence-specific 1H assignments are available for the polypeptide backbone and the amino acid side chains. The automated method assigned 70-90% of all NOESY cross peaks, which is on average 10% less than with the interactive approach, and only between 0.8% and 2.4% of the automatically assigned peaks had a different assignment than in the corresponding manually assigned peak lists. The structures obtained with NOAH/DIANA are in close agreement with those from manually assigned peak lists, and with both approaches the residual constraint violations correspond to high-quality NMR structure determinations. Systematic comparisons of the bundles of conformers that represent corresponding automatically and interactively determined structures document the absence of significant bias in either approach, indicating that an important step has been made towards automation of structure determination from NMR spectra.

Automated sequence-specific NMR assignment of homologous proteins using the program GARANT.

The program GARANT (General Algorithm for Resonance AssignmeNT) for automated sequence-specific NMR assignment of proteins is based on the mapping of peaks predicted from the amino acid sequence onto the peaks observed in multidimensional spectra [C. Bartels, P. Güntert, M. Billeter and K. Wüthrich (1996) J. Comput. Chem., manuscript submitted for publication]. In this paper we demonstrate the potential of GARANT for the assignment of homologous proteins when either the three-dimensional structure or the chemical shifts of the parent protein are known. In these applications, GARANT utilizes supplementary information either in the form of interatomic distances derived from the three-dimensional structure, in order to add nuclear Overhauser effects reflecting the tertiary structure to the list of expected peaks, or in the form of the chemical shifts of the parent protein, in order to obtain a better estimate of the positions of the expected peaks. The procedure is illustrated with three different proteins: (i) a mutant form of Tendamistat (74 residues), using homonuclear 2D (1)H NMR spectra and either the three-dimensional structure or the chemical shifts of the wild-type protein; (ii) the mutant Antp(C39S, W56S) homeodomain (68 residues), using homonuclear 2D (1)H NMR spectra and the three-dimensional structure of the Antp(C39S) homeodomain; and (iii) free cyclophilin A (165 residues), using heteronuclear 3D NMR spectra and the three-dimensional structure of a cyclophilin A-cyclosporin A complex. In these three systems nearly complete assignment of the polypeptide backbone resonances and assignment of over 80% of the amino acid side-chain resonances was obtained without manual intervention.

Hydration and DNA recognition by homeodomains.

A 2-nanosecond molecular dynamics (MD) simulation of an Antennapedia homeodomain-DNA complex in explicit solvent water at ambient temperature and pressure was performed to supplement experimental nuclear magnetic resonance (NMR) data on the structure and dynamics of this complex. In addition to direct protein-DNA contacts, the MD trajectory attributes an essential role for specific DNA recognition to hydration water molecules that mediate intermolecular contacts. The simulation provides a detailed description of the pathways of hydration water molecules exchanging in and out of the protein-DNA interface and indicates that the residence times of these "interior" waters are on the nanosecond time scale, near the lower end of the range determined by NMR.

Complete relaxation matrix refinement of NMR structures of proteins using analytically calculated dihedral angle derivatives of NOE intensities.

A new method for refining three-dimensional (3D) NMR structures of proteins is described, which takes account of the complete relaxation pathways. Derivatives of the NOE intensities with respect to the dihedral angles are analytically calculated, and efficiently evaluated with the use of a filter technique for identifying the dominant terms of these derivatives. This new method was implemented in the distance geometry program DIANA. As an initial test, we refined 30 rigid distorted helical structures, using a simulated data set of NOE distance constraints for a rigid standard alpha-helix. The final root-mean-square deviations of the refined structures relative to the standard helix were less than 0.1 A, and the R-factors dropped from values between 7% and 32% to values of less than 0.5% in all cases, which compares favorably with the results from distance geometry calculations. In particular, because spin diffusion was not explicitly considered in the evaluation of 'exact' 1H-1H distances corresponding to the simulated NOE intensities, a group of nearly identical distance geometry structures was obtained which had about 0.5 A root-mean-square deviation from the standard alpha-helix. Further test calculations using an experimental NOE data set recorded for the protein trypsin inhibitor K showed that the complete relaxation matrix refinement procedure in the DIANA program is functional also with systems of practical interest.