Prediction of protein structure: the problem of fold multiplicity.

Three-dimensional (3D) models of four CASP3 targets were calculated using a simple modeling procedure that includes prediction of regular secondary structure, analysis of possible beta-sheet topologies, assembly of amphiphilic helices and beta-sheets to bury their nonpolar surfaces, and adjustment of side-chain conformers and loops to provide close packing and saturation of the "hydrogen bond potential" (exposure of all polar groups to water or their involvement in intramolecular hydrogen bonds). It has been found that this approach allows construction of 3D models that, in some cases, properly reproduce the structural class of the protein (such as beta-barrel or beta-sandwich of definite shape and size) and details of tertiary structure (such as pairing of beta-strands), although all four models were more or less incorrect. Remarkably, some models had fewer water-exposed nonpolar side-chains, more hydrogen bonds, and smaller holes than the corresponding native structures (although the models had a larger water-accessible nonpolar surface). The results obtained indicate that hydrophobicity patterns do not unequivocally determine protein folds, and that any ab initio or fold recognition methods that operate with imprecise potential energy functions, or use crude geometrical approximations of the peptide chain, will probably produce many different nonnative structures.

Structural organization of G-protein-coupled receptors.

Atomic-resolution structures of the transmembrane 7-alpha-helical domains of 26 G-protein-coupled receptors (GPCRs) (including opsins, cationic amine, melatonin, purine, chemokine, opioid, and glycoprotein hormone receptors and two related proteins, retinochrome and Duffy erythrocyte antigen) were calculated by distance geometry using interhelical hydrogen bonds formed by various proteins from the family and collectively applied as distance constraints, as described previously [Pogozheva et al., Biophys. J., 70 (1997) 1963]. The main structural features of the calculated GPCR models are described and illustrated by examples. Some of the features reflect physical interactions that are responsible for the structural stability of the transmembrane alpha-bundle: the formation of extensive networks of interhelical H-bonds and sulfur-aromatic clusters that are spatially organized as 'polarity gradients'; the close packing of side-chains throughout the transmembrane domain; and the formation of interhelical disulfide bonds in some receptors and a plausible Zn2+ binding center in retinochrome. Other features of the models are related to biological function and evolution of GPCRs: the formation of a common 'minicore' of 43 evolutionarily conserved residues; a multitude of correlated replacements throughout the transmembrane domain; an Na(+)-binding site in some receptors, and excellent complementarity of receptor binding pockets to many structurally dissimilar, conformationally constrained ligands, such as retinal, cyclic opioid peptides, and cationic amine ligands. The calculated models are in good agreement with numerous experimental data.

Opioid receptor three-dimensional structures from distance geometry calculations with hydrogen bonding constraints.

Three-dimensional structures of the transmembrane, seven alpha-helical domains and extracellular loops of delta, mu, and kappa opioid receptors, were calculated using the distance geometry algorithm, with hydrogen bonding constraints based on the previously developed general model of the transmembrane alpha-bundle for rhodopsin-like G-protein coupled receptors (Biophys. J. 1997. 70:1963). Each calculated opioid receptor structure has an extensive network of interhelical hydrogen bonds and a ligand-binding crevice that is partially covered by a beta-hairpin formed by the second extracellular loop. The binding cavities consist of an inner "conserved region" composed of 18 residues that are identical in delta, mu, and kappa opioid receptors, and a peripheral "variable region," composed of 19 residues that are different in delta, mu, and kappa subtypes and are responsible for the subtype specificity of various ligands. Sixteen delta-, mu-, or kappa-selective, conformationally constrained peptide and nonpeptide opioid agonists and antagonists and affinity labels were fit into the binding pockets of the opioid receptors. All ligands considered have a similar spatial arrangement in the receptors, with the tyramine moiety of alkaloids or Tyr1 of opioid peptides interacting with conserved residues in the bottom of the pocket and the tyramine N+ and OH groups forming ionic interactions or H-bonds with a conserved aspartate from helix III and a conserved histidine from helix VI, respectively. The central, conformationally constrained fragments of the opioids (the disulfide-bridged cycles of the peptides and various ring structures in the nonpeptide ligands) are oriented approximately perpendicular to the tyramine and directed toward the extracellular surface. The results obtained are qualitatively consistent with ligand affinities, cross-linking studies, and mutagenesis data.

Thermodynamic model of secondary structure for alpha-helical peptides and proteins.

