Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.

Investigations of DNA-mediated protein oxidation

DNA charge transport (CT) involves the efficient transfer of electrons or electron holes through the DNA π-stack over long molecular distances of at least 100 base-pairs. Despite this shallow distance dependence, DNA CT is sensitive to mismatches or lesions that disrupt π-stacking and is critically dependent on proper electronic coupling of the donor and acceptor moieties into the base stack. Favorable DNA CT is very rapid, occurring on the picosecond timescale. Because of this speed, electron holes equilibrate along the DNA π-stack, forming a characteristic pattern of DNA damage at low oxidation potential guanine multiplets. Furthermore, DNA CT may be used in a biological context. DNA processing enzymes with 4Fe4S clusters can perform DNA-mediated electron transfer (ET) self-exchange reactions with other 4Fe4S cluster proteins, even if the proteins are quite dissimilar, as long as the DNA-bound [4Fe4S]^3+/2+ redox potentials are conserved. This mechanism would allow low copy number DNA repair proteins to find their lesions efficiently within the cell. DNA CT may also be used biologically for the long-range, selective activation of redox-active transcription factors. Within this work, we pursue other proteins that may utilize DNA CT within the cell and further elucidate aspects of the DNA-mediated ET self-exchange reaction of 4Fe4S cluster proteins.

Dps proteins, bacterial mini-ferritins that protect DNA from oxidative stress, are implicated in the survival and virulence of pathogenic bacteria. One aspect of their protection involves ferroxidase activity, whereby ferrous iron is bound and oxidized selectively by hydrogen peroxide, thereby preventing formation of damaging hydroxyl radicals via Fenton chemistry. Understanding the specific mechanism by which Dps proteins protect the bacterial genome could inform the development of new antibiotics. We investigate whether DNA-binding E. coli Dps can utilize DNA CT to protect the genome from a distance. An intercalating ruthenium photooxidant was employed to generate oxidative DNA damage via the flash-quench technique, which localizes to a low potential guanine triplet. We find that Dps loaded with ferrous iron, in contrast to Apo-Dps and ferric iron-loaded Dps which lack available reducing equivalents, significantly attenuates the yield of oxidative DNA damage at the guanine triplet. These data demonstrate that ferrous iron-loaded Dps is selectively oxidized to fill guanine radical holes, thereby restoring the integrity of the DNA. Luminescence studies indicate no direct interaction between the ruthenium photooxidant and Dps, supporting the DNA-mediated oxidation of ferrous iron-loaded Dps. Thus DNA CT may be a mechanism by which Dps efficiently protects the genome of pathogenic bacteria from a distance.

Further work focused on spectroscopic characterization of the DNA-mediated oxidation of ferrous iron-loaded Dps. X-band EPR was used to monitor the oxidation of DNA-bound Dps after DNA photooxidation via the flash-quench technique. Upon irradiation with poly(dGdC)₂, a signal arises with g = 4.3, consistent with the formation of mononuclear high-spin Fe(III) sites of low symmetry, the expected oxidation product of Dps with one iron bound at each ferroxidase site. When poly(dGdC)₂ is substituted with poly(dAdT)₂, the yield of Dps oxidation is decreased significantly, indicating that guanine radicals facilitate Dps oxidation. The more favorable oxidation of Dps by guanine radicals supports the feasibility of a long-distance protection mechanism via DNA CT where Dps is oxidized to fill guanine radical holes in the bacterial genome produced by reactive oxygen species.

We have also explored possible electron transfer intermediates in the DNA-mediated oxidation of ferrous iron-loaded Dps. Dps proteins contain a conserved tryptophan residue in close proximity to the ferroxidase site (W52 in E. coli Dps). In comparison to WT Dps, in EPR studies of the oxidation of ferrous iron-loaded Dps following DNA photooxidation, W52Y and W52A mutants were deficient in forming the characteristic EPR signal at g = 4.3, with a larger deficiency for W52A compared to W52Y. In addition to EPR, we also probed the role of W52 Dps in cells using a hydrogen peroxide survival assay. Bacteria containing W52Y Dps survived the hydrogen peroxide challenge more similarly to those containing WT Dps, whereas cells with W52A Dps died off as quickly as cells without Dps. Overall, these results suggest the possibility of W52 as a CT hopping intermediate.

