Synthetic molecular machines for active self-assembly : prototype algorithms, designs, and experimental study

Computer science and electrical engineering have been the great success story of the twentieth century. The neat modularity and mapping of a language onto circuits has led to robots on Mars, desktop computers and smartphones. But these devices are not yet able to do some of the things that life takes for granted: repair a scratch, reproduce, regenerate, or grow exponentially fast–all while remaining functional.

This thesis explores and develops algorithms, molecular implementations, and theoretical proofs in the context of “active self-assembly” of molecular systems. The long-term vision of active self-assembly is the theoretical and physical implementation of materials that are composed of reconfigurable units with the programmability and adaptability of biology’s numerous molecular machines. En route to this goal, we must first find a way to overcome the memory limitations of molecular systems, and to discover the limits of complexity that can be achieved with individual molecules.

One of the main thrusts in molecular programming is to use computer science as a tool for figuring out what can be achieved. While molecular systems that are Turing-complete have been demonstrated [Winfree, 1996], these systems still cannot achieve some of the feats biology has achieved.

One might think that because a system is Turing-complete, capable of computing “anything,” that it can do any arbitrary task. But while it can simulate any digital computational problem, there are many behaviors that are not “computations” in a classical sense, and cannot be directly implemented. Examples include exponential growth and molecular motion relative to a surface.

Passive self-assembly systems cannot implement these behaviors because (a) molecular motion relative to a surface requires a source of fuel that is external to the system, and (b) passive systems are too slow to assemble exponentially-fast-growing structures. We call these behaviors “energetically incomplete” programmable behaviors. This class of behaviors includes any behavior where a passive physical system simply does not have enough physical energy to perform the specified tasks in the requisite amount of time.

As we will demonstrate and prove, a sufficiently expressive implementation of an “active” molecular self-assembly approach can achieve these behaviors. Using an external source of fuel solves part of the the problem, so the system is not “energetically incomplete.” But the programmable system also needs to have sufficient expressive power to achieve the specified behaviors. Perhaps surprisingly, some of these systems do not even require Turing completeness to be sufficiently expressive.

Building on a large variety of work by other scientists in the fields of DNA nanotechnology, chemistry and reconfigurable robotics, this thesis introduces several research contributions in the context of active self-assembly.

We show that simple primitives such as insertion and deletion are able to generate complex and interesting results such as the growth of a linear polymer in logarithmic time and the ability of a linear polymer to treadmill. To this end we developed a formal model for active-self assembly that is directly implementable with DNA molecules. We show that this model is computationally equivalent to a machine capable of producing strings that are stronger than regular languages and, at most, as strong as context-free grammars. This is a great advance in the theory of active self- assembly as prior models were either entirely theoretical or only implementable in the context of macro-scale robotics.

We developed a chain reaction method for the autonomous exponential growth of a linear DNA polymer. Our method is based on the insertion of molecules into the assembly, which generates two new insertion sites for every initial one employed. The building of a line in logarithmic time is a first step toward building a shape in logarithmic time. We demonstrate the first construction of a synthetic linear polymer that grows exponentially fast via insertion. We show that monomer molecules are converted into the polymer in logarithmic time via spectrofluorimetry and gel electrophoresis experiments. We also demonstrate the division of these polymers via the addition of a single DNA complex that competes with the insertion mechanism. This shows the growth of a population of polymers in logarithmic time. We characterize the DNA insertion mechanism that we utilize in Chapter 4. We experimentally demonstrate that we can control the kinetics of this re- action over at least seven orders of magnitude, by programming the sequences of DNA that initiate the reaction.

In addition, we review co-authored work on programming molecular robots using prescriptive landscapes of DNA origami; this was the first microscopic demonstration of programming a molec- ular robot to walk on a 2-dimensional surface. We developed a snapshot method for imaging these random walking molecular robots and a CAPTCHA-like analysis method for difficult-to-interpret imaging data.

Pharmacology and pore-forming domains of the cystic fibrosis transmembrane conductance regulator

The cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel member of the ATP-binding cassette (ABC) superfamily of membrane proteins. CFTR has two homologous halves, each consisting of six transmembrane spanning domains (TM) followed by a nucleotide binding fold, connected by a regulatory (R) domain. This thesis addresses the question of which domains are responsible for Cl^- selectivity, i.e., which domains line the channel pore.

