Studies on Hydrogen Bonds Part V-Hydrogen Bonding in Energy Minimization Studies of Peptides

Abstract

An energy term, representing the N—H…O type of hydrogen bond, which is a function of the hydrogen bond length (R) and angle (θ) has been introduced in an energy minimization program, taking into consideration its interpolation with the non-bonded energy for borderline values of R and θ. The details of the mathematical formulation of the derivatives of the hydrogen bond function as applicable to the energy minimization have been given. The minimization technique has been applied to hydrogen bonded two and three linked peptide units (γ-turns and β-turns), and having Gly, Ala and Pro side chains. Some of the conformational highlights of the resulting minimum energy conformations are a) the occurrence of the expected 4?1 hydrogen bond in all of the β-turn tripeptide sequences and b) the presence of an additional 3?1 hydrogen bond in some of the type I and II tripeptides with the hydrogen bonding scheme in such type I β-turns occurring in a bifurcated form. These and other conformational features have been discussed in the light of experimental evidence and theoretical predictions of other workers.     

Hydrogen-bonded turns in proteins: the case for a recount

Beta-turns are sites at which proteins change their overall chain direction, and they occur with high frequency in globular proteins. The Protein Data Bank has many instances of conformations that resemble beta-turns but lack the characteristic N-H(i) --> O=C(i - 3) hydrogen bond of an authentic beta-turn. Here, we identify potential hydrogen-bonded beta-turns in the coil library, a Web-accessible database utility comprised of all residues not in repetitive secondary structure, neither alpha-helix nor beta-sheet (http://www.roselab.jhu.edu/coil). In particular, candidate turns were identified as four-residue segments satisfying highly relaxed geometric criteria but lacking a strictly defined hydrogen bond. Such candidates were then subjected to a minimization protocol to determine whether slight changes in torsion angles are sufficient to shift the conformation into reference-quality geometry without deviating significantly from the original structure. This approach of applying constrained minimization to known structures reveals a substantial population of previously unidentified, stringently defined, hydrogen-bonded beta-turns. In particular, 33% of coil library residues were classified as beta-turns prior to minimization. After minimization, 45% of such residues could be classified as beta-turns, with another 8% in 3(10) helixes (which closely resemble type III beta-turns). Of the remaining coil library residues, 37% have backbone dihedral angles in left-handed polyproline II structure.     

Study of beta-turns in globular proteins

The formation of beta-turns in globular proteins has been studied by the method of molecular mechanics. Statistical method of discriminant analysis was applied to calculate energy components and sequences of oligopeptide segments, and after this prediction of I type beta-turns has been drawn. The accuracy of true positive prediction is 65%. Components of conformational energy considerably affecting beta-turn formation were delineated. There are torsional energy, energy of hydrogen bonds, and van der Waals energy.     

Analysis and prediction of the different types of beta-turn in proteins

beta-Turns have been extracted from 59 non-identical proteins (resolution 2 A) using the standard criterion that the distance between C alpha (i) and C alpha (i + 3) is less than 7 A (1 A = 0.1 nm). The beta-turns have been classified, using phi, psi angles, into seven conventional turn types (I, I', II, II', IV, VIa, VIb) and a new class of beta-turn, designated type VIII, in which the central residues (i + 1, i + 2) adopt an alpha R beta conformation. Most beta-turn types are found in various topological environments, with the exception of I' and II' beta-turns, where 83% and 50%, respectively, are found in beta-hairpins. Sufficient data have been gathered to enable, for the first time, the separate statistical analysis of type I and II beta-turns. The two turn types have been shown to be strikingly different in their sequence preferences. Type I turns favour Asp, Asn, Ser and Cys at i; Asp, Ser, Thr and Pro at i + 1; Asp, Ser, Asn and Arg at i + 2; Gly, Trp and Met at i + 3, whilst type II turns prefer Pro at i + 1; Gly and Asn at i + 2; Gln and Arg at i + 3. These preferences have been explained by the specific side-chain interactions observed within the X-ray structures. The positional trends for type I and II beta-turns have been incorporated into the simple empirical predictive algorithm originally developed by P.N. Lewis et al. The program has improved the positional prediction of beta-turns, and has enhanced and extended the method by predicting the type of beta-turn. Since the observed preferences reflect local interactions these predictions are applicable not only to proteins, but also to peptides, many of which are thought to contain beta-turns.     

Mimicry by asx- and ST-turns of the four main types of beta-turn in proteins

Hydrogen-bonded beta-turns in proteins occur in four categories: type I (the most common), type II, type II', and type I'. Asx-turns resemble beta-turns, in that both have an NH. . .OC hydrogen bond forming a ring of 10 atoms. Serine and threonine side chains also commonly form hydrogen-bonded turns, here called ST-turns. Asx-turns and ST-turns can be categorized into four classes, based on side chain rotamers and the conformation of the central turn residue, which are geometrically equivalent to the four types of beta-turns. We propose asx- and ST-turns be named using the type I, II, I', and II' beta-turn nomenclature. Using this, the frequency of occurrence of both asx- and ST-turns is: type II' > type I > type II > type I', whereas for beta-turns it is type I > type II > type I' > type II'. Almost all type II asx-turns occur as a recently described three residue feature named an asx-nest.     

Conformational analysis of cyclolinopeptide A, a cyclic nonapeptide: nuclear Overhauser effect and energy minimization studies

The conformation of cyclolinopeptide A [cyclo(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val)], a naturally occurring cyclic nonapeptide has been investigated in dimethylsulfoxide solution by 270 MHz 1H-nmr. A complete assignment of all C alpha H and NH resonances has been accomplished using two-dimensional correlated spectroscopy and nuclear Overhauser effects (NOEs). Analysis of interresidue NOEs and JHNC alpha H values permit construction of a molecular model for the cyclic peptide backbone. The crude model derived from nmr has been used as a starting point for energy minimization, which yields a refined structure largely compatible with nmr observations. The major features of the conformation of cyclolinopeptide A are a Type VI beta-turn centered at Pro(1)-Pro(2), with a cis peptide bond between these residues and a gamma-turn (C7 structure) centered at Ile(6). Two intramolecular hydrogen bonds Val(9) CO--Phe(3)NH (4----1) and Leu(5) CO--Ile(7)NH (3----1) are observed in the low-energy conformation. The limited solvent accessibility observed for the Val(9) and Leu(5) NH groups in the nmr studies are rationalized in terms of steric shielding.     

Depsipeptide analogues of elastin repeating sequences: synthesis

Depsipeptide analogues of peptide sequences can help in elucidating the role of specific hydrogen bonds in determining the conformation in peptides. The repeating pentapeptide and hexapeptide sequences of elastin have been suggested to contain a type II beta-turn with a 4----1 hydrogen bond. Depsipeptide analogues of the repeating sequences of elastin in which this 4----1 hydrogen bond cannot exist were synthesized. A fragment condensation approach was employed in which the depsipeptide ester bond was introduced early in the synthesis. This approach proved to be effective, although the increased lability of the depsipeptide ester bond resulted in side products and low yields in some reactions.     

gβ- and gγ-turns in proteins revisited: A new set of amino acid turn-type dependent positional preferences and potentials

The number of beta-turns in a representative set of 426 protein three-dimensional crystal structures selected from the recent Protein Data Bank has nearly doubled and the number of gamma-turns in a representative set of 320 proteins has increased over seven times since the previous analysis. Beta-turns (7153) and gamma-turns (911) extracted from these proteins were used to derive a revised set of type-dependent amino acid positional preferences and potentials. Compared with previous results, the preference for proline, methionine and tryptophan has increased and the preference for glutamine, valine, glutamic acid and alanine has decreased for beta-turns. Certain new amino acid preferences were observed for both turn types and individual amino acids showed turn-type dependent positional preferences. The rationale for new amino acid preferences are discussed in the light of hydrogen bonds and other interactions involving the turns. Where main-chain hydrogen bonds of the type NH(i + 3) --> CO(i) were not observed for some beta-turns, other main-chain hydrogen bonds or solvent interactions were observed that possibly stabilize such beta-turns. A number of unexpected isolated beta-turns with proline at i + 2 position were also observed. The NH(i + 2) --> CO(i) hydrogen bond was observed for almost all gamma-turns. Nearly 20% classic gamma-turns and 43% inverse gamma-turns are isolated turns.     

Structure characterization of the central repetitive domain of high molecular weight gluten proteins. II. Characterization in solution and in the dry state.

The structure of the central repetitive domain of high molecular weight HMW) wheat gluten proteins was characterized in solution and in the dry state using HMW proteins Bx6 and Bx7 and a subcloned, bacterially expressed part of the repetitive domain of HMW Dx5. Model studies of the HMW consensus peptides PGQGQQ and GYYPTSPQQ formed the basis for the data analysis (van Dijk AA et al., 1997, Protein Sci 6:637-648). In solution, the repetitive domain contained a continuous nonoverlapping series of both type I and type II II beta-turns at positions predicted from the model studies; type II beta-turns occurred at QPGQ and QQGY sequences and type I beta-turns at YPTS and SPQQ. The subcloned part of the HMW Dx5 repetitive domain sometimes migrated as two bands on SDS-PAGE; we present evidence that this may be caused by a single amino acid insertion that disturbs the regular structure of beta-turns. The type I beta-turns are lost when the protein is dried on a solid surface, probably by conversion to type II beta-turns. The homogeneous type II beta-turn distribution is compatible with the formation of a beta-spiral structure, which provides the protein with elastic properties. The beta-turns and thus the beta-spiral are stabilized by hydrogen bonds within and between turns. Reformation of this hydrogen bonding network after, e.g., mechanical disruption may be important for the elastic properties of gluten proteins.     

