Mapping tertiary interactions in protein folding reactions: a novel mass spectrometry- and chemical synthesis-based approach

A novel mass spectrometry- and chemical synthesis-based approach for studying protein folding reactions is described, and its initial application to study the folding/unfolding reaction of a homo-hexameric enzyme 4-oxalocrotonate (4OT) is reported. This new approach involves the application of total chemical synthesis to prepare protein analogues that contain a photoreactive amino acid site-specifically incorporated into their primary amino acid sequence. To this end, a photoreactive amino acid-containing analogue of 4OT in which Pro-1 was replaced with p-benzoyl-l-phenylalanine (Bpa) was prepared. This analogue can be used to map structurally specific protein-protein interactions in 4OT’s native folded state. These photocrosslinking studies and peptide mapping results with (P1Bpa)4OT indicate that this construct is potentially useful for probing the structural properties of equilibrium and kinetic intermediates in 4OT’s folding reaction.     

A biosynthetic route to photoclick chemistry on proteins

Light-induced chemical reactions exist in nature, regulating many important cellular and organismal functions, e.g., photosensing in prokaryotes and vision formation in mammals. Here, we report the genetic incorporation of a photoreactive unnatural amino acid, p-(2-tetrazole)phenylalanine (p-Tpa), into myoglobin site-specifically in E. coli by evolving an orthogonal tRNA/aminoacyl-tRNA synthetase pair and the use of p-Tpa as a bioorthogonal chemical "handle" for fluorescent labeling of p-Tpa-encoded myoglobin via the photoclick reaction. Moreover, we elucidated the structural basis for the biosynthetic incorporation of p-Tpa into proteins by solving the X-ray structure of p-Tpa-specific aminoacyl-tRNA synthetase in complex with p-Tpa. The genetic encoding of this photoreactive amino acid should make it possible in the future to photoregulate protein function in living systems.     

A synthesis strategy yielding skeletally diverse small molecules combinatorially

The efficient synthesis of small molecules having many molecular skeletons is an unsolved problem in diversity-oriented synthesis (DOS). We describe the development and application of a synthesis strategy that uses common reaction conditions to transform a collection of similar substrates into a collection of products having distinct molecular skeletons. The substrates have different appendages that pre-encode skeletal information, called sigma-elements. This approach is analogous to the natural process of protein folding in which different primary sequences of amino acids are transformed into macromolecules having distinct three-dimensional structures under common folding conditions. Like sigma-elements, the amino acid sequences pre-encode structural information. An advantage of using folding processes to generate skeletal diversity in DOS is that skeletal information can be pre-encoded into substrates in a combinatorial fashion, similar to the way protein structural information is pre-encoded combinatorially in polypeptide sequences, thus making it possible to generate skeletal diversity in an efficient manner. This efficiency was realized in the context of a fully encoded, split-pool synthesis of approximately 1260 compounds potentially representing all possible combinations of building block, stereochemical, and skeletal diversity elements.     

Photolabile and paramagnetic reagents for the investigation of transmembrane signaling events.

To investigate the dynamics of membrane processes that may be integral components of specific transmembrane signaling events we have synthesized several novel paramagnetic probes and their photoreactive counterparts. The structure of these probes was designed to (1) restrict "flipping" across the membrane bilayer; (2) contain paramagnetic or photoreactive moieties that could be placed at specific depths within the bilayer; (3) provide information about membrane structure as well as dynamics of protein movement; and (4) in the case of the photoreactive probes, be of high specific radioactivity. The molecules described in this paper consist of amino acid, dipeptide, or carbohydrate groups attached to arylazide- or nitroxide-bearing fatty acids. The synthesis and initial characterization of these membrane probes is described.     

Total chemical synthesis of crambin

