Effect of protein backbone folding on the stability of protein-ligand complexes

The role played by the degree of folding of protein backbones in explaining the binding energetics of protein-ligand interactions has been studied. We analyzed the protein/peptide interactions in the RNase-S system in which amino acids at two positions of the peptide S have been mutated. The global degree of folding of the protein S correlates in a significant way with the free energy and enthalpy of the protein-peptide interactions. A much better correlation is found with the local contribution to the degree of folding of one amino acid residue: Thr36. This residue is shown to have a destabilizing interaction with Lys41, which interacts directly with peptide S. Another system, consisting of the interactions of small organic molecules with HIV-1 protease was also studied. In this case, the global change in the degree of folding of the protease backbone does not explain the binding energetics of protein-ligand interactions. However, a significant correlation is observed between the free energy of binding and the contribution of two amino acid residues in the HVI-1 protease: Gly49 and Ile66. In general, it was observed that the changes in the degree of folding are not restricted to the binding site of the protein chain but are distributed along the whole protein backbone. This study provides a basis for further consideration of the degree of folding as a parameter for empirical structural parametrizations of the binding energetics of protein folding and binding.     

Dissecting the energetics of protein alpha-helix C-cap termination through chemical protein synthesis

The alpha-helix is a fundamental protein structural motif and is frequently terminated by a glycine residue. Explanations for the predominance of glycine at the C-cap terminal portions of alpha-helices have invoked uniquely favorable energetics of this residue in a left-handed conformation or enhanced solvation of the peptide backbone because of the absence of a side chain. Attempts to quantify the contributions of these two effects have been made previously, but the issue remains unresolved. Here we have used chemical protein synthesis to dissect the energetic basis of alpha-helix termination by comparing a series of ubiquitin variants containing an L-amino acid or the corresponding D-amino acid at the C-cap Gly35 position. D-Amino acids can adopt a left-handed conformation without energetic penalty, so the contributions of conformational strain and backbone solvation can thus be separated. Analysis of the thermodynamic data revealed that the preference for glycine at the C' position of a helix is predominantly a conformational effect.     

NMR study and comparison of the antigenic properties of a peptide recognized by two HIV-1 neutralizing antibodies

Fab–peptide complexes formed between a 15 residue peptide derived from the HIV-1 gp120 V3 loop and two of its cognate monoclonal antibodies, 5023A and 5025A, were studied using isotope-edited solution nuclear magnetic resonance (NMR) techniques. Since these antibodies neutralize HIV-1 virus with different strain specificities, this study was conducted to better understand the nature of these differences. The amide proton and nitrogen NMR resonances of specific residues were used to monitor the backbone of this peptide in these complexes. Three central residues of this peptide (‘RAF’) were found to be strongly affected by binding to both antibodies. Several other peptide residues were affected by binding to antibody 5023A but not 5025A. The antibody epitopes mapped by NMR are similar to those obtained previously via PEPSCAN at higher pH. One main difference between the PEPSCAN and NMR determined epitopes for 5023A involved two glycine residues of the peptide. By NMR, one of these glycines was more dramatically affected by antibody binding than predicted by PEPSCAN, while the other was much less so. © 1998 John Wiley & Sons, Ltd.     

Proline cis-trans isomers in calbindin D9k observed by X-ray crystallography.

In a structure of recombinant bovine calbindin D9k, determined crystallographically to 1.6 A resolution, a proline in mixed, approximately equally populated, cis and trans conformation is observed. Isomers of this kind have not been reported in structure determinations of calbindin D9k to 2.3 A resolution or in any other crystallographically determined protein structure. The cis-trans isomerization occurs at the peptide bond between Gly42 and Pro43, which is in agreement with results from two-dimensional 1H nuclear magnetic resonance spectroscopy experiments on solutions of calbindin D9k. Alternative backbone stretches have been modeled and refined by stereochemical restrained least-squares refinement for the segment Lys41 to Pro43. The final R-value was 0.188. The structural perturbations accompanying the cis-trans isomerization are found to be very localized. The largest positional differences are observed at residue Gly42, in which the alternative positions of the oxygen atom are 3.6 A apart.     

