Comparing four different approaches for the determination of inter‐residue interactions provides insight for the structure prediction of helical membrane proteins

Studying inter‐residue interactions provides insight into the folding and stability of both soluble and membrane proteins and is essential for developing computational tools for protein structure prediction. As the first step, various approaches for elucidating such interactions within protein structures have been proposed and proven useful. Since different approaches may grasp different aspects of protein structural folds, it is of interest to systematically compare them. In this work, we applied four approaches for determining inter‐residue interactions to the analysis of three distinct structure datasets of helical membrane proteins and compared their correlation to the three individual quality measures of structures in these datasets. These datasets included one of 35 structures of rhodopsin receptors and bacterial rhodopsins determined at various resolutions, one derived from the HOMEP benchmark dataset previously reported, and one comprising of 139 homology models. It was found that the correlation between the average number of inter‐residue interactions obtained by applying the four approaches and the available structure quality measures varied quite significantly among them. The best correlation was achieved by the approach focusing exclusively on favorable inter‐residue interactions. These results provide interesting insight for the development of objective quality measure for the structure prediction of helical membrane proteins. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 547–556, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Solution structure of a novel T‐cell adhesion inhibitor derived from the fragment of ICAM‐1 receptor: Cyclo(1,8)‐Cys‐Pro‐Arg‐Gly‐Gly‐Ser‐Val‐Cys

This study is aimed at elucidating the structure of a novel T‐cell adhesion inhibitor, cyclo(1,8)‐CPRGGSVC using one‐ and two‐dimensional (2D) ¹H NMR and molecular dynamics (MD) simulation. The peptide is derived from the sequence of its parent peptide cIBR (cyclo(1,12)‐PenPRGGSVLVTGC), which is a fragment of intercellular adhesion molecule‐1 (ICAM‐1). Our previous results show that the cyclo(1,8)‐CPRGGSVC peptide binds to the LFA‐1 I‐domain and inhibits heterotypic T‐cell adhesion, presumably by blocking the LFA‐1/ICAM‐1 interactions. The structure of the peptide was determined using NMR and MD simulation in aqueous solution. Our results indicate that the peptide adopts type‐I β‐turn conformation at the Pro2‐Arg3‐Gly4‐Gly5 (PRGG) sequence. The β‐turn structure at the PRGG motif is well conserved in cIBR peptide and ICAM‐1 receptor, which suggests the importance of the PRGG motif for the biological activity of cyclo(1,8)‐CPRGGSVC peptide. Meanwhile, the Gly5‐Ser6‐Val7‐Cys8‐Cys1 (GSVCC) sequence forms a “turn‐like” random coil structure that does not belong to any structured motif. Therefore, cyclo(1,8)‐CPRGGSVC peptide has only one structured region at the PRGG sequence, which may play an important role in the binding of the peptide to the LFA‐1 I‐domain. The conserved β‐turn conformation of the PRGG motif in ICAM‐1, cIBR, and cyclo(1,8)‐CPRGGSVC peptides can potentially be used to design peptidomimetics. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 633–641, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Peptidic modulators of protein‐protein interactions: Progress and challenges in computational design

With the decline in productivity of drug‐development efforts, novel approaches to rational drug design are being introduced and developed. Naturally occurring and synthetic peptides are emerging as novel promising compounds that can specifically and efficiently modulate signaling pathways in vitro and in vivo. We describe sequence‐based approaches that use peptides to mimic proteins in order to inhibit the interaction of the mimicked protein with its partners. We then discuss a structure‐based approach, in which protein‐peptide complex structures are used to rationally design and optimize peptidic inhibitors. We survey flexible peptide docking techniques and discuss current challenges and future directions in the rational design of peptidic inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 505–513, 2009. This article was originally published online as an accepted preprint. The “Published Online”date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Ion channel activity of transmembrane segment 6 of Escherichia coli proton‐dependent manganese transporter

