β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

The structure of a shellfish specific GST class glutathione S‐transferase from antarctic bivalve Laternula elliptica reveals novel active site architecture

Glutathione‐S‐transferases have been identified in all the living species examined so far, yet little is known about their function in marine organisms. In a previous report, the recently identified GST from Antarctic bivalve Laternula elliptica (LeGST) was classified into the rho class GST, but there are several unique features of LeGST that may justify reclassification, which could represent specific shellfish GSTs. Here, we determined the crystal structure of LeGST, which is a shellfish specific class of GST. The structural analysis showed that the relatively open and wide hydrophobic H‐site of the LeGST allows this GST to accommodate various substrates. These results suggest that the H‐site of LeGST may be the result of adaptation to their environments as sedentary organisms. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Temperature‐induced partially unfolded state of hUBF HMG Box‐5: Conformational and dynamic investigations of the Box‐5 thermal intermediate ensemble

Human upstream binding factor (hUBF) HMG Box‐5 is a highly conserved protein domain, containing 84 amino acids and belonging to the family of the nonspecific DNA‐binding HMG boxes. Its native structure adopts a twisted L shape, which consists of three α‐helices and two hydrophobic cores: the major wing and the minor wing. In this article, we report a reversible three‐state thermal unfolding equilibrium of hUBF HMG Box‐5, which is investigated by differential scanning calorimetry (DSC), circular dichroism spectroscopy, fluorescence spectroscopy, and NMR spectroscopy. DSC data show that Box‐5 unfolds reversibly in two separate stages. Spectroscopic analyses suggest that different structural elements exhibit noncooperative transitions during the unfolding process and that the major form of the Box‐5 thermal intermediate ensemble at 55°C shows partially unfolded characteristics. Compared with previous thermal stability studies of other boxes, it appears that Box‐5 possesses a more stable major wing and two well separated subdomains. NMR chemical shift index and sequential ¹H^N_i‐¹H^N_i+1 NOE analyses indicate that helices 1 and 2 are native‐like in the thermal intermediate ensemble, while helix 3 is partially unfolded. Detailed NMR relaxation dynamics are compared between the native state and the intermediate ensemble. Our results implicate a fluid helix‐turn‐helix folding model of Box‐5, where helices 1 and 2 potentially form the helix 1‐turn‐helix 2 motif in the intermediate, while helix 3 is consolidated only as two hydrophobic cores form to stabilize the native structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural basis of urea-induced unfolding of Fasciola gigantica glutathione S-transferase

Glutathione S-transferases (GSTs) are enzymes that are involved in the detoxification of harmful electrophilic endogenous and exogenous compounds by conjugating with glutathione (GSH). The liver fluke GSTs have multifunctional roles in the host–parasite interaction, such as general detoxification and bile acid sequestration to synthase activity. The GSTs have been highlighted as vaccine candidates towards parasitic flukes. In this study, we have thoroughly examined the urea-induced unfolding of a mu-class Fasciola gigantica GST1 ( FgGST1) using spectroscopic techniques and molecular dynamic simulations. FgGST1 is a highly cooperative molecule, because during urea-induced equilibrium unfolding, a concurrent unfolding of the protein without stabilization of any folded intermediate was observed. The protein was stabilized with conformational free energy of about ~12.36 kcal/mol. The protein loses its activity with increasing urea concentration, as the GSH molecule is not able to bind to the protein. We also studied the fluorescence quenching of Trp residues and the obtained K _SV data that provided additional information on the unfolding of FgGST1. Molecular dynamic trajectories simulated in different urea concentrations and temperatures indicated that urea destabilizes FgGST1 structure by weakening hydrophobic interactions and the hydrogen bond network. We observed a precise correlation between the in vitro and in silico studies.     

Molecular Origin of the Stability Difference in Four Shark IgNAR Constant Domains

