Crossover isomer bias is the primary sequence-dependent property of immobilized Holliday junctions

Recombination of genes is essential to the evolution of genetic diversity, the segregation of chromosomes during cell division, and certain DNA repair processes. The Holliday junction, a four-arm, four-strand branched DNA crossover structure, is formed as a transient intermediate during genetic recombination and repair processes in the cell. The recognition and subsequent resolution of Holliday junctions into parental or recombined products appear to be critically dependent on their three-dimensional structure. Complementary NMR and time-resolved fluorescence resonance energy transfer experiments on immobilized four-arm DNA junctions reported here indicate that the Holliday junction cannot be viewed as a static structure but rather as an equilibrium mixture of two conformational isomers. Furthermore, the distribution between the two possible crossover isomers was found to depend on the sequence in a manner that was not anticipated on the basis of previous low-resolution experiments.

Sequence dependent hypermutation of the immunoglobulin heavy chain in cultured B cells

The variable (V) regions of immunoglobulin heavy and light chains undergo high rates of somatic mutation during the immune response. Although point mutations accumulate throughout the V regions and their immediate flanking sequences, analysis of large numbers of mutations that have arisen in vivo reveal that the triplet AGC appears to be most susceptible to mutation. We have stably transfected B cell lines with γ2a heavy chain constructs containing TAG nonsense codons in their V regions that are part of either a putative (T)AGC hot spot or a (T)AGA non-hot spot motif. Using an ELISA spot assay to detect revertants and fluctuation analysis to determine rates of mutation, the rate of reversion of the TAG nonsense codon has been determined for different motifs in different parts of the V region. In the NSO plasma cell line, the (T)AGC hot spot motif mutates at rates of ≈6 × 10−4/bp per generation and ≈3 × 10−5/bp per generation at residues 38 and 94 in the V region. At each of these locations, the (T)AGC hot spot motif is 20–30 times more likely to undergo mutation than the (T)AGA non-hot spot motif. Moreover, the AGA non-hot spot motif mutates at as high a rate as the hot spot motif when it is located adjacent to hot spot motifs, suggesting that more extended sequences influence susceptibility to mutation.

Transplantation of a 17-amino acid alpha-helical DNA-binding domain into an antibody molecule confers sequence-dependent DNA recognition.

Recombinant antibodies capable of sequence-specific interactions with nucleic acids represent a class of DNA- and RNA-binding proteins with potential for broad application in basic research and medicine. We describe the rational design of a DNA-binding antibody, Fab-Ebox, by replacing a variable segment of the immunoglobulin heavy chain with a 17-amino acid domain derived from TFEB, a class B basic helix-loop-helix protein. DNA-binding activity was studied by electrophoretic mobility-shift assays in which Fab-Ebox was shown to form a specific complex with DNA containing the TFEB recognition motif (CACGTG). Similarities were found in the abilities of TFEB and Fab-Ebox to discriminate between oligodeoxyribonucleotides containing altered recognition sequences. Comparable interference of binding by methylation of cytosine residues indicated that Fab-Ebox and TFEB both contact DNA through interactions along the major groove of double-stranded DNA. The results of this study indicate that DNA-binding antibodies of high specificity can be developed by using the modular nature of both immunoglobulins and transcription factors.

The Role of the Membrane-spanning Domain Sequence in Glycoprotein-mediated Membrane Fusion

The role of glycoprotein membrane-spanning domains in the process of membrane fusion is poorly understood. It has been demonstrated that replacing all or part of the membrane-spanning domain of a viral fusion protein with sequences that encode signals for glycosylphosphatidylinositol linkage attachment abrogates membrane fusion activity. It has been suggested, however, that the actual amino acid sequence of the membrane-spanning domain is not critical for the activity of viral fusion proteins. We have examined the function of Moloney murine leukemia virus envelope proteins with substitutions in the membrane-spanning domain. Envelope proteins bearing substitutions for proline 617 are processed and incorporated into virus particles normally and bind to the viral receptor. However, they possess greatly reduced or undetectable capacities for the promotion of membrane fusion and infectious virus particle formation. Our results imply a direct role for the residues in the membrane-spanning domain of the murine leukemia virus envelope protein in membrane fusion and its regulation. They also support the thesis that membrane-spanning domains possess a sequence-dependent function in other protein-mediated membrane fusion events.

The effects of sequence context on DNA curvature.

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.

Architectural specificity in chromatin structure at the TATA box in vivo: Nucleosome displacement upon β-phaseolin gene activation

Extensive studies of the β-phaseolin (phas) gene in transgenic tobacco have shown that it is highly active during seed embryogenesis but is completely silent in leaf and other vegetative tissues. In vivo footprinting revealed that the lack of even basal transcriptional activity in vegetative tissues is associated with the presence of a nucleosome that is rotationally positioned with base pair precision over three phased TATA boxes present in the phas promoter. Positioning is sequence-dependent because an identical rotational setting is obtained upon nucleosome reconstitution in vitro. A comparison of DNase I and dimethyl sulfate footprints in vivo and in vitro strongly suggests that this repressive chromatin architecture is remodeled concomitant with gene activation in the developing seed. This leads to the disruption of histone-mediated DNA wrapping and the assembly of the TATA boxes into a transcriptionally competent nucleoprotein complex.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Antitumor activity and pharmacokinetics of a mixed-backbone antisense oligonucleotide targeted to the RIα subunit of protein kinase A after oral administration

