Backbone dynamics of sequence specific recognition and binding by the yeast Pho4 bHLH domain probed by NMR

Backbone dynamics of the basic/helix-loop-helix domain of Pho4 from Saccharomyces cerevisae have been probed by NMR techniques, in the absence of DNA, nonspecifically bound to DNA and bound to cognate DNA. Alpha proton chemical shift indices and nuclear Overhauser effect patterns were used to elucidate the secondary structure in these states. These secondary structures are compared to the co-crystal complex of Pho4 bound to a cognate DNA sequence (Shimizu T. Toumoto A, Ihara K, Shimizu M, Kyogou Y, Ogawa N, Oshima Y, Hakoshima T, 1997, EMBO J 15: 4689-4697). The dynamic information provides insight into the nature of this DNA binding domain as it progresses from free in solution to a specifically bound DNA complex. Relative to the unbound form, we show that formation of either the nonspecific and cognate DNA bound complexes involves a large change in conformation and backbone dynamics of the basic region. The nonspecific and cognate complexes, however, have nearly identical secondary structure and backbone dynamics. We also present evidence for conformational flexibility at a highly conserved glutamate basic region residue. These results are discussed in relation to the mechanism of sequence specific recognition and binding.     

Mechanical measurement of single-molecule binding rates: kinetics of DNA helix-destabilization by T4 gene 32 protein

Bacteriophage T4 gene 32 protein (gp32) is a single-stranded DNA (ssDNA) binding protein, and is essential for DNA replication, recombination and repair. While gp32 binds preferentially and cooperatively to ssDNA, it has not been observed to lower the thermal melting temperature of natural double-stranded DNA (dsDNA). However, in single-molecule stretching experiments, gp32 significantly destabilizes lambda DNA. In this study, we develop a theory of the effect of the protein on single dsDNA stretching curves, and apply it to the measured dependence of the DNA overstretching force on pulling rate in the presence of the full-length and two truncated forms of the protein. This allows us to calculate the rate of cooperative growth of single clusters of protein along ssDNA that are formed as the dsDNA molecule is stretched, as well as determine the site size of the protein binding to ssDNA. The rate of cooperative binding (ka) of both gp32 and of its proteolytic fragment *I (which lacks 48 residues from the C terminus) varies non-linearly with protein concentration, and appears to exceed the diffusion limit. We develop a model of protein association with the ends of growing clusters of cooperatively bound protein enhanced by 1-D diffusion along dsDNA, under the condition of protein excess. Upon globally fitting ka versus protein concentration, we determine the binding site size and the non-cooperative binding constants to dsDNA for gp32 and I. Our experiment mimics the growth of clusters of gp32 that likely exist at the DNA replication fork in vivo, and explains the origin of the "kinetic block" to dsDNA melting by gene 32 protein observed in thermal melting experiments.     

Analysis of the human replication protein A:Rad52 complex: evidence for crosstalk between RPA32, RPA70, Rad52 and DNA

The eukaryotic single-stranded DNA-binding protein, replication protein A (RPA), is essential for DNA replication, and plays important roles in DNA repair and DNA recombination. Rad52 and RPA, along with other members of the Rad52 epistasis group of genes, repair double-stranded DNA breaks (DSBs). Two repair pathways involve RPA and Rad52, homologous recombination and single-strand annealing. Two binding sites for Rad52 have been identified on RPA. They include the previously identified C-terminal domain (CTD) of RPA32 (residues 224-271) and the newly identified domain containing residues 169-326 of RPA70. A region on Rad52, which includes residues 218-303, binds RPA70 as well as RPA32. The N-terminal region of RPA32 does not appear to play a role in the formation of the RPA:Rad52 complex. It appears that the RPA32CTD can substitute for RPA70 in binding Rad52. Sequence homology between RPA32 and RPA70 was used to identify a putative Rad52-binding site on RPA70 that is located near DNA-binding domains A and B. Rad52 binding to RPA increases ssDNA affinity significantly. Mutations in DBD-D on RPA32 show that this domain is primarily responsible for the ssDNA binding enhancement. RPA binding to Rad52 inhibits the higher-order self-association of Rad52 rings. Implications for these results for the "hand-off" mechanism between protein-protein partners, including Rad51, in homologous recombination and single-strand annealing are discussed.     

Dynamic elements of replication protein A at the crossroads of DNA replication,recombination, and repair

Abstract

