Functional cloning, heterologous expression, and purification of two different N-deoxyribosyltransferases from Lactobacillus helveticus

Lactobacillus helveticus contains two types of N-deoxyribosyltransferases: DRTase I catalyzes the transfer of 2'-deoxyribose between purine bases exclusively whereas DRTase II is able to transfer the 2'-deoxyribose between two pyrimidine or between pyrimidine and purine bases. An Escherichia coli strain, auxotrophic for guanine and unable to use deoxyguanosine as source of guanine, was constructed to clone the corresponding genes. By screening a genomic bank for the production of guanine, the L. helveticus ptd and ntd genes coding for DRTase I and II, respectively, were isolated. Although the two genes have no sequence similarity, the two deduced polypeptides display 25.6% identity, with most of the residues involved in substrate binding and the active site nucleophile Glu-98 being conserved. Overexpression and purification of the two proteins shows that DRTase I is specific for purines with a preference for deoxyinosine (dI) > deoxyadenosine > deoxyguanosine as donor substrates whereas DRTase II has a strong preference for pyrimidines as donor substrates and purines as base acceptors. Purine analogues were substrates as acceptor bases for both enzymes. Comparison of DRTase I and DRTase II activities with dI as donor or hypoxanthine as acceptor and colocalization of the ptd and add genes suggest a specific role for DRTase I in the metabolism of dI.     

A code governing specific binding of regulatory proteins to DNA and structure of stereospecific sites of regulatory proteins]

A model is proposed for the structure of stereospecific sites in regulatory proteins. On its basis a possible code is suggested that governs the binding of regulatory proteins at specific control sites on DNA. Stereospecific sites of regulatory proteins are assumed to contain pairs of antiparallel polypeptide chain segments which form a right-hand twisted antiparallel beta-sheet, with single-stranded regions at the ends of the beta-structure. The model predicts that binding reaction between a regulatory protein and double-helical DNA is a cooperative phenomenon and is accompanied by significant structural alteration at the stereospecific site of the protein. Half of hydrogen bonds normally existing in beta-structure are broken upon complex formation with DNA and a new set of hydrogen bonds is formed between polypeptide amide groups and DNA base pairs. In a stereospecific site, one chain (t-chain) is attached through hydrogen bonds to the carbonyl oxygens of pyramides and N3 adenines lying in one DNA strand, while the second polypeptide chain (g chain) is hydrogen bonded to the 2-amino groups of guanine residues lying in the opposite DNA strand. The amide groups serve as specific reaction sites being hydrogen bond acceptors in g-chain and hydrogen bond donors in t-chain. The single-stranded portions of t- and g-chains lying in neighbouring subunits of regulatory protein interact with each other forming deformed beta-sheets. The recognition of regulatory sequences by proteins is based on the structural complementarity between stereospecific sites of regulatory proteins and base pairs sequences at the control sites. An essential feature of these sequences is the asymmetrical distribution of guanine residues between the two DNA strands. The code predicts that there are six fundamental amino acid residues (serine, threonine, asparagine, histidine, glutamine and cysteine) whose sequence in stereospecific site determines the base pair sequence to which a given regulatory protein would bind preferentially. The code states a correspondence between four amino acid residues at the stereospecific site of regulatory protein with the two residues being in t- and g-segments, respectively, and AT(GC) base pair at the control site. It is thus possible to determine which amino acid residues in the repressor and which base pairs in the operator DNA are involved in specific interactions with each other, as exemplified by lac repressor binding to lac operator.     

Environment of tryptophan side chains in proteins

Although relatively rare, the tryptophan residue (Trp), with its large hydrophobic surface, has a unique role in the folded structure and the binding site of many proteins, and its fluorescence properties make it very useful in studying the structures and dynamics of protein molecules in solution. An analysis has been made of its environment and the geometry of its interaction with neighbors using 719 Trp residues in 180 different protein structures. The distribution of the number of partners interacting with the Trp aromatic ring shows a peak at 6 (considering protein residues only) and 8 (including water and substrate molecules also). The means of the solvent-accessible surface areas of the ring show an exponential decrease with the increase in the number of partners; this relationship can be used to assess the efficiency of packing of residues around Trp. Various residues exhibit different propensities of binding the Trp side chain. The aromatic residues, Met and Pro have high values, whereas the smaller and polar-chain residues have weaker propensities. Most of the interactions are with residues far away in sequence, indicating the importance of Trp in stabilizing the tertiary structure. Of all the ring atoms NE1 shows the highest number of interactions, both along the edge (hydrogen bonding) as well as along the face. Various weak but specific interactions, engendering stability to the protein structure, have been identified.     

