Local DNA Structure Research Articles

The importance of DNA sequence for nucleosome positioning in transcriptional regulation.

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Single-molecule analysis of DNA base-stacking energetics using patterned DNA nanostructures

The DNA double helix structure is stabilized by base-pairing and base-stacking interactions. However, a comprehensive understanding of dinucleotide base-stacking energetics is lacking. Here we combined multiplexed DNA-based point accumulation in nanoscale topography (DNA-PAINT) imaging with designer DNA nanostructures and measured the free energy of dinucleotide base stacking at the single-molecule level. Multiplexed imaging enabled us to extract the binding kinetics of an imager strand with and without additional dinucleotide stacking interactions. The DNA-PAINT data showed that a single additional dinucleotide base stacking results in up to 250-fold stabilization for the DNA duplex nanostructure. We found that the dinucleotide base-stacking energies vary from −0.95 ± 0.12 kcal mol−1 to −3.22 ± 0.04 kcal mol−1 for C|T and A|C base-stackings, respectively. We demonstrate the application of base-stacking energetics in designing DNA-PAINT probes for multiplexed super-resolution imaging, and efficient assembly of higher-order DNA nanostructures. Our results will aid in designing functional DNA nanostructures, and DNA and RNA aptamers, and facilitate better predictions of the local DNA structure.

SNP rs10420324 in the AMPA receptor auxiliary subunit TARP \u03b3-8 regulates the susceptibility to antisocial personality disorder

In the brain, AMPA receptors mediate fast excitatory neurotransmission, the dysfunction of which leads to neuropsychiatric disorders. Synaptic function of AMPA receptors is tightly controlled by a protein group called transmembrane AMPAR regulatory proteins (TARPs). TARP γ-8 (also known as CACNG8) preferentially expresses in the hippocampus, cortex and subcortical regions that are critical for emotion generation indicating its association with psychiatric disorders. Here, we identified rs10420324 (T/G), a SNP located in the human CACNG8 gene, regulated reporter gene expression in vitro and TARP γ-8 expression in the human brain. A guanine at the locus (rs10420324G) suppressed transcription likely through modulation of a local G-quadruplex DNA structure. Consistent with these observations, the frequency of rs10420324G was higher in patients with anti-social personality disorder (ASPD) than in controls, indicating that rs10420324G in CACNG8 is more voluntary for ASPD. We then characterized the behavior of TARP γ-8 knockout and heterozygous mice and found that consistent with ASPD patients who often exhibit impulsivity, aggression, risk taking, irresponsibility and callousness, a decreased γ-8 expression in mice displayed similar behaviors. Furthermore, we found that a decrease in TARP γ-8 expression impaired synaptic AMPAR functions in layer 2–3 pyramidal neurons of the prefrontal cortex, a brain region that inhibition leads to aggression, thus explaining, at least partially, the neuronal basis for the behavioral abnormality. Taken together, our study indicates that TARP γ-8 expression level is associated with ASPD, and that the TARP γ-8 knockout mouse is a valuable animal model for studying this psychiatric disease.

Removal of heat-sensitive clustered damaged DNA sites is independent of double-strand break repair.

DNA double-strand breaks (DSBs) are the most deleterious lesions that can arise in cells after ionizing radiation or radiometric drug treatment. In addition to prompt DSBs, DSBs may also be produced during repair, evolving from a clustered DNA damaged site, which is composed of two or more distinct lesions that are located within two helical turns. A specific type of cluster damage is the heat-sensitive clustered site (HSCS), which transforms into DSBs upon treatment at elevated temperatures. The actual lesions or mechanisms that mediate the HSCS transformation into DSBs are unknown. However, there are two possibilities; either these lesions are transformed into DSBs due to DNA lesion instability, e.g., transfer of HSCS into single-strand breaks (SSBs), or they are formed due to local DNA structure instability, e.g., DNA melting, where two SSBs on opposite strands meet and transform into a DSB. The importance of these processes in living cells is not understood, but they significantly affect estimates of DSB repair capacity. In this study, we show that HSCS removal in human cells is not affected by defects in DSB repair or inhibition of DSB repair. Under conditions where rejoining of prompt DSBs was almost completely inhibited, heat-sensitive DSBs were successfully rejoined, without resulting in increased DSB levels, indicating that HSCS do not transfer into DSB in cells under physiological conditions. Furthermore, analysis by atomic force microscopy suggests that prolonged heating of chromosomal DNA can induce structural changes that facilitate transformation of HSCS into DSB. In conclusion, the HSCS do not generate additional DSBs at physiological temperatures in human cells, and the repair of HSCS is independent of DSB repair.

