A protein-DNA binding mechanism proceeds through multi-state or two-state parallel pathways

Kinetics of DNA strand transfer between polymerase and proofreading exonuclease active sites regulates error correction during high-fidelity replication

We show that T7 DNA polymerase (pol) and exonuclease (exo) domains contribute to selective error correction during DNA replication by regulating bidirectional strand transfer between the two active sites. To explore the kinetic basis for selective removal of mismatches, we used a fluorescent cytosine analog (1,3-diaza-2-oxophenoxazine) to monitor the kinetics of DNA transfer between the exo and pol sites. We globally fit stopped-flow fluorescence and base excision kinetic data and compared results obtained with ssDNA versus duplex DNA to resolve how DNA transfer governs exo specificity. We performed parallel studies using hydrolysis-resistant phosphorothioate oligonucleotides to monitor DNA transfer to the exo site without hydrolysis. ssDNA binds to the exo site at the diffusion limit (10⁹ M^-1 s^-1, K_d = 40 nM) followed by fast hydrolysis of the 3'-terminal nucleotide (>5000 s^-1). Analysis using duplex DNA with a 3'-terminal mismatch or a buried mismatch exposed a unique intermediate state between pol and exo active sites and revealed that transfer via the intermediate to the exo site is stimulated by free nucleoside triphosphates. Transfer from the exo site back to the pol site after cleavage is fast and efficient. We propose a model to explain why buried mismatches are removed faster than single 3'-terminal mismatches and thereby provide an additional opportunity for error correction. Our data provide the first comprehensive model to explain how DNA transfer from pol to exo active sites and back again after base excision allow efficient selective mismatch removal during DNA replication to improve fidelity by more than 1000-fold.

Biophysical Properties of Self-Assembled Immune Signals Impact Signal Processing and the Nature of Regulatory Immune Function

Outcomes during immunotherapy are impacted not only by the specific therapeutic signals and pharmacodynamics, but also by the biophysical forms in which signals are delivered. This integration is determinative in autoimmunity because the disease is caused by immune dysregulation and inflammation. Unfortunately, the links between nanomaterial design, biophysical properties, and immune regulation are poorly defined. Here we designed cationic peptide antigens with defined charge distributions and then used electrostatics to assemble these peptides into complexes with anionic regulatory cues. We first show complexes induce antigen-specific tolerance during myelin-driven autoimmunity. We next show the affinity between these immune cues is controlled by charge balance and that affinity confers distinct biophysical properties important in immunological processing, including antigen availability. The underlying binding affinities between the self-assembled signals influences inflammatory gene expression in dendritic cells and antigen-specific regulatory outcomes in self-reactive transgenic T cells. This granular understanding of nanomaterial-immune interactions contributes to a more rational immunotherapy design.

Kinetics of DNA-protein association and dissociation by stopped-flow spectroscopy

The ability of proteins to recognize a specific nucleotide sequence or structural feature in the genome is central to DNA replication and repair as well as gene transcription. To understand how DNA-binding proteins accomplish this feat it is important to characterize how they efficiently recognize specific DNA sequences or various structural features in the milieu of a crowded nucleus containing 50mM basepairs of DNA. This minimally entails the determination of thermodynamic (equilibrium) parameters that characterize the relative affinities of DNA-binding proteins with their cognate and non-cognate DNA sequences of interest. Additionally, elucidating the kinetic parameters that govern DNA-protein association and dissociation can provide essential mechanistic insights into the underlying mechanisms, and can inform on how these processes can be controlled under varying cellular conditions. Unfortunately, the acquisition of kinetic data can appear intimidating and are usually not included in the characterization of most DNA binding proteins. Here we provide detailed theoretical considerations and practical, easy-to-follow protocols for using stopped-flow spectroscopy to study the kinetics of DNA-protein association and dissociation, using poly-(ADP ribose) polymerase 1 as a practical example.

Poly(ADP-ribose) polymerase 1 searches DNA via a 'monkey bar' mechanism

Poly(ADP-ribose) polymerase 1 (PARP1) is both a first responder to DNA damage and a chromatin architectural protein. How PARP1 rapidly finds DNA damage sites in the context of a nucleus filled with undamaged DNA, to which it also binds, is an unresolved question. Here, we show that PARP1 association with DNA is diffusion-limited, and release of PARP1 from DNA is promoted by binding of an additional DNA molecule that facilitates a ‘monkey bar’ mechanism, also known as intersegment transfer. The WGR-domain of PARP1 is essential to this mechanism, and a point mutation (W589A) recapitulates the altered kinetics of the domain deletion. Demonstrating the physiological importance of the monkey bar mechanism for PARP1 function, the W589A mutant accumulates at sites of DNA damage more slowly following laser micro-irradiation than wild-type PARP1. Clinically relevant inhibitors of PARP1 did not alter the rate or mechanism of the release of PARP1 from DNA.

