Structural characterization of the N-terminal mineral modification domains from the molluscan crystal-modulating biomineralization proteins, AP7 and AP24

The AP7 and AP24 proteins represent a class of mineral-interaction polypeptides that are found in the aragonite-containing nacre layer of mollusk shell (H. rufescens). These proteins have been shown to preferentially interfere with calcium carbonate mineral growth in vitro. It is believed that both proteins play an important role in aragonite polymorph selection in the mollusk shell. Previously, we demonstrated the 1-30 amino acid (AA) N-terminal sequences of AP7 and AP24 represent mineral interaction/modification domains in both proteins, as evidenced by their ability to frustrate calcium carbonate crystal growth at step edge regions. In this present report, using free N-terminal, C(alpha)-amide "capped" synthetic polypeptides representing the 1-30 AA regions of AP7 (AP7-1 polypeptide) and AP24 (AP24-1 polypeptide) and NMR spectroscopy, we confirm that both N-terminal sequences possess putative Ca (II) interaction polyanionic sequence regions (2 x -DD- in AP7-1, -DDDED- in AP24-1) that are random coil-like in structure. However, with regard to the remaining sequences regions, each polypeptide features unique structural differences. AP7-1 possesses an extended beta-strand or polyproline type II-like structure within the A11-M10, S12-V13, and S28-I27 sequence regions, with the remaining sequence regions adopting a random-coil-like structure, a trait common to other polyelectrolyte mineral-associated polypeptide sequences. Conversely, AP24-1 possesses random coil-like structure within A1-S9 and Q14-N16 sequence regions, and evidence for turn-like, bend, or loop conformation within the G10-N13, Q17-N24, and M29-F30 sequence regions, similar to the structures identified within the putative elastomeric proteins Lustrin A and sea urchin spicule matrix proteins. The similarities and differences in AP7 and AP24 N-terminal domain structure are discussed with regard to joint AP7-AP24 protein modification of calcium carbonate growth.     

Modeling AFM-induced PEVK extension and the reversible unfolding of Ig/FNIII domains in single and multiple titin molecules

Molecular elasticity is associated with a select number of polypeptides and proteins, such as titin, Lustrin A, silk fibroin, and spider silk dragline protein. In the case of titin, the globular (Ig) and non-globular (PEVK) regions act as extensible springs under stretch; however, their unfolding behavior and force extension characteristics are different. Using our time-dependent macroscopic method for simulating AFM-induced titin Ig domain unfolding and refolding, we simulate the extension and relaxation of hypothetical titin chains containing Ig domains and a PEVK region. Two different models are explored: 1) a series-linked WLC expression that treats the PEVK region as a distinct entropic spring, and 2) a summation of N single WLC expressions that simulates the extension and release of a discrete number of parallel titin chains containing constant or variable amounts of PEVK. In addition to these simulations, we also modeled the extension of a hypothetical PEVK domain using a linear Hooke's spring model to account for "enthalpic" contributions to PEVK elasticity. We find that the modified WLC simulations feature chain length compensation, Ig domain unfolding/refolding, and force-extension behavior that more closely approximate AFM, laser tweezer, and immunolocalization experimental data. In addition, our simulations reveal the following: 1) PEVK extension overlaps with the onset of Ig domain unfolding, and 2) variations in PEVK content within a titin chain ensemble lead to elastic diversity within that ensemble.     

A kinetic molecular model of the reversible unfolding and refolding of titin under force extension.