Expression of an antimicrobial peptide identified in the male reproductive system of rats

Bin1b is a β-defensins-like molecule originally isolated from the rat epididymis. Owing to its bactericidal activity, Bin1b may have therapeutic properties suitable for the treatment of sexually transmitted diseases. The amino terminus of the mature Bin1b peptide contains a conserved myristoylated Gly residue. We studied the requirement of the terminal myristoylated Gly residue in the bactericidal activity of Bin1b and found that the terminal myristoylated Gly residue is not essential for the bactericidal activity. In addition, we expressed the tandem repeats of Bin1b in Escherichia coli and found that two tandem repeats of Bin1b protein were successfully expressed. The bacterially expressed tandem Bin1b repeats may be used in a diverse array of biochemical and cell biological studies.     

Engineering an unnatural Nalpha-anchored disulfide into BPTI by total chemical synthesis: structural and functional consequences

A disulfide-engineered analogue of bovine pancreatic trypsin inhibitor (BPTI), ((N(alpha)-(CH2)2S-)Gly38)BPTI, has been prepared using a thioester-mediated auxiliary functional group chemical ligation of a N(alpha)-ethanethiol-containing peptide segment with a peptide-alphaCOSR segment. In this study, Nalpha-(ethanethiol)Gly38 replaces the native Cys38, providing the sulfhydryl group required for ligation and folding. Comparisons between ((Nalpha-(CH2)2SH)Gly38)BPTI, synthetic native BPTI and reference BPTI purchased from Sigma were made using mass spectroscopy, enzyme inhibitor association constant determination (K(a)) and 1H-nuclear magnetic resonance total correlated spectroscopy (1H-NMR TOCSY) measurements. The K(a) value for ((Nalpha-(CH2)2SH)Gly38)BPTI was approximately 20-fold lower than synthetic and reference BPTI, which was attributed to perturbations in the binding loop of the protein (near Cys14). This hypothesis was confirmed by two-dimensional (2D) 1H-NMR TOCSY experiments. The data reported here demonstrate that total chemical synthesis by auxiliary functional group chemical ligation is a practical method for the synthesis of a novel class of biologically active protein analogues containing additional functional groups linked to the protein backbone.     

Role of the conserved acidic residue Asp21 in the structure of phosphatidylinositol 3-kinase Src homology 3 domain: circular dichroism and nuclear magnetic resonance studies

To elucidate a role of the Src homology 3 (SH3)-conserved acidic residue Asp21 of the phosphatidylinositol 3-kinase (PI3K) SH3 domain, structural changes induced by the D21N mutation (Asp21 --> Asn) were examined by circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopies. In the previous study, we demonstrated that environmental alterations occurred at the side chains of Trp55 and some Tyr residues from the comparison of the near-UV CD spectra of the PI3K SH3 domain with or without a D21N mutation [Okishio, N., et al. (2000) Biopolymers 57, 208-217]. In this work, the affected Tyr residues were identified as Tyr14 and Tyr73 by the CD analysis of a series of mutants, in which every single Tyr residue was replaced by a Phe residue with or without a D21N mutation. The (1)H and (15)N resonance assignments of the PI3K SH3 domain and its D21N mutant revealed that significant chemical shift changes occurred to the aromatic side-chain protons of Trp55 and Tyr14 upon the D21N mutation. All these aromatic residues are implicated in ligand recognition. In addition, the NMR analysis showed that the backbone conformations of Lys15-Asp23, Gly54-Trp55, Asn57-Gly58, and Gly67-Pro70 were affected by the D21N mutation. Furthermore, the (15)N[(1)H] nuclear Overhauser effect values of Tyr14, Glu19, and Glu20 were remarkably changed by the mutation. These results show that the D21N mutation causes structural deformation of more than half of the ligand binding cleft of the domain and provide evidence that Asp21 plays an important role in forming a well-ordered ligand binding cleft in cooperation with the RT loop (Lys15-Glu20).     

Transformation of the mechanism of triple-helix peptide folding in the absence of a C-terminal nucleation domain and its implications for mutations in collagen disorders

