Model of the alphaLbeta2 integrin I-domain/ICAM-1 DI interface suggests that subtle changes in loop orientation determine ligand specificity

The interaction of the alphaLbeta2 integrin with its cellular ligand the intercellular adhesion molecule-1 (ICAM-1) is critical for the tight binding interaction between most leukocytes and the vascular endothelium before transendothelial migration to the sites of inflammation. In this article we have modeled the alphaL subunit I-domain in its active form, which was computationally docked with the D1 domain of the ICAM-1 to probe potential protein-protein interactions. The experimentally observed key interaction between the carboxylate of Glu 34 in the ICAM-1 D1 domain and the metal ion-dependent adhesion site (MIDAS) in the open alphaL I-domain was consistently reproduced by our calculations. The calculations reveal the nature of the alphaLbeta2/ICAM-1 interaction and suggest an explanation for the increased ligand-binding affinity in the "open" versus the "closed" conformation of the alphaL I-domain. A mechanism for substrate selectivity among alphaL, alphaM, and alpha2 I-domains is suggested whereby the orientation of the loops within the I-domain is critical in mediating the interaction of the Glu 34 carboxylate of ICAM-1 D1 with the MIDAS.     

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.     

Vibrational infrared conformational studies of model peptides representing the semicrystalline domains of Bombyx mori silk fibroin

The structural organization of Bombyx mori silk fibroin was investigated by infrared (IR) spectroscopy. To this aim, (AG)15 and other model peptides of varying chain length, containing tyrosine (Y), valine (V), and serine (S) in the basic (AG)n sequence were synthesized by the solid phase method and their spectroscopic properties were determined. Both the position and the relative content of Y, V, and S residues in the (AG)n model system appeared critical in determining the preferred conformation, i.e., silk I, silk II, and unordered structures. Curve fitting analysis in the amide I range showed that the model peptides with prevailing silk II structure displayed different beta-sheet content, which was dependent on the degree of interruption of the (AG)n sequence. In this regard, the bands at about 1000 and 980 cm(-1), specifically assigned to the AG sequence of the B. mori silk fibroin chain, were identified as marker of the degree of interruption of the (AG)n sequence.A stable silk I structure was observed only when the Y residue was located near the chain terminus, while a silk I --> silk II conformational transition occurred when it was positioned in the central region of the peptide.Analysis of the second-derivative spectra in the amide I range allowed us to identify a band at 1639 cm(-1) (4 --> 1 hydrogen-bonded type II beta-turns), which is characteristic of the silk I conformation.     

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.     

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.     

Probing structural requirements of fMLP receptor: on the size of the hydrophobic pocket corresponding to residue 2 of the tripeptide.

The conformationally constrained f-L-Met-Ac(n)c-L-Phe-OMe (n = 4,9-12) tripeptides, analogues of the chemoattractant f-L-Met-L-Leu-L-Phe-OH, were synthesized in solution by classical methods and fully characterized. These compounds and the published f-L-Met-Xxx-L-Phe-OMe (Xxx = Aib and Ac(n)c where n = 3, 5-8) analogues were compared to determine the combined effect of backbone preferred conformation and side-chain bulkiness at position 2 on the relation of 3D-structure to biological activity. A conformational study of all the analogues was performed in solution by FT-IR absorption and 1H-NMR techniques. In parallel, each peptide was tested for its ability to induce chemotaxis, superoxide anion production and lysozyme secretion from human neutrophils. The biological and conformational data are discussed in relation to the proposed model of the chemotactic receptor on neutrophils, in particular of the hydrophobic pocket accommodating residue 2 of the tripeptide.     

Preferred conformations of RGDX tetrapeptides to inhibit the binding of fibrinogen to platelets