Genome Instability Is Promoted by the Chromatin-Binding Protein Spn1 in Saccharomyces cerevisiae.

Cells expend a large amount of energy to maintain their DNA sequence. DNA repair pathways, cell cycle checkpoint activation, proofreading polymerases, and chromatin structure are ways in which the cell minimizes changes to the genome. During replication, the DNA-damage tolerance pathway allows the replication forks to bypass damage on the template strand. This avoids prolonged replication fork stalling, which can contribute to genome instability. The DNA-damage tolerance pathway includes two subpathways: translesion synthesis and template switch. Post-translational modification of PCNA and the histone tails, cell cycle phase, and local DNA structure have all been shown to influence subpathway choice. Chromatin architecture contributes to maintaining genome stability by providing physical protection of the DNA and by regulating DNA-processing pathways. As such, chromatin-binding factors have been implicated in maintaining genome stability. Using Saccharomyces cerevisiae, we examined the role of Spn1 (Suppresses postrecruitment gene number 1), a chromatin-binding and transcription elongation factor, in DNA-damage tolerance. Expression of a mutant allele of SPN1 results in increased resistance to the DNA-damaging agent methyl methanesulfonate, lower spontaneous and damage-induced mutation rates, along with increased chronological life span. We attribute these effects to an increased usage of the template switch branch of the DNA-damage tolerance pathway in the spn1 strain. This provides evidence for a role of wild-type Spn1 in promoting genome instability, as well as having ties to overcoming replication stress and contributing to chronological aging.

Systematic prediction of DNA shape changes due to CpG methylation explains epigenetic effects on protein\u2013DNA binding

BackgroundDNA shape analysis has demonstrated the potential to reveal structure-based mechanisms of protein–DNA binding. However, information about the influence of chemical modification of DNA is limited. Cytosine methylation, the most frequent modification, represents the addition of a methyl group at the major groove edge of the cytosine base. In mammalian genomes, cytosine methylation most frequently occurs at CpG dinucleotides. In addition to changing the chemical signature of C/G base pairs, cytosine methylation can affect DNA structure. Since the original discovery of DNA methylation, major efforts have been made to understand its effect from a sequence perspective. Compared to unmethylated DNA, however, little structural information is available for methylated DNA, due to the limited number of experimentally determined structures. To achieve a better mechanistic understanding of the effect of CpG methylation on local DNA structure, we developed a high-throughput method, methyl-DNAshape, for predicting the effect of cytosine methylation on DNA shape.ResultsUsing our new method, we found that CpG methylation significantly altered local DNA shape. Four DNA shape features—helix twist, minor groove width, propeller twist, and roll—were considered in this analysis. Distinct distributions of effect size were observed for different features. Roll and propeller twist were the DNA shape features most strongly affected by CpG methylation with an effect size depending on the local sequence context. Methylation-induced changes in DNA shape were predictive of the measured rate of cleavage by DNase I and suggest a possible mechanism for some of the methylation sensitivities that were recently observed for human Pbx-Hox complexes.ConclusionsCpG methylation is an important epigenetic mark in the mammalian genome. Understanding its role in protein–DNA recognition can further our knowledge of gene regulation. Our high-throughput methyl-DNAshape method can be used to predict the effect of cytosine methylation on DNA shape and its subsequent influence on protein–DNA interactions. This approach overcomes the limited availability of experimental DNA structures that contain 5-methylcytosine.

Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism.

Torsional restraints on DNA change in time and space during the life of the cell and are an integral part of processes such as gene expression, DNA repair and packaging. The mechanical behavior of DNA under torsional stress has been studied on a mesoscopic scale, but little is known concerning its response at the level of individual base pairs and the effects of base pair composition. To answer this question, we have developed a geometrical restraint that can accurately control the total twist of a DNA segment during all-atom molecular dynamics simulations. By applying this restraint to four different DNA oligomers, we are able to show that DNA responds to both under- and overtwisting in a very heterogeneous manner. Certain base pair steps, in specific sequence environments, are able to absorb most of the torsional stress, leaving other steps close to their relaxed conformation. This heterogeneity also affects the local torsional modulus of DNA. These findings suggest that modifying torsional stress on DNA could act as a modulator for protein binding via the heterogeneous changes in local DNA structure.

