Aromatic-aromatic interactions in the zinc finger motif. Analysis of the two-dimensional nuclear magnetic resonance structure of a mutant domain

Temperature effect on the structure and conformational fluctuations in two zinc knuckles from the mouse mammary tumor virus

Zinc fingers are small protein domains in which zinc plays a structural role, contributing to the stability of the zinc-peptide complex. Zinc fingers are structurally diverse and are present in proteins that perform a broad range of functions in various cellular processes, such as replication and repair, transcription and translation, metabolism and signaling, cell proliferation, and apoptosis. Zinc fingers typically function as interaction modules and bind to a wide variety of compounds, such as nucleic acids, proteins, and small molecules. In this study, we investigated the structural properties, in solution, of the proximal and distal zinc knuckles of the nucleocapsid (NC) protein from the mouse mammary tumor virus (MMTV) (MMTV NC). For this purpose, we performed a series of molecular dynamics simulations in aqueous solution at 300 K, 333 K, and 348 K. The temperature effect was evaluated in terms of root mean square deviation of the backbone atoms and root mean square fluctuation of the coordinating residue atoms. The stability of the zinc coordination sphere was analyzed based upon the time profile of the interatomic distances between the zinc ions and the chelator atoms. The results indicate that the hydrophobic character of the proximal zinc finger is dominant at 333 K. The low mobility of the coordinating residues suggests that the strong electrostatic effect exerted by the zinc ion on its coordinating residues is not influenced by the increase in temperature. The evolution of the structural parameters of the coordination sphere of the distal zinc finger at 300 K gives us a reasonable picture of the unfolding pathway, as proposed by Bombarda and coworkers (Bombarda et al., 2005), which can predict the binding order of the four conserved ligand-binding residues. Our results support the conclusion that the structural features can vary significantly between the two zinc knuckles of MMTV NC.

Stability and folding behavior analysis of zinc-finger using simple models

Zinc-fingers play crucial roles in regulating gene expression and mediating protein-protein interactions. In this article, two different proteins (Sp1f2 and FSD-1) are investigated using the Gaussian network model and anisotropy elastic network model. By using these simple coarse-grained methods, we analyze the structural stabilization and establish the unfolding pathway of the two different proteins, in good agreement with related experimental and molecular dynamics simulation data. From the analysis, it is also found that the folding process of the zinc-finger motif is predominated by several factors. Both the zinc ion and C-terminal loop affect the folding pathway of the zinc-finger motif. Knowledge about the stability and folding behavior of zinc-fingers may help in understanding the folding mechanisms of the zinc-finger motif and in designing new zinc-fingers. Meanwhile, these simple coarse-grained analyses can be used as a general and quick method for mechanistic studies of metalloproteins.

Metal-coupled folding of Cys2His2 zinc-finger

Zinc-fingers, which widely exist in eukaryotic cell and play crucial roles in life processes, depend on the binding of zinc ion for their proper folding. To computationally study the zinc-coupled folding of the zinc-fingers, charge transfer and metal induced protonation/deprotonation effects have to be considered. Here, by attempting to implicitly account for such effects in classical molecular dynamics and performing intensive simulations with explicit solvent for the peptides with and without zinc binding, we investigate the folding of the Cys2His2-type zinc-finger motif and the coupling between the peptide folding and zinc binding. We find that zinc ion not only stabilizes the native structure but also participates in the whole folding process. It binds to the peptide at an early stage of folding and directs or modulates the folding and stabilizations of the component beta-hairpin and alpha-helix. Such a crucial role of zinc binding is mediated by the packing of the conserved hydrophobic residues. We also find that the packing of the hydrophobic residues and the coordination of the native ligands are coupled. Meanwhile, the processes of zinc binding, mis-ligation, ligand exchange, and zinc induced secondary structure conversion as well as the water behavior due to the involvement of zinc ion are characterized. Our results are in good agreement with related experimental observations and provide significant insight into the general mechanisms of the metal cofactor dependent protein folding and other metal-induced conformational changes of biological importance.