The conformational study on Arg-Gly-Asp (RGD)-containing tetrapeptides in the unhydrated and hydrated states has been carried out using the force field ECEPP/3 and the hydration shell model. The tetrapeptides studied here are H-RGDX-OH (X = Trp, Tyr, Phe, Leu, Val, Cys, Gln, and Ser), which show the inhibitory activity for binding of fibrinogen to platelets in the order of RGDW approximately equal to RGDY approximately equal to RGDF approximately equal to RGDL > RGDV > or = RGDC > or = RGDQ > or = RGDS. The backbone conformations with two C(7) backbone-to-backbone hydrogen bonds between Asp and Arg residues and between Xaa and Gly residues are in common most probable for the RGD sequence of RGDX tetrapeptides in the hydrated state. The dominant beta-turns for RGDX are found to be the types V' and IV at Gly-Asp and Asp-Xaa sequences, respectively, which are quite similar to the types II' and I (or II), respectively. However, it cannot be ruled out that the extended conformations are also remarkably feasible for RGDX tetrapeptides in water by peering the distributions of backbone conformations. These calculated results are consistent with the experimental results on RGD-containing proteins and conformationally constrained RGD-containing peptides. The reason why the RGDX becomes more potent as the side chain of the X residue is more hydrophobic may be ascribed to that the more hydrophobic is the residue X, the more populated are beta-turn structures for the Gly-Asp sequence. The hydrophobic side chain of X residue exposed to water is likely to interact with the hydrophobic region of receptor easily.     

Structural motifs in the maturation process of peptide hormones. The somatostatin precursor. I. A CD conformational study.

Synthetic peptides reproducing both the native domain around the dibasic cleavage site of prosomatostatin, and mutated sequences there of, previously assayed in site-directed mutagenesis experiments, have been studied by CD in different solvent systems, such as water, TFE/H2O, MeCN/H2O and aqueous SDS, in order to ascertain the ability of each solvent to stabilize secondary structural motifs. A combination of deconvolution methods and empirical calculations, that allow subtraction of the contributions due to unordered structures from the spectra, suggests that mainly two distinct families of ordered conformers containing alpha-helix and/or structurally different beta-turns are present in solution, the relative stability of the different conformers depending on the nature of the solvent. The presence of beta-turns is in line with a previous NMR study in DMSO and DMSO/H2O. Comparison of the CD spectra in aqueous SDS of peptides undergoing processing with a sequence not processed in vivo shows that only the latter possesses a stable and detectable alpha-helix population. This observation suggests that the structuration involving beta-turns but no alpha-helix, which was observed by CD both in SDS and organic solvent/H2O mixtures at high water contents, might be of biological significance. The similarity of this structuration to molecular models obtained from NMR data in DMSO and DMSO/H2O is discussed.     

Conformational analysis of endomorphin-2 by molecular dynamics methods

Endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH(2)) is a highly potent and selective mu-opioid receptor agonist. A conformational analysis of EM2 was carried out by simulated annealing (SA) and molecular dynamics (MD) methods. Molecular modeling was conducted on both neutral (N-terminal NH(2)) and charged (N-terminal NH(3) (+)) molecules. Based on the results of NMR investigations showing an equilibrium mixture of cis and trans Tyr(1)-Pro(2) peptide bonds for EM2 in solution, simulations were performed with restrained cis-Pro and trans-Pro peptide bonds, too. A separate SA study with unrestrained Pro peptide bonds was also conducted. Preferred conformational states are presented in Ramachandran plots. The g(+), g(-), and trans populations of the aromatic amino acid residue side chains were determined in chi(1) space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic rings and the Pro(2) ring, were determined. The energy distribution of the structures obtained by SA was calculated. The preferred secondary structural elements were different kinds of beta-turns, an inverse gamma-turn located in the N-terminal region, and regular and inverse gamma-turns located in the C-terminal region. These turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

Structural characterization of human elastin derived peptides containing the GXXP sequence

The degradation of elastin, the insoluble biopolymer of tropoelastin, can lead to the production of small peptides. These elastin-derived peptides (EDPs) are playing a key role in cellular behavior within the extracellular matrix, showing a great variety of biological effects such as chemotaxis, stimulation of cell proliferation, ion flux modifications, vasorelaxation, and inflammatory enzymes secretion. It has also been demonstrated recently that EDPs containing the GXXPG motif could induce pro-MMP1 and pro-MMP3 upregulation. Elastolysis could then cause collagen degradation and play an important role in the aging process. Many experimental studies have been devoted to EDPs, but their structure/activity relationships are not well elucidated yet. However, the assumption that their active conformation is a type VIII beta-turn on GXXP was highly suggested on the basis of predictive statistical calculations. Investigation of the EDPs three-dimensional (3D) structure would provide useful information for drug-design strategies to propose specific inhibitors. The work presented here reports theoretical results obtained from molecular dynamics simulations performed over 128 human EDPs containing the GXXP motif. We show that all the peptides, for which the central residues are not glycines, adopt a canonical (or very close to) type VIII beta-turn structure on the GXXP sequence. Amino acids surrounding this motif are also important for the structural behavior. Any residue located before the GXXP motif (XGXXP) increases the beta-turn stabilization, whereas the residue located after GXXP (GXXPX) has no significant structural effect. Moreover, we show their biological activity can be correlated with their ability to exhibit a type VIII beta-turn conformation.     