Precise mapping of 17 deletion breakpoints within the central hotspot deletion region (introns 50 and 51) of the DMD gene.

Exon deletions in the human DMD gene, which encodes the dystrophin protein, are the molecular defect in 50-70% of cases of Duchenne/Becker muscular dystrophies. Deletions are preferentially clustered in the 5' (exons 2-20) and the central (exons 45-53) region of DMD, likely because local DNA structure predisposes to specific breakage or recombination events. Notably, innovative therapeutic strategies may rescue dystrophin function by homology-based specific targeting of sequences within the central DMD hot spot deletion region. To further study molecular mechanisms that generate such frequent genome variations and to identify residual intronic sequences, we sequenced 17 deletion breakpoints within introns 50 and 51 of DMD and analyzed the surrounding genomic architecture. There was no breakpoint clustering within the introns nor extensive homology between sequences adjacent to each junction. However, at or near the breakpoint, we found microhomology, short tandem repeats, interspersed repeat elements and short sequence stretches that predispose to DNA deletion or bending. Identification of such structural elements contributes to elucidate general mechanisms generating deletion within the DMD gene. Moreover, precise mapping of deletion breakpoints and localization of repeated elements are of interest, because residual intronic sequences may be targeted by therapeutic strategies based on genome editing correction.

TRANSCRIPTION. Allosteric transcriptional regulation via changes in the overall topology of the core promoter.

Many transcriptional activators act at a distance from core promoter elements and work by recruiting RNA polymerase through protein-protein interactions. We show here how the prokaryotic regulatory protein CueR both represses and activates transcription by differentially modulating local DNA structure within the promoter. Structural studies reveal that the repressor state slightly bends the promoter DNA, precluding optimal RNA polymerase-promoter recognition. Upon binding a metal ion in the allosteric site, CueR switches into an activator conformation. It maintains all protein-DNA contacts but introduces torsional stresses that kink and undertwist the promoter, stabilizing an A-form DNA-like conformation. These factors switch on and off transcription by exerting dynamic control of DNA stereochemistry, reshaping the core promoter and making it a better or worse substrate for polymerase.

A DNA shape-based regulatory score improves position-weight matrix-based recognition of transcription factor binding sites.

The position-weight matrix (PWM) is a useful representation of a transcription factor binding site (TFBS) sequence pattern because the PWM can be estimated from a small number of representative TFBS sequences. However, because the PWM probability model assumes independence between individual nucleotide positions, the PWMs for some TFs poorly discriminate binding sites from non-binding-sites that have similar sequence content. Since the local three-dimensional DNA structure ('shape') is a determinant of TF binding specificity and since DNA shape has a significant sequence-dependence, we combined DNA shape-derived features into a TF-generalized regulatory score and tested whether the score could improve PWM-based discrimination of TFBS from non-binding-sites. We compared a traditional PWM model to a model that combines the PWM with a DNA shape feature-based regulatory potential score, for accuracy in detecting binding sites for 75 vertebrate transcription factors. The PWM+shape model was more accurate than the PWM-only model, for 45% of TFs tested, with no significant loss of accuracy for the remaining TFs. The shape-based model is available as an open-source R package at that is archived on the GitHub software repository at https://github.com/ramseylab/regshape/. stephen.ramsey@oregonstate.edu Supplementary data are available at Bioinformatics online.

Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes.

