Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.     

Dynamics may significantly influence the estimation of interatomic distances in biomolecular X-ray structures

Atomic positions obtained by X-ray crystallography are time and space averages over many molecules in the crystal. Importantly, interatomic distances, calculated between such average positions and frequently used in structural and mechanistic analyses, can be substantially different from the more appropriate time-average and ensemble-average interatomic distances. Using crystallographic B-factors, one can deduce corrections, which have so far been applied exclusively to small molecules, to obtain correct average distances as a function of the type of atomic motion. Here, using 4774 high-quality protein X-ray structures, we study the significance of such corrections for different types of atomic motion. Importantly, we show that for distances shorter than 5 Å, corrections greater than 0.5 Å may apply, especially for noncorrelated or anticorrelated motion. For example, 14% of the studied structures have at least one pair of atoms with a correction of ≥ 0.5 Å in the case of noncorrelated motion. Using molecular dynamics simulations of villin headpiece, ubiquitin, and SH3 domain unit cells, we demonstrate that the majority of average interatomic distances in these proteins agree with noncorrelated corrections, suggesting that such deviations may be truly relevant. Importantly, we demonstrate that the corrections do not significantly affect stereochemistry and the overall quality of final refined X-ray structures, but can provide marked improvements in starting unrefined models obtained from low-resolution X-ray data. Finally, we illustrate the potential mechanistic and biological significance of the calculated corrections for KcsA ion channel and show that they provide indirect evidence that motions in its selectivity filter are highly correlated.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

Role of the Adenine Ligand on the Stabilization of the Secondary and Tertiary Interactions in the Adenine Riboswitch

Riboswitches are RNA-based genetic control elements that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket. To facilitate an atomic detail interpretation of experimental investigations on the role of the adenine ligand on the conformational properties and kinetics of folding of the add adenine riboswitch, we performed molecular dynamics simulations in both the presence and the absence of the ligand. In the absence of ligand, structural deviations were observed in the J23 junction and the P1 stem. Destabilization of the P1 stem in the absence of ligand involves the loss of direct stabilizing interactions with the ligand, with additional contributions from the J23 junction region. The J23 junction of the riboswitch is found to be more flexible, and the tertiary contacts among the junction regions are altered in the absence of the adenine ligand; results suggest that the adenine ligand associates and dissociates from the riboswitch in the vicinity of J23. Good agreement was obtained with the experimental data with the results indicating dynamic behavior of the adenine ligand on the nanosecond time scale to be associated with the dynamic behavior of hydrogen bonding with the riboswitch. Results also predict that direct interactions of the adenine ligand with U74 of the riboswitch are not essential for stable binding although it is crucial for its recognition. The possibility of methodological artifacts and force-field inaccuracies impacting the present observations was checked by additional molecular dynamics simulations in the presence of 2,6-diaminopurine and in the crystal environment.     

NMR Structure of the Amino-Terminal Domain of the Lambda Integrase Protein in Complex with DNA: Immobilization of a Flexible Tail Facilitates Beta-Sheet Recognition of the Major Groove

The integrase protein (Int) from bacteriophage lambda is the archetypal member of the tyrosine recombinase family, a large group of enzymes that rearrange DNA in all domains of life. Int catalyzes the insertion and excision of the viral genome into and out of the Escherichia coli chromosome. Recombination transpires within higher-order nucleoprotein complexes that form when its amino-terminal domain binds to arm-type DNA sequences that are located distal to the site of strand exchange. Arm-site binding by Int is essential for catalysis, as it promotes Int-mediated bridge structures that stabilize the recombination machinery. We have elucidated how Int is able to sequence specifically recognize the arm-type site sequence by determining the solution structure of its amino-terminal domain (Int^N, residues Met1 to Leu64) in complex with its P′2 DNA binding site. Previous studies have shown that Int^N adopts a rare monomeric DNA binding fold that consists of a three-stranded antiparallel beta-sheet that is packed against a carboxy-terminal alpha helix. A low-resolution crystal structure of the full-length protein also revealed that the sheet is inserted into the major groove of the arm-type site. The solution structure presented here reveals how Int^N specifically recognizes the arm-type site sequence. A novel feature of the new solution structure is the use of an 11-residue tail that is located at the amino terminus. DNA binding induces the folding of a 3₁₀ helix in the tail that projects the amino terminus of the protein deep into the minor groove for stabilizing DNA contacts. This finding reveals the structural basis for the observation that the “unstructured” amino terminus is required for recombination.     