Peptides containing 4-amino-1,2-dithiolane-4-carboxylic acid (Adt): conformation of Boc-Adt-Adt-NHMe and NH...S interactions.

To study the conformational preferences induced by the insertion of the 4-amino-1,2-dithiolane-4-carboxylic acid (Adt) residue into a peptide backbone, the achiral N-protected dipeptide methylamide Boc-Adt-Adt-NHMe (1) was synthesized and its crystal state and solution conformation studied and compared with that exhibited by its carba-analogue Boc-Ac5c-Ac5c-NHMe containing two residues of 1-aminocyclopentane-1-carboxylic acid (Ac5c). Compound 1 in the crystal adopts a type-III beta-turn conformation and an analogous structure is that preferred in chloroform solution as established by 1H-NMR and NOE information. In the crystal packing three different Adt rings form a cavity and the involved sulphur atoms give rise to unusual multiple interactions with one NH group. The chemical nature of these intermolecular and intramolecular main-chain...side-chain NH...S interactions are discussed in terms of quantum chemical calculations.     

Crystal-state 3D-structural characterization of novel, Aib-based, turn and helical peptides.

The crystal-state conformations of the hexapeptide amide Pht-(Aib)(6)-NH-C(CH(3))(2)-O-OtBu (7), the hexapeptide Ac-L-aIle-(Aib)(5)-OtBu (6), the pentapeptide Z-(Aib)(3)-L-Glu(OtBu)-Aib-O-(CH(2))(2)-(1)Nap (5), the tetrapeptides Z-(Aib)(2)-L-His(N(tau)-Trt)-Aib-OMe (4 I) and Z-(Aib)(2)-L-Nva-Aib-OtBu (4 II), the tripeptide Pyr-(Aib)(3)-OtBu (3 I), the dipeptide amides Pyr-(Aib)(2)-(4)NH-TEMPO (3 II) and Piv-(Aib)(2)-NH-C(CH(3))(2)-O-OtBu (3 III), and the dipeptides Pht-Aib-betaAc(6)c-OtBu (2 I), Pht-Aib-NH-C(CH(3))(2)-O-OtBu (2 II) and Boc-gGly-mAib-OH (2 III) have been determined by X-ray diffraction analyses. All peptides investigated are characterized by one or more turn/helix forming Aib residues. Except the three short dipeptides, all are folded into C==O...H--N intramolecularly H-bonded 3(10)-helices, or into various types of beta-turns. In the structure of 6, two independent molecules of opposite screw sense were observed in the asymmetric unit, generating diastereomeric 3(10)-helices.     

Metal ion-binding ability of tetrapeptides containing alpha-aminoisobutyric acid.

alpha-Aminoisobutyric acid (Aib), one of the Calpha,alpha-disubstituted glycines, is a sterically hindered amino acid that acts as a conformational constraint in peptides. However, studies for the application of the ability of Aib to control conformation are quite few. The paper focuses on the molecular recognition ability of acyclic oligopeptides containing Aib. Liquid-liquid extraction of nine kinds of metal ions from aqueous layers to nonpolar organic layers with acyclic tetrapeptides, X-Trp-Xaa2-Gly-Xaa4-NH-Ar (X = H or C6H5CH2OCO (Z), Xaa2 = Aib or Gly, Xaa4 = Leu or Ala, Ar = phenyl or 3,5-dimethylphenyl) was examined using picrate as the anion of ion pairs. The extraction behaviour of the metal ions with the tetrapeptides was investigated in the pH range from 3 to 9. In the case of basic pH regions, Cu(II) and Ag(I) were effectively extracted with Trp-Aib-Gly-Leu-NH-Ar. Pd(II) was specifically extracted with Trp-Aib-Gly-Leu-NH-Ar in acidic pH regions. The extraction percent (%E) of the peptide host, which has a 3,5-dimethylphenyl group, was even larger than that of the host, which has a phenyl group. Moreover, Pd(II) was extracted with a peptide host which has Leu and a 3,5-dimethylphenyl group in the absence of picrate as the anion of ion pairs. The free alpha-amino group, the turn conformation and the hydrophobicity of peptide molecules were important factors for the extraction of the metals.     

