Structural segments and residue propensities in protein-RNA interfaces: comparison with protein-protein and protein-DNA complexes

The interface of a protein molecule that is involved in binding another protein, DNA or RNA has been characterized in terms of the number of unique secondary structural segments (SSSs), made up of stretches of helix, strand and non-regular (NR) regions. On average 10-11 segments define the protein interface in protein-protein (PP) and protein-DNA (PD) complexes, while the number is higher (14) for protein-RNA (PR) complexes. While the length of helical segments in PP interaction increases with the interface area, this is not the case in PD and PR complexes. The propensities of residues to occur in the three types of secondary structural elements (SSEs) in the interface relative to the corresponding elements in the protein tertiary structures have been calculated. Arg, Lys, Asn, Tyr, His and Gln are preferred residues in PR complexes; in addition, Ser and Thr are also favoured in PD interfaces.     

Thermodynamic and electrostatic properties of ternary Oxytricha nova TEBP-DNA complex

Telomeres constitute the nucleoprotein ends of eukaryotic chromosomes which are essential for their proper function. Telomere end binding protein (TEBP) from Oxytricha nova was among the first telomeric proteins, which were well characterized biologically. TEBP consists of two protein subunits (alpha, beta) and forms a ternary complex with single stranded telomeric DNA containing tandem repeats TTTTGGGG. This work presents the characterization of the thermodynamic and electrostatic properties of this complex by computational chemistry methods (continuum Poisson-Boltzmann and solvent accessible surface calculations). Our calculations give a new insight into molecular properties of studied system. Based on the thermodynamic analysis we provide a rationale for the experimental observation that alpha and ssDNA forms a binary complex and the beta subunit joins alpha:ssDNA complex only after the latter is formed. Calculations of distribution of the molecular electrostatic potential for protein subunits alone and for all possible binary complexes revealed the important role of the "guiding funnel" potential generated by alpha:ssDNA complex. This potential may help the beta subunit to dock to the already formed alpha:DNA intermediate in highly steric and electrostatic favorable manner. Our pK(a) calculations of TEBP are able to explain the experimental mobility shifts of the complex in electrophoretic non-denaturating gels.     

Differential affinity and cooperativity functions of the amino-terminal 70 residues of lambda integrase

The site-specific recombinase (Int) of bacteriophage lambda is a heterobivalent DNA-binding protein that binds two different classes of DNA-binding sites within its recombination target sites. The several functions of Int are apportioned between a large carboxy-terminal domain that cleaves and ligates DNA at each of its four "core-type" DNA-binding sites and a small amino-terminal domain, whose primary function is binding to each of its five "arm-type" DNA sites, which are distant from the core region. Int bridges between the two classes of binding sites are facilitated by accessory DNA-bending proteins that along with Int comprise higher-order recombinogenic complexes. We show here that although the 64 amino-terminal residues of Int bind efficiently to a single arm site, this protein cannot form doubly bound complexes on adjacent arm sites. However, 1-70 Int does show the same cooperative binding to adjacent arm sites as the full length protein. We also found that 1-70 Int specifies cooperative interactions with the accessory protein Xis when the two are bound to their adjacent cognate sites P2 and X1, respectively. To complement the finding that these two amino-terminal domain functions (along with arm DNA binding) are all specified by residues 1-70, we determined that Thr75 is the first residue of the minimal carboxy-terminal domain, thereby identifying a specific interdomain linker region. We have measured the affinity constants for Int binding to each of the five arm sites and the cooperativity factors for Int binding to the two pairs of adjacent arm sites, and we have identified several DNA structural features that contribute to the observed patterns of Int binding to arm sites. Taken together, the results highlight several interesting features of arm DNA binding that invite speculation about additional levels of complexity in the regulation of lambda site-specific recombination.     

Nucleotide variability and linkage disequilibrium patterns at the porcine FABP5 gene