Secondary structure of the Arg-Gly-Asp recognition site in proteins involved in cell-surface adhesion. Evidence for the occurrence of nested beta-bends in the model hexapeptide GRGDSP

The primary sequence Arg-Gly-Asp has been found in a number of proteins which bind to cell surface receptors. Studies with synthetic peptides have shown that the presence of charged side chains alone is not sufficient to confer binding activity. Application of folding algorithms to proteins and peptides having similar sequences indicates that binding activity is strongly correlated with the presence of two or more closely spaced residues that each have a high probability of initiating a beta-bend. Circular dichroic studies on the hexapeptide GRGDSP, whose sequence is contained in fibronectin and which also shows binding activity, demonstrate that it adopts an unusual conformation in aqueous solution. 1H-NMR spectra of the peptide in aqueous solution show that the two amide hydrogens of Asp4 and Ser5 exchange very slowly. Computer-assisted modeling using restrained molecular dynamics and energy minimization results in conformations that include two beta-bends of type III-III or III-I (hydrogen bonds 4----1 and 5----2), fully consistent with constraints imposed by 1H- and 13C-NMR data. It is suggested that this unusual secondary structure provides an additional specificity determinant.     

A study of core domains, and the core domain-domain interaction of cytochrome c fragment complex.

To gain insight into the folding mechanism of the cytochrome c complex, we prepared a complete set of homologous and hybrid two-fragment ferric complexes of four different types and related complexes from horse, tuna, yeast iso-l, and Candida cytochromes c. The complexes were characterized for structural properties. Apparent equilibrium constants of the complexes were determined to calculate delta G0 for binding. The results have allowed us to assign four core domains of the complex. A core domain is a structural region containing a hydrophobic core and the surrounding shell which folds and unfolds as a unit. Core domain 1 folds by itself and consists essentially of the right channel structure, found by R. E. Dickerson and colleagues, and a part of the heme. Core domains 2, 3, and 4, respectively, are assigned based on the cores located on the left (the Fe-S bond) and right sides and at the bottom of heme. Evidence of the core domain-domain interaction to stabilize the Fe-S bond, combined with the kinetic studies by G. R. Parr and H. Taniuchi, has led to a model of two alternative folding orders of the core domains for the horse type I complex: domain 1----3----2----4 or 1----2----3----4. Furthermore, delta G0 variation between the complexes has shown non-additive behavior, indicating the existence of a residue-residue interaction between the heme- and apofragments in the complex. Evidence suggests that this interaction in most cases occurs within or through the core groups of the ordered interface between the heme- and the apo-fragments formed by folding of core domains 1, 2, and 3. Evidence also suggests that such core group interaction manifests itself in the interaction to stabilize the Fe-S bond and may be manifested in the core domain-domain interaction.     