DNA-modified electrodes have become an essential tool for the study of the redox chemistry of DNA processing enzymes with 4Fe4S clusters. In many cases, it is necessary to investigate different complex samples and substrates in parallel in order to elucidate this chemistry. Therefore, we optimized and characterized a multiplexed electrochemical platform with the 4Fe4S cluster base excision repair glycosylase Endonuclease III (EndoIII). Closely packed DNA films, where the protein has limited surface accessibility, produce EndoIII electrochemical signals sensitive to an intervening mismatch, indicating a DNA-mediated process. Multiplexed analysis allowed more robust characterization of the CT-deficient Y82A EndoIII mutant, as well as comparison of a new family of mutations altering the electrostatics surrounding the 4Fe4S cluster in an effort to shift the reduction potential of the cluster. While little change in the DNA-bound midpoint potential was found for this family of mutants, likely indicating the dominant effect of DNA-binding on establishing the protein redox potential, significant variations in the efficiency of DNA-mediated electron transfer were apparent. On the basis of the stability of these proteins, examined by circular dichroism, we proposed that the electron transfer pathway in EndoIII can be perturbed not only by the removal of aromatic residues but also through changes in solvation near the cluster.

While the 4Fe4S cluster of EndoIII is relatively insensitive to oxidation and reduction in solution, we have found that upon DNA binding, the reduction potential of the [4Fe4S]^3+/2+ couple shifts negatively by approximately 200 mV, bringing this couple into a physiologically relevant range. Demonstrated using electrochemistry experiments in the presence and absence of DNA, these studies do not provide direct molecular evidence for the species being observed. Sulfur K-edge X-ray absorbance spectroscopy (XAS) can be used to probe directly the covalency of iron-sulfur clusters, which is correlated to their reduction potential. We have shown that the Fe-S covalency of the 4Fe4S cluster of EndoIII increases upon DNA binding, stabilizing the oxidized [4Fe4S]³⁺ cluster, consistent with a negative shift in reduction potential. The 7% increase in Fe-S covalency corresponds to an approximately 150 mV shift, remarkably similar to DNA electrochemistry results. Therefore we have obtained direct molecular evidence for the shift in 4Fe4S reduction potential of EndoIII upon DNA binding, supporting the feasibility of our model whereby these proteins can utilize DNA CT to cooperate in order to efficiently find DNA lesions inside cells.

In conclusion, in this work we have explored the biological applications of DNA CT. We discovered that the DNA-binding bacterial ferritin Dps can protect the bacterial genome from a distance via DNA CT, perhaps contributing to pathogen survival and virulence. Furthermore, we optimized a multiplexed electrochemical platform for the study of the redox chemistry of DNA-bound 4Fe4S cluster proteins. Finally, we have used sulfur K-edge XAS to obtain direct molecular evidence for the negative shift in 4Fe4S cluster reduction potential of EndoIII upon DNA binding. These studies contribute to the understanding of DNA-mediated protein oxidation within cells.

Investigations of Pyrrole-Imidazole Polyamide Effects on DNA Replication

Pyrrole–Imidazole polyamides are programmable, cell-permeable small molecules that bind in the minor groove of double-stranded DNA sequence-specifically. Polyamide binding has been shown to alter the local helical structure of DNA, disrupt protein-DNA interactions, and modulate endogenous gene expression. Py–Im polyamides targeted to the androgen receptor-DNA interface have been observed to decrease expression of androgen-regulated genes, upregulate p53, and induce apoptosis in a hormone-sensitive prostate cancer cell line. Here we report that androgen response element (ARE)-targeted polyamides induced DNA replication stress in a hormone-insensitive prostate cancer cell line. The ATR checkpoint kinase was activated in response to this stress, causing phosphorylation of MCM2, and FANCD2 was monoubiquitinated. Surprisingly, little single-stranded DNA was exhibited, and the ATR targets RPA2 and Chk1 were not phosphorylated. We conclude that polyamide induces relatively low level replication stress, and suggest inhibition of the replicative helicase as a putative mechanism based on in vitro assays. We also demonstrate polyamide-induced inhibition of DNA replication in cell free extracts from X. laevis oocytes. In this system, inhibition of chromatin decondensation is observed, preventing DNA replication initiation. Finally, we show that Py-Im polyamides targeted to the ARE and ETS binding sequence downregulate AR- and ERG-driven signaling in a prostate cancer cell line harboring the TMPRSS2-ERG fusion. In a mouse xenograft model, ARE-targeted polyamide treatment reduced growth of the tumor.

On the interaction of actinomycin and DNA

Experiments have been accomplished that (a) further define the nature of the strong, G-containing DNA binding sites for actinomycin D (AMD), and (b) quantitate the in vitro inhibition of E. coli RNA polymerase activity by T7 DNA-bound AMD.

Twenty-five to forty percent of the G's of crab dAT are disallowed as strong AMD binding sites. The G's are measured to be randomly distributed, and, therefore, this datum cannot be explained on the basis of steric interference alone. Poly dAC:TG binds as much AMD and as strongly as any natural DNA, so the hypothesis that the unique strong AMD binding sites are G and a neighboring purine is incorrect. The datum can be explained on the basis of both steric interference and the fact that TGA is a disallowed sequence for strong AMD binding.