Classification of pseudo pairs between nucleotide bases and amino acids by analysis of nucleotide-protein complexes

Nucleotide bases are recognized by amino acid residues in a variety of DNA/RNA binding and nucleotide binding proteins. In this study, a total of 446 crystal structures of nucleotide-protein complexes are analyzed manually and pseudo pairs together with single and bifurcated hydrogen bonds observed between bases and amino acids are classified and annotated. Only 5 of the 20 usual amino acid residues, Asn, Gln, Asp, Glu and Arg, are able to orient in a coplanar fashion in order to form pseudo pairs with nucleotide bases through two hydrogen bonds. The peptide backbone can also form pseudo pairs with nucleotide bases and presents a strong bias for binding to the adenine base. The Watson-Crick side of the nucleotide bases is the major interaction edge participating in such pseudo pairs. Pseudo pairs between the Watson-Crick edge of guanine and Asp are frequently observed. The Hoogsteen edge of the purine bases is a good discriminatory element in recognition of nucleotide bases by protein side chains through the pseudo pairing: the Hoogsteen edge of adenine is recognized by various amino acids while the Hoogsteen edge of guanine is only recognized by Arg. The sugar edge is rarely recognized by either the side-chain or peptide backbone of amino acid residues.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Binding specificity of the two major DNA-binding proteins in human serum.

The two major DNA-binding proteins of human serum (DNA-binding protein 1 and DNA-binding protein 2) were shown to bind preferentially to single-stranded polynucleotides rich in guanine residues. Equilibrium competition experiments using a nitrocellulose filter assay system containing labeled human lymphocyte DNA and various competing natural and synthetic polynucleotides indicated that both proteins recognized sequences of bases containing a keto group in either position 6 (purines) or 4 (pyrimidines) and that these keto groups must be readily accessible for effective binding to occur. Guanine was shown to be the preferred nucleotide through inhibition experiments using a series of synthetic homopolymers and a series of bacterial DNAs of differing G + C content. The relationship between protein affinity and G + C content was shown to be directly proportional. The equilibrium constants for the binding of the human lymphocyte DNA by both proteins were on the order of 10(-6) M, and the length of the nucleotide sequence necessary for effective binding was found to be 12 to 18 bases using a series of oligomers of poly(dG).     

Learning about protein hydrogen bonding by minimizing contrastive divergence

Defining the strength and geometry of hydrogen bonds in protein structures has been a challenging task since early days of structural biology. In this article, we apply a novel statistical machine learning technique, known as contrastive divergence, to efficiently estimate both the hydrogen bond strength and the geometric characteristics of strong interpeptide backbone hydrogen bonds, from a dataset of structures representing a variety of different protein folds. Despite the simplifying assumptions of the interatomic energy terms used, we determine the strength of these hydrogen bonds to be between 1.1 and 1.5 kcal/mol, in good agreement with earlier experimental estimates. The geometry of these strong backbone hydrogen bonds features an almost linear arrangement of all four atoms involved in hydrogen bond formation. We estimate that about a quarter of all hydrogen bond donors and acceptors participate in these strong interpeptide hydrogen bonds.     

Insights into activation and RNA binding of trp RNA-binding attenuation protein (TRAP) through all-atom simulations