The heterotrimeric eukaryotic Replication protein A (RPA) is a master regulator of numerous DNA metabolic processes. For a long time, it has been viewed as an inert protector of ssDNA and a platform for assembly of various genome maintenance and signaling machines. Later, the modular organization of the RPA DNA binding domains suggested a possibility for dynamic interaction with ssDNA. This modular organization has inspired several models for the RPA-ssDNA interaction that aimed to explain how RPA, the high-affinity ssDNA binding protein, is replaced by the downstream players in DNA replication, recombination, and repair that bind ssDNA with much lower affinity. Recent studies, and in particular single-molecule observations of RPA-ssDNA interactions, led to the development of a new model for the ssDNA handoff from RPA to a specific downstream factor where not only stability and structural rearrangements but also RPA conformational dynamics guide the ssDNA handoff. Here we will review the current knowledge of the RPA structure, its dynamic interaction with ssDNA, and how RPA conformational dynamics may be influenced by posttranslational modification and proteins that interact with RPA, as well as how RPA dynamics may be harnessed in cellular decision making.     

Unraveling protein dynamics through fast spectral density mapping

Spectral density mapping at multiple NMR field strengths is probably the best method to describe the dynamical behavior of a protein in solution through the analysis of ¹⁵N heteronuclear relaxation parameters. Nevertheless, such analyses are scarcely reported in the literature, probably because this method is excessively demanding in spectrometer measuring time. Indeed, when using n different magnetic fields and assuming the validity of the high frequency approximation, the discrete sampling of the spectral density function with 2n + 1 points needs the measurement of 3n ¹⁵N heteronuclear relaxation measurements (n R ₁, n R ₂, and n¹⁵N{¹H}NOEs). Based on further approximations, we proposed a new strategy that allows us to describe the spectral density with n + 2 points, with the measurement of a total of n + 2 heteronuclear relaxation parameters. Applied to the dynamics analysis of the protein p13 ^MTCP1 at three different NMR fields, this approach allowed us to divide by nearly a factor of two the total measuring time, without altering further results obtained by the “model free” analysis of the resulting spectral densities. Furthermore, simulations have shown that this strategy remains applicable to any low isotropically tumbling protein ( ns), and is valid for the types of motion generally envisaged for proteins.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

Antibody variable region binding by Staphylococcal protein A: thermodynamic analysis and location of the Fv binding site on E-domain.

Immunoglobulins of human heavy chain subgroup III have a binding site for Staphylococcal protein A on the heavy chain variable domain (V(H)), in addition to the well-known binding site on the Fc portion of the antibody. Thermodynamic characterization of this binding event and localization of the Fv-binding site on a domain of protein A is described. Isothermal titration calorimetry (ITC) was used to characterize the interaction between protein A or fragments of protein A and variants of the hu4D5 antibody Fab fragment. Analysis of binding isotherms obtained for titration of hu4D5 Fab with intact protein A suggests that 3-4 of the five immunoglobulin binding domains of full length protein A can bind simultaneously to Fab with a Ka of 5.5+/-0.5 x 10(5) M(-1). A synthetic single immunoglobulin binding domain, Z-domain, does not bind appreciably to hu4D5 Fab, but both the E and D domains are functional for hu4D5 Fab binding. Thermodynamic parameters for titration of the E-domain with hu4D5 Fab are n = 1.0+/-0.1, Ka = 2.0+/-0.3 x 10(5) M(-1), and deltaH = -7.1+/-0.4 kcal mol(-1). Similar binding thermodynamics are obtained for titration of the isolated V(H) domain with E-domain indicating that the E-domain binding site on Fab resides within V(H). E-domain binding to an IgG1 Fc yields a higher affinity interaction with thermodynamic parameters n = 2.2+/-0.1, Ka > 1.0 x 10(7) M(-1), and deltaH = -24.6+/-0.6 kcal mol(-1). Fc does not compete with Fab for binding to E-domain indicating that the two antibody fragments bind to different sites. Amide 1H and 15N resonances that undergo large changes in NMR chemical shift upon Fv binding map to a surface defined by helix-2 and helix-3 of E-domain, distinct from the Fc-binding site observed in the crystal structure of the B-domain/Fc complex. The Fv-binding region contains negatively charged residues and a small hydrophobic patch which complements the basic surface of the region of the V(H) domain implicated previously in protein A binding.     

Phosphorylation-induced changes in backbone dynamics of the dematin headpiece C-terminal domain

Dematin is an actin-binding protein abundant in red blood cells and other tissues. It contains a villin-type ‘headpiece’ F-actin-binding domain at its extreme C-terminus. The isolated dematin headpiece domain (DHP) undergoes a significant conformational change upon phosphorylation. The mutation of Ser74 to Glu closely mimics the phosphorylation of DHP. We investigated motions in the backbone of DHP and its mutant DHPS74E using several complementary NMR relaxation techniques: laboratory frame ¹⁵N NMR relaxation, which is sensitive primarily to the ps–ns time scale, cross-correlated chemical shift modulation NMR relaxation detecting correlated μs–ms time scale motions of neighboring ¹³C′ and ¹⁵N nuclei, and cross-correlated relaxation of two ¹⁵N–¹H dipole–dipole interactions detecting slow motions of backbone NH vectors in successive amino acid residues. The results indicate a reduction in mobility upon the mutation in several regions of the protein. The additional salt bridge formed in DHPS74E that links the N- and C-terminal subdomains is likely to be responsible for these changes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Kinetic regulation of single DNA molecule denaturation by T4 gene 32 protein structural domains

