beta-sheet propensity and its correlation with parameters based on conformation

Amino acid propensities for secondary structures and its variation across protein structures using exhaustive PDB data

Amino acid propensities for protein secondary structures are vital for protein structure prediction, understanding folding, and design, and have been studied using various theoretical and experimental methods. Traditional assessments of average propensities using statistical methods have been done on relatively smaller dataset for only a few secondary structures. They also involve averaging out the environmental factors and lack insights into consistency of preferences across diverse protein structures. While a few studies have explored variations in propensities across protein structural classes and folds, exploration of such variations across protein structures remains to be carried out. In this work, we have revised the average propensities for all six different secondary structures, namely α-helix, β-strand, 310-helix, π-helix, turn and coil, analyzing the most exhaustive dataset available till date using two robust secondary structure assignment algorithms, DSSP and STRIDE. The propensities evaluated here can serve as a standard reference. Moreover, we present here, for the first time, the propensities within individual protein structures and investigated how the preferences of residues and more interestingly, of their groups formed based on their structural features, vary across different unique structures. We devised a novel approach- the minimal set analysis, based on the propensity distribution of residues, which along with the group propensities led us to the conclusion that a residue's preference for a specific secondary structure is primarily dictated by its side chain's structural features. The findings in this study provide a more insightful picture of residues propensities and can be useful in protein folding and design studies.

Structure-Based Rational and General Strategy for Stabilizing Single-Chain T-Cell Receptors to Enhance Affinity

The T-cell receptor (TCR) is a crucial molecule in cellular immunity. The single-chain T-cell receptor (scTCR) is a potential format in TCR therapeutics because it eliminates the possibility of αβ-TCR mispairing. However, its poor stability and solubility impede the in vitro study and manufacturing of therapeutic applications. In this study, some conserved structural motifs are identified in variable domains regardless of germlines and species. Theoretical analysis helps to identify those unfavored factors and leads to a general strategy for stabilizing scTCRs by substituting residues at exact IMGT positions with beneficial propensities on the consensus sequence of germlines. Several representative scTCRs are displayed to achieve stability optimization and retain comparable binding affinities with the corresponding αβ-TCRs in the range of μM to pM. These results demonstrate that our strategies for scTCR engineering are capable of providing the affinity-enhanced and specificity-retained format, which are of great value in facilitating the development of TCR-related therapeutics.

Structure Relaxation Approximation (SRA) for Elucidation of Protein Structures from Ion Mobility Measurements (II). Protein Complexes

Characterizing structures of protein complexes and their disease-related aberrations is essential to understanding molecular mechanisms of many biological processes. Electrospray ionization coupled with hybrid ion mobility/mass spectrometry (ESI-IM/MS) methods offer sufficient sensitivity, sample throughput, and dynamic range to enable systematic structural characterization of proteomes. However, because ESI-IM/MS characterizes ionized protein systems in the gas phase, it generally remains unclear to what extent the protein ions characterized by IM/MS have retained their solution structures. Here, we discuss the first application of our computational structure relaxation approximation [Bleiholder, C.; et al. J. Phys. Chem. B 2019, 123 (13), 2756-2769] to assign structures of protein complexes in the range from ∼16 to ∼60 kDa from their "native" IM/MS spectra. Our analysis shows that the computed IM/MS spectra agree with the experimental spectra within the errors of the methods. The structure relaxation approximation (SRA) indicates that native backbone contacts appear largely retained in the absence of solvent for the investigated protein complexes and charge states. Native contacts between polypeptide chains of the protein complex appear to be retained to a comparable extent as contacts within a folded polypeptide chain. Our computations also indicate that the hallmark "compaction" often observed for protein systems in native IM/MS measurements appears to be a poor indicator of the extent to which native residue-residue interactions are lost in the absence of solvent. Further, the SRA indicates that structural reorganization of the protein systems in IM/MS measurements appears driven largely by remodeling of the protein surface that increases its hydrophobic content by approximately 10%. For the systems studied here, this remodeling of the protein surface appears to occur mainly by structural reorganization of surface-associated hydrophilic amino acid residues not associated with β-strand secondary structure elements. Properties related to the internal protein structure, as assessed by void volume or packing density, appear unaffected by remodeling of the surface. Taken together, the structural reorganization of the protein surface appears to be generic in nature and to sufficiently stabilize protein structures to render them metastable on the time scale of IM/MS measurements.