A Point Mutation in NEMO Associated with Anhidrotic Ectodermal Dysplasia with Immunodeficiency Pathology Results in Destabilization of the Oligomer and Reduces Lipopolysaccharide- and Tumor Necrosis Factor-mediated NF-κB Activation

The NEMO (NF-kappaB essential modulator) protein plays a crucial role in the canonical NF-kappaB pathway as the regulatory component of the IKK (IkappaB kinase) complex. The human disease anhidrotic ectodermal dysplasia with immunodeficiency (EDA-ID) has been recently linked to mutations in NEMO. We investigated the effect of an alanine to glycine substitution found in the NEMO polypeptide of an EDA-ID patient. This pathogenic mutation is located within the minimal oligomerization domain of the protein, which is required for the IKK activation in response to diverse stimuli. The mutation does not dramatically change the native-like state of the trimer, but temperature-induced unfolding studied by circular dichroism showed that it leads to an important loss in the oligomer stability. Furthermore, fluorescence studies showed that the tyrosine located in the adjacent zinc finger domain, which is possibly required for NEMO ubiquitination, exhibits an alteration in its spectral properties. This is probably due to a conformational change of this domain, providing evidence for a close interaction between the oligomerization domain and the zinc finger. In addition, functional complementation assays using NEMO-deficient pre-B and T lymphocytes showed that the pathogenic mutation reduced TNF-alpha and LPS-induced NF-kappaB activation by altering the assembly of the IKK complex. Altogether, our findings provide understanding as to how a single point mutation in NEMO leads to the observed EDA-ID phenotype in relation to the NEMO-dependent mechanism of IKK activation.

Contribution of individual zinc ligands to metal binding and peptide folding of zinc finger peptides

Little is known about the contribution of individual zinc-ligating amino acid residues for coupling between zinc binding and protein folding in zinc finger domains. To understand such roles of each zinc ligand, four zinc finger mutant peptides corresponding to the second zinc finger domain of Sp1 were synthesized. In the mutant peptides, glycine was substituted for one of four zinc ligands. Their metal binding and folding properties were spectroscopically characterized and compared to those of the native zinc finger peptide. In particular, the electronic charge-transfer and d-d bands of the Co(II)-substituted peptide complexes were used to examine the metal coordination number and geometry. Fluorescence emission studies revealed that the mutant peptides are capable of binding zinc despite removing one ligand. Circular dichroism results clearly showed the induction of an alpha-helix by zinc binding. In addition, the structures of certain mutant zinc finger peptides were simulated by molecular dynamics calculation. The information indicates that His23 and the hydrophobic core formed between the alpha-helix and the beta-sheet play an essential role in alpha-helix induction. This report demonstrates that each ligand does not contribute equally to alpha-helix formation and coordination geometry in the zinc finger peptide.

Metal-dependent folding and stability of nuclear hormone receptor DNA-binding domains

The nuclear/hormone receptors are an extensive family of ligand-activated transcription factors that recognise DNA targets through a highly conserved, structurally autonomous DNA-binding domain. The compact structure of the DNA-binding domain is supported by two zinc ions, each of which is co-ordinated by the tetrahedral arrangement of thiol groups from four cysteine residues. Metal binding is expected to be linked with deprotonation of the co-ordinating thiol groups and folding of the polypeptide. Using a variety of biophysical approaches, we characterise these linked equilibria for the isolated DNA-binding domains (DBD) of the receptors for estrogen and glucocorticoid. Mass spectrometry and equilibrium denaturation indicate that, near neutral pH, approximately four of the eight co-ordinating thiol groups release protons with zinc uptake, in agreement with the expected pK(a) change for the -SH group in the presence of the metal. Mass spectrometry reveals that the protein charge distribution changes with the uptake of zinc and that metal binding is co-operative. The co-operativity is consistent with observations from equilibrium denaturation, which indicate that the folding event is a two-state process. A crucial residue that stabilises the equilibrium structure of the DBD fold itself is a cysteine residue situated in the hydrophobic core of all known nuclear hormone receptors (but not involved in metal binding): it appears to be conserved absolutely for its unique combination of size and hydrophobicity. Stabilisation of the DBDs could be achieved by truncating the flexible, basic termini, suggesting that like-charge clusters may have deleterious effects on protein folds. While the metal-free apo protein and the chemically denatured state have little defined secondary structure, these states were expanded only partially in comparison with the native structure, according to data from small-angle X-ray scattering. The comparatively compact shapes of the denatured and apo forms may explain, in part, the marginal stability of the native fold.