Single base substitutions (SBSs) and insertions/deletions are critical for generating population diversity and can lead both to inherited disease and cancer. Whereas on a genome-wide scale SBSs are influenced by cellular factors, on a fine scale SBSs are influenced by the local DNA sequence-context, although the role of flanking sequence is often unclear. Herein, we used bioinformatics, molecular dynamics and hybrid quantum mechanics/molecular mechanics to analyze sequence context-dependent mutagenesis at mononucleotide repeats (A-tracts and G-tracts) in human population variation and in cancer genomes. SBSs and insertions/deletions occur predominantly at the first and last base-pairs of A-tracts, whereas they are concentrated at the second and third base-pairs in G-tracts. These positions correspond to the most flexible sites along A-tracts, and to sites where a ‘hole’, generated by the loss of an electron through oxidation, is most likely to be localized in G-tracts. For A-tracts, most SBSs occur in the direction of the base-pair flanking the tracts. We conclude that intrinsic features of local DNA structure, i.e. base-pair flexibility and charge transfer, render specific nucleotides along mononucleotide runs susceptible to base modification, which then yields mutations. Thus, local DNA dynamics contributes to phenotypic variation and disease in the human population.

A repetitive DNA element regulates expression of the Helicobacter pylori sialic acid binding adhesin by a rheostat-like mechanism.

During persistent infection, optimal expression of bacterial factors is required to match the ever-changing host environment. The gastric pathogen Helicobacter pylori has a large set of simple sequence repeats (SSR), which constitute contingency loci. Through a slipped strand mispairing mechanism, the SSRs generate heterogeneous populations that facilitate adaptation. Here, we present a model that explains, in molecular terms, how an intergenically located T-tract, via slipped strand mispairing, operates with a rheostat-like function, to fine-tune activity of the promoter that drives expression of the sialic acid binding adhesin, SabA. Using T-tract variants, in an isogenic strain background, we show that the length of the T-tract generates multiphasic output from the sabA promoter. Consequently, this alters the H. pylori binding to sialyl-Lewis x receptors on gastric mucosa. Fragment length analysis of post-infection isolated clones shows that the T-tract length is a highly variable feature in H. pylori. This mirrors the host-pathogen interplay, where the bacterium generates a set of clones from which the best-fit phenotypes are selected in the host. In silico and functional in vitro analyzes revealed that the length of the T-tract affects the local DNA structure and thereby binding of the RNA polymerase, through shifting of the axial alignment between the core promoter and UP-like elements. We identified additional genes in H. pylori, with T- or A-tracts positioned similar to that of sabA, and show that variations in the tract length likewise acted as rheostats to modulate cognate promoter output. Thus, we propose that this generally applicable mechanism, mediated by promoter-proximal SSRs, provides an alternative mechanism for transcriptional regulation in bacteria, such as H. pylori, which possesses a limited repertoire of classical trans-acting regulatory factors.

Preferential binding of p53 tumor suppressor to p21 promoter sites that contain inverted repeats capable of forming cruciform structure

p53 Is one of the most critical proteins involved in protecting organisms from malignancies and its gene is frequently mutated in these diseases. p53 Functions as a transcription factor and its role in the cell is mediated by sequence-specific DNA binding. Although the genome contains many p53-binding sequences, the p53 protein binds only a subset of these sequences with high affinity. One likely mechanism of how p53 binds DNA effectively underlies its ability to recognize selective local DNA structure. We analyzed the possibility of cruciform structure formation within different regions of the p21 gene promoter. p53 protein remarkably activates the transcription of p21 gene after genotoxic treatment. In silico analysis showed that p21 gene promoter contains numerous p53 target sequences, some of which have inverted repeats capable of forming cruciform structures. Using chromatin immunoprecipitation, we demonstrated that p53 protein binds preferentially to sequences that not only contain inverted repeats but also have the ability to create local cruciform structures. Gel retardation assay also revealed strong preference of the p53 protein for response element in superhelical state, with cruciform structure in the DNA sequence. Taken together, our results suggest that p53 response element’s potential for cruciform structure formation could be an additional determinant in p53 DNA-binding machinery.

Tandem multiplication of the IS26-flanked amplicon with theblaSHV-5gene within plasmid p1658/97