Synthetic peptides corresponding to the sixth transmembrane segment (TMS6) of secondary‐active transporter MntH (Proton‐dependent Manganese Transporter) from Escherichia coli and its two mutations in the functionally important conserved histidine residue were used as a model for structure–function study of MntH. The secondary structure of the peptides was estimated in different environments using circular dichroism spectroscopy. These peptides interacted with and adopted helical conformations in lipid membranes. Electrophysiological experiments demonstrated that TMS6 was able to form multi‐state ion channels in model biological membranes. Electrophysiological properties of these weakly cation‐selective ion channels were strongly dependent on the surrounding pH. Manganese ion, as a physiological substrate of MntH, enhanced the conductivity of TMS6 channels, influenced the transition between closed and open states, and affected the peptide conformations. Moreover, functional properties of peptides carrying two different mutations of His²¹¹ were analogous to in vivo functional characteristics of Nramp/MntH proteins mutated at homologous residues. Hence, a single functionally important TMS can retain some of the functional properties of the full‐length protein. These findings could contribute to understanding the structure–function relationship at the molecular level. However it remains unclear to what extent the peptide‐specific channel activity represents a functional aspect of the full‐length membrane carrier protein. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 718–726, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Sequence‐dependent folding and unfolding of ligand‐bound purine riboswitches

Riboswitch regulation of gene expression requires ligand‐mediated RNA folding. From the fluorescence lifetime distribution of bound 2‐aminopurine ligand, we resolve three RNA conformers (C_o, C_i, C_c) of the liganded G‐ and A‐sensing riboswitches from Bacillus subtilis. The ligand binding affinities, and sensitivity to Mg²⁺, together with results from mutagenesis, suggest that C_o and C_i are partially unfolded species compromised in key loop‐loop interactions present in the fully folded C_c. These data verify that the ligand‐bound riboswitches may dynamically fold and unfold in solution, and reveal differences in the distribution of folded states between two structurally homologous purine riboswitches: Ligand‐mediated folding of the G‐sensing riboswitch is more effective, less dependent on Mg²⁺, and less debilitated by mutation, than the A‐sensing riboswitch, which remains more unfolded in its liganded state. We propose that these sequence‐dependent RNA dynamics, which adjust the balance of ligand‐mediated folding and unfolding, enable different degrees of kinetic discrimination in ligand binding, and fine‐tuning of gene regulatory mechanisms. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 953–965, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Chaperone‐like activities of different molecular forms of β‐casein. Importance of polarity of N‐terminal hydrophilic domain

As a member of intrinsically unstructured protein family, β‐casein (β‐CN) contains relatively high amount of prolyl residues, adopts noncompact and flexible structure and exhibits chaperone‐like activity in vitro. Like many chaperones, native β‐CN does not contain cysteinyl residues and exhibits strong tendencies for self‐association. The chaperone‐like activities of three recombinant β‐CNs wild type (WT) β‐CN, C4 β‐CN (with cysteinyl residue in position 4) and C208 β‐CN (with cysteinyl residue in position 208), expressed and purified from E. coli, which, consequently, lack the phosphorylated residues, were examined and compared with that of native β‐CN using insulin and alcohol dehydrogenase as target/substrate proteins. The dimers (β‐CND) of C4‐β‐CN and C208 β‐CN were also studied and their chaperone‐like activities were compared with those of their monomeric forms. Lacking phosphorylation, WT β‐CN, C208 β‐CN, C4 β‐CN and C4 β‐CND exhibited significantly lower chaperone‐like activities than native β‐CN. Dimerization of C208 β‐CN with two distal hydrophilic domains considerably improved its chaperone‐like activity in comparison with its monomeric form. The obtained results demonstrate the significant role played by the polar contributions of phosphorylated residues and N‐terminal hydrophilic domain as important functional elements in enhancing the chaperone‐like activity of native β‐CN. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 623–632, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

A proposed model for dragline spider silk self‐assembly: Insights from the effect of the repetitive domain size on fiber properties