Using carefully defined in vitro conditions, there is one RNA synthesized per T7 DNA by E. coli RNA polymerase. The rate of the RNA polymerase-catalyzed reaction conforms to the equation 1/rate = 1/k_A(ATP) + 1/K_G(GTP) + 1/K_C(CTP) + 1/K_U(UTP) T7 DNA-bound AMD has only modest effects on initiation and termination of the polymerase-catalyzed reaction, but a large inhibitory effect on propagation. In the presence of bound AMD, k_G and k_C are decreased, whereas k_A and k_U are unaffected. These facts are interpreted to mean that on the microscopic level, on the average, the rates of incorporation of ATP and UTP are the same in the absence or presence of bound AMD, but that the rates of incorporation of GTP and CTP are decreased in the presence of AMD.

The interactions of basic proteins and DNA

I. ELECTROPHORESIS OF THE NUCLEIC ACIDS

A zone electrophoresis apparatus using ultraviolet optics has been constructed to study nucleic acids at concentrations less than 0.004%. Native DNA has a mobility about 15% higher than denatured DNA over a range of conditions. Otherwise, the electrophoretic mobility is independent of molecular weight, base composition or source. DNA mobilities change in the expected way with pH but the fractional change in mobility is less than the calculated change in charge. A small decrease in mobility accompanies an increase in ionic strength. RNA’s from various sources have mobilities slightly lower than denatured DNA except for s-RNA which travels slightly faster. The important considerations governing the mobility of nucleic acids appear to be the nature of the hydrodynamic segment, and the binding of counterions. The differences between electrophoresis and sedimentation stem from the fact that all random coil polyelectrolytes are fundamentally free draining in electrophoresis.

II. THE CYTOCHROME C/DNA COMPLEX

The basic protein, cytochrome c, has been complexed to DNA. Up to a cytochrome:DNA mass ratio of 2, a single type of complex is formed. Dissociation of this complex occurs between 0.05F and 0.1F NaCl. The complexing of cytochrome to DNA causes a slight increase in the melting temperature of the DNA, and a reduction of the electrophoretic mobility proportional to the decrease in net charge. Above a cytochrome:DNA mass ratio of 2.5, a different type of complex is formed. The results suggest that complexes such as are formed in the Kleinschmidt technique of electron microscopy would not exist in bulk solution and are exclusively film phenomena.

III. STUDIES OF THE ELECTROPHORESIS AND MELTING BEHAVIOUR OF NUCLEOHISTONES

Electrophoresis studies on reconstituted nucleohistones indicate that the electrophoretic mobility for these complexes is a function of the net charge of the complex. The mobility is therefore dependent on the charge density of the histone complexing the DNA, as well as on the histone/DNA ratio. It is found that the different histones affect the transition from native to denatured DNA in different ways. It appears that histone I is exchanging quite rapidly between DNA molecules in 0.01 F salt, while histone II is irreversibly bound. Histone III-IV enhances the capacity of non-strand separated denatured DNA to reanneal. Studies on native nucleoproteins indicate that there are no gene-sized uncomplexed DNA regions in any preparations studied.

IV. THE DISSOCIATION OF HISTONE FROM CALF THYMUS CROMATIN

Calf thymus nucleoprotein was treated with varying concentrations of NaCl. The identity of the histones associated and dissociated from the DNA at each salt concentration was determined by gel electrophoresis. It was found that there is no appreciable histone dissociation below 0.4 F NaCl. The lysine rich histones dissociate between 0.4 and 0.5 F NaCl. Their dissociation is accompanies by a marked increase in the solubility of the chromatin. The moderately lysine rich histones dissociate mainly between 0.8 and 1.1 F NaCl. There are two arginine rich histone components: the first dissociates between 0.8 F and 1.1 F NaCl, but the second class is the very last to be dissociated from the DNA (dissociation beginning at 1.0 F NaCl). By 2.0 F NaCl, essentially all the histones are dissociated.

The properties of the extracted nucleoprotein were studied. The electrophoretic mobility increases and the melting temperature decreases as more histones are dissociated from the DNA. A comparison with the dissociation of histones from DNA in NaClO₄ shows that to dissociate the same class of histones, the concentration of NaCl required is twice that of NaClO₄.

Differenzierung von Pangasius (Pangasius hypophthalmus) und Asiatischem Rotflossenwels (Hemibagrus wyckioides) durch Protein- und DNA-Analyse