Crambin is a small (46 amino acids) protein isolated from the seeds of the plant Crambe abyssinica. Crambin has been extensively used as a model protein for the development of advanced crystallography and NMR techniques and for computational folding studies. We set out to establish synthetic access to crambin. Initially, we synthesized the 46 amino acid polypeptide by native chemical ligation of two distinct sets of peptide segments (15 + 31 and 31 + 15 residues). The synthetic polypeptide chain folded in good yield to give native crambin containing three disulfide bonds. The chemically synthesized crambin was characterized by LC-MS and by 2D-NMR. However, the 31-residue peptide segments were difficult to purify, and this caused an overall low yield for the synthesis. To overcome this problem, we synthesized crambin by the native chemical ligation of three segments (15 + 16 + 15 residues). Total synthesis using the ligation of three segments gave more than a 10-fold increase in yield and a protein product of exceptionally high purity. This work demonstrates the efficacy of chemical protein synthesis by the native chemical ligation of three segments and establishes efficient synthetic access to the important model protein crambin for experimental studies of protein folding and stability.     

A statistical model for predicting protein folding rates from amino acid sequence with structural class information

Prediction of protein folding rates from amino acid sequences is one of the most important challenges in molecular biology. In this work, I have related the protein folding rates with physical-chemical, energetic and conformational properties of amino acid residues. I found that the classification of proteins into different structural classes shows an excellent correlation between amino acid properties and folding rates of two- and three-state proteins, indicating the importance of native state topology in determining the protein folding rates. I have formulated a simple linear regression model for predicting the protein folding rates from amino acid sequences along with structural class information and obtained an excellent agreement between predicted and experimentally observed folding rates of proteins; the correlation coefficients are 0.99, 0.96 and 0.95, respectively, for all-alpha, all-beta and mixed class proteins. This is the first available method, which is capable of predicting the protein folding rates just from the amino acid sequence with the aid of generic amino acid properties and structural class information.     

Predicting protein folding rates using the concept of Chou's pseudo amino acid composition

One of the most important challenges in computational and molecular biology is to understand the relationship between amino acid sequences and the folding rates of proteins. Recent works suggest that topological parameters, amino acid properties, chain length and the composition index relate well with protein folding rates, however, sequence order information has seldom been considered as a property for predicting protein folding rates. In this study, amino acid sequence order was used to derive an effective method, based on an extended version of the pseudo-amino acid composition, for predicting protein folding rates without any explicit structural information. Using the jackknife cross validation test, the method was demonstrated on the largest dataset (99 proteins) reported. The method was found to provide a good correlation between the predicted and experimental folding rates. The correlation coefficient is 0.81 (with a highly significant level) and the standard error is 2.46. The reported algorithm was found to perform better than several representative sequence-based approaches using the same dataset. The results indicate that sequence order information is an important determinant of protein folding rates.     

Exploration of two-dimensional hydrophobic-polar lattice model by combining local search with elastic net algorithm

Determining the functional conformation of a protein from its amino acid sequence remains a central problem in computational biology. In this paper, we establish the mathematical optimal model of protein folding problem (PFP) on two-dimensional space based on the minimal energy principle. A novel hybrid of elastic net algorithm and local search method (ENLS) is applied successfully to simulations of protein folding on two-dimensional hydrophobic-polar (HP) lattice model. Eight HP benchmark instances with up to 64 amino acids are tested to verify the effectiveness of proposed approach and model. In several cases, the ENLS method finds new lower energy states. The numerical results show that it is drastically superior to other methods in finding the ground state of a protein.     

Oxidative Deselenization of Selenocysteine: Applications for Programmed Ligation at Serine

Despite the unique chemical properties of selenocysteine (Sec), ligation at Sec is an under‐utilized methodology for protein synthesis. We describe herein an unprecedented protocol for the conversion of Sec to serine (Ser) in a single, high‐yielding step. When coupled with ligation at Sec, this transformation provides a new approach to programmed ligations at Ser residues. This new reaction is compatible with a wide range of functionality, including the presence of unprotected amino acid side chains and appended glycans. The utility of the methodology is demonstrated in the rapid synthesis of complex glycopeptide fragments of the epithelial glycoproteins MUC5AC and MUC4 and through the total synthesis of the structured, cysteine (Cys)‐free protein eglin C.     

Flexible Folding: Disulfide-Containing Peptides and Proteins Choose the Pathway Depending on the Environments