Amino acids in SRS1 and SRS6 are critical for furanocoumarin metabolism by CYP6B1v1, a cytochrome P450 monooxygenase

CYP6B1v1 is the principal cytochrome P450 monooxygenase (P450) that detoxifies dietary furanocoumarins in the guts of Papilio polyxenes, the black swallowtail caterpillar. Sequence alignments and structure comparisons of CYP6B1v1 with the mouse CYP2A5 and bacterial CYP102 proteins, which are also capable of metabolizing the linear furanocoumarin xanthotoxin (8-methoxypsoralen), suggested that Phe116, His117, Val368 and Phe484 might be active site residues. In a homology model developed for CYP6B1v1, the side chains of Phe116 and His117 located in the B'-C loop of SRS1 are predicted to be positioned above the haem plane, while the side chain of Phe484 located in SRS6 is predicted near the entrance of the catalytic pocket. Site-directed mutagenesis of residues Phe116, His117 and Phe484 indicated that these residues represent several of those that determine this protein's stability and substrate specificity. Whereas all aromatic mutants of Phe116 and Phe484 generated CO-difference spectra with maxima at 450 nm indicative of correctly configured monooxygenases, aromatic mutants of Phe116 exhibited reduced reactivities toward some furanocoumarins and aromatic mutants of Phe484 eliminated all reactivities toward furanocoumarins. All single and double aliphatic mutants of Phe116, His117 and Phe484 and aromatic mutants of His117 generated carbon monoxide (CO) difference spectra with maxima at 420 nm (P420) indicative of incorrectly configured monooxygenases. These studies define residues Phe116, His117 and Phe484 as determinants of this insect P450's catalytic site integrity and residues Phe116 and Phe484 as determinants of its substrate specificity. Conservation of Phe116 and His117 in an array of lepidopteran CYP6B proteins implies that these amino acids serve a similar function in other monooxygenases of the insect CYP6B subfamily.

The hidden thermodynamics of a zinc finger

The Zn finger provides a model for studies of protein structure and stability. Its core contains a conserved phenylalanine residue adjoining three architectural elements: a beta-hairpin, an alpha-helix and a tetrahedral Zn(2+)-binding site. Here, we demonstrate that the consensus Phe is not required for high-affinity Zn(2+) binding but contributes to the specification of a precise DNA-binding surface. Substitution of Phe by leucine in a ZFY peptide permits Zn(2+)-dependent folding. Although a native-like structure is retained, structural fluctuations lead to attenuation of selected nuclear Overhauser enhancements and accelerated amide proton exchange. Surprisingly, wild-type Zn affinity is maintained by entropy-enthalpy compensation (EEC): a hidden entropy penalty (TDeltaDeltaS 7kcal/mol) is balanced by enhanced enthalpy of association (DeltaDeltaH -7kcal/mol) at 25 degrees C. Because the variant is less well ordered than the Phe-anchored domain, the net change in entropy is opposite to the apparent change in configurational entropy. By analogy to the thermodynamics of organometallic complexation, we propose that EEC arises from differences in solvent reorganization. Exclusion of Leu among biological sequences suggests an evolutionary constraint on the dynamics of a Zn finger.

Compound effects of point mutations causing campomelic dysplasia/autosomal sex reversal upon SOX9 structure, nuclear transport, DNA binding, and transcriptional activation