Bacteriophage T4 gene 32 protein (gp32) specifically binds single-stranded DNA, a property essential for its role in DNA replication, recombination, and repair. Although on a thermodynamic basis, single-stranded DNA binding proteins should lower the thermal melting temperature of double-stranded DNA (dsDNA), gp32 does not. Using single molecule force spectroscopy, we show for the first time that gp32 is capable of slowly destabilizing natural dsDNA. Direct measurements of single DNA molecule denaturation and renaturation kinetics in the presence of gp32 and its proteolytic fragments reveal three types of kinetic behavior, attributable to specific protein structural domains, which regulate gp32's helix-destabilizing capabilities. Whereas the full-length protein exhibits very slow denaturation kinetics, a truncate lacking the acidic C-domain exhibits much faster kinetics. This may reflect a steric blockage of the DNA binding site and/or a conformational change associated with this domain. Additional removal of the N-domain, which is needed for binding cooperativity, further increases the DNA denaturation rate, suggesting that both of these domains are critical to the regulation of gp32's helix-destabilization capabilities. This regulation is potentially biologically significant because uncontrolled helix-destabilization would be lethal to the cell. We also obtain equilibrium measurements of the helix-coil transition free energy in the presence of these proteins for the first time.     

Determinants of replication protein A subunit interactions revealed using a phosphomimetic peptide

Replication protein A (RPA) is a eukaryotic ssDNA-binding protein and contains three subunits: RPA70, RPA32, and RPA14. Phosphorylation of the N-terminal region of the RPA32 subunit plays an essential role in DNA metabolism in processes such as replication and damage response. Phosphorylated RPA32 (pRPA32) binds to RPA70 and possibly regulates the transient RPA70-Bloom syndrome helicase (BLM) interaction to inhibit DNA resection. However, the structural details and determinants of the phosphorylated RPA32–RPA70 interaction are still unknown. In this study, we provide molecular details of the interaction between RPA70 and a mimic of phosphorylated RPA32 (pmRPA32) using fluorescence polarization and NMR analysis. We show that the N-terminal domain of RPA70 (RPA70N) specifically participates in pmRPA32 binding, whereas the unphosphorylated RPA32 does not bind to RPA70N. Our NMR data revealed that RPA70N binds pmRPA32 using a basic cleft region. We also show that at least 6 negatively charged residues of pmRPA32 are required for RPA70N binding. By introducing alanine mutations into hydrophobic positions of pmRPA32, we found potential points of contact between RPA70N and the N-terminal half of pmRPA32. We used this information to guide docking simulations that suggest the orientation of pmRPA32 in complex with RPA70N. Our study demonstrates detailed features of the domain-domain interaction between RPA70 and RPA32 upon phosphorylation. This result provides insight into how phosphorylation tunes transient bindings between RPA and its partners in DNA resection.     

Fast detection of quadruplex structure in DNA by the intrinsic fluorescence of a single-stranded DNA binding protein

Single-stranded guanine-rich (G-rich) DNA can fold into a four-stranded G-quadruplex structure and such structures are implicated in important biological processes and therapeutic applications. So far, bioinformatic analysis has identified up to several hundred thousand of putative quadruplex sequences in the genome of human and other animal. Given such a large number of sequences, a fast assay would be desired to experimentally verify the structure of these sequences. Here we describe a method that identifies the quadruplex structure by a single-stranded DNA binding protein from a thermoautotrophic archaeon. This protein binds single-stranded DNA in the unfolded, but not in the folded form. Upon binding to DNA, its fluorescence can be quenched by up to 70%. Formation of quadruplex greatly reduces fluorescence quenching in a K+-dependent manner. This structure-dependent quenching provides simple and fast detection of quadruplex in DNA at low concentration without DNA labelling.     

Different activities of the largest subunit of replication protein A cooperate during SV40 DNA replication

Replication protein A (RPA) is a stable heterotrimeric complex consisting of p70, p32 and p14 subunits. The protein plays a crucial role in SV40 minichromosome replication. Peptides of p70 representing interaction sites for the smaller two subunits, DNA as well as the viral initiator protein large T-antigen (Tag) and the cellular DNA polymerase alpha-primase (Pol) all interfered with the replication process indicating the importance of the different p70 activities in this process. Inhibition by the peptide disrupting protein-protein interactions was observed only during the pre-initiation stage prior to primer synthesis, suggesting the formation of a stable initiation complex between RPA, Tag and Pol at the primer end.     

