Backbone Hydrogen Bonds Research Articles

Mutation of Charged Residues to Neutral Ones Accelerates Urea Denaturation of HP-35

Following the studies of urea denaturation of β-hairpins using molecular dynamics, in this paper, molecular dynamics simulations of two peptides, a 35 residue three helix bundle villin headpiece protein HP-35 and its doubly norleucine-substituent mutant (Lys24Nle/Lys29Nle) HP-35 NleNle, were undertaken in urea solutions to understand the molecular mechanism of urea denaturation of α-helices. The mutant HP-35 NleNle was found to denature more easily than the wild type. During the expansion of the small hydrophobic core, water penetration occurs first, followed by that of urea molecules. It was also found that the initial hydration of the peptide backbone is achieved through water hydrogen bonding with the backbone CO groups during the denaturation of both polypeptides. The mutation of the two charged lysine residues to apolar norleucine enhances the accumulation of urea near the hydrophobic core and facilitates the denaturation process. Urea also interacts directly with the peptide backbone as well as side chains, thereby stabilizing nonnative conformations. The mechanism revealed here is consistent with the previous study on secondary structure of β-hairpin polypeptide, GB1, PEPTIDE 1, and TRPZIP4, suggesting that there is a general mechanism in the denaturation of protein backbone hydrogen bonds by urea.

Molecular Basis for the Structural Stability of an Enclosed β-Barrel Loop

We present molecular dynamics simulation studies of the structural stability of an enclosed loop in the β domain of the Escherichia coli O157:H7 autotransporter EspP. Our investigation revealed that, in addition to its excellent resistance to thermal perturbations, EspP loop 5 (L5) also has remarkable mechanical stability against pulling forces along the membrane norm. These findings are consistent with the experimental report that EspP L5 helps to maintain the permeability barrier in the outer membrane. In contrast to the major secondary structure elements of globular proteins such as ubiquitin, whose resistance to thermal and mechanical perturbations depends mainly on backbone hydrogen bonds and hydrophobic interactions, the structural stability of EspP L5 can be attributed mainly to geometric constraints and side-chain interactions dominated by hydrogen bonds. Examination of B-factors from available high-resolution structures of membrane-embedded β barrels indicates that most of the enclosed loops have stable structures. This finding suggests that loops stabilized by geometric constraints and side-chain interactions might be used more generally to restrict β-barrel channels for various functional purposes.

Protein β-Sheet Nucleation Is Driven by Local Modular Formation

Despite its central role in the protein folding process, the specific mechanism(s) behind beta-sheet formation has yet to be determined. For example, whether the nucleation of beta-sheets, often containing strands separated in sequence by many residues, is local or not remains hotly debated. Here, we investigate the initial nucleation step of beta-sheet formation by performing an analysis of the smallest beta-sheets in a non-redundant dataset on the grounds that the smallest sheets, having undergone little growth after nucleation, will be enriched for nucleating characteristics. We find that the residue propensities are similar for small and large beta-sheets as are their interstrand pairing preferences, suggesting that nucleation is not primarily driven by specific residues or interacting pairs. Instead, an examination of the structural environments of the two-stranded sheets shows that virtually all of them are contained in single, compact structural modules, or when multiple modules are present, one or both of the chain termini are involved. We, therefore, find that beta-nucleation is a local phenomenon resulting either from sequential or topological proximity. We propose that beta-nucleation is a result of two opposite factors; that is, the relative rigidity of an associated folding module that holds two stretches of coil close together in topology coupled with sufficient chain flexibility that enables the stretches of coil to bring their backbones in close proximity. Our findings lend support to the hydrophobic zipper model of protein folding (Dill, K. A., Fiebig, K. M., and Chan, H. S. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 1942-1946). Implications for protein folding are discussed.

Probing the Strand Orientation and Registry Alignment in the Propagation of Amyloid Fibrils