Molecular elasticity is a physicomechanical property that is associated with a select number of polypeptides and proteins, such as the giant muscle protein, titin, and the extracellular matrix protein, tenascin. Both proteins have been the subject of atomic force microscopy (AFM), laser tweezer, and other in vitro methods for examining the effects of force extension on the globular (FNIII/Ig-like) domains that comprise each protein. In this report we present a time-dependent method for simulating AFM force extension and its effect on FNIII/Ig domain unfolding and refolding. This method treats the unfolding and refolding process as a standard three-state protein folding model (U right arrow over left arrow T right arrow over left arrow F, where U is the unfolded state, T is the transition or intermediate state, and F is the fully folded state), and integrates this approach within the wormlike chain (WLC) concept. We simulated the effect of AFM tip extension on a hypothetical titin molecule comprised of 30 globular domains (Ig or FNIII) and 25% Pro-Glu-Val-Lys (PEVK) content, and analyzed the unfolding and refolding processes as a function of AFM tip extension, extension rate, and variation in PEVK content. In general, we find that the use of a three-state protein-folding kinetic-based model and the implicit inclusion of PEVK domains can accurately reproduce the experimental force-extension curves observed for both titin and tenascin proteins. Furthermore, our simulation data indicate that PEVK domains exhibit extensibility behavior, assist in the unfolding and refolding of FNIII/Ig domains in the titin molecule, and act as a force "buffer" for the FNIII/Ig domains, particularly at low and moderate extension forces.     

In vitro unfolding,refolding, and polymerization of human gammaD crystallin,a protein involved in cataract formation

Human gammaD crystallin (HgammaD-Crys), a major protein of the human eye lens, is a primary component of cataracts. This 174-residue primarily beta-sheet protein is made up of four Greek keys separated into two domains. Mutations in the human gene sequence encoding HgammaD-Crys are implicated in early-onset cataracts in children, and the mutant protein expressed in Escherichia coli exhibits properties that reflect the in vivo pathology. We have characterized the unfolding, refolding, and competing aggregation of human wild-type HgammaD-Crys as a function of guanidinium hydrochloride (GuHCl) concentration at neutral pH and 37 degrees C, using intrinsic tryptophan fluorescence to monitor in vitro folding. Wild-type HgammaD-Crys exhibited reversible refolding above 1.0 M GuHCl. The GuHCl unfolded protein was more fluorescent than its native counterpart despite the absence of metal or ion-tryptophan interactions. Aggregation of refolding intermediates of HgammaD-Crys was observed in both equilibrium and kinetic refolding processes. The aggregation pathway competed with productive refolding at denaturant concentrations below 1.0 M GuHCl, beyond the major conformational transition region. Atomic force microscopy of samples under aggregating conditions revealed the sequential appearance of small nuclei, thin protofibrils, and fiber bundles. The HgammaD-Crys fibrous aggregate species bound bisANS appreciably, indicating the presence of exposed hydrophobic pockets. The mechanism of HgammaD-Crys aggregation may provide clues to understanding age-onset cataract formation in vivo.     

Effect of single amino acid replacements on the folding and stability of dihydrofolate reductase from Escherichia coli

The role of the secondary structure in the folding mechanism of dihydrofolate reductase from Escherichia coli was probed by studying the effects of amino acid replacements in two alpha helices and two strands of the central beta sheet on the folding and stability. The effects on stability could be qualitatively understood in terms of the X-ray structure for the wild-type protein by invoking electrostatic, hydrophobic, or hydrogen-bonding interactions. Kinetic studies focused on the two slow reactions that are thought to reflect the unfolding/refolding of two stable native conformers to/from their respective folding intermediates [Touchette, N. A., Perry, K. M., & Matthews, C. R. (1986) Biochemistry 25, 5445-5452]. Replacements at three different positions in helix alpha B selectively alter the relaxation time for unfolding while a single replacement in helix alpha C selectively alters the relaxation time for refolding. This behavior is characteristic of mutations that change the stability of the protein but do not affect the rate-limiting step. In striking contrast, replacements in strands beta F and beta G can affect both unfolding and refolding relaxation times. This behavior shows that these mutations alter the rate-limiting step in these native-to-intermediate folding reactions. It is proposed that the intermediates have an incorrectly formed beta sheet whose maturation to the structure found in the native conformation is one of the slow steps in folding.     

Mechanical stability and differentially conserved physical–chemical properties of titin Ig‐domains

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single‐molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I‐band of a sarcomere is composed of about 40 Ig‐domains and is exposed to force under normal physiological conditions and connects the free‐hanging ends of the myosin filaments to the Z‐disc. Recent single‐molecule force spectroscopy data show a mechanical hierarchy in the I‐band domains. Domains near the C‐terminus in this region unfold at forces two to three times greater than domains near the beginning of the I‐band. Though all of these Ig‐domains are thought to share a fold and topology common to members of the Ig‐like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I‐band Ig domain to specific, conserved physical–chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into “strong” and “weak” families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the β‐sandwich fold, and force sensitive residues are not only confined to the A′‐G region. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Denaturation of cobalt-carbonic anhydrase B from bovine erythrocyte.

