Synthesis and DNA-binding ability of Spl protein zinc finger domain and its peptidomimetics

The second zinc finger fragment of Sp1 (Spl-ZF2), its mutant (Spl-ZF2/HT. E20→H, R23→T), and two mimic analogues (ZF20 and ZF15) were synthesized by stepwise solid phase technique. The CD spectra and UV-visible spectrum with CoCl2 indicated that the formation of zinc finger structure was affected not only by the hy-drophobic amino acids but also by the change of the distance between Cys and His. Gel-retardation electrophoresis as-says indicated that the Grlu and Arg residues are very important for recognition. A single zinc finger like Spl-ZF2 isable to bind DNA sequence specifically.     

Functional analyses of the Dof domain, a zinc finger DNA-binding domain, in a pumpkin DNA-binding protein AOBP

AOBP, a DNA-binding protein in pumpkin, contains a Dof domain that is composed of 52 amino acid residues and is highly conserved in several DNA-binding proteins of higher plants. The Dof domain has a significant resemblance to Cys2/Cys2 zinc finger DNA-binding domains of steroid hormone receptors and GATA1, but has a longer putative loop where an extra Cys residue is conserved. We show that the Dof domain in AOBP functions as a zinc finger DNA-binding domain and suggest that the Cys residue uniquely conserved in the putative loop might negatively regulate the binding to DNA.     

Identification of amino acids important for target recognition by the DNA:m5C methyltransferase M.NgoPII by alanine-scanning mutagenesis of residues at the protein-DNA interface

DNA:m(5)C MTases comprise a catalytic domain with conserved residues of the active site and a strongly diverged TRD with variable residues involved in DNA recognition and binding. To date, crystal structures of 2 DNA:m(5)C MTases complexed with the substrate DNA have been obtained; however, for none of these enzymes has the importance of the whole set of DNA-binding residues been comprehensively studied. We built a comparative model of M.NgoPII, a close homologue and isomethylomer of M.HaeIII, and systematically analyzed the effect of alanine substitutions for the complete set of amino acid residues from its TRD predicted to be important for DNA binding and target recognition. Our data demonstrate that only 1 Arg residue is indispensable for the MTase activity in vivo and in vitro, and that mutations of only a few other residues cause significant reduction of the activity in vitro, with little effect on the activity in vivo. The identification of dispensable protein-DNA contacts in the wild-type MTase will serve as a platform for exhaustive combinatorial mutagenesis aimed at the design of new contacts, and thus construction of enzyme variants that retain the activity but exhibit potentially new substrate preferences.     

Hydration of protein-protein interfaces

We present an analysis of the water molecules immobilized at the protein-protein interfaces of 115 homodimeric proteins and 46 protein-protein complexes, and compare them with 173 large crystal packing interfaces representing nonspecific interactions. With an average of 15 waters per 1000 A2 of interface area, the crystal packing interfaces are more hydrated than the specific interfaces of homodimers and complexes, which have 10-11 waters per 1000 A2, reflecting the more hydrophilic composition of crystal packing interfaces. Very different patterns of hydration are observed: Water molecules may form a ring around interfaces that remain "dry," or they may permeate "wet" interfaces. A majority of the specific interfaces are dry and most of the crystal packing interfaces are wet, but counterexamples exist in both categories. Water molecules at interfaces form hydrogen bonds with protein groups, with a preference for the main-chain carbonyl and the charged side-chains of Glu, Asp, and Arg. These interactions are essentially the same in specific and nonspecific interfaces, and very similar to those observed elsewhere on the protein surface. Water-mediated polar interactions are as abundant at the interfaces as direct protein-protein hydrogen bonds, and they may contribute to the stability of the assembly.     

Characterization of HMBP-2, a DNA-binding protein that binds to HIV-1 LTR when only one of the three Sp1 sites is methylated

HIV-1 methylation binding protein-1 (HMBP-1) (formerly called HMBP) is a protein found in human cell nuclei that binds with enhanced affinity to a fragment of the HIV-1 long terminal repeat sequence (LTR) containing three Sp1 sites when all three sites are methylated. HMBP-2 is another protein present in the nuclei of human T helper lymphocytes and HeLa cells that binds to the HIV-1 LTR. HMBP-2 binds preferentially to the same region of the HIV-1 LTR as does HMBP-1, but HMBP-2 binds best when only one of the three Sp1 sites is methylated. HMBP-2 can be separated from HMBP-1 chromatographically, and dimethyl sulfate (DMS) methylation interference analysis indicates that their binding sites are not identical. HMBP-2 represents a novel protein factor capable of binding to a partially methylated region of the HIV-1 LTR.     