A thermodynamic model describing formation of alpha-helices by peptides and proteins in the absence of specific tertiary interactions has been developed. The model combines free energy terms defining alpha-helix stability in aqueous solution and terms describing immersion of every helix or fragment of coil into a micelle or a nonpolar droplet created by the rest of protein to calculate averaged or lowest energy partitioning of the peptide chain into helical and coil fragments. The alpha-helix energy in water was calculated with parameters derived from peptide substitution and protein engineering data and using estimates of nonpolar contact areas between side chains. The energy of nonspecific hydrophobic interactions was estimated considering each alpha-helix or fragment of coil as freely floating in the spherical micelle or droplet, and using water/cyclohexane (for micelles) or adjustable (for proteins) side-chain transfer energies. The model was verified for 96 and 36 peptides studied by 1H-nmr spectroscopy in aqueous solution and in the presence of micelles, respectively ([set 1] and [set 2]) and for 30 mostly alpha-helical globular proteins ([set 3]). For peptides, the experimental helix locations were identified from the published medium-range nuclear Overhauser effects detected by 1H-nmr spectroscopy. For sets 1, 2, and 3, respectively, 93, 100, and 97% of helices were identified with average errors in calculation of helix boundaries of 1.3, 2.0, and 4.1 residues per helix and an average percentage of correctly calculated helix-coil states of 93, 89, and 81%, respectively. Analysis of adjustable parameters of the model (the entropy and enthalpy of the helix-coil transition, the transfer energy of the helix backbone, and parameters of the bound coil), determined by minimization of the average helix boundary deviation for each set of peptides or proteins, demonstrates that, unlike micelles, the interior of the effective protein droplet has solubility characteristics different from that for cyclohexane, does not bind fragments of coil, and lacks interfacial area.

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints.

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.

Development of a model for the delta-opioid receptor pharmacophore. 4. Residue 3 dehydrophenylalanine analogues of Tyr-c[D-Cys-Phe-D-Pen]OH (JOM-13) confirm required gauche orientation of aromatic side chain.

We have previously proposed a model for the delta-opioid receptor binding conformation of the high affinity tetrapeptide Tyr-c[D-Cys-Phe-D-Pen]OH (JOM-13) based on experimental and theoretical conformational analysis of this peptide and a correlation of conformational preferences of further conformationally restricted analogues of this tetrapeptide with their receptor binding affinities. A key element of this model is the requirement that the Phe3 side chain exist in the chi 1 = -60 degrees conformation. Conformational calculations on the residue 3 dehydrophenylalanine analogues of JOM-13 suggest that while the dehydro (Z) phenylalanine analogue can be superimposed easily with the proposed binding conformer of JOM-13, the dehydro(E)phenylalanine analogue cannot. These results lead to the prediction that the dehydro(Z)phenylalanine analogue should display similar delta-receptor binding affinity as JOM-13 while the dehydro(E)phenylalanine analogue is expected to bind less avidly. Synthesis and subsequent opioid receptor binding analysis of the dehydrophenylalanine analogues of JOM-13 confirm these predictions, lending support to the delta-pharmacophore model.

Development of a model for the delta-opioid receptor pharmacophore: 3. Comparison of the cyclic tetrapeptide, Tyr-c[D-Cys-Phe-D-Pen]OH with other conformationally constrained delta-receptor selective ligands.

We have previously proposed a model of the delta-opioid receptor bound conformation for the cyclic tetrapeptide, Tyr-c[D-Cys-Phe-D-Pen]OH (JOM-13) based on its conformational analysis and from conformation-affinity relationships observed for its analogues with modified first and third residues. To further verify the model, it is compared here with results of conformational and structure-activity studies for other known conformationally constrained delta-selective ligands: the cyclic pentapeptide agonist, Tyr-c[D-Pen-Gly-Phe-D-Phe]OH (DPDPE): the peptide antagonist, Tyr-Tic-Phe-PheOH (TIPP); the alkaloid agonist, 7-spiroindanyloxymorphone (SIOM); and the related alkaloid antagonist, oxymorphindole (OMI). A candidate delta-bound conformer is identified for DPDPE that provides spatial overlap of the functionally important N-terminal NH3+ and C-terminal COO- groups and the aromatic rings of the Tyr and Phe residues in both cyclic peptides. It is shown that all delta-selective ligands considered have similar arrangements of their pharmacophoric elements, i.e., the tyramine moiety and a second aromatic ring (i.e., the rings of Phe3, Phe4, and Tic2 residues in JOM-13, DPDPE, and TIPP, respectively; the indole ring system in OMI, and the indanyl ring system in SIOM). The second aromatic rings, while occupying similar regions of space throughout the analogues considered, have different orientations in agonists and antagonists, but identical orientations in peptide and alkaloid ligands with the same agonistic or antagonistic properties. These results agree with the previously proposed binding model for JOM-13, are consistent with the view that delta-opioid agonists and antagonists share the same binding site, and support the hypothesis of a similar mode of binding for opioid peptides and alkaloids.