In the last few decades, development of novel experimental techniques, such as new types of disulfide (SS)-forming reagents and genetic and chemical technologies for synthesizing designed artificial proteins, is opening a new realm of the oxidative folding study where peptides and proteins can be folded under physiologically more relevant conditions. In this review, after a brief overview of the historical and physicochemical background of oxidative protein folding study, recently revealed folding pathways of several representative peptides and proteins are summarized, including those having two, three, or four SS bonds in the native state, as well as those with odd Cys residues or consisting of two peptide chains. Comparison of the updated pathways with those reported in the early years has revealed the flexible nature of the protein folding pathways. The significantly different pathways characterized for hen-egg white lysozyme and bovine milk α-lactalbumin, which belong to the same protein superfamily, suggest that the information of protein folding pathways, not only the native folded structure, is encoded in the amino acid sequence. The application of the flexible pathways of peptides and proteins to the engineering of folded three-dimensional structures is an interesting and important issue in the new realm of the current oxidative protein folding study.     

Fmoc-POAC: [(9-fluorenylmethyloxycarbonyl)-2,2,5,5-tetramethylpyrrolidine-N-oxyl-3-amino-4-carboxylic acid]: a novel protected spin labeled beta-amino acid for peptide and protein chemistry

The stable free radical 2,2,6,6-tetramethylpiperidine-N-oxyl-4-amino-4-carboxylic acid (TOAC) is the only spin labeled amino acid that has been used to date to successfully label peptide sequences for structural studies. However, severe difficulty in coupling the subsequent amino acid has been the most serious shortcoming of this paramagnetic marker. This problem stems from the low nucleophilicity of TOAC's amine group towards the acylation reaction during peptide chain elongation. The present report introduces the alternative beta-amino acid 2,2,5,5-tetramethylpyrrolidine-N-oxyl-3-amino-4-carboxylic acid (POAC), potentially useful in peptide and protein chemistry. Investigations aimed at addressing the stereochemistry of this cyclic molecule through X-ray diffraction measurements of crystalline and bulk samples revealed that it consists only of the trans conformer. The 9-fluorenylmethyloxycarbonyl group (Fmoc) was chosen for temporary protection of the POAC amine function, allowing insertion of the probe at any position in a peptide sequence. The vasoactive octapeptide angiotensin II (All, DRVYIHPF) was synthesized by replacing Pro7 with POAC. The reaction of Fmoc-POAC with the peptidyl-resin occurred smoothly, and the coupling of the subsequent amino acid showed a much faster reaction when compared with TOAC. POAC7-AII was obtained in good yield, demonstrating that, in addition to TOAC, POAC is a convenient amino acid for the synthesis of spin labeled peptide analogues. The present findings open the possibility of a wide range of chemical and biological applications for this novel beta-amino acid derivative, including structural investigations involving its differentiated bend-inducing characteristics.     

Synthesis of peptides and proteins without cysteine residues by native chemical ligation combined with desulfurization.

The highly chemoselective reaction between unprotected peptides bearing an N-terminal Cys residue and a C-terminal thioester enables the total and semi-synthesis of complex polypeptides. Here we extend the utility of this native chemical ligation approach to non-cysteine containing peptides. Since alanine is a common amino acid in proteins, ligation at this residue would be of great utility. To achieve this goal, a specific alanine residue in the parent protein is replaced with cysteine to facilitate synthesis by native chemical ligation. Following ligation, selective desulfurization of the resulting unprotected polypeptide product with H(2)/metal reagents converts the cysteine residue to alanine. This approach, which provides a general method to prepare alanyl proteins from their cysteinyl forms, can be used to chemically synthesize a variety of polypeptides, as demonstrated by the total chemical syntheses of the cyclic antibiotic microcin J25, the 56-amino acid streptococcal protein G B1 domain, and a variant of the 110-amino acid ribonuclease, barnase.     

Foldable subunits of helix protein

Detection of foldable subunits in proteins is an important approach to understand their evolutions and find building motifs for de novo protein design. Using united-residue model, we simulated the folding of a six-helix protein with a length of 120 amino acids (C-terminal domain of Ku86). The folding behaviors, structural topology and sequence repetition of this protein all suggest that it may have a two-fold quasi-repetition or symmetry in its sequence and structure. Therefore, we simulated the folding of its two halves (1–60 and 61–120 amino acids) and find that they can fold into native conformations independently. It is also found that their folding behaviors are very similar to other three-helix bundles. This suggests that this protein may be divided into two foldable halves.     

