Site-specific conformational determination in thermal unfolding studies of helical peptides using vibrational circular dichroism with isotopic substitution

Understanding the detailed mechanism of protein folding requires dynamic, site-specific stereochemical information. The short time response of vibrational spectroscopies allows evaluation of the distribution of populations in rapid equilibrium as the peptide unfolds. Spectral shifts associated with isotopic labels along with local stereochemical sensitivity of vibrational circular dichroism (VCD) allow determination of the segment sequence of unfolding. For a series of alanine-rich peptides that form α-helices in aqueous solution, we used isotopic labeling and VCD to demonstrate that the α-helix noncooperatively unwinds from the ends with increasing temperature. For these blocked peptides, the C-terminal is frayed at 5°C. Ab initio level theoretical simulations of the IR and VCD band shapes are used to analyze the spectra and to confirm the conformation of the labeled components. The VCD signals associated with the labeled residues are amplified by coupling to the nonlabeled parts of the molecule. Thus small labeled segments are detectable and stereochemically defined in moderately large peptides in this report of site-specific peptide VCD conformational analysis.

Observation via one-dimensional 13Calpha NMR of local conformational substates in thermal unfolding equilibria of a synthetic analog of the GCN4 leucine zipper.

Synthesis of a 33-residue, capped leucine zipper analogous to that in GCN4 is reported. Histidine and arginine residues are mutated to lysine to reduce the unfolding temperature. CD and ultracentrifugation studies indicate that the molecule is a two-stranded coiled coil under benign conditions. Versions of the same peptide are made with 99% 13Calpha at selected sites. One-dimensional 13C NMR spectra are assigned by inspection and used to study thermal unfolding equilibria over the entire transition from 8 to 73 degrees C. Spectra at the temperature extremes establish the approximate chemical shifts for folded and unfolded forms at each labeled site. Resonances for each amino acid appear at both locations at intermediate T, indicating that folded and unfolded forms interconvert slowly (> >2 ms) on the NMR time scale. Moreover, near room temperature, the structured form's resonance is double at several, but not all, sites, indicating at least two slowly interconverting, structured, local conformational substates. Analysis of the dynamics is possible. For example, near room temperature at the Val-9, Ala-24, and Gly-31 positions, the equilibrium constant for interconversion of the two structured forms is near unity and the time scale is > or= 10-20 ms.

Thermal unfolding of human high-density apolipoprotein A-1: implications for a lipid-free molten globular state.

Apolipoprotein A-1 (apoA-1) in complex with high-density lipoprotein is critically involved in the transport and metabolism of cholesterol and in the pathogenesis of atherosclerosis. We reexamined the thermal unfolding of lipid-free apoA-1 in low-salt solution at pH approximately 7, by using differential scanning calorimetry and circular dichroism. At protein concentrations <5 mg/ml, thermal unfolding of apoA-1 is resolved as an extended peak (25 degrees C-90 degrees C) that can be largely accounted for by a single reversible non-two-state transition with midpoint Tm 57 +/- 1 degree C, calorimetric enthalpy deltaH(Tm)= 200 +/- 20 kcal/mol (1 kcal = 4.18 kJ), van't Hoff enthalpy deltaHv(Tm) approximately 32.5 kcal/mol, and cooperativity deltaHv(Tm)/deltaH(Tm) approximately 0.16. The enthalpy deltaH(Tm) can be accounted for by melting of the alpha-helical structure that is inferred by CD to constitute approximately 60% of apoA-1 amino acids. Farand near-UV CD spectra reveal noncoincident melting of the secondary and tertiary structural elements and indicate a well-defined secondary structure but a largely melted tertiary structure for apoA-1 at approximately 37 degrees C and pH 7. This suggests a molten globular-like state for lipid-free apoA-1 under near-physiological conditions. Our results suggest that in vivo lipid binding by apoA-1 may be mediated via the molten globular apolipoprotein state in plasma.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

A thymidine triphosphate shape analog lacking Watson–Crick pairing ability is replicated with high sequence selectivity

Compound 1 (F), a nonpolar nucleoside analog that is isosteric with thymidine, has been proposed as a probe for the importance of hydrogen bonds in biological systems. Consistent with its lack of strong H-bond donors or acceptors, F is shown here by thermal denaturation studies to pair very poorly and with no significant selectivity among natural bases in DNA oligonucleotides. We report the synthesis of the 5′-triphosphate derivative of 1 and the study of its ability to be inserted into replicating DNA strands by the Klenow fragment (KF, exo− mutant) of Escherichia coli DNA polymerase I. We find that this nucleotide derivative (dFTP) is a surprisingly good substrate for KF; steady-state measurements indicate it is inserted into a template opposite adenine with efficiency (Vmax/Km) only 40-fold lower than dTTP. Moreover, it is inserted opposite A (relative to C, G, or T) with selectivity nearly as high as that observed for dTTP. Elongation of the strand past F in an F–A pair is associated with a brief pause, whereas that beyond A in the inverted A–F pair is not. Combined with data from studies with F in the template strand, the results show that KF can efficiently replicate a base pair (A–F/F–A) that is inherently very unstable, and the replication occurs with very high fidelity despite a lack of inherent base-pairing selectivity. The results suggest that hydrogen bonds may be less important in the fidelity of replication than commonly believed and that nucleotide/template shape complementarity may play a more important role than previously believed.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