In the last years farmed Pangasius (Tra-Pangasius, Pangasius hypophthalmus) from Vietnam has reached a considerable market share, whereas aquaculture of Asian Redtail Catfish (Hemibagrus wyckioides) is in its infancy. Recently it has been detected by food control authorities in Hamburg, that Pangasius fillets have been mislabelled and sold as fillets produced from Asian Redtail catfish. The necessity to improve the analytical methods for differentiation of Pangasius and Redtail Catfish prompted us to evaluate the suitability of isoelectric focusing (IEF) and DNA-analysis for identification of the two species. IEF of water soluble proteins was found to be a fast, reliable and economical method for differentiation of raw fillets of Pangasius and Redtail Catfish, as long as reference material is available. PCR-based DNA analysis was performed as follows: (i) amplification of a 464 bp segment of the cytochrome b gene; (ii) sequencing of the PCR product; (iii) comparison of the sequence with entries in GenBank using BLAST. The sequences of both species differed considerably, allowing the unequivocal differentiation between P. hypophthalmus and H. wyckioides. Kurzfassung Pangasius (Schlankwels, Tra-Pangasius, Pangasius hypophthalmus) hat sich innerhalb weniger Jahre zu einem bedeutenden Zuchtfisch entwickelt, während die Aquakultur des Asiatischen Rotflossenwelses (Hemibagrus wyckioides) in Vietnam noch in einem relativ kleinen Maßstab stattfindet. Kürzlich wurde von der Lebensmittelüberwachung in Hamburg nachgewiesen, dass im Handel erhältliche Filets mit der Deklaration „Rotflossenwels“ aus Pangasius hergestellt worden waren. Vor diesem Hintergrund wurden zwei Methoden auf ihre Eignung zur Differenzierung von Pangasius und Rotflossenwels geprüft. Es zeigte sich, dass sowohl die isoelektrische Fokussierung (IEF) wasserlöslicher Proteine als auch die PCR-basierte DNA-Analyse zur Unterscheidung beider Arten gut geeignet ist. Die IEF stellt eine schnelle und kostengünstige Untersuchungsmethode dar, die allerdings Referenzmaterial benötigt. Mit Hilfe der PCR (Polymerase-Kettenreaktion) wurde ein Abschnitt des Cytochrom b-Gens vervielfältigt und sequenziert. Die Sequenzen von P. hypophthalmus und H. wyckioides wiesen beträchtliche Unterschiede auf. Es wird diskutiert, wie sich durch Vergleich dieser Sequenzen mit Einträgen in Gendatenbanken unbekannte Proben beider Arten sicher zuordnen lassen.

Bases fisiológicas del crecimiento diferencial en especies de almeja del género Ruditapes: estudio de alometrías

La coexistencia de las especies R. philippinarum y R. decussatus en la costa atlántica, empleando los mismos recursos y con tasas de crecimiento diferenciadas, parece indicio de estrategias fisiológicas diferentes. Dichas estrategias pueden identificarse mediante la cuantificación de parámetros fisiológicos que intervienen en el crecimiento, el cual es resultante del balance energético. En el presente estudio, se ha hecho una aproximación a ello, eligiendo como parámetros a analizar, la tasa de aclaramiento (componente de ganancia energética) y la tasa metabólica (componente de pérdida energética). El crecimiento diferencial viene determinado por diferentes combinaciones de comportamientos fisiológicos, recogidos en modelos fisiológicos de crecimiento (modelo de adquisición de energía, de limitación de costes basales y de eficiencia metabólica). Con el objetivo de realizar una comparación inter-específica entre ambas especies, se ha estudiado además la relación alométrica de estos parámetros con el tamaño, ya que la distribución de tamaños es una característica de la población y constituye la base sobre la que se construyen curvas de crecimiento. De acuerdo a informes emitidos por la FAO, R. philippinarum presenta mayores tasas de crecimiento por lo que cabría esperar mayor tasa de aclaramiento y menor tasa metabólica en esta especie. Sin embargo los resultados no se corresponden con ello, siendo la tasa de aclaramiento mayor para R. decussatus y no encontrándose diferencias en la tasa metabólica. Las relaciones alométricas tampoco difieren entre especies. No se puede asociar la mayor tasa de crecimiento de R. philippinarum a ninguno de los modelos de crecimiento, ya que las diferencias inter-específicas en los parámetros estudiados son nulas o contrarias a lo esperado. La falta de disponibilidad de semillas de la especie R. decussatus para este estudio limita en parte estas conclusiones, ya que los patrones de crecimiento diferencial pueden diferir profundamente entre estadíos juvenil y adulto. Alternativamente las anomalías comentadas podrían explicarse por la presencia de híbridos de ambas especies entre los ejemplares estudiados y donde la heterosis del crecimiento podría disipar las diferencias entre especies.

Studies of dinucleoside monophosphates and monomer-polynucleotide interactions by proton magnetic resonance

The nature of the intra- and intermolecular base-stacking interactions involving several dinucleoside monophosphates in aqueous solution have been investigated by proton magnetic resonance spectrosocopy, and this method has been applied to a study of the interaction of polyuridylic acid with purine and adenosine monomers.