Our cells constantly withstand damage that can lead to breaks in the strands of our DNA. These cuts need to be fixed for the cell to stay healthy. When a break happens, one of the first responders to the scene is a protein known as PARP1. It binds to the ruptured strand (or strands) and then it recruits other repair agents to that location. But first, PARP1 needs to scan for cuts and notches amongst an overwhelming amount of DNA that is still intact. This is a complicated task, especially since the protein tends to bind both broken and unbroken DNA. How does it not stay ‘stuck’ on an undamaged portion of the genome?

Here, Rudolph et al. use a combination of biochemical techniques and cell biology to show that PARP1 travels through our genome by swinging from one DNA location to another, the way a child grabs onto monkey bars. One of the DNA-binding domains of PARP1, known as the WGR-domain, acts like an arm and initiates the movement by gripping onto a new segment of DNA. In fact, chopping off the WGR-domain or disabling it through mutations makes PARP1 worse at finding DNA breakages in the cell.

Unfixed DNA damage can lead to a cell becoming cancerous; ultimately, if the breakages keep accumulating, the cell does not survive. This makes PARP1 an important target for cancer treatment. Indeed, certain drugs already rely on trapping the protein so that tumor cells die. Understanding how cells cope with DNA damage and exactly how PARP1 works could help in the fight against cancer.

Direct observation of a transient ternary complex during IκBα-mediated dissociation of NF-κB from DNA

We previously demonstrated that IκBα markedly increases the dissociation rate of DNA from NF-κB. The mechanism of this process remained a puzzle because no ternary complex was observed, and structures show that the DNA and IκBα binding sites on NF-κB are overlapping. The kinetics of interaction of IκBα with NF-κB and its complex with DNA were analyzed by using stopped-flow experiments in which fluorescence changes in pyrene-labeled DNA or the native tryptophan in IκBα were monitored. Rate constants governing the individual steps in the reaction were obtained from analysis of the measured rate vs. concentration profiles. The NF-κB association with DNA is extremely rapid with a rate constant of 1.5 × 10(8) M(-1)⋅s(-1). The NF-κB-DNA complex dissociates with a rate constant of 0.41 s(-1), yielding a KD of 2.8 nM. When IκBα is added to the NF-κB-DNA complex, we observe the formation of a transient ternary complex in the first few milliseconds of the fluorescence trace, which rapidly rearranges to release DNA. The rate constant of this IκBα-mediated dissociation is nearly equal to the rate constant of association of IκBα with the NF-κB-DNA complex, showing that IκBα is optimized to repress transcription. The rate constants for the individual steps of a more folded mutant IκBα were also measured. This mutant associates with NF-κB more rapidly than wild-type IκBα, but it associates with the NF-κB-DNA complex more slowly and also is less efficient at mediating dissociation of the NF-κB-DNA complex.

Minute time scale prolyl isomerization governs antibody recognition of an intrinsically disordered immunodominant epitope

Conformational rearrangements in antibody·antigen recognition are essential events where kinetic discrimination of isomers expands the universe of combinations. We investigated the interaction mechanism of a monoclonal antibody, M1, raised against E7 from human papillomavirus, a prototypic viral oncoprotein and a model intrinsically disordered protein. The mapped 12-amino acid immunodominant epitope lies within a "hinge" region between the N-terminal intrinsically disordered and the C-terminal globular domains. Kinetic experiments show that despite being within an intrinsically disordered region, the hinge E7 epitope has at least two populations separated by a high energy barrier. Nuclear magnetic resonance traced the origin of this barrier to a very slow (t(1/2)∼4 min) trans-cis prolyl isomerization event involving changes in secondary structure. The less populated (10%) cis isomer is the binding-competent species, thus requiring the 90% of molecules in the trans configuration to isomerize before binding. The association rate for the cis isomer approaches 6 × 10(7) M(-1) s(-1), a ceiling for antigen-antibody interactions. Mutagenesis experiments showed that Pro-41 in E7Ep was required for both binding and isomerization. After a slow postbinding unimolecular rearrangement, a consolidated complex with K(D) = 1.2 × 10(-7) M is reached. Our results suggest that presentation of this viral epitope by the antigen-presenting cells would have to be "locked" in the cis conformation, in opposition to the most populated trans isomer, in order to select the specific antibody clone that goes through affinity and kinetic maturation.

Classes of fast and specific search mechanisms for proteins on DNA

Problems of search and recognition appear over different scales in biological systems. In this review we focus on the challenges posed by interactions between proteins, in particular transcription factors, and DNA and possible mechanisms which allow for fast and selective target location. Initially we argue that DNA-binding proteins can be classified, broadly, into three distinct classes which we illustrate using experimental data. Each class calls for a different search process and we discuss the possible application of different search mechanisms proposed over the years to each class. The main thrust of this review is a new mechanism which is based on barrier discrimination. We introduce the model and analyze in detail its consequences. It is shown that this mechanism applies to all classes of transcription factors and can lead to a fast and specific search. Moreover, it is shown that the mechanism has interesting transient features which allow for stability at the target despite rapid binding and unbinding of the transcription factor from the target.