The IncF plasmid p1658/97 (c. 125 kb) from Escherichia coli isolates recovered during a clonal outbreak in a hospital in Warsaw, Poland, in 1997 contains the extended-spectrum β-lactamase (ESBL) gene bla(SHV-5), originated from the Klebsiella pneumoniae chromosome. A region containing the bla(SHV-5) gene is flanked by two IS26 copies and its copy number multiplies spontaneously within p1658/97 and RecA-deficient E. coli strains. Here, we demonstrate that the amplified IS26-bla(SHV-5) units were arranged in tandems, containing up to more than 10 units, which could raise ceftazidime MICs for host strains from 4 μg mL(-1) to more than 128 μg mL(-1). Successive deletions within p1658/97, located outside the amplifiable module and encompassing even as little as c. 15% of the plasmid, blocked the amplification. Moreover, the complementing re-introduction of the deleted fragments in trans did not restore the process. Similarly, insertions of a 1-kb DNA fragment into the amplicon inhibited its self-multiplication ability. The module was able to transmit into another IS26-containing plasmid by recombination. The results prompted us to speculate that local DNA structure, especially favorable in p1658/97, might have been responsible for the IS26-bla(SHV-5) multiplication ability.

Functional Implications of Local DNA Structures in Regulatory Motifs

The three-dimensional structure of DNA has been proposed to be a major determinant for functional transcription factors (TFs) and DNA interaction. Here, we use hydroxyl radical cleavage pattern as a measure of local DNA structure. We compared the conservation between DNA sequence and structure in terms of information content and attempted to assess the functional implications of DNA structures in regulatory motifs. We used statistical methods to evaluate the structural divergence of substituting a single position within a binding site and applied them to a collection of putative regulatory motifs. The following are our major observations: (i) we observed more information in structural alignment than in the corresponding sequence alignment for most of the transcriptional factors; (ii) for each TF, majority of positions have more information in the structural alignment as compared to the sequence alignment; (iii) we further defined a DNA structural divergence score (SD score) for each wild-type and mutant pair that is distinguished by single-base mutation. The SD score for benign mutations is significantly lower than that of switch mutations. This indicates structural conservation is also important for TFBS to be functional and DNA structures will provide previously unappreciated information for TF to realize the binding specificity.

Temporal behavior of DNA thermal stability in the presence of platinum compounds. Role of monofunctional and bifunctional adducts

Penetrating into cell nuclei, antitumor drug cisplatin sequentially forms various intermediate and final adducts destroying local DNA structure. The demonstrated disappearance of the fine structure of melting curve of long DNAs along with a strong decrease in melting enthalpy conforms to the structural impact. However, the negative thermal effect (δTm) caused by cisplatin is relatively small if neutral medium is used in melting experiments. Cisplatin's inactive analogs transplatin and diethylenetriaminechloroplatinum {Pt[(dien)Cl]Cl} also distort DNA structure but their thermal effect is even positive. We have found that the use of alkaline medium in melting experiments strengthens the negative thermal effect for cisplatin. For transplatin and Pt[(dien)Cl]Cl, the thermal effect becomes negative that makes it qualitatively consistent with structural distortions. Those changes are explained by elimination of nonspecific electrostatic stabilization of DNA under platination. Additionally, alkaline medium fixes intermediate states of DNA platination and makes them stable against heating. These results allowed us to monitor δTm under binding of platinum compounds to DNA and their further transformation. The kinetic and thermal characteristics of monofunctional and bifunctional adducts were evaluated. It has been demonstrated that monofunctional adducts of cisplatin, transplatin and Pt[(dien)Cl]Cl produce approximately the same thermal destabilization. Cisplatin intrastrand crosslinks cause a two-fold stronger thermal destabilization than its monofunctional adducts. The value of δTm for cisplatin's final adducts is ten times larger than for transplatin. This difference mainly comes from the much stronger thermal destabilizing power of cisplatin's intrastrand crosslinks, which are responsible for antitumor activity of this compound.

Enzymatic Excision of Uracil Residues in Nucleosomes Depends on the Local DNA Structure and Dynamics

The excision of uracil bases from DNA is accomplished by the enzyme uracil DNA glycosylase (UNG). Recognition of uracil bases in free DNA is facilitated by uracil base pair dynamics, but it is not known whether this same mechanistic feature is relevant for detection and excision of uracil residues embedded in nucleosomes. Here we investigate this question using nucleosome core particles (NCPs) generated from Xenopus laevis histones and the high-affinity "Widom 601" positioning sequence. The reactivity of uracil residues in NCPs under steady-state multiple-turnover conditions was generally decreased compared to that of free 601 DNA, mostly because of anticipated steric effects of histones. However, some sites in NCPs had equal or even greater reactivity than free DNA, and the observed reactivities were not readily explained by simple steric considerations or by global DNA unwrapping models for nucleosome invasion. In particular, some reactive uracils were found in occluded positions, while some unreactive uracils were found in exposed positions. One feature of many exposed reactive sites is a wide DNA minor groove, which allows penetration of a key active site loop of the enzyme. In single-turnover kinetic measurements, multiphasic reaction kinetics were observed for several uracil sites, where each kinetic transient was independent of the UNG concentration. These kinetic measurements, and supporting structural analyses, support a mechanism in which some uracils are transiently exposed to UNG by local, rate-limiting nucleosome conformational dynamics, followed by rapid trapping of the exposed state by the enzyme. We present structural models and plausible reaction mechanisms for the reaction of UNG at three distinct uracil sites in the NCP.