Backbone dynamics in an intramolecular prolylpeptide-SH3 complex from the diphtheria toxin repressor, DtxR

The diphtheria toxin repressor contains an SH3-like domain that forms an intramolecular complex with a proline-rich (Pr) peptide segment and stabilizes the inactive state of the repressor. Upon activation of diphtheria toxin repressor (DtxR) by transition metals, this intramolecular complex must dissociate as the SH3 domain and Pr segment form different interactions in the active repressor. Here we investigate the dynamics of this intramolecular complex using backbone amide nuclear spin relaxation rates determined using NMR spectroscopy and molecular dynamics trajectories. The SH3 domain in the unbound and bound states showed typical dynamics in that the secondary structures were fairly ordered with high generalized order parameters and low effective correlation times, while residues in the loops connecting β-strands exhibited reduced generalized order parameters and required additional motional terms to adequately model the relaxation rates. Residues forming the Pr segment exhibited low-order parameters with internal rotational correlation times on the order of 0.6 ns-1 ns. Further analysis showed that the SH3 domain was rich in millisecond time scale motions while the Pr segment exhibited motions on the 100 μs time scale. Molecular dynamics simulations indicated structural rearrangements that may contribute to the observed relaxation rates and, together with the observed relaxation rate data, suggested that the Pr segment exhibits a binding ↔ unbinding equilibrium. The results here provide new insights into the nature of the intramolecular complex and provide a better understanding of the biological role of the SH3 domain in regulating DtxR activity.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

A new structural domain in the Escherichia coli RcsC hybrid sensor kinase connects histidine kinase and phosphoreceiver domains

The Rcs signalling pathway controls a variety of physiological functions like capsule synthesis, cell division or motility in prokaryotes. The Rcs regulation cascade, involving a multi-step phosphorelay between the two membrane-bound hybrid sensor kinases RcsC and RcsD and the global regulator RcsB, is, up to now, one of the most complicated regulatory systems in bacteria. To understand the structural basis of Rcs signal transduction, NMR spectroscopy was employed to determine the solution structure of the RcsC C terminus, possessing a phosphoreceiver domain (RcsC-PR), and a region previously described as a long linker between the histidine kinase domain of RcsC (RcsC-HK) and the RcsC-PR. We have found that the linker region comprises an independent structural domain of a new alpha/beta organization, which we named RcsC-ABL domain (Alpha/Beta/Loop). The ABL domain appears to be a conserved and unique structural element of RcsC-like kinases with no significant sequence homology to other proteins. The second domain of the C terminus, the RcsC-PR domain, represents a well-folded CheY-like phosphoreceiver domain with the central parallel beta-sheet covered with two alpha-helical layers on both sides. We have mapped the interaction of RcsC-ABL and RcsC-PR with the histidine phosphotransfer domain (HPt) of RcsD. In addition we have characterized the interaction with and the conformational effects of Mg2+ and the phosphorylation mimetic BeF(-)(3) on RcsC-ABL and RcsC-PR.     

Horse Prion Protein NMR Structure and Comparisons with Related Variants of the Mouse Prion Protein

The NMR structure of the horse (Equus caballus) cellular prion protein at 25 °C exhibits the typical PrP^C [cellular form of prion protein (PrP)] global architecture, but in contrast to most other mammalian PrP^Cs, it contains a well-structured loop connecting the β2 strand with the α2 helix. Comparison with designed variants of the mouse prion protein resulted in the identification of a single amino acid exchange within the loop, D167S, which correlates with the high structural order of this loop in the solution structure at 25 °C and is unique to the PrP sequences of equine species. The β2-α2 loop and the α3 helix form a protein surface epitope that has been proposed to be the recognition area for a hypothetical chaperone, “protein X,” which would promote conversion of PrP^C into the disease-related scrapie form and thus mediate intermolecular interactions related to the transmission barrier for transmissible spongiform encephalopathies (TSEs) between different species. The present results are evaluated in light of recent indications from in vivo experiments that the local β2-α2 loop structure affects the susceptibility of transgenic mice to TSEs and the fact that there are no reports on TSE in horses.     

Lactose Binding to Galectin-1 Modulates Structural Dynamics, Increases Conformational Entropy, and Occurs with Apparent Negative Cooperativity