Improving the stability of antibodies for manufacture and shelf life is one of the main focuses of antibody engineering. One stabilization strategy is to perform specific mutations in human antibodies based on highly stable antibodies in other species. To identify the key residues for mutagenesis, it is necessary to understand the roles of these residues in stabilizing the antibody. Here, we use molecular dynamics simulations to study the molecular origin of the four shark immunoglobulin new antigen receptors constant domains (C1–C4). According to the unfolding pathways and the conformational free energy surfaces in 8 M urea at 380 K, the C2 domain is the most stable, followed by C4, C1, and C3, which agrees with the experimental findings. The C1 and C3 domains follow a common unfolding pathway in which the unfolding starts from the edge strands, particularly strand g, and then gradually progresses to the inner strands. Detailed structural analysis of the C2 domain reveals a “sandwich-like” R339-E322-R341 salt-bridge cluster on strand g, which grants ultrahigh stability to the C2 domain. We further design two sets of mutations by mutating E322 to alanine or setting all atomic charges in E322 to zero to break the salt-bridge cluster in the C2 domain, which confirms the importance of the salt bridges in stability. In the C4 domain, the D80-K104 salt bridge on strand g also strengthens the stability. On the other hand, in the C1 and C3 domains, there is no salt bridge on strand g. In addition to the salt bridges, the overall hydrophobicity score of the hydrophobic core is also positively correlated with the domain stability. Our findings provide a detailed microscopic picture of the molecular origin of the four shark immunoglobulin new antigen receptors constant domains that not only explains the differences in their structural stability but also provides important insights into future antibody design.     

Differential backbone dynamics of companion helices in the extended helical coiled‐coil domain of a bacterial chemoreceptor

Cytoplasmic domains of transmembrane bacterial chemoreceptors are largely extended four‐helix coiled coils. Previous observations suggested the domain was structurally dynamic. We probed directly backbone dynamics of this domain of the transmembrane chemoreceptor Tar from Escherichia coli using site‐directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. Spin labels were positioned on solvent‐exposed helical faces because EPR spectra for such positions reflect primarily polypeptide backbone movements. We acquired spectra for spin‐labeled, intact receptor homodimers solubilized in detergent or inserted into native E. coli lipid bilayers in Nanodiscs, characterizing 16 positions distributed throughout the cytoplasmic domain and on both helices of its helical hairpins, one amino terminal to the membrane‐distal tight turn (N‐helix), and the other carboxyl terminal (C‐helix). Detergent solubilization increased backbone dynamics for much of the domain, suggesting that loss of receptor activities upon solubilization reflects wide‐spread destabilization. For receptors in either condition, we observed an unanticipated difference between the N‐ and C‐helices. For bilayer‐inserted receptors, EPR spectra from sites in the membrane‐distal protein‐interaction region and throughout the C‐helix were typical of well‐structured helices. In contrast, for approximately two‐thirds of the N‐helix, from its origin as the AS‐2 helix of the membrane‐proximal HAMP domain to the beginning of the membrane‐distal protein‐interaction region, spectra had a significantly mobile component, estimated by spectral deconvolution to average approximately 15%. Differential helical dynamics suggests a four‐helix bundle organization with a pair of core scaffold helices and two more dynamic partner helices. This newly observed feature of chemoreceptor structure could be involved in receptor function.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

Forced Unfolding of Apocytochrome <Emphasis Type="Italic">b</Emphasis><Subscript>5</Subscript> by Steered Molecular Dynamics Simulation

Apocytochrome b5 (apocyt b5), a small b-type cytochrome with heme prosthetic group removal, has been subjected to steered molecular dynamics (SMD) simulations for investigating the consequences of mechanical force-induced unfolding. Both constant velocity (0.5 and 1.0 A/ps) and constant force (500, 750 and 1000 pN) stretching have been employed to model forced unfolding of apocyt b5. The results of SMD simulations elucidate that apocyt b5 is protected against external stress mainly through the interstrand hydrogen bonding between its beta1-beta2 and beta2-beta3 strands, highlighting the importance of hydrophobic core 2 in stabilization of apocyt b5. The existence of intermediate states manifested by current simulations in the forced unfolding pathway of apocyt b5 is different from the observations in pervious thermal or chemical unfolding studies in the absence of force. The present study could thus provide insights into the relationship between the two cooperative functional modules of apocyt b5 and also guide the rational molecular design of heme proteins.     

Solution structures and backbone dynamics of the ribosomal protein S6 and its permutant P54‐55