The pmr spectra of adenylyl (3' → 5') cytidine (ApC) and cytidylyl (3' → 5') adenosine (CpA) have been studied as a function of concentration and temperature. The results of these studies indicate that the intramolecular base-stacking interactions between the adenine and cytosine bases of these dinucleoside monophosphates are rather strong, and that the stacking tendencies are comparable for the two sequence isomers. The chemical shifts of the cytosine H₅ and adenine H₂ protons, and their variations with temperature, were shown to be consistent with stacked conformations in which both bases of the dinucleoside monophosphates are preferentially oriented in the anti conformation as in similar dApdC, and dCpdA (dA = deoxyadenosine; dC = deoxycytidine) segments in double helical DNA. The intramolecular stacking interaction was found to have a pronounced effect on the conformations of the ribose moieties, and these conformational changes are discussed. The concentration studies indicate extensive self-association of these dinucleoside monophosphates, and analysis of the concentration data facilitated determination of the dimerization constant for the association process as well as the nature of the intermolecular complexes.

The dependence of the ribose conformation upon the extent of intramolecular base-stacking was used to demonstrate that the base-base interaction in cytidylyl (3' → 5') cytidine (CpC) is rather strong, while there appears to be little interaction between the two uracil bases of uridylyl (3' → 5') uridine (UpU).

Studies of the binding of purine to several ribose and deoxyribose dinucleoside monophosphates show that the mode of interaction is base-stacking, and evidence for the formation of a purine-dinucleoside monophosphate intercalated complex is presented. The purine proton resonances are markedly broadened in this complex, and estimates of the purine linewidths in the complex and the equilibrium constant for purine intercalation are obtained.

A study of the interaction of unsubstitued purine with polyuridylic acid at 29°C by pmr indicated that purine binds to the uracil bases of the polymer by base-stacking. The severe broadening of the purine proton resonances observed provides strong evidence for the intercalation of purine between adjacent uracil bases of poly U. This interaction does not result in a more rigid or ordered structure for the polymer.

Investigation of the interaction between adenosine and polyuridylic acid revealed two modes of interaction between the monomer and the polymer, depending on the temperature. At temperatures above 26°C or so, monomeric adenosine binds to poly U by noncooperative A-U base stacking. Below this temperature, a rigid triple-stranded 1A:2U complex is formed, presumably via cooperative hydrogen-bonding as has previously been reported.

These results clearly illustrate the importance of base-stacking in non-specific interactions between bases, nucleosides and nucleotides, and also reveal the important role of the base-stacking interactions in cooperatively for med structures involving specific base-pairing where both types of interaction are possible.

Isolation and characteristics of a bacteriophage for bacillus stearothermophilus

A bacteriophage (TØ3) which infects the thermophilic bacterium Bacillus stearothermophilus ATCC 8005 was isolated and characterized. Infection of the bacterium by the bacteriophage was carried out at 60°C, the optimum growth temperature of the host. At 60°C the phage has a latent period of 18 minutes and a burst size of about 200. The phage is comparatively thermostable in broth. The half life of the phage is 400 minutes at 60°C, 120 minutes at 65°C, 40 minutes at 70°C and 12 minutes at 75°C. The activation energy for the heat inactivation of TØ3 is 56,000 cal. The buoyant density of TØ3 in a cesium chloride density gradient is 1.526.

Electron micrographs of TØ3 indicate that the phage has a regular hexagonal shaped head 57 mμ long. The morphology of the head is compatible with icosahedral symmetry. Each edge of the head is 29 mμ long, and there are 6 or 7 subunits along each edge. The tail of TØ3 is 125 mμ long and 10 mμ wide. There are about 30 cross striations that are spaced at 3.9 mμ intervals along the tail.

The DNA of phage TØ3 has a melting temperature of 88.5°C. Heat denatured TØ3 DNA can be extensively annealed in a high ionic strength environment. The buoyant density of TØ3 DNA in a cesium chloride density gradient is 1.695. TØ3 DNA contains: 42.7% guanine plus cytosine, as determined from the melting temperature; 43% guanine plus cytosine, as determined from the buoyant density; and 40.2% guanine plus cytosine, as determined by chromatographic separation and spectrophotometric estimation of the bases. The molecular weight of TØ3 DNA is 16.7 X 10⁶ as determined from the band width of the TØ3 DNA concentration distribution in a cesium chloride density gradient. Electron microscopy of TØ3 DNA revealed a single linear molecule that is 11.7 μ long. This corresponds to a molecular weight of 22.5 X 10⁶.

