Comparison of the backbone dynamics of a natural and a consensus designed 3-TPR domain

The tetratricopeptide repeat (TPR) is a 34-amino acid helix-turn-helix motif that occurs in tandem arrays in numerous proteins. Here we compare the backbone dynamics of a natural 3-repeat TPR domain, from the protein UBP, with the behavior of a designed protein CTPR3, which consists of three identical consensus TPR units. Although the three tandem TPR repeats in both CTPR3 and UBP behave as a single unit, with no evidence of independent repeat motions, the data indicate that certain positions in UBP are significantly more flexible than are the corresponding positions in CTPR3. Most of the dynamical changes occur at or adjacent to positions that are involved in intra-repeat packing interactions. These observations lead us to suggest that the three-TPR domain of UBP does not incorporate optimized packing, compared to that seen in the idealized CTPR. The natural TPR domain is not only less stable overall than CTPR3, but also presents increased local flexibility at the positions where the sequences differs from the conserved consensus.     

Mapping the energy landscape of repeat proteins using NMR-detected hydrogen exchange

Repeat proteins contain tandem arrays of a simple structural motif. In contrast to globular proteins, repeat proteins are stabilized only by interactions between residues that are relatively close together in the sequence, with no ”long-range” interactions. Our work focuses on the tetratricopeptide repeat (TPR), a 34 amino acid helix-turn-helix motif found in tandem arrays in many natural proteins. Earlier, we reported the design and characterization of a series of consensus TPR (CTPR) proteins, which are built as arrays of multiple tandem copies of a 34 amino acid consensus sequence. Here, we present the results of extensive hydrogen exchange (HX) studies of the folding-unfolding behavior of two CTPR proteins (CTPR2 and CTPR3). We used HX to detect and characterize partially folded species that are populated at low frequency in the nominally folded state. We show that for both proteins the equilibrium folding-unfolding transition is non-two-state, but sequential, with the outermost helices showing a significantly higher probability than inner helices of being unfolded. We show that the experimentally observed unfolding behavior is consistent with the predictions of a simple Ising model, in which individual helices are treated as ”spin-equivalents”. The results that we present have general implications for our understanding of the thermodynamic properties of repeat proteins.     

Self‐association of TPR domains: Lessons learned from a designed,consensus‐based TPR oligomer

The tetratricopeptide repeat (TPR) motif is a protein–protein interaction module that acts as an organizing centre for complexes regulating a multitude of biological processes. Despite accumulating evidence for the formation of TPR oligomers as an additional level of regulation there is a lack of structural and solution data explaining TPR self‐association. In the present work we characterize the trimeric TPR‐containing protein YbgF, which is linked to the Tol system in Gram‐negative bacteria. By subtracting previously identified TPR consensus residues required for stability of the fold from residues conserved across YbgF homologs, we identified residues involved in oligomerization of the C‐terminal YbgF TPR domain. Crafting these residues, which are located in loop regions between TPR motifs, onto the monomeric consensus TPR protein CTPR3 induced the formation of oligomers. The crystal structure of this engineered oligomer shows an asymmetric trimer where stacking interactions between the introduced tyrosines and displacement of the C‐terminal hydrophilic capping helix, present in most TPR domains, are key to oligomerization. Asymmetric trimerization of the YbgF TPR domain and CTPR3Y3 leads to the formation of higher order oligomers both in the crystal and in solution. However, such open‐ended self‐association does not occur in full‐length YbgF suggesting that the protein's N‐terminal coiled‐coil domain restricts further oligomerization. This interpretation is borne out in experiments where the coiled‐coil domain of YbgF was engineered onto the N‐terminus of CTPR3Y3 and shown to block self‐association beyond trimerization. Our study lays the foundations for understanding the structural basis for TPR domain self‐association and how such self‐association can be regulated in TPR domain‐containing proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Relationships between protein structure and dynamics from a database of NMR-derived backbone order parameters