Rational Remodeling of Atypical Scaffolds for the Design of Photoswitchable Cannabinoid Receptor Tools

Azobenzene-embedded photoswitchable ligands are the widely used chemical tools in photopharmacological studies. Current approaches to azobenzene introduction rely mainly on the isosteric replacement of typical azologable groups. However, atypical scaffolds may offer more opportunities for photoswitch remodeling, which are chemically in an overwhelming majority. Herein, we investigate the rational remodeling of atypical scaffolds for azobenzene introduction, as exemplified in the development of photoswitchable ligands for the cannabinoid receptor 2 (CB2). Based on the analysis of residue-type clusters surrounding the binding pocket, we conclude that among the three representative atypical arms of the CB2 antagonist, AM10257, the adamantyl arm is the most appropriate for azobenzene remodeling. The optimizing spacer length and attachment position revealed AzoLig 9 with excellent thermal bistability, decent photopharmacological switchability between its two configurations, and high subtype selectivity. This structure-guided approach gave new impetus in the extension of new chemical spaces for tool customization for increasingly diversified photo-pharmacological studies and beyond.

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

Supramolecular control of self-assembling terthiophene-peptide conjugates through the amino acid side chain

The self-assembly of oligothiophene-peptide conjugates can be directed through the systematic variation of the peptide sequence into different nanostructures, including flat spicules, nanotubes, spiral sheets, and giant, flat sheets. Furthermore, the assembly of these molecules is not controlled by steric interactions between the amino acid side chains.

A hierarchical order within protein structures underlies large separations between strands in β-sheets

Protein β-sheets often involve nonlocal interactions between parts of the polypeptide chain that are separated by hundreds of residues, raising the question of how these nonlocal contacts form. A recent study of the smallest β-sheets found that their formation was not driven by signals hidden in the primary sequence. Instead, the strands in these sheets were either local in sequence, or, when separated by large sequential distances, the intervening residues were found to fold into compact modules that anchored distant parts of the chain in close spatial proximity. Here, we examine larger β-sheets to investigate the extensibility of this principle. From an analysis of the β-sheets in a nonredundant protein dataset, we find that a highly ordered hierarchical relationship exists in the intervening structure between nonlocal β-strands. This observation is almost universal: virtually all β-sheets, no matter their complexity, appear to adopt an antiparallel model to manage the nonlocal aspects of their assembly, one where the chain, having left the vicinity of an unfinished β-sheet, retraces its steps via the same route to complete the initial sheet. Exceptions typically involve unstructured regions at chain termini. Moreover, an analysis of the residues involved in nonlocal crossstrand interactions did not produce any evidence of a signal hidden in the sequence that might direct long-range interactions. These results build on those reported for the smallest sheets, suggesting that sheet formation is either local in sequence or local in space following prior folding events that anchor disparate parts of the chain in close proximity.

Dependence of α-helical and β-sheet amino acid propensities on the overall protein fold type

Background

A large number of studies have been carried out to obtain amino acid propensities for α-helices and β-sheets. The obtained propensities for α-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the β-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects β-sheet propensity.

Results

We calculated amino acid propensities for α-helices and β-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that α-helix propensities do not differ significantly by fold, but β-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for α-helix, but such is not the case for the β-sheet propensities. We also found some fold dependence on amino acid frequency in β-strands. Folds with a high Ser, Thr and Asn content at exposed sites in β-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = −0.90) and to have flat β-sheets. At buried sites in β-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = −0.93). "All-β" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "α/β" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The α-helix propensities are similar for all folds and for exposed and buried residues. However, β-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in β-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a β-strand or association between two β sheets.