Folding abnormalities of the triple helix have been demonstrated in collagen diseases such as osteogenesis imperfecta in which the mutation leads to the substitution of a single Gly in the (Gly-X-Y)n sequence pattern by a larger residue. Model peptides can be used to clarify the details of normal collagen folding and the consequences of the interruption of that folding by a Gly substitution. NMR and CD studies show that placement of a (GPO)4 nucleation domain at the N terminus rather than the C terminus of a native collagen sequence allows the formation of a stable triple helix but alters the folding mechanism. Although C- to N-terminal directional folding occurs when the nucleation domain is at the C terminus, there is no preferential folding direction when the nucleation domain is at the N terminus. The lack of zipper-like directional folding does not interfere with triple-helix formation, and when a Gly residue is replaced by Ser to model an osteogenesis imperfecta mutation, the peptide with the N-terminal (GPO)4 domain can still form a good triple helix N-terminal to the mutation site. These peptide studies raise the possibility that mutant collagen could fold in a C to N direction in a zipper-like manner up to the mutation site and that completion of the triple helix N-terminal to the mutation would involve an alternative mechanism.     

Solution NMR evidence for a cis Tyr-Ala peptide group in the structure of [Pro93Ala] bovine pancreatic ribonuclease A

Proline peptide group isomerization can result in kinetic barriers in protein folding. In particular, the cis proline peptide conformation at Tyr92-Pro93 of bovine pancreatic ribonuclease A (RNase A) has been proposed to be crucial for chain folding initiation. Mutation of this proline-93 to alanine results in an RNase A molecule, P93A, that exhibits unfolding/refolding kinetics consistent with a cis Tyr92-Ala93 peptide group conformation in the folded structure (Dodge RW, Scheraga HA, 1996, Biochemistry 35:1548-1559). Here, we describe the analysis of backbone proton resonance assignments for P93A together with nuclear Overhauser effect data that provide spectroscopic evidence for a type VI beta-bend conformation with a cis Tyr92-Ala93 peptide group in the folded structure. This is in contrast to the reported X-ray crystal structure of [Pro93Gly]-RNase A (Schultz LW, Hargraves SR, Klink TA, Raines RT, 1998, Protein Sci 7:1620-1625), in which Tyr92-Gly93 forms a type-II beta-bend with a trans peptide group conformation. While a glycine residue at position 93 accommodates a type-II bend (with a positive value of phi93), RNase A molecules with either proline or alanine residues at this position appear to require a cis peptide group with a type-VI beta-bend for proper folding. These results support the view that a cis Pro93 conformation is crucial for proper folding of wild-type RNase A.     

Osteogenesis imperfecta missense mutations in collagen: structural consequences of a glycine to alanine replacement at a highly charged site

Glycine is required as every third residue in the collagen triple helix, and a missense mutation leading to the replacement of even one Gly in the repeating (Gly-Xaa-Yaa)(n) sequence with a larger residue leads to a pathological condition. Gly to Ala missense mutations are highly underrepresented in osteogenesis imperfecta (OI) and other collagen diseases, suggesting that the smallest replacement residue, Ala, might cause the least structural perturbation and mildest clinical consequences. The relatively small number of Gly to Ala mutation sites that do lead to OI must have some unusual features, such as greater structural disruption because of local sequence environment or location at a biologically important site. Here, peptides are used to model a severe OI case in which a Gly to Ala mutation is found within a highly stabilizing Lys-Gly-Asp sequence environment. Nuclear magnetic resonance, circular dichroism, and differential scanning calorimetry studies indicate this Gly to Ala replacement leads to a substantial loss of triple-helix stability and nonequivalence of the Ala residues in the three chains such that only one of the three Ala residues is capable of forming a good backbone hydrogen bond. Examination of reported OI Gly to Ala mutations suggests their preferential location at known collagen binding sites, and we propose that structural defects caused by Ala replacements may lead to pathology when they interfere with interactions.     

Secondary structure in the solution conformation of the proteinase inhibitor IIA from bull seminal plasma by nuclear magnetic resonance