A general method to design dominant negatives to B-HLHZip proteins that abolish DNA binding

We describe a method to design dominant-negative proteins (D-N) to the basic helix–loop–helix–leucine zipper (B-HLHZip) family of sequence-specific DNA binding transcription factors. The D-Ns specifically heterodimerize with the B-HLHZip dimerization domain of the transcription factors and abolish DNA binding in an equimolar competition. Thermal denaturation studies indicate that a heterodimer between a Myc B-HLHZip domain and a D-N consisting of a 12-amino acid sequence appended onto the Max dimerization domain (A-Max) is −6.3 kcal⋅mol−1 more stable than the Myc:Max heterodimer. One molar equivalent of A-Max can totally abolish the DNA binding activity of a Myc:Max heterodimer. This acidic extension also has been appended onto the dimerization domain of the B-HLHZip protein Mitf, a member of the transcription factor enhancer binding subfamily, to produce A-Mitf. The heterodimer between A-Mitf and the B-HLHZip domain of Mitf is −3.7 kcal⋅mol−1 more stable than the Mitf homodimer. Cell culture studies show that A-Mitf can inhibit Mitf-dependent transactivation both in acidic extension and in a dimerization-dependent manner. A-Max can inhibit Myc-dependent foci formation twice as well as the Max dimerization domain (HLHZip). This strategy of producing D-Ns may be applicable to other B-HLHZip or B-HLH proteins because it provides a method to inhibit the DNA binding of these transcription factors in a dimerization-specific manner.

Antigenic peptides containing large PEG loops designed to extend out of the HLA-A2 binding site form stable complexes with class I major histocompatibility complex molecules.

Recognition of peptides bound to class I major histocompatibility complex (MHC) molecules by specific receptors on T cells regulates the development and activity of the cellular immune system. We have designed and synthesized de novo cyclic peptides that incorporate PEG in the ring structure for binding to class I MHC molecules. The large PEG loops are positioned to extend out of the peptide binding site, thus creating steric effects aimed at preventing the recognition of class I MHC complexes by T-cell receptors. Peptides were synthesized and cyclized on polymer support using high molecular weight symmetrical PEG dicarboxylic acids to link the side chains of lysine residues substituted at positions 4 and 8 in the sequence of the HLA-A2-restricted human T-lymphotrophic virus type I Tax peptide. Cyclic peptides promoted the in vitro folding and assembly of HLA-A2 complexes. Thermal denaturation studies using circular dichroism spectroscopy showed that these complexes are as stable as complexes formed with antigenic peptides.

Partial molar volume, surface area, and hydration changes for equilibrium unfolding and formation of aggregation transition state: High-pressure and cosolute studies on recombinant human IFN-γ

The equilibrium dissociation of recombinant human IFN-γ was monitored as a function of pressure and sucrose concentration. The partial molar volume change for dissociation was −209 ± 13 ml/mol of dimer. The specific molar surface area change for dissociation was 12.7 ± 1.6 nm2/molecule of dimer. The first-order aggregation rate of recombinant human IFN-γ in 0.45 M guanidine hydrochloride was studied as a function of sucrose concentration and pressure. Aggregation proceeded through a transition-state species, N*. Sucrose reduced aggregation rate by shifting the equilibrium between native state (N) and N* toward the more compact N. Pressure increased aggregation rate through increased solvation of the protein, which exposes more surface area, thus shifting the equilibrium away from N toward N*. The changes in partial molar volume and specific molar surface area between the N* and N were −41 ± 9 ml/mol of dimer and 3.5 ± 0.2 nm2/molecule, respectively. Thus, the structural change required for the formation of the transition state for aggregation is small relative to the difference between N and the dissociated state. Changes in waters of hydration were estimated from both specific molar surface area and partial molar volume data. From partial molar volume data, estimates were 25 and 128 mol H2O/mol dimer for formation of the aggregation transition state and for dissociation, respectively. From surface area data, estimates were 27 and 98 mol H2O/mol dimer. Osmotic stress theory yielded values ≈4-fold larger for both transitions.

Adjustment of conformational flexibility is a key event in the thermal adaptation of proteins