Nuclear magnetic resonance spectroscopic and computer-stimulated structural analyses of a heptapeptide sequence found around the N-glycosylation site of a proline-rich glycoprotein from human parotid saliva.

The proline-rich glycoprotein from human parotid saliva has a common heptapeptide sequence around four of six N-glycosylation sites (Maeda, N., H. S. Kim, E. A. Azen, and O. J. Smithies, 1985, J. Biol. Chem., 20:11123-11130). A synthetic model of the heptamer protein sequence, NH2-Q(1)-G(2)-G(3)-N(4)-Q(5)-S(6)-Q(7)-CONH2, was examined by nuclear magnetic resonance (NMR) spectroscopy and the ECEPP/2-VAO4A (Empirical Conformation Energy Program for Peptides) energy minimization computer algorithm (Scheraga, H. A., 1982, Quantum Chemistry Program Exchange, 454; Powell, M. J. D., 1964, Quantum Chemistry Program Exchange, 60). The NMR spectrum was almost completely assigned in dimethylsulfoxide-d6 (DMSO), and the amide chemical shift temperature dependence, phi dihedral angles, and chi 1 rotamer populations elucidated. These data indicated that a significant population of the heptamer could exist as a type I beta-turn [4----1 between Q(5) and G(2)] and/or a type II' beta-turn [4----1 between (Q)5 and G(2) and/or a gamma-turn [3----1 between Q(5) and G(3)] with the amino acid chi 1 torsion angles weighted toward the gauche- conformation. Starting from these three possible conformations, the ECEPP/2-VAO4A rigid geometry energy minimization program was used to find the localized predominant in vacuo structures of this heptapeptide sequence. The type II' beta-turn conformation best fits the data based on internuclear hydrogen-bonding distances, minimum potential energy considerations, and the NMR parameters.     