Protein flexibility directs DNA recognition by the papillomavirus E2 proteins

Although DNA flexibility is known to play an important role in DNA–protein interactions, the importance of protein flexibility is less well understood. Here, we show that protein dynamics are important in DNA recognition using the well-characterized human papillomavirus (HPV) type 6 E2 protein as a model system. We have compared the DNA binding properties of the HPV 6 E2 DNA binding domain (DBD) and a mutant lacking two C-terminal leucine residues that form part of the hydrophobic core of the protein. Deletion of these residues results in increased specific and non-specific DNA binding and an overall decrease in DNA binding specificity. Using ¹⁵N NMR relaxation and hydrogen/deuterium exchange, we demonstrate that the mutation results in increased flexibility within the hydrophobic core and loop regions that orient the DNA binding helices. Stopped-flow kinetic studies indicate that increased flexibility alters DNA binding by increasing initial interactions with DNA but has little or no effect on the structural rearrangements that follow this step. Taken together these data demonstrate that subtle changes in protein dynamics have a major influence on protein–DNA interactions.

Thermodynamics of cooperative DNA recognition at a replication origin and transcription regulatory site

Binding cooperativity guides the formation of protein−nucleic acid complexes, in particular those that are highly regulated such as replication origins and transcription sites. Using the DNA binding domain of the origin binding and transcriptional regulator protein E2 from human papillomavirus type 16 as model, and through isothermal titration calorimetry analysis, we determined a positive, entropy-driven cooperativity upon binding of the protein to its cognate tandem double E2 site. This cooperativity is associated with a change in DNA structure, where the overall B conformation is maintained. Two homologous E2 domains, those of HPV18 and HPV11, showed that the enthalpic−entropic components of the reaction and DNA deformation can diverge. Because the DNA binding helix is almost identical in the three domains, the differences must lie dispersed throughout this unique dimeric β-barrel fold. This is in surprising agreement with previous results for this domain, which revealed a strong coupling between global dynamics and DNA recognition.

Mutational analysis of kinetic partitioning in protein folding and protein-DNA binding

Kinetic partitioning between competing routes is present in many biological processes. Here, we propose a methodology to characterize kinetic partitioning through site-directed mutagenesis and apply it to parallel routes for unfolding of the TI I27 protein and for recognition of its target DNA by the human papillomavirus E2 protein. The balance between the two competing reaction routes can be quantified by the partitioning constant K(p). K(p) is easily modulated by point mutations, opening the way for the rational design of kinetic partitioning. Conserved wild-type residues strongly favor one of the two competing reactions, suggesting that in these systems there is an evolutionary pressure to shift partitioning towards a certain route. The mutations with the largest effects on partitioning cluster together in space, defining the protein regions most relevant for the modulation of partitioning. Such regions are neither fully coincident with nor strictly segregated from the regions that are important from each competing reaction. We dissected the mutational effects on partitioning into the contributions from each competing route using a new parameter called pi-value. The results suggest how the design of kinetic partitioning may be approached in each case.

Experimental snapshots of a protein-DNA binding landscape

Protein recognition of DNA sites is a primary event for gene function. Its ultimate mechanistic understanding requires an integrated structural, dynamic, kinetic, and thermodynamic dissection that is currently limited considering the hundreds of structures of protein-DNA complexes available. We describe a protein-DNA-binding pathway in which an initial, diffuse, transition state ensemble with some nonnative contacts is followed by formation of extensive nonnative interactions that drive the system into a kinetic trap. Finally, nonnative contacts are slowly rearranged into native-like interactions with the DNA backbone. Dissimilar protein-DNA interfaces that populate along the DNA-binding route are explained by a temporary degeneracy of protein-DNA interactions, centered on "dual-role" residues. The nonnative species slow down the reaction allowing for extended functionality.

Searching fast for a target on DNA without falling to traps

Genomic expression depends critically on both the ability of regulatory proteins to locate specific target sites on DNA within seconds and on the formation of long-lived (many minutes) complexes between these proteins and the DNA. Equilibrium experiments show that indeed regulatory proteins bind tightly to their target site. However, they also find strong binding to other nonspecific sites which act as traps that can dramatically increase the time needed to locate the target. This gives rise to a conflict between the speed and stability requirements. Here we suggest a simple mechanism which can resolve this long-standing paradox.

A strained DNA binding helix is conserved for site recognition, folding nucleation, and conformational modulation