A residue in helical conformation in the native state adopts a β-strand conformation in the folding transition state despite its high and canonical Φ-value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.

STITCHER: Dynamic assembly of likely amyloid and prion β-structures from secondary structure predictions

The supersecondary structure of amyloids and prions, proteins of intense clinical and biological interest, are difficult to determine by standard experimental or computational means. In addition, significant conformational heterogeneity is known or suspected to exist in many amyloid fibrils. Previous work has demonstrated that probability-based prediction of discrete β-strand pairs can offer insight into these structures. Here, we devise a system of energetic rules that can be used to dynamically assemble these discrete β-strand pairs into complete amyloid β-structures. The STITCHER algorithm progressively 'stitches' strand-pairs into full β-sheets based on a novel free-energy model, incorporating experimentally observed amino-acid side-chain stacking contributions, entropic estimates, and steric restrictions for amyloidal parallel β-sheet construction. A dynamic program computes the top 50 structures and returns both the highest scoring structure and a consensus structure taken by polling this list for common discrete elements. Putative structural heterogeneity can be inferred from sequence regions that compose poorly. Predictions show agreement with experimental models of Alzheimer's amyloid beta peptide and the Podospora anserina Het-s prion. Predictions of the HET-s homolog HET-S also reflect experimental observations of poor amyloid formation. We put forward predicted structures for the yeast prion Sup35, suggesting N-terminal structural stability enabled by tyrosine ladders, and C-terminal heterogeneity. Predictions for the Rnq1 prion and alpha-synuclein are also given, identifying a similar mix of homogenous and heterogeneous secondary structure elements. STITCHER provides novel insight into the energetic basis of amyloid structure, provides accurate structure predictions, and can help guide future experimental studies.

Position-specific propensities of amino acids in the β-strand

Background

Despite the importance of β-strands as main building blocks in proteins, the propensity of amino acid in β-strands is not well-understood as it has been more difficult to determine experimentally compared to α-helices. Recent studies have shown that most of the amino acids have significantly high or low propensity towards both ends of β-strands. However, a comprehensive analysis of the sequence dependent amino acid propensities at positions between the ends of the β-strand has not been investigated.

Results

The propensities of the amino acids calculated from a large non-redundant database of proteins are found to be highly position-specific and vary continuously throughout the length of the β-strand. They follow an unexpected characteristic periodic pattern in inner positions with respect to the cap residues in both termini of β-strands; this periodic nature is markedly different from that of the α-helices with respect to the strength and pattern in periodicity. This periodicity is not only different for different amino acids but it also varies considerably for the amino acids belonging to the same physico-chemical group. Average hydrophobicity is also found to be periodic with respect to the positions from both termini of β-strands.

Conclusions

The results contradict the earlier perception of isotropic nature of amino acid propensities in the middle region of β-strands. These position-specific propensities should be of immense help in understanding the factors responsible for β-strand design and efficient prediction of β-strand structure in unknown proteins.

Protein beta-sheet nucleation is driven by local modular formation

Despite its central role in the protein folding process, the specific mechanism(s) behind beta-sheet formation has yet to be determined. For example, whether the nucleation of beta-sheets, often containing strands separated in sequence by many residues, is local or not remains hotly debated. Here, we investigate the initial nucleation step of beta-sheet formation by performing an analysis of the smallest beta-sheets in a non-redundant dataset on the grounds that the smallest sheets, having undergone little growth after nucleation, will be enriched for nucleating characteristics. We find that the residue propensities are similar for small and large beta-sheets as are their interstrand pairing preferences, suggesting that nucleation is not primarily driven by specific residues or interacting pairs. Instead, an examination of the structural environments of the two-stranded sheets shows that virtually all of them are contained in single, compact structural modules, or when multiple modules are present, one or both of the chain termini are involved. We, therefore, find that beta-nucleation is a local phenomenon resulting either from sequential or topological proximity. We propose that beta-nucleation is a result of two opposite factors; that is, the relative rigidity of an associated folding module that holds two stretches of coil close together in topology coupled with sufficient chain flexibility that enables the stretches of coil to bring their backbones in close proximity. Our findings lend support to the hydrophobic zipper model of protein folding (Dill, K. A., Fiebig, K. M., and Chan, H. S. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 1942-1946). Implications for protein folding are discussed.