The amplitude of protein backbone NH group motions on a time-scale faster than molecular tumbling may be determined by analysis of (15)N NMR relaxation data according to the Lipari-Szabo model free formalism. An internet-accessible database has been compiled containing 1855 order parameters from 20 independent NMR relaxation studies on proteins whose three-dimensional structures are known. A series of statistical analyses has been performed to identify relationships between the structural features and backbone dynamics of these proteins. Comparison of average order parameters for different amino acid types indicates that amino acids with small side-chains tend to have greater backbone flexibility than those with large side-chains. In addition, the motions of a given NH group are also related to the sizes of the neighboring amino acids in the primary sequence. The secondary structural environment appears to influence backbone dynamics relatively weakly, with only subtle differences between the order parameter distributions of loop structures and regular hydrogen bonded secondary structure elements. However, NH groups near helix termini are more mobile on average than those in the central regions of helices. Tertiary structure influences are also relatively weak but in the expected direction, with more exposed residues being more flexible on average than residues that are relatively inaccessible to solvent.     

Determination of the molecular dynamics of alamethicin using 13C NMR: implications for the mechanism of gating of a voltage-dependent channel.

Alamethicin is a channel-forming peptide antibiotic that produces a highly voltage-dependent conductance in planar bilayers. To provide insight into the mechanisms for its voltage dependence, the dynamics of the peptide were examined in solution using nuclear magnetic resonance. Natural-abundance 13C spin-lattice relaxation rates and 13C-1H nuclear Overhauser effects of alamethicin were measured at two magnetic field strengths in methanol. This information was interpreted using a model-free approach to obtain values for the overall correlation times as well as the rates and amplitudes of the internal motions of the peptide. The picosecond, internal motions of alamethicin are highly restricted along the peptide backbone and indicate that it behaves as a rigid helical rod in solution. The side chain carbons exhibit increased segmental motion as their distance from the peptide backbone is increased; however, these motions are not unrestricted. Methyl group dynamics are also consistent with the restricted motions observed for the backbone carbons. There is no evidence from these dynamics measurements for a hinged motion of the peptide about proline-14. Alamethicin appears to be slightly less structured in methanol than in the membrane; as a result, alamethicin is also expected to behave as a rigid helix in the membrane. This suggests that the gating of this peptide involves changes in the orientation of the entire helix, rather than the movement of a segment of the peptide backbone.     

Calorimetric study of a series of designed repeat proteins: modular structure and modular folding

Repeat proteins comprise tandem arrays of a small structural motif. Their structure is defined and stabilized by interactions between residues that are close in the primary sequence. Several studies have investigated whether their structural modularity translates into modular thermodynamic properties. Tetratricopeptide repeat proteins (TPRs) are a class in which the repeated unit is a 34 amino acid helix-turn-helix motif. In this work, we use differential scanning calorimetry (DSC) to study the equilibrium stability of a series of TPR proteins with different numbers of an identical consensus repeat, from 2 to 20, CTPRa2 to CTPRa20. The DSC data provides direct evidence that the folding/unfolding transition of CTPR proteins does not fit a two-state folding model. Our results confirm and expand earlier studies on TPR proteins, which showed that apparent two-state unfolding curves are better fit by linear statistical mechanics models: 1D Ising models in which each repeat is treated as an independent folding unit.     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

Thermal stability of chicken brain α-spectrin repeat 17: a spectroscopic study