Fatty acid binding protein 5 (FABP5) is a major positional and physiological candidate gene for the porcine FAT1 QTL on SSC4. Here we characterize the nucleotide polymorphism and haplotype variability of FABP5 and we compare it with that of FABP4, given their close physical location and similar metabolic roles. DNA resequencing of the FABP5 gene region in 29 pigs from 14 breeds and in European and Japanese wild boars revealed 36 polymorphisms in 5.2 kb, and a nucleotide diversity of 0.19%, comparable to values reported in other domestic species but sixfold lower than that previously found for FABP4. Remarkably, both the nucleotide variability and the haplotype structure of FABP5 and FABP4 were dramatically different, and the Hudson-Kreitman-Aguadé test was highly significant. Nevertheless, both genes also had similarities. The neighbour-joining trees of their haplotypes did not show a geographical arrangement for any of the genes. Besides, both genes presented a similar extent and pattern of linkage disequilibrium. Haplotype blocks did not extend for large stretches ( approximately 1 kb in both genes), and the number of tag SNPs required to capture all variability was higher than previously expected. Our findings indicate that FABP4 and FABP5 have undergone different selective or evolutive processes. The fact that haplotype blocks were so small may require us to increase the number of SNPs in prospective whole-genome association studies in the pig.     

Gene-sized macronuclear DNA molecules are clustered in micronuclear chromosomes of the ciliate Oxytricha nova.

Following the sexual phase of its life cycle, the hypotrichous ciliate Oxytricha nova transforms a copy of its chromosomal micronucleus into a macronucleus containing short, linear DNA molecules with an average size of 2.2 kilobase pairs. In addition, more than 90% of the DNA sequences in the micronuclear genome are eliminated during this process. We have examined the organization of macronuclear DNA molecules in the micronuclear chromosomes. Macronuclear DNA molecules were found to be clustered and separated by less than 550 base pairs in two cloned segments of micronuclear DNA. Recombinant clones of two macronuclear DNA molecules that are adjacent in the micronucleus were also isolated and examined by DNA sequencing. The two macronuclear DNA molecules were found to be separated by only 90 base pairs in the micronuclear genome.     

Close linkage between the albumin and Gc loci in the horse

Evidence for close linkage between the structural loci for albumin and Gc protein in the horse was presented. A recombination frequency (c) of 0.009 ± 0.006 (95 % confidence limits: 0.001 < c < 0.032) was estimated. These results were based on a study of a large sire family comprising 223 offspring from informative matings. No evidence of linkage disequilibrium was observed in one horse population studied.     

Thermodynamic characterization of binding Oxytricha nova single strand telomere DNA with the alpha protein N-terminal domain

The Oxytricha nova telemere binding protein alpha subunit binds single strand DNA and participates in a nucleoprotein complex that protects the very ends of chromosomes. To understand how the N-terminal, DNA binding domain of alpha interacts with DNA we measured the stoichiometry, enthalpy (DeltaH), entropy (DeltaS), and dissociation constant (K(D-DNA)) for binding telomere DNA fragments at different temperatures and salt concentrations using native gel electrophoresis and isothermal titration calorimetry (ITC). About 85% of the total free energy of binding corresponded with non-electrostatic interactions for all DNAs. Telomere DNA fragments d(T(2)G(4)), d(T(4)G(4)), d(G(3)T(4)G(4)), and d(G(4)T(4)G(4)) each formed monovalent protein complexes. In the case of d(T(4)G(4)T(4)G(4)), which has two tandemly repeated d(TTTTTGGGG) telomere motifs, two binding sites were observed. The high-affinity "A site" has a dissociation constant, K(D-DNA(A)) = 13(+/-4) nM, while the low-affinity "B site" is characterized by K(D-DNA(B)) = 5600(+/-600) nM at 25 degrees C. Nucleotide substitution variants verified that the A site corresponds principally with the 3'-terminal portion of d(T(4)G(4)T(4)G(4)). The relative contributions of entropy (DeltaS) and enthalpy (DeltaH) for binding reactions were DNA length-dependent as was heat capacity (DeltaCp). These trends with respect to DNA length likely reflect structural transitions in the DNA molecule that are coupled with DNA-protein association. Results presented here are important for understanding early intermediates and subsequent stages in the assembly of the full telomere nucleoprotein complex and how binding events can prepare the telomere DNA for extension by telomerase, a critical event in telomere biology.     