The ribosomal protein S6 from Thermus thermophilus has served as a model system for the study of protein folding, especially for understanding the effects of circular permutations of secondary structure elements. This study presents the structure of a permutant protein, the 96‐residue P^54‐55, and the structure of its 101‐residue parent protein S6^wt in solution. The data also characterizes the effects of circular permutation on the backbone dynamics of S6. Consistent with crystallographic data on S6^wt, the overall solution structures of both P^54‐55 and S6^wt show a β‐sheet of four antiparallel β‐strands with two α‐helices packed on one side of the sheet. In clear contrast to the crystal data, however, the solution structure of S6^wt reveals a disordered loop in the region between β‐strands 2 and 3 (Leu43‐Phe60) instead of a well‐ordered stretch and associated hydrophobic mini‐core observed in the crystal structure. Moreover, the data for P^54‐55 show that the joined wild‐type N‐ and C‐terminals form a dynamically robust stretch with a hairpin structure that complies with the in silico design. Taken together, the results explain why the loop region of the S6^wt structure is relatively insensitive to mutational perturbations, and why P^54‐55 is more stable than S6^wt: the permutant incision at Lys54‐Asp55 is energetically neutral by being located in an already disordered loop whereas the new hairpin between the wild‐type N‐ and C‐termini is stabilizing.     

Sequential unfolding of the hemolysin two‐partner secretion domain from Proteus mirabilis

Protein secretion is a major contributor to Gram‐negative bacterial virulence. Type Vb or two‐partner secretion (TPS) pathways utilize a membrane bound β‐barrel B component (TpsB) to translocate large and predominantly virulent exoproteins (TpsA) through a nucleotide independent mechanism. We focused our studies on a truncated TpsA member termed hemolysin A (HpmA265), a structurally and functionally characterized TPS domain from Proteus mirabilis. Contrary to the expectation that the TPS domain of HpmA265 would denature in a single cooperative transition, we found that the unfolding follows a sequential model with three distinct transitions linking four states. The solvent inaccessible core of HpmA265 can be divided into two different regions. The C‐proximal region contains nonpolar residues and forms a prototypical hydrophobic core as found in globular proteins. The N‐proximal region of the solvent inaccessible core, however, contains polar residues. To understand the contributions of the hydrophobic and polar interiors to overall TPS domain stability, we conducted unfolding studies on HpmA265 and site‐specific mutants of HpmA265. By correlating the effect of individual site‐specific mutations with the sequential unfolding results we were able to divide the HpmA265 TPS domain into polar core, nonpolar core, and C‐terminal subdomains. Moreover, the unfolding studies provide quantitative evidence that the folding free energy for the polar core subdomain is more favorable than for the nonpolar core and C‐terminal subdomains. This study implicates the hydrogen bonds shared among these conserved internal residues as a primary means for stabilizing the N‐proximal polar core subdomain.     

Molecular dynamics simulations elucidate the mode of protein recognition by Skp1 and the F‐box domain in the SCF complex

Polyubiquitination of the target protein by a ubiquitin transferring machinery is key to various cellular processes. E3 ligase Skp1‐Cul1‐F‐box (SCF) is one such complex which plays crucial role in substrate recognition and transfer of the ubiquitin molecule. Previous computational studies have focused on S‐phase kinase‐associated protein 2 (Skp2), cullin, and RING‐finger proteins of this complex, but the roles of the adapter protein Skp1 and F‐box domain of Skp2 have not been determined. Using sub‐microsecond molecular dynamics simulations of full‐length Skp1, unbound Skp2, Skp2‐Cks1 (Cks1: Cyclin‐dependent kinases regulatory subunit 1), Skp1‐Skp2, and Skp1‐Skp2‐Cks1 complexes, we have elucidated the function of Skp1 and the F‐box domain of Skp2. We found that the L₁₆ loop of Skp1_, which was deleted in previous X‐ray crystallography studies, can offer additional stability to the ternary complex via its interactions with the C‐terminal tail of Skp2. Moreover, Skp1 helices H6, H7, and H8 display vivid conformational flexibility when not bound to Skp2, suggesting that these helices can recognize and lock the F‐box proteins. Furthermore, we observed that the F‐box domain could rotate (5°–129°), and that the binding partner determined the degree of conformational flexibility. Finally, Skp1 and Skp2 were found to execute a domain motion in Skp1‐Skp2 and Skp1‐Skp2‐Cks1 complexes that could decrease the distance between ubiquitination site of the substrate and the ubiquitin molecule by 3 nm. Thus, we propose that both the F‐box domain of Skp2 and Skp1‐Skp2 domain motions displaying preferential conformational control can together facilitate polyubiquitination of a wide variety of substrates. Proteins 2016; 84:159–171. © 2015 Wiley Periodicals, Inc.     

