Probing contacts to the DNA backbone in the trp repressor-operator sequence-specific protein-nucleic acid complex using diastereomeric methylphosphonate analogues

Corepressor effect on androgen receptor activity varies with the length of the CAG encoded polyglutamine repeat and is dependent on receptor/corepressor ratio in prostate cancer cells

The response of prostate cells to androgens reflects a combination of androgen receptor (AR) transactivation and transrepression, but how these two processes differ mechanistically and influence prostate cancer risk and disease outcome remain elusive. Given recent interest in targeting AR transrepressive processes, a better understanding of AR/corepressor interaction and responses is warranted. Here, we used transactivation and interaction assays with wild-type and mutant ARs, and deletion AR fragments, to dissect the relationship between AR and the corepressor, silencing mediator for retinoic acid and thyroid hormone receptors (SMRT). We additionally tested how these processes are influenced by AR agonist and antagonist ligands, as well as by variation in the polyglutamine tract in the AR amino terminal domain (NTD), which is encoded by a polymorphic CAG repeat in the gene. SMRT was recruited to the AR ligand binding domain by agonist ligand, and as determined by the effect of strategic mutations in activation function 2 (AF-2), requires a precise conformation of that domain. A distinct region of SMRT also mediated interaction with the AR-NTD via the transactivation unit 5 (TAU5; residues 315-538) region. The degree to which SMRT was able to repress AR increased from 17% to 56% as the AR polyglutamine repeat length was increased from 9 to 42 residues, but critically this effect could be abolished by increasing the SMRT:AR molar ratio. These data suggest that the extent to which the CAG encoded polyglutamine repeat influences AR activity represents a balance between corepressor and coactivator occupancy of the same ligand-dependent and independent AR interaction surfaces. Changes in the homeostatic relationship of AR to these molecules, including SMRT, may explain the variable penetrance of the CAG repeat and the loss of AR signaling flexibility in prostate cancer progression.

Evolution of acceptor stem tRNA recognition by class II prolyl-tRNA synthetase

Aminoacyl-tRNA synthetases (AARS) are an essential family of enzymes that catalyze the attachment of amino acids to specific tRNAs during translation. Previously, we showed that base-specific recognition of the tRNA(Pro) acceptor stem is critical for recognition by Escherichia coli prolyl-tRNA synthetase (ProRS), but not for human ProRS. To further delineate species-specific differences in acceptor stem recognition, atomic group mutagenesis was used to probe the role of sugar-phosphate backbone interactions in recognition of human tRNA(Pro). Incorporation of site-specific 2'-deoxynucleotides, as well as phosphorothioate and methylphosphonate modifications within the tRNA acceptor stem revealed an extensive network of interactions with specific functional groups proximal to the first base pair and the discriminator base. Backbone functional groups located at the base of the acceptor stem, especially the 2'-hydroxyl of A66, are also critical for aminoacylation catalytic efficiency by human ProRS. Therefore, in contrast to the bacterial system, backbone-specific interactions contribute significantly more to tRNA recognition by the human enzyme than base-specific interactions. Taken together with previous studies, these data show that ProRS-tRNA acceptor stem interactions have co-adapted through evolution from a mechanism involving 'direct readout' of nucleotide bases to one relying primarily on backbone-specific 'indirect readout'.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.

Mechanisms of coupling between DNA recognition specificity and catalysis in EcoRI endonuclease

Proteins that bind to specific sites on DNA often do so in order to carry out catalysis or specific protein-protein interaction while bound to the recognition site. Functional specificity is enhanced if this second function is coupled to correct DNA site recognition. To analyze the structural and energetic basis of coupling between recognition and catalysis in EcoRI endonuclease, we have studied stereospecific phosphorothioate (PS) or methylphosphonate (PMe) substitutions at the scissile phosphate GpAATTC or at the adjacent phosphate GApATTC in combination with molecular-dynamics simulations of the catalytic center with bound Mg2+. The results show the roles in catalysis of individual phosphoryl oxygens and of DNA distortion and suggest that a "crosstalk ring" in the complex couples recognition to catalysis and couples the two catalytic sites to each other.

Modeling of flap endonuclease interactions with DNA substrate