Surface properties and conformation of Nephila clavipes spider recombinant silk proteins at the air-water interface

The dragline fiber of spiders is composed of two proteins, the major ampullate spidroins I and II (MaSpI and MaSpII). To better understand the assembly mechanism and the properties of these proteins, the adsorption behavior of the recombinant proteins of the spider Nephila clavipes produced by Nexia Biotechnologies Inc. has been studied at the air-water interface using ellipsometry, surface pressure, rheological, and infrared measurements. The results show that the adsorption is more rapid and more molecules are present at the interface for MaSpII than for MaSpI. MaSpII has thus a higher affinity for the interface than MaSpI, which is consistent with its higher aggregation propensity in water. The films formed at the interface consist of networks containing a high content of intermolecular beta-sheets as revealed by the in situ polarization modulation infrared absorption reflection spectra. The infrared results further demonstrate that, for MaSpI, the beta-sheets are formed as soon as the proteins adsorb to the interface while for MaSpII the beta-sheet formation occurs more slowly. The amount of beta-sheets is lower for MaSpII than for MaSpI, most likely due to the presence of proline residues in its sequence. Both proteins form elastic films, but they are heterogeneous for MaSpI and homogeneous for MaSpII most probably as a result of a more ordered and slower aggregation process for MaSpII. This difference in their mechanism of assembly and interfacial behaviors does not seem to arise from their overall hydrophobicity or from a specific pattern of hydrophobicity, but rather from the longer polyalanine motifs, lower glycine content, and higher proline content of MaSpII. The propensity of both spidroins to form beta-sheets, especially the polyalanine blocks, suggests the participation of both proteins in the silk's beta-sheet crystallites.

Understanding the mechanism of beta-sheet folding from a chemical and biological perspective

Perturbing the structure of the Pin1 WW domain, a 34-residue protein comprised of three beta-strands and two intervening loops has provided significant insight into the structural and energetic basis of beta-sheet folding. We will review our current perspective on how structure acquisition is influenced by the sequence, which determines local conformational propensities and mediates the hydrophobic effect, hydrogen bonding, and analogous intramolecular interactions. We have utilized both traditional site-directed mutagenesis and backbone mutagenesis approaches to alter the primary structure of this beta-sheet protein. Traditional site-directed mutagenesis experiments are excellent for altering side-chain structure, whereas amide-to-ester backbone mutagenesis experiments modify backbone-backbone hydrogen bonding capacity. The transition state structure associated with the folding of the Pin1 WW domain features a partially H-bonded, near-native reverse turn secondary structure in loop 1 that has little influence on thermodynamic stability. The thermodynamic stability of the Pin1 WW domain is largely determined by the formation of a small hydrophobic core and by the formation of desolvated backbone-backbone H-bonds enveloped by this hydrophobic core. Loop 1 engineering to the consensus five-residue beta-bulge-turn found in most WW domains or a four-residue beta-turn found in most beta-hairpins accelerates folding substantially relative to the six-residue turn found in the wild type Pin1 WW domain. Furthermore, the more efficient five- and four-residue reverse turns now contribute to the stability of the three-stranded beta-sheet. These insights have allowed the design of Pin1 WW domains that fold at rates that approach the theoretical speed limit of folding.

Denatured-state energy landscapes of a protein structural database reveal the energetic determinants of a framework model for folding