Binding energetics of phosphorus-containing inhibitors of thermolysin

The importance of a specific hydrogen bond between thermolysin and a phosphonamidate inhibitor, Z-NHCH2-PO(O-)-Leu-Leu (1) [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington D.C.) 235, 569-571], has been reevaluated. We have determined the inhibition constants (binding free energies) for thermolysin of phosphonamidate n-hexyl-P(O)(O-)-Leu-Trp-NHMe (4), phosphonate n-hexyl-P-(O)(O-)OCH(iBu)CO-Trp-NHMe (5), and phosphinates n-hexyl-P(O)(O-)CH2CH(iBu)CO-Trp-NHMe (6) and Z-NHCH2PO(O-)CH2CH(iBu)CO-Leu (3). Replacement of the P-NH group by P-CH2 (1----3 and 4----6) weakens the overall binding free energy by about 1.5 kcal/mol. A negligible difference in solvation energy has been measured for these pairs, and the basicity of the P-O- ligand for zinc in each pair remains nearly unchanged as determined by pH titration of their 31P NMR resonances. Therefore, this value of 1.5 kcal/mol can be assigned to the specific hydrogen bond known to exist between the P-NH of 1 and thermolysin [Tronrud, D. E., Holden, H. M., & Matthews, B. W. (1987) Science (Washington, D.C.) 235, 871-574] and inferred to exist between 4 and the enzyme. Substitution of P-O for P-NH (1----2 [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington, D.C.) 235, 569-571] and 4----5) weakens the overall binding free energy by 4.1 kcal/mol for each pair as the basicity of the P-O- ligand decreases by about 1.6 pH units. The measured solvation energy difference between 4 and 5 (and by inference between 1 and 2) is negligible.(ABSTRACT TRUNCATED AT 250 WORDS)     

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.     

