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.

Relay lens-free birefringence-customized stackable beam-steering modules for optical interconnects

A previously suggested birefringence-customized modular optical interconnect technique is extended for lens-free relay operation. Various lens-free relay imaging models are developed. We claim that the lens-free relay system is important in simplifying an optical interconnect system whenever the imaging conditions permit. To verify the validity of various proposed concepts, we experimentally implemented some 8 x 8 optical permutation modules. High-power efficiency and low channel cross talk were experimentally observed. In general, the larger the channel spacing, the less the cross talk. A quantitative cross-talk measurement of the lens-free relay system shows that, for a fixed channel width of 0.5 mm and channel spacings of 0.5, 1, and 2 mm, a less than -20-dB cross-talk performance can be guaranteed for lens-free relay distances of 40, 280, and 430 mm, respectively. (C) 1998 Optical Society of America.

An in vitro biomolecular breadboard for prototyping synthetic biological circuits

Biomolecular circuit engineering is critical for implementing complex functions in vivo, and is a baseline method in the synthetic biology space. However, current methods for conducting biomolecular circuit engineering are time-consuming and tedious. A complete design-build-test cycle typically takes weeks' to months' time due to the lack of an intermediary between design ex vivo and testing in vivo. In this work, we explore the development and application of a "biomolecular breadboard" composed of an in-vitro transcription-translation (TX-TL) lysate to rapidly speed up the engineering design-build-test cycle. We first developed protocols for creating and using lysates for conducting biological circuit design. By doing so we simplified the existing technology to an affordable ($0.03/uL) and easy to use three-tube reagent system. We then developed tools to accelerate circuit design by allowing for linear DNA use in lieu of plasmid DNA, and by utilizing principles of modular assembly. This allowed the design-build-test cycle to be reduced to under a business day. We then characterized protein degradation dynamics in the breadboard to aid to implementing complex circuits. Finally, we demonstrated that the breadboard could be applied to engineer complex synthetic circuits in vitro and in vivo. Specifically, we utilized our understanding of linear DNA prototyping, modular assembly, and protein degradation dynamics to characterize the repressilator oscillator and to prototype novel three- and five-node negative feedback oscillators both in vitro and in vivo. We therefore believe the biomolecular breadboard has wide application for acting as an intermediary for biological circuit engineering.

Proton-Coupled Reduction of N2 Facilitated by Molecular Fe Complexes

The activation of Fe-coordinated N₂ via the formal addition of hydrogen atom equivalents is explored in this thesis. These reactions may occur in nitrogenase enzymes during the biological conversion of N₂ to NH₃. To understand these reactions, the N₂ reactivity of a series of molecular Fe(N₂) platforms is investigated. A trigonal pyramidal, carbon-ligated Fe^I complex was prepared that displays a similar geometry to that of the resting state 'belt' Fe atoms of nitrogenase. Upon reduction, this species was shown to coordinate N₂, concomitant with significant weakening of the C-Fe interaction. This hemilability of the axial ligand may play a critical role in mediating the interconversion of Fe(N_xH_y) species during N₂ conversion to NH₃. In fact, a trigonal pyramidal borane-ligated Fe complex was shown to catalyze this transformation, generating up to 8.49 equivalents of NH₃. To shed light on the mechanistic details of this reaction, protonation of a borane-ligated Fe(N₂) complex was investigated and found to give rise to a mixture of species that contains an iron hydrazido(2-) [Fe(NNH₂)] complex. The identification of this species is suggestive of an early N-N bond cleavage event en route to NH₃ production, but the highly-reactive nature of this complex frustrated direct attempts to probe this possibility. A structurally-analogous silyl-ligated Fe(N₂) complex was found to react productively with hydrogen atom equivalents, giving rise to an isolable Fe(NNH₂) species. Spectroscopic and crystallographic studies benefited from the enhanced stability of this complex relative to the borane analogue. One-electron reduction of this species initiates a spontaneous disproportionation reaction with an iron hydrazine [Fe(NH₂NH₂)] complex as the predominant reaction product. This transformation provides support for an Fe-mediated N₂ activation mechanism that proceeds via a late N-N bond cleavage. In hopes of gaining more fundamental insight into these reactions, a series of Fe(CN) complexes were prepared and reacted with hydrogen-atom equivalents. Significant quantities of CH₄ and NH₃ are generated in these reactions as a result of complete C-N bond activation. A series of Fe(CNH_x) were found to be exceptionally stable and may be intermediates in these reactions. The stability of these compounds permitted collection of thermodynamic parameters pertinent to the unique N-H bonds. This data is comparatively discussed with the theoretically-predicted data of the N₂-derived Fe(NNH_x) species. Exceptionally-weak N-H bond enthalpies are found for many of these compounds, and sheds light on their short-lived nature and tendency to evolve H₂. As a whole, these works both establish and provide a means to understand Fe-mediated N₂ activation via the addition of hydrogen atom equivalents.

Metallogenic Provinces of the Southwestern United States and Northern Mexico

Some of the metallogenic provinces of the southwestern United States and northern Mexico are defined by the geographic distribution of trace elements in the primary sulfide minerals chalcopyrite and sphalerite. The elements investigated include antimony, arsenic, bismuth, cadmium, cobalt, gallium, germanium, indium, manganese, molybdenum, nickel, silver, tellurium, thallium, and tin. Of these elements, cobalt, gallium, germanium, indium, nickel, silver, and tin exhibit the best defined geographic distribution.