Nuclear magnetic resonance data on the protease inhibitor IIA from bull seminal plasma were used to determine the secondary structure elements in the solution conformation of the protein. The experimental data were obtained from analyses of two-dimensional 1H nuclear magnetic resonance spectra at 500 and 360 MHz and include details of inter-residue nuclear Overhauser enhancements, vicinal spin-spin coupling constants and the sequence location of slowly exchanging amide protons. Accurate measurement of coupling constants and reliable assignments of nuclear Overhauser enhancements were facilitated by the use of absorption mode two-dimensional spectroscopy and large data matrices. It is shown that the peptide backbone is extended from residues 4 to 7, followed by a poorly defined helical region from residues 8 to 13 with a marked change of direction at residue Phe10. Residues 15 to 19 are extended and there is a kink at residue Glu20. Residues 22 to 27 form the central strand of a triple-stranded antiparallel beta-sheet, of which the other two strands are residues 29 to 33 and 49 to 53. Residues 34 to 46 form a helix. The tight turn in the beta-sheet is of type I geometry, and there is a beta-bulge at residue His53.     

Conformation of the metastasis-inhibiting laminin pentapeptide

The binding of cancer cells to the basement membrane glycoprotein laminin appears to be a critical step in the metastatic process. This binding can be inhibited competitively by a specific pentapeptide sequence (Tyr-Ile-Gly-Ser-Arg) of the laminin B1 chain, and this peptide can prevent metastasis formationin vivo. However, other similar pentapeptide sequences (e.g., Tyr-Ile-Gly-Ser-Glu) have been found to be much less active in metastasis inhibition, raising the possibility that such amino acid substitutions produce structural changes responsible for altering binding to the laminin receptor. In this study, conformational energy analysis has been used to determine the three-dimensional structures of these peptides. The results indicate that the substitution of Glu for the terminal Arg produces a significant conformational change in the peptide backbone at the middle Gly residue. These results have important implications for the design of drugs that may be useful in preventing metastasis formation and tumor spread.     

A conformational analysis of leucine enkephalin as a function of pH

The conformations of Leu enkephalin in aqueous solution have been investigated as a function of pH using molecular dynamics simulations. The simulations suggest the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each) at neutral pH, but is always unfolded at low or high pH. The folded form at neutral pH possesses a 2 --> 5 hydrogen bond and a close head to tail separation. No significant intramolecular hydrogen bonding of the carbonyl oxygens was observed in either the folded or unfolded forms of the peptide. Analysis of the Gly carbonyl oxygens and terminal groups indicated that, while the conformational population distribution of Leu enkephalin did vary noticeably as a function of pH, their hydration was essentially independent of pH and in agreement with the available NMR data. Further study indicated that the unfolded state of the peptide was not random in nature and consisted of one major unfolded backbone arrangement stabilized by a persistent hydrophobic interaction between the side chains of Tyr and Leu.     

The effect of a side chain-backbone swap on protein stability.

To assess the relative importance of backbone hydrogen bonding (H-bonding) vs. side chain hydrophobicity in protein structural formation, a method called side chain-backbone swap is proposed. Such a method swaps the side chain and backbone portions of certain amino acid residues, such as Asp, Glu, Asn, Gln, Lys, and Arg. Such a swap retains the sequence of a polypeptide and preserves the identity of the backbone linkage. On the other hand, the swap disrupts backbone H-bonding geometry because of the introduction of extra methylene groups into the peptide backbone. In this project, we chose the two-stranded alpha-helical coiled-coil to implement side chain-backbone swap. A pair of 36-residue peptides was designed. The two peptides have identical sequence with four residues in each heptad repeat occupied by glutamyl residues. Each glutamic acid was incorporated either as alpha-glutamyl residue (the peptide is denoted as alpha-Glu-36) or as gamma-glutamyl residue (the peptide is denoted as gamma-Glu-36). The inter-conversion between the two peptides constitutes a side chain-backbone swap. Residues constituting the hydrophobic core of the coiled-coil, however, are left unchanged. The peptide pair was characterized by circular dichroism spectroscopy, reversed-phase liquid chromatography (RPLC), and two-dimensional nuclear magnetic resonance (NMR). The results indicate that alpha-Glu-36 is a two-stranded alpha-helical coiled-coil while gamma-Glu-36 lacks stable structural elements. It is concluded that, at least for coiled-coils where hydrophobic interactions are predominantly long-range, local backbone H-bonding is a required for structural formation, consistent with a hierarchic folding mechanism. The methodological implication of side chain-backbone swap is also discussed.     