Detailed knowledge of the structure and growth mechanism of amyloid fibrils is important for understanding the disease process. Recently, solid-state NMR and other spectroscopic data have revealed the equilibrium organization of the tertiary structure of fibrils formed by various segments of beta-amyloid peptides. A three-step "dock-and-lock" mechanism for fibril growth has been proposed on the basis of the kinetic data. Here we use all-atom replica-exchange molecular dynamics simulations in generalized-Born implicit solvent to probe the mechanism of tertiary structure propagation in fibrils of Abeta(16-22) modeled as an oligomer consisting of two beta-sheets each having four strands. The data show that following association with the oligomer, but before being fully locked onto the existing beta-sheet, the added monomer predominantly samples states with the antiparallel strand orientation, but both in- and one-residue shifted backbone hydrogen bond alignments. The in-register state, which is the experimentally observed equilibrium alignment, is marginally more stable than the registry-shifted one. These results suggest that, following the fast docking step, the added monomer dynamically slides in the backbone registry, and stabilization of the preferential alignment must occur in the second locking step as the monomer becomes fully integrated with the fibril. We also delineate the electrostatic and hydrophobic effects in directing the registry alignment during monomer addition. Surprisingly, the in-register alignment provides both increased cross-strand electrostatic attraction and hydrophobic surface burial. Finally, our data support the notion that side chain hydrophobic burial is a major driving force for beta-sheet assembly.

Golden Rule for Buttressing Vulnerable Soluble Proteins

Local weaknesses in the structure of soluble proteins have received little attention. The structure may be inherently weak at sites where hydration of the protein backbone is locally hampered by formation of an intramolecular hydrogen bond which in turn is not fully stabilized through burial within a hydrophobic environment. The result is insufficient compensation for the thermodynamic cost of dehydrating the backbone polar groups. This work shows that these structural deficiencies, the unburied backbone hydrogen bonds, are compensated in natural proteins by disulfide bonds that are needed to maintain the structural integrity. Examination of all PDB-reported soluble structures reveals that, after suitable normalization, the number of disulfide bonds, X, correlates tightly with the number of unburied backbone hydrogen bonds, Y, beyond the baseline level Y = 20, revealing a simple balance relation: Y = 5X + 20. This equation introduces a 1:5 ratio associated with the buttressing of soluble proteins with structural deficiencies. The results are justified on thermodynamic grounds and have implications for biomolecular engineering as they introduce two constants of universal applicability determining the architecture of soluble proteins.

Zipper-Like Unfolding of β-Sheets Accessed by Pioneer Water Molecules: Atomic Resolution of Forced Unfold Reveals Different Mechanisms for Parallel and Antiparallel Motifs

In this study, quantum mechanical calculations were used for an atomic level investigation of the β-sheet unfolding mechanism aided by pioneer water molecules accessing the structural motif. Results indicate that there is a qualitatively different forced unfold mechanism for parallel and antiparallel β-sheets. In the case of parallel β-sheets, the presence of only a single water molecule could already be enough to stimulate rupture of consecutive backbone hydrogen bonds by stepping from one H-bond to the next one, similarly as a slider opens up a zipper. The extension curves and energetics obtained at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level of theory correlate well and may explain the hyperfine resolution of experimentally observed sawtooth patterns in single molecule studies where external pulling force was applied.

Computational Backbone Mutagenesis of Aβ Peptides: Probing the Role of Backbone Hydrogen Bonds in Aggregation

Using replica exchange molecular dynamics (REMD) and united atom implicit solvent model we examine the role of backbone hydrogen bonds (HBs) in Abeta aggregation. The importance of HBs appears to depend on the aggregation stage. The backbone HBs have little effect on the stability of Abeta dimers or on their aggregation interface. The HBs also do not play a critical role in initial binding of Abeta peptides to the amyloid fibril. Their elimination does not change the continuous character of Abeta binding nor its temperature. However, cancellation of HBs forming between incoming Abeta peptides and the fibril disrupts the locked fibril-like states in the bound peptides. Without the support of HBs, bound Abeta peptides form few long beta-strands on the fibril edge. As a result, the deletion of peptide-fibril HBs is expected to impede fibril growth. As for the peptides bound to Abeta fibril the deletion of interpeptide HBs reduces the beta propensity in the dimers making them less competent for amyloid assembly. These simulation findings together with the backbone mutagenesis experiments suggest that a viable strategy for arresting fibril growth is the disruption of interpeptide HBs.