Tryptophan biosynthesis in Bacillus stearothermophilus is regulated by a trp RNA binding attenuation protein (TRAP). It is a ring-shaped 11-mer of identical 74 residue subunits. Tryptophan binding pockets are located between adjacent subunits, and tryptophan binding activates TRAP to bind RNA. Here, we report results from all-atom molecular dynamics simulations of the system, complementing existing extensive experimental studies. We focus on two questions. First, we look at the activation mechanism, of which relatively little is known experimentally. We find that the absence of tryptophan allows larger motions close to the tryptophan binding site, and we see indication of a conformational change in the BC loop. However, complete deactivation seems to occur on much longer time scales than the 40 ns studied here. Second, we study the TRAP-RNA interactions. We look at the relative flexibilities of the different bases in the complex and analyze the hydrogen bonds between the protein and RNA. We also study the role of Lys37, Lys56, and Arg58, which have been experimentally identified as essential for RNA binding. Hydrophobic stacking of Lys37 with the nearby RNA base is confirmed, but we do not see direct hydrogen bonding between RNA and the other two residues, in contrast to the crystal structure. Rather, these residues seem to stabilize the RNA-binding surface, and their positive charge may also play a role in RNA binding. Simulations also indicate that TRAP is able to attract RNA nonspecifically, and the interactions are quantified in more detail using binding energy calculations. The formation of the final binding complex is a very slow process: within the simulation time scale of 40 ns, only two guanine bases become bound (and no others), indicating that the binding initiates at these positions. In general, our results are in good agreement with experimental studies, and provide atomic-scale insights into the processes.     

Adaptive recognition in RNA complexes with peptides and protein modules

Recently, progress has been made towards the structural characterization of the novel folds of RNA-bound arginine-rich peptides and the architecture of their peptide-binding RNA pockets in viral and phage systems. These studies are based on an approach whereby the peptide and RNA components are minimalist modular domains that undergo adaptive structural transitions upon complex formation. Such complexes are characterized by recognition alignments in which the tertiary fold of the RNA generates binding pockets with the potential to envelop minimal elements of protein secondary structure. Strikingly, the peptides fold as isolated alpha-helical or beta-hairpin folds within their RNA major-groove targets, without the necessity of additional appendages for anchorage within the binding pocket. The RNA peptide-binding pocket architectures are sculptured through precisely positioned mismatches, triples and looped-out bases, which accommodate amino acid sidechains through hydrophobic, hydrogen bonding and ionic intermolecular contacts. By contrast, protein modules associated with the HIV-1 nucleocapsid and MS2 phage coat target their RNA binding sites through the insertion of specificity-determining RNA base residues within conserved hydrophobic pockets and crevices on the protein surface, with the bases anchored through hydrogen bonding interactions. These alternative strategies of RNA recognition at the peptide and protein module level provide novel insights into the principles, patterns and diversity of the adaptive transitions associated with the recognition process.     

Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

On the Information Content of Protein Sequences

Abstract

TATA-box binding protein (TBP) in a monomelic form and the complexes it forms with DNA have been elucidated with molecular dynamics simulations. Large TBP domain motions (bend and twist) are detected in the monomer as well as in the DNA complexes; these motions can be important for TBP binding of DNA. TBP interacts with guanine bases flanking the TATA element in the simulations of the complex; these interactions may explain the preference for guanine observed at these DNA positions. Side chains of some TBP residues at the binding interface display significant dynamic flexibility that results in ‘flipflop’ contacts involving multiple base pairs of the DNA. We discuss the possible functional significance of these observations.     

Molecular Complexes of acridine orange and nucleosides

Summary The addition of various nucleosides to the aqueous AO solution brings about the red shift of the absorption band of AO monomer and the enhancement of the AO fluorescence emission. These phenomena are attributable to the formation of a kind of molecular complex between AO monomer and nucleoside.The absorption and fluorescence characteristics and their thermal behaviours enable us to determine the association constant and the binding energy. The binding energy of AO with purines is larger than that with pyrimidines, and the association constant between AO and the deoxyribonucleoside is larger than that between AO and the ribonucleoside. For the molecular complexes dealt with, the face-to-face arrangement of AO and nucleoside, linking of AO with sugar by hydrogen bridge, may be more preferential than the side-by-side arrangement with the direct linkage between AO and nucleic acid base by hydrogen bonding. The Van der Waals-London interactions may be one of the essential factors for the binding in these molecular complexes.The association constant and the binding energy for AO-DNA and -RNA systems were also determined. The magnitude of these quantities seems to reflect the difference in the structure of these nucleic acids. The rather open structure of RNA compared with DNA is in favour of affinity with AO and gives the larger association constant than that for AO-DNA. The binding energy is somewhat larger for AO-DNA complex than for AO-RNA complex, probably due to the structural difference of the base arrangement; the stacked base pairs for DNA and the stacked bases for RNA.This may explain the selective degradation [12] of guanine photo-sensitized by dyes either in the free state or when incorporated to DNA or RNA, as discussed elsewhere.     