Human mutations in the transcription factor SOX9 cause campomelic dysplasia/autosomal sex reversal. Here we identify and characterize two novel heterozygous mutations, F154L and A158T, that substitute conserved "hydrophobic core" amino acids of the high mobility group domain at positions thought to stabilize SOX9 conformation. Circular dichroism studies indicated that both mutations disrupt alpha-helicity within their high mobility group domain, whereas tertiary structure is essentially maintained as judged by fluorescence spectroscopy. In cultured cells, strictly nuclear localization was observed for wild type SOX9 and the F154L mutant; however, the A158T mutant showed a 2-fold reduction in nuclear import efficiency. Importin-beta was demonstrated to be the nuclear transport receptor recognized by SOX9, with both mutant proteins binding importin-beta with wild type affinity. Whereas DNA bending was unaffected, DNA binding was drastically reduced in both mutants (to 5% of wild type activity in F154L, 17% in A158T). Despite this large effect, transcriptional activation in cultured cells was only reduced to 26% in F154L and 62% in A158T of wild type activity, suggesting that a small loss of SOX9 transactivation activity could be sufficient to disrupt proper regulation of target genes during bone and testis formation. Thus, clinically relevant mutations of SOX9 affect protein structure leading to compound effects of reduced nuclear import and reduced DNA binding, the net effect being loss of transcriptional activation.

Aromatic Ring Currents at a Protein Surface: Use of (1)H-NMR Chemical Shifts to Refine the Structure of a Naked β Sheet

Abstract The naked β sheet, a newly recognized motif of protein structure, exhibits ordered surfaces in the absence of a conventional hydrophobic core. A model is provided by an archaeal Zn ribbon homologous to eukaryotic RNA polymerase II subunit 9 (RPB9). This subunit, which regulates transcriptional start-site selection and downstream pausing, contains Zn(2+)-binding motifs similar to those of general transcription factors TFIIB and TFIIS. Interestingly, distance-geometry yields two models of the archaeal Zn ribbon differing in the orientation of a conserved tyrosine side chain on the well-ordered surface of the naked β-sheet. The models are equally consistent with conventional restraints and otherwise contain indistinguishable structural features, including a tetrahedral Cys(4) Zn(2+)-binding sites, four antiparallel β-strands, and disordered loop. Due to the change in tyrosine orientation and correlated changes in the configuration of neighboring side chains, the two models predict inequivalent patterns of aromatic ring-current shifts. The observed secondary shifts of adjoining resonances are shown to be consistent with one model but not the other. In the consistent model the surface of the β-sheet contains successive aromatic edge-to-face contacts in accord with semi-classical and ab initio potentials. We speculate that the aromatic-rich surface of the hyperthermophilic RPB9 domain contributes its thermodynamic stability and provides a nucleic-acid-binding site in the eukaryotic and archaeal transcriptional machinery. The present study demonstrates how the reduced dimensionality of a surface can lead to ambiguities in the interpretation of nuclear Overhauser enhancements. The results illustrate the utility of chemical shifts at such a surface in the cross-validation of a high-resolution solution structure.

High-resolution structure of an archaeal zinc ribbon defines a general architectural motif in eukaryotic RNA polymerases

BACKGROUND

Transcriptional initiation and elongation provide control points in gene expression. Eukaryotic RNA polymerase II subunit 9 (RPB9) regulates start-site selection and elongational arrest. RPB9 contains Cys4 Zn(2+)-binding motifs which are conserved in archaea and homologous to those of the general transcription factors TFIIB and TFIIS.

RESULTS

The structure of an RPB9 domain from the hyperthermophilic archaeon Thermococcus celer was determined at high resolution by NMR spectroscopy. The structure consists of an apical tetrahedral Zn(2+)-binding site, central beta sheet and disordered loop. Although the structure lacks a globular hydrophobic core, the two surfaces of the beta sheet each contain well ordered aromatic rings engaged in serial edge-to-face interactions. Basic sidechains are clustered near the Zn(2+)-binding site. The disordered loop contains sidechains conserved in TFIIS, including acidic residues essential for the stimulation of transcriptional elongation.

CONCLUSIONS

The planar architecture of the RPB9 zinc ribbon-distinct from that of a conventional globular domain-can accommodate significant differences in the alignment of polar, non-polar and charged sidechains. Such divergence is associated with local and non-local changes in structure. The RPB9 structure is distinguished by a fourth beta strand (extending the central beta sheet) in a well ordered N-terminal segment and also differs from TFIIS (but not TFIIB) in the orientation of its apical Zn(2+)-binding site. Cys4 Zn(2+)-binding sites with distinct patterns of polar, non-polar and charged residues are conserved among unrelated RNAP subunits and predicted to form variant zinc ribbons.

