Predicting the effects of basepair mutations in DNA-protein complexes by thermodynamic integration

AbstractThermodynamically rigorous free energy methods in principle allow the exact computation of binding free energies in biological systems. Here, we use thermodynamic integration together with molecular dynamics simulations of a DNA-protein complex to compute relative binding free energies of a series of mutants of a protein-binding DNA operator sequence. A guanine-cytosine basepair that interacts strongly with the DNA-binding protein is mutated into adenine-thymine, cytosine-guanine, and thymine-adenine. It is shown that basepair mutations can be performed using a conservative protocol that gives error estimates of ∼10% of the change in free energy of binding. Despite the high CPU-time requirements, this work opens the exciting opportunity of being able to perform basepair scans to investigate protein-DNA binding specificity in great detail computationally

Structural and mutagenic analysis of the RM controller protein C.Esp1396I

Bacterial restriction-modification (RM) systems are comprised of two complementary enzymatic activities that prevent the establishment of foreign DNA in a bacterial cell: DNA methylation and DNA restriction. These two activities are tightly regulated to prevent over-methylation or auto-restriction. Many Type II RM systems employ a controller (C) protein as a transcriptional regulator for the endonuclease gene (and in some cases, the methyltransferase gene also). All high-resolution structures of C-protein/DNA-protein complexes solved to date relate to C.Esp1396I, from which the interactions of specific amino acid residues with DNA bases and/or the phosphate backbone could be observed. Here we present both structural and DNA binding data for a series of mutations to the key DNA binding residues of C.Esp1396I. Our results indicate that mutations to the backbone binding residues (Y37, S52) had a lesser affect on DNA binding affinity than mutations to those residues that bind directly to the bases (T36, R46), and the contributions of each side chain to the binding energies are compared. High-resolution X-ray crystal structures of the mutant and native proteins showed that the fold of the proteins was unaffected by the mutations, but also revealed variation in the flexible loop conformations associated with DNA sequence recognition. Since the tyrosine residue Y37 contributes to DNA bending in the native complex, we have solved the structure of the Y37F mutant protein/DNA complex by X-ray crystallography to allow us to directly compare the structure of the DNA in the mutant and native complexes

Radiation damage to nucleoprotein complexes in macromolecular crystallography

Significant progress has been made in macromolecular crystallography over recent years in both the understanding and mitigation of X-ray induced radiation damage when collecting diffraction data from crystalline proteins. In contrast, despite the large field that is productively engaged in the study of radiation chemistry of nucleic acids, particularly of DNA, there are currently very few X-ray crystallographic studies on radiation damage mechanisms in nucleic acids. Quantitative comparison of damage to protein and DNA crystals separately is challenging, but many of the issues are circumvented by studying pre-formed biological nucleoprotein complexes where direct comparison of each component can be made under the same controlled conditions. Here a model protein-DNA complex C.Esp1396I is employed to investigate specific damage mechanisms for protein and DNA in a biologically relevant complex over a large dose range (2.07-44.63 MGy). In order to allow a quantitative analysis of radiation damage sites from a complex series of macromolecular diffraction data, a computational method has been developed that is generally applicable to the field. Typical specific damage was observed for both the protein on particular amino acids and for the DNA on, for example, the cleavage of base-sugar N1-C and sugar-phosphate C-O bonds. Strikingly the DNA component was determined to be far more resistant to specific damage than the protein for the investigated dose range. At low doses the protein was observed to be susceptible to radiation damage while the DNA was far more resistant, damage only being observed at significantly higher doses

Identification of a structural and functional domain in xNAP1 involved in protein–protein interactions

xNAP1 (Xenopus nucleosome assembly protein) belongs to the family of nucleosome assembly proteins (NAPs) and shares 92% identity with human and mouse NAP1. NAPs have been reported to have a role in nucleosome assembly, cell cycle regulation, cell proliferation and transcriptional control, although the precise function of NAP1 is still not clear. Here we report the identification of a putative domain of xNAP1 by limited proteolysis. This domain has been mapped in the xNAP1 protein sequence to residues 38–282 and thus lacks the acidic sequences at the N- and C-termini. We have studied this domain and related fragments in vitro and by a functional assay involving over-expression of the protein in Xenopus laevis embryos. Analytical ultracentrifugation shows that removal of the acidic N- and C-terminal regions does not prevent the formation of larger multimers, which are predominantly hexadecamers. Injection of mRNA encoding the full-length xNAP1 or the putative domain and other related constructs into Xenopus embryos gave identical phenotypes. These results are discussed in relation to protein–protein interactions between NAP1 octamers and a possible ‘squelching’ mechanism

