Dihedral Restraints Research Articles

High-resolution solution structure of the EGF-like domain of heregulin-alpha.

The solution structure of the 63-residue heregulin-alpha (HRG-alpha) epidermal growth factor (EGF)-like domain, corresponding to residues 177-239 of HRG-alpha, has been determined to high resolution using data from two-dimensional and three-dimensional homo- and heteronuclear NMR spectroscopy. The structure is based on a total of 887 internuclear distance and dihedral restraints derived from data obtained using unlabeled and uniformly 15N-labeled protein samples, at pH 4.5, 20 degrees C. A total of 20 structures were calculated using a hybrid distance geometry-simulated annealing approach with the program DGII, followed by restrained molecular dynamics using the program DISCOVER. The average maximum violations are 0.12 +/- 0.01 angstroms and 1.4 +/- 0.3 degrees for distance and dihedral restraints, respectively. The backbone (N,C(alpha),C) atomic rms distribution about the mean coordinates for residues 3-23 and 31-49 is 0.29 +/- 0/07 angstroms. The N-and C-terminal residues (1-2 and 50-63) and 24-30 are disordered. Comparison of the HRG-alpha EGF-like domain structure with the previously determined structure of human EGF [Hommel et al. (1992) J. Mol. Biol. 227, 271-282] reveals a high degree of structural similarity; excluding the N-terminal region (residues 1-13), the disordered phi-loop region (residues 24-30) that contains a three-residue insertion in HRG-alpha relative to hEGF, and the disordered C-terminal region (residues 50-63), the C(alpha) alignment between the HRG-alpha and hEGF minimized mean structures has a rms difference of approximately 1 angstrom. In HRG-alpha the N-terminal residues 2-6 form a well-defined beta strand rather than being disordered as found for hEGF. This structural difference correlates with functional data which suggest that the N-terminal region of the HRG-alpha EGF-like domain is responsible for the observed receptor specificity differences between HRG-alpha and EGF.

Structures of a series of 6-kDa trypsin inhibitors isolated from the stigma of Nicotiana alata.

The three-dimensional structures of a series of 6-kDa trypsin inhibitors isolated from the stigma of the ornamental tobacco Nicotiana alata have been determined by 1H NMR spectroscopy combined with simulated annealing calculations. The proteins, T1-T4, are proteolytically cleaved from a 40.3-kDa precursor protein, NA-proPI, together with a chymotrypsin inhibitor, C1, the structure of which was reported recently [Nielsen, K.J., Heath, R.L., Anderson, M.A., & Craik, D.J. (1994) J. Mol. Biol. 242, 231-243]. Each of the proteinase inhibitors comprises 53 amino acids, including 8 cysteine residues which are linked to form 4 disulfide bridges. The proteins have a high degree of sequence identity and differ mainly in residues around the putative reactive sites. The structure of T1 was determined using a set of 533 interproton distance restraints derived from NOESY spectra, combined with 33 dihedral restraints derived from 3JNH-H alpha coupling constants and 16 hydrogen bonds. The structures of the remaining inhibitors (T2-T4) were deduced to be almost identical to T1, on the basis of their similar chemical shifts and 2D spectra. The current study demonstrates that the structures of the trypsin inhibitors (T1-T4) are similar to that previously found for the chymotrypsin inhibitor, C1. Despite differences in sequence, there is conservation in backbone geometry between the reactive site loops of the two classes of inhibitors. From this, it is clear that the nature of the side chain on the primary binding residue, rather than the backbone fold, is the main determinant of the enzyme specificities of these proteinase inhibitors.

A genetic algorithm-based protocol for docking ensembles of small ligands using experimental restraints

A genetic algorithm (GA) based method for docking ensembles of small, flexible ligands to receptor proteins using NMR-derived constraints is described. In this method, three translations and rotations of the ligand and the dihedral angles of the ligand are represented by binary strings and evolve under the genetic operators of cross-over, mutation, migration and selection. The fitness function for the selection process includes distance and dihedral restraints and a repulsive van der Waals term. The GA was applied to a three-atom model system as well as to the streptavidin-biotin complex using simulated intermolecular distance restraints. In both systems, the GA was able to obtain low-energy conformations when only a single binding site was simulated. Calculations were also performed using distance restraints from two distinct binding sites. In these simulations, the GA was able to obtain low-energy conformations corresponding to ligand molecules in each of the two sites. The inclusion of additional ligands in the ensemble did not result in an energetic benefit, confirming that only two ligand conformations were necessary to fulfill the distance restraints. This method allows for a direct investigation of the minimum number of ligand orientations necessary to fulfill experimental distance restraints, and simultaneously yields detailed structural information about each site.

Resonance assignments and solution structure of the second RNA-binding domain of sex-lethal determined by multidimensional heteronuclear magnetic resonance.

The RNA-binding protein Sex-lethal (Sxl) is a critical regulator of sexual differentiation and dosage compensation in Drosophila. This regulatory activity is a consequence of the ability of Sxl to bind uridine-rich RNA tracts involved in pre-mRNA splicing. Sxl contains two RNP consensus-type RNA-binding domains (RBDs). A structural study of a portion of Sxl (amino acids 199-294) containing the second RNA-binding domain (RBD-2) using multidimensional heteronuclear NMR is presented here. Nearly complete 1H, 13C, and 15N resonance assignments have been obtained from 15N- and 13C/15N-uniformly labeled protein. These assignments were used to analyze 3D 15N-separated NOESY and 13C/13C-separated 4D NOESY spectra which produced 494 total and 169 long-range NOE-derived distance restraints. Along with 41 backbone dihedral restraints, these distance restraints were employed to generate an intermediate-resolution family of calculated structures, which exhibits the beta alpha beta-beta alpha beta tertiary fold found in other RBD-containing proteins. The RMSD to the average structure for the backbone atoms of residues 11-93 is 1.55 +/- 0.30 A, while the RMSD for backbone atoms involved in secondary structure is 0.76 +/- 0.14 A. A capping box [Harper, E.T., & Rose, G.D. (1993) Biochemistry 32, 7605-7609] was identified at the N-terminus of the first helix and has been characterized by short- and medium-range NOEs. Finally, significant structural similarities and differences between Sxl RBD-2 and other RBD-containing proteins are discussed.