Dragline spider silk has been intensively studied for its superior qualities as a biomaterial. In previous studies, we made use of the baculovirus mediated expression system for the production of a recombinant Araneus diadematus spider silk dragline ADF4 protein and its self‐assembly into intricate fibers in host insect cells. In this study, our aim was to explore the function of the major repetitive domain of the dragline spider silk. Thus, we generated an array of synthetic proteins, each containing a different number of identical repeats up to the largest recombinantly expressed spider silk to date. Study of the self‐assembly properties of these proteins showed that depending on the increasing number of repeats they give rise to different assembly phenotypes, from a fully soluble protein to bona fide fibers with superior qualities. The different assembly forms, the corresponding chemical resistance properties obtained as well as ultrastructural studies, revealed novel insights concerning the structure and intermolecular interactions of the repetitive and nonrepetitive domains. Based on these observations and current knowledge in the field, we hereby present a comprehensive hypothetical model for the mechanism of dragline silk self‐assembly and fiber formation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 458–468, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Kinetics of the self‐aggregation and film formation of poly‐L‐proline at high temperatures explored by circular dichroism spectroscopy

Poly‐L ‐proline has been used as a model system for various purposes over a period of more than 60 years. Its relevance among the protein/peptide community stems from its use as a reference system for determining the conformational distributions of unfolded peptides and proteins, its use as a molecular ruler, and from the pivotal role of proline residues in conformational transitions and protein–protein interactions. While several studies indicate that polyproline can aggregate and precipitate in aqueous solution, a systematic study of the aggregation process is still outstanding. We found, by means of UV‐circular dichroism and IR measurements, that polyproline is predominately monomeric at room temperature at millimolar concentrations. Upon heating, the polypeptide stays in its monomeric state until the temperature reaches a threshold of ca. 60°C. At higher temperatures, the peptide aggregates as a film on the inside surface of the employed cuvette. The process proceeds on a time scale of 10³ s and can best be described by a bi‐exponential relaxation function. The respective CD and IR spectra are qualitatively different from the canonical spectra of polyproline in aqueous solution, and are indicative of a highly packed state. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 451–457, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Conformational stability of the full‐atom hexameric model of the ClpB chaperone from Escherichia coli

The Escherichia coli heat shock protein ClpB, a member of the Hsp100 family, plays a crucial role in cellular thermotolerance. In co‐operation with the Hsp70 chaperone system, it is able to solubilize proteins aggregated by heat shock conditions and refold them into the native state in an ATP‐dependent way. It was established that the mechanism of ClpB action depends on the formation of a ring‐shaped hexameric structure and the translocation of a protein substrate through an axial channel. The structural aspects of this process are not fully known. By means of homology modeling and protein–protein docking, we obtained a model of the hexameric arrangement of the full‐length ClpB protein complexed with ATP. A molecular dynamics simulation of this model was performed to assess its flexibility and conformational stability. The high mobility of the “linker” M‐domain, essential for the renaturing activity of ClpB, was demonstrated, and the size and shape of central channel were analyzed. In this model, we propose the coordinates for a loop between b4 and B6 structural elements, not defined in previous structural research, which faces the inside of the channel and may therefore play a role in substrate translocation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 47–60, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Structural transition from dimeric to tetrameric i‐motif,caused by the presence of TAA at the 3′‐end of human telomeric C‐rich sequence

Widely dispersed in genomic DNA, the tandem C‐rich repetitive stretches may fold below physiological pH, into i‐motif structures, stabilized by C·C⁺ pairing. Herein, structural status of a 9‐mer stretch d(CCCTAACCC), [the truncated double repeat of human telomeric sequence], and its extended version, comprising of additional ? TAA segment at the 3′‐end, representing the complete double repeat d(CCCTAACCCTAA), has been investigated. The pH dependent monophasic UV‐melting, Gel and CD data suggested that while the truncated version adopts a bimolecular i‐motif structure, its complete double repeat (12‐mer) sequence exists in two (bimolecular and tetramolecular) forms. A model is proposed for the tetramolecular i‐motif with conventional C · C⁺ base pairs, additionally stabilized by asymmetric A · A base pairs at the ?3′ TAA flanking ends and Watson–Crick A · T hydrogen bonding between intervening bases on antiparallel strands. Expanding the known topologies of DNA i‐motifs, such atypical geometries of i‐motifs may have implications in their recognition by proteins. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 150–160, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Glycerol‐induced folding of unstructured disulfide‐deficient lysozyme into a native‐like conformation