Structure validation by Calpha geometry: phi,psi and Cbeta deviation

Geometrical validation around the Calpha is described, with a new Cbeta measure and updated Ramachandran plot. Deviation of the observed Cbeta atom from ideal position provides a single measure encapsulating the major structure-validation information contained in bond angle distortions. Cbeta deviation is sensitive to incompatibilities between sidechain and backbone caused by misfit conformations or inappropriate refinement restraints. A new phi,psi plot using density-dependent smoothing for 81,234 non-Gly, non-Pro, and non-prePro residues with B < 30 from 500 high-resolution proteins shows sharp boundaries at critical edges and clear delineation between large empty areas and regions that are allowed but disfavored. One such region is the gamma-turn conformation near +75 degrees,-60 degrees, counted as forbidden by common structure-validation programs; however, it occurs in well-ordered parts of good structures, it is overrepresented near functional sites, and strain is partly compensated by the gamma-turn H-bond. Favored and allowed phi,psi regions are also defined for Pro, pre-Pro, and Gly (important because Gly phi,psi angles are more permissive but less accurately determined). Details of these accurate empirical distributions are poorly predicted by previous theoretical calculations, including a region left of alpha-helix, which rates as favorable in energy yet rarely occurs. A proposed factor explaining this discrepancy is that crowding of the two-peptide NHs permits donating only a single H-bond. New calculations by Hu et al. [Proteins 2002 (this issue)] for Ala and Gly dipeptides, using mixed quantum mechanics and molecular mechanics, fit our nonrepetitive data in excellent detail. To run our geometrical evaluations on a user-uploaded file, see MOLPROBITY (http://kinemage.biochem.duke.edu) or RAMPAGE (http://www-cryst.bioc.cam.ac.uk/rampage).     

Quantitative analysis of cyclic beta-turn models.

The beta-turn is a frequently found structural unit in the conformation of globular proteins. Although the circular dichroism (CD) spectra of the alpha-helix and beta-pleated sheet are well defined, there remains some ambiguity concerning the pure component CD spectra of the different types of beta-turns. Recently, it has been reported (Hollósi, M., Kövér, K.E., Holly, S., Radics, L., & Fasman, G.D., 1987, Biopolymers 26, 1527-1572; Perczel, A., Hollósi, M., Foxman, B.M., & Fasman, G.D., 1991a, J. Am. Chem. Soc. 113, 9772-9784) that some pseudohexapeptides (e.g., the cyclo[(delta)Ava-Gly-Pro-Aaa-Gly] where Aaa = Ser, Ser(OtBu), or Gly) in many solvents adopt a conformational mixture of type I and the type II beta-turns, although the X-ray-determined conformation was an ideal type I beta-turn. In addition to these pseudohexapeptides, conformational analysis was also carried out on three pseudotetrapeptides and three pseudooctapeptides. The target of the conformation analysis reported herein was to determine whether the ring stress of the above beta-turn models has an influence on their conformational properties. Quantitative nuclear Overhauser effect (NOE) measurements yielded interproton distances. The conformational average distances so obtained were interpreted utilizing molecular dynamics (MD) simulations to yield the conformational percentages. These conformational ratios were correlated with the conformational weights obtained by quantitative CD analysis of the same compounds. The pure component CD curves of type I and type II beta-turns were also obtained, using a recently developed algorithm (Perczel, A., Tusnády, G., Hollósi, M., & Fasman, G.D., 1991b, Protein Eng. 4(6), 669-679). For the first time the results of a CD deconvolution, based on the CD spectra of 14 beta-turn models, were assigned by quantitative NOE results. The NOE experiments confirmed the ratios of the component curves found for the two major beta-turns by CD analysis. These results can now be used to enhance the conformational determination of globular proteins on the basis of their CD spectra.     

Structure and hydration of the amylopectin trisaccharide building blocks--Synthesis, NMR, and molecular dynamics

To gain insight into the molecular details and hydration of amylopectin, the five constituting trisaccharides have been chemically synthesized as their methyl alpha-glycosides. All five trisaccharides were subjected to 950 MHz NMR spectroscopy for complete assignment and nanosecond molecular dynamics trajectories were calculated to study the structure and dynamics of the trisaccharides in aqueous solution. Systematic analysis of the simulation data revealed several examples of bridging water molecules playing an important role in the stabilization of specific amylopectin conformations, which was also supported by the experimental NMR data such as interresidue NOE's and heteronuclear scalar couplings between nuclei from neighboring residues. Although alpha-maltotriose, alpha-iso-maltotriose, alpha-panose and alpha-isopanose are relatively well characterized structures, the study also includes one less characterized trisaccharide with the structure alphaGlcp(1-->4)alphaGlcp(1-->6)alphaGlcp. This trisaccharide, tentatively labelled alpha-forkose, is located at the branch point of amylopectin, forking the amylopectin into two strands that align into double-helical segments. The results show that the conformation of alpha-forkose takes a natural bend form which fits well into the structure of the double-helical segment of amylopectin. As the only trisaccharide in this study the structure of alpha-forkose is not significantly influenced by the hydration. In contrast, alpha-isopanose takes a restricted, but rather extended form due to an exceptionally strong localized water density. The two homo-linkage oligomers, alpha-maltotriose and alpha-iso-maltotriose, showed to be the most extended and the most flexible trimers, respectively, providing regular structure for crystalline domains and maximum linker flexibility for amorphous domains.     

Mechanism-based protein design: attempted "nucleation-condensation" approach to a possible minimal helix-bundle protein

In an intended mechanism-based de novo approach, a 22-mer peptide was so designed as to make it both a stereochemically nucleatable and hydrophobically condensable minimal globular protein. Framework-like nucleation of a triple-helix bundle was targeted by employing as folding nucleators composite beta-turns that could both nucleate helices and place them in close juxtaposition for possible interhelical interaction. To promote the targeted triple-helix bundle to condense as a globular protein, an amphipathic sequence pattern was adopted for possible hydrophobic interhelical interaction. A predominantly helicogenic 22-mer amphipathic peptide was thus designed, punctuating it with composite type II'-III and type II-Asx type beta-turns as the helix nucleators cum chain reversal elements. The peptide made by solid-phase synthesis was shown by NMR and CD to be a nascent and distorted triple-helix bundle in a trifluoroethanol (TFE)-water mixture, but more or less a random coil in water. A fold nucleation effect is evident in the TFE-water mixture, but apparently the hydrophobic effect cannot sustain the peptide conformational order in water. A lack of synergy between folding nucleation and hydrophobic condensation of the peptide is possible. Indeed, a mismatch between the sequential H,P pattern of the peptide and its nascent-type globular fold in a TFE-water mixture is evident based on a simulated annealing study guided by NMR.     

Solution structure of a D,L-alternating oligonorleucine as a model of double-stranded antiparallel beta-helix

Conformational characteristics of alternating D,L linear peptides are of particular interest because of their capacity to form transmembrane channels with different transport properties, as some natural antibiotics do. Single- and double-stranded beta-helical structures are common for alternating D,L peptides. The stability of the beta-helix depends on several structural factors, such as the backbone peptide length, type and position of side chains, and nature of terminal groups. The NMR and molecular dynamics solution conformation of a synthetic alternating D,L-oligopeptide with 15 norleucines (XVMe) has been used as a model to get insight in to the conformational features of double-stranded beta-helix structures. The NH chemical shift values (delta(NH)) and long-range nuclear Overhauser effects (NOE) cross peaks, in particular interstrand connectivities, clearly point to an antiparallel double-stranded beta-helix for the XVMe major conformation in solution. An extensive set of distances (from NOE cross peaks) and H-bonds (from delta(NH)) has been included in the molecular dynamics calculations. The experimental NMR data and theoretical calculations clearly indicate that the most probable conformation of XVMe in solution is a double-strand antiparallel beta(5.6) increasing decreasing-helix structure.