Conformational features of a natural break in the type IV collagen Gly-X-Y repeat

Fibrillar collagens have an absolute requirement for Gly as every 3rd residue, whereas breaks in the Gly-X-Y repeating pattern are found normally in the triple helix domains of non-fibrillar collagens, such as type IV collagen in basement membranes. In this study, a model 30-mer peptide is designed to include the interruption GPOGAAVMGPOGPO found in the alpha5 chain of type IV collagen. The GAAVM peptide forms a stable triple helix, with Tm= 29 degrees C. When compared with a control peptide with Gly as every 3rd residue, the GAAVM peptide has a marked decrease in the 225 nm maximum of its CD spectrum and a 10 degrees C drop in stability. A 50% decrease in calorimetric enthalpy is observed, which may result from disruption of ordered water structure anchored by regularly placed backbone carbonyls. NMR studies on specific 15N-labeled residues within the GAAVM peptide indicate a normal triple helical structure for Gly-Pro-Hyp residues flanking the break. The sequence within the break is not disordered but shows altered hydrogen exchange rates and an abnormal Val chemical shift. It was previously reported that a peptide designed to model a similar kind of interruption in the peptide (Pro-Hyp-Gly)10, (GPOGPOPOGPO), is unable to form a stable triple helix, and replacement of GAA by GPO or VM by PO within the GAAVM break decreases the stability. Thus, rigid imino acids are unfavorable within a break, despite their favorable stabilization of the triple helix itself. These results suggest some non-random structure typical of this category of breaks in the Gly-X-Y repeat of the triple helix.     

Evidence for a folded conformation of methionine- and leucine-enkephalin in a membrane environment

Transfer of an aqueous-soluble peptide hormone or neurotransmitter such as [Met]- or [Leu]enkephalin (Tyr1-Gly2-Gly3-Phe4-Met5(Leu5)), to the lipid-rich environment of its membrane-embedded receptor protein may convert the peptide into a ("bioactive") conformation required for eliciting biological activity. We have examined by high-resolution nuclear magnetic resonance (NMR) spectroscopy the conformational parameters of free enkephalin in aqueous solution versus those of enkephalin bound to lysophosphatidylcholine micelles using two approaches: 1) exchange rates, line broadening, coupling constants, and chemical shift changes of enkephalin backbone peptide N-H protons were measured for free and membrane-bound peptide in H2O (360 MHz, pH 5.6, 20 degrees C). A selective upfield shift observed for the Met5(Leu5) N-H proton upon lipid binding was interpreted in terms of its incorporation into an intramolecular H-bond. 2) 13C chemical shift changes induced by the shift reagent praseodymium nitrate (Pr(NO3)3) were compared in the presence and absence of lipid micelles. Significant changes occurring in Gly2 carbon atoms in membrane-bound enkephalin suggested the relative proximity of this residue to the Pr3+ atom (bound to the Met5(Leu5) COOH-terminal carboxylate 4 residues away). These combined results, in conjunction with studies on the specific interactions of enkephalin substituents with the micelles (Deber, C. M., and Behnam, B. A., (1984) Proc. Natl. Acad. Sci. U. S. A. 81, 61-65) suggest that enkephalin folds into an intramolecularly H-bonded beta-turn structure (with an H-bond between Gly2 C = O and Met5 NH) in the lipid environment. Such folding could facilitate the positioning of strategic residues in vivo as the hormone diffuses toward its receptor.     

Site-specific NMR monitoring of cis-trans isomerization in the folding of the proline-rich collagen triple helix