Structural determinants of binding and specificity in transforming growth factor-receptor interactions.

Transforming growth factor (TGF-beta) protein families are cytokines that occur as a large number of homologous proteins. Three major subgroups of these proteins with marked specificities for their receptors have been found-TGF-beta, activin/inhibin, and bone morphogenic protein. Although structural information is available for some members of the TGF-beta family of ligands and receptors, very little is known about the way these growth factors interact with the extracellular domains of their cell surface receptors, especially the type II receptor. In addition, the elements that are the determinants of binding and specificity of the ligands are poorly understood. The structure of the extracellular domain of the receptor is a three-finger fold similar to some toxin structures. Amino acid exchanges between multiply aligned homologous sequences of type II receptors point to a residue at the surface, specifically finger 1, as the determinant of ligand specificity and complex formation. The "knuckle" epitope of ligands was predicted to be the surface that interacts with the type II receptor. The residues on strands beta2, beta3, beta7, beta8 and the loop region joining beta2 and beta3 and joining beta7 and beta8 of the ligands were identified as determinants of binding and specificity. These results are supported by studies on the docking of the type II receptor to the ligand dimer-type I receptor complex.     

On the nature of cavities on protein surfaces: application to the identification of drug-binding sites

In this article we introduce a new method for the identification and the accurate characterization of protein surface cavities. The method is encoded in the program SCREEN (Surface Cavity REcognition and EvaluatioN). As a first test of the utility of our approach we used SCREEN to locate and analyze the surface cavities of a nonredundant set of 99 proteins cocrystallized with drugs. We find that this set of proteins has on average about 14 distinct cavities per protein. In all cases, a drug is bound at one (and sometimes more than one) of these cavities. Using cavity size alone as a criterion for predicting drug-binding sites yields a high balanced error rate of 15.7%, with only 71.7% coverage. Here we characterize each surface cavity by computing a comprehensive set of 408 physicochemical, structural, and geometric attributes. By applying modern machine learning techniques (Random Forests) we were able to develop a classifier that can identify drug-binding cavities with a balanced error rate of 7.2% and coverage of 88.9%. Only 18 of the 408 cavity attributes had a statistically significant role in the prediction. Of these 18 important attributes, almost all involved size and shape rather than physicochemical properties of the surface cavity. The implications of these results are discussed. A SCREEN Web server is available at http://interface.bioc.columbia.edu/screen.     

Single-molecule analysis of chromatin: changing the view of genomes one molecule at a time

Wrapping DNA into chromatin provides a wealth of regulatory mechanisms that ensure normal growth and development in eukaryotes. Our understanding of chromatin structure, including nucleosomes and non-histone protein-DNA interactions, has benefited immensely from nuclease and chemical digestion techniques. DNA-bound proteins, such as histones or site-specific factors, protect DNA against nuclease cleavage and generate large nucleosomal or small regulatory factor footprints. Chromatin subject to distinct modes of regulation often coincides with sites of nuclease hypersensitivity or nucleosome positioning. An inherent limitation of cleavage-based analyses has been the inability to reliably analyze regions of interest when levels of digestion depart from single-hit kinetics. Moreover, cleavage-based techniques provide views that are averaged over all the molecules in a sample population. Therefore, in cases of occupancy of multiple regulatory elements by factors, one cannot define whether the factors are bound to the same or different molecules in the population. The recent development of DNA methyltransferase-based, single-molecule MAP-IT technology overcomes limitations of ensemble approaches and has opened numerous new avenues in chromatin research. Here, we review the strengths, limitations, applications and future prospects of MAP-IT ranging from structural issues to mechanistic questions in eukaryotic chromatin regulation.     

Distance-constraint approach to protein folding. I. Statistical analysis of protein conformations in terms of distances between residues