Development of a model for the delta opioid receptor pharmacophore. 1. Conformationally restricted Tyr1 replacements in the cyclic delta receptor selective tetrapeptide Tyr-c[D-Cys-Phe-D-Pen]OH (JOM-13).

A series of analogues of the conformationally restricted delta opioid receptor selective tetrapeptide Tyr-c[D-Cys-Phe-D-Pen]OH (JOM 13) was prepared in which the conformationally labile Tyr residue was replaced with several less flexible tyrosine analogues. Among these tyrosine analogues were the bicyclic structures 1,2,3,4-tetrahydro-7-hydroxyisoquinoline-3-carboxylic acid (HO-Tic), 2-amino-6-hydroxytetralin-2-carboxylic acid (Hat), and 2-amino-5-hydroxyindan-2-carboxylic acid (Hai) in which rotations about the C alpha-C beta and C beta-C gamma bonds are restricted due to cyclization of the side chain to the backbone. Also examined were analogues in which tyrosine was replaced with either trans-3-(4'-hydroxyphenyl)proline (t-Hpp) or cis-3-(4'-hydroxyphenyl)proline (c-Hpp), residues in which rotations about C alpha-C beta, but not C beta-C gamma, are restricted. Both the t-Hpp1 and c-Hpp1 analogues displayed delta receptor binding affinity similar to the parent Tyr1-containing peptide, while the D-Hat1, L-Hat1, and L-Hai1 analogues exhibited somewhat lower affinity. The results observed for the t-Hpp1 and c-Hpp1 analogues are particularly significant since these two residues have little accessible conformational space in common. Since the binding conformation of residue 1 must be included in this limited conformational intersection, its elucidation is facilitated. Bioassay results from guinea pig ileum and mouse vas deferens preparations are in general agreement with the binding results; however some potency discrepancies are observed. These discrepancies may reflect different selectivities among delta receptor subtypes for the analogues or may represent differing efficacies among these conformationally restricted peptides. The conformational properties of the parent tetrapeptide and the residue 1-modified analogues were studied by molecular mechanics computations. All these peptides share a common rigid tripeptide cycle with a single energetically preferred backbone conformation and three different conformers of the D-Cys, D-Pen disulfide bridge, two of which are observed in the solid state and in aqueous solution, as previously determined from X-ray crystallography and 1H NMR spectroscopy data (Lomize, A; et al. J. Am. Chem. Soc. 1994, 116, 429-436). All the peptides have similar sets of low-energy conformations of their common flexible elements, the Phe3 side chain and the peptide group between the first residue and the rigid tripeptide cycle. However, possible conformations of the first residue differ and depend on the covalent constraints incorporated into the side chain.(ABSTRACT TRUNCATED AT 400 WORDS)

Two-dimensional 1H-NMR study of the spatial structure of neurotoxin II from Naja naja oxiana.

The spatial structure of neurotoxin II from the venom of the central Asian cobra Naja naja oxiana was determined by two-dimensional 1H-NMR techniques and computational analysis. Nearly complete proton resonance assignments for 61 amino acid residues have been made using two-dimensional (2D) homonuclear total correlated spectroscopy, 2D homonuclear double-quantum-filtered correlated spectroscopy and 2D homonuclear NOE spectroscopy (NOESY) experiments. The cross-peak volumes in NOESY spectra spin-spin coupling constants of vicinal protons NH-C alpha H and C alpha H-C beta H and the observation of slow deuterium exchange of amide protons were used to define local structure and a set of constraints for distance geometry program DIANA. The average root-mean-square deviations are 53 pm for backbone heavy atoms and 118 pm for all heavy atoms of 19 final neurotoxin II conformations. The spatial structure is characterized by a short double-stranded (residues 1-5 and 13-17) and a triple-stranded (residues 22-30, 33-41 and 50-54) antiparallel beta-sheets.

Spatial structure of (34-65)bacterioopsin polypeptide in SDS micelles determined from nuclear magnetic resonance data.