2SS[6‐127,64‐80] variant of lysozyme which has two disulfide bridges, Cys6‐Cys127 and Cys64‐Cys80, and lacks the other two disulfide bridges, Cys30‐Cys115 and Cys76‐Cys94, was quite unstructured in water, but a part of the polypeptide chain was gradually frozen into a native‐like conformation with increasing glycerol concentration. It was monitored from the protection factors of amide hydrogens against H/D exchange. In solution containing various concentrations of glycerol, H/D exchange reactions were carried out at pH* 3.0 and 4°C. Then, ¹H‐¹⁵N‐HSQC spectra of partially deuterated protein were measured in a quenching buffer for H/D exchange (95% DMSO/5% D₂O mixture at pH* 5.5 adjusted with dichloroacetate). In a solution of 10% glycerol, the protection factors were nearly equal to 10 at most of residues. With increasing glycerol concentration, some selected regions were further protected, and their protection factors reached about a 1000 in 30% glycerol solution. The highly protected residues were included in A‐, B‐, and C‐helices and β3‐strand, and especially centered on Ile 55 and Leu 56. In 2SS[6‐127,64‐80], long‐range interactions were recovered due to the preferential hydration by glycerol in the hydrophobic box of the α‐domain. Glycerol‐induced recovering of the native‐like structure is discussed from the viewpoint of molten globules growing with the protein folding. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 665–675, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Study of the interactions between a proline‐rich protein and a flavan‐3‐ol by NMR: Residual structures in the natively unfolded protein provides anchorage points for the ligands

Astringency is one of the major organoleptic properties of food and beverages that are made from plants, such as tea, chocolate, beer, or red wine. This sensation is thought to be due to interactions between tannins and salivary proline‐rich proteins, which are natively unfolded proteins. A human salivary proline‐rich protein, namely IB‐5, was produced by the recombinant method. Its interactions with a model tannin, epigallocatechin gallate (EGCG), the major flavan‐3‐ol in green tea, were studied here. Circular dichroism experiments showed that IB‐5 presents residual structures (PPII helices) when the ionic strength is close to that in saliva. In the presence of these residual structures, IB‐5 undergoes an increase in structural content upon binding to EGCG. NMR data corroborated the presence of preformed structural elements within the protein prior to binding and a partial assignment was proposed, showing partial structuration. TOCSY experiments showed that amino acids that are involved in PPII helices are more likely to interact with EGCG than those in random coil regions, as if they were anchorage points for the ligand. The signal from IB‐5 in the DOSY NMR spectrum revealed an increase in polydispersity upon addition of EGCG while the mean hydrodynamic radius remained unchanged. This strongly suggests the formation of IB‐5/EGCG aggregates. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 745–756, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Field‐based comparison of ligand and coactivator binding sites of nuclear receptors

A structure‐based comparison of the ligand‐binding domains of 35 nuclear receptors from five different subfamilies is presented. Their ligand and coactivator binding sites are characterized using knowledge‐based contact preference fields for hydrophobic and hydrophilic interactions implemented in the MOE modeling environment. Additionally, for polar knowledge‐based field points the preference for negative or positive electrostatic interactions is estimated using the Poisson‐Boltzmann equation. These molecular‐interaction fields are used to cluster the nuclear receptor family based on similarities of their binding sites. By analyzing the similarities and differences of hydrophobic and polar fields in binding pockets of related receptors it is possible to identify conserved interactions in ligand and coactivator binding pockets, which support e.g. design of specific ligands during lead optimization or virtual screening as docking filter. Examples of remarkable similarities between ligand binding sites of members from phylogenetically different nuclear receptor families (RXR, RAR, HNF4, NR5) and differences between closely related subtypes (LXR, RAR, TR) are discussed in more detail. Significant similarities and differences of coactivator binding sites are shown for NR3Cs, LXRs and PPARs. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 884–894, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Hsp90 and co‐chaperones twist the functions of diverse client proteins