Structure-specific 5' nucleases play an important role in DNA replication and repair uniquely recognizing an overlap flap DNA substrate and processing it into a DNA nick. However, in the absence of a high-resolution structure of the enzyme/DNA complex, the mechanism underlying this recognition and substrate specificity, which is key to the enzyme's function, remains unclear. Here, we propose a three-dimensional model of the structure-specific 5' flap endonuclease from Pyrococcus furiosus in its complex with DNA. The model is based on the known X-ray structure of the enzyme and a variety of biochemical and molecular dynamics (MD) data utilized in the form of distance restraints between the enzyme and the DNA. Contacts between the 5' flap endonuclease and the sugar-phosphate backbone of the overlap flap substrate were identified using enzyme activity assays on substrates with methylphosphonate or 2'-O-methyl substitutions. The enzyme footprint extends two to four base-pairs upstream and eight to nine base-pairs downstream of the cleavage site, thus covering 10-13 base-pairs of duplex DNA. The footprint data are consistent with a model in which the substrate is bound in the DNA-binding groove such that the downstream duplex interacts with the helix-hairpin-helix motif of the enzyme. MD simulations to identify the substrate orientation in this model are consistent with the results of the enzyme activity assays on the methylphosphonate and 2'-O-methyl-modified substrates. To further refine the model, 5' flap endonuclease variants with alanine point substitutions at amino acid residues expected to contact phosphates in the substrate and one deletion mutant were tested in enzyme activity assays on the methylphosphonate-modified substrates. Changes in the enzyme footprint observed for two point mutants, R64A and R94A, and for the deletion mutant in the enzyme's beta(A)/beta(B) region, were interpreted as being the result of specific interactions in the enzyme/DNA complex and were used as distance restraints in MD simulations. The final structure suggests that the substrate's 5' flap interacts with the enzyme's helical arch and that the helix-hairpin-helix motif interacts with the template strand in the downstream duplex eight base-pairs from the cleavage site. This model suggests specific interactions between the 3' end of the upstream oligonucleotide and the enzyme. The proposed structure presents the first detailed description of substrate recognition by structure-specific 5' nucleases.

Exocyclic groups in the minor groove influence the backbone conformation of DNA

Exocyclic groups in the minor groove of DNA modulate the affinity and positioning of nucleic acids to the histone protein. The addition of exocyclic groups decreases the formation of this protein-DNA complex, while their removal increases nucleosome formation. On the other hand, recent theoretical results show a strong correlation between the B(I)/B(II) phosphate backbone conformation and the hydration of the grooves of the DNA. We performed a simulation of the d(CGCGAATTCGCG)2 Drew Dickerson dodecamer and one simulation of the d(CGCIAATTCGCG)2 dodecamer in order to investigate the influence of the exocyclic amino group of guanine. The removal of the amino group introduces a higher intrinsic flexibility to DNA supporting the suggestions that make the enhanced flexibility responsible for the enlarged histone complexation affinity. This effect is attributed to changes in the destacking interactions of both strands of the DNA. The differences in the hydration of the minor groove could be the explanation of this flexibility. The changed hydration of the minor groove also leads to a different B(I)/B(II) substate pattern. Due to the fact that the histone preferentially builds contacts with the backbone of the DNA, we propose an influence of these B(I)/B(II) changes on the nucleosome formation process. Thus, we provide an additional explanation for the enhanced affinity to the histone due to removal of exocyclic groups. In terms of B(I)/B(II) we are also able to explain how minor groove binding ligands could affect the nucleosome assembly without disrupting the structure of DNA.

Studies of the Escherichia coli Trp repressor binding to its five operators and to variant operator sequences

The Escherichia coli Trp repressor binds to promoters of very different sequence and intrinsic activity. Its mode of binding to trp operator DNA has been studied extensively yet remains highly controversial. In order to examine the selectivity of the protein for DNA, we have used electromobility shift assays (EMSAs) to study its binding to synthetic DNA containing the core sequences of each of its five operators and of operator variants. Our results for DNA containing sequences of two of the operators, trpEDCBA and aroH are similar to those of previous studies. Up to three bands of lower mobility than the free DNA are obtained which are assigned to complexes of stoichiometry 1 : 1, 2 : 1 and 3 : 1 Trp repressor dimer to DNA. The mtr and aroL operators have not been studied previously in vitro. For DNA containing these sequences, we observe predominantly one retarded band in EMSA with mobility corresponding to 2 : 1 complexes. We have also obtained retardation of DNA containing the trpR operator sequence, which has only been previously obtained with super-repressor Trp mutants. This gives bands with mobilities corresponding to 1 : 1 and 2 : 1 complexes. In contrast, DNA containing containing a symmetrized trpR operator sequence, trpRs, gives a single retarded band with mobility corresponding solely to a 1 : 1 protein dimer-DNA complex. Using trpR operator variants, we show that a change in a single base pair in the core 20 base pairs can alter the number of retarded DNA bands in EMSA and the length of the DNase I footprint observed. This shows that the binding of the second dimer is sequence selective. We propose that the broad selectivity of Trp repressor coupled to tandem 2 : 1 binding, which we have observed with all five operator sequences, enables the Trp repressor to bind to a limited number of sites with diverse sequences. This allows it to co-ordinately control promoters of different intrinsic strength. This mechanism may be of importance in a number of promoters that bind multiple effector molecules.

Sequence specific hybridization properties of methylphosphonate oligodeoxynucleotides

Methylphosphonate oligodeoxynucleotides (MPO's) with isomerically pure Rp-configurated methylphosphonates (MP's) were synthesized by block coupling of ApT and TpA dinucleoside methylphosphonates (DMP's, p indicating MP-linkage). Oligonucleotide duplexes (20 mers) with these Rp-MP's showed almost the same melting temperatures (Tm) as those with phosphorodiester bonds. Further a dependence of the duplex stability from the nucleosides (bases) adjacent to the MP moiety was observed. For the first time thermodynamic parameters for the duplex to coil transition of isomerically pure MP's were determined from the concentration dependence of the Tm. CD-spectra of the duplexes show structural changes which can be associated with the transition to a compact helix with higher helix winding angle.