To address this question, novel blockers of CFTR were characterized. CFTR was heterologously expressed in Xenopus oocytes to study the mechanism of block by two closely related arylaminobenzoates, diphenylamine-2-carboxylic acid (DPC) and flufenamic acid (FFA). Block by both is voltage-dependent, with a binding site ≈ 40% through the electric field of the membrane. DPC and FFA can both reach their binding site from either side of the membrane to produce a flickering block of CFTR single channels. In addition, DPC block is influenced by Cl^- concentration, and DPC blocks with a bimolecular forward binding rate and a unimolecular dissociation rate. Therefore, DPC and FFA are open-channel blockers of CFTR, and a residue of CFTR whose mutation affects their binding must line the pore.

Screening of site-directed mutants for altered DPC binding affinity reveals that TM-6 and TM-12 line the pore. Mutation of residue 5341 in TM-6 abolishes most DPC block, greatly reduces single-channel conductance, and alters the direction of current rectification. Additional residues are found in TM-6 (K335) and TM-12 (T1134) whose mutations weaken or strengthen DPC block; other mutations move the DPC binding site from TM-6 to TM-12. The strengthened block and lower conductance due to mutation T1134F is quantitated at the single-channel level. The geometry of DPC and of the residues mutated suggest α-helical structures for TM-6 and TM-12. Evidence is presented that the effects of the mutations are due to direct side-chain interaction, and not to allosteric effects propagated through the protein. Mutations are also made in TM-11, including mutation S1118F, which gives voltage-dependent current relaxations. The results may guide future studies on permeation through ABC transporters and through other Cl^- channels.

Equilibrium in dynamic economic models

The main theme running through these three chapters is that economic agents are often forced to respond to events that are not a direct result of their actions or other agents actions. The optimal response to these shocks will necessarily depend on agents' understanding of how these shocks arise. The economic environment in the first two chapters is analogous to the classic chain store game. In this setting, the addition of unintended trembles by the agents creates an environment better suited to reputation building. The third chapter considers the competitive equilibrium price dynamics in an overlapping generations environment when there are supply and demand shocks.

The first chapter is a game theoretic investigation of a reputation building game. A sequential equilibrium model, called the "error prone agents" model, is developed. In this model, agents believe that all actions are potentially subjected to an error process. Inclusion of this belief into the equilibrium calculation provides for a richer class of reputation building possibilities than when perfect implementation is assumed.

In the second chapter, maximum likelihood estimation is employed to test the consistency of this new model and other models with data from experiments run by other researchers that served as the basis for prominent papers in this field. The alternate models considered are essentially modifications to the standard sequential equilibrium. While some models perform quite well in that the nature of the modification seems to explain deviations from the sequential equilibrium quite well, the degree to which these modifications must be applied shows no consistency across different experimental designs.

The third chapter is a study of price dynamics in an overlapping generations model. It establishes the existence of a unique perfect-foresight competitive equilibrium price path in a pure exchange economy with a finite time horizon when there are arbitrarily many shocks to supply or demand. One main reason for the interest in this equilibrium is that overlapping generations environments are very fruitful for the study of price dynamics, especially in experimental settings. The perfect foresight assumption is an important place to start when examining these environments because it will produce the ex post socially efficient allocation of goods. This characteristic makes this a natural baseline to which other models of price dynamics could be compared.

Signal transduction with hybridization chain reactions

Some of the most exciting developments in the field of nucleic acid engineering include the utilization of synthetic nucleic acid molecular devices as gene regulators, as disease marker detectors, and most recently, as therapeutic agents. The common thread between these technologies is their reliance on the detection of specific nucleic acid input markers to generate some desirable output, such as a change in the copy number of an mRNA (for gene regulation), a change in the emitted light intensity (for some diagnostics), and a change in cell state within an organism (for therapeutics). The research presented in this thesis likewise focuses on engineering molecular tools that detect specific nucleic acid inputs, and respond with useful outputs.