Hsp90 molecular chaperones are required for the stability and activity of a diverse range of client proteins that have critical roles in signal transduction, cellular trafficking, chromatin remodeling, cell growth, differentiation, and reproduction. Mammalian cells contain three types of Hsp90s: cytosolic Hsp90, mitochondrial Trap‐1, and Grp94 of the endoplasmic reticulum. Each of the Hsp90s, as well as the bacterial homolog, HtpG, hydrolyze ATP and undergo similar conformational changes. Unlike the other forms of Hsp90, cytosolic Hsp90 function is dependent on a battery of co‐chaperone proteins that regulate the ATPase activity of Hsp90 or direct Hsp90 to interact with specific client proteins. This review will summarize what is known about Hsp90's ability to mediate the folding and activation of diverse client proteins that contribute to human diseases, such as cancer and fungal and viral infections. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 211–217, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

The role of the Val57 amino‐acid residue in the hinge loop of the human cystatin C. Conformational studies of the beta2‐L1‐beta3 segments of wild‐type human cystatin C and its mutants

Human cystatin C (HCC) is one of the amyloidogenic proteins to be shown to oligomerize via a three‐dimensional domain swapping mechanism. This process precedes the formation of a stable dimer and proceeds particularly easily in the case of the L68Q mutant. According to the proposed mechanism, dimerization of the HCC precedes conformational changes within the β2 and β3 strands. In this article, we present conformational studies, using circular dichroism and MD methods, of the β2‐L1‐β3 (His43‐Thr72) fragment of the HCC involved in HCC dimer formation. We also carried out studies of the β2‐L1‐β3 peptide, in which the Val57 residue was replaced by residues promoting β‐turn structure formation (Asp, Asn, or Pro). The present study established that point mutation could modify the structure of the L1 loop in the β‐hairpin peptide. Our results showed that the L1 loop in the peptide excised from human cystatin C is broader than that in cystatin C. In the HCC protein, broadening of the L1 loop together with the unfavorable L68Q mutation in the hydrophobic pocket could be a force sufficient to cause the partial unfolding and then the opening of HCC or its L68Q mutant structure for further dimerization. We presume further that the Asp57 and Asn57 mutations in the L1 loop of HCC could stabilize the closed form of HCC, whereas the Pro57 mutation could lead to the opening of the HCC structure and then to dimer/oligomer formation. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 373–383, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Stable intermediates determine proteins' primary unfolding sites in the presence of surfactants

Despite detailed knowledge of the overall structural changes and stoichiometries of surfactant binding, little is known about which protein regions constitute the preferred sites of attack for initial unfolding. Here we have exposed three proteins to limited proteolysis at anionic (SDS) and cationic (DTAC) surfactant concentrations corresponding to specific conformational transitions, using the surfactant‐robust broad‐specificity proteases Savinase and Alcalase. Cleavage sites are identified by SDS‐PAGE and N‐terminal sequencing. We observe well‐defined cleavage fragments, which suggest that flexibility is limited to certain regions of the protein. Cleavage sites for α‐lactalbumin and myoglobin correspond to regions identified in other studies as partially unfolded at low pH or in the presence of organic solvents. For Tnfn3, which does not form partially folded structures under other conditions, cleavage sites can be rationalized from the structure of the protein's folding transition state and the position of loops in the native state. Nevertheless, they are more sensitive to choice of surfactant and protease, probably reflecting a heterogeneous and fluctuating ensemble of partially unfolded structures. Thus, for proteins accumulating stable intermediates on the folding pathway, surfactants encourage the formation of these states, while the situation is more complex for proteins that do not form these intermediates. © 2008 Wiley Periodicals, Inc. Biopolymers 91: 221–231, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com