Characterizing molecular interactions in different bacteriorhodopsin assemblies by single-molecule force spectroscopy

Using single-molecule force spectroscopy we characterized inter- and intramolecular interactions stabilizing structural segments of individual bacteriorhodopsin (BR) molecules assembled into trimers and dimers, and monomers. While the assembly of BR did not vary the location of these structural segments, their intrinsic stability could change up to 70% increasing from monomer to dimer to trimer. Since each stable structural segment established one unfolding barrier, we conclude that the locations of unfolding barriers were determined by intramolecular interactions but that their strengths were strongly influenced by intermolecular interactions. Subtracting the unfolding forces of the BR trimer from that of monomer allowed us to calculate the contribution of inter- and intramolecular interactions to the membrane protein stabilization. Statistical analyses showed that the unfolding pathways of differently assembled BR molecules did not differ in their appearance but in their population. This suggests that in our experiments the membrane protein assembly does not necessarily change the location of unfolding barriers within the protein, but certainly their strengths, and thus alters the probability of a protein to choose certain unfolding pathways.     

Conservation of molecular interactions stabilizing bovine and mouse rhodopsin

Rhodopsin is the light receptor that initiates phototransduction in rod photoreceptor cells. The structure and function of rhodopsin are tightly linked to molecular interactions that stabilize and determine the receptor's functional state. Single-molecule force spectroscopy (SMFS) was used to localize and quantify molecular interactions that structurally stabilize bovine and mouse rhodopsin from native disk membranes of rod photoreceptor cells. The mechanical unfolding of bovine and mouse rhodopsin revealed nine major unfolding intermediates, each intermediate defining a structurally stable segment in the receptor. These stable structural segments had similar localization and occurrence in both bovine and mouse samples. For each structural segment, parameters describing their unfolding energy barrier were determined by dynamic SMFS. No major differences were observed between bovine and mouse rhodopsin, thereby implying that the structures of both rhodopsins are largely stabilized by similar molecular interactions.     

Molecular dynamics simulation of the unfolding of individual bacteriorhodopsin helices in sodium dodecyl sulfate micelles

We report molecular dynamics simulations of the trends in the changes in secondary structure of the seven individual helices of bacteriorhodopsin when inserted into sodium dodecyl sulfate (SDS) micelles, and their dependence on the amino acid sequence. The results indicate that the partitioning of the helices in the micelles and their stability are dependent on the hydrophobicity of the transmembrane segments. Helices A, B, and E are stable and retain their initial secondary structure throughout the 100 ns simulation time. In contrast, helices C, D, F, and G show structural perturbations within the first 10 ns. The instabilities are localized near charged residues within the transmembrane segments. The overall structural instability of the helix is correlated with its partitioning to the surface of the micelle and its interaction with polar groups there. The in silico experiments were performed to complement the in vitro experiments that examined the partial denaturation of bacteriorhodopsin in SDS described in the preceding article (DOI 10.1021/bi201769z ). The simulations are consistent with the trends revealed by the experimental results but strongly underestimate the extent of helix to extended coil transformation. The reason may be either that the sampling time was not sufficiently long or, more interestingly, that interhelix residue interactions play a role in the unfolding of the helices.     

Molecular mechanism of alpha 1(I)-osteogenesis imperfecta/Ehlers-Danlos syndrome: unfolding of an N-anchor domain at the N-terminal end of the type I collagen triple helix

We demonstrate that 85 N-terminal amino acids of the alpha1(I) chain participate in a highly stable folding domain, acting as the stabilizing anchor for the amino end of the type I collagen triple helix. This anchor region is bordered by a microunfolding region, 15 amino acids in each chain, which include no proline or hydroxyproline residues and contain a chymotrypsin cleavage site. Glycine substitutions and amino acid deletions within the N-anchor domain induce its reversible unfolding above 34 degrees C. The overall triple helix denaturation temperature is reduced by 5-6 degrees C, similar to complete N-anchor removal. N-propeptide partially restores the stability of mutant procollagen but not sufficiently to prevent N-anchor unfolding and a conformational change at the N-propeptide cleavage site. The ensuing failure of N-proteinase to cleave at the misfolded site leads to incorporation of pN-collagen into fibrils. Similar, but weaker, effects are caused by G88E substitution in the adjacent triplet, which appears to alter N-anchor structure as well. As in Ehlers-Danlos syndrome (EDS) VIIA/B, fibrils containing pN-collagen are thinner and weaker causing EDS-like laxity of large and small joints and paraspinal ligaments. However, distinct structural consequences of N-anchor destabilization result in a distinct alpha1(I)-osteogenesis imperfecta (OI)/EDS phenotype.     