Heat denatured TØ3 DNA forms two bands in a cesium chloride density gradient, one at a density of 1.707 and the other at a density of 1.715. After the separated bands are mixed and annealed in the centrifuge cell, the renatured TØ3 DNA forms a single band at a density of 1.699. These results indicate that the two complementary strands of TØ3 DNA have different buoyant densities in cesium chloride, presumably because they have different base compositions.

The characteristics of TØ3 are compared with those of other phages. A hypothesis is presented for a relationship between the base composition of one strand of TØ3 DNA and the amino acid composition of the proteins of TØ3.

I. The structure and properties of closed circular duplex DNA. II. Chemically interacting systems at equilibrium in a buoyant density gradient

I. The binding of the intercalating dye ethidium bromide to closed circular SV 40 DNA causes an unwinding of the duplex structure and a simultaneous and quantitatively equivalent unwinding of the superhelices. The buoyant densities and sedimentation velocities of both intact (I) and singly nicked (II) SV 40 DNAs were measured as a function of free dye concentration. The buoyant density data were used to determine the binding isotherms over a dye concentration range extending from 0 to 600 µg/m1 in 5.8 M CsCl. At high dye concentrations all of the binding sites in II, but not in I, are saturated. At free dye concentrations less than 5.4 µg/ml, I has a greater affinity for dye than II. At a critical amount of dye bound I and II have equal affinities, and at higher dye concentration I has a lower affinity than II. The number of superhelical turns, τ, present in I is calculated at each dye concentration using Fuller and Waring's (1964) estimate of the angle of duplex unwinding per intercalation. The results reveal that SV 40 DNA I contains about -13 superhelical turns in concentrated salt solutions.

The free energy of superhelix formation is calculated as a function of τ from a consideration of the effect of the superhelical turns upon the binding isotherm of ethidium bromide to SV 40 DNA I. The value of the free energy is about 100 kcal/mole DNA in the native molecule. The free energy estimates are used to calculate the pitch and radius of the superhelix as a function of the number of superhelical turns. The pitch and radius of the native I superhelix are 430 Å and 135 Å, respectively.

A buoyant density method for the isolation and detection of closed circular DNA is described. The method is based upon the reduced binding of the intercalating dye, ethidium bromide, by closed circular DNA. In an application of this method it is found that HeLa cells contain in addition to closed circular mitochondrial DNA of mean length 4.81 microns, a heterogeneous group of smaller DNA molecules which vary in size from 0.2 to 3.5 microns and a paucidisperse group of multiples of the mitochondrial length.

II. The general theory is presented for the sedimentation equilibrium of a macromolecule in a concentrated binary solvent in the presence of an additional reacting small molecule. Equations are derived for the calculation of the buoyant density of the complex and for the determination of the binding isotherm of the reagent to the macrospecies. The standard buoyant density, a thermodynamic function, is defined and the density gradients which characterize the four component system are derived. The theory is applied to the specific cases of the binding of ethidium bromide to SV 40 DNA and of the binding of mercury and silver to DNA.

Proton magnetic resonance studies of ribonucleic acid complexes. I. Complexes of biological bases and oligonucleotides with RNA. II. Template recognition and the degeneracy of the genetic code

Part I. Complexes of Biological Bases and Oligonucleotides with RNA

The physical nature of complexes of several biological bases and oligonucleotides with single-stranded ribonucleic acids have been studied by high resolution proton magnetic resonance spectroscopy. The importance of various forces in the stabilization of these complexes is also discussed.

Previous work has shown that purine forms an intercalated complex with single-stranded nucleic acids. This complex formation led to severe and stereospecific broadening of the purine resonances. From the field dependence of the linewidths, T₁ measurements of the purine protons and nuclear Overhauser enhancement experiments, the mechanism for the line broadening was ascertained to be dipole-dipole interactions between the purine protons and the ribose protons of the nucleic acid.

The interactions of ethidium bromide (EB) with several RNA residues have been studied. EB forms vertically stacked aggregates with itself as well as with uridine, 3'-uridine monophosphate and 5'-uridine monophosphate and forms an intercalated complex with uridylyl (3' → 5') uridine and polyuridylic acid (poly U). The geometry of EB in the intercalated complex has also been determined.

The effect of chain length of oligo-A-nucleotides on their mode of interaction with poly U in D₂0 at neutral pD have also been studied. Below room temperatures, ApA and ApApA form a rigid triple-stranded complex involving a stoichiometry of one adenine to two uracil bases, presumably via specific adenine-uracil base pairing and cooperative base stacking of the adenine bases. While no evidence was obtained for the interaction of ApA with poly U above room temperature, ApApA exhibited complex formation of a 1:1 nature with poly U by forming Watson-Crick base pairs. The thermodynamics of these systems are discussed.

Part II. Template Recognition and the Degeneracy of the Genetic Code