A simple physical model for the prediction and design of protein-DNA interactions

Protein-DNA interactions are crucial for many biological processes. Attempts to model these interactions have generally taken the form of amino acid-base recognition codes or purely sequence-based profile methods, which depend on the availability of extensive sequence and structural information for specific structural families, neglect side-chain conformational variability, and lack generality beyond the structural family used to train the model. Here, we take advantage of recent advances in rotamer-based protein design and the large number of structurally characterized protein-DNA complexes to develop and parameterize a simple physical model for protein-DNA interactions. The model shows considerable promise for redesigning amino acids at protein-DNA interfaces, as design calculations recover the amino acid residue identities and conformations at these interfaces with accuracies comparable to sequence recovery in globular proteins. The model shows promise also for predicting DNA-binding specificity for fixed protein sequences: native DNA sequences are selected correctly from pools of competing DNA substrates; however, incorporation of backbone movement will likely be required to improve performance in homology modeling applications. Interestingly, optimization of zinc finger protein amino acid sequences for high-affinity binding to specific DNA sequences results in proteins with little or no predicted specificity, suggesting that naturally occurring DNA-binding proteins are optimized for specificity rather than affinity. When combined with algorithms that optimize specificity directly, the simple computational model developed here should be useful for the engineering of proteins with novel DNA-binding specificities.     

Ecotin: lessons on survival in a protease-filled world.

Ecotin, an Escherichia coli periplasmic protein of 142 amino acids, has been shown to be a potent inhibitor of a group of homologous serine proteases with widely differing substrate recognition. It is highly effective against a number of enzymes, including both pancreatic and neutrophil-derived elastases, chymotrypsin, trypsin, factor Xa, and kallikrein. Recent structural and functional studies on ecotin and its interactions with different serine proteases have clarified these initial observations and revealed the remarkable features of this protein in inhibiting a strikingly large subset of the chymotrypsin family of serine proteases. The structures of the ecotin:serine protease complexes provide the first examples of protein-protein recognition where the concept of specificity of interactions needs to be reexamined. The binding sites show a fluidity of protein contacts derived from ecotin's innate flexibility in fitting itself to proteases while strongly interfering with their function.     

Identification of the ubiquitin interfacial residues in a ubiquitin-E2 covalent complex

One of the key intermediates formed during the protein ubiquitination cycle is a covalent complex between ubiquitin (Ub) and the conjugation enzyme, UBC1. In order to probe the interface between these two proteins we have formed the covalent complex in situ (in the NMR tube) using Ub, the catalytic domain of UBC1, UBC1450, an activation enzyme, E1, and Mg²⁺-ATP. The size of the Ub-UBC1450 complex (25 kDa) and its relatively short lifetime ( 4 h) makes assignment of the backbone resonances in the covalent species difficult. In order to monitor the formation and identify the interface in the complex we have used fast ¹H-¹⁵N HSQC spectra to monitor the decay of ¹H-¹⁵N correlations as a function of time until the complex formed reached about 90%. The residual peak intensities were used to probe the surface of interaction between Ub and UBC1450 and provided a clear surface of interaction on Ub.     

Capping complex formation at the slow-growing end of the actin filament

Actin filaments are polar; their barbed (fast-growing) and pointed (slow-growing) ends differ in structure and dynamic properties. The slow-growing end is regulated by tropomodulins, a family of capping proteins that require tropomyosins for optimal function. There are four tropomodulin isoforms; their distributions vary depending on tissue type and change during development. The C-terminal half of tropomodulin contains one compact domain represented by alternating α-helices and β-structures. The tropomyosin-independent actin-capping site is located at the C-terminus. The N-terminal half has no regular structure; however, it contains a tropomyosin-dependent actin-capping site and two tropomyosin-binding sites. One tropomodulin molecule can bind two tropomyosin molecules. Effectiveness of tropomodulin binding to tropomyosin depends on the tropomyosin isoform. Regulation of tropomodulin binding at the pointed end as well as capping effectiveness in the presence of specific tropomyosins may affect formation of local cytoskeleton and dynamics of actin filaments in cells.     