Selenium and Selenocysteine in Protein Chemistry

Selenocysteine, the selenium‐containing analogue of cysteine, is the twenty‐first proteinogenic amino acid. Since its discovery almost fifty years ago, it has been exploited in unnatural systems even more often than in natural systems. Selenocysteine chemistry has attracted the attention of many chemists in the field of chemical biology owing to its high reactivity and resulting potential for various applications such as chemical modification, chemical protein (semi)synthesis, and protein folding, to name a few. In this Minireview, we will focus on the chemistry of selenium and selenocysteine and their utility in protein chemistry.     

Chiral Photoresponsive Tetrathiazoles That Provide Snapshots of Folding States

Herein, we designed chiral photoresponsive tetra(2‐phenylthiazole)s, which induce a diastereoselective 6π‐electrocyclization reaction in a helically folded structure to freeze the conformational interconversions. The folding conformation with one helical turn of tetra(2‐phenylthiazole)s was supported by multiple intramolecular noncovalent interactions including vicinal S???N interheteroatom interactions and CH–π and π–π stacking interactions between nonadjacent units, as found in X‐ray crystal structures as well as quantum chemical calculations. The introduction of a chiral group at both ends of tetra(2‐phenylthiazole) dictates the preferential folding into a one‐handed helix conformation by the simultaneous operation of S???O and multiple CH–π interactions that involve the chiral end groups. Since the tetra(2‐phenylthiazole)s possess two equivalent photoreactive 6π‐electron systems and the folded conformation is suitable for photoinduced electrocyclization reaction, they undergo a photocyclization reaction in a stereoselective manner to memorize the chirality of the helix in a resulting diastereomeric closed form.     

Importance of native-state topology for determining the folding rate of two-state proteins

Understanding the relationship between amino acid sequences and folding rate of proteins is a challenging task similar to protein folding problem. In this work, we have analyzed the relative importance of protein sequence and structure for predicting the protein folding rates in terms of amino acid properties and contact distances, respectively. We found that the parameters derived with protein sequence (physical-chemical, energetic, and conformational properties of amino acid residues) show very weak correlation (|r| < 0.39) with folding rates of 28 two-state proteins, indicating that the sequence information alone is not sufficient to understand the folding rates of two-state proteins. However, the maximum positive correlation obtained for the properties, number of medium-range contacts, and alpha-helical tendency reveals the importance of local interactions to initiate protein folding. On the other hand, a remarkable correlation (r varies from -0.74 to -0.88) has been obtained between structural parameters (contact order, long-range order, and total contact distance) and protein folding rates. Further, we found that the secondary structure content and solvent accessibility play a marginal role in determining the folding rates of two-state proteins. Multiple regression analysis carried out with the combination of three properties, beta-strand tendency, enthalpy change, and total contact distance improved the correlation to 0.92 with protein folding rates. The relative importance of existing methods along with multiple-regression model proposed in this work will be discussed. Our results demonstrate that the native-state topology is the major determinant for the folding rates of two-state proteins.     

Structure prediction of an S-layer protein by the mean force method

S-layer proteins have a wide range of application potential due to their characteristic features concerning self-assembling, assembling on various surfaces, and forming of isoporous structures with functional groups located on the surface in an identical position and orientation. Although considerable knowledge has been experimentally accumulated on the structure, biochemistry, assemble characteristics, and genetics of S-layer proteins, no structural model at atomic resolution has been available so far. Therefore, neither the overall folding of the S-layer proteins-their tertiary structure-nor the exact amino acid or domain allocations in the lattices are known. In this paper, we describe the tertiary structure prediction for the S-layer protein SbsB from Geobacillus stearothermophilus PV72/p2. This calculation was based on its amino acid sequence using the mean force method (MF method) achieved by performing molecular dynamic simulations. This method includes mainly the thermodynamic aspects of protein folding as well as steric constraints of the amino acids and is therefore independent of experimental structure analysis problems resulting from biochemical properties of the S-layer proteins. Molecular dynamic simulations were performed in vacuum using the simulation software NAMD. The obtained tertiary structure of SbsB was systematically analyzed by using the mean force method, whereas the verification of the structure is based on calculating the global free energy minimum of the whole system. This corresponds to the potential of mean force, which is the thermodynamically most favorable conformation of the protein. Finally, an S-layer lattice was modeled graphically using CINEMA4D and compared with scanning force microscopy data down to a resolution of 1 nm. The results show that this approach leads to a thermodynamically favorable atomic model of the tertiary structure of the protein, which could be verified by both the MF Method and the lattice model.     