Galectins are a family of lectins with a conserved carbohydrate recognition domain that interacts with β-galactosides. By binding cell surface glycoconjugates, galectin-1 (gal-1) is involved in cell adhesion and migration processes and is an important regulator of tumor angiogenesis. Here, we used heteronuclear NMR spectroscopy and molecular modeling to investigate lactose binding to gal-1 and to derive solution NMR structures of gal-1 in the lactose-bound and unbound states. Structure analysis shows that the β-strands and loops around the lactose binding site, which are more open and dynamic in the unbound state, fold in around the bound lactose molecule, dampening internal motions at that site and increasing motions elsewhere throughout the protein to contribute entropically to the binding free energy. CD data support the view of an overall more open structure in the lactose-bound state. Analysis of heteronuclear single quantum coherence titration binding data indicates that lactose binds the two carbohydrate recognition domains of the gal-1 dimer with negative cooperativity, in that the first lactose molecule binds more strongly (K₁ = 21 ± 6 × 10³ M^− 1) than the second (K₂ = 4 ± 2 × 10³ M^− 1). Isothermal calorimetry data fit using a sequential binding model present a similar picture, yielding K₁ = 20 ± 10 × 10³ M^− 1 and K₂ = 1.67 ± 0.07 × 10³ M^− 1. Molecular dynamics simulations provide insight into structural dynamics of the half-loaded lactose state and, together with NMR data, suggest that lactose binding at one site transmits a signal through the β-sandwich and loops to the second binding site. Overall, our results provide new insight into gal-1 structure-function relationships and to protein-carbohydrate interactions in general.     

Myosin binding protein C positioned to play a key role in regulation of muscle contraction: structure and interactions of domain C1

Myosin binding protein C (MyBP-C) is a thick filament protein involved in the regulation of muscle contraction. Mutations in the gene for MyBP-C are the second most frequent cause of hypertrophic cardiomyopathy. MyBP-C binds to myosin with two binding sites, one at its C-terminus and another at its N-terminus. The N-terminal binding site, consisting of immunoglobulin domains C1 and C2 connected by a flexible linker, interacts with the S2 segment of myosin in a phosphorylation-regulated manner. It is assumed that the function of MyBP-C is to act as a tether that fixes the S1 heads in a resting position and that phosphorylation releases the S1 heads into an active state. Here, we report the structure and binding properties of domain C1. Using a combination of site-directed mutagenesis and NMR interaction experiments, we identified the binding site of domain C1 in the immediate vicinity of the S1-S2 hinge, very close to the light chains. In addition, we identified a zinc binding site on domain C1 in close proximity to the S2 binding site. Its zinc binding affinity (K_d of approximately 10-20 μM) might not be sufficient for a physiological effect. However, the familial hypertrophic cardiomyopathy-related mutation of one of the zinc ligands, glutamine 210 to histidine, will significantly increase the binding affinity, suggesting that this mutation may affect S2 binding. The close proximity of the C1 binding site to the hinge, the light chains and the S1 heads also provides an explanation for recent observations that (a) shorter fragments of MyBP-C unable to act as a tether still have an effect on the actomyosin ATPase and (b) as to why the myosin head positions in phosphorylated wild-type mice and MyBP-C knockout mice are so different: Domain C1 bound to the S1-S2 hinge is able to manipulate S1 head positions, thus influencing force generation without tether. The potentially extensive extra interactions of C1 are expected to keep it in place, while phosphorylation dislodges the C1-C2 linker and domain C2. As a result, the myosin heads would always be attached to a tether that has phosphorylation-dependent length regulation.     

Predicting the effect of a point mutation on a protein fold: the villin and advillin headpieces and their Pro62Ala mutants

Homology modeling of unknown proteins is based on the assumption that highly similar sequences are likely to share the same fold. However, this does not provide any information on the stability of a given fold, which is ultimately determined by the subtle interplay of enthalpic and entropic contributions. Herein it is shown that ab initio atomistic simulations can be used to predict the effect of point mutations on the stability of a protein fold. The calculations indicate that the fold stabilities of two proteins of similar sequence and identical fold, the villin and advillin C-terminal headpiece fragments, are different and that the same P62A point mutation has a dramatic effect on the fold of villin but a minor one on that of advillin. These predictions were subsequently validated by NMR and CD experiments.     