Position-specific denatured-state thermodynamics were determined for a database of human proteins by use of an ensemble-based model of protein structure. The results of modeling denatured protein in this manner reveal important sequence-dependent thermodynamic properties in the denatured ensembles as well as fundamental differences between the denatured and native ensembles in overall thermodynamic character. The generality and robustness of these results were validated by performing fold-recognition experiments, whereby sequences were matched with their respective folds based on amino acid propensities for the different energetic environments in the protein, as determined through cluster analysis. Correlation analysis between structure and energetic information revealed that sequence segments destined for beta-sheet in the final native fold are energetically more predisposed to a broader repertoire of states than are sequence segments destined for alpha-helix. These results suggest that within the subensemble of mostly unstructured states, the energy landscapes are dominated by states in which parts of helices adopt structure, whereas structure formation for sequences destined for beta-strand is far less probable. These results support a framework model of folding, which suggests that, in general, the denatured state has evolutionarily evolved to avoid low-energy conformations in sequences that ultimately adopt beta-strand. Instead, the denatured state evolved so that sequence segments that ultimately adopt alpha-helix and coil will have a high intrinsic structure formation capability, thus serving as potential nucleation sites.

Optimization of linear disorder predictors yields tight association between crystallographic disorder and hydrophobicity

X-ray crystallographic protein structures often contain disordered regions that are observed as missing electron density. Diffraction data may give little or no direct evidence as to the specific nature of disordered regions. We have developed a weighted window-based disorder predictor optimized using crystallographic data. Performance of a predictor is strongly influenced by chain termini. Optimized score adjustment values for amino- and carboxy-terminal positions demonstrate a simple, monotonic relationship between disorder and residue distance from termini. This optimized disorder predictor performs similarly to DISOPRED2 on crystallographically disordered regions. Data-optimized residue disorder propensities show strong linear correlation with experimentally determined amino acid transfer energies between water and hydrogen-bonding organic solvents, which primarily reflect residue hydrophobicity (exemplified by the Nozaki-Tanford hydrophobicity scale). Disorder propensities do not correlate as well with transfer energies between water and apolar solvents, which primarily reflect a different hydropathic property: residue hydrophilicity (also reflected by the Kyte-Doolittle hydropathy scale). Our results suggest that while hydrophobic side-chain interactions are primarily involved in determining stability of the folded conformation, hydrogen bonding, and similar polar interactions are primarily involved in conformational and interaction specificity.

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.

Synthesis and structural model of an alpha(2,6)-sialyl-t glycosylated MUC1 eicosapeptide under physiological conditions

To study the effect of O-glycosylation on the conformational propensities of a peptide backbone, a 20-residue peptide (GSTAPPAHGVTSAPDTRPAP) representing the full length tandem repeat sequence of the human mucin MUC1 and its analogue glycosylated with the (2,6)-sialyl-T antigen on Thr11, were prepared and investigated by NMR and molecular modeling. The peptides contain both the GVTSAP sequence, which is an effective substrate for GalNAc transferases, and the PDTRP fragment, a known epitope recognized by several anti-MUC1 monoclonal antibodies. It has been shown that glycosylation of threonine in the GVTSAP sequence is a prerequisite for subsequent glycosylation of the serine at GVTSAP. Furthermore, carbohydrates serve as additional epitopes for MUC1 antibodies. Investigation of the solution structure of the sialyl-T glycoeicosapeptide in a H(2)O/D(2)O mixture (9:1) under physiological conditions (25 degrees C and pH 6.5) revealed that the attachment of the saccharide side-chain affects the conformational equilibrium of the peptide backbone near the glycosylated Thr11 residue. For the GVTSA region, an extended, rod-like secondary structure was found by restrained molecular dynamics simulation. The APDTR region formed a turn structure which is more flexibly organized. Taken together, the joined sequence GVTSAPDTR represents the largest structural model of MUC1 derived glycopeptides analyzed so far.

