A structure-based method for identifying DNA-binding proteins and their sites of DNA-interaction

A classification model of a DNA-binding protein chain was created based on identification of alpha helices within the chain likely to bind to DNA. Using the model, all chains in the Protein Data Bank were classified. For many of the chains classified with high confidence, previous documentation for DNA-binding was found, yet no sequence homology to the structures used to train the model was detected. The result indicates that the chain model can be used to supplement sequence based methods for annotating the function of DNA-binding. Four new candidates for DNA-binding were found, including two structures solved through structural genomics efforts. For each of the candidate structures, possible sites of DNA-binding are indicated by listing the residue ranges of alpha helices likely to interact with DNA.     

Osteogenesis imperfecta murine: interaction between type I collagen homotrimers

Types I, II, and III collagens are believed to have evolved from the same homotrimer ancestor and they have substantial sequence homology, but type I molecules are alpha1(I)(2)alpha2(I) heterotrimers, unlike homotrimeric types II and III. It is believed that the alpha2(I) chain first appeared in lower vertebrates and that it plays a particularly important role in bone formation. For instance, spontaneous mutations resulting in non- functional alpha2 chains and formation of type I homotrimers cause severe bone pathology (osteogenesis imperfecta) in humans and in animals. However, the exact role of the alpha2 chain is not known. Here, we report measurements of intermolecular forces between collagen helices in native and reconstituted fibers composed of type I homotrimers, heterotrimers and their mix. For comparison, we report forces between type II homotrimers in reconstituted fibers. In agreement with previous studies, we find that the absence of the alpha2 chain reduces temperature-favored attraction between collagen helices, either because of the difference in amino acid sequence of the alpha1 and alpha2 chains or because of more extensive post-translational modification of homotrimers. We find that forces between helices in fibers from type I (as well as type II) homotrimers are not sensitive to pH between pH 6 and 7.5, in contrast to type I heterotrimers. Apparently, the effect of pH is related to extra histidine residues present on alpha2 chains but not on alpha1 chains. Finally, our measurements indicate that the alpha2 chain is responsible for binding some soluble compound(s), possibly glycosaminoglycans, whose displacement results, e.g., in the loss of tendon crystallinity. The ability of the alpha2 chain to bind non-collagen matrix components may be particularly important for bone matrix formation and mineralization.     

Structural and sequence characteristics of long alpha helices in globular proteins.

Elucidation of the detailed structural features and sequence requirements for alpha helices of various lengths could be very important in understanding secondary structure formation in proteins and, hence, in the protein folding mechanism. An algorithm to characterize the geometry of an alpha helix from its C(alpha) coordinates has been developed and used to analyze the structures of long alpha helices (number of residues > or = 25) found in globular proteins, the crystal structure coordinates of which are available from the Brookhaven Protein Data Bank. All long alpha helices can be unambiguously characterized as belonging to one of three classes: linear, curved, or kinked, with a majority being curved. Analysis of the sequences of these helices reveals that the long alpha helices have unique sequence characteristics that distinguish them from the short alpha helices in globular proteins. The distribution and statistical propensities of individual amino acids to occur in long alpha helices are different from those found in short alpha helices, with amino acids having longer side chains and/or having a greater number of functional groups occurring more frequently in these helices. The sequences of the long alpha helices can be correlated with their gross structural features, i.e., whether they are curved, linear, or kinked, and in case of the curved helices, with their curvature.     

Mapping backbone dynamics in solution with site-directed spin labeling: GCN4-58 bZip free and bound to DNA

In site-directed spin labeling, a nitroxide-containing side chain is introduced at selected sites in a protein. The EPR spectrum of the labeled protein encodes information about the motion of the nitroxide on the nanosecond time scale, which has contributions from the rotary diffusion of the protein, from internal motions in the side chain, and from backbone fluctuations. In the simplest model for the motion of noninteracting (surface) side chains, the contribution from the internal motion is sequence independent, as is that from protein rotary diffusion. Hence, differences in backbone motions should be revealed by comparing the sequence-dependent motions of nitroxides at structurally homologous sites. To examine this model, nitroxide side chains were introduced, one at a time, along the GCN4-58 bZip sequence, for which NMR (15)N relaxation experiments have identified a striking gradient of backbone mobility along the DNA-binding region [Bracken et al. (1999) J. Mol. Biol. 285, 2133]. Spectral simulation techniques and a simple line width measure were used to extract dynamical parameters from the EPR spectra, and the results reveal a mobility gradient similar to that observed in NMR relaxation, indicating that side chain motions mirror backbone motions. In addition, the sequence-dependent side chain dynamics were analyzed in the DNA/protein complex, which has not been previously investigated by NMR relaxation methods. As anticipated, the backbone motions are damped in the DNA-bound state, although a gradient of motion persists with residues at the DNA-binding site being the most highly ordered, similar to those of helices on globular proteins.     