The interaction of ApApG and poly U was studied as a model system for the codon-anticodon interaction of tRNA and mRNA in vivo. ApApG was shown to interact with poly U below ~20°C. The interaction was of a 1:1 nature which exhibited the Hoogsteen bonding scheme. The three bases of ApApG are in an anti conformation and the guanosine base appears to be in the lactim tautomeric form in the complex.

Due to the inadequacies of previous models for the degeneracy of the genetic code in explaining the observed interactions of ApApG with poly U, the "tautomeric doublet" model is proposed as a possible explanation of the degenerate interactions of tRNA with mRNA during protein synthesis in vivo.

Studies on T-even bacteriophage DNA

Part I. The regions of sequence homology and non-homology between the DNA molecules of T2, T4, and T6 have been mapped by the electron microscopic heteroduplex method. The heteroduplex maps have been oriented with respect to the T4 genetic map. They show characteristic, reproducible patterns of substitution and deletion loops. All heteroduplex molecules show more than 85% homology. Some of the loop patterns in T2/T4 heteroduplexes are similar to those in T4/T6.

We find that the rII, the lysozyme and ac genes, the D region, and gene 52 are homologous in T2, T4, and T6. Genes 43 and 47 are probably homologous between T2 and T4. The region of greatest homology is that bearing the late genes. The host range region, which comprises a part of gene 37 and all of gene 38, is heterologous in T2, T4, and T6. The remainder of gene 37 is partially homologous in the T2/T4 heteroduplex (Beckendorf, Kim and Lielausis, 1972) but it is heterologous in T4/T6 and in T2/T6. Some of the tRNA genes are homologous and some are not. The internal protein genes in general seem to be non-homologous.

The molecular lengths of the T-even DNAs are the same within the limit of experimental error; their calculated molecular weights are correspondingly different due to unequal glucosylation. The size of the T2 genome is smaller than that of T4 or T6, but the terminally repetitious region in T2 is larger. There is a length distribution of the terminal repetition for any one phage DNA, indicating a variability in length of the DNA molecules packaged within the phage.

Part II. E. coli cells infected with phage strains carrying extensive deletions encompassing the gene for the phage ser-tRNA are missing the phage tRNAs normally present in wild type infected cells. By DNA-RNA hybridization we have demonstrated that the DNA complementary to the missing tRNAs is also absent in such deletion mutants. Thus the genes for these tRNAs must be clustered in the same region of the genome as the ser-tRNA gene. Physical mapping of several deletions of the ser-tRNA and lysozyme genes, by examination of heteroduplex DNA in the electron microscope, has enabled us to locate the cluster, to define its maximum size, and to order a few of the tRNA genes within it. That such deletions can be isolated indicates that the phage-specific tRNAs from this cluster are dispensable.

Part III. Genes 37 and 38 between closely related phages T2 and T4 have been compared by genetic, biochemical, and hetero-duplex studies. Homologous, partially homologous and non-homologous regions of the gene 37 have been mapped. The host range determinant which interacts with the gene 38 product is identified.

Part IV. A population of double-stranded ØX-RF DNA molecules carrying a deletion of about 9% of the wild-type DNA has been discovered in a sample cultivated under conditions where the phage lysozyme gene is nonessential. The structures of deleted monomers, dimers, and trimers have been studied by the electron microscope heteroduplex method. The dimers and trimers are shown to be head-to-tail repeats of the deleted monomers. Some interesting examples of the dynamical phenomenon of branch migration in vitro have been observed in heteroduplexes of deleted dimer and trimer strands with undeleted wild-type monomer viral strands.

Desenvolvimento e validação de sistema multiplex X-STR para identificação humana por DNA