A statistical analysis of protein conformations in terms of the distance between residues, represented by their C atoms, is presented. We consider four factors that contribute to the determination of the distanced _i,i+k between a given pair ofith and(i+k)th residues in the native conformation of a globular protein: (1) the distancek along the chain, (2) the size of the protein, (3) the conformational states of theith to(i+k)th residues, and (4) the amino acid types of the and(i+k)th residues. In order to account for the dependence on the distancek along the chain, the statistics are taken for three ranges, viz., short, medium, and long ranges (k8; 9k20; andk21; respectively). In the statistics of short-range distances, a mean distanceD _k and its standard deviationS _k are calculated for each value ofk, with and without taking into account the conformational states of all residues fromi toi+k (factors 1 and 3). As an Appendix, the relations for converting from the distances between residues into other conformational parameters are discussed. In the statistics of long-range distances, a reduced distanced* _ij (the actual distance divided by the radius of gyration) is used to scale the data so that they become independent of protein size, and then a mean reduced distanceD _l (a, a) and its standard deviation _l (a, a) are calculated for each amino acid pair (a, a) (factors 2 and 4). The effect of the neighboring residues along the chain on the value of the distanced* _ij is explored by a linear regression analysis between the actual reduced distanced* _ij and the mean value over theD _l for all possible pairs of residues in the two segments of the (i–2)th to the (i+2)th and the (j–2)th to the (j+2)th residues. The effect is assessed in terms of the tangentA _l (a, a) of the calculated regression line for each amino acid pair (a, a). In the statistics of medium-range distances, only factors 1 and 4 are considered, to simplify the analysis. The scaled distanced _i,i+k =(d _i,i+k -D _k )/S _k is used to eliminate the dependence onk, the distance along the chain. The propertiesD _m (a, a), _m (a, a) andA _m (a, a) corresponding toD _l (a, a), _l (a, a), andA _l (a, a), and also calculated for each amino acid pair (a, a). The results are interpreted as follows: the smaller values ofD _l (a, a) andD _m (a, a) indicate a preference of the pair (a, a) for a contact (e.g., pairs between hydrophobic amino acids, and pairs of Cys with aromatic amino acids), and the larger values of these quantities indicate a preference for distant mutual location (e.g., pairs between strong hydrophilic amino acids); the smaller values of _l (a, a) and _m (a, a) indicate a strong preference for either contact or noncontact (e.g., pairs between hydrophobic amino acids, and pairs between strong hydrophobic and hydrophilic amino acids, respectively), and the larger values of these quantities indicate the ambivalent/neutral nature of the preference for contact and noncontact (e.g., pairs containing Ser or Thr); the smaller values ofA _l (a, a) andA _m (a, a) indicate that the distance of an (a, a) pair is determined independently of the amino acid character of the neighboring residues along the chain (e.g., some pairs of Cys or Met with other amino acids) and the larger values of these quantities indicare that such amino acid character contributes strongly to the determination of the distance (e.g., pairs containing Ser or Thr, and pairs between amino acids with small side chains). The difference between the statistics for the long- and medium-range distances is also discussed; the former reflect the difference between the hydrophobic and hydrophilic character of the residues, but the latter cannot be easily interpretable only in terms of hydrophobicity and hydrophilicity. The data analyzed here are used in the optimization of an object function to compute protein conformation in a subsequent paper.     

Dynamic states of the DNA repair enzyme AlkB regulate product release

The 2-oxoglutarate (2OG)- and Fe(2+)-dependent dioxygenase AlkB couples the demethylation of modified DNA to the decarboxylation of 2OG. Extensive crystallographic analyses have shown no evidence of significant structural differences between complexes binding either 2OG or succinate. By using nuclear magnetic resonance spectroscopy, we have shown that the AlkB-succinate and AlkB-2OG complexes have significantly different dynamic properties in solution. 2OG makes the necessary contacts between the metal site and the large beta-sheet to maintain a fully folded conformation. Oxidative decarboxylation of 2OG to succinate leads to weakening of a main contact with the large beta-sheet, resulting in an enhanced dynamic state. These conformational fluctuations allow for the replacement of succinate in the central core of the protein and probably contribute to the effective release of unmethylated DNA. We also propose that the inherent dynamics of the co-product complex and the subsequent increased molecular ordering of the co-substrate complex have a role in DNA damage recognition.     

TRIADs: a new class of proteins with a novel cysteine-rich signature.