We have attempted to elucidate the mechanism of the acquisition of active-site conformation in enzymes by studying the unfolding and refolding of Co(II) carbonic anhydrase. This enzyme possesses characteristic absorption and CD spectra in the visible region, which reflect primarily the local conformational environment of the active site where the Co(II) is bound. The conformation of the cobalt–protein is found to be almost the same as the native zinc–protein. It has lower thermal stability and unfolds at a lower concentration of guanidinium chloride, even though a single-phase transition is observed in these denaturations. The unfolding to the random-coiled protein was followed by various spectroscopic methods which reflect different conformational changes. The results show that separable stages of unfolding are observable for the cobalt–protein. It appears that the collapse of the active-site conformation, as followed by changes in visible absorption and CD parameters, occurs immediately after the loosening of the molecule. This is followed by the exposure of buried aromatic residues and the unraveling of the backbone secondary structure. However, renaturation of the cobalt–protein results in an inactive protein which has spectroscopic properties significantly different from those of the native cobalt–protein.     

Folding myoglobin within a sol-gel glass: protein folding constrained to a small volume

The unfolding and refolding reaction of myoglobin was examined in solution and within a porous silica sol-gel glass. The sol-gel pores constrain the protein to a volume that is the same size and shape as the folded native state accompanied by a few layers of water solvation. Denaturants such as low pH buffers can be diffused through the gel pores to the protein to initiate unfolding and refolding. Acid-induced unfolding was hindered by the steric constraints imposed by the gel pores such that more denaturing conditions were required within the gel than in solution to create the unfolded state. No new folding intermediates were observed. Refolding of myoglobin was not complete in millimolar pH 7 buffer alone. Addition of 25% glycerol to the pH 7 buffer resulted in nearly complete refolding, and the use of 1 M phosphate buffer resulted in complete refolding. The role of this cosolvent and salt in disrupting the ordered water surrounding the protein within the gel is discussed in light of the Hofmeister series and entropic trapping via a diminished hydrophobic effect within the gel. These results are consistent with the premises of folding models in which secondary and tertiary structures are considered to form within a compact conformation of the protein backbone.     

Differential scanning calorimetric,circular dichroism,and Fourier transform infrared spectroscopic characterization of the thermal unfolding of xylanase A from Streptomyces lividans

The thermal unfolding of xylanase A from Streptomyces lividans, and of its isolated substrate binding and catalytic domains, was studied by differential scanning calorimetry and Fourier transform infrared and circular dichroism spectroscopy. Our calorimetric studies show that the thermal denaturation of the intact enzyme is a complex process consisting of two endothermic events centered near 57 and 64 degrees C and an exothermic event centered near 75 degrees C, all of which overlap slightly on the temperature scale. A comparison of the data obtained with the intact enzyme and isolated substrate binding and catalytic domains indicate that the lower- and higher-temperature endothermic events are attributable to the thermal unfolding of the xylan binding and catalytic domains, respectively, whereas the higher-temperature exothermic event arises from the aggregation and precipitation of the denatured catalytic domain. Moreover, the thermal unfolding of the two domains of the native enzyme are thermodynamically independent and differentially sensitive to pH. The unfolding of the substrate binding domain is a reversible two-state process and, under appropriate conditions, the refolding of this domain to its native conformation can occur. In contrast, the unfolding of the catalytic domain is a more complex process in which two subdomains unfold independently over a similar temperature range. Also, the unfolding of the catalytic domain leads to aggregation and precipitation, which effectively precludes the refolding of the protein to its native conformation. These observations are compatible with the results of our spectroscopic studies, which show that the catalytic and substrate binding domains of the enzyme are structurally dissimilar and that their native conformations are unaffected by their association in the intact enzyme. Thus, the calorimetric and spectroscopic data demonstrate that the S. lividans xylanase A consists of structurally dissimilar catalytic and substrate binding domains that, although covalently linked, undergo essentially independent thermal denaturation. These observations provide valuable new insights into the structure and thermal stability of this enzyme and should assist our efforts at engineering xylanases that are more thermally robust and otherwise better suited for industrial applications.     