The Conformational Properties of the Glc3Man Unit Suggest Conformational Biasing within the Chaperone-assisted Glycoprotein Folding Pathway

A major puzzle is: are all glycoproteins routed through the ER calnexin pathway irrespective of whether this is required for their correct folding? Calnexin recognizes the terminal Glcα1-3Manα linkage, formed by trimming of the Glcα1-2Glcα1-3Glcα1-3Manα (Glc₃Man) unit in Glc₃Man₉GlcNAc₂. Different conformations of this unit have been reported. We have addressed this problem by studying the conformation of a series of N-glycans; i.e. Glc₃ManOMe, Glc₃Man₄,₅,₇GlcNAc₂ and Glc₁Man₉GlcNAc₂ using 2D NMR NOESY, ROESY, T-ROESY and residual dipolar coupling experiments in a range of solvents, along with solution molecular dynamics simulations of Glc₃ManOMe. Our results show a single conformation for the Glcα1-2Glcα and Glcα1-3Glcα linkages, and a major (65%) and a minor (30%) conformer for the Glcα1-3Manα linkage. Modeling of the binding of Glc₁Man₉GlcNAc₂ to calnexin suggests that it is the minor conformer that is recognized by calnexin. This may be one of the mechanisms for controlling the rate of recruitment of proteins into the calnexin/calreticulin chaperone system and enabling proteins that do not require such assistance for folding to bypass the system. This is the first time evidence has been presented on glycoprotein folding that suggests the process may be optimized to balance the chaperone-assisted and chaperone-independent pathways.     

Structure of the Flexible Amino-Terminal Domain of Prion Protein Bound to a Sulfated Glycan

The intrinsically disordered amino-proximal domain of hamster prion protein (PrP) contains four copies of a highly conserved octapeptide sequence, PHGGGWGQ, that is flanked by two polycationic residue clusters. This N-terminal domain mediates the binding of sulfated glycans, which can profoundly influence the conversion of PrP to pathological forms and the progression of prion disease. To investigate the structural consequences of sulfated glycan binding, we performed multidimensional heteronuclear (¹H, ¹³C, ¹⁵N) NMR (nuclear magnetic resonance), circular dichroism (CD), and fluorescence studies on hamster PrP residues 23-106 (PrP 23-106) and fragments thereof when bound to pentosan polysulfate (PPS). While the majority of PrP 23-106 remain disordered upon PPS binding, the octarepeat region adopts a repeating loop-turn structure that we have determined by NMR. The β-like turns within the repeats are corroborated by CD data demonstrating that these turns are also present, although less pronounced, without PPS. Binding to PPS exposes a hydrophobic surface composed of aligned tryptophan side chains, the spacing and orientation of which are consistent with a self-association or ligand binding site. The unique tryptophan motif was probed by intrinsic tryptophan fluorescence, which displayed enhanced fluorescence of PrP 23-106 when bound to PPS, consistent with the alignment of tryptophan side chains. Chemical-shift mapping identified binding sites on PrP 23-106 for PPS, which include the octarepeat histidine and an N-terminal basic cluster previously linked to sulfated glycan binding. These data may in part explain how sulfated glycans modulate PrP conformational conversions and oligomerizations.     

NMR Structure of the N-Terminal Domain of Capsid Protein from the Mason-Pfizer Monkey Virus

The high-resolution structure of the N-terminal domain (NTD) of the retroviral capsid protein (CA) of Mason-Pfizer monkey virus (M-PMV), a member of the betaretrovirus family, has been determined by NMR. The M-PMV NTD CA structure is similar to the other retroviral capsid structures and is characterized by a six α-helix bundle and an N-terminal β-hairpin, stabilized by an interaction of highly conserved residues, Pro1 and Asp57. Since the role of the β-hairpin has been shown to be critical for formation of infectious viral core, we also investigated the functional role of M-PMV β-hairpin in two mutants (i.e., ΔP1NTDCA and D57ANTDCA) where the salt bridge stabilizing the wild-type structure was disrupted. NMR data obtained for these mutants were compared with those obtained for the wild type. The main structural changes were observed within the β-hairpin structure; within helices 2, 3, and 5; and in the loop connecting helices 2 and 3. This observation is supported by biochemical data showing different cleavage patterns of the wild-type and the mutated capsid-nucleocapsid fusion protein (CANC) by M-PMV protease. Despite these structural changes, the mutants with disrupted salt bridge are still able to assemble into immature, spherical particles. This confirms that the mutual interaction and topology within the β-hairpin and helix 3 might correlate with the changes in interaction between immature and mature lattices.     