Predicting helical segments in proteins by a helix-coil transition theory with parameters derived from a structural database of proteins

A novel helix-coil transition theory has been developed. This new theory contains more types of interactions than similar theories developed earlier. The parameters of the models were obtained from a database of 351 nonhomologous proteins. No manual adjustment of the parameters was performed. The interaction parameters obtained in this manner were found to be physically meaningful, consistent with current understanding of helix stabilizing/destabilizing interactions. Novel insights into helix stabilizing/destabilizing interactions have also emerged from this analysis. The theory developed here worked well in sorting out helical residues from amino acid sequences. If the theory was forced to make prediction on every residue of a given amino acid sequence, its performance was the best among ten other secondary structural prediction algorithms in distinguishing helical residues from nonhelical ones. The theory worked even better if one only required it to make prediction on residues that were “predictable” (identifiable by the theory); >90% predictive reliability could be achieved. The helical residues or segments identified by the helix-coil transition theory can be used as secondary structural contraints to speed up the prediction of the three-dimensional structure of a protein by reducing the dimension of a computational protein folding problem. Possible further improvements of this helix-coil transition theory are also discussed. Proteins 28:344–359, 1997. © 1997 Wiley-Liss, Inc.     

Unfolding of apomyoglobin helices by synchrotron radiolysis and mass spectrometry.

The synchrotron X-ray protein radiolysis technique is based on a quantitative determination of the extent and the site of millisecond radiolytic oxidation of amino-acid side chains by mass spectrometry. The amino acids most susceptible to radiolytic oxidation are cysteine, methionine, phenylalanine, tyrosine, tryptophan, proline, histidine, and leucine. These residues serve as reactive markers within a protein structure that can be used to monitor changes in solvent accessibility during folding or as part of macromolecular interactions. To monitor the unfolding, the extent of radiolytic products of side chains of reactive amino acids is quantitatively measured by mass spectrometry as a function of the denaturant concentration following proteolysis. This approach provides site-specific unfolding isotherms for various segments of a protein without the use of mutation or labeling techniques. Application of this technique to the equilibrium urea unfolding of apomyoglobin at pH 7.8 has demonstrated the cooperative unfolding of helices A to C consistent with midpoints, DeltaG, and m values derived from fluorescence data. The G helix, in contrast, showed a local unfolding behavior. The similarity of the thermodynamic data derived by this synchrotron-based method for helix A (containing two oxidizable tryptophan residues) to that of the fluorescence data indicates that the limited oxidation of proteins by exposure to X-rays on millisecond timescales does not alter the structure of apomyglobin. This supports the viability of the method for the study of protein folding and the mapping of protein interaction sites.     

Detecting molecular interactions that stabilize, activate and guide ligand-binding of the sodium/proton antiporter MjNhaP1 from Methanococcus jannaschii

Integral membrane proteins are involved in virtually every cellular process. Precisely regulating these machineries would allow controlling many human and vertebrate diseases. Embedded into cellular membranes, membrane proteins establish molecular interactions that sensitively react to environmental changes and to molecular compounds, such as ligands or inhibitors. We applied atomic force microscopy (AFM) to image the Na(+)/H(+) antiporter MjNhaP1 from Methanococcus jannaschii, and single-molecule force spectroscopy (SMFS) to probe molecular interactions that drive the protein structure-function relationship. High-resolution AFM topographs showed the dimeric assembly of MjNhaP1 being reconstituted into a lipid bilayer. SMFS of MjNhaP1 unraveled molecular interactions stabilizing individual structural domains. Transmembrane domains exhibited certain probabilities to unfold individually or cooperatively with other domains resulting in different unfolding pathways. Helices VIII and X established pH sensitive interactions altering significantly upon MjNhaP1 activation, while removal of the ligand (Na(+)) destabilized the entire antiporter except helix VIII. It is assumed that Asp234/235 of helix VIII are involved in the ligand-binding site and that helix X plays a functional role in the activation of the transporter.     

Stability of bacteriorhodopsin alpha-helices and loops analyzed by single-molecule force spectroscopy