Spackling the crack: stabilizing human fibroblast growth factor-1 by targeting the N and C terminus beta-strand interactions

The beta-trefoil protein human fibroblast growth factor-1 (FGF-1) is made up of a six-stranded antiparallel beta-barrel closed off on one end by three beta-hairpins, thus exhibiting a 3-fold axis of structural symmetry. The N and C terminus beta-strands hydrogen bond to each other and their interaction is postulated from both NMR and X-ray structure data to be important in folding and stability. Specific mutations within the adjacent N and C terminus beta-strands of FGF-1 are shown to provide a substantial increase in stability. This increase is largely correlated with an increased folding rate constant, and with a smaller but significant decrease in the unfolding rate constant. A series of stabilizing mutations are subsequently combined and result in a doubling of the DeltaG value of unfolding. When taken in the context of previous studies of stabilizing mutations, the results indicate that although FGF-1 is known for generally poor thermal stability, the beta-trefoil architecture appears capable of substantial thermal stability. Targeting stabilizing mutations within the N and C terminus beta-strand interactions of a beta-barrel architecture may be a generally useful approach to increase protein stability. Such stabilized mutations of FGF-1 are shown to exhibit significant increases in effective mitogenic potency, and may prove useful as "second generation" forms of FGF-1 for application in angiogenic therapy.

Phosphorylation-Induced Conformational Changes in Short Peptides Probed by Fluorescence Resonance Energy Transfer in the 10 Å Domain

Phosphorylation-induced conformational changes in short polypeptides were probed by a fluorescence resonance energy transfer (FRET) method by employing a short-distance FRET pair (R(0) approximately 10 A) based on tryptophan as natural donor and a 2,3-diazabicyclo[2.2.2]oct-2-ene-labeled asparagine (Dbo) as synthetic acceptor. Two substrates for kinases, LeuArgArgTrpSerLeuGly-Dbo (peptide I) and TrpLysArgThrLeuArgArg-Dbo (peptide II), were investigated, with serine and threonine, respectively, as phosphorylation sites. Steady-state and time-resolved fluorescence experiments in H(2)O revealed a decrease in FRET efficiency for peptide I and an increase for peptide II; this suggested that the effective distances between donor and acceptor increased and decreased, respectively. The same trends and similar absolute variations in effective donor-acceptor distances were observed in propylene glycol, a less polar and highly viscous solvent; this suggested that the variations are due to intrinsic structural preferences. Fitting of the time-resolved decay traces according to a distribution function model (Gaussian distribution) provided the mean donor-acceptor distances, which showed an increase upon phosphorylation for peptide I (from 9.7 to 10.5 A) and a decrease for peptide II (from 10.9 to 9.3 A) in H(2)O. The broadness (half-width) of the distributions, which provides a measure of the rigidity of the peptides, remained similar upon phosphorylation of peptide I (3.0 versus 3.1 A), but decreased for peptide II (from 3.1 to 0.73 A in H(2)O); this suggests a more compact, structured conformation upon phosphorylation of the latter peptide. The elongation of the peptide backbone (by ca. 0.7 A) for peptide I is attributed to an increase in steric demand upon phosphorylation, which favors an extended conformation. The contraction (by ca. 1.4 A) and structural rigidification of peptide II is attributed to attractive Coulombic interactions and hydrogen bonding between the phosphate group and the arginine residues.

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.

Targeting the gatekeeper residue in phosphoinositide 3-kinases