Seven-helix bundles: molecular modeling via restrained molecular dynamics.

Simulated annealing via restrained molecular dynamics (SA/MD) has been used to model compact bundles of seven approximately (anti)parallel alpha-helices. Seven such helix bundles occur, e.g., in bacteriorhodopsin, in rhodopsin, and in the channel-forming N-terminal domain of Bacillus thuringiensis delta-endotoxin. Two classes of model are considered: (a) those consisting of seven Ala20 peptide chains; and (b) those containing a single polypeptide chain, made up of seven Ala20 helices linked by GlyN interhelix loops (where N = 5 or 10). Three different starting C alpha templates for SA/MD are used, in which the seven helices are arranged (a) on a left-handed circular template, (b) on a bacteriorhodopsin-like template, or (c) on a zig-zag template. The ensembles of models generated by SA/MD are analyzed in terms of their geometry and energetics, and the most stable structures from each ensemble are examined in greater detail. Structures resembling bacteriorhodopsin and structures resembling delta-endotoxin are both represented among the most stable structures. delta-Endotoxin-like structures arise from both circular and bacteriorhodopsin-like C alpha templates. A third helix-packing mode occurs several times among the stable structures, regardless of the C alpha template and of the presence or absence of interhelix loops. It is characterized by a "4 + 1" core, in which four helices form a distorted left-handed supercoil around a central, buried helix. The remaining two helices pack onto the outside of the core. This packing mode is comparable with that proposed for rhodopsin on the basis of two-dimensional electron crystallographic and sequence analysis studies.     

A hydrogen bonded chain in bacteriorhodopsin by computer modelling approach

The seven alpha-helical segments of Bacteriorhodopsin (BR) passing through the membrane are investigated for a continuous Hydrogen Bonded Chain (HBC). The study is carried out by computer modelling approach. It is assumed that the seven helices are placed as (AGFEDCB), which has been accepted as the best model by several groups. Helices A, D, E and G are considered to be present in right handed alpha-helical conformation. The inter-orientation of these helices are represented by Eulerian angles alpha, beta and gamma. For the helices B, C and F which contain Proline in the middle, several conformational possibilities were considered. In these cases apart from the Eulerian angles alpha, beta and gamma, the dihedral angles phi p-1 and psi p-1 of the residues that are succeeded by Proline residue in the helical regions were also used in fixing the position of the helices with respect to each other. All these parameters were varied to fit with the top, middle and bottom distances reported by electron diffraction studies. Good fit was obtained for all right handed alpha-helical conformations and also for helices B, C and F with a left handed turn at the residue preceeding proline. Hence two structures were analysed for continuous HBC. Structure I which contained all the seven helices in right handed alpha-helical conformation and Structure II, which had the helices A, D, E and G in right handed conformation and the helices B, C and F in right handed alpha-helical conformation with a left handed turn at the residue preceeding proline. All possible staggered conformations were considered for the side chains and the inter atomic distances were analysed for Hydrogen bonds. It was possible to obtain a continuous chain in both the structures which includes most of the residues found to be important by the experiments. However Lys-216 has to be considered in two different conformations to connect the cytoplasmic side with the extra cellular side. The overall height spanned by HBC is about 25A. The chains obtained by both the structures I and II are analysed in terms of the conformational parameters. It has also been possible to place the retinal in the region as predicted by the experiments. The Tryptophan residues which affect the spectral characteristics can be aligned on either side of the retinal.     

Comparison between TRF2 and TRF1 of their telomeric DNA-bound structures and DNA-binding activities

Mammalian telomeres consist of long tandem arrays of double-stranded telomeric TTAGGG repeats packaged by the telomeric DNA-binding proteins TRF1 and TRF2. Both contain a similar C-terminal Myb domain that mediates sequence-specific binding to telomeric DNA. In a DNA complex of TRF1, only the single Myb-like domain consisting of three helices can bind specifically to double-stranded telomeric DNA. TRF2 also binds to double-stranded telomeric DNA. Although the DNA binding mode of TRF2 is likely identical to that of TRF1, TRF2 plays an important role in the t-loop formation that protects the ends of telomeres. Here, to clarify the details of the double-stranded telomeric DNA-binding modes of TRF1 and TRF2, we determined the solution structure of the DNA-binding domain of human TRF2 bound to telomeric DNA; it consists of three helices, and like TRF1, the third helix recognizes TAGGG sequence in the major groove of DNA with the N-terminal arm locating in the minor groove. However, small but significant differences are observed; in contrast to the minor groove recognition of TRF1, in which an arginine residue recognizes the TT sequence, a lysine residue of TRF2 interacts with the TT part. We examined the telomeric DNA-binding activities of both DNA-binding domains of TRF1 and TRF2 and found that TRF1 binds more strongly than TRF2. Based on the structural differences of both domains, we created several mutants of the DNA-binding domain of TRF2 with stronger binding activities compared to the wild-type TRF2.     