Understanding the folding of the proline-rich collagen triple helix requires consideration of the effects of proline cis-trans isomerization and may shed light on the misfolding of collagen in connective tissue diseases. Folding was monitored in real time by heteronuclear 2D NMR spectroscopy for the (15)N labeled positions in the triple-helical peptide T1-892 [GPAGPAGPVGPAGARGPAGPOGPOGPOGPOGV]. In the equilibrium unfolded monomer form, each labeled residue showed multiple peaks with interconversion rates consistent with cis-trans isomerization of Gly-Pro and Pro-Hyp bonds. Real-time NMR studies on the folding of T1-892 showed slow decay of monomer peaks and a concomitant increase in trimer peaks. Gly25 in the C-terminal rich (Gly-Pro-Hyp)(4) domain folds first, consistent with its being a nucleation domain. Analysis of the kinetics indicates that the folding of Gly25 is biphasic and the slower step represents cis-trans isomerization of imino acids. This illustrates that nucleation is limited by cis-trans isomerization. Monitoring Gly6, Gly10, Ala12, and Gly13 monomer and trimer peaks captures the C- to N-terminal propagation of the triple helix, which is also limited by Gly-Pro cis-trans isomerization events. The zipper-like nature of the propagation process is confirmed by the slower rate of folding of Ala6 compared to Gly13, reflecting the larger number of isomerization events encountered by the more N-terminal Ala6. The cis-trans isomerization events at multiple proline residues is a complex statistical process which can be visualized by these NMR studies.     

Application of peptide-mediated ring current shifts to the study of neurophysin-peptide interactions: a partial model of the neurophysin-peptide complex

Perdeuteriated peptides were synthesized that are capable of binding to the hormone binding site of neurophysin but that differ in the position of aromatic residues. The binding of these peptides to bovine neurophysin I and its des-1-8 derivative was studied by proton nuclear magnetic resonance spectroscopy in order to identify protein residues near the binding site through the observation of differential ring current effects on assignable protein resonances. Phenylalanine in position 3 of bound peptides was shown to induce significant ring current shifts in several resonances assignable to the 1-8 sequence, including those of Leu-3 and/or Leu-5, but was without effect on Tyr-49 ring protons. The magnitude of these shifts was dependent on the identity of peptide residue 1. By contrast, the sole demonstrable direct effect of an aromatic residue in position 1 was a downfield shift in Tyr-49 ring protons. Study of peptide binding to des-1-8-neurophysin demonstrated similar conformations of native and des-1-8 complexes except for the environment of Tyr-49, confirmed the peptide-induced ring current shift assignments in native neurophysin, and indicated an effect of binding on Thr-9. These observations are integrated with other results to provide a partial model of neurophysin-peptide complexes that places the ring of Tyr-49 at a distance 5-10 A from residue 1 of bound peptide and that places both the 1-8 sequence and the protein backbone region containing Tyr-49 proximal to each other and to peptide residue 3.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular cloning of crustacean pigment dispersing hormone precursor.

The cDNA encoding the precursor of the pigment dispersing hormone (PDH) of the shore crab, Carcinus maenas, was isolated and sequenced. The precursor consists of a putative 22 amino acid signal peptide, a putative 33 residue peptide of unknown function, and the 18 amino acid mature PDH, followed by a Gly residue which serves as a possible amide donor. The deduced mature PDH amino acid sequence is identical to those of Uca pugilator and Cancer magister, previously determined by Edman degradation.     

Backbone dynamics of green fluorescent protein and the effect of histidine 148 substitution

Green fluorescent protein (GFP) and its mutants have become valuable tools in molecular biology. GFP has been regarded as a very stable and rigid protein with the beta-barrel shielding the chromophore from the solvent. Here, we report the 15N nuclear magnetic resonance (NMR) studies on the green fluorescent protein (GFPuv) and its mutant His148Gly. 15N NMR relaxation studies of GFPuv show that most of the beta-barrel of GFP is rigid on the picosecond to nanosecond time scale. For several regions, including the first alpha-helix and beta-sheets 3, 7, 8, and 10, increased hydrogen-deuterium exchange rates suggest a substantial conformational flexibility on the microsecond to millisecond time scales. Mutation of residue 148 located in beta-sheet 7 is known to have a strong impact on the fluorescence properties of GFPs. UV absorption and fluorescence spectra in combination with 1H-15N NMR spectra indicate that the His148Gly mutation not only reduces the absorption of the anionic chromophore state but also affects the conformational stability, leading to the appearance of doubled backbone amide resonances for a number of residues. This suggests the presence of two conformations in slow exchange on the NMR time scale in this mutant.