A single residue in the ATP binding pocket of protein kinases-termed the gatekeeper-has been shown to control sensitivity to a wide range of small molecule inhibitors (Chem. Biol.2004, 11, 691; Chem. Biol.1999, 6, 671). Kinases that possess a small side chain at this position (Thr, Ala, or Gly) are readily targeted by structurally diverse classes of inhibitors, whereas kinases that possess a larger residue at this position are broadly resistant. Recently, lipid kinases of the phosphoinositide 3-kinase (PI3-K) family have become the focus of intense research interest as potential drug targets (Chem. Biol.2003, 10, 207; Curr. Opin. Pharmacol.2003, 3, 426). In this study, we identify the residue that corresponds structurally to the gatekeeper in PI3-Ks, and explore its importance in controlling enzyme activity and small molecule sensitivity. Isoleucine 848 of p110alpha was mutated to alanine and glycine, but the mutated kinase was found to have severely impaired enzymatic activity. A structural bioinformatic comparison of this kinase with its yeast orthologs identified second site mutations that rescued the enzymatic activity of the I848A kinase. To probe the dimensions of the gatekeeper pocket, a focused panel of analogs of the PI3-K inhibitor LY294002 was synthesized and its activity against gatekeeper mutated and wild-type p110alpha was assessed.

Dramatic acceleration of protein folding by stabilization of a nonnative backbone conformation

Through a mutagenic investigation of Gly-48, a highly conserved position in the Src homology 3 domain, we have discovered a series of amino acid substitutions that are highly destabilizing, yet dramatically accelerate protein folding, some up to 10-fold compared with the wild-type rate. The unique folding properties of these mutants allowed for accurate measurement of their folding and unfolding rates in water with no denaturant by using an NMR spin relaxation dispersion technique. A strong correlation was found between beta-sheet propensity and the folding rates of the Gly-48 mutants, even though Gly-48 lies in an unusual non-beta-strand backbone conformation in the native state. This finding indicates that the accelerated folding rates of the Gly-48 mutants are the result of stabilization of a nonnative beta-strand conformation in the transition-state structure at this position, thus providing the first, to our knowledge, experimentally elucidated example of a mechanism by which folding can occur fastest through a nonnative conformation. We also demonstrate that residues that are most stabilizing in the transition-state structure are most destabilizing in the native state, and also cause the greatest reductions in in vitro functional activity. These data indicate that the unusual native conformation of the Gly-48 position is important for function, and that evolutionary selection for function can result in a domain that folds at a rate far below the maximum possible.

Probing the conformational and dynamical effects of O-glycosylation within the immunodominant region of a MUC1 peptide tumor antigen

MUC1 mucin is a large transmembrane glycoprotein, the extracellular domain of which is formed by a repeating 20 amino acid sequence, GVTSAPDTRPAPGSTAPPAH. In normal breast epithelial cells, the extracellular domain is densely covered with highly branched complex carbohydrate structures. However, in neoplastic breast tissue, the extracellular domain is under-glycosylated, resulting in the exposure of a highly immunogenic core peptide epitope (PDTRP in bold above), as well as in the exposure of normally cryptic core Tn (GalNAc), STn (sialyl alpha2-6 GalNAc) and TF (Gal beta1-3 GalNAc) carbohydrates. Here, we report the results of 1H NMR structural studies, natural abundance 13C NMR relaxation measurements and distance-restrained MD simulations designed to probe the structural and dynamical effects of Tn-glycosylation within the PDTRP core peptide epitope. Two synthetic peptides were studied: a nine-residue MUC1 peptide of the sequence, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6-Arg7-Pro8-Ala9, and a Tn-glycosylated version of this peptide, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6(alphaGalNAc)-Arg7-Pro8-Ala9. The results of these studies show that a type I beta-turn conformation is adopted by residues PDTR within the PDTRP region of the unglycosylated MUC1 sequence. The existence of a similar beta-turn within the PDTRP core peptide epitope of the under-glycosylated cancer-associated MUC1 mucin protein might explain the immunodominance of this region in vivo, as the presence of defined secondary structure within peptide epitope regions has been correlated with increased immunogenicity in other systems. Our results have also shown that Tn glycosylation at the central threonine within the PDTRP core epitope region shifts the conformational equilibrium away from the type I beta-turn conformation and toward a more rigid and extended state. The significance of these results are discussed in relation to the possible roles that peptide epitope secondary structure and glycosylation state may play in MUC1 tumor immunogenicity.