Four contributions to the field of nucleic acid engineering are presented: (1) the construction of a single nucleotide polymorphism (SNP) detector based on the mechanism of hybridization chain reaction (HCR); (2) the utilization of a single-stranded oligonucleotide molecular Scavenger as a means of enhancing HCR selectivity; (3) the implementation of Quenched HCR, a technique that facilitates transduction of a nucleic acid chemical input into an optical (light) output, and (4) the engineering of conditional probes that function as sequence transducers, receiving target signal as input and providing a sequence of choice as output. These programmable molecular systems are conceptually well-suited for performing wash-free, highly selective rapid genotyping and expression profiling in vitro, in situ, and potentially in living cells.

Chain-forming zintl antimonides as novel thermoelectric materials

Zintl phases, a subset of intermetallic compounds characterized by covalently-bonded "sub-structures," surrounded by highly electropositive cations, exhibit precisely the characteristics desired for thermoelectric applications. The requirement that Zintl compounds satisfy the valence of anions through the formation of covalent substructures leads to many unique, complex crystal structures. Such complexity often leads to exceptionally low lattice thermal conductivity due to the containment of heat in low velocity optical modes in the phonon dispersion. To date, excellent thermoelectric properties have been demonstrated in several Zintl compounds. However, compared with the large number of known Zintl phases, very few have been investigated as thermoelectric materials.

From this pool of uninvestigated compounds, we selected a class of Zintl antimonides that share a common structural motif: anionic moieties resembling infinite chains of linked MSb₄ tetrahedra, where $M$ is a triel element. The compounds discussed in this thesis (A₅M₂Sb₆ and A₃MSb₃, where A = Ca or Sr and M = Al, Ga and In) crystallize as four distinct, but closely related "chain-forming" structure types. This thesis describes the thermoelectric characterization and optimization of these phases, and explores the influence of their chemistry and structure on the thermal and electronic transport properties. Due to their large unit cells, each compound exhibits exceptionally low lattice thermal conductivity (0.4 - 0.6 W/mK at 1000 K), approaching the predicted glassy minimum at high temperatures. A combination of Density Functional calculations and classical transport models were used to explain the experimentally observed electronic transport properties of each compound. Consistent with the Zintl electron counting formalism, A₅M₂Sb₆ and A₃MSb₃ phases were found to have filled valence bands and exhibit intrinsic electronic properties. Doping with divalent transition metals (Zn²⁺ and Mn²⁺) on the M³⁺ site, or Na¹⁺ on the A³⁺ site allowed for rational control of the carrier concentration and a transition towards degenerate semiconducting behavior. In optimally-doped samples, promising peak zT values between 0.4 and 0.9 were obtained, highlighting the value of continued investigations of complex Zintl phases.

Synthesis and potency of long end-associative polymers for mist control

Long linear polymers that are end-functionalized with associative groups were studied as additives to hydrocarbon fluids to mitigate the fire hazard associated with the presence of mist in a crash scenario. These polymers were molecularly designed to overcome both the shear-degradation of long polymer chains in turbulent flows, and the chain collapse induced by the random placement of associative groups along polymer backbones. Architectures of associative groups on the polymer chain ends that were tested included clusters of self-associative carboxyl groups and pairs of hetero-complementary associative units.

Linear polymers with clusters of discrete numbers of carboxyl groups on their chain ends were investigated first: an innovative synthetic strategy was devised to achieve unprecedented backbone lengths and precise control of the number of carboxyl groups on chain ends (N). We found that a very narrow range of N allows the co-existence of sufficient end-association strength and polymer solubility in apolar media. Subsequent steady-flow rheological study on solution behavior of such soluble polymers in apolar media revealed that the end-association of very long chains in apolar media leads to the formation of flower-like micelles interconnected by bridging chains, which trap significant fraction of polymer chains into looped structures with low contribution to mist-control. The efficacy of very long 1,4-polybutadiene chains end-functionalized with clusters of four carboxyl groups as mist-control additives for jet fuel was further tested. In addition to being shear-resistant, the polymer was found capable of providing fire-protection to jet fuel at concentrations as low as 0.3wt%. We also found that this polymer has excellent solubility in jet fuel over a wide range of temperature (-30 to +70°C) and negligible interference with dewatering of jet fuel. It does not cause an adverse increase in viscosity at concentrations where mist-control efficacy exists.

Four pairs of hetero-complementary associative end-groups of varying strengths were subsequently investigated, in the hopes of achieving supramolecular aggregates with both mist-control ability and better utilization of polymer building blocks. Rheological study of solutions of the corresponding complementary associative polymer pairs in apolar media revealed the strength of complementary end-association required to achieve supramolecular aggregates capable of modulating rheological properties of the solution.