The combination of high-resolution atomic force microscopy imaging and single-molecule force spectroscopy allows the identification, selection, and mechanical investigation of individual proteins. In a recent paper we had used this technique to unfold and extract single bacteriorhodopsins (BRs) from native purple membrane patches. We show that subsets of the unfolding spectra can be classified and grouped to reveal detailed insight into the individualism of the unfolding pathways. We have further developed this technique and analysis to report here on the influence of pH effects and local mutations on the stability of individual structural elements of BR against mechanical unfolding. We found that, although the seven transmembrane alpha-helices predominantly unfold in pairs, each of the helices may also unfold individually and in some cases even only partially. Additionally, intermittent states in the unfolding process were found, which are associated with the stretching of the extracellular loops connecting the alpha-helices. This suggests that polypeptide loops potentially act as a barrier to unfolding and contribute significantly to the structural stability of BR. Chemical removal of the Schiff base, the covalent linkage of the photoactive retinal to the helix G, resulted in a predominantly two-step unfolding of this helix. It is concluded that the covalent linkage of the retinal to helix G stabilizes the structure of BR. Trapping mutant D96N in the M state of the proton pumping photocycle did not affect the unfolding barriers of BR.     

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.     

On the thermal unfolding character of globular proteins

A theoretical model is presented to study the stepwise thermal unfolding of globular proteins using the stabilizing/destabilizing characters of amino acid residues in protein crystals. A multiple regression relation connecting the melting temperature and the amounts of stabilizing and destabilizing groups of residues in a protein, when used for the thermal behavior of peptide segments, provides reliable results on the stepwise unfolding nature of the protein. In ribonuclease A, the shell residues 16–22 are predicted to unfold earlier in the temperature range 30–45°C; the -sheet structures undergo thermal denaturation as a single cooperative unit and there is evidence indicating the segment 106–118 as a nucleation site. In ribonuclease S, the S-peptide unfolds earlier than S-protein. The predicted average and the range of melting temperatures, and the folding pathways of a set of globular proteins, agree very well with the experimental results. The results obtained in the present study indicate that (i) most of the nucleation parts possess high relative thermal stability, (ii) the unfolded state retains some residual structure, and (iii) some segments undergo gradual and overlapping thermal denaturation.     

Local secondary structure in ribonuclease A denatured by guanidine · HCl near 1 °C

Recent work has shown that a short α-helix can be stable in water near 1 °C when stabilized by specific interactions between side-chains, while earlier “host-guest” results with random copolymers have shown that a short α-helix is unstable in water at all temperatures in the absence of stabilizing side-chain interactions. As regards the mechanism of protein folding, it is now reasonable on energetic grounds to consider isolated α-helices and β-sheets as the first intermediates on the pathway of protein folding. Proton nuclear magnetic resonance is used here to detect isolated secondary structures in ribonuclease A denatured by guanidine · HCl (GuHCl). Temperatures near 1 °C are used because the low-temperature stability of the C-peptide helix may be a general property of isolated secondary structures in water.Our procedure is to titrate with GuHCl the C2H resonance lines of the four histidine residues of denatured ribonuclease A. Studies of model peptides (C-peptide (lactone) and C-peptide carboxylate, residues 1 to 13 of ribonuclease A; S-peptide, residues 1 to 20) show linear titration curves for the C2H resonance of His12 above 0.5 M-GuHCl, once helix unfolding is complete. Deviations from this line are used to monitor helix formation. The GuHCl titration curves of the other three histidine residues are also linear, once unfolding is complete. The results show that the helix found in C-peptide and S-peptide is also found in denatured ribonuclease A, where it behaves as an isolated helix not stabilized significantly by interactions with other chain segments. Studies of denatured S-protein show that the remaining three His residues, His48, His105 and His119, are involved in structure only below 1 m-GuHCl at 9 °C, pH 1.9. The nature of this structure is not known. The main conclusion from this work is that the His12 helix can be observed as a stable, isolated helix in denatured ribonuclease A near 1 °C, and that none of the other three His residues is involved in a comparably stable local structure. In native ribonuclease A, His12 is within an α-helix and the other three His residues are involved in a 3-stranded β-sheet structure.The helix-coil transition of C-peptide has also been studied for other side-chain resonances by GuHCl titration. Typically, but not always, the titration curves are linear after helix unfolding takes place and resonance lines from different residues of the same amino acid type can be resolved in GuHCl solutions. This is true of the four histidine residues of ribonuclease A although their pK values in 5 m-GuHCl are nearly the same. In C-peptide, the βCH₃ resonance of Ala6 is affected strongly by GuHCl while the lines of Ala4 and Ala5 are shifted only weakly by GuHCl. Evidently the interactions between GuHCl and side-chains in an unfolded peptide depend upon neighboring groups.     

