Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K_d values that vary between 10^− 4 and 10^− 14 M, characteristics that are explained by a ‘dual-recognition’ mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential ΔΔGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental ΔΔG data, particularly for interactions not mediated by interfacial water molecules. ΔΔG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Amino Acid Insertion Reveals a Necessary Three-Helical Intermediate in the Folding Pathway of the Colicin E7 Immunity Protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ?-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Probing the folding intermediate of Rd-apocyt b562 by protein engineering and infrared T-jump

Small proteins often fold in an apparent two-state manner with the absence of detectable early-folding intermediates. Recently, using native-state hydrogen exchange, intermediates that exist after the rate-limiting transition state have been identified for several proteins. However, little is known about the folding kinetics from these post-transition intermediates to their corresponding native states. Herein, we have used protein engineering and a laser-induced temperature-jump (T-jump) technique to investigate this issue and have applied it to Rd-apocyt b(562) , a four-helix bundle protein. Previously, it has been shown that Rd-apocyt b(562) folds via an on-pathway hidden intermediate, which has only the N-terminal helix unfolded. In the present study, a double mutation (V16G/I17A) in the N-terminal helix of Rd-apocyt b(562) was made to further increase the relative population of this intermediate state at high temperature by selectively destabilizing the native state. In the circular dichroism thermal melting experiment, this mutant showed apparent two-state folding behavior. However, in the T-jump experiment, two kinetic phases were observed. Therefore, these results are in agreement with the idea that a folding intermediate is populated on the folding pathway of Rd-apocyt b(562) . Moreover, it was found that the exponential growth rate of the native state from this intermediate state is roughly (25 microsec)(-1) at 65 degrees C.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Cooperative folding in a multi-domain protein

Most protein domains are found in multi-domain proteins, yet most studies of protein folding have concentrated on small, single-domain proteins or on isolated domains from larger proteins. Spectrin domains are small (106 amino acid residues), independently folding domains consisting of three long alpha-helices. They are found in multi-domain proteins with a number of spectrin domains in tandem array. Structural studies have shown that in these arrays the last helix of one domain forms a continuous helix with the first helix of the following domain. It has been demonstrated that a number of spectrin domains are stabilised by their neighbours. Here we investigate the molecular basis for cooperativity between adjacent spectrin domains 16 and 17 from chicken brain alpha-spectrin (R16 and R17). We show that whereas the proteins unfold as a single cooperative unit at 25 degrees C, cooperativity is lost at higher temperatures and in the presence of stabilising salts. Mutations in the linker region also cause the cooperativity to be lost. However, the cooperativity does not rely on specific interactions in the linker region alone. Most mutations in the R17 domain cause a decrease in cooperativity, whereas proteins with mutations in the R16 domain still fold cooperatively. We propose a mechanism for this behaviour.     

Branching in the sequential folding pathway of cytochrome c

Previous results indicate that the folding pathways of cytochrome c and other proteins progressively build the target native protein in a predetermined stepwise manner by the sequential formation and association of native-like foldon units. The present work used native state hydrogen exchange methods to investigate a structural anomaly in cytochrome c results that suggested the concerted folding of two segments that have little structural relationship in the native protein. The results show that the two segments, an 18-residue omega loop and a 10-residue helix, are able to unfold and refold independently, which allows a branch point in the folding pathway. The pathway that emerges assembles native-like foldon units in a linear sequential manner when prior native-like structure can template a single subsequent foldon, and optional pathway branching is seen when prior structure is able to support the folding of two different foldons.     

Microsecond hydrophobic collapse in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.     

Topology and sequence in the folding of a TIM barrel protein: global analysis highlights partitioning between transient off-pathway and stable on-pathway folding intermediates in the complex folding mechanism of a (betaalpha)8 barrel of unknown function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.     

The stability of transmembrane helix interactions measured in a biological membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.     

The Determinants of Stability and Folding in Evolutionarily Diverged Cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by F?rster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.     

Model of the outer membrane potential generation by the inner membrane of mitochondria.

Voltage-dependent anion channels in the outer mitochondrial membrane are strongly regulated by electrical potential. In this work, one of the possible mechanisms of the outer membrane potential generation is proposed. We suggest that the inner membrane potential may be divided on two resistances in series, the resistance of the contact sites between the inner and outer membranes and the resistance of the voltage-dependent anion channels localized beyond the contacts in the outer membrane. The main principle of the proposed mechanism is illustrated by simplified electric and kinetic models. Computational behavior of the kinetic model shows a restriction of the steady-state metabolite flux through the mitochondrial membranes at relatively high concentration of the external ADP. The flux restriction was caused by a decrease of the voltage across the contact sites and by an increase in the outer membrane potential (up to +60 mV) leading to the closure of the voltage-dependent anion channels localized beyond the contact sites. This mechanism suggests that the outer membrane potential may arrest ATP release through the outer membrane beyond the contact sites, thus tightly coordinating mitochondrial metabolism and aerobic glycolysis in tumor and normal proliferating cells.     

Purification,crystal structure determination and functional characterization of type III antifreeze proteins from the European eelpout Zoarces viviparus

Antifreeze proteins (AFPs) are essential components of many organisms adaptation to cold temperatures. Fish type III AFPs are divided into two groups, SP isoforms being much less active than QAE1 isoforms. Two type III AFPs from Zoarces viviparus, a QAE1 (ZvAFP13) and an SP (ZvAFP6) isoform, are here characterized and their crystal structures determined. We conclude that the higher activity of the QAE1 isoforms cannot be attributed to single residues, but rather a combination of structural effects. Furthermore both ZvAFP6 and ZvAFP13 crystal structures have water molecules around T18 equivalent to the tetrahedral-like waters previously identified in a neutron crystal structure. Interestingly, ZvAFP6 forms dimers in the crystal, with a significant dimer interface. The presence of ZvAFP6 dimers was confirmed in solution by native electrophoresis and gel filtration. To our knowledge this is the first report of dimerization of AFP type III proteins.