Differential stability of the bovine prion protein upon urea unfolding

Prion diseases, or transmissible spongiform encephalopathies, are a group of infectious neurological diseases associated with the structural conversion of an endogenous protein (PrP) in the central nervous system. There are two major forms of this protein: the native and noninfectious cellular form, PrP^C; and the misfolded, infectious, and proteinase K‐resistant form, PrP^Sc. The C‐terminal domain of PrP^C is mainly α‐helical in structure, whereas PrP^Sc in known to aggregate into an assembly of β‐sheets, forming amyloid fibrils. To identify the regions of PrP^C potentially involved in the initial steps of the conversion to the infectious conformation, we have used high‐resolution NMR spectroscopy to characterize the stability and structure of bovine recombinant PrP^C (residues 121 to 230) during unfolding with the denaturant urea. Analysis of the 800 MHz ¹H NMR spectra reveals region‐specific information about the structural changes occurring upon unfolding. Our data suggest that the dissociation of the native β‐sheet of PrP^C is a primary step in the urea‐induced unfolding process, while strong hydrophobic interactions between helices α1 and α3, and between α2 and α3, stabilize these regions even at very high concentrations of urea.     

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.     

The stability of cylindrin β‐barrel amyloid oligomer models—A molecular dynamics study

Small‐soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out‐off register antiparallel β‐strands in the shape of a β‐barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β‐barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out‐off‐register interface) in close proximity, which in turn stabilizes the side‐chain and main‐chain hydrogen bonds, and the salt bridges on the outer surface along the weak out‐of‐register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β‐barrel.Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Antiapoptotic Bcl‐2 homolog CED‐9 in Caenorhabditis elegans: Dynamics of BH3 and CED‐4 binding regions and comparison with mammalian antiapoptotic Bcl‐2 proteins

Proteins belonging to Bcl‐2 family regulate intrinsic cell death pathway. Although mammalian antiapoptotic Bcl‐2 members interact with multiple proapoptotic proteins, the Caenorhabditis elegans Bcl‐2 homolog CED‐9 is known to have only two proapoptotic partners. The BH3‐motif of proapoptotic proteins bind to the hydrophobic groove of prosurvival proteins formed by the Bcl‐2 helical fold. CED‐9 is also known to interact with CED‐4, a homolog of the human cell death activator Apaf1. We have performed molecular dynamics simulations of CED‐9 in two forms and compared the results with those of mammalian counterparts Bcl‐X_L, Bcl‐w, and Bcl‐2. Our studies demonstrate that the region forming the hydrophobic cleft is more flexible compared with the CED‐4‐binding region, and this is generally true for all antiapoptotic Bcl‐2 proteins studied. CED‐9 is the most stable protein during simulations and its hydrophobic pocket is relatively rigid explaining the absence of functional redundancy in CED‐9. The BH3‐binding region of Bcl‐2 is less flexible among the mammalian proteins and this lends support to the studies that Bcl‐2 binds to less number of BH3 peptides with high affinity. The C‐terminal helix of CED‐9 lost its helical character because of a large number of charged residues. We speculate that this region probably plays a role in intracellular localization of CED‐9. The BH4‐motif accessibility in CED‐9 and Bcl‐w is controlled by the loop connecting the first two helices. Although CED‐9 adopts the same Bcl‐2 fold, our studies highlight important differences in the dynamic behavior of CED‐9 and mammalian antiapoptotic homologs. Proteins 2014; 82:1035–1047. © 2013 Wiley Periodicals, Inc.     

Crystal structure of a super leucine zipper,an extended two‐stranded super long coiled coil

Coiled coil is a ubiquitous structural motif in proteins, with two to seven alpha helices coiled together like the strands of a rope, and coiled coil folding and assembly is not completely understood. A GCN4 leucine zipper mutant with four mutations of K3A, D7A, Y17W, and H18N has been designed, and the crystal structure has been determined at 1.6 Å resolution. The peptide monomer shows a helix trunk with short curved N‐ and C‐termini. In the crystal, two monomers cross in 35° and form an X‐shaped dimer, and each X‐shaped dimer is welded into the next one through sticky hydrophobic ends, thus forming an extended two‐stranded, parallel, super long coiled coil rather than a discrete, two‐helix coiled coil of the wild‐type GCN4 leucine zipper. Leucine residues appear at every seventh position in the super long coiled coil, suggesting that it is an extended super leucine zipper. Compared to the wild‐type leucine zipper, the N‐terminus of the mutant has a dramatic conformational change and the C‐terminus has one more residue Glu 32 determined. The mutant X‐shaped dimer has a large crossing angle of 35° instead of 18° in the wild‐type dimer. The results show a novel assembly mode and oligomeric state of coiled coil, and demonstrate that mutations may affect folding and assembly of the overall coiled coil. Analysis of the formation mechanism of the super long coiled coil may help understand and design self‐assembling protein fibers.     