Light chain determines the binding property of human anti-dsDNA IgG autoantibodies

We have previously prepared human anti-double-stranded (ds) DNA IgG Fab clones using phage-display technology. Nucleotide sequence analysis of genes of immunoglobulin (Ig) heavy and light chain variable regions in these Fab clones suggested that the DNA-binding activity of the clones depended on light chain usage. To confirm the role of the light chain in antibody binding to DNA, we constructed in the present study's new recombined Fab clones by heavy and light chain shuffling between the original anti-dsDNA Fab clones. Clones constructed by pairing Fdgamma fragments with the light chain from a high DNA-binding clone showed high DNA-binding activities, whereas other constructed clones using light chains from low DNA-binding clones showed low DNA-binding activities. Our results indicate that light chains in anti-dsDNA antibodies can determine the DNA-binding activity of the antibodies. Ig chain shuffling of phage-display antibodies may be useful for investigating the molecular mechanisms for antigen-antibody binding of human autoantibodies.     

Where is the weak linkage in the globin chain?

Hemoglobinopathies are important inherited disorders with high prevalence in many tropical countries. Prediction of protein nanostructure and function is a great challenge in proteomics and structural genomics. Identifying the point vulnerable to mutation is a new trend in research on disorders at the genomic and proteomic level. A bioinformatics analysis was performed to determine the positions that tend to correspond with peptide motifs in the amino acid sequence of alpha and beta globin chains. To identify the weak linkage in alpha globin and beta globin chains, a new bioinformatics tool, GlobPlot, was used. For the alpha globin chain, 22 positions were identified: the disorders were found at positions 3-8, 38-42, 46-51, and 75-79. For the beta globin chain, 46 positions were identified: the disorders were found at positions 61-146. The study showed that weak linkages in alpha globin and beta globin chains can be identified and provide good information for predicting possible new mutations that could lead to new hemoglobinopathies.     

A structural model of the acetylcholine receptor channel based on partition energy and helix packing calculations.

A structural model of the transmembrane portion of the acetylcholine receptor was developed from sequences of all its subunits by using transfer energy calculations to locate transmembrane alpha-helices and to calculate which helical side chains should be in contact with water inside the channel, with portions of other transmembrane helices, or with lipid hydrocarbon chains. "Knobs-into-holes" side chain packing calculations were used with other factors to stack the transmembrane alpha-helices together. In the model each subunit has the following structures in order along the sequence from the NH2 terminus: a large extracellular domain of undetermined structure, a short apolar alpha-helix that lies on the extracellular lipid surface of the membrane; three apolar transmembrane alpha-helices (I, II, and III), a cytoplasmic domain of undetermined structure, an amphipathic transmembrane alpha-helix (L) that forms the channel lining, a short extracellular alpha-helix, another apolar transmembrane alpha-helix (IV), and a small cytoplasmic domain formed by the COOH-terminal end of the chain. Three concentric layers form the pore. A bundle of five amphipathic L helices forms the channel lining. This bundle is surrounded by a bundle of 10 alternating II and III helices. Helices I and IV cover portions of the outer surface of the bundle formed by helices II and III. Positions of disulfide bridges are predicted and a mechanism for opening and closing conformational changes is proposed that requires tilting transmembrane helices and possibly a thiol-disulfide interchange reaction.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Family of G protein alpha chains: amphipathic analysis and predicted structure of functional domains