Further characterization of a depurinated DNA-purine base insertion activity from cultured human fibroblasts.

The purification from cultured human fibroblasts of a protein that binds specifically to partially depurinated DNA and inserts purines into those sites is described. The purine insertion, but not the binding, requires K+. The DNA binding can be saturated with increasing apurinic sites and is weakened by the presence of adenine or guanine. Base insertion into depurinated DNA is specific for adenine or guanine; none is observed with dATP or dGTP. When the depurinated DNA substrate is specifically cleaved with apurinic endonuclease, no purine insertion occurs. Guanine insertion does not occur into tRNA or depyrimidinated DNA, and thymine is not inserted into either depyrimidinated DNA or depurinated DNA. Purine insertion activity follows Michaelis-Menten kinetics with respect to purintes; the apparent Km values for both adenine and guanine are 5 microM. The enzyme binds the purine bases very tightly. Adenine binding saturates at less than 1 microM adenine, perhaps reflecting the low intracellular adenine concentration. The binding protein specific for UV-irradiated DNA (Feldberg, R.S., and Grossman, L. (1976) Biochemistry 15, 2402-2408) had no detectable purine or pyrimidine base insertion activity with depurinated or depyrimidinated DNAs.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Relationship between intramolecular hydrogen bonding and solvent accessibility of side-chain donors and acceptors in proteins

This study shows that intramolecular hydrogen bonding in proteins depends on the accessibility of donors and acceptors to water molecules. The frequency of occurrence of H-bonded side chains in proteins is inversely proportional to the solvent accessibility of their donors and acceptors. Estimates of the notional free energy of hydrogen bonding suggest that intramolecular hydrogen-bonding interactions of buried and half-buried donors and acceptors can contribute favorably to the stability of a protein, whereas those of solvent-exposed polar atoms become less favorable or unfavorable.     

Host-guest scale of left-handed polyproline II helix formation

Despite the clear importance of the left-handed polyproline II (PPII) helical conformation in many physiologically important processes as well as its potential significance in protein unfolded states, little is known about the physical determinants of this conformation. We present here a scale of relative PPII helix-forming propensities measured for all residues, except tyrosine and tryptophan, in a proline-based host peptide system. Proline has the highest measured propensity in this system, a result of strong steric interactions that occur between adjacent prolyl rings. The other measured propensities are consistent with backbone solvation being an important component in PPII helix formation. Side chain to backbone hydrogen bonding may also play a role in stabilizing this conformation. The PPII helix-forming propensity scale will prove useful in future studies of the conformational properties of proline-rich sequences as well as provide insights into the prevalence of PPII helices in protein unfolded states.     

The 2.2 A crystal structure of Hsp33: a heat shock protein with redox-regulated chaperone activity

BACKGROUND: One strategy that cells employ to respond to environmental stresses (temperature, oxidation, and pathogens) is to increase the expression of heat shock proteins necessary to maintain viability. Several heat shock proteins function as molecular chaperones by binding unfolded polypeptides and preventing their irreversible aggregation. Hsp33, a highly conserved bacterial heat shock protein, is a redox-regulated molecular chaperone that appears to protect cells against the lethal effects of oxidative stress. RESULTS: The 2.2 A crystal structure of a truncated E. coli Hsp33 (residues 1-255) reveals a domain-swapped dimer. The core domain of each monomer (1-178) folds with a central helix that is sandwiched between two beta sheets. The carboxyl-terminal region (179-235), which lacks the intact Zn binding domain of Hsp33, folds into three helices that pack on the other subunit. The interface between the two core domains is comprised of conserved residues, including a rare Glu-Glu hydrogen bond across the dyad axis. Two potential polypeptide binding sites that span the dimer are observed: a long groove containing pockets of conserved and hydrophobic residues, and an intersubunit 10-stranded beta sheet "saddle" with a largely uncharged or hydrophobic surface. CONCLUSIONS: Hsp33 is a dimer in the crystal structure. Solution studies confirmed that this dimer reflects the structural changes that occur upon activation of Hsp33 as a molecular chaperone. Patterns of conserved residues and surface charges suggest that two grooves might be potential binding sites for protein folding intermediates.     

Molecular modeling study of the norfloxacin-DNA complex

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5'-d(ATACGTAT)(2). Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.