Pressure as a denaturing agent in studies of single-point mutants of an amyloidogenic protein human cystatin c

Recently, we presented a convenient method combining a deuterium‐hydrogen exchange and electrospray mass spectrometry for studying high‐pressure denaturation of proteins (Stefanowicz et al., Biosci Rep 2009; 30:91–99). Here, we present results of pressure‐induced denaturation studies of an amyloidogenic protein—the wild‐type human cystatin C (hCC) and its single‐point mutants, in which Val57 residue from the hinge region was substituted by Asn, Asp or Pro, respectively. The place of mutation and the substituting residues were chosen mainly on a basis of theoretical calculations. Observation of H/D isotopic exchange proceeding during pressure induced unfolding and subsequent refolding allowed us to detect differences in the proteins stability and folding dynamics. On the basis of the obtained results we can conclude that proline residue at the hinge region makes cystatin C structure more flexible and dynamic, what probably facilitates the dimerization process of this hCC variant. Polar asparagine does not influence stability of hCC conformation significantly, whereas charged aspartic acid in 57 position makes the protein structure slightly more prone to unfolding. Our experiments also point out pressure denaturation as a valuable supplementary method in denaturation studies of mutated proteins. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Refolding studies on the tetrameric loop deletion mutant RM6 of ROP protein.

Previous DSC and X-ray studies on RM6, a loop deletion mutant of wtROP protein, have shown that removal of five amino acids from the loop causes a dramatic reorganization of the wild-type structure. The new tetrameric molecule exhibits a significantly higher stability (Lassalle, M.W. et al., J. Mol. Biol., 1998, 279, 987-1000) and unfolds in a second order reaction (Lassalle, M.W. and Hinz, H.-J., Biochemistry, 1998, 37, 8465-8472). In the present investigation we report extensive refolding studies of RM6 at different temperatures and GdnHCl concentrations monitored by CD and fluorescence to probe for changes in secondary and tertiary structure, respectively. The measurements permitted us to determine activation parameters as a function of denaturant concentration. The results demonstrate convincingly that the variation with GdnHCl concentration of the activation parameters deltaH#, deltaS# and deltaG# is very similar for unfolding and refolding. For both processes the activation properties approach a maximum in the vicinity of the denaturant concentration, c(K=1), where the equilibrium constant equals 1, i.e. deltaG0 equals zero. CD and fluorescence refolding kinetics are described by identical constants suggesting that the formation of secondary and tertiary structure occurs simultaneously. Refolding is, however, characterized by a more complex mechanism than unfolding. Although the general pattern is dominated by the sequence monomers to dimers to tetramers, parallel side reactions involving dimers and monomers have to be envisaged in the initial folding phase, supporting the view that the native state of RM6 can be reached by several rather than a single pathway.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long‐range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51–56 [IG(51–56)] and 50–57 [IG(50–57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the β‐turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283–323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the β‐hairpin structure of peptides corresponding to the C‐terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long‐range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the β‐hairpin structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

The study of protein mechanics with the atomic force microscope.

The unfolding and folding of single protein molecules can be studied with an atomic force microscope (AFM). Many proteins with mechanical functions contain multiple, individually folded domains with similar structures. Protein engineering techniques have enabled the construction and expression of recombinant proteins that contain multiple copies of identical domains. Thus, the AFM in combination with protein engineering has enabled the kinetic analysis of the force-induced unfolding and refolding of individual domains as well as the study of the determinants of mechanical stability.     

A kinetic study of the folding of staphylococcal nuclease using size-exclusion chromatography.