Experimental comparison of energy landscape features of ubiquitin family proteins

Small ubiquitin-related modifiers (SUMO1 and SUMO2) are ubiquitin family proteins, structurally similar to ubiquitin, differing in terms of their amino acid sequence and functions. Therefore, they provide a great platform for investigating sequence-structure-stability-function relationship. Here, we used chemical denaturation in comparing the folding-unfolding pathways of the SUMO proteins with their structural homologue ubiquitin (UF45W-pseudo wild-type [WT] tryptophan variant) with structurally analogous tryptophan mutations (SUMO1 [S1F66W], SUMO2 [S2F62W]). Equilibrium denaturation studies report that ubiquitin is the most stable protein among the three. The observed denaturant-dependent folding rates of SUMOs are much lower than ubiquitin and primarily exhibit a two-state folding pathway unlike ubiquitin, which has a kinetic folding intermediate. We hypothesize that, as SUMO proteins start off as slow folders, they avoid stabilizing their folding intermediates and the presence of which might further slow-down their folding rates. The denaturant-dependent unfolding of ubiquitin is the fastest, followed by SUMO2, and slowest for SUMO1. However, the spontaneous unfolding rate constant is the lowest for ubiquitin (~40 times), and similar for SUMOs. This correlation between thermodynamic stability and kinetic stability is achieved by having different unfolding transition state positions with respect to the solvent-accessible surface area, as quantified by the Tanford β _u values: ubiquitin (0.42) > SUMO2 (0.20) > SUMO1 (0.16). The results presented here highlight the unique energy landscape features which help in optimizing the folding-unfolding rates within a structurally homologous protein family.     

Helix capping in the GCN4 leucine zipper.

Capping interactions associated with specific sequences at or near the ends of alpha-helices are important determinants of the stability of protein secondary and tertiary structure. We investigate here the role of the helix-capping motif Ser-X-X-Glu, a sequence that occurs frequently at the N termini of alpha helices in proteins, on the conformation and stability of the GCN4 leucine zipper. The 1.8 A resolution crystal structure of the capped molecule reveals distinct conformations, packing geometries and hydrogen-bonding networks at the amino terminus of the two helices in the leucine zipper dimer. The free energy of helix stabilization associated with the hydrogen-bonding and hydrophobic interactions in this capping structure is -1.2 kcal/mol, evaluated from thermal unfolding experiments. A single cap thus contributes appreciably to stabilizing the terminated helix and thereby the native state. These results suggest that helix capping plays a further role in protein folding, providing a sensitive connector linking alpha-helix formation to the developing tertiary structure of a protein.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

How Ala-->Gly mutations in different helices affect the stability of the apomyoglobin molten globule

The apomyoglobin molten globule has a complex, partly folded structure with a folded A[B]GH subdomain; the factors determining its stability are not yet known in detail. Ala-->Gly mutations, made at solvent-exposed positions, are used to probe the role of helix propensity of individual helices in stabilizing the molten globule. Molten globule stability is measured by reversible urea unfolding, monitored both by circular dichroism and by tryptophan fluorescence. Two-state unfolding is tested by superposition of these two unfolding curves, and stability data are reported only for variants which satisfy the superposition test. Results for sites Q8 in the A helix and E109 in the G helix confirm that the helix propensities of the A and G helices both strongly affect molten globule stability, in contrast to results for the G65A/G73A double mutant which show that changing the helix propensity of the E-helix sequence has no significant stabilizing effect. Changing the helix propensity of the B-helix sequence with the G23A/G25A double mutant affects molten globule stability to an intermediate extent, confirming an earlier report that this mutant has increased stability. These results are consistent with the bipartite structure for the molten globule in which the A, G, and H helices are stably folded, while the long E helix is unfolded and the B helix has intermediate stability. Some differences are found in the shapes of the unfolding curves of different mutants even though they satisfy the superposition test for two-state unfolding, and possible explanations are discussed.     

Inter- and intrasubunit interactions during the formation of RNA polymerase assembly intermediate