Spectrin is a rod-like multi-modular protein that is mainly composed of triple-helical repeats. These repeats show very similar 3D-structures but variable conformational and thermodynamical stabilities, which may be of great importance for the flexibility and dynamic behaviour of spectrin in the cell. For instance, repeat 17 (R17) of the chicken brain spectrin α-chain is four times less stable than neighbouring repeat 16 (R16) in terms of ?G. The structure of spectrin repeats has mainly been investigated by X-ray crystallography, but the structures of a few repeats, e.g. R16, have also been determined by NMR spectroscopy. Here, we undertook a detailed characterization of the neighbouring R17 by NMR spectroscopy. We assigned most backbone resonances and observed NOE restraints, relaxation values and coupling constants that all indicated that the fold of R17 is highly similar to that of R16, in agreement with previous X-ray analysis of a tandem repeat of the two domains. However, (15)N heteronuclear NMR spectra measured at different temperatures revealed particular features of the R17 domain that might contribute to its lower stability. Conformational exchange appeared to alter the linker connecting R17 to R16 as well as the BC-loop in close proximity. In addition, heat-induced splitting was observed for backbone resonances of a few spatially related residues including V99 of helix C, which in R16 is replaced by the larger hydrophobic tryptophan residue that is relatively conserved among other spectrin repeats. These data support the view that the substitution of tryptophan by valine at this position may contribute to the lower stability of R17.     

TPR重复和ELTR排型:重复的长度变化是一种功能进化的机制

TPR重复最初定义为一个34个氨基酸的蛋白结构重复。本文用PATTINPROT软件对基因库中TPR重复长度多样性进行了分析,发现40和42个氨基酸TPR重复大量存在。含40和42个氨基酸TPR重复序列蛋白的功能分析说明,这些长的TPR重复可以赋予蛋白新的功能。根据以上分析,提出重复序列的长度变化可能是功能进化的一种发生机制。     

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Collective motions on ns-µs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from eight different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second b strand (S2) and the second a helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.     

Removal of a consensus proline is not sufficient to allow tetratricopeptide repeat oligomerization

Tetratricopeptide repeat (TPR) domains are ubiquitous protein interaction domains that adopt a modular antiparallel array of α‐helices. The TPR fold typically adopts a monomeric state, and consensus TPRs sequences successfully fold into the expected monomeric topology. The versatility of the TPR fold also supports different quaternary structures, which may function as regulatory switches. One example is yeast mitochondrial fission 1 (Fis1) that appears to interconvert between monomer and dimer states in regulating division of peroxisomes and mitochondria. Whether human Fis1 can also interconvert like the yeast molecule is unknown. A TPR consensus proline residue present in human Fis1 is absent in the yeast molecule and, when added, prevents yeast Fis1 dimerization suggesting that the TPR consensus proline might have persisted to prevent TPR oligomerization. Here, we address this question with human Fis1 and the consensus TPR protein CTPR3. We demonstrate that human Fis1 does not form a noncovalent dimer via its TPR domain, despite conditions that favor dimerization of the yeast protein. We also show that the presence of the consensus proline is not sufficient to forbid TPR dimerization. Lastly, an analysis of all available TPR protein structures (22 nonredundant structures, totaling 64 TPRs—42 with the consensus proline and 22 without) revealed that the consensus proline is not necessary for turn formation, but does favor shorter turns. This work suggests the TPR consensus proline is not to prevent oligomerization, but to favor tight turns between repeats.     

Backbone motions in a crystalline protein from field-dependent 2H-NMR relaxation and line-shape analysis

We have used 2H-nmr to study backbone dynamics of the 2H-labeled, slowly exchanging amide sites of fully hydrated, crystalline hen egg white lysozyme. Order parameters are determined from the residual quadrupole coupling and values increase from S2 = 0.85 at 290 K to S2 = 0.94 at 200 K. Dynamical rates are determined from spin-lattice relaxation at three nmr frequencies (38.8, 61.5, and 76.7 MHz). The approach used here is thus distinct from solution nmr studies where dynamical amplitudes and rates are both determined from relaxation measurements. At temperatures below 250 K, relaxation is independent of the nmr frequency indicating that backbone motions are fast compared to the nmr frequencies. However, as the temperature is increased above 250 K, relaxation is significantly more efficient at the lowest frequency, which shows, in addition, the presence of motions that are slow compared to the nmr frequencies. Using the values of S2 determined from the residual quadrupole coupling and a model-free relaxation formalism that allows for fast and slow internal motions, we conclude that these slow motions have correlation times in the range of 0.1 to 1.0 microsecond and are effectively frozen out at 250 K where fast motions of the amide planes with approximately 15 ps effective correlation times and 9 degrees rms amplitudes dominate relaxation. The fast internal motions increase slightly in amplitude as the temperature rises toward 290 K, but the correlation time, as is also observed in solution nmr studies of RNase H, is approximately constant. These findings are consistent with hypotheses of dynamic glass transitions in hydrated proteins arising from temperature-dependent damping of harmonic modes of motion above the transition point.     

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.     