The G proteins transduce hormonal and other signals into regulation of enzymes such as adenylyl cyclase and retinal cGMP phosphodiesterase. Each G protein contains an alpha subunit that binds and hydrolyzes guanine nucleotides and interacts with beta gamma subunits and specific receptor and effector proteins. Amphipathic and secondary structure analysis of the primary sequences of five different alpha chains (bovine alpha s, alpha t1 and alpha t2, mouse alpha i, and rat alpha o) predicted the secondary structure of a composite alpha chain (alpha avg). The alpha chains contain four short regions of sequence homologous to regions in the GDP binding domain of bacterial elongation factor Tu (EF-Tu). Similarities between the predicted secondary structures of these regions in alpha avg and the known secondary structure of EF-Tu allowed us to construct a three-dimensional model of the GDP binding domain of alpha avg. Identification of the GDP binding domain of alpha avg defined three additional domains in the composite polypeptide. The first includes the amino terminal 41 residues of alpha avg, with a predicted amphipathic alpha helical structure; this domain may control binding of the alpha chains to the beta gamma complex. The second domain, containing predicted beta strands and alpha helices, several of which are strongly amphipathic, probably contains sequences responsible for interaction of alpha chains with effector enzymes. The predicted structure of the third domain, containing the carboxy terminal 100 amino acids, is predominantly beta sheet with an amphipathic alpha helix at the carboxy terminus. We propose that this domain is responsible for receptor binding.(ABSTRACT TRUNCATED AT 250 WORDS)     

Analysis of lipooligosaccharide biosynthesis in the Neisseriaceae

Neisserial lipooligosaccharide (LOS) contains three oligosaccharide chains, termed the alpha, beta, and gamma chains. We used Southern hybridization experiments on DNA isolated from various Neisseria spp. to determine if strains considered to be nonpathogenic possessed DNA sequences homologous with genes involved in the biosynthesis of these oligosaccharide chains. The presence or absence of specific genes was compared to the LOS profiles expressed by each strain, as characterized by their mobilities on sodium dodecyl sulfate-polyacrylamide gel electrophoresis gel and their reactivities with various LOS-specific monoclonal antibodies. A great deal of heterogeneity was seen with respect to the presence of genes encoding glycosyltransferases in Neisseria. All pathogenic species were found to possess DNA sequences homologous with the lgt gene cluster, a group of genes needed for the synthesis of the alpha chain. Some of these genes were also found to be present in strains considered to be nonpathogenic, such as Neisseria lactamica, N. subflava, and N. sicca. Some nonpathogenic Neisseria spp. were able to express high-molecular-mass LOS structures, even though they lacked the DNA sequences homologous with rfaF, a gene whose product must act before gonococcal and meningococcal LOS can be elongated. Using a PCR amplification strategy, in combination with DNA sequencing, we demonstrated that N. subflava 44 possessed lgtA, lgtB, and lgtE genes. The predicted amino acid sequence encoded by each of these genes suggested that they encoded functional proteins; however, structural analysis of LOS isolated from this strain indicated that the bulk of its LOS was not modified by these gene products. This suggests the existence of an additional regulatory mechanism that is responsible for the limited expression of these genes in this strain.     

Computer simulations of membrane protein folding: structure and dynamics

A lattice model of membrane proteins with a composite energy function is proposed to study their folding dynamics and native structures using Monte Carlo simulations. This model successfully predicts the seven helix bundle structure of sensory rhodopsin I by practicing a three-stage folding. Folding dynamics of a transmembrane segment into a helix is further investigated by varying the cooperativity in the formation of alpha helices for both random folding and assisted folding. The chain length dependence of the folding time of a hydrophobic segment to a helical state is studied for both free and anchored chains. An unusual length dependence in the folding time of anchored chains is observed.     

NMR solution structure of the archaebacterial chromosomal protein MC1 reveals a new protein fold

The three-dimensional structure of methanogen chromosomal protein 1 (MC1), a chromosomal protein extracted from the archaebacterium Methanosarcina sp. CHTI55, has been solved using (1)H NMR spectroscopy. The small basic protein MC1 contains 93 amino acids (24 basic residues against 12 acidic residues). The main elements of secondary structures are an alpha helix and five beta strands, arranged as two antiparallel beta sheets (a double one and a triple one) packed in an orthogonal manner forming a barrel. The protein displays a largely hydrophilic surface and a very compact hydrophobic core made up by side chains at the interface of the two beta sheets and the helix side facing the interior of the protein. The MC1 solution structure shows a globular protein with overall dimensions in the range of 34-40 A, which potentially corresponds to a DNA-binding site of 10-12 base pairs. The presumed DNA-binding site is located on the sequence comprising residues K62-P82, which is formed by a part of strands II2 and II3 belonging to the triple-stranded antiparallel beta sheet and a loop flanked by prolines P68 and P76. The tryptophan W74 that is expected to play a key role in the DNA-binding according to photocross-linking experiments was found completely exposed to the solvent, in a good position to interact with DNA. The overall fold of MC1, characterized by its linking beta-beta-alpha-beta-beta-loop-beta, is different from other known DNA-binding proteins. Its structure suggests a different DNA-binding mode than those of the histone-like proteins HU or HMGB. Thus, MC1 may be classified as a member of a new family.