Triad1 was recently identified as a nuclear RING finger protein, which is up-regulated during retinoic acid induced granulocytic differentiation of acute leukemia cells. Here we show that a cysteine-rich domain (C6HC), present in Triad1, is conserved in at least 24 proteins encoded by various eukaryotes. The C6HC consensus pattern C-x(4)-C-x(14-30)-C-x(1-4)-C-x(4)-C-x(2)-C-x(4)-H-x(4)-C defines this structure as the fourth family member of the zinc-binding RING, LIM, and LAP/PHD fingers. Strikingly, in 22 of 24 proteins the C6HC domain is flanked by two RING finger structures. We have termed the novel C6HC motif DRIL (double RING finger linked). The strong conservation of the larger tripartite TRIAD (two RING fingers and DRIL) structure indicates that the three subdomains are functionally linked and identifies a novel class of proteins.     

Interaction interfaces of protein domains are not topologically equivalent across families within superfamilies: Implications for metabolic and signaling pathways

Using a data set of aligned protein domain superfamilies of known three-dimensional structure, we compared the location of interdomain interfaces on the tertiary folds between members of distantly related protein domain superfamilies. The data set analyzed is comprised of interdomain interfaces, with domains occurring within a polypeptide chain and those between two polypeptide chains. We observe that, in general, the interfaces between protein domains are formed entirely in different locations on the tertiary folds in such pairs. This variation in the location of interface happens in protein domains involved in a wide range of functions, such as enzymes, adapters, and domains that bind protein ligands, or cofactors. While basic biochemical functionality is preserved at the domain superfamily level, the effect of biochemical function on protein assemblies is different in these protein domains related by superfamily. The divergence between proteins, in most cases, is coupled with domain recruitment, with different modes of interaction with the recruited domain. This is in complete contrast to the observation that in closely related homologous protein domains, almost always the interaction interfaces are topologically equivalent. In a small subset of interacting domains within proteins related by remote homology, we observe that the relative positioning of domains with respect to one another is preserved. Based on the analysis of multidomain proteins of known or unknown structure, we suggest that variation in protein-protein interactions in members within a superfamily could serve as diverging points in otherwise parallel metabolic or signaling pathways. We discuss a few representative cases of diverging pathways involving domains in a superfamily.     

Structure and dynamics of the DNA-binding protein HU of B. stearothermophilus investigated by Raman and ultraviolet-resonance Raman spectroscopy

The histone-like protein HU of Bacillus stearothermophilus (HUBst) is a 90-residue homodimer that binds nonspecifically to B DNA. Although the structure of the HUBst:DNA complex is not known, the proposed DNA-binding surface consists of extended arms that project from an alpha-helical platform. Here, we report Raman and ultraviolet-resonance Raman (UVRR) spectra diagnostic of subunit secondary structures and indicative of key side-chains lining the proposed DNA-binding surface. Raman conformation markers show that the DNA-binding arms of the dimer contain beta-stranded structure in excess (eight +/- two residues per subunit) of that reported previously. Important among side-chain markers are Met (701 cm(-1)), Ala (908 cm(-1)), Arg (1082 cm(-1)), and Pro (1457 cm(-1)). The Ala marker undergoes a substantial shift (908 --> 893 cm(-1)) on deuteration of alanyl peptide sites, indicating a coupled side-chain/main-chain mode of diagnostic value in the identification of exchange-protected alanines. A large subset of alanines (67%) in the alpha-helical core exhibits robust resistance to exchange. A quantitative study of NH --> ND exchange exploiting newly identified amide II' markers of helical (1440 cm(-1)) and nonhelical (1472 cm(-1)) conformations of HUBst indicates unexpected flexibility at the dimer interface, which is manifested in rapid exchange of 80% of peptide sites. The results establish a basis for subsequent Raman and UVRR investigations of HUBst:DNA complexes and provide a framework for applications to other DNA-binding architectural proteins.     

Molecular dynamics simulations of the cytochrome c3-rubredoxin complex from Desulfovibrio vulgaris.

Molecular dynamics simulations have been carried out on the complex formed between the tetraheme cytochrome c3 and the iron protein rubredoxin from the sulfate-reducing bacterium Desulfovibrio vulgaris. These simulations were performed both with explicit solvent water molecules included, and without solvent molecules using a distance-dependent dielectric constant to approximate the screening effects of solvent. The results of both simulations are strikingly different, indicating that the representation of environmental effects is important in such simulations. For example, a striking adaptation of the two proteins seen in the nonsolvated simulation is not seen when explicit solvent water is included; in fact, the complex appears to become weaker in the solvated simulation. Nonetheless, the iron-iron distance decreases more significantly in the solvated simulation than in the nonsolvated simulation. It was found that in both cases molecular dynamics optimized the structures further than energy minimization alone.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.