Hydrophobic Side Chain Dynamics of a Glutamate Receptor Ligand Binding Domain

Ionotropic glutamate receptors are ligand-gated ion channels that mediate much of the fast excitatory neurotransmission in the central nervous system. The extracellular ligand binding core (S1S2) of the GluR2 subtype of ionotropic glutamate receptors can be produced as a soluble protein with properties essentially identical to the corresponding domain in the intact, membrane-bound protein. Using a variety of biophysical techniques, much has been learned about the structure and dynamics of S1S2 and the relationship between its ligand-induced conformational changes and the function of the receptor. It is clear that dynamic processes are essential to the function of ionotropic glutamate receptors. We have isotopically labeled side chain methyls of GluR2 S1S2 and used NMR spectroscopy to study their dynamics on the ps-ns and mus-ms time scales. Increased disorder is seen in regions that are part of the key dimer interface in the intact protein. When glutamate is bound, the degree of ps-ns motion is less than that observed with other ligands, suggesting that the physiological agonist binds to a preformed binding site. At the slower time scales, the degree of S1S2 flexibility induced by ligand binding is greatest for willardiine partial agonists, least for antagonists, and intermediate for full agonists. Notable differences among bound ligands are in the region of the protein that forms a hinge between two lobes that close upon agonist binding, and along the beta-sheet in Lobe 2. These motions provide clues as to the functional properties of partial agonists and to the conformational changes associated with lobe closure and channel activation.

Synthesis and structure–activity relationship studies of peptidomimetic PKB/Akt inhibitors: The significance of backbone interactions

Elevated levels of activated Protein Kinase B (PKB/Akt) have been detected in many types of human cancer. In contrast to ATP site inhibitors, substrate-based inhibitors are more likely to be selective because of extensive interactions with the specific substrate binding site. Unfortunately, peptide-based inhibitors lack important pharmacological properties that are required of drug candidates. Chemical modifications of potent peptide inhibitors, such as peptoids and N α-methylated amino acids, may overcome these drawbacks, while maintaining potency. We present a structure-activity relationship study of a potent, peptide-based PKB/Akt inhibitor, PTR6154. The study was designed to evaluate backbone modifications on the inhibitory activity of PTR6154. Two peptidomimetic libraries, peptoid and N α-methylation, based on PTR6154, were synthesized and evaluated for in vitro PKB/Akt inhibition efficiency. All the peptoid analogs reduced potency significantly, as well as most of the members of the N-methyl library, suggesting that the backbone conformation and/or hydrogen bond interactions of PTR6154 derivatives are essential for inhibition activity. Two N-terminal members of the N-methyl library did not decrease potency and can be used as future drug leads.

Side-Chain χ1 Conformations in Urea-Denatured Ubiquitin and Protein G from 3J Coupling Constants and Residual Dipolar Couplings