Solid-state 23Na NMR determination of the number and coordination of sodium cations bound to Oxytricha nova telomere repeat d(G4T4G4)

We report a solid-state (23)Na NMR study of the bound sodium cations in a G-quadruplex formed by Oxytricha nova telomere DNA repeat, d(G(4)T(4)G(4)) (Oxy-1.5). Using a 2D multiple-quantum magic-angle spinning (23)Na NMR method, we observed three sodium cations residing inside the quadruplex channel of the Na(+) form of Oxy-1.5. Each of these sodium cations is sandwiched between two G-quartets. We found no evidence for sodium cations in the T(4) loop region. For comparison, solid-state (15)N MAS NMR spectra were also obtained for the (15)NH(4)(+) form of Oxy-1.5. The insufficient resolution in the (15)N MAS NMR spectra did not permit determination of the number of NH(4)(+) ions inside the quadruplex channel. The solid-state (23)Na and (15)N NMR spectra for Oxy-1.5 were also compared with those obtained for guanosine 5'-monophosphate.     

DNA G-quartets in a 1.86 A resolution structure of an Oxytricha nova telomeric protein-DNA complex.

The Oxytricha nova telomere end binding protein (OnTEBP) recognizes, binds and protects the single-stranded 3'-terminal DNA extension found at the ends of macronuclear chromosomes. The structure of this complex shows that the single strand GGGGTTTTGGGG DNA binds in a deep cleft between the two protein subunits of OnTEBP, adopting a non-helical and irregular conformation. In extending the resolution limit of this structure to 1.86 A, we were surprised to find a G-quartet linked dimer of the GGGGTTTTGGGG DNA also packing within the crystal lattice and interacting with the telomere end binding protein. The G-quartet DNA exhibits the same structure and topology as previously observed in solution by NMR with diagonally crossing d(TTTT) loops at either end of the four-stranded helix. Additionally, the crystal structure reveals clearly visible Na(+), and specific patterns of bound water molecules in the four non-equivalent grooves. Although the G-quartet:protein contact surfaces are modest and might simply represent crystal packing interactions, it is interesting to speculate that the two types of telomeric DNA-protein interactions observed here might both be important in telomere biology.     

Electrostatic potential distribution of the gene V protein from Ff phage facilitates cooperative DNA binding: a model of the GVP-ssDNA complex.

The crystal structure of the gene V protein (GVP) from the Ff filamentous phages (M13, fl, fd) has been solved for the wild-type and two mutant (Y41F and Y41H) proteins at high resolution. The Y41H mutant crystal structure revealed crystal packing interactions, which suggested a plausible scheme for constructing the polymeric protein shell of the GVP-single-stranded DNA (ssDNA) complex (Guan Y, et al., 1994, Biochemistry 33:7768-7778). The electrostatic potentials of the isolated and the cooperatively formed protein shell have been calculated using the program GRASP and they revealed a highly asymmetric pattern of the electrostatic charge distribution. The inner surface of the putative DNA-binding channel is positively charged, whereas the opposite outer surface is nearly neutral. The electrostatic calculation further demonstrated that the formation of the helical protein shell enhanced the asymmetry of the electrostatic distribution. A model of the GVP-ssDNA complex with the n = 4 DNA-binding mode could be built with only minor conformational perturbation to the GVP protein shell. The model is consistent with existing biochemical and biophysical data and provides clues to the properties of GVP, including the high cooperatively of the protein binding to ssDNA. The two antiparallel ssDNA strands form a helical ribbon with the sugar-phosphate backbones at the middle and the bases pointing away from each other. The bases are stacked and the Phe 73 residue is intercalated between two bases. The optimum binding to a tetranucleotide unit requires the participation of four GVP dimers, which may explain the cooperativity of the GVP binding to DNA.     