Sequential assignment and structure determination of spider toxin .omega.-Aga-IVB

The solution structure of a peptide toxin isolated from funnel web spider venom, omega-Aga-IVB, was determined by 2D NMR methods. omega-Aga-IVB is a high-affinity specific blocker of P-type voltage-dependent calcium channels. Nearly all of the proton resonances of this 48-residue protein were assigned using conventional 2D homonuclear NMR experiments. The three-dimensional structure of the molecule was determined by simulated annealing. The distance and dihedral restraints used in the structure calculations were derived from NOESY and COSY-type experiments, respectively. Mass spectrometric analysis of omega-Aga-IVB suggests that the protein contains four disulfide bonds. In the absence of chemical data to identify the pattern of cysteine pairing, the disulfide bonds of the toxin are proposed from the NMR data and subsequent structural calculations. The structure of the toxin can be described as a three-stranded anti-parallel beta sheet connected by flexible loops. A striking feature of the structure is that the C-terminal 10 residues of this protein adopt random coil conformations. Several positively charged amino acid side chains are found localized on one face of the molecule, in close proximity to the C-terminal tail. This observation has led us to propose a speculative model of the toxins blockade mechanism.

Human epidermal growth factor. High resolution solution structure and comparison with human transforming growth factor alpha.

The solution structure of the 53 amino acid peptide hormone, human epidermal growth factor (hEGF), has been determined to high resolution from nuclear magnetic resonance (n.m.r.) data. A large number of internuclear distance and dihedral restraints was obtained, including data from uniformly 15N-labelled hEGF. Dynamical simulated annealing methods using the program XPLOR were used for structure calculation. An improved protocol was developed combining efficient conformational searching at a reduced computational cost. The general fold of the calculated structures compared well with that of a derivative of the carboxy-terminally truncated hEGF determined previously. A group of 44 structures were calculated with no violations greater than 0.3 A and 3 degrees for distance and dihedral restraints, respectively. The average pairwise root mean square (r.m.s.) deviation of all backbone atoms for these structures was 2.25 A for all 53 residues, 0.92 A for the bulk of the protein, and 0.23 A for the functionally important carboxy-terminal domain. Two new helical segments containing highly conserved amino acids have been identified; one between cysteines 6 and 14 and a second at the end of the carboxy-terminal domain. New insight into the molecular architecture of the site of putative receptor binding was provided by comparing the structure of hEGF with its biologically equipotent analogue, human transforming growth factor alpha. This comparison revealed a close structural relationship between the two growth factors and provides an improved understanding of the structure/function relationships in EGF.

Two-dimensional 1H NMR study of recombinant insect defensin A in water: resonance assignments, secondary structure and global folding.

A 500 MHz 2D 1H NMR study of recombinant insect defensin A is reported. This defense protein of 40 residues contains 3 disulfide bridges, is positively charged and exhibits antibacterial properties. 2D NMR maps of recombinant defensin A were fully assigned and secondary structure elements were localized. The set of NOE connectivities, 3JNH-alpha H coupling constants as well as 1H/2H exchange rates and delta delta/delta T temperature coefficients of NH protons strongly support the existence of an alpha-helix (residues 14-24) and of an antiparallel beta-sheet (residues 27-40). Models of the backbone folding were generated by using the DISMAN program and energy refined by using the AMBER program. This was done on the basis of: (i) 133 selected NOEs, (ii) 21 dihedral restraints from 3JNH-alpha H coupling constants, (iii) 12 hydrogen bonds mostly deduced from 1H/2H exchange rates or temperature coefficients, in addition to 9 initial disulfide bridge covalent constraints. The two secondary structure elements and the two bends connecting them involve approximately 70% of the total number of residues, which impose some stability in the C-terminal part of the molecule. The remaining N-terminal fragment forms a less well defined loop. This spatial organization, in which a beta-sheet is linked to an alpha-helix by two disulfide bridges and to a large loop by a third disulfide bridge, is rather similar to that found in scorpion charybdotoxin and seems to be partly present in several invertebrate toxins.

Conformation of secretin in dimethyl sulfoxide solution. NMR studies and restrained molecular dynamics.

The solution conformation of the 27-residue polypeptide hormone secretin in dimethyl sulfoxide has been determined on the basis of 1H-NMR measurements. The experimental data set used in the structure determination consisted of 98 nuclear-Overhauser-enhancement-derived interproton and dihedral angle restraints from coupling constants. The NH-NH and H alpha-NH NOEs were determined from build-up rates, while the remaining distances were classified in a qualitative manner. The structure calculations consisted of two phases. First, dynamical simulated annealing calculations were carried out to find conformations of the peptide which satisfy NOE and phi dihedral restraints. The convergence of ten calculated structures was good except for those regions of the molecule where NOE data were not unambiguous. From the calculated set another initial structure was built which was again minimized in several 5-ps calculations now employing the full empirical energy function. The resulting structures of secretin reveal conformationally well-defined regions, but not a single uniform secondary structure. The structure is different from the calculated structure from trifluoroethanol/water measurements.