Current NMR information on side-chain conformations of unfolded protein states is sparse due to the poor dispersion particularly of side-chain proton resonances. We present here optimized schemes for the detection of (3)J(HalphaHbeta), (3)J(NHbeta), and (3)J(C'Hbeta) scalar and (1)D(CbetaHbeta) residual dipolar couplings (RDCs) in unfolded proteins. For urea-denatured ubiquitin and protein G, up to six (3)J-couplings to (1)H(beta) are detected, which define the chi(1) angle at very high precision. Interpretation of the (3)J couplings by a model of mixed staggered chi(1) rotamers yields excellent agreement and also provides stereoassignments for (1)H(beta) methylene protons. For all observed amino acids with the exception of leucine, the chemical shift of (1)H(beta3) protons was found downfield from (1)H(beta2). For most residues, the precision of individual chi(1) rotamer populations is better than 2%. The experimental chi(1) rotamer populations are in the vicinity of averages obtained from coil regions in folded protein structures. However, individual variations from these averages of up to 40% are highly significant and indicate sequence- and residue-specific interactions. Particularly strong deviations from the coil average are found for serine and threonine residues, an effect that may be explained by a weakening of side-chain to backbone hydrogen bonds in the urea-denatured state. The measured (1)D(CbetaHbeta) RDCs correlate well with predicted RDCs that were calculated from a sterically aligned coil model ensemble and the (3)J-derived chi(1) rotamer populations. This agreement supports the coil model as a good first approximation of the unfolded state. Deviations between measured and predicted values at certain sequence locations indicate that the description of the local backbone conformations can be improved by incorporation of the RDC information. The ease of detection of a large number of highly precise side-chain RDCs opens the possibility for a more rigorous characterization of both side-chain and backbone conformations in unfolded proteins.

Extended Morse Function Model for Angle-Dependent Hydrogen Bond in Protein−Protein Interactions

Backbone-backbone, backbone-asparagine, and serine-backbone hydrogen bonds (HBs) are the most abundant interactions at the interface of protein-protein complex. Here, we propose an angle-dependent potential energy function for these HBs constructed by the product of the radial and the angular Morse functions whose various parameters are optimized with high-level density functional theory (DFT) calculations. The new angular variables, the interatomic distance between the donor and the acceptor atoms (R(theta)) and that between the hydrogen and the base atom of the acceptor (R(phi)), are employed to define the angular Morse functions. The angular part in the new potential function is found to be comparable in the magnitude of energy values to the radial one, which is consistent with the significant angular dependence of HBs. The HB binding energies calculated with the new potential function compare well with those obtained by high-level DFT calculations with the associated squared correlation coefficients ranging from 0.82 to 0.85. This agreement indicates the suitability of the new energy functions as a potential function for HB in modeling the protein-protein interactions.

A Molecular Dynamics Study of Human Defensins HBD-1 and HNP-3 in Water

Mammalian defensins are crucial components of the innate immune system. They are characterized by three disulfide bridges and exhibit broad spectrum antibacterial activity. The spacing between the cysteines and disulfide connectivities in the two classes of defensins, the α- and β-forms, are different. The structural motif of 3 β-strands appears to be conserved in α and β-defensins despite differences in disulfide connectivities and spacing between cysteines. In this study, Molecular Dynamics Simulations (MDS) have been carried out to study the conformational behavior of α- and β-defensins with and without disulfide bridges. Our results indicate that β-strands in the C-terminal region of HBD-1 and HNP-3 do not unfold during the course of MDS. The segment adopting α-helix in HBD-1 unfolds early during the simulations. The backbone hydrogen bonds in HBD-1 and HNP-3 are broken during MDS. When the disulfide bonds are absent, the N-terminal β-strand unfolds by 20 ns but β-strands are observed in the C-terminal region of HNP-3. HBD-1, without disulfide bridges, unfolds to a greater extent during the course of the MDS. Examination of distances between sulfur atoms of cysteines without disulfide bridges during the simulations indicate that there is no specific preference for native disulfide bridges, which could be the reason for the experimental observation of non-native disulfide bridge formation during chemical synthesis of human α- and β-defensins. Since defensins with non-native disulfide bridges are biologically active, the exact three dimensional structures observed for native HBD-1 and HNP-3 does not appear to be essential for exhibiting antibacterial activity.

Mining protein dynamics from sets of crystal structures using “consensus structures”

In this work, we describe two novel approaches to utilize the dynamic structure information implicitly contained in large crystal structure data sets. The first approach visualizes both consistent as well as variable ligand-induced changes in ligand-bound compared with apo protein crystal structures. For this purpose, information was mined from B-factors and ligand-induced residue displacements in multiple crystal structures, minimizing experimental error and noise. With this approach, the mechanism of action of non-nucleoside reverse transcriptase inhibitors (NNRTIs) as an inseparable combination of distortion of protein dynamics and conformational changes of HIV-1 reverse transcriptase was corroborated (a combination of the previously proposed "molecular arthritis" and "distorted site" mechanisms). The second approach presented here uses "consensus structures" to map common binding features that are present in a set of structures of NNRTI-bound HIV-1 reverse transcriptase. Consensus structures are based on different levels of structural overlap of multiple crystal structures and are used to analyze protein-ligand interactions. The structures are shown to yield information about conserved hydrogen bonding interactions as well as binding-pocket flexibility, shape, and volume. From the consensus structures, a common wild type NNRTI binding pocket emerges. Furthermore, we were able to identify a conserved backbone hydrogen bond acceptor at P236 and a novel hydrophobic subpocket, which are not yet utilized by current drugs. Our methods introduced here reinterpret the atom information and make use of the data variability by using multiple structures, complementing classical 3D structural information of single structures.

Structure of the N-terminal Regulatory Domain of a Plant NADPH Oxidase and Its Functional Implications

Plant NADPH oxidases (Rboh, for respiratory burst oxidase homolog) produce reactive oxygen species that are key regulators of various cellular events including plant innate immunity. Rbohs possess a highly conserved cytoplasmic N-terminal region containing two EF-hand motifs that regulate Rboh activity. Rice (Oryza sativa) RbohB (OsRbohB) is regulated by the direct binding of a small GTPase (Rac1) to this regulatory region as well as by Ca(2+) binding to the EF-hands. Here, we present the atomic structure of the N-terminal region of OsRbohB. The structure reveals that OsRbohB forms a unique dimer stabilized by swapping the EF-hand motifs. We identified two additional EF-hand-like motifs that were not predicted from sequence data so far. These EF-hand-like motifs together with the swapped EF-hands form a structure similar to that found in calcineurin B. We observed conformational changes mediated by Ca(2+) binding to only one EF-hand. Structure-based in vitro pulldown assays and NMR titration experiments defined the OsRac1 binding interface within the coiled-coil region created by swapping the EF-hands. In addition, we demonstrate a direct intramolecular interaction between the N and C terminus, and that the complete N-terminal cytoplasmic region is required for this interaction. The structural features and intramolecular interactions characterized here might be common elements shared by Rbohs that contribute to the regulation of reactive oxygen species production.

Importance of Protein Context on the Unfolding Pathways of β-hairpins

Small β-hairpin peptides as well as three-stranded WW domains have been used as models for the folding of β-hairpins in larger proteins. Previous studies of the FBP28 WW domain proposed that side chain contacts between residues in the strands and not the precise order of backbone hydrogen bond formation guide β-hairpin folding. But how applicable is the folding of model systems, such as FBP28 WW domain, to the folding of β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding molecular dynamics simulations of three proteins that share a double hairpin motif structurally similar to WW domains: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain.

RNA structural motifs that entail hydrogen bonds involving sugar–phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

FoldAmyloid: a method of prediction of amyloidogenic regions from protein sequence

Amyloidogenic regions in polypeptide chains are very important because such regions are responsible for amyloid formation and aggregation. It is useful to be able to predict positions of amyloidogenic regions in protein chains. Two characteristics (expected probability of hydrogen bonds formation and expected packing density of residues) have been introduced by us to detect amyloidogenic regions in a protein sequence. We demonstrate that regions with high expected probability of the formation of backbone-backbone hydrogen bonds as well as regions with high expected packing density are mostly responsible for the formation of amyloid fibrils. Our method (FoldAmyloid) has been tested on a dataset of 407 peptides (144 amyloidogenic and 263 non-amyloidogenic peptides) and has shown good performance in predicting a peptide status: amyloidogenic or non-amyloidogenic. The prediction based on the expected packing density classified correctly 75% of amyloidogenic peptides and 74% of non-amyloidogenic ones. Two variants (averaging by donors and by acceptors) of prediction based on the probability of formation of backbone-backbone hydrogen bonds gave a comparable efficiency. With a hybrid-scale constructed by merging the above three scales, our method is correct for 80% of amyloidogenic peptides and for 72% of non-amyloidogenic ones. Prediction of amyloidogenic regions in proteins where positions of amyloidogenic regions are known from experimental data has also been done. In the proteins, our method correctly finds 10 out of 11 amyloidogenic regions. The FoldAmyloid server is available at http://antares.protres.ru/fold-amyloid/.

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.

Structure of Escherichia coli Hfq bound to polyriboadenylate RNA

Hfq is a small, highly abundant hexameric protein that is found in many bacteria and plays a critical role in mRNA expression and RNA stability. As an "RNA chaperone," Hfq binds AU-rich sequences and facilitates the trans annealing of small RNAs (sRNAs) to their target mRNAs, typically resulting in the down-regulation of gene expression. Hfq also plays a key role in bacterial RNA decay by binding tightly to polyadenylate [poly(A)] tracts. The structural mechanism by which Hfq recognizes and binds poly(A) is unknown. Here, we report the crystal structure of Escherichia coli Hfq bound to the poly(A) RNA, A(15). The structure reveals a unique RNA binding mechanism. Unlike uridine-containing sequences, which bind to the "proximal" face, the poly(A) tract binds to the "distal" face of Hfq using 6 tripartite binding motifs. Each motif consists of an adenosine specificity site (A site), which is effected by peptide backbone hydrogen bonds, a purine nucleotide selectivity site (R site), and a sequence-nondiscriminating RNA entrance/exit site (E site). The resulting implication that Hfq can bind poly(A-R-N) triplets, where R is a purine nucleotide and N is any nucleotide, was confirmed by binding studies. Indeed, Hfq bound to the oligoribonucleotides (AGG)(8), (AGC)(8), and the shorter (A-R-N)(4) sequence, AACAACAAGAAG, with nanomolar affinities. The abundance of (A-R-N)(4) and (A-R-N)(5) triplet repeats in the E. coli genome suggests additional RNA targets for Hfq. Further, the structure provides insight into Hfq-mediated sRNA-mRNA annealing and the role of Hfq in RNA decay.

Fibronectin forms the most extensible biological fibers displaying switchable force-exposed cryptic binding sites

Rather than maximizing toughness, as needed for silk and muscle titin fibers to withstand external impact, the much softer extracellular matrix fibers made from fibronectin (Fn) can be stretched by cell generated forces and display extraordinary extensibility. We show that Fn fibers can be extended more than 8-fold (>700% strain) before 50% of the fibers break. The Young's modulus of single fibers, given by the highly nonlinear slope of the stress-strain curve, changes orders of magnitude, up to MPa. Although many other materials plastically deform before they rupture, evidence is provided that the reversible breakage of force-bearing backbone hydrogen bonds enables the large strain. When tension is released, the nano-sized Fn domains first contract in the crowded environment of fibers within seconds into random coil conformations (molten globule states), before the force-bearing hydrogen bond networks that stabilize the domain's secondary structures are reestablished within minutes (double exponential). The exposure of cryptic binding sites on Fn type III modules increases steeply upon stretching. Thus fiber extension steadily up-regulates fiber rigidity and cryptic epitope exposure, both of which are known to differentially alter cell behavior. Finally, since stress-strain relationships cannot directly be measured in native extracellular matrix (ECM), the stress-strain curves were correlated with stretch-induced alterations of intramolecular fluorescence resonance energy transfer (FRET) obtained from trace amounts of Fn probes (mechanical strain sensors) that can be incorporated into native ECM. Physiological implications of the extraordinary extensibility of Fn fibers and contraction kinetics are discussed.