The spatial structure of a synthetic 32-residue polypeptide, an analog of the membrane-spanning segment B (residues 34-65) of bacterioopsin of Halobacterium halobium, incorporated into perdeuterated sodium dodecyl sulfate micelles, was determined from 1H NMR data. The structure determination included the following steps: (1) local structure analysis; (2) structure calculations using the distance geometry program DIANA; (3) systematic search for energetically allowed side-chain rotamers consistent with NOESY cross-peak volumes; (4) random generation of peptide conformations in allowed conformational space. The obtained structure has a right-handed alpha-helical region from Lys41 to Leu62 with a kink of 27 degrees at Pro50. The C-cap Gly63 adopts a conformation with phi = 87 +/- 6 degrees, psi = 43 +/- 10 degrees typical to a left-handed helix. The N-terminal part (residues 34-40) is exposed to the aqueous phase and lacks an ordered conformation. The secondary structure of segment B in micelles is consistent with the high-resolution electron cryomicroscopy model of bacteriorhodopsin (Henderson et al. (1990) J. Mol. Biol., 213, 899-929).

Three-dimensional structure of proteolytic fragment 163-231 of bacterioopsin determined from nuclear magnetic resonance data in solution.

546 NOESY cross-peak volumes were measured in the two-dimensional NOESY spectrum of proteolytic fragment 163-231 of bacterioopsin in organic solution. These data and 42 detected hydrogen bonds were applied for determining the peptide spatial structure. The fold of the polypeptide chain was determined by local structure analysis, a distance geometry approach and systematic search for energetically allowed side-chain rotamers which are consistent with experimental NOESY cross-peak volumes. The effective rotational correlation time of 6 ns for the molecule was evaluated from optimization of the local structure to meet NOE data and from the dependence on mixing time of the NiH/Ci alpha H cross-peak volumes of the residues in alpha-helical conformation. The resulting structure has two well defined alpha-helical regions, 168-191 and 198-227, with root-mean-square deviation 44 pm and 69 pm, respectively, between the backbone atoms in 14 final energy refined conformations. The alpha-helices correspond to transmembrane segments F and G of bacteriorhodopsin. The segment F contains proline 186, which introduces a kink of about 25 degrees with a disruption of the hydrogen bond with the NH group of the following residue. The segments are connected by a flexible loop region 192-197. Torsion angles chi 1 are unequivocally defined for 62% of side chains in the alpha-helices but half of them differ from electron cryo-microscopy (ECM) model of bacteriorhodopsin, apparently because of the low resolution of ECM. Nevertheless, the F and G segments can be packed as in the ECM model and with side-chain conformations consistent with all NMR data in solution.

[Refinement of the spatial structure of the gramicidin A ion channel]. / Utochnenie prostranstvennoi struktury ionnogo kanala gramitsidina A.

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.

[Determination of the spatial structure of insectotoxin 15A from Buthus erpeus by (1)H-NMR spectroscopy data]. / Opredelenie prostranstvennoi struktury insektotoksina 15A Buthus eupeus po dannym spektroskopii (1)H-IaMR.

The solution structure of insectotoxin 15A (35 residues) from scorpion Buthus eupeus was determined on the basis of 386 interproton distance restraints 12 hydrogen-bonding restraints and 113 dihedral angle restraints derived from 1H NMR experiments. A group of 20 structures was calculated with the distance geometry program DIANA followed by the restrained energy minimization with the program CHARMM. The atomic RMS distribution about the mean coordinate position is 0.64 +/- 0.11 A for the backbone atoms and 1.35 +/- 0.20 A for all atoms. The structure contains an alpha-helix (residues 10-20) and a three-stranded antiparallel beta-sheet (residues 2-5, 24-28 and 29-33). A pairing of the eight cysteine residues of insectotoxin 15A was established basing on NMR data. Three disulfide bridges (residues 2-19, 16-31 and 20-33) connect the alpha-helix with the beta-sheet, and the fourth one (5-26) joins beta-strands together. The spatial fold of secondary structure elements (the alpha-helix and the beta-sheet) of the insectotoxin 15A is very similar to those of the other short and long scorpion toxins in spite of a low (about 20%) sequence homology.

[Conformational analysis of a segment in bacterioopsin by two-dimensional (1)H-NMR spectroscopy]. / Konformatsionnyi analiz segmenta v bakterioopsina po dannym dvumernoispektroskopii (1)H-IaMR.