Overexpression of the RIα subunit of cAMP-dependent protein kinase (PKA) has been demonstrated in various human cancers. PKA has been suggested as a potential target for cancer therapy. The goal of the present study was to evaluate an anti-PKA antisense oligonucleotide (mixed-backbone oligonucleotide) as a therapeutic approach to human cancer treatment. The identified oligonucleotide inhibited the growth of cell lines of human colon cancer (LS174T, DLD-1), leukemia (HL-60), breast cancer (MCF-7, MDA-MB-468), and lung cancer (A549) in a time-, concentration-, and sequence-dependent manner. In a dose-dependent manner, the oligonucleotide displayed in vivo antitumor activity in severe combined immunodeficient and nude mice bearing xenografts of human cancers of the colon (LS174T), breast (MDA-MB-468), and lung (A549). The routes of drug administration were intraperitoneal and oral. Synergistic effects were found when the antisense oligonucleotide was used in combination with the cancer chemotherapeutic agent cisplatin. The pharmacokinetics of the oligonucleotide after oral administration of 35S-labeled oligonucleotide into tumor-bearing mice indicated an accumulation and retention of the oligonucleotide in tumor tissue. This study further provides a basis for clinical studies of the antisense oligonucleotide targeted to the RIα subunit of PKA (GEM 231) as a cancer therapeutic agent used alone or in combination with conventional chemotherapy.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Mapping the intrinsic curvature and flexibility along the DNA chain

The energy of DNA deformation plays a crucial and active role in its packaging and its function in the cell. Considerable effort has gone into developing methodologies capable of evaluating the local sequence-directed curvature and flexibility of a DNA chain. These studies thus far have focused on DNA constructs expressly tailored either with anomalous flexibility or curvature tracts. Here we demonstrate that these two structural properties can be mapped also along the chain of a “natural” DNA with any sequence on the basis of its scanning force microscope (SFM) images. To know the orientation of the sequence of the investigated DNA molecules in their SFM images, we prepared a palindromic dimer of the long DNA molecule under study. The palindromic symmetry also acted as an internal gauge of the statistical significance of the analysis carried out on the SFM images of the dimer molecules. It was found that although the curvature modulus is not efficient in separating static and dynamic contributions to the curvature of the population of molecules, the curvature taken with its direction (its sign in two dimensions) permits the direct separation of the intrinsic curvature from the flexibility contributions. The sequence-dependent flexibility seems to vary monotonically with the chain's intrinsic curvature; the chain rigidity was found to modulate as its local thermodynamic stability and does not correlate with the dinucleotide chain rigidities evaluation made from x-ray data by other authors.

Identification of the yeast cytidine deaminase CDD1 as an orphan C→U RNA editase

Yeast co-expressing rat APOBEC-1 and a fragment of human apolipoprotein B (apoB) mRNA assembled functional editosomes and deaminated C6666 to U in a mooring sequence-dependent fashion. The occurrence of APOBEC-1-complementing proteins suggested a naturally occurring mRNA editing mechanism in yeast. Previously, a hidden Markov model identified seven yeast genes encoding proteins possessing putative zinc-dependent deaminase motifs. Here, only CDD1, a cytidine deaminase, is shown to have the capacity to carry out C→U editing on a reporter mRNA. This is only the second report of a cytidine deaminase that can use mRNA as a substrate. CDD1-dependent editing was growth phase regulated and demonstrated mooring sequence-dependent editing activity. Candidate yeast mRNA substrates were identified based on their homology with the mooring sequence-containing tripartite motif at the editing site of apoB mRNA and their ability to be edited by ectopically expressed APOBEC-1. Naturally occurring yeast mRNAs edited to a significant extent by CDD1 were, however, not detected. We propose that CDD1 be designated an orphan C→U editase until its native RNA substrate, if any, can be identified and that it be added to the CDAR (cytidine deaminase acting on RNA) family of editing enzymes.

Prediction of protein deamidation rates from primary and three-dimensional structure

A method for the quantitative estimation of instability with respect to deamidation of the asparaginyl (Asn) residues in proteins is described. The procedure involves the observation of several simple aspects of the three-dimensional environment of each Asn residue in the protein and a calculation that includes these observations, the primary amino acid residue sequence, and the previously reported complete set of sequence-dependent rates of deamidation for Asn pentapeptides. This method is demonstrated and evaluated for 23 proteins in which 31 unstable and 167 stable Asn residues have been reported and for 7 unstable and 63 stable Asn residues that have been reported in 61 human hemoglobin variants. The relative importance of primary structure and three-dimensional structure in Asn deamidation is estimated.