NMR analysis of the conformational properties of the trapped on-pathway folding intermediate of the bacterial immunity protein Im7

Previous work shows that the transiently populated, on-pathway intermediate in Im7 folding contains three of the four native alpha-helices docked around a core stabilised by native and non-native interactions. To determine the structure and dynamic properties of this species in more detail, we have used protein engineering to trap the intermediate at equilibrium and analysed the resulting proteins using NMR spectroscopy and small angle X-ray scattering. Four variants were created. In L53AI54A, two hydrophobic residues within helix III are truncated, preventing helix III from docking stably onto the developing hydrophobic core. In two other variants, the six residues encompassing the native helix III were replaced with three (H3G3) or six (H3G6) glycine residues. In the fourth variant, YY, two native tyrosine residues (Tyr55 and Tyr56) were re-introduced into H3G6 to examine their role in determining the properties of the intermediate ensemble. All four variants show variable peak intensities and broad peak widths, consistent with these proteins being conformationally dynamic. Chemical shift analyses demonstrated that L53AI54A and YY contain native-like secondary structure in helices I and IV, while helix II is partly formed and helix III is absent. Lack of NOEs and rapid NH exchange for L53AI54A, combined with detailed analysis of the backbone dynamics, indicated that the hydrophobic core of this variant is not uniquely structured, but fluctuates on the NMR timescale. The results demonstrate that though much of the native-like secondary structure of Im7 is present in the variants, their hydrophobic cores remain relatively fluid. The comparison of H3G3/H3G6 and L53AI54A/YY suggests that Tyr55 and/or Tyr56 interact with the three-helix core, leading other residues in this region of the protein to dock with the core as folding progresses. In this respect, the three-helix bundle acts as a template for formation of helix III and the creation of the native fold.     

Structural and Functional Characterization of the NHR1 Domain of the Drosophila Neuralized E3 Ligase in the Notch Signaling Pathway

The Notch signaling pathway is critical for many developmental processes and requires complex trafficking of both Notch receptor and its ligands, Delta and Serrate. In Drosophila melanogaster, the endocytosis of Delta in the signal-sending cell is essential for Notch receptor activation. The Neuralized protein from D. melanogaster (Neur) is a ubiquitin E3 ligase, which binds to Delta through its first neuralized homology repeat 1 (NHR1) domain and mediates the ubiquitination of Delta for endocytosis. Tom, a Bearded protein family member, inhibits the Neur-mediated endocytosis through interactions with the NHR1 domain. We have identified the domain boundaries of the novel NHR1 domain, using a screening system based on our cell-free protein synthesis method, and demonstrated that the identified Neur NHR1 domain had binding activity to the 20-residue peptide corresponding to motif 2 of Tom by isothermal titration calorimetry experiments. We also determined the solution structure of the Neur NHR1 domain by heteronuclear NMR methods, using a ¹⁵N/¹³C-labeled sample. The Neur NHR1 domain adopts a characteristic β-sandwich fold, consisting of a concave five-stranded antiparallel β-sheet and a convex seven-stranded antiparallel β-sheet. The long loop (L6) between the β6 and β7 strands covers the hydrophobic patch on the concave β-sheet surface, and the Neur NHR1 domain forms a compact globular fold. Intriguingly, in spite of the slight, but distinct, differences in the topology of the secondary structure elements, the structure of the Neur NHR1 domain is quite similar to those of the B30.2/SPRY domains, which are known to mediate specific protein-protein interactions. Further NMR titration experiments of the Neur NHR1 domain with the 20-residue Tom peptide revealed that the resonances originating from the bottom area of the β-sandwich (the L3, L5, and L11 loops, as well as the tip of the L6 loop) were affected. In addition, a structural comparison of the Neur NHR1 domain with the first NHR domain of the human KIAA1787 protein, which is from another NHR subfamily and does not bind to the 20-residue Tom peptide, suggested the critical amino acid residues for the interactions between the Neur NHR1 domain and the Tom peptide. The present structural study will shed light on the role of the Neur NHR1 domain in the Notch signaling pathway.     

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 β-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 β-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 β-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 β-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.