Structure of benzyloxycarbonyl-L-alanyl-alpha-aminoisobutyl-alpha-aminoisobutylic acid

The x-ray diffraction study of the C-terminally unprotected tripeptide benzyloxy-carbonyl-L-alanyl-alpha-aminoisobutyl-alpha- aminoisobutylic acid (Z-L-Ala-Aib-Aib-OH) has shown that the molecule adopts a consecutive type III beta-turn, which characterizes a right-handed 3(10) helix. A very weak 4----1 intramolecular hydrogen bond with the long N...O distance of 3.32 A, and a unique "oxy analogue" of the 4----1 hydrogen bond wih the O...O distance of 2.77 A, were observed. The stability of the observed conformation with the asymmetric Aib residues was theoretically evaluated by a conformation-energy calculation. The stereochemical characteristics of Aib and Ala residues were made clear by a comparison of the conformations of the short peptides containing both Aib and Ala residues.     

Expanded turn conformations: characterization and sequence-structure correspondence in alpha-turns with implications in helix folding

Like the beta-turns, which are characterized by a limiting distance between residues two positions apart (i, i+3), a distance criterion (involving residues at positions i and i+4) is used here to identify alpha-turns from a database of known protein structures. At least 15 classes of alpha-turns have been enumerated based on the location in the phi,psi space of the three central residues (i+1 to i+3)-one of the major being the class AAA, where the residues occupy the conventional helical backbone torsion angles. However, moving towards the C-terminal end of the turn, there is a shift in the phi,psi angles towards more negative phi, such that the electrostatic repulsion between two consecutive carbonyl oxygen atoms is reduced. Except for the last position (i+4), there is not much similarity in residue composition at different positions of hydrogen and non-hydrogen bonded AAA turns. The presence or absence of Pro at i+1 position of alpha- and beta-turns has a bearing on whether the turn is hydrogen-bonded or without a hydrogen bond. In the tertiary structure, alpha-turns are more likely to be found in beta-hairpin loops. The residue composition at the beginning of the hydrogen bonded AAA alpha-turn has similarity with type I beta-turn and N-terminal positions of helices, but the last position matches with the C-terminal capping position of helices, suggesting that the existence of a "helix cap signal" at i+4 position prevents alpha-turns from growing into helices. Our results also provide new insights into alpha-helix nucleation and folding.     

Solid state and solution conformation of Boc-L-Met-Aib-L-Phe-OMe. Beta-turn conformation of a sequence related to an active chemotactic peptide analog

The conformation of the peptide Boc-L-Met-Aib-L-Phe-OMe has been studied in the solid state and solution by X-ray diffraction and 1H n.m.r., respectively. The peptide differs only in the N-terminal protecting group from the biologically active chemotactic peptide analog formyl-L-Met-Aib-L-Phe-OMe. The molecules adopt a type-II beta-turn in the solid state with Met and Aib as the corner residues (phi Met = -51.8 degrees, psi Met = 139.5 degrees, phi Aib = 58.1 degrees, psi Aib = 37.0 degrees). A single, weak 4----1 intramolecular hydrogen bond is observed between the Boc CO and Phe NH groups (N---O 3.25 A, N-H---O 128.4 degrees). 1H n.m.r. studies, using solvent and temperature dependencies of NH chemical shifts and paramagnetic radical induced line broadening of NH resonances, suggest that the Phe NH is solvent shielded in CDCl3 and (CD3)2SO. Nuclear Overhauser effects observed between Met C alpha H and Aib NH protons provide evidence of the occurrence of Met-Aib type-II beta-turns in these solvents.     

Waardenburg syndrome (WS): the analysis of a single family with a WS1 mutation showing linkage to RFLP markers on human chromosome 2q.