The kinetics of the hydrodynamic volume change accompanying the reversible unfolding of staphylococcal nuclease have been observed by size-exclusion chromatography at 4 degrees C and pH 7.0 using the denaturant guanidine hydrochloride. The observed chromatographic profiles have been simulated by a six-component unfolding/refolding mechanism using a consistent set of equilibrium and kinetic parameters. The native protein is an equilibrium mixture of the cis and trans isomers of the peptide bond preceding proline-117. The native conformation containing the cis isomer dominates the equilibrium mixture, is more stable, and unfolds more slowly at its transition midpoint. The denatured protein is an equilibrium mixture of at least four components, the cis/trans isomers of proline-117 and one of the five remaining prolines. The dominant refolding pathway is initiated from the denatured component containing the trans isomer of proline-117. The six-component mechanism is consistent with tryptophan fluorescence kinetic measurements of the wild-type protein and with chromatographic measurements of a mutant P117G protein.     

Distinct dynamics and interaction patterns in H‐ and K‐Ras oncogenic P‐loop mutants

Despite years of study, the structural or dynamical basis for the differential reactivity and oncogenicity of Ras isoforms and mutants remains unclear. In this study, we investigated the effects of amino acid variations on the structure and dynamics of wild type and oncogenic mutants G12D, G12V, and G13D of H‐ and K‐Ras proteins. Based on data from µs‐scale molecular dynamics simulations, we show that the overall structure of the proteins remains similar but there are important differences in dynamics and interaction networks. We identified differences in residue interaction patterns around the canonical switch and distal loop regions, and persistent sodium ion binding near the GTP particularly in the G13D mutants. Our results also suggest that different Ras variants have distinct local structural features and interactions with the GTP, variations that have the potential to affect GTP release and hydrolysis. Furthermore, we found that H‐Ras proteins and particularly the G12V and G13D variants are significantly more flexible than their K‐Ras counterparts. Finally, while most of the simulated proteins sampled the effector‐interacting state 2 conformational state, G12V and G13D H‐Ras adopted an open switch state 1 conformation that is defective in effector interaction. These differences have implications for Ras GTPase activity, effector or exchange factor binding, dimerization and membrane interaction. Proteins 2017; 85:1618–1632. © 2017 Wiley Periodicals, Inc.     

Reversible and irreversible unfolding of multi-domain proteins

In contrast to single-domain proteins unfolding of larger multi-domain proteins is often irreversible. In a comparative case study on three different multi-domain proteins (phosphoglycerate kinase: PGK and two homologous alpha-amylases: TAKA and BLA) we investigated properties of unfolded states and their ability to fold back into the native state. For this purpose guanidine hydrochloride, alkaline pH, and thermal unfolded states were characterized. Structural alterations upon unfolding and refolding transitions were monitored using fluorescence and CD spectroscopy. Static and dynamic light scattering was employed to follow aggregation processes. Furthermore, proper refolding was also investigated by enzyme activity measurements. While for PGK at least partial reversible unfolding transitions were observed in most cases, we found reversible unfolding for TAKA in the case of alkaline pH and GndHCl induced unfolding. BLA exhibits reversible unfolding only under conditions with high concentrations of protecting osmolytes (glycerol), indicating that aggregation of the unfolded state is the main obstacle to achieve proper refolding for this protein. Structural properties, such as number and size of domains, secondary structure contents and compositions within domains, and domain topology were analyzed and considered in the interpretation of differences in refolding behavior of the investigated proteins.     

Unfolding study of a trimeric membrane protein AcrB

The folding of a multi‐domain trimeric α‐helical membrane protein, Escherichia coli inner membrane protein AcrB, was investigated. AcrB contains both a transmembrane domain and a large periplasmic domain. Protein unfolding in sodium dodecyl sulfate (SDS) and urea was monitored using the intrinsic fluorescence and circular dichroism spectroscopy. The SDS denaturation curve displayed a sigmoidal profile, which could be fitted with a two‐state unfolding model. To investigate the unfolding of separate domains, a triple mutant was created, in which all three Trp residues in the transmembrane domain were replaced with Phe. The SDS unfolding profile of the mutant was comparable to that of the wild type AcrB, suggesting that the observed signal change was largely originated from the unfolding of the soluble domain. Strengthening of trimer association through the introduction of an inter‐subunit disulfide bond had little effect on the unfolding profile, suggesting that trimer dissociation was not the rate‐limiting step in unfolding monitored by fluorescence emission. Under our experimental condition, AcrB unfolding was not reversible. Furthermore, we experimented with the refolding of a monomeric mutant, AcrB_Δloop, from the SDS unfolded state. The CD spectrum of the refolded AcrB_Δloop superimposed well onto the spectra of the original folded protein, while the fluorescence spectrum was not fully recovered. In summary, our results suggested that the unfolding of the trimeric AcrB started with a local structural rearrangement. While the refolding of secondary structure in individual monomers could be achieved, the re‐association of the trimer might be the limiting factor to obtain folded wild‐type AcrB.     