Hierarchy of structure loss in MD simulations of src SH3 domain unfolding.

To complement experimental studies of the src SH3 domain folding, we studied 30 independent, high-temperature, molecular dynamics simulations of src SH3 domain unfolding. These trajectories were observed to differ widely from each other. Thus, rather than analyzing individual trajectories, we sought to identify the recurrent features of the high-temperature unfolding process. The conformations from all simulations were combined and then divided into groups based on the number of native contacts. Average occupancies of each side-chain hydrophobic contact and hydrogen bond in the protein were then determined. In the symmetric funnel limit, the occupancies of all contacts should decrease in concert with the loss in total number of native contacts. If there is a lack of symmetry or hierarchy to the unfolding process, the occupancies of some contacts should decrease more slowly, and others more rapidly. Despite the heterogeneity of the individual trajectories, the ensemble averaging revealed an order to the unfolding process: contacts between the N and C-terminal strands are the first to disappear, whereas contacts within the distal beta-hairpin and a hydrogen-bonding network involving the distal loop beta-turn and the diverging turn persist well after the majority of the native contacts are lost. This hierarchy of events resembles but is somewhat less pronounced than that observed in our experimental studies of the folding of src SH3 domain.     

A rapid procedure to isolate isotopically labeled peptides for NMR studies: application to the Disabled‐2 sulfatide‐binding motif

A procedure for obtaining isotopically labeled peptides, by combining affinity chromatography, urea‐equilibrated gel filtration, and hydrophobic chromatography procedures, is presented using the Disabled‐2 (Dab2) sulfatide‐binding motif (SBM) as a proof of concept. The protocol is designed to isolate unstructured, membrane‐binding, recombinant peptides that co‐purify with bacterial proteins (e.g., chaperones). Dab2 SBM is overexpressed in bacteria as an isotopically labeled glutathione S‐transferase (GST) fusion protein using minimal media containing [¹⁵N] ammonium chloride as the nitrogen source. The fusion protein is purified using glutathione beads, and Dab2 SBM is released from GST using a specific protease. It is then dried, resuspended in urea to release the bound bacterial protein, and subjected to urea‐equilibrated gel filtration. Urea and buffer reagents are removed using an octadecyl column. The peptide is eluted with acetonitrile, dried, and stored at ?80 °C. Purification of Dab2 SBM can be accomplished in 6 days with a yield of ~2 mg/l of culture. The properties of Dab2 SBM can be studied in the presence of detergents using NMR spectroscopy. Although this method also allows for the purification of unlabeled peptides that co‐purify with bacterial proteins, the procedure is more relevant to isotopically labeled peptides, thus alleviating the cost of peptide production. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

The key event in force-induced unfolding of Titin's immunoglobulin domains

Steered molecular dynamics simulation of force-induced titin immunoglobulin domain I27 unfolding led to the discovery of a significant potential energy barrier at an extension of approximately 14 A on the unfolding pathway that protects the domain against stretching. Previous simulations showed that this barrier is due to the concurrent breaking of six interstrand hydrogen bonds (H-bonds) between beta-strands A' and G that is preceded by the breaking of two to three hydrogen bonds between strands A and B, the latter leading to an unfolding intermediate. The simulation results are supported by Angstrom-resolution atomic force microscopy data. Here we perform a structural and energetic analysis of the H-bonds breaking. It is confirmed that H-bonds between strands A and B break rapidly. However, the breaking of the H-bond between strands A' and G needs to be assisted by fluctuations of water molecules. In nanosecond simulations, water molecules are found to repeatedly interact with the protein backbone atoms, weakening individual interstrand H-bonds until all six A'-G H-bonds break simultaneously under the influence of external stretching forces. Only when those bonds are broken can the generic unfolding take place, which involves hydrophobic interactions of the protein core and exerts weaker resistance against stretching than the key event.