Identification of a structural and functional domain in xNAP1 involved in protein–protein interactions

xNAP1 (Xenopus nucleosome assembly protein) belongs to the family of nucleosome assembly proteins (NAPs) and shares 92% identity with human and mouse NAP1. NAPs have been reported to have a role in nucleosome assembly, cell cycle regulation, cell proliferation and transcriptional control, although the precise function of NAP1 is still not clear. Here we report the identification of a putative domain of xNAP1 by limited proteolysis. This domain has been mapped in the xNAP1 protein sequence to residues 38–282 and thus lacks the acidic sequences at the N- and C-termini. We have studied this domain and related fragments in vitro and by a functional assay involving over-expression of the protein in Xenopus laevis embryos. Analytical ultracentrifugation shows that removal of the acidic N- and C-terminal regions does not prevent the formation of larger multimers, which are predominantly hexadecamers. Injection of mRNA encoding the full-length xNAP1 or the putative domain and other related constructs into Xenopus embryos gave identical phenotypes. These results are discussed in relation to protein–protein interactions between NAP1 octamers and a possible ‘squelching’ mechanism

Transcription regulation of the type II restriction-modification system AhdI

The Restriction-modification system AhdI contains two convergent transcription units, one with genes encoding methyltransferase subunits M and S and another with genes encoding the controller (C) protein and the restriction endonuclease (R). We show that AhdI transcription is controlled by two independent regulatory loops that are well-optimized to ensure successful establishment in a naïve bacterial host. Transcription from the strong MS promoter is attenuated by methylation of an AhdI site overlapping the -10 element of the promoter. Transcription from the weak CR promoter is regulated by the C protein interaction with two DNA-binding sites. The interaction with the promoter-distal high-affinity site activates transcription, while interaction with the weaker promoter-proximal site represses it. Because of high levels of cooperativity, both C protein-binding sites are always occupied in the absence of RNA polymerase, raising a question how activated transcription is achieved. We develop a mathematical model that is in quantitative agreement with the experiment and indicates that RNA polymerase outcompetes C protein from the promoter-proximal-binding site. Such an unusual mechanism leads to a very inefficient activation of the R gene transcription, which presumably helps control the level of the endonuclease in the cell

Transcription regulation of the type II restriction-modification system AhdI

The Restriction-modification system AhdI contains two convergent transcription units, one with genes encoding methyltransferase subunits M and S and another with genes encoding the controller (C) protein and the restriction endonuclease (R). We show that AhdI transcription is controlled by two independent regulatory loops that are well-optimized to ensure successful establishment in a naïve bacterial host. Transcription from the strong MS promoter is attenuated by methylation of an AhdI site overlapping the -10 element of the promoter. Transcription from the weak CR promoter is regulated by the C protein interaction with two DNA-binding sites. The interaction with the promoter-distal high-affinity site activates transcription, while interaction with the weaker promoter-proximal site represses it. Because of high levels of cooperativity, both C protein-binding sites are always occupied in the absence of RNA polymerase, raising a question how activated transcription is achieved. We develop a mathematical model that is in quantitative agreement with the experiment and indicates that RNA polymerase outcompetes C protein from the promoter-proximal-binding site. Such an unusual mechanism leads to a very inefficient activation of the R gene transcription, which presumably helps control the level of the endonuclease in the cell

Transcription regulation of restriction-modification system Esp1396I

The convergently transcribed restriction (R) and methylase (M) genes of the Restriction–Modification system Esp1396I are tightly regulated by a controller (C) protein that forms part of the CR operon. We have mapped the transcriptional start sites from each promoter and examined the regulatory role of C.Esp1396I in vivo and in vitro. C-protein binding at the CR and M promoters was analyzed by DNA footprinting and a range of biophysical techniques. The distal and proximal C-protein binding sites at the CR promoter are responsible for activation and repression, respectively. In contrast, a C-protein dimer binds to a single site at the M-promoter to repress the gene, with an affinity much greater than for the CR promoter. Thus, during establishment of the system in a naïve host, the activity of the M promoter is turned off early, preventing excessive synthesis of methylase. Mutational analysis of promoter binding sites reveals that the tetranucleotide inverted repeats long believed to be important for C-protein binding to DNA are less significant than previously thought. Instead, symmetry-related elements outside of these repeats appear to be critical for the interaction and are discussed in terms of the recent crystal structure of C.Esp139I bound to the CR promoter