The data indicate that chalcopyrite is the preferred host for tin and perhaps molybdenum; sphalerite is the preferred host for cadmium, gallium, germanium, indium, and manganese; galena is the preferred host for antimony, bismuth, silver, tellurium, and thallium; and pyrite is the preferred host for cobalt, nickel, and perhaps arsenic. With respect to the two minerals chalcopyrite and sphalerite, antimony, arsenic, molybdenum, nickel, silver, and tin prefer chalcopyrite; and bismuth, cadmium, cobalt, gallium, germanium, indium, manganese, and thallium prefer sphalerite. This distribution probably is the result of the interaction of several factors, among which are these: the various radii of the elements, the association due to chemical similarities of the major and trace elements, and the degree of ionic versus covalent and metallic character of the metal-sulfur bonds in chalcopyrite and sphalerite. The type of deposit, according to a temperature classification, appears to be of minor importance in determining the trace element content of chalcopyrite and sphalerite.

A preliminary investigation of large single crystals of sphalerite and chalcopyrite indicates that the distribution within a single crystal of some elements such as cadmium in sphalerite and indium and silver in chalcopyrite is relatively uniform, whereas the distribution of some other elements such as cobalt and manganese in sphalerite is somewhat less uniform and the distribution of tin in sphalerite is extremely erratic. The variations in trace element content probably are due largely to variations in the composition of the fluids during the growth of the crystals, but the erratic behavior of tin in sphalerite perhaps is related to the presence of numerous cavities and inclusions in the crystal studied.

Maps of the geographic distribution of trace elements in chalcopyrite and sphalerite exhibit three main belts of greater than average trace element content, which are called the Eastern, Central, and Western belts. These belts are consistent in trend and position with a beltlike distribution of copper, gold, lead, zinc, silver, and tungsten deposits and with most of the major tectonic features. However, there appear to be no definite time relationships, for as many as four metallogenic epochs, from Precambrian to late Tertiary, are represented by ore deposits within the Central belt.

The evidence suggests that the beltlike features have a deep seated origin, perhaps in the sub-crust or outer parts of the mantle, and that the deposits within each belt might be genetically related through a beltlike compositional heterogeneity in the source regions of the ores. Hence, the belts are regarded as metallogenic provinces.

Combinatorial properties of finite geometric lattices

Let L be a finite geometric lattice of dimension n, and let w(k) denote the number of elements in L of rank k. Two theorems about the numbers w(k) are proved: first, w(k) ≥ w(1) for k = 2, 3, ..., n-1. Second, w(k) = w(1) if and only if k = n-1 and L is modular. Several corollaries concerning the "matching" of points and dual points are derived from these theorems.

Both theorems can be regarded as a generalization of a theorem of de Bruijn and Erdös concerning ʎ= 1 designs. The second can also be considered as the converse to a special case of Dilworth's theorem on finite modular lattices.

These results are related to two conjectures due to G. -C. Rota. The "unimodality" conjecture states that the w(k)'s form a unimodal sequence. The "Sperner" conjecture states that a set of non-comparable elements in L has cardinality at most max/k {w(k)}. In this thesis, a counterexample to the Sperner conjecture is exhibited.

Modelagem, simulação e otimização operacional de unidades de síntese de metanol

A indústria de processos químicos tem sofrido consideráveis transformações devido ao acirramento da competitividade. Importantes progressos tecnológicos têm sido atingidos através de técnicas de modelagem, simulação e otimização visando o aumento da lucratividade e melhoria contínua nos processos industriais. Neste contexto, as plantas de metanol, um dos mais importantes produtos petroquímicos, podem ser destacadas. Atualmente, a principal matéria-prima para obtenção de metanol é o gás natural. A produção do metanol é caracterizada por três etapas: geração de gás de síntese, conversão do gás de síntese em metanol (unidade de síntese ou loop de síntese) e purificação do produto na especificação requerida. Os custos fixos e variáveis da unidade de síntese são fortemente dependentes das variáveis operacionais, como temperatura, pressão, razão de reciclo e composição da carga. Desta forma, foi desenvolvido um conjunto de modelos e algoritmos computacionais para representar matematicamente unidades de síntese de metanol. O modelo apresenta operações unitárias associadas aos seguintes equipamentos: divisores de correntes, misturadores de correntes, compressores, trocadores de calor, vasos de flash e reatores. Inicialmente, foi proposto um simulador estacionário, que serviu como base para um pseudo-estacionário, o qual contempla a desativação do catalisador por sinterização térmica. Os simuladores foram criados segundo uma arquitetura seqüencial modular e empregou-se o método de substituição sucessiva para a convergência dos reciclos. O estudo envolveu dois fluxogramas típicos, um constituído por reatores adiabáticos em série, enquanto o outro constituído por um reator tipo quench. Uma análise do efeito das principais variáveis operacionais foi realizada para o reator e para o loop de síntese. Estudou-se também o efeito da desativação do catalisador ao longo do tempo. Uma ferramenta de otimização operacional foi empregada para alcançar a máxima produção manipulando as injeções de carga fria na entrada dos leitos catalíticos. Formulou-se também um problema de maximização do lucro em função da otimização de paradas da unidade para troca do catalisador. Os resultados obtidos apontam que a ferramenta desenvolvida é promissora para a compreensão e otimização da unidade objeto deste estudo

Decoding cosets of first-order Reed-Muller code

Proper encoding of transmitted information can improve the performance of a communication system. To recover the information at the receiver it is necessary to decode the received signal. For many codes the complexity and slowness of the decoder is so severe that the code is not feasible for practical use. This thesis considers the decoding problem for one such class of codes, the comma-free codes related to the first-order Reed-Muller codes.

A factorization of the code matrix is found which leads to a simple, fast, minimum memory, decoder. The decoder is modular and only n modules are needed to decode a code of length 2ⁿ. The relevant factorization is extended to any code defined by a sequence of Kronecker products.

The problem of monitoring the correct synchronization position is also considered. A general answer seems to depend upon more detailed knowledge of the structure of comma-free codes. However, a technique is presented which gives useful results in many specific cases.