Left-Handed Intercalated DNA Double Helix: Rendezvous of Ethidium and Actinomycin D in the Z-Helical Conformation Space

Abstract

It is now very well recognized that the DNA double helix is conformationally pluralistic and that this flexibility is derived from internal motions due to backbone torsions. But what is less apparent is that such internal motions can occur in a correlated fashion and express themselves in a wide variety of structural motifs and phenomena. For example, flexibility inherent in the DNA molecule can lead to a family of Z-DNA, LZ1 and LZ2 being the two extremes and correlated internal motion can cause LZ1?LZ2 transition. More interestingly, such motions manifest themselves as breathing modes on the DNA lattice resulting in the sequence specific intercalation sites. Following a detailed stereochemical analyses we observed that the intercalation site for ethidium is located at the dCpdG sequence of the intercalated LZ1 helix (LZ1*) while that for actinomycin D is located at the dGpdC sequence of the intercalated LZ2 helix (LZ2*). From the stereochemistry of the drug binding we make experimentally testable predictions which are in fact supported by a few recent experimental studies. These studies also show that a left-handed intercalated B-DNA model is a viable intermediate in the Z to B transition which can hold the drug with binding energy comparable to that of the intercalated right-handed B-DNA.     

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Collective motions on ns-µs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from eight different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second b strand (S2) and the second a helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.Key words: prion proteins structural stability, molecular dynamics simulation, essential collective dynamics, protein dynamic domains, biomolecular NMR, rigid loop     

Dynamics of the Mrf-2 DNA-binding domain free and in complex with DNA

Mrf-2 is a member of a new class of DNA-binding proteins known as the AT-rich interaction domain family or ARID. Chemical shift indices and characteristic NOE values indicate that the three-dimensional structure of the Mrf-2 ARID in complex with DNA is nearly identical to that of the free protein. The backbone dynamics of the Mrf-2 domain free and in complex with DNA have been characterized by (15)N NMR relaxation measurements and model-free analysis. Chemical shift perturbations and dynamic studies suggest that two flexible interhelical loops, the flexible C-terminal tail, and one alpha-helix are involved in DNA recognition, indicating the importance of protein dynamics in DNA binding. Some well-structured regions, in particular the putative DNA-contacting helix, in Mrf-2 show a decrease in the order parameters (S(2)) upon complex formation. The less well-structured loops and the unstructured C-terminus show reduced flexibility upon DNA binding. In addition, the model-free analysis indicates motions on the picosecond to nanosecond and micro- to millisecond time scales at the DNA-binding surface of the bound Mrf-2 ARID, suggesting a model where interactions between the protein and DNA are highly dynamic.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Pathways and intermediates in forced unfolding of spectrin repeats

Spectrin repeats are triple-helical coiled-coil domains found in many proteins that are regularly subjected to mechanical stress. We used atomic force microscopy technique and steered molecular dynamics simulations to study the behavior of a wild-type spectrin repeat and two mutants. The experiments indicate that spectrin repeats can form stable unfolding intermediates when subjected to external forces. In the simulations the unfolding proceeded via a variety of pathways. Stable intermediates were associated to kinking of the central helix close to a proline residue. A mutant stabilizing the central helix showed no intermediates in experiments, in agreement with simulation. Spectrin repeats may thus function as elastic elements, extendable to intermediate states at various lengths.