Ignacio E. Sánchez

Evidence for Sequential Barriers and Obligatory Intermediates in Apparent Two-state Protein Folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.

Hammond behavior versus ground state effects in protein folding: evidence for narrow free energy barriers and residual structure in unfolded states.

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in phi-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.

Non-linear rate-equilibrium free energy relationships and Hammond behavior in protein folding

Non-linear rate-equilibrium relationships upon mutation or changes in solvent conditions are frequently observed in protein folding reactions and are usually interpreted in terms of Hammond behavior. Here we first give a general overview over the concept of transition state movements in chemical reactions and discuss its application to protein folding. We then show examples for genuine Hammond behavior and for apparent transition state movements caused by other effects like changes in the rate-limiting step of the folding reaction or ground state effects, i.e. structural changes in either the native state or the unfolded state. These examples show that apparent transition state movements can easily be mistaken for Hammond behavior. We describe experimental tests using self- and cross-interaction parameters to distinguish between structural changes in a single transition state following Hammond behavior and apparent transition state movements caused by other effects.

Targeting mechanism of the retinoblastoma tumor suppressor by a prototypical viral oncoprotein. Structural modularity, intrinsic disorder and phosphorylation of human papillomavirus E7.

DNA tumor viruses ensure genome amplification by hijacking the cellular replication machinery and forcing infected cells to enter the S phase. The retinoblastoma (Rb) protein controls the G1/S checkpoint, and is targeted by several viral oncoproteins, among these the E7 protein from human papillomaviruses (HPVs). A quantitative investigation of the interaction mechanism between the HPV16 E7 protein and the RbAB domain in solution revealed that 90% of the binding energy is determined by the LxCxE motif, with an additional binding determinant (1.0u2003kcal·mol−1) located in the C‐terminal domain of E7, establishing a dual‐contact mode. The stoichiometry and subnanomolar affinity of E7 indicated that it can bind RbAB as a monomer. The low‐risk HPV11 E7 protein bound 2.0 kcal·mol−1 more weakly than the high‐risk HPV16 and HPV18 type counterparts, but the modularity and binding mode were conserved. Phosphorylation at a conserved casein kinase II site in the natively unfolded N‐terminal domain of E7 affected the local conformation by increasing the polyproline II content and stabilizing an extended conformation, which allowed for a tighter interaction with the Rb protein. Thus, the E7–RbAB interaction involves multiple motifs within the N‐terminal domain of E7 and at least two conserved interaction surfaces in RbAB. We discussed a mechanistic model of the interaction of the Rb protein with a viral target in solution, integrated with structural data and the analysis of other cellular and viral proteins, which provided information about the balance of interactions involving the Rb protein and how these determine the progression into either the normal cell cycle or transformation.

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.

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.

Detecting repetitions and periodicities in proteins by tiling the structural space.

The notion of energy landscapes provides conceptual tools for understanding the complexities of protein folding and function. Energy landscape theory indicates that it is much easier to find sequences that satisfy the “Principle of Minimal Frustration” when the folded structure is symmetric (Wolynes, P. G. Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci. U.S.A.1996, 93, 14249–14255). Similarly, repeats and structural mosaics may be fundamentally related to landscapes with multiple embedded funnels. Here we present analytical tools to detect and compare structural repetitions in protein molecules. By an exhaustive analysis of the distribution of structural repeats using a robust metric, we define those portions of a protein molecule that best describe the overall structure as a tessellation of basic units. The patterns produced by such tessellations provide intuitive representations of the repeating regions and their association toward higher order arrangements. We find that some protein architectures can be described as nearly periodic, while in others clear separations between repetitions exist. Since the method is independent of amino acid sequence information, we can identify structural units that can be encoded by a variety of distinct amino acid sequences.

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.

Indirect DNA Readout on the Protein Side: Coupling between Histidine Protonation, Global Structural Cooperativity, Dynamics, and DNA Binding of the Human Papillomavirus Type 16 E2C Domain

DNA sequence recognition by the homodimeric C-terminal domain of the human papillomavirus type 16 E2 protein (E2C) is known to involve both direct readout and DNA-dependent indirect readout mechanisms, while protein-dependent indirect readout has been deduced but not directly observed. We have investigated coupling between specific DNA binding and the dynamics of the unusual E2C fold, using pH as an external variable. Nuclear magnetic resonance and isothermal titration calorimetry show that pH titration of His318 in the complex interface and His288 in the core of the domain is coupled to both binding and the dynamics of the beta-barrel core of E2C, with a tradeoff between dimer stability and function. Specific DNA binding is, in turn, coupled to the slow dynamics and amide hydrogen exchange in the entire beta-barrel, reaching residues far apart from the DNA recognition elements but not affecting the two helices of each monomer. The changes are largest in the dimerization interface, suggesting that the E2C beta-barrel acts as a hinge that regulates the relative position of the DNA recognition helices. In conclusion, the cooperative dynamics of the human papillomavirus type 16 E2C beta-barrel is coupled to sequence recognition in a protein-dependent indirect readout mechanism. The patterns of residue substitution in genital papillomaviruses support the importance of the protonation states of His288 and His318 and suggest that protein-dependent indirect readout and histidine pH titration may regulate DNA binding in the cell.

Cysteine-rich positions outside the structural zinc motif of human papillomavirus E7 provide conformational modulation and suggest functional redox roles

The E7 protein from high-risk human papillomavirus is essential for cell transformation in cervical, oropharyngeal, and other HPV-related cancers, mainly through the inactivation of the retinoblastoma (Rb) tumor suppressor. Its high cysteine content (~7%) and the observation that HPV-transformed cells are under oxidative stress prompted us to investigate the redox properties of the HPV16 E7 protein under biologically compatible oxidative conditions. The seven cysteines in HPV16 E7 remain reduced in conditions resembling the basal reduced state of a cell. However, under oxidative stress, a stable disulfide bridge forms between cysteines 59 and 68. Residue 59 has a protective effect on the other cysteines, and its mutation leads to an overall increase in the oxidation propensity of E7, including cysteine 24 central to the Rb binding motif. Gluthationylation of Cys 24 abolishes Rb binding, which is reversibly recovered upon reduction. Cysteines 59 and 68 are located 18.6 Å apart, and the formation of the disulfide bridge leads to a large structural rearrangement while retaining strong Zn association. These conformational and covalent changes are fully reversible upon restoration of the reductive environment. In addition, this is the first evidence of an interaction between the N-terminal intrinsically disordered and the C-terminal globular domains, known to be highly and separately conserved among human papillomaviruses. The significant conservation of such noncanonical cysteines in HPV E7 proteins leads us to propose a functional redox activity. Such an activity adds to the previously discovered chaperone activity of E7 and supports the picture of a moonlighting pathological role of this paradigmatic viral oncoprotein.