The solution structure of motilin from NMR distance constraints, distance geometry, molecular dynamics, and an iterative full relaxation matrix refinement

Maximum entropy reconstruction of joint phi, psi-distribution with a coil-library prior: the backbone conformation of the peptide hormone motilin in aqueous solution from phi and psi-dependent J-couplings

In this paper, we present a new method for structure determination of flexible "random-coil" peptides. A numerical method is described, where the experimentally measured 3J(H(alpha)Nalpha) and [3J(H(alpha)Nalpha+1 couplings, which depend on the phi and psi dihedral angles, are analyzed jointly with the information from a coil-library through a maximum entropy approach. The coil-library is the distribution of dihedral angles found outside the elements of the secondary structure in the high-resolution protein structures. The method results in residue specific joint phi,psi-distribution functions, which are in agreement with the experimental J-couplings and minimally committal to the information in the coil-library. The 22-residue human peptide hormone motilin, uniformly 15N-labeled was studied. The 3J(H(alpha)-N(i+1)) were measured from the E.COSY pattern in the sequential NOESY cross-peaks. By employing homodecoupling and an in-phase/anti-phase filter, sharp H(alpha)-resonances (about 5 Hz) were obtained enabling accurate determination of the coupling with minimal spectral overlap. Clear trends in the resulting phi,psi-distribution functions along the sequence are observed, with a nascent helical structure in the central part of the peptide and more extended conformations of the receptor binding N-terminus as the most prominent characteristics. From the phi,psi-distribution functions, the contribution from each residue to the thermodynamic entropy, i.e., the segmental entropies, are calculated and compared to segmental entropies estimated from 15N-relaxation data. Remarkable agreement between the relaxation and J-couplings based methods is found. Residues belonging to the nascent helix and the C-terminus show segmental entropies, of approximately -20 J K(-1) mol(-1) and -12 J K(-1) mol(-1), respectively, in both series. The agreement between the two estimates of the segmental entropy, the agreement with the observed J-couplings, the agreement with the CD experiments, and the assignment of population to sterically allowed conformations show that the phi,psi-distribution functions are indeed meaningful and useful descriptions of the conformational preferences for each residue in this flexible peptide.

Design and synthesis of motilin antagonists derived from the [1-4] fragment of porcine motilin

A series of cyclic peptides having the general structure H-Phe-c[-N(epsilon)-Lys-X-NH-(CH(2))(n)-CO-] were designed on the basis of structure-activity relationship studies of motilin. All were motilin antagonists. The cyclic peptides, in which X is a 3-tert-butyl-substituted tyrosine residue (H-Phe-c[-N(epsilon)-Lys-Tyr(3-tBu)-beta Ala-] (3), H-Phe-c[-N(epsilon)-Lys-Tyr(3-tBu)-Gly-] (6), H-Phe-c[-N(epsilon)-Lys-Tyr(3-tBu)-Abu-] (7), and H-Phe-c[-N(epsilon)-Lys-Tyr(3-tBu)-Ahx-] (8)) showed potent motilin receptor antagonistic activity in the rabbit smooth muscle (pA(2) > 7). The 3-tert-butyl Tyr was found to be the moiety responsible for enhanced binding to the motilin receptor, while the size of the ring had little importance.

Design and synthesis of N-terminal cyclic motilin partial peptides: a novel pure motilin antagonist

Motilin antagonist was designed and synthesized on the basis of the structure-activity relationship analysis of porcine motilin that we reported recently. The drug design was performed on a specific concept to reduce a flexibility of peptide conformation of porcine motilin partial peptide by its cyclization. The cyclic peptide was synthesized using Boc (tert-butyloxycarbonyl) solid phase methodology, followed by cyclization using the azide procedure, and tested for the binding activity to motilin receptor and smooth muscle contractile activity. The cyclic peptides 3, 4, and 5 showed antagonistic property on contraction assay (pA2 [the negative logarithm of molar concentration of antagonist causing a 2-hold shift to the right of the concentration-response curve for motilin]: 4.5, 4.34, and 4.04, respectively, in rabbit duodenum) and no contractile activity even at high concentration.

Analysis of error propagation from NMR-derived internuclear distances into molecular structure of cyclo-pro-gly

Analytical expressions have been derived that translate uncertainties in distance constraints (obtained from NMR investigations) into uncertainties in atom positions in the maximum likelihood (ML) structure consistent with these inputs. As a test of this approach, a comparison was made between test structures reconstructed by the new ML approach, which yields a single structure and a covariance matrix for coordinates, and those reconstructed by metric matrix distance-geometry (MMDG), which yields a family of structures that sample uncertainty space. The test structures used were 560 polyhedra, with edges of arbitrary length containing up to 50 vertices, and one polyhedron, with 100 vertices; randomized distance constraints generated from these structures were used in reconstructing the polyhedra. The uncertainties derived from the two methods showed excellent agreement, and the correlation improved, as expected, with increasingly larger numbers of MMDG structures. This agreement supports the validity of the rapid analytical ML approach, which requires the calculation of only a single structure. As a second test of the ML method, the approach was applied to the determination of uncertainties in the structure of a cyclic dipeptide, cyclo(DL-Pro-Gly) (cPG), derived from NMR cross-relaxation data. The input data were interproton distances calculated from NOEs measured for a solution of the peptide in 2:1 DMSO:H2O at -40 degreesC (so as to yield large negative NOEs). In order to evaluate effects of the quality of the input spectral parameters on the precision of the resulting NMR structure, information from the covalent geometry of cPG was not used in the structure calculations. Results obtained from the analytical ML approach compared favorably with those from the much slower random-walk variant of the Monte Carlo method applied to the same input data. As a third test, the ML approach was used with synthetic structural constraints for a small protein; the results indicate that it will be feasible to use this rapid method to translate uncertainties associated with a given set of distance restraints into uncertainties in atom positions in larger molecules.

A study of melittin, motilin and galanin in reversed micellar environments, using circular dichroism spectroscopy

Circular dichroism spectroscopy has been used to study the behaviour of the cytolytic peptide melittin, the intestinal peptide hormone motilin (porcine) and the neuropeptide galanin (porcine) in various reversed micellar systems. The micellar systems used contained sodium dodecyl sulphate, bis(2-ethylhexyl) sulfosuccinate, n-dodecyltrimethylammonium chloride or polyoxyethylene(7) lauryl ether. Various structural changes of the peptides, induced either by varying the water content or the surface charge of the reversed micelles, could be monitored. Melittin has in all micellar systems a large amount of alpha-helix, and is almost unaffected by both water content and the surface charge of the reversed micelles. Motilin on the other hand attains an alpha-helical structure at low water content only. The surface charges seem to be of importance for the association between motilin and the hydrated reversed micellar surface. Galanin has the most complicated behaviour with a large dependence on surface charge and with a water content dependence which varies with the surfactant used. Stabilization of alpha-helical secondary structures was only seen in negatively charged reversed micelles. These observations indicate a specific interaction between galanin and surfactant, probably of electrostatic nature.

A complete relaxation matrix refinement of the solution structure of a paramagnetic metalloprotein: reduced HiPIP I from Ectothiorhodospira halophila

We have accounted for the effect of paramagnetism on the intensities of NOEs in a 73-residue paramagnetic metalloprotein, the reduced high-potential iron sulfur protein ISO I from Ectothiorhodospira halophila, whose solution structure had been recently solved by us. The paramagnetic effects were dealt with through a suitably modified complete relaxation matrix approach. We have then recalculated the structure through a distance geometry program by minimizing the difference between the sixth roots of the calculated and experimental NOEs. The average RMSD, calculated on residues 4-71, within the structures constituting the two families decreased from 0.67 to 0.46 angstrom for backbone atoms and from 1.23 to 1.06 angstroms for all heavy atoms. The structures in the new family are for the most part within the indetermination of the previous, less resolved, family. A few specific differences are detected and related to the presence of non-negligible paramagnetic effects, which are now properly evaluated.

Structural effects of the selective reduction of amide carbonyl groups in motilin 1-12 as determined by nuclear magnetic resonance

Motilin is a 22-residue peptide stimulating stomach and intestinal motility. The motilin 1-12 fragment displays biological effects similar to the native peptide. Selective reduction of the amide carbonyl groups to form CH2NH analogs leads to a significant reduction in activity for the first two N-terminal positions and to a complete loss of activity for all other positions. The structures of motilin 1-12 and ten reduced analogs were investigated using the temperature dependence of the amide NH chemical shifts. In all the analogs, the structure of the N-terminal region (residues 1-5) was different from the structure of motilin 1-12, which is characterized by hydrogen bonding between Phe1 and Ile4. The structure of the C-terminal region of analogs was similar to the structure of motilin 1-12 for the first two reduction positions only (1-2 and 2-3), indicating that the C-terminal portion of motilin 1-12 is more critical for biological activity. Complete structural characterizations of motilin 1-12, [CH2NH]1-2, and [CH2NH]4-5-motilin 1-12 were performed by two-dimensional NMR spectroscopy and molecular modeling. The structural features observed confirm the differences based on the temperature dependence of the amide NH chemical shifts. These results demonstrate that conservation of the amide bond rigidity is essential for the activity of non-hydrolyzable analogs.

Three-dimensional structure of the highly conserved seventh transmembrane domain of G-protein-coupled receptors

The S/T-X1-X2-N-P-X3-X4-Y highly conserved sequence of the seventh transmembrane (TM VII) segment of G-protein-coupled receptors is not present in the photon receptor bacteriorhodopsin TM VII domain. Despite this noticeable discrepancy in sequence, the X-ray structure of bacteriorhodopsin is generally used as the key structure for modelling all G-protein-coupled receptors. Thus, a kinked trans Pro-helix is usually accepted for the TM VII three-dimensional structure of G-protein-coupled receptors, although Asn-Pro dipeptide mainly induces a type I/III beta-turn conformation in both model peptides and proteins. NMR studies in various solvents and molecular calculations were undertaken in order to gain insight into the conformational behaviour of a 15-residue peptide from the tachykinin NK-1 TM VII domain incorporating this common sequence. The low solubility of this membrane-embedded peptide precludes methanol or micellar systems mimicking membrane environment; thus only dimethylsulfoxide (Me2SO) or chloroform/Me2SO mixture could be used. We also found that perfluoro-tert-butanol, which has not been previously used for NMR studies, constitutes an excellent alternative solvent for the analysis of hydrophobic peptides. The postulated kinked trans-Pro helix was only present as a minor conformer in Me2SO and an equilibrium between helical and extended structures existed. From NOE data a type I/III beta-structure, centered around Pro9-Ile10, probably stabilized by an Asx turn, may be postulated. Addition of chloroform in Me2SO increased the percentage of folded structures but no preferential conformation could be proposed. In perfluoro-tert-butanol/CD3OD (9:1) the N- and C-terminal regions presented an alpha-helical structure, and these two domains were linked by a hinge around Asn-Pro with a gamma-turn for the preceding residue Tyr7 and either a type I/III beta-turn around Pro9-Ile10 or alpha R orientations for these residues, which are both stabilized by an Asx turn. As determined by energy calculations, these structures were equally as stable as the kinked trans-Pro helix and could constitute key structures for analysing the conformational changes and/or the dynamics of TM VII segment induced by the ligand when interacting with the receptor.

Induction of secondary structure in the peptide hormone motilin by interaction with phospholipid vesicles

Motilin is an intestinal peptide hormone that binds to a membrane bound receptor located in the gut tissue. Circular dichroism (CD) was used to study the interaction between either porcine or rabbit motilin or a 1-16 fragment of porcine motilin, with model systems of lipid membranes: sodium dodecyl sulphate (SDS), 1,2-dioleoyl-sn-glycero-3-phosphoglycerol (DOPG) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The CD measurements show significant induction of secondary structure in both motilins and the fragment when negatively charged vesicles (DOPG) or negatively charged micelles (SDS) were present. In contrast, neutral DOPC vesicles did not induce any change in the secondary structure compared to water, in which a random-like secondary structure dominates. The induced secondary structure in the presence of DOPG vesicles is very close to that induced by a mixed aqueous solution containing 30% hexafluoroisopropanol, in which previous NMR-studies have resulted in a three-dimensional solution structure of porcine motilin. In both porcine and rabbit motilin the alpha-helix content is about 50%. This is in agreement with the presence of an amphipathic helix in the C-terminal half of motilin interacting with phospholipid membranes. The interaction appears to be mainly electrostatic in nature, and does not induce any significant alterations in the vesicle, as monitored by EPR studies of spin labels located at the fifth carbon atom of the backbone in a stearic acid molecule. In the 1-16 fragment the alpha-helical content induced by DOPG and SDS is only about 20%.

NMR and circular dichroic studies of the solution structure of conformationally constrained antigenic peptides

Circular dichroic and nuclear magnetic resonance spectroscopies were used to evaluate the conformational properties in solution of a series of 20-amino-acid peptides derived from the primary structure of an antigen from Echinococcus granulosus. The linear peptide corresponding to the sequence 93-112 in the antigen was found to populate in a significant proportion the alpha-helix conformational state. In the presence of 2,2,2-trifluoroethanol, a cosolvent known to stabilize peptide secondary structure, the helical population, estimated from circular dichroic spectra, increases up to 60-70%. Two-dimensional nuclear magnetic resonance studies under these conditions showed that the segment K96-K108 meets all the criteria of an alpha-helix at 281 K and 298 K. Three different variants were synthesized with the same or similar primary structure but containing a lactam-bridged (>) side chain: D107 > K110, D97 > K100 and K94 > E98. Generally, the observed helical content in these variants was lower than in the parent molecule and the stability of the helical conformation decreased in the order D107, K110, K94, E98, D97, K100. Analysis of chemical shift and nuclear Overhauser enhancement data suggested that the lactam rings induce significant distortions of the local features of helix secondary structure. The possible factors of helix destabilization induced by lactam bridges, observed in the studied peptides are discussed in relation to the stabilizing effect of ion pairs in model compounds.

A refined three-dimensional solution structure of a carboxy terminal fragment of apolipoprotein CII

The three-dimensional structure of a synthetic fragment of human apolipoprotein CII (apo-CII) in 35%, 1,1,1,3,3,3-hexafluoro-2-propanol (HFP) has been determined on the basis of distance and intensity constraints derived from two-dimensional proton nuclear magnetic resonance measurements. The NOE crosspeak build-up rates were converted to distance constraints which were used in the distance geometry program DI-ANA. A set of one hundred structures were generated and of these ten structures were used in molecular dynamics simulations using the program XPLOR. This program enabled a direct minimization between the difference of the two-dimensional NOE intensities and those calculated from the full relaxation matrix. In this way spin diffusion is fully taken into account, which can be seen from the considerable improvement of the R-factor after the relaxation matrix refinement. These calculations show that this fragment, which corresponds to the carboxy terminal 30 amino acids of intact apo-CII and which retains its ability to activate lipoprotein lipase, is essentially flexible, but has three defined secondary structural elements. The most significant one is an alpha-helix between residues 67 and 74. The following three residues adopt a turn-like structure. Another turn of alpha-helix is seen between residues 56 and 59. The effect of the solvent system on the secondary structure was studied by circular dichroism spectroscopy. The results show that the mixed aqueous 35% HFP solvent induces secondary structure of a very similar nature to the one induced by sodium dodecyl sulphate.

Structures of DNA-binding mutant zinc finger domains: implications for DNA binding

Studies of Cys2-His2 zinc finger domains have revealed that the structures of individual finger domains in solution determined by NMR spectroscopy are strikingly similar to the structure of fingers bound to DNA determined by X-ray diffraction. Therefore, detailed structural analyses of single finger domains that contain amino acid substitutions known to affect DNA binding in the whole protein can yield information concerning the structural ramifications of such mutations. We have used this approach to study two mutants in the N-terminal finger domain of ADR1, a yeast transcription factor that contains two Cys2-His2 zinc finger sequences spanning residues 102-159. Two point mutants at position 118 in the N-terminal zinc finger (ADR1b: 102-130) that adversely affect the DNA-binding activity of ADR1 have previously been identified: H118A and H118Y. The structures of wild-type ADR1b and the two mutant zinc finger domains were determined using two-dimensional nuclear magnetic resonance spectroscopy and distance geometry and were refined using a complete relaxation matrix method approach (REPENT) to improve agreement between the models and the nuclear Overhauser effect spectroscopy data from which they were generated. The molecular architecture of the refined wild-type ADR1b domain is presented in detail. Comparisons of wild-type ADR1b and the two mutants revealed that neither mutation causes a significant structural perturbation. The structures indicate that the DNA binding properties of the His 118 mutants are dependent on the identity of the side chain at position 118, which has been postulated to make a direct DNA contact in the wild-type ADR1 protein. The results suggest that the identity of the side chain at the middle DNA contact position in Cys2-His2 zinc fingers may be changed with impunity regarding the domain structure and can affect the affinity of the protein-DNA interaction.

Structure and dynamics of motilin. Time-resolved fluorescence of peptide hormone with single tyrosine residue

Time-resolved fluorescence and CD spectroscopy were used to characterize the structure and dynamics of the peptide hormone motilin with a single tyrosine residue among its 22 amino acids. CD spectroscopy showed that secondary structure is independent of concentration in the range 1 x 10(-5)-2.6 x 10(-4) M, and of the presence of DOPC lipid vesicles, but is strongly induced by addition of hexafluoroisopropanol. The fluorescence studies with tyrosine as the intrinsic fluorophore, performed at the MAX synchrotron laboratory at Lund, showed that three fluorescence lifetimes (0.4 ns, 1.7 ns and 3.6 ns at 20 degrees C) and two rotational correlation times (0.4 ns and 5 ns at 20 degrees C) were needed to account for the data. The different decay times are interpreted as representing ground-state rotamers interconverting slowly on the ns time scale. The rotational correlation times are ascribed to local angular motion of the tyrosyl ring, and global motion of the whole peptide, respectively.

D-amino acid and alanine scans of the bioactive portion of porcine motilin

A recent systematic study of porcine motilin fragments has clearly shown that biological activity resides in the amino-terminal end. The amino-terminal tetradecapeptide retains more than 90% of the potency of the full molecule. We now examined the effect of replacement of residues 1 through 11 by either their D-isomer or by alanine in [Leu13]pMOT(1-14). Peptides were synthesized using Fmoc solid phase methodology, purified by HPLC, and assayed for their ability to displace bound motilin (rabbit antral smooth muscle homogenate) and to induce contractions (isolated rabbit duodenal segments). The negative logarithm of the concentration displacing 50% of the tracer (pIC50), or producing 50% of the maximal contractile response (pEC50), was determined. All compounds were still full agonists. A reduction in potency of more than two log units was seen for the compounds in which residues 1 (Phe), 4 (Ile), and 7 (Tyr) were replaced by Ala and residues 3 (Pro), 4 (Ile), and 6 (Thr) by their D-isomer. The largest drop was noted for the analogs substituted at position 4. For all compounds there was an almost perfect correlation between the pIC50 and the pEC50 values (r = 0.96), although the pEC50 was consistently smaller. These results show that the biological activity of motilin is mainly determined by the first seven residues. The pharmacophore consists of the aromatic rings from Phe1 and Tyr7 and the aliphatic side chains from Val2 and Ile4. Pro3, Phe5, and Thr6 may stabilize the bioactive conformation.