Radioprobing the conformation of DNA in a p53-DNA complex

Purpose: The frequency of DNA strand breaks produced by the decay of Auger electron-emitting radionuclides is inversely proportional to the distance of DNA nucleotides from the decay site; and thus is very sensitive to changes in the local conformation of the DNA. Analysis of the frequency of DNA breaks, or radioprobing, gives valuable information about the local DNA structure. More than 10 years ago, we demonstrated the feasibility of radioprobing using a DNA-repressor complex with a known structure. Herein, we used radioprobing to study the conformation of DNA in complex with the tumor suppressor protein 53 (p53). Several structures of p53-DNA complexes have been solved by X-ray crystallography. These structures, obtained with the p53 DNA binding domain, a truncated form, laid the groundwork for understanding p53-DNA interactions and their relation to p53 functions. However, whether all observed stereochemical details are relevant to the native p53-DNA complex remains unclear. A common theme of the crystallographic structures is the lack of significant bending in the central part of the DNA response element. In contrast, gel electrophoresis and electron microscopy data showed strong DNA bending and overtwisting upon binding to the native p53 tetramer.Methods: To analyze DNA in complex with p53, we incorporated 125I-dCTP in two different positions of synthetic duplexes containing the consensus p53-binding site.Results: The most significant changes in the break frequency distributions were detected close to the center of the binding site, which is consistent with an increase in DNA twisting in this region and local DNA bending and sliding.Conclusions: Our data confirm the main results of the studies made in solution and lay a foundation for systematic examination of interactions between DNA and native p53 using 125I radioprobing.

A porphycene-DNA hybrid and its DNA-templated interactions with a porphyrin

A porphycene with a hydroxyethyl side chain was coupled to the 3′-terminus of an oligodeoxynucleotide via solid-phase synthesis. The resulting porphycene-DNA hybrid binds to a complementary region of a DNA target strand with greater affinity than the unmodified control oligonucleotide, resulting in an increase of the UV-melting point of 12.7 °C. Duplex formation is accompanied by an increase in porphycene fluorescence at 640 nm by 18%. When a tetrakis(p-hydroxyphenyl)porphyrin appended to the 5′-terminus of another DNA-strand is brought into close proximity of the porphycene by hybridizing it to the downstream-region of the template strand, 94% of the porphycene fluorescence is quenched. By quenching each others fluorescence to different degrees, porphycene and porphyrin, together, report on local DNA structure in a fashion reminiscent of that of molecular beacons. The introduction of porphycenes and the porphycene-porphyrin "two hybrid system'' to DNA-based structuring may open up new avenues to designed functional materials.

Distribution of MGEs and their insertion sites in the Macaca mulatta genome.

Mobile genetic elements (MGEs) are fragments of DNA that can move around within the genome through retrotransposition. These are responsible for various important events such as gene inactivation, transduction, regulation of gene expression and genome expansion. The present work involves the identification and study of the distribution of Alu and L1 retrotransposons in the genome of Macaca mulatta, an extensively used organism in biomedical studies. We also make comparisons with MGE distributions in other primate genomes and study the physicochemical properties of the local DNA structure around the transposon insertion site using ELAN. The present work also includes computational testing of the pre-insertion loci in order to detect unique features based on DNA structure, thermodynamic considerations and protein interaction measures. Although there is significant sequence divergence between the elements of M. mulatta and H. sapiens, their genome wide distribution is very similar; comparing the distribution of L1’s in all available X chromosome sequences suggests a common mechanism behind the spread of MGE’s in primate genomes.