The spatial structure of a synthetic peptide, an analogue of the membrane spanning segment B (residues 34-65) of bacterioopsin from Halobacterium halobium, has been refined. Backbone torsion angles were derived from intensities of short-range interproton NOEs. These, together with a complete set of the NOEs integral intensities formed the basis for the three-dimensional structure refinement by the energy minimization with consideration of NOE penalty functions. Analysis indicates the right-handed alpha-helical conformation of segment B extending from Asp-38 to Tyr-64 with a kink of the helical axis (27 degrees) at Pro-50. The most stable region with an average root-mean-square deviation of 0.43 A between the backbone atoms includes residues 42-60 in six energy refined structures. The N-terminal part of segment B (residues 34-37) has no ordered conformation. The inferred structure is in close agreement with the electron cryomicroscopy structure of bacteriorhodopsin, differing from it in conformations of most of the side chains.

[Determination of the local structure of the protein insectotoxin I5A from the scorpion Buthus eupeus from 1H-NMR spectroscopy data]. / Opredelenie lokal'noi struktury belkovogo insektotoksina I5A skorpiona Buthus eupeus po dannym spektroskopii 1H-IaMR.

The local structure (torsion angles phi, psi and chi 1 of amino acid residues) of insectotoxin I5A (35 residues) of scorpion Buthus eupeus has been determined from cross-peak integral intensities in two-dimensional nuclear Overhauser enhancement (NOESY) spectra and spin coupling constants of vicinal H--NC alpha--H and H--C alpha C beta--H protons. The local structure determination was carried out by fitting complete relaxation matrix of peptide unit protons (protons of a given residue and NH proton of the next residue in the amino acid sequence) with experimental NOESY cross-peak intensities. The obtained intervals of backbone torsional angles phi and psi consistent with NMR data were determined for all but Gly residues. The predominant C alpha--C beta rotamer of the side chain has been unambiguously determined for 42% of the insectotoxin amino acid residues whereas for another 46% residues experimental data are fitted equally well with two rotamers. Stereospecific assignments were obtained for 38% of beta-methylene groups. The determined torsional angles phi, psi and chi 1 correspond to the sterically allowed conformations of the amino acid residues and agree with the insectotoxin secondary structure established earlier by 1H NMR spectroscopy.

[Determination of local conformation of proteins from 1H-NMR spectroscopy data]. / Metodika opredeleniia lokal'noi konformatsii belkov po dannym spektroskopii 1H-IaMR.

A method is proposed to determine conformations of amino acid residues of the protein and effective correlation time tau c from cross-peak intensities in two-dimensional nuclear Overhauser enhancement (NOESY) spectra. The method consists in fitting complete relaxation matrix of dipeptide unit protons to experimental cross-peak intensities by varying phi, psi, chi torsional angles and tau c. To verify the method, NOESY spectra of basic pancreatic trypsin inhibitor (BPTI) were theoretically generated at mixing times tau m = 25-300 ms and tau c = 4 ns and used for local structure determination. The method works well with optimum for measurement of NOE intensities tau m 100-200 ms. As a result, the backbone phi, psi torsion angles were unambiguously determined at tau m = 100 ms for all but Gly residues of BPTI, and chi 1 angles were determined for the majority of side chains. The obtained dipeptide unit conformations are very close to the BPTI crystallographic structure: root mean square deviation (RMSD) of interproton distances within dipeptide units, on the average, is 0.08 A (maximal deviation 0.44 A), and RMSD of phi and psi angles are 18 and 9 degrees, respectively (maximal deviations are 44 and 22 degrees).

1H-NMR study of gramicidin A transmembrane ion channel. Head-to-head right-handed, single-stranded helices.

The structure of [Val1]gramicidin A incorporated into sodium dodecyl-d25 sulphate micelles has been studied by two-dimensional proton NMR spectroscopy. Analysis of nuclear Overhauser effects, spin-spin couplings and solvent accessibility of NH groups show that the conformation of the Na+ complex of gramicidin A in detergent micelles, which in many ways mimic the phospholipid bilayer of biomembranes, is an N-terminal to N-terminal (head-to-head) dimer (Formula: see text) formed by two right-handed, single-stranded beta 6.3 helices with 6.3 residues per turn, differing from Urry's structure by handedness of the helices.

[Theoretical conformation analysis of MCD peptide]. / Teoreticheskii konformatsionnyi analiz MCD-peptida.

The spatial structure of the MCD-peptide from bee venom has been calculated basing on the known sequence of 22 amino acid. The a priori calculations produce a system of two disulfide bonds, identical to that observed in the native structure. The calculated structure of MCD-peptide is close to that proposed earlier for the homologues peptide tertiapin and is confirmed by NMR and CD data.