Lessons from zinc-binding peptides

Zinc-finger domains are small metal-binding modules that are found in a wide range of gene regulatory proteins. Peptides corresponding to these domains have provided valuable model systems for examining a number of biophysical parameters entirely unrelated to their nucleic acid binding properties. These include the chemical basis for metal-ion affinity and selectivity, thermodynamic properties related to hydrophobic packing and beta-sheet propensities, and constraints on the generation of ligand-binding and potential catalytic sites. These studies have laid the foundation for applications such as the generation of optically detected zinc probes and the design of metal-binding peptides and proteins with desired spectroscopic and chemical properties.

Probing weakly polar interactions in cytochrome c

Theoretical, statistical, and model studies suggest that proteins are stabilized by weakly polar attractions between sulfur atoms and properly oriented aromatic rings. The two sulfur-containing amino acids, methionine and cysteine, occur frequently among functional alleles in random mutant libraries of Saccharomyces cerevisiae iso-1-cytochrome c genes at positions that form a weakly polar aromatic-aromatic interaction, the wild-type protein. To determine if a weakly polar sulfur-aromatic interaction replaced the aromatic-aromatic interaction, the structure and stability of two variants were examined. Phenylalanine 10, which interacts with tyrosine 97, was replaced by methionine and cysteine. The cysteine was modified to form the methionine and cysteine analog, S-methyl cysteine (CysSMe). Proton NMR studies indicate that changing Phe 10 to Met or CysSMe affects only local structure and that the structures of sulfur-containing variants are nearly identical. Analysis of chemical shifts and nuclear Overhauser effect data indicates that both sulfur-containing side chains are in position to form a weakly polar interaction with Tyr 97. The F10M and F10CSMe variants are 2-3 kcal mol-1 less stable than iso-1-cytochrome c at 300 K. Comparison of the stabilities of the F10M and F10CSMe variants allows evaluation of the potential weakly polar interaction between the additional sulfur atom of F10CSMe and the aromatic moiety of Tyr 97. The F10CSMe;C102T variant is 0.7 +/- 0.3 kcal mol-1 more stable than the F10M;C102T protein. The increased stability is explained by the difference in hydrophobicity of the sulfur-containing side chains. We conclude that any weakly polar interaction between the additional sulfur and the aromatic ring is too weak to detect or is masked by destabilizing contributions to the free energy of denaturation.

A 1H-NMR study of the transcription factor 1 from Bacillus subtilis phage SPO1 by selective 2H-labeling. Complete assignment and structural analysis of the aromatic resonances for a 22-kDa homodimer

1H-NMR experiments have been performed on transcription factor 1 (TF1) encoded by Bacillus subtilis phage SPO1. To study this 22-kDa homodimeric DNA-binding protein, a selective 2H-labeling strategy has been employed. Complete sequence-specific assignments of all the resonances from the five aromatic residues were determined by a modified standard sequential-assignment procedure. The reduced contribution of spin diffusion upon the long-mixing-time nuclear-Overhauser-enhancement spectroscopy for the selectively 2H-labeled variants, as opposed to the fully 1H-containing protein, has allowed for the identification of the spin systems and of the long-range dipolar contacts between Phe28 and Phe47 protons in the protein core and between Phe61 and Phe97 protons. The latter suggests an interaction between the proposed beta-ribbon DNA-binding arm and the carboxy terminus of the paired monomer. A previously proposed TF1 structural model [Geiduschek, E. P., Schneider, G. J. & Sayre, M. H. (1990) J. Struct. Biol. 104, 84-90)] has been modified using constrained-energy-minimization calculations incorporating the experimentally determined set of aromatic-to-aromatic contacts. This new model has been analyzed with regard to the relative mobility and the relative solvent accessibility of the aromatic residues which have been measured by the nonselective T1 relaxation times of the aromatic resonances for the fully 1H-containing protein and the relaxation time enhancements upon selective 2H-labeling, respectively.