Both self-associating and complementary associating polymers have therefore been found to resist shear degradation. The successful strategy of building soluble, end-associative polymers with either self-associative or complementary associative groups will guide the next generation of mist-control technology.

New strategies for the total synthesis of aza-propellane natural products

The propellane alkaloids comprise a large class of natural products that possess varying degrees of structural complexity and biological activity. The earliest of these to be isolated was acutumine, a chlorinated alkaloid that has been shown to exhibit selective T-cell cytotoxicity and antiamnesic properties. Alternatively, the hasubanan family of natural products has garnered considerable attention from the synthetic community in part due to its structural similarities to morphine. While these alkaloids have been the subject of numerous synthetic studies over the last forty years, very few enantioselective total syntheses have been reported to date.

As part of a research program directed towards the synthesis of various alkaloid natural products, we have developed a unified strategy for the preparation of the hasubanan and acutumine alkaloids. Specifically, a highly diastereoselective 1,2-addition of organometallic reagents to benzoquinone-derived tert-butanesulfinimines was established, which provides access to enantioenriched 4-aminocyclohexadienone products. This methodology enabled the enantioselective construction of functionalized dihydroindolones, which were found to undergo intramolecular Friedel-Crafts conjugate additions to furnish the propellane cores of several hasubanan alkaloids. As a result of these studies, the first enantioselective total syntheses of 8-demethoxyrunanine and cepharatines A, C, and D were accomplished in 9-11 steps from commercially available starting materials.

More recent efforts have focused on applying the sulfinimine methodology to the synthesis of a more structurally complex propellane alkaloid, acutumine. Extensive studies have determined that a properly functionalized dihydroindolone undergoes a photochemical [2+2] cycloaddition followed by a lactone fragmentation/Dieckmann cyclization to establish the carbocyclic framework of the natural product. The preparation of more appropriately oxidized propellane intermediates is currently under investigation, and is anticipated to facilitate our synthetic endeavors toward acutumine.

Investigations of the origin of stereocontrol in syndiospecific Ziegler-Natta polymerizations

In order to expand our understanding of the mechanism of stereocontrol in syndiospecific α-olefin polymerization, a family of Cs-symmetric, ansa-group 3 metallocenes was targeted as polymerization catalysts. The syntheses of new ansa-yttrocene and scandocene derivatives that employ the doubly [SiMe₂]- bridged ligand array (1,2-SiMe₂)₂{C₅H-3,5-(CHMe₂)₂} (where R = t- butyl, tBuThp; where R = i-propyl, iPrThp) are described. The structures of tBuThpY(µ-Cl)₂K(THF)₂, tBuThpSc(µ-Cl)₂K(Et₂O)₂, tBuThpYCH(SiMe₃)₂, Y₂{µ₂-(tBuThp)₂}(µ₂-H)₂, and tBuThpSc(µ-CH₃)₂ have been examined by single crystal X-ray diffraction methods. Ansa-yttrocenes and scandocenes that incorporate the singly [CPh₂]-bridged ligand array (CPh₂)(C₅H₄)(C₁₃H₈)(where C₅H₄ = Cp, cyclopentadienyl; where C₁₃H₈ = Flu, fluourenyl) have also been prepared. Select meallocene alkyl complexes are active single component catalysts for homopolymerization of propylene and 1-pentene. The scandocene tetramethylaluminate complexes generate polymers with the highes molecular weights of the series. Under all conditions examined atactic polymer microstructures are observed, suggesting a chain-end mechanism for stereocontrol.

A series of ansa-tantalocenes have been prepared as models for Ziegler-Natta polymerization catalysts. A singly bridged ansa-tantalocene trimethyl complex, Me₂Si(η⁵-C₅H₄)₂TaMe₃, has been prepared and used for the synthesis of a tantalocene ethylene-methyl complex. Addition of propylene to this ethylene-methyl adduct results in olefin exchange to give a mixture of endo and exo propylene isomers. Doubly-silylene bridged ansa-tantalocene complexes have been prepared with the tBuThp ligand; a tantalocene trimethyl complex and a tantalocene methylidene-methyl complex have been synthesized and characterized by X-ray diffraction. Thermolysis of the methylidene-methyl complex affords the corresponding ethylene-hydride complex. Addition of either propylene or styrene to this ethylene-hydride compound results in olefin exchange. In both cases, only one product isomer is observed. Studies of olefin exchange with ansa-tantalocene olefin-hydride and olefin-methyl complexes have provided information about the important steric influences for olefin coordination in Ziegler-Natta polymerization.