3-Isopropylmalate dehydrogenase (IPMDH, E.C. 1.1.1.85) from the thermophilic bacterium Thermus thermophilus HB8 is homologous to IPMDH from the mesophilic Escherichia coli, but has an approximately 17°C higher melting temperature. Its temperature optimum is 22–25°C higher than that of the E. coli enzyme; however, it is hardly active at room temperature. The increased conformational rigidity required to stabilize the thermophilic enzyme against heat denaturation might explain its different temperature-activity profile. Hydrogen/deuterium exchange studies were performed on this thermophilic-mesophilic enzyme pair to compare their conformational flexibilities. It was found that Th. thermophilus IPMDH is significantly more rigid at room temperature than E. coli IPMDH, whereas the enzymes have nearly identical flexibilities under their respective optimal working conditions, suggesting that evolutionary adaptation tends to maintain a “corresponding state” regarding conformational flexibility. These observations confirm that conformational fluctuations necessary for catalytic function are restricted at room temperature in the thermophilic enzyme, suggesting a close relationship between conformational flexibility and enzyme function.

Solid-state synthesis and mechanical unfolding of polymers of T4 lysozyme

Recent advances in single molecule manipulation methods offer a novel approach to investigating the protein folding problem. These studies usually are done on molecules that are naturally organized as linear arrays of globular domains. To extend these techniques to study proteins that normally exist as monomers, we have developed a method of synthesizing polymers of protein molecules in the solid state. By introducing cysteines at locations where bacteriophage T4 lysozyme molecules contact each other in a crystal and taking advantage of the alignment provided by the lattice, we have obtained polymers of defined polarity up to 25 molecules long that retain enzymatic activity. These polymers then were manipulated mechanically by using a modified scanning force microscope to characterize the force-induced reversible unfolding of the individual lysozyme molecules. This approach should be general and adaptable to many other proteins with known crystal structures. For T4 lysozyme, the force required to unfold the monomers was 64 ± 16 pN at the pulling speed used. Refolding occurred within 1 sec of relaxation with an efficiency close to 100%. Analysis of the force versus extension curves suggests that the mechanical unfolding transition follows a two-state model. The unfolding forces determined in 1 M guanidine hydrochloride indicate that in these conditions the activation barrier for unfolding is reduced by 2 kcal/mol.

A general two-process model describes the hydrogen exchange behavior of RNase A in unfolding conditions.

When NMR hydrogen exchange was used previously to monitor the kinetics of RNase A unfolding, some peptide NH protons were found to show EX2 exchange (detected by base catalysis) in addition to the expected EX1 exchange, whose rate is limited by the kinetic unfolding process. In earlier work, two groups showed independently that a restricted two-process model successfully fits published hydrogen exchange rates of native RNase A in the range 0-0.7 M guanidinium chloride. We find that this model predicts properties that are very different from the observed properties of the EX2 exchange reactions of RNase A in conditions where guanidine-induced unfolding takes place. The model predicts that EX2 exchange should be too fast to measure by the technique used, whereas it is readily measurable. Possible explanations for the contradiction are considered here, and we show that removing the restriction from the earlier two-process model is sufficient to resolve the contradiction; instead of specifying that exchange caused by global unfolding occurs by the EX2 mechanism, we allow it to occur by the general mechanism, which includes both the EX1 and EX2 cases. It is logical to remove this restriction because global unfolding of RNase A is known to give rise to EX1 exchange in these unfolding conditions. Resolving the contradiction makes it possible to determine whether populated unfolding intermediates contribute to the EX2 exchange, and this question is considered elsewhere. The results and simulations indicate that moderate or high denaturant concentrations readily give rise to EX1 exchange in native proteins. Earlier studies showed that hydrogen exchange in native proteins typically occurs by the EX2 mechanism but that high temperatures or pH values above 7 may give rise to EX1 exchange. High denaturant concentrations should be added to the list of variables likely to cause EX1 exchange.

Binding studies of cationic thymidyl deoxyribonucleic guanidine to RNA homopolynucleotides.

Deoxyribonucleic guanidine is a potential antisense agent that is generated via the replacement of the negative phosphodiester linkages of DNA [--O--(PO2-)--O--] with positively-charged guanidinium (g) linkages [--NH--C(==NH2+)--NH--]. A pentameric thymidyl deoxyribonucleic guanidine molecule [d(Tg)4T-azido] has been shown to base pair specifically to poly(rA) with an unprecedented affinity. Both double and triple strands consisting of one and two equivalents of d(Tg)4T-azido paired with one equivalent of poly(rA) are indicated by thermal denaturation experiments. At an ionic strength of 0.22, the five bases of d(Tg)4T-azido are estimated to dissociate from a double helix with poly(rA) at > 100 degrees C! The effect of ionic strength on thermal denaturation is very pronounced, with stability greatest at low ionic strengths. The method of continuous variation indicates that there is an equilibrium complex with a molar ratio of d(Tg) to r(Ap) or d(Ap) of 2:1. Based on this evidence, models of the structures of d(Tg)9T-azido bound to r(Ap)9A are proposed.