Novas metodologias de análise molecular voltadas para estudos populacionais, clínicos, evolutivos, da biodiversidade e identificação forense foram desenvolvidas com base em marcadores microssátelites ou STR Short Tandem Repeats. Os marcadores STR, que estão amplamente espalhados nos genomas e se caracterizam por apresentar alto grau de polimorfismo, podem ser analisados a partir da amplificação por PCR (Reação em Cadeia da polimerase). A análise foi facilitada a partir do desenvolvimento de sistemas de amplificação simultânea de múltiplos STR (multiplex STR) e com a detecção automatizada dos produtos de amplificação marcados por fluorescência. Recentemente, o uso de marcadores STR do cromossomo X (X-STR) tornou-se significativo na prática forense. Devido ao seu modo de transmissão, os X-STR são úteis em situações particulares de investigação de relações de parentesco, apresentando vantagens sobre o uso de STR autossômicos. Este estudo teve como principal objetivo o desenvolvimento e validação de sistema multiplex, denominado LDD (X-STR) Decaplex, capaz de amplificar dez loci X-STR (DXS7133, DXS7424, DXS8378, DXS6807, DXS7132, DXS10074, DXS7423, DXS8377, GATA172D05 e DXS10101) para aplicação em genética populacional, identificação e análises forenses. Utilizando o LDD (X-STR) Decaplex 170 indivíduos autodenominados afrodescendentes, não aparentados geneticamente, foram genotipados. As freqüências alélicas e genotípicas não apresentaram desvio do equilíbrio de Hardy-Weinberg e estão em concordância com aquelas observadas em outros estudos. Os haplótipos observados foram únicos em indivíduos de amostra masculina. A análise de desequilíbrio de ligação não revelou associação entre os marcadores X-STR. A diversidade genética foi elevada, variando entre 0,6218 para o locus DXS7133 a 0,9327 para o locus DXS8377. Os parâmetros de Probabilidade de Vinculação (PV), Índice de Vinculação (IV), Poder de Exclusão (PE), Poder de Discriminação e Razão de Verossimilhança foram também elevados, demonstraram que os dez X-STRs são altamente polimórficos e discriminativos na população estudada. A concentração mínima de DNA para a amplificação dos loci do LDD (X-STR) Decaplex é de 0,5 ng e verificamos que amplificação por PCR pode ser afetada quando são adicionados mais de 5 ng de DNA nas reações. Os percentuais de bandas stutter foram elevados para os loci DXS7132 e DXS8377. No teste de reprodutibilidade observamos consistência entre as tipagem de diferentes amostras biológicas, incluindo as de restos mortais. No teste de mistura a proporção limite em que observamos a coexistência de duas espécies biológicas foi de 2,5:1ng (feminino-masculino). Os resultados evidenciaram que os loci do LDD (X-STR) Decaplex são altamente informativos, consistindo, em conjunto, uma ferramenta importante em estudos de identificação humana e de relações de parentesco.

Avaliação da suscetibilidade antimicrobiana e do perfil de macrorrestrição do DNA de isolados humanos de Salmonella serovar Typhimurium fagotipo 193, no período de 1970 a 2008

Para analisar cepas de Salmonella ser. Typhimurium isoladas de processos entéricos e extraintestinais humanos ocorridos no período de 1970 a 2008 de diferentes regiões do país foram selecionadas, com base nos registros contidos no banco de dados do Laboratório de Enterobactérias do IOC/FIOCRUZ, RJ, amostras do fagotipo prevalente 193, visando precipuamente o reconhecimento de clones epidêmicos. Foram selecionadas 553 cepas de Salmonella ser. Typhimurium fagotipo 193 representadas por 91, 65, 70 e 327 amostras referentes as décadas de 70, 80, 90 e ao período de 2000 a 2008, respectivamente. Na análise global da sensibilidade destas cepas, 52% apresentaram um ou mais marcadores de resistência a antibióticos incluídos no perfil ACSSuT. Este perfil de resistência completo foi verificado em 20,9% dos isolados, sendo os 21,9% restantes, sensíveis a todas as drogas testadas, especialmente no período de 2000 a 2008, representadas por 121 amostras (37,0%) em relação as 327 culturas dessa época. O maior percentual de resistência foi observado nas amostras da década de 70 (99%) sendo o perfil ACSSuT detectado em 35,2% dos isolados, ressaltando-se que todas as amostras foram isoladas de processos gastroentéricos ocorridos na cidade de São Paulo. Ao longo das quatro décadas de estudo, descreve-se um ponto de ruptura entre a prevalência de resistência e a suscetibilidade na transição entre as décadas de 80 e 90. Embora o número de isolados de Salmonella ser. Typhimurium fagotipo 193 tenha aumentado no último período considerado, o percentual de mono e multirresistência aos antimicrobianos se situou em nível elevado (63,0%). A análise do polimorfismo obtido após macrorrestrição com a enzima XbaI revelou que cepas isoladas na década de 90 apresentaram elevado percentual de similaridade (≥85%) com cepas isoladas recentemente (período de 2000-2008), sendo agrupadas nos mesmos subclusters. Por outro lado, as cepas da década de 70 inserem-se em subclusters independentes, embora o percentual de similaridade entre tais subclusters e os demais seja ≥70%; o mesmo sendo observado para as cepas isoladas durante a década de 80. Em conclusão, este estudo mostrou que a prevalência de isolados humanos de Salmonella ser. Typhimurium fagotipo 193 no Brasil vem progredindo desde a década de 1990, enquanto a detecção do modelo R (ACSSuT) está diminuindo e a avaliação através da PFGE indicou a presença de multiplicidade de perfis de macrorrestrição no fagotipo 193, entretanto com elevados percentuais de similaridade entre si, sugerindo alguma clonalidade, tendo em vista o período entre o isolamento e a análise