We used yeast two-hybrid and in vitro co-immobilization assays to study the interaction between the Escherichia coli RNA polymerase (RNAP) alpha and beta subunits during the formation of alpha(2)beta, a physiological RNAP assembly intermediate. We show that a 430-amino acid-long fragment containing beta conserved segments F, G, H, and a short part of segment I forms a minimal domain capable of specific interaction with alpha. The alpha-interacting domain is held together by protein-protein interactions between beta segments F and I. Residues in catalytically important beta segments H and I directly participate in alpha binding; substitutions of strictly conserved segment H Asp(1084) and segment I Gly(1215) abolish alpha(2)beta formation in vitro and are lethal in vivo. The importance of these beta amino acids in alpha binding is fully supported by the structural model of the Thermus aquaticus RNAP core enzyme. We also demonstrate that determinants of RNAP assembly are conserved, and that a homologue of beta Asp(1084) in A135, the beta-like subunit of yeast RNAP I, is responsible for interaction with AC40, the largest alpha-like subunit. However, the A135-AC40 interaction is weak compared with the E. coli alpha-beta interaction, and A135 mutation that abolishes the interaction is phenotypically silent. The results suggest that in eukaryotes additional RNAP subunits orchestrate the enzyme assembly by stabilizing weak, but specific interactions of core subunits.     

A heuristic approach to predicting the tertiary structure of bovine somatotropin

A combination of a heuristic approach and energy minimization was used to predict the three-dimensional structure of bovine somatotropin (bSt), also known as bovine growth hormone, a protein of 191 amino acids. The starting points for energy minimizations were generated from the following two types of inputs: (a) the amino acid sequence and (b) the heuristic inputs, which were derived according to physical, chemical, and biological principles by piecing together all useful information available. The predicted 3-D structure of the bSt molecule has all the features observed in four-helix bundle proteins. The four alpha-helices in bSt are intimately packed to form an assembly with an approximately square cross section. All the adjacent alpha-helices are antiparallel, with a somewhat tilted angle between each of the adjacent pairs so that the assembly of the four helices looks like a left-handed twisted bundle. There are two disulfide bonds in the bSt structure: one "hooking" the middle of a long loop with helix 4 so as to pull the long loop onto the surface of the helix bundle and the other "hooking" the C-terminal segment with the same helix so as to force the C-terminal segment to bend toward the helix bundle. As a consequence, a considerable part of the surface of the four-helix bundle is closely packed or intimately embraced by the loop segments. The predicted bSt structure has a hydrophobic core and a hydrophilic exterior surface. The energetic analysis of the predicted bSt structure indicates that the interaction between helices and loops plays a dominant role in stabilizing the four-helix bundle structure from the viewpoint of both electrostatic and nonbonded interactions. A technique called FOLD was meanwhile developed, by which one can fold a polypeptide chain into any shape as desired. This tool proved to be very useful during the heuristic model-building process.     

Deoxymyoglobin studied by the conformational normal mode analysis. I. Dynamics of globin and the heme-globin interaction

Dynamic properties of deoxymyoglobin are studied theoretically by the analysis of conformational fluctuations. Root-mean-square atomic fluctuations and distance fluctuations between different segments reveal the mechanical construction of the molecule. Eight alpha-helices behave as relatively rigid bodies and corner regions are more flexible, showing larger fluctuations. More particularly, corner regions EF and GH are specific in that flanking alpha-helices extend their rigidity up to a point in the corner region and the two rigid segments are connected flexibly at that point. The FG corner is exceptional. A segment from the F helix to the beginning of the G helix, in which the FG corner is included, becomes relatively rigid by means of strong interactions with the heme group. The whole myoglobin molecule is divided into two large units of motion, one extending from the B to the E helix, and the other from the F to the H helix. These two units are connected covalently by the EF corner. However, dynamic interactions between these two units take place mainly through contacts between helices B and G and not through the EF corner. From correlation coefficients between fluctuational motions of residues and the heme group, 55 residues are identified as having strong dynamic interactions with the heme moiety. Among them, 18 residues in the three segments, one consisting of residues from the C helix to the CD corner, a second consisting of the E helix, and a third from the F helix to the beginning of the G helix, are in close contact with the heme group. Twenty-two of the 55 residues are within four residues of the 18 residues in their sequential residue number and are more than 3 A away from the heme group. The other 15 residues are located further in the sequential residue number and are all found in helices A and H. They are more than 6 A away from the heme group. By the use of correlation coefficients of fluctuations between residues, it is found that dynamic interaction with the heme group is transmitted to the A helix and the beginning of the H helix in the direction Leu(E15)----[Val(All) and Trp(A12)]. The transmission to the C-terminal end of the H helix is mediated by a long segment, from the end of the EF corner to the beginning of the G helix, that lies on the heme group and has close contacts over a wide range.(ABSTRACT TRUNCATED AT 400 WORDS)