Specific association of a small protein with the telomeric DNA-protein complex during the onset of leaf senescence in Arabidopsis thaliana

Telomeres and their changes in length throughout the life span of cells have been intensively investigated in different organisms. Telomere length is assumed to control replicative senescence in mammalian cells. However, only very few data are available on the developmental dynamics of plant telomeres. Here, changes of telomere length and DNA-protein structure of Arabidopsis thaliana telomeres were analysed in different stages of development, with the main focus resting on the transition from pre-senescent to senescent leaves. The lengths of the telomeres, ranging from ca. 2.0 to 6.5 kb, do not significantly change during plant development indicating that telomere length is not involved in differentiation and replicative senescence nor in post-mitotic senescence of A. thaliana. In dedifferentiated cultured cells a slight increase in length can be determined. The nucleoprotein structure of the telomeric DNA was investigated by gel mobility shift assays, with synthetic oligonucleotides and nuclear protein extracts derived from four defined stages of post-mitotic leaf senescence. In all four stages, a highly salt-resistant DNA-protein complex was formed with the double-stranded as well as with the single-stranded G-rich telomeric DNA. An additional DNA-protein complex was identified in nuclear protein extracts isolated from plants in the transition stage from pre-senescence to senescence. The protein components of the DNA-protein complexes were analysed on native PAGE and SDS-PAGE gels. A protein of 67 kDa (ATBP1) bound to the telomeric DNA in all developmental stages. An additional protein of merely 22 kDa (ATBP2) was associated via protein-protein interaction with ATBP1 to form a higher-order complex exclusively during the onset of senescence. DNA interaction of this higher-order protein complex seems to be restricted to double-stranded telomeric DNA. The defined period of ATBP1/ATBP2 complex formation with the telomeric DNA probably indicates that ATBP2 is involved in the onset of post-mitotic leaf senescence by either disturbing an established or establishing an additional function exhibited by the telomeres in the interphase nuclei.     

Electrostatic interactions in the association of proteins: an analysis of the thrombin-hirudin complex.

The role of electrostatic interactions in stabilization of the thrombin-hirudin complex has been investigated by means of two macroscopic approaches: the modified Tanford-Kirkwood model and the finite-difference method for numerical solution of the Poisson-Boltzmann equations. The electrostatic potentials around the thrombin and hirudin molecules were asymmetric and complementary, and it is suggested that these fields influence the initial orientation in the process of the complex formation. The change of the electrostatic binding energy due to mutation of acidic residues in hirudin has been calculated and compared with experimentally determined changes in binding energy. In general, the change in electrostatic binding energy for a particular mutation calculated by the modified Tanford-Kirkwood approach agreed well with the experimentally observed change. The finite-difference approach tended to overestimate changes in binding energy when the mutated residues were involved in short-range electrostatic interactions. Decreases in binding energy caused by mutations of amino acids that do not make any direct ionic interactions (e.g., Glu 61 and Glu 62 of hirudin) can be explained in terms of the interaction of these charges with the positive electrostatic potential of thrombin. Differences between the calculated and observed changes in binding energy are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

Telomere terminal transferase activity in the hypotrichous ciliate Oxytricha nova and a model for replication of the ends of linear DNA molecules.

We have found abundant telomere-specific terminal transferase activity in crude macronuclear extracts from vegetatively growing cells of the hypotrichous ciliate Oxytricha nova. This activity adds two to seven tandem repeats of the sequence GGGGTTTT (the Oxytricha telomeric repeat) to the 3' end of oligonucleotide primers ending in repeats of G4T4 and always adds the repeats in the proper phase. The activity requires the presence of micromolar amounts of dGTP and dTTP as well as single-stranded oligomer primers ending 3' with repeats of the Oxytricha telomeric sequence. A nuclease activity is present in the extracts which is closely balanced with telomere terminal transferase activity. We propose a simple model for replication of the ends of linear DNA molecules based on the telomere terminal transferase.