Waardenburg syndrome type I (WS1; MIM 19350) is caused by a pleiotropic, autosomal dominant mutation with variable penetrance and expressivity. Of individuals with this mutation, 20%-25% are hearing impaired. A multilocus linkage analysis of RFLP data from a single WS1 family with 11 affected individuals indicates that the WS1 mutation in this family is linked to the following four marker loci located on the long arm of chromosome 2: ALPP (alkaline phosphatase, placental), FN1 (fibronectin 1), D2S3 (a unique-copy DNA segment), and COL6A3 (collagen VI, alpha 3). For the RFLP marker loci, a multilocus linkage analysis using MLINK produced a peak lod (Z) of 3.23 for the following linkage relationships and recombination fractions (theta i): (ALPP----.000----FN1)----.122----D2S3----.267----CO L6A3. A similar analysis produced a Z of 6.67 for the following linkage relationships and theta i values among the markers and WS1: (FN1----.000----WS1----.000----ALPP)----.123----D2S 3----.246----COL6A3. The data confirm the conclusion of Foy et al. that at least some WS1 mutations map to chromosome 2q.     

Conformation of the polar headgroup of sphingomyelin and its analogues

The conformation of the polar headgroup of synthetic D-erythro-stearoylsphingomyelin (1), its L-threo-isomer (2) and phosphorothioyl analogues of 1 (3 and 4) has been studied in detail by high-resolution NMR spectroscopy. In both monomeric and aggregated states the phosphocholine function of 1 adopts the synclinal conformation (alpha 5 torsional angle), in analogy with phosphatidylcholine (Hauser, H., Guyer, W., Pascher, I., Skrabal, P. and Sundell, S. (1980) Biochemistry 19, 366-373). The conformation about the C1-C2 bond (theta 1 angle) of the sphingosine backbone is predominantly -synclinal, analogously to the conformation of the crystalline galactosyl cerebroside (Pascher, I. and Sundell, S. (1977) Chem. Phys. Lipids 20, 175-191). In contrast, the L-threo-isomer displays unrestricted rotation about C1-C2 bond. The possibility of the existence of a hydrogen bond between the 3-hydroxyl function and the bridged oxygen atom of sphingosine responsible for the different conformation of 1 and 2 is discussed. The modification of the phosphate function in 1 with sulfur has no significant effect on the conformation of the resulting analogues. The conformation of all studied compounds about the C-O phosphoester bonds (alpha 1 and alpha 4 torsion angles) is mainly antiperiplanar. Similar to other double-chain phospholipids, sphingomyelin shows a preference towards the antiperiplanar conformation about the C2-C3 bond.     

Proline-containing beta-turns. IV. Crystal and solution conformations of tert.-butyloxycarbonyl-L-prolyl-D-alanine and tert.-butyloxycarbonyl-L-prolyl-D-alanyl-L-alanine

The conformations of the dipeptide t-Boc-Pro-DAla-OH and the tripeptide t-Boc-Pro-DAla-Ala-OH have been determined in the crystalline state by X-ray diffraction and in solution by CD, n.m.r. and i.r. techniques. The unit cell of the dipeptide crystal contains two independent molecules connected by intermolecular hydrogen bonds. The urethane-proline peptide bond is in the cis orientation in both the molecular forms while the peptide bond between Pro and DAla is in the trans orientation. The single dipeptide molecule exhibits a "bent" structure which approximates to a partial beta-turn. The tripeptide adopts the 4----1 hydrogen-bonded type II beta-turn with all trans peptide bonds. In solution, the CD and i.r. data on the dipeptide indicate an ordered conformation with an intramolecular hydrogen bond. N.m.r. data indicate a significant proportion of the conformer with a trans orientation at the urethane-proline peptide bond. The temperature coefficient of the amide proton of this conformer in DMSO-d6 points to a 3----1 intramolecular hydrogen bond. Taken together, the data on the dipeptide in solution indicate the presence (in addition to the cis conformer) of a C7 conformation which is absent in the crystalline state. The spectral data on the tripeptide indicate the presence of the type II beta-turn in solution in addition to the nonhydrogen-bonded conformer with the cis peptide bond between the urethane and proline residues. The relevance of these data to studies on the substrate specificity of collagen prolylhydroxylase is pointed out.