Search

A collection of programs is presented to analyze the thermodynamics of folding of linear repeat proteins using a 1D Ising model to determine intrinsic folding and interfacial coupling free energies. Expressions for folding transitions are generated for a series of constructs with different repeat numbers and are globally fitted to transitions for these constructs. These programs are designed to analyze Ising parameters for capped homopolymeric consensus repeat constructs as well as heteropolymeric constructs that contain point substitutions, providing a rigorous framework for analysis of the effects of mutation on intrinsic and directional (i.e., N- vs. C-terminal) interfacial coupling free-energies. A bootstrap analysis is provided to estimate parameter uncertainty as well as correlations among fitted parameters. Rigorous statistical analysis is essential for interpreting fits using the complex models required for Ising analysis of repeat proteins, especially heteropolymeric repeat proteins. Programs described here are available at https://github.com/barricklab-at-jhu/Ising_programs.

Transcription activator-like effector proteins (TALEs) contain large numbers of repeats that bind double-stranded DNA, wrapping around DNA to form a continuous superhelix. Since unbound TALEs retain superhelical structure, it seems likely that DNA binding requires a significant structural distortion or partial unfolding. In this study, we use nearest-neighbor “Ising” analysis of consensus TALE (cTALE) repeat unfolding to quantify intrinsic folding free energies, coupling energies between repeats, and the free energy distribution of partly unfolded states, and to determine how those energies depend on the sequence that determines DNA-specificity (called the “RVD”). We find a moderate level of cooperativity for both the HD and NS RVD sequences (stabilizing interfaces combined with unstable repeats), as has been seen in other linear repeat proteins. Surprisingly, RVD sequence identity influences both the overall stability and the balance of intrinsic repeat stability and interfacial coupling energy. Using parameters from the Ising analysis, we have analyzed the distribution of partly folded states as a function of cTALE length and RVD sequence. We find partly unfolded states where one or more repeats are unfolded to be energetically accessible. Mixing repeats with different RVD sequences increases the population of partially folded states. Local folding free energies plateau for central repeats, suggesting that TALEs access partially folded states where a single internal repeat is unfolded while adjacent repeats remain folded. This breakage should allow TALEs to access superhelically-broken states, and may facilitate DNA binding.

Transcription activator-like effectors (TALEs) bind DNA through an array of tandem 34-residue repeats. How TALE repeat domains wrap around DNA, often extending more than 1.5 helical turns, without using external energy is not well understood. Here, we examine the kinetics of DNA binding of TALE arrays with varying numbers of identical repeats. Single molecule fluorescence analysis and deterministic modeling reveal conformational heterogeneity in both the free- and DNA-bound TALE arrays. Our findings, combined with previously identified partly folded states, indicate a TALE instability that is functionally important for DNA binding. For TALEs forming less than one superhelical turn around DNA, partly folded states inhibit DNA binding. In contrast, for TALEs forming more than one turn, partly folded states facilitate DNA binding, demonstrating a mode of ‘functional instability’ that facilitates macromolecular assembly. Increasing repeat number slows down interconversion between the various DNA-free and DNA-bound states., eLife digest The DNA contains all the information needed to build an organism. It is made up of two strands that wind around each other like a twisted ladder to form the double helix. The strands consist of sugar and phosphate molecules, which attach to one of for bases. Genes are built from DNA, and contain specific sequences of these bases. Being able to modify DNA by deleting, inserting or changing certain sequences allows researchers to engineer tissues or even organisms for therapeutical and practical applications. One of these gene editing tools is the so-called transcription activator-like effector protein (or TALE for short). TALE proteins are derived from bacteria and are built from simple repeating units that can be linked to form a string-like structure. They have been found to be unstable proteins. To bind to DNA, TALES need to follow the shape of the double helix, adopting a spiral structure, but how exactly TALE proteins thread their way around the DNA is not clear. To investigate this, Geiger-Schuller et al. monitored single TALE units using fluorescent microscopy. This way, they could exactly measure the time it takes for single TALE proteins to bind and release DNA. The results showed that some TALE proteins bind DNA quickly, whereas others do this slowly. Using a computer model to analyze the different speeds of binding suggested that the fast binding comes from partly unfolded proteins that quickly associate with DNA, and that the slow binding comes from rigid, folded TALE proteins, which have a harder time wrapping around DNA. This suggest that the unstable nature of TALEs, helps these proteins to bind to DNA and turn on genes. These findings will help to design future TALE-based gene editing tools and also provide more insight into how large molecules can assemble into complex structures. A next step will be to identify TALE repeats with unstable states and to test TALE gene editing tools that have intentionally placed unstable units.

Designed helical repeats (DHRs) are modular helix–loop–helix–loop protein structures that are tandemly repeated to form a superhelical array. Structures combining tandem DHRs demonstrate a wide range of molecular geometries, many of which are not observed in nature. Understanding cooperativity of DHR proteins provides insight into the molecular origins of Rosetta-based protein design hyperstability and facilitates comparison of energy distributions in artificial and naturally occurring protein folds. Here, we use a nearest-neighbor Ising model to quantify the intrinsic and interfacial free energies of four different DHRs. We measure the folding free energies of constructs with varying numbers of internal and terminal capping repeats for four different DHR folds, using guanidine-HCl and glycerol as destabilizing and solubilizing cosolvents. One-dimensional Ising analysis of these series reveals that, although interrepeat coupling energies are within the range seen for naturally occurring repeat proteins, the individual repeats of DHR proteins are intrinsically stable. This favorable intrinsic stability, which has not been observed for naturally occurring repeat proteins, adds to stabilizing interfaces, resulting in extraordinarily high stability. Stable repeats also impart a downhill shape to the energy landscape for DHR folding. These intrinsic stability differences suggest that part of the success of Rosetta-based design results from capturing favorable local interactions.

Rad6 is a yeast E2 ubiquitin conjugating enzyme that monoubiquitinates histone H2B in conjunction with the E3, Bre1, but can non-specifically modify histones on its own. We determined the crystal structure of a Rad6∼Ub thioester mimic, which revealed a network of interactions in the crystal in which the ubiquitin in one conjugate contacts Rad6 in another. The region of Rad6 contacted is located on the distal face of Rad6 opposite the active site, but differs from the canonical E2 backside that mediates free ubiquitin binding and polyubiquitination activity in other E2 enzymes. We find that free ubiquitin interacts weakly with both non-canonical and canonical backside residues of Rad6 and that mutations of non-canonical residues have deleterious effects on Rad6 activity comparable to those observed to mutations in the canonical E2 backside. The effect of non-canonical backside mutations is similar in the presence and absence of Bre1, indicating that contacts with non-canonical backside residues govern the intrinsic activity of Rad6. Our findings shed light on the determinants of intrinsic Rad6 activity and reveal new ways in which contacts with an E2 backside can regulate ubiquitin conjugating activity.