Ancillary ligand effects : from zirconium(IV)-catalyzed homogeneous propylene polymerization to platinum(II)-mediated C-H bond activation

A series of C_s- and C₁-symmetric doubly-linked ansa-metallocenes of the general formula {1,1'-SiMe₂-2,2'-E-('ƞ⁵-C₅H₂-4-R¹)-(ƞ⁵-C₅H-3',5'-(CHMe₂)₂)}ZrC₂ (E = SiMe₂ (1), SiPh₂ (2), SiMe₂ -SiMe₂ (3); R¹ = H, CHMe₂, C₅H₉, C₆H₁₁, C₆H₅) has been prepared. When activated by methylaluminoxane, these are active propylene polymerization catalysts. 1 and 2 produce syndiotactic polypropylenes, and 3 produces isotactic polypropylenes. Site epimerization is the major pathway for stereoerror formation for 1 and 2. In addition, the polymer chain has slightly stronger steric interaction with the diphenylsilylene linker than with the dimethylsilylene linker. This results in more frequent site epimerization and reduced syndiospecificity for 2 compared to 1.

C₁-Symmetric ansa-zirconocenes [1,1 '-SiMe₂-(C₅H₄)-(3-R-C₅H₃)]ZrCl₂ (4), [1,1 '-SiMe₂-(C₅H₄)-(2,4-R₂-C₅H₂)]ZrCl₂ (5) and [1,1 '-SiMe₂-2,2 '-(SiMe₂-SiMe₂)-(C₅H₃)-( 4-R-C₅H₂)]ZrCl₂ (6) have been prepared to probe the origin of isospecificity in 3. While 4 and 3 produce polymers with similar isospecificity, 5 and 6 give mostly hemi-isotactic-like polymers. It is proposed that the facile site epimerization via an associative pathway allows rapid equilibration of the polymer chain between the isospecific and aspecific insertion sites. This results in more frequent insertion from the isospecific site, which has a lower kinetic barrier for chain propagation. On the other hand, site epimerization for 5 and 6 is slow. This leads to mostly alternating insertion from the isospecific and aspecific sites, and consequently, a hemi-isotactic-like polymers. In comparison, site epimerization is even slower for 3, but enchainment from the aspecific site has an extremely high kinetic barrier for monomer coordination. Therefore, enchainment occurs preferentially from the isospecific site to produce isotactic polymers.

A series of cationic complexes [(ArN=CR-CR=NAr)PtMe(L)]⁺[BF₄]⁺ (Ar = aryl; R = H, CH₃; L = water, trifluoroethanol) has been prepared. They react smoothly with benzene at approximately room temperature in trifluoroethanol solvent to yield methane and the corresponding phenyl Pt(II) cations, via Pt(IV)-methyl-phenyl-hydride intermediates. The reaction products of methyl-substituted benzenes suggest an inherent reactivity preference for aromatic over benzylic C-H bond activation, which can however be overridden by steric effects. For the reaction of benzene with cationic Pt(II) complexes, in which the diimine ligands bear 3,5-disubstituted aryl groups at the nitrogen atoms, the rate-determining step is C-H bond activation. For the more sterically crowded analogs with 2,6-dimethyl-substituted aryl groups, benzene coordination becomes rate-determining. The more electron-rich the ligand, as reflected by the CO stretching frequency in the IR spectrum of the corresponding cationic carbonyl complex, the faster the rate of C-H bond activation. This finding, however, does not reflect the actual C-H bond activation process, but rather reflects only the relative ease of solvent molecules displacing water molecules to initiate the reaction. That is, the change in rates is mostly due to a ground state effect. Several lines of evidence suggest that associative substitution pathways operate to get the hydrocarbon substrate into, and out of, the coordination sphere; i.e., that benzene substitution proceeds by a solvent- (TFE-) assisted associative pathway.

The development of an asymmetric palladium-catalyzed conjugate addition and its application toward the total syntheses of taiwaniaquinoid natural products

The asymmetric synthesis of quaternary stereocenters remains a challenging problem in organic synthesis. Past work from the Stoltz laboratory has resulted in methodology to install quaternary stereocenters α- or γ- to carbonyl compounds. Thus, the asymmetric synthesis of β-quaternary stereocenters was a desirable objective, and was accomplished by engineering the palladium-catalyzed addition of arylmetal organometallic reagents to α,β-unsaturated conjugate acceptors.

Herein, we described the rational design of a palladium-catalyzed conjugate addition reactions utilizing a catalyst derived from palladium(II) trifluoroacetate and pyridinooxazole ligands. This reaction is highly tolerant of protic solvents and oxygen atmosphere, making it a practical and operationally simple reaction. The mild conditions facilitate a remarkably high functional group tolerance, including carbonyls, halogens, and fluorinated functional groups. Furthermore, the reaction catalyzed conjugate additions with high enantioselectivity with conjugate acceptors of 5-, 6-, and 7-membered ring sizes. Extension of the methodology toward the asymmetric synthesis of flavanone products is presented, as well.

A computational and experimental investigation into the reaction mechanism provided a stereochemical model for enantioinduction, whereby the α-methylene protons adjacent the enone carbonyl clashes with the tert-butyl groups of the chiral ligand. Additionally, it was found that the addition of water and ammonium hexafluorophosphate significantly increases the reaction rate without sacrificing enantioselectivity. The synergistic effects of these additives allowed for the reaction to proceed at a lower temperature, and thus facilitated expansion of the substrate scope to sensitive functional groups such as protic amides and aryl bromides. Investigations into a scale-up synthesis of the chiral ligand (S)-tert-butylPyOx are also presented. This three-step synthetic route allowed for synthesis of the target compound of greater than 10 g scale.

Finally, the application of the newly developed conjugate addition reaction toward the synthesis of the taiwaniaquinoid class of terpenoid natural products is discussed. The conjugate addition reaction formed the key benzylic quaternary stereocenter in high enantioselectivity, joining together the majority of the carbons in the taiwaniaquinoid scaffold. Efforts toward the synthesis of the B-ring are presented.

Formal methods for control synthesis in partially observed environments : application to autonomous robotic manipulation

Modern robots are increasingly expected to function in uncertain and dynamically challenging environments, often in proximity with humans. In addition, wide scale adoption of robots requires on-the-fly adaptability of software for diverse application. These requirements strongly suggest the need to adopt formal representations of high level goals and safety specifications, especially as temporal logic formulas. This approach allows for the use of formal verification techniques for controller synthesis that can give guarantees for safety and performance. Robots operating in unstructured environments also face limited sensing capability. Correctly inferring a robot's progress toward high level goal can be challenging.

This thesis develops new algorithms for synthesizing discrete controllers in partially known environments under specifications represented as linear temporal logic (LTL) formulas. It is inspired by recent developments in finite abstraction techniques for hybrid systems and motion planning problems. The robot and its environment is assumed to have a finite abstraction as a Partially Observable Markov Decision Process (POMDP), which is a powerful model class capable of representing a wide variety of problems. However, synthesizing controllers that satisfy LTL goals over POMDPs is a challenging problem which has received only limited attention.

This thesis proposes tractable, approximate algorithms for the control synthesis problem using Finite State Controllers (FSCs). The use of FSCs to control finite POMDPs allows for the closed system to be analyzed as finite global Markov chain. The thesis explicitly shows how transient and steady state behavior of the global Markov chains can be related to two different criteria with respect to satisfaction of LTL formulas. First, the maximization of the probability of LTL satisfaction is related to an optimization problem over a parametrization of the FSC. Analytic computation of gradients are derived which allows the use of first order optimization techniques.

The second criterion encourages rapid and frequent visits to a restricted set of states over infinite executions. It is formulated as a constrained optimization problem with a discounted long term reward objective by the novel utilization of a fundamental equation for Markov chains - the Poisson equation. A new constrained policy iteration technique is proposed to solve the resulting dynamic program, which also provides a way to escape local maxima.

The algorithms proposed in the thesis are applied to the task planning and execution challenges faced during the DARPA Autonomous Robotic Manipulation - Software challenge.