Identification of a residue critical for maintaining the functional conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.     

A Single Residue Determines the Cooperative Binding Property of a Primosomal DNA Replication Protein,PriB, to Single-Stranded DNA

PriB is a primosomal protein required for re-initiation of replication in bacteria. We characterized and compared the DNA-binding properties of PriB from Salmonella enterica serovar Typhimurium LT2 (StPriB) and Escherichia coli (EcPriB). Only one residue of EcPriB, V6, was different in StPriB (replaced by A6). Previous structural information revealed that this residue is located on the putative dimer-dimer interface of PriB and is not involved in single-stranded DNA (ssDNA) binding. The cooperative binding mechanism of StPriB to DNA is, however, very different from that of EcPriB. Unlike EcPriB, which forms a single complex with ssDNAs of various lengths, StPriB forms two or more distinct complexes. Based on these results, as well as information on structure, binding modes for forming a stable complex of PriB with ssDNA of 25 nucleotides (nt), (EcPriB)₂₅, and (StPriB)₂₅ are proposed.     

Backbone dynamics of the human MIA protein studied by (15)N NMR relaxation: implications for extended interactions of SH3 domains

The melanoma inhibitory activity (MIA) protein is a clinically valuable marker in patients with malignant melanoma as enhanced values diagnose metastatic melanoma stages III and IV. Here, we report the backbone dynamics of human MIA studied by (15)N NMR relaxation experiments. The folded core of human MIA is found to be rigid, but several loops connecting beta-sheets, such as the RT-loop for example, display increased mobility on picosecond to nanosecond time scales. One of the most important dynamic features is the pronounced flexibility of the distal loop, comprising residues Asp 68 to Ala 75, where motions on time scales up to milliseconds occur. Further, significant exchange contributions are observed for residues of the canonical binding site of SH3 domains including the RT-loop, the n-Src loop, for the loop comprising residues 13 to 19, which we refer to as the"disulfide loop", in part for the distal loop, and the carboxyl terminus of human MIA. The functional importance of this dynamic behavior is discussed with respect to the biological activity of several point mutations of human MIA. The results of this study suggest that the MIA protein and the recently identified highly homologous fibrocyte-derived protein (FDP)/MIA-like (MIAL) constitute a new family of secreted proteins that adopt an SH3 domain-like fold in solution with expanded ligand interactions.     

Structure and dynamics of the N-terminal half of hepatitis C virus core protein: An intrinsically unstructured protein

Hepatitis C virus core protein plays an important role in the assembly and packaging of the viral genome. We have studied the structure of the N-terminal half of the core protein (C82) which was shown to be sufficient for the formation of nucleocapsid-like particle (NLP) in vitro and in yeast. Structural bioinformatics analysis of C82 suggests that it is mostly unstructured. Circular dichroism and structural NMR data indicate that C82 lacks secondary structure. Moreover, NMR relaxation data shows that C82 is highly disordered. These results indicate that the N-terminal half of the HCV core protein belongs to the growing family of intrinsically unstructured proteins (IUP). This explains the tendency of the hepatitis C virus core protein to interact with several host proteins, a well-documented characteristic of IUPs.     

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its “essential space,” are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the β-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion. Proteins 31:370–382, 1998. © 1998 Wiley-Liss, Inc.     

Comparison of (13)C(alpha)H and (15)NH backbone dynamics in protein GB1

This study presents a site-resolved experimental view of backbone C(alpha)H and NH internal motions in the 56-residue immunoglobulin-binding domain of streptococcal protein G, GB1. Using (13)C(alpha)H and (15)NH NMR relaxation data [T(1), T(2), and NOE] acquired at three resonance frequencies ((1)H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(omega) = 2omegaJ(omega) to provide a model-independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(omega) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all C(alpha)H and NH groups. Deconvolution of F(omega) curves effectively separates overall tumbling and internal motional correlation time distributions to yield more accurate order parameters than determined by using standard model free approaches. Compared to NH groups, C(alpha)H internal motions are more broadly distributed on the nanosecond time scale, and larger C(alpha)H order parameters are related to correlated bond rotations for C(alpha)H fluctuations. Motional parameters for NH groups are more structurally correlated, with NH order parameters, for example, being larger for residues in more structured regions of beta-sheet and helix and generally smaller for residues in the loop and turns. This is most likely related to the observation that NH order parameters are correlated to hydrogen bonding. This study contributes to the general understanding of protein dynamics and exemplifies an alternative and easier way to analyze NMR relaxation data.