Nucleic acid recognition is often mediated by alpha-helices or disordered regions that fold into alpha-helix on binding. A peptide bearing the DNA recognition helix of HPV16 E2 displays type II polyproline (PII) structure as judged by pH, temperature, and solvent effects on the CD spectra. NMR experiments indicate that the canonical alpha-helix is stabilized at the N-terminus, while the PII forms at the C-terminus half of the peptide. Re-examination of the dihedral angles of the DNA binding helix in the crystal structure and analysis of the NMR chemical shift indexes confirm that the N-terminus half is a canonical alpha-helix, while the C-terminal half adopts a 3(10) helix structure. These regions precisely match two locally driven folding nucleii, which partake in the native hydrophobic core and modulate a conformational switch in the DNA binding helix. The peptide shows only weak and unspecific residual DNA binding, 10(4)-fold lower affinity, and 500-fold lower discrimination capacity compared with the domain. Thus, the precise side chain conformation required for modulated and tight physiological binding by HPV E2 is largely determined by the noncanonical strained alpha-helix conformation, "presented" by this unique architecture. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 432-443, 2009.

The interaction of hyperthermophilic TATA-box binding protein with single-stranded DNA is entropically favorable and exhibits a large negative heat capacity change at high salt concentration

We have investigated the thermodynamics of the interaction between the TATA-box-binding protein from Pyrococcus horikoshii (PhoTBP) and its target DNA (TATA-1). The interaction between PhoTBP and double-stranded DNA (dsDNA) is entropically favorable and enthalpically unfavorable. The thermodynamic parameters for TATA-1 duplex formation in the presence of PhoTBP, that is, ternary PhoTBP-dsDNA complexation, are similar to those for TATA-1 duplex formation, which is enthalpically favorable. Surface plasmon resonance analysis indicates that the interaction between PhoTBP and single-stranded DNA (ssDNA) of TATA-1 is entropy driven and has a large negative heat capacity change (-1.19 kcal mol(-1) K(-1)) at high salt concentration (800 mM NaCl). These results suggest that the favorable entropic effect corresponding to the interaction between PhoTBP and dsDNA is due not to ternary complexation but to the interaction between PhoTBP and ssDNA. This report is the first to describe the thermodynamics of the interaction between TBP and ssDNA.

Transition state for protein-DNA recognition

We describe the formation of protein-DNA contacts in the two-state route for DNA sequence recognition by a transcriptional regulator. Surprisingly, direct sequence readout establishes in the transition state and constitutes the bottleneck of complex formation. Although a few nonspecific ionic interactions are formed at this early stage, they mainly play a stabilizing role in the final consolidated complex. The interface is fairly plastic in the transition state, likely because of a high level of hydration. The overall picture of this two-state route largely agrees with a smooth energy landscape for binding that speeds up DNA recognition. This "direct" two-state route differs from the parallel multistep pathway described for this system, which involves nonspecific contacts and at least two intermediate species that must involve substantial conformational rearrangement in either or both macromolecules.

Kinetically driven refolding of the hyperstable EBNA1 origin DNA-binding dimeric beta-barrel domain into amyloid-like spherical oligomers

The Epstein-Barr nuclear antigen 1 (EBNA1) is essential for DNA replication and episome segregation of the viral genome, and participates in other gene regulatory processes of the Epstein-Barr virus in benign and malignant diseases related to this virus. Despite the participation of other regions of the protein in evading immune response, its DNA binding, dimeric beta-barrel domain (residues 452-641) is necessary and sufficient for the main functions. This domain has an unusual topology only shared by another viral origin binding protein (OBP), the E2 DNA binding domain of papillomaviruses. Both the amino acid and DNA target sequences are completely different for these two proteins, indicating a link between fold conservation and function. In this work we investigated the folding and stability of the DNA binding domain of EBNA1 OBP and found it is extremely resistant to chemical, temperature, and pH denaturation. The thiocyanate salt of guanidine is required for obtaining a complete transition to a monomeric unfolded state. The unfolding reaction is extremely slow and shows a marked uncoupling between tertiary and secondary structure, indicating the presence of intermediate species. The Gdm.SCN unfolded protein refolds to fully soluble and spherical oligomeric species of 1.2 MDa molecular weight, with identical fluorescence centre of spectral mass but different intensity and different secondary structure. The refolded spherical oligomers are substantially less stable than the native recombinant dimer. In keeping with the substantial structural rearrangement in the oligomers, the spherical oligomers do not bind DNA, indicating that the DNA binding site is either disrupted or participates in the oligomerization interface. The puzzling extreme stability of a dimeric DNA binding domain from a protein from a human infecting virus in addition to a remarkable kinetically driven folding where all molecules do not return to the most stable original species suggests a co-translational and directional folding of EBNA1 in vivo, possibly assisted by folding accessory proteins. Finally, the oligomers bind Congo red and thioflavin-T, both characteristic of repetitive beta-sheet elements of structure found in amyloids and their soluble precursors. The stable nature of the "kinetically trapped" oligomers suggest their value as models for understanding amyloid intermediates, their toxic nature, and the progress to amyloid fibers in misfolding diseases. The possible role of the EBNA1 spherical oligomers in the virus biology is discussed.

Antibody recognition of a flexible epitope at the DNA binding site of the human papillomavirus transcriptional regulator E2