Effect of an alternative disulfide bond on the structure, stability, and folding of human lysozyme

Human lysozyme has four disulfide bonds, one of which, Cys65-Cys81, is included in a long loop of the beta-domain. A cysteine-scanning mutagenesis in which the position of Cys65 was shifted within a continuous segment from positions 61 to 67, with fixed Cys81, has previously shown that only the mutant W64CC65A, which has a nonnative Cys64-Cys81 disulfide, can be correctly folded and secreted by yeast. Here, using the W64CC65A mutant, we investigated the effects of an alternative disulfide bond on the structure, stability, and folding of human lysozyme using circular dichroism (CD) and fluorescence spectroscopy combined with a stopped-flow technique. Although the mutant is expected to have a different main-chain structure from that of the wild-type protein around the loop region, far- and near-UV CD spectra show that the native state of the mutant has tightly packed side chains and secondary structure similar to that of the wild-type. Guanidine hydrochloride-induced equilibrium unfolding transition of the mutant is reversible, showing high stability and cooperativity of folding. In the kinetic folding reaction, both proteins accumulate a similar burst-phase intermediate having pronounced secondary structure within the dead time of the measurement and fold into the native structure by means of a similar folding mechanism. Both the kinetic refolding and unfolding reactions of the mutant protein are faster than those of the wild-type, but the increase in the unfolding rate is larger than that of the refolding rate. The Gibbs' free-energy diagrams obtained from the kinetic analysis suggest that the structure around the loop region in the beta-domain of human lysozyme is formed after the transition state of folding, and thus, the effect of the alternative disulfide bond on the structure, stability, and folding of human lysozyme appears mainly in the native state.     

600 ps Molecular dynamics reveals stable substructures and flexible hinge points in cAMP dependent protein kinase

Molecular dynamics simulations of the catalytic subunit of cAMP dependent protein kinase (cAPK) have been performed in an aqueous environment. The relations among the protein hydrogen‐bonding network, secondary structural elements, and the internal motions of rigid domains were examined. The values of fluctuations of protein dihedral angles during dynamics show quite distinct maxima in the regions of loops and minima in the regions of α‐helices and β‐strands. Analyses of conformation snapshots throughout the run show stable subdomains and indicate that these rigid domains are constrained during the dynamics by a stable network of hydrogen bonds. The most stable subdomain during the dynamics was in the small lobe including part of the carboxy‐terminal tail. The most significant flexible region was the highly conserved glycine‐rich loop between β strands 1 and 2 in the small lobe. Many of the main chain dihedral angle changes measured in a comparison of the crystallographic structures of “open” and “closed” conformations of cAPK correspond to the highly flexible residues found during dynamics. © 1999 John Wiley & Sons, Inc. Biopoly 50: 513–524, 1999     

Simulated refolding of stretched titin immunoglobulin domains

Steered molecular dynamics (SMD) is used to investigate forced unfolding and spontaneous refolding of immunoglobulin I27, a domain of the muscle protein titin. Previous SMD simulations revealed the events leading to stretch-induced unfolding of I27, the rupture of hydrogen bonds bridging beta-strands A and B, and those bridging beta-strands A' and G, the latter rupture occurring at an extension of approximately 15 A and preceding the complete unfolding. Simulations are now used to study the refolding of partially unfolded I27 domains. The results reveal that stretched domains with ruptured interstrand hydrogen bonds shrink along the extension direction. Two types of refolding patterns are recognized: for separated beta-strands A' and G, in most simulations five of the six hydrogen bonds between A' and G stably reformed in 2 ns, whereas for separated beta-strands A and B hydrogen bonds seldom reformed in eight 2-ns simulations. The mechanical stability of the partially refolded intermediates has been tested by re-stretching.