Future directions of structural mass spectrometry using hydroxyl radical footprinting

Hydroxyl radical protein footprinting coupled to mass spectrometry has been developed over the last decade and has matured to a powerful method for analyzing protein structure and dynamics. It has been successfully applied in the analysis of protein structure, protein folding, protein dynamics, and protein–protein and protein–DNA interactions. Using synchrotron radiolysis, exposure of proteins to a ‘white’ X‐ray beam for milliseconds provides sufficient oxidative modification to surface amino acid side chains, which can be easily detected and quantified by mass spectrometry. Thus, conformational changes in proteins or protein complexes can be examined using a time‐resolved approach, which would be a valuable method for the study of macromolecular dynamics. In this review, we describe a new application of hydroxyl radical protein footprinting to probe the time evolution of the calcium‐dependent conformational changes of gelsolin on the millisecond timescale. The data suggest a cooperative transition as multiple sites in different molecular subdomains have similar rates of conformational change. These findings demonstrate that time‐resolved protein footprinting is suitable for studies of protein dynamics that occur over periods ranging from milliseconds to seconds. In this review, we also show how the structural resolution and sensitivity of the technology can be improved as well. The hydroxyl radical varies in its reactivity to different side chains by over two orders of magnitude, thus oxidation of amino acid side chains of lower reactivity are more rarely observed in such experiments. Here we demonstrate that the selected reaction monitoring (SRM)‐based method can be utilized for quantification of oxidized species, improving the signal‐to‐noise ratio. This expansion of the set of oxidized residues of lower reactivity will improve the overall structural resolution of the technique. This approach is also suggested as a basis for developing hypothesis‐driven structural mass spectrometry experiments. Copyright © 2010 John Wiley & Sons, Ltd.     

Ribonuclease T1: Structure,Function, and Stability

Proteins carry out the most important and difficult tasks in all living organisms. To do so, they must often interact specifically with other small and large molecules. This requires that they fold to a globular conformation with a unique active site that is used for the specific interaction. Consequently, protein folding can be regarded as the “secret of life”. Biochemists and chemists have a great interest in elucidating the mechanism by which proteins fold and in predicting the folded conformation and its stability given just the amino acid sequence. This challenge is sometimes called the “protein folding problem”. The ability to construct proteins differing in sequence by one or more amino acids and to analyze their three-dimensional structures by X-ray crystallography and NMR spectroscopy is a powerful tool for investigating the conformational stability and folding of proteins. Several proteins are now under intensive study by this approach. One of these is ribonuclease T1.     

Transient enzyme-substrate recognition monitored by real-time NMR

Slow protein folding processes during which kinetic folding intermediates occur for an extended time can lead to aggregation and dysfunction in living cells. Therefore, protein folding helpers have evolved, which prevent proteins from aggregation and/or speed up folding processes. In this study, we present the structural characterization of a long-living transient folding intermediate of RNase T1 (S54G/P55N) by time-resolved NMR spectroscopy. NMR resonances of this kinetic folding intermediate could be assigned mainly by a real-time 3D BEST-HNCA. These assignments were the basis to investigate the interaction sites between the protein folding helper enzyme SlyD(1-165) (SlyD*) from Escherichia coli (E. coli) and this kinetic intermediate at a residue resolution. Thus, we investigated the Michaelis-Menten complex of this enzyme reaction, because the NMR data acquisition was performed during the actual catalysis. The interaction surface of the transient folding intermediate is restricted to a region around the peptidyl-prolyl bond (Y38-P39), whose isomerization is catalyzed by SlyD*. The interaction surface regarding SlyD* extends from specific amino acids of the FKBP domain forming the peptidyl-prolyl cis/trans-isomerase active site to almost the entire IF domain. This illustrates an effective interplay between the two functional domains of SlyD* to facilitate protein folding catalysis.