We have obtained a monoclonal antibody (ED15) against the C-terminal DNA-binding domain of the high-risk human papillomavirus strain-16 E2 protein that strongly interferes with its DNA-binding activity. We here characterize the recognition mechanism of this antibody and find that the ED15-E2 interaction has a strong electrostatic component, which correlates with the high proportion of acidic residues found in the antibody combining site. Further circular dichroism experiments in the presence of phosphate show that, in addition to electrostatic screening of key potential interactions, ionic strength affects the conformation of the epitope. In addition, the interaction is strongly modulated by pH, which correlates with the local flexibility of the epitope rather than the presence of pH sensitive residues at the interface. Noticeably, this finding is well correlated with the strong entropic component of the interaction. Site directed mutagenesis indicates that the ED15 epitope involves at least part of the DNA-binding helix of E2, explaining the mAb inhibitory activity. At physiological salt concentrations, the equilibrium dissociation constant of the E2-ED15 interaction is 10(-7) M and the association rate is 10(4) M-1 s-1, at least 1 order of magnitude slower than those generally reported in the most extensively described "nonflexible" antibody-protein interactions, indicating the presence of a slow conformational rearrangement on the antigen as the rate-limiting step. The crucial role of antigen flexibility in antibody-protein recognition is discussed.

Comprehensive comparison of the interaction of the E2 master regulator with its cognate target DNA sites in 73 human papillomavirus types by sequence statistics

Mucosal human papillomaviruses (HPVs) are etiological agents of oral, anal and genital cancer. Properties of high- and low-risk HPV types cannot be reduced to discrete molecular traits. The E2 protein regulates viral replication and transcription through a finely tuned interaction with four sites at the upstream regulatory region of the genome. A computational study of the E2–DNA interaction in all 73 types within the alpha papillomavirus genus, including all known mucosal types, indicates that E2 proteins have similar DNA discrimination properties. Differences in E2–DNA interaction among HPV types lie mostly in the target DNA sequence, as opposed to the amino acid sequence of the conserved DNA-binding alpha helix of E2. Sequence logos of natural and in vitro selected sites show an asymmetric pattern of conservation arising from indirect readout, and reveal evolutionary pressure for a putative methylation site. Based on DNA sequences only, we could predict differences in binding energies with a standard deviation of 0.64 kcal/mol. These energies cluster into six discrete affinity hierarchies and uncovered a fifth E2-binding site in the genome of six HPV types. Finally, certain distances between sites, affinity hierarchies and their eventual changes upon methylation, are statistically associated with high-risk types.

Increased stability and DNA site discrimination of "single chain" variants of the dimeric beta-barrel DNA binding domain of the human papillomavirus E2 transcriptional regulator

Human papillomavirus infects millions of people worldwide and is a causal agent of cervical cancer in women. The HPV E2 protein controls the expression of all viral genes through binding of its dimeric C-terminal domain (E2C) to its target DNA site. We engineered monomeric versions of the HPV16 E2C, in order to probe the link of the dimeric beta-barrel fold to stability, dimerization, and DNA binding. Two single-chain variants, with 6 and 12 residue linkers (scE2C-6 and scE2C-12), were purified and characterized. Spectroscopy and crystallography show that the native structure is unperturbed in scE2C-12. The single chain variants are stabilized with respect to E2C, with effective concentrations of 0.6 to 6 mM. The early folding events of the E2C dimer and scE2C-12 are very similar and include formation of a compact species in the submillisecond time scale and a non-native monomeric intermediate with a half-life of 25 ms. However, monomerization changes the unfolding mechanism of the linked species from two-state to three-state, with a high-energy intermediate. Binding to the specific target site is up to 5-fold tighter in the single chain variants. Nonspecific DNA binding is up to 7-fold weaker in the single chain variants, leading to an overall 10-fold increased site discrimination capacity, the largest described so far for linked DNA binding domains. Titration calorimetric binding analysis, however, shows almost identical behavior for dimer and single-chain species, suggesting very subtle changes behind the increased specificity. Global analysis of the mechanisms probed suggests that the dynamics of the E2C domain, rather than the structure, are responsible for the differential properties. Thus, the plastic and dimeric nature of the domain did not evolve for a maximum affinity, specificity, and stability of the quaternary structure, likely because of regulatory reasons and for roles other than DNA binding played by partly folded dimeric or monomeric conformers.

An analytical rate expression for the kinetics of gene transcription mediated by dimeric transcription factors

To model gene transcription kinetics, empirical fitting with the Hill function or S-system is often used. In this study, we derived an analytical expression for gene transcription rates in a manner similar to that developed for enzyme kinetics to describe the kinetics of gene transcription mediated by dimeric transcription factors (TFs) such as Gcn4p, a Saccharomyces cerevisiae master gene regulator. We showed that the analytical rate expression and its parameters estimated from several sets of experimental data could accurately reproduce the experimentally measured promoter-binding activity of Gcn4p. Furthermore, the analytical rate expression allowed us to derive analytically, rather than fit empirically, the parameters of the Hill function and S-system for use in modelling transcription kinetics. We found that a plot of gene transcription rate against Gcn4p concentration gave a sigmoidal dose-response curve with a positive co-operativity Hill coefficient (approximately 1.25), in accordance with previous experimental findings on the promoter binding of dimeric TFs. The characteristics of the dose-response curve around the estimated cellular Gcn4p concentration suggest that transcription regulation is efficiently controlled under physiological conditions. This work is a useful initial step towards analytically modelling and simulating complicated gene transcription networks.

A quasi-spontaneous amyloid route in a DNA binding gene regulatory domain: The papillomavirus HPV16 E2 protein

The DNA binding domain of papillomavirus E2 proteins is at the center of the regulation of gene transcription and replication of the virus. Its unique fold consists of a beta-barrel domain that combines an eight-stranded dimeric beta-barrel core interface with two symmetrical DNA binding alpha-helices and other two helices, packed against the central barrel. Treatment with low amounts of trifluoroethanol readily leads to a mostly beta-sheet oligomeric species, with a loss of near-UV circular dichroism signal and increase in its ANS binding capacity, indicating that buried hydrophobic surfaces become accessible to the solvent. This species subsequently undergoes a slow transition into amyloid aggregates as determined by light scattering and Congo red and thioflavin T binding. Electron microscopy shows short amyloid fibers with a curly aspect as the end product. The amyloid route is completely prevented by addition of stoichiometrical amounts of specific DNA, strongly suggesting that unfolding of the DNA binding alpha-helix is required for the formation of the intermediate. The slow nature of this expanded beta-oligomeric species and the availability of several different conformational probes make it an excellent model for investigating amyloid mechanisms. The mild perturbation required for entering an amyloid route is indicative of a preexisting equilibrium. Oligomerization processes are required for the assembly of transcription initiation and DNA replication machineries, where proteins from different viruses must come together with host cell proteins. The E2 protein is a virus-encoded multifunctional master regulator that may exert one of its multiple functions through its ability to oligomerize.

Stepwise binding and bending of DNA by Escherichia coli integration host factor

Integration host factor (IHF) is a prokaryotic protein required for the integration of lambda phage DNA into its host genome. An x-ray crystal structure of the complex shows that IHF binds to the minor groove of DNA and bends the double helix by 160 degrees [Rice PA, Yang S, Mizuuchi K, Nash HA (1996) Cell 87:1295-1306]. We sought to dissect the complex formation process into its component binding and bending reaction steps, using stopped-flow fluorimetry to observe changes in resonance energy transfer between DNA-bound dyes, which in turn reflect distance changes upon bending. Different DNA substrates that are likely to increase or decrease the DNA bending rate were studied, including one with a nick in a critical kink position, and a substrate with longer DNA ends to increase hydrodynamic friction during bending. Kinetic experiments were carried out under pseudofirst-order conditions, in which the protein concentration is in substantial excess over DNA. At lower concentrations, the reaction rate rises linearly with protein concentration, implying rate limitation by the bimolecular reaction step. At high concentrations the rate reaches a plateau value, which strongly depends on temperature and the nature of the DNA substrate. We ascribe this reaction limit to the DNA bending rate and propose that complex formation is sequential at high concentration: IHF binds rapidly to DNA, followed by slower DNA bending. Our observations on the bending step kinetics are in agreement with results using the temperature-jump kinetic method.

Mechanism of DNA Recognition at a Viral Replication Origin

Recognition of the DNA origin by the Epstein-Barr nuclear antigen 1 (EBNA1) protein is the primary event in latentphase genome replication of the Epstein-Barr virus, a model for replication initiation in eukaryotes. We carried out an extensive thermodynamic and kinetic characterization of the binding mechanism of the DNA binding domain of EBNA1, EBNA1452-641, to a DNA fragment containing a single specific origin site. The interaction displays a binding energy of 12.7 kcal mol-1, with 11.9 kcal mol-1 coming from the enthalpic change with a minimal entropic contribution. Formation of the EBNA1452-641.DNA complex is accompanied by a heat capacity change of -1.22 kcal mol-1 K-1, a very large value considering the surface area buried, which we assign to an unusually apolar protein-DNA interface. Kinetic dissociation experiments, including fluorescence anisotropy and a continuous native electrophoretic mobility shift assay, confirmed that two EBNA1.DNA complex conformers are in slow equilibrium; one dissociates slowly (t1/2 approximately 41 min) through an undissociated intermediate species and the other corresponds to a fast twostep dissociation route (t1/2 approximately 0.8 min). In line with this, at least two parallel association events from two populations of protein conformers are observed, with on-rates of 0.25-1.6x10(8) m-1 s-1, which occur differentially either in excess protein or DNA molecules. Both parallel complexes undergo subsequent firstorder rearrangements of approximately 2.0 s-1 to yield two consolidated complexes. These parallel association and dissociation routes likely allow additional flexible regulatory events for site recognition depending on site availability according to nucleus environmental conditions, which may lock a final recognition event, dissociate and re-bind, or slide along the DNA.

The recognition of local DNA conformation by the human papillomavirus type 6 E2 protein

The E2 proteins are transcription/replication factors from papillomaviruses. Human papillomaviruses (HPVs) can be broadly divided in two groups; low-risk HPV subtypes cause benign warts while high-risk HPVs give rise to cervical cancer. Although a range of crystal structures of E2 DNA-binding domains (DBD) from both high- and low-risk HPV subtypes have been reported previously, structures of E2 DBD:DNA complexes have only been available for high-risk HPV18 and bovine papillomavirus (BPV1). In the present study we report the unliganded and DNA complex structures of the E2 DBD from the low-risk HPV6. As in the previous E2-DNA structures, complex formation results in considerable bending of the DNA, which is facilitated by sequences with A:T-rich spacers that adopt a pre-bent conformation. The low-risk HPV6 E2-DNA complex differs from the earlier structures in that minimal deformation of the protein accompanies complex formation. Stopped-flow kinetic studies confirm that both high- and low-risk E2 proteins adapt their structures on binding to DNA, although this is achieved more readily for HPV6 E2. It therefore appears that the higher selectivity of the HPV6 E2 protein may arise from its limited molecular adaptability, a property that might distinguish the behaviour of E2 proteins from high- and low-risk HPV subtypes.

Structural and thermodynamic basis for the enhanced transcriptional control by the human papillomavirus strain-16 E2 protein

Strain 16 of the human papillomavirus is responsible for the largest number of cases of cervical cancers linked to this virus, and the E2 protein is the transcriptional regulator of all viral genes. We present the first structure for the DNA binding domain of HPV16 E2 bound to DNA, and in particular, to a natural cognate sequence. The NMR structure of the protein backbone reveals that the overall conformation remains virtually unchanged, and chemical shift analysis of the protein bound to a shorter DNA duplex uncovered a contact out of the minimal E2 DNA binding site, made by lysine 349. This contact was confirmed by titration calorimetry and mutagenesis, with a contribution of 1.0 kcal mol(-)(1) to binding energy. HPV16 E2 has the highest DNA binding affinity and exerts a strict transcriptional control, translated into the repression of the E6 and E7 oncogenes. These novel features provide the structural and thermodynamic basis for this tight transcriptional control, the loss of which correlates with carcinogenesis.

WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

Water is essential for life in many ways, and without it biomolecules might no longer truly be biomolecules. In particular, water is important to the structure, stability, dynamics, and function of biological macromolecules. In protein folding, water mediates the collapse of the chain and the search for the native topology through a funneled energy landscape. Water actively participates in molecular recognition by mediating the interactions between binding partners and contributes to either enthalpic or entropic stabilization. Accordingly, water must be included in recognition and structure prediction codes to capture specificity. Thus water should not be treated as an inert environment, but rather as an integral and active component of biomolecular systems, where it has both dynamic and structural roles. Focusing on water sheds light on the physics and function of biological machinery and self-assembly and may advance our understanding of the natural design of proteins and nucleic acids.

Specificity in DNA recognition by a peptide from papillomavirus E2 protein

The E2 proteins of papillomavirus specifically bind to double-stranded DNA containing the consensus sequence ACCG-N4-CGGT, where N is any nucleotide. Here, we show the binding and recognition of dissimilar DNA sequences by an 18 amino-acid peptide (alpha1E2), which corresponds to the DNA-recognition helix, alpha-helix-1. Isothermal DNA binding assays performed with the DNA consensus sequence show saturable curves with alpha1E2 peptide, and the alpha1E2 peptide is converted to an ordered conformation upon complexation. Measurements performed with non-specific DNA sequence fail to saturate, a behavior characteristic of non-specific binding. Binding of the alpha1E2 peptide to these DNA sequences display a different counter-ion dependence, indicating a dissimilar, sequence-dependent mechanism of interaction. Quantitative stoichiometric measurements revealed the specificity in alpha1E2 peptide recognition of the ACCG half-site, demonstrating capacity for discrimination of nucleic acid bases sequences without the need of a whole protein architecture.

Electrostatic potential of aminoacyl-tRNA synthetase navigates tRNA on its pathway to the binding site

In the first stage of a diffusion-controlled enzymatic reaction, aminoacyl-tRNA synthetases (aaRSs) interact with cognate tRNAs forming non-specific encounters. The aaRSs catalyzing the same overall aminoacylation reaction vary greatly in subunit organization, structural domain composition and amino acid sequence. The diffusional association of aaRS and tRNA was found to be governed by long-range electrostatic interactions when the homogeneous negative potential of tRNA fits to the patches of positive potential produced by aaRS; one patch for each tRNA substrate molecule. Considering aaRS as a molecule with anisotropic reactivity and on the basis of continuum electrostatics and Smoluchowski's theory, the reaction conditions for tRNA-aaRS diffusional encounters were formulated. The domains, categorized as enzymatically relevant, appeared to be non-essential for field sculpturing at long distances. On the other hand, a set of complementary domains exerts primary control on the aaRS isopotential surface formation. Subdividing the aaRS charged residues into native, conservative and non-conservative subsets, we evaluated the contribution of each group to long-range electrostatic potential. Surprisingly, the electrostatic potential landscapes generated by native and non-conservative subsets are fairly similar, thus suggesting the non-conservative subset is developed specifically for efficient tRNA attraction.

Free energy contributions to direct readout of a DNA sequence

The energetic contributions of individual DNA-contacting side chains to specific DNA recognition in the human papillomavirus 16 E2C-DNA complex is small (less than 1.0 kcal mol(-1)), independent of the physical and chemical nature of the interaction, and is strictly additive. The sum of the individual contributions differs 1.0 kcal mol(-1) from the binding energy of the wild-type protein. This difference corresponds to the contribution from the deformability of the DNA, known as "indirect readout." Thus, we can dissect the energetic contribution to DNA binding into 90% direct and 10% indirect readout components. The lack of high energy interactions indicates the absence of "hot spots," such as those found in protein-protein interfaces. These results are compatible with a highly dynamic and "wet" protein-DNA interface, yet highly specific and tight, where individual interactions are constantly being formed and broken.

Positive contribution of hydration on DNA binding by E2c protein from papillomavirus

Protein-nucleic acid interactions are responsible for the regulation of key biological events such as genomic transcription and recombination and viral replication. However, the recognition mechanisms involved in these processes are not completely understood. Here, we investigate the dominant forces involved in protein-protein and protein-DNA interactions for the 80-amino-acid C-terminal domain of the E2 protein (E2c) from human papillomavirus (HPV-16). The E2c protein is a homodimer that specifically binds to double-stranded DNA containing the consensus sequence ACCG-N(4)-CGGT, where N is any nucleotide. DNA binding affinity is reduced by lowering water chemical potential, accompanied by an increase in cooperativity. Wyman linkage relations between affinity and water chemical potential indicate that 11 additional water molecules are bound in the formation of the complex between E2c and DNA. Salt dissociation isotherms showed that 10 counterions are released upon association, even at low water activity, indicating that this latter variable does not change the electrostatic component of the interaction. Further analysis demonstrates a strong dependence of cooperativity of binding on the protein concentration. Altogether, these results reveal a novel binding pathway in which the consolidated complex may achieve its final form via a monomer-DNA intermediate, which favors the binding of a second monomer. This molecular mechanism reveals the contributions of multiple conformers in a tight virus genome modulation that seems to be important in the cell infection scenario.

Monomeric and dimeric bZIP transcription factor GCN4 bind at the same rate to their target DNA site

Basic leucine zipper (bZIP) transcription factors are dimeric proteins that recognize dyadic and mostly palindromic DNA sites. Dimerization of bZIP transcription factor GCN4 is linked to the folding of its C-terminal leucine zipper domain. However, monomeric GCN4, lacking a folded leucine zipper, also recognizes the DNA site with dimerization taking place on the DNA. Here we report the kinetics of DNA recognition by unfolded monomeric and folded dimeric derivatives of GCN4 using a 19 bp double-stranded DNA containing a palindromic CRE site. The rate of DNA binding of both monomeric and dimeric GCN4 has a bimolecular rate constant of 3-5 x 10(8) M(-1) s(-1), which is near the diffusion limit. Because the rate of dimerization of GCN4 is slower (1.7 x 10(7) M(-1) s(-1)) than the rate of DNA association, the formation of the dimeric GCN4-DNA complex through consecutive binding of two monomers (monomer pathway) is faster when starting from free monomers. Thus, the results presented here support facilitated and rapid target recognition by the monomeric transcription factor. However, DNA binding of preformed folded dimeric GCN4 is as rapid as complex formation through the monomer pathway. Therefore, the monomer and dimer pathways are kinetically equivalent if monomeric and dimeric GCN4 are at equilibrium. Hence, the dimer pathway may also have a role under in vivo conditions. The lower affinity of GCN4 in which two DNA contacting residues have been mutated is due exclusively to the faster dissociation of the mutant protein-DNA complex and not to slower complex formation.

Predicting indirect readout effects in protein-DNA interactions

Recognition of DNA by proteins relies on direct interactions with specific DNA-functional groups, along with indirect effects that reflect variable energetics in the response of DNA sequences to twisting and bending distortions induced by proteins. Predicting indirect readout requires knowledge of the variations in DNA curvature and flexibility in the affected region, which we have determined for a series of DNA-binding sites for the E2 regulatory protein by using the cyclization kinetics method. We examined 16 sites containing different noncontacted spacer sequences, which vary by more than three orders of magnitude in binding affinity. For 15 of these sites, the variation in affinity was predicted within a factor of 3, by using experimental curvature and flexibility values and a statistical mechanical theory. The sole exception was traced to differential magnesium ion binding.