Unnatural amino acids in the synthesis and semisynthesis of metalloprotein motifs

The design, synthesis, and characterization of two novel metalloprotein motifs is presented. The first project involved the design and construction of a protein motif which was programmed to form a tetradentate metal complex upon the addition of metal cations. The overall structure of the motif was based on a ββ super-secondary structure consisting of a flexible peptide sequence flanked by metal binding regions located at the carboxy and amino termini. The metal binding region near the amino terminus was constructed from a reverse turn motif with two metal ligating residues, (2R, 3R)-β-methyl-cysteine and histidine. Selection of the peptide sequence for this region was based on the conformational analysis of a series of tetrapeptides designed to form reverse turns in solution.

The stereospecific syntheses of a series of novel bipyridyl- and phenanthrolylsubstituted amino acids was carried out to provide ligands for the carboxy terminus metal binding region. These residues were incorporated into peptide sequences using solid phase peptide synthesis protocols, and metal binding studies indicated that the metal binding properties of these ligands was dictated by the specific regioisomer of the heteroaromatic ring and the peptide primary sequence.

Finally, a peptide containing optimized components for the metal binding regions was prepared to test the ability of the compound to form the desired intramolecular peptide:metal cation complexes. Metal binding studies demonstrated that the peptide formed monomeric complexes with very high metal cation binding affinities and that the two metal binding regions act cooperatively in the metal binding process. The use of these systems in the design of proteins capable of regulating naturally occurring proteins is discussed.

The second project involved the semisynthesis of two horse heart cytochrome c mutants incorporating the bipyridyl-amino acids at position 72 of the protein sequence. Structural studies on the proteins indicated that the bipyridyl amino acids had a neglible effect on the protein structure. One of the mutants was modified with Ru(bpy)_2^(+2) to form a redox-active protein, and the modified protein was found to have enhanced electron transfer properties between the heme and the introduced metal site.

Structural and mechanistic motifs in membrane proteins: the three-dimensional modelling of rhodopsin, band 3, and the nicotinic acetylcholine receptor

Because so little is known about the structure of membrane proteins, an attempt has been made in this work to develop techniques by which to model them in three dimensions. The procedures devised rely heavily upon the availability of several sequences of a given protein. The modelling procedure is composed of two parts. The first identifies transmembrane regions within the protein sequence on the basis of hydrophobicity, β-turn potential, and the presence of certain amino acid types, specifically, proline and basic residues. The second part of the procedure arranges these transmembrane helices within the bilayer based upon the evolutionary conservation of their residues. Conserved residues are oriented toward other helices and variable residues are positioned to face the surrounding lipids. Available structural information concerning the protein's helical arrangement, including the lengths of interhelical loops, is also taken into account. Rhodopsin, band 3, and the nicotinic acetylcholine receptor have all been modelled using this methodology, and mechanisms of action could be proposed based upon the resulting structures.

Specific residues in the rhodopsin and iodopsin sequences were identified, which may regulate the proteins' wavelength selectivities. A hinge-like motion of helices M3, M4, and M5 with respect to the rest of the protein was proposed to result in the activation of transducin, the G-protein associated with rhodopsin. A similar mechanism is also proposed for signal transduction by the muscarinic acetylcholine and β-adrenergic receptors.

The nicotinic acetylcholine receptor was modelled with four trans-membrane helices per subunit and with the five homologous M2 helices forming the cation channel. Putative channel-lining residues were identified and a mechanism of channel-opening based upon the concerted, tangential rotation of the M2 helices was proposed.

Band 3, the anion exchange protein found in the erythrocyte membrane, was modelled with 14 transmembrane helices. In general the pathway of anion transport can be viewed as a channel composed of six helices that contains a single hydrophobic restriction. This hydrophobic region will not allow the passage of charged species, unless they are part of an ion-pair. An arginine residue located near this restriction is proposed to be responsible for anion transport. When ion-paired with a transportable anion it rotates across the barrier and releases the anion on the other side of the membrane. A similar process returns it to its original position. This proposed mechanism, based on the three-dimensional model, can account for the passive, electroneutral, anion exchange observed for band 3. Dianions can be transported through a similar mechanism with the additional participation of a histidine residue. Both residues are located on M10.

A biochemical and structural characterization of Drosophila neuroglian

A number of cell-cell interactions in the nervous system are mediated by immunoglobulin gene superfamily members. For example, neuroglian, a homophilic neural cell adhesion molecule in Drosophila, has an extracellular portion comprising six C- 2 type immunoglobulin-like domains followed by five fibronectin type III (FnIII) repeats. Neuroglian shares this domain organization and significant sequence identity with Ll, a murine neural adhesion molecule that could be a functional homologue. Here I report the crystal structure of a proteolytic fragment containing the first two FnIII repeats of neuroglian (NgFn 1,2) at 2.0Å. The interpretation of photomicrographs of rotary shadowed Ng, the entire extracellular portion of neuroglian, and NgFnl-5, the five neuroglian Fn III domains, is also discussed.

The structure of NgFn 1,2 consists of two roughly cylindrical β-barrel structural motifs arranged in a head-to-tail fashion with the domains meeting at an angle of ~120, as defined by the cylinder axes. The folding topology of each domain is identical to that previously observed for single FnIII domains from tenascin and fibronectin. The domains of NgFn1,2 are related by an approximate two fold screw axis that is nearly parallel to the longest dimension of the fragment. Assuming this relative orientation is a general property of tandem FnIII repeats, the multiple tandem FnIII domains in neuroglian and other proteins are modeled as thin straight rods with two domain zig-zag repeats. When combined with the dimensions of pairs of tandem immunoglobulin-like domains from CD4 and CD2, this model suggests that neuroglian is a long narrow molecule (20 - 30 Å in diameter) that extends up to 370Å from the cell surface.

In photomicrographs, rotary shadowed Ng and NgFn1-5 appear to be highly flexible rod-like molecules. NgFn 1-5 is observed to bend in at least two positions and has a mean total length consistent with models generated from the NgFn1,2 structure. Ng molecules have up to four bends and a mean total length of 392 Å, consistent with a head-to-tail packing of neuroglian's C2-type domains.

Enantioselective synthesis of pyrroloindolines and tryptophan derivatives by an asymmetric protonation reaction

Pyrroloindoline and unnatural tryptophan motifs are important targets for synthesis based on their incorporation into a diverse array of biologically active natural products. Both types of alkaloids have also found applications as chiral catalysts and tryptophan derivatives are commonly employed as biological probes. On account of their applications, these frameworks have inspired the development of numerous enantioselective, catalytic reactions. In particular, the past few years have witnessed an impressive number of novel approaches for pyrroloindoline formation.

The first project described herein involves our contribution to pyrroloindoline research. We have developed an (R)-BINOL•SnCl₄-catalyzed formal (3 + 2) cycloaddition reaction between 3-substituted indoles and 2-amidoacrylates that affords pyrroloindoline-2-carboxylates bearing an all-carbon quaternary center. Mechanistic studies have elucidated that the enantiodetermining step is a highly face-selective catalyst-controlled protonation reaction. The subsequent application of this asymmetric protonation strategy to the synthesis of a variety of enantioenriched tryptophan derivatives is also discussed.

Sterochemically modified polyamides for recognition in the minor groove of DNA

The design of synthetic molecules that recognize specific sequences of DNA is an ongoing challenge in molecular medicine. Cell-permeable small molecules targeting predetermined DNA sequences offer a potential approach for offsetting the abnormal effects of misregulated gene-expression. Over the past twenty years, Professor Peter B. Dervan has developed a set of pairing rules for the rational design of minor groove binding polyamides containing pyrrole (Py), imidazole (Im), and hydroxypyrrole (Hp). Polyamides have illustrated the capability to permeate cells and inhibit transcription of specific genes in vivo. This provides impetus to identify structural elements that expand the repetoire of polyamide motifs with recognition properties comparable to naturally occurring DNA binding proteins. Through the introduction of chiral amino acids, we have developed chiral polyamides with stereochemically regulated binding characteristics. In addition, chiral substituents have facilitated the development of new polyamide motifs that broaden binding site sizes targetable by this class of ligands.

Tunable heparan sulfate glycomimetics for modulating chemokine activity

Heparan sulfate (HS) glycosaminoglycans participate in critical biological processes by modulating the activity of a diverse set of protein binding partners. Such proteins include all known members of the chemokine superfamily, which are thought to guide the migration of distinct subsets of immune cells through their interactions with HS proteoglycans on endothelial cell surfaces. Animal-derived heparin polysaccharides have been shown to reduce inflammation levels through the inhibition of HS-chemokine interactions; however, the clinical usage of heparin as an anti-inflammatory drug is hampered by its anticoagulant activity and potential risk for side effects, such as heparin-induced thrombocytopenia (HIT).

Here, we describe an expedient, divergent synthesis to prepare defined glycomimetics of HS that recapitulate the macromolecular structure and biological activity of natural HS glycosaminoglycans. Our synthetic approach uses a core disaccharide precursor to generate a library of four differentially sulfated polymers. We show that a trisulfated glycopolymer antagonizes the chemotactic activities of pro-inflammatory chemokine RANTES with similar potency as heparin polysaccharide, without potentiating the anticoagulant activities of antithrombin III. The same glycopolymer also inhibited the homeostatic chemokine SDF-1 with significantly more efficacy than heparin. Our work offers a general strategy for modulating chemokines and dissecting the pleiotropic functions of HS/heparin through the presentation of defined sulfation motifs within multivalent polymeric scaffolds.

Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

Structure-, stereochemistry-, and metal-regulated DNA binding/cleaving molecules

In the cell, the binding of proteins to specific sequences of double helical DNA is essential for controlling the processes of protein synthesis (at the level of DNA transcription) and cell proliferation (at the level of DNA replication). In the laboratory, the sequence-specific DNA binding/cleaving properties of restriction endonuclease enzymes (secreted by microorganisms to protect them from foreign DNA molecules) have helped to fuel a revolution in molecular biology. The strength and specificity of a protein:DNA interaction depend upon structural features inherent to the protein and DNA sequences, but it is now appreciated that these features (and therefore protein:DNA complexation) may be altered (regulated) by other protein:DNA complexes, or by environmental factors such as temperature or the presence of specific organic molecules or inorganic ions. It is also now appreciated that molecules much smaller than proteins (including antibiotics of molecular weight less than 2000 and oligonucleotides) can bind to double-helical DNA in sequence-specific fashion. Elucidation of structural motifs and microscopic interactions responsible for the specific molecular recognition of DNA leads to greater understanding of natural processes and provides a basis for the design of novel sequence-specific DNA binding molecules. This thesis describes the synthesis and DNA binding/cleaving characteristics of molecules designed to probe structural, stereochemical, and environmental factors that regulate sequence-specific DNA recognition.

Chapter One introduces the DNA minor groove binding antibiotics Netropsin and Distamycin A, which are di- and tri(N-methylpyrrolecarboxamide) peptides, respectively. The method of DNA affinity cleaving, which has been employed to determine DNA binding properties of designed synthetic molecules is described. The design and synthesis of a series of Netropsin dimers linked in tail-to-tail fashion (by oxalic, malonic, succinic, or fumaric acid), or in head-to-tail fashion (by glycine, β-alanine, and γ-aminobutanoic acid (Gaba)) are presented. These Bis(Netropsin)s were appended with the iron-chelating functionality EDTA in order to make use of the technique of DNA affinity cleaving. Bis(Netropsin)-EDTA compounds are analogs of penta(N-methylpyrrolecarboxamide)-EDTA (P5E), which may be considered a head-to-tail Netropsin dimer linked by Nmethylpyrrolecarboxamide. Low- and high-resolution analysis of pBR322 DNA affinity cleaving by the iron complexes of these molecules indicated that small changes in the length and nature of the linker had significant effects on DNA binding/cleaving efficiency (a measure of DNA binding affinity). DNA binding/cleaving efficiency was found to decrease with changes in the linker in the order β-alanine > succinamide > fumaramide > N-methylpyrrolecarboxamide > malonamide >glycine, γ-aminobutanamide > oxalamide. In general, the Bis(Netropsin)-EDTA:Fe compounds retained the specificity for seven contiguous A:T base pairs characteristic of P5E:Fe binding. However, Bis(Netropsin)Oxalamide- EDTA:Fe exhibited decreased specificity for A:T base pairs, and Bis(Netropsin)-Gaba-EDT A:Fe exhibited some DNA binding sites of less than seven base pairs. Bis(Netropsin)s linked with diacids have C₂-symmmetrical DNA binding subunits and exhibited little DNA binding orientation preference. Bis(Netropsin)s linked with amino acids lack C₂-symmetrical DNA binding subunits and exhibited higher orientation preferences. A model for the high DNA binding orientation preferences observed with head-to-tail DNA minor groove binding molecules is presented.

Chapter Two describes the design, synthesis, and DNA binding properties of a series of chiral molecules: Bis(Netropsin)-EDTA compounds with linkers derived from (R,R)-, (S,S)-, and (RS,SR)-tartaric acids, (R,R)-, (S,S)-, and (RS,SR)-tartaric acid acetonides, (R)- and (S)-malic acids, N ,N-dimethylaminoaspartic acid, and (R)- and (S)-alanine, as well as three constitutional isomers in which an N-methylpyrrolecarboxamide (P1) subunit and a tri(N-methylpyrrolecarboxamide)-EDTA (P3-EDTA) subunit were linked by succinic acid, (R ,R)-, and (S ,S)-tartaric acids. DNA binding/cleaving efficiencies among this series of molecules and the Bis(Netropsin)s described in Chapter One were found to decrease with changes in the linker in the order β-alanine > succinamide > P1-succinamide-P3 > fumaramide > (S)-malicamide > N-methylpyrrolecarboxamide > (R)-malicamide > malonamide > N ,N-dimethylaminoaspanamide > glycine = Gaba = (S,S)-tartaramide = P1-(S,S)-tanaramide-P3 > oxalamide > (RS,SR)-tartaramide = P1- (R,R)-tanaramide-P3 > (R,R)-tartaramide (no sequence-specific DNA binding was detected for Bis(Netropsin)s linked by (R)- or (S)-alanine or by tartaric acid acetonides). The chiral molecules retained DNA binding specificity for seven contiguous A:T base pairs. From the DNA affinity cleaving data it could be determined that: 1) Addition of one or two substituents to the linker of Bis(Netropsin)-Succinamide resulted in stepwise decreases in DNA binding affinity; 2) molecules with single hydroxyl substituents bound DNA more strongly than molecules with single dimethylamino substituents; 3) hydroxyl-substituted molecules of (S) configuration bound more strongly to DNA than molecules of (R) configuration. This stereochemical regulation of DNA binding is proposed to arise from the inherent right-handed twist of (S)-enantiomeric Bis(Netropsin)s versus the inherent lefthanded twist of (R)-enantiomeric Bis(Netropsin)s, which makes the (S)-enantiomers more complementary to the right-handed twist of B form DNA.

Chapter Three describes the design and synthesis of molecules for the study of metalloregulated DNA binding phenomena. Among a series of Bis(Netropsin)-EDTA compounds linked by homologous tethers bearing four, five, or six oxygen atoms, the Bis(Netropsin) linked by a pentaether tether exhibited strongly enhanced DNA binding/cleaving in the presence of strontium or barium cations. The observed metallospecificity was consistent with the known affinities of metal cations for the cyclic hexaether 18-crown-6 in water. High-resolution DNA affinity cleaving analysis indicated that DNA binding by this molecule in the presence of strontium or barium was not only stronger but of different sequence-specificity than the (weak) binding observed in the absence of metal cations. The metalloregulated binding sites were consistent with A:T binding by the Netropsin subunits and G:C binding by a strontium or barium:pentaether complex. A model for the observed positive metalloregulation and novel sequence-specificity is presented. The effects of 44 different cations on DNA affinity cleaving by P5E:Fe were examined. A series of Bis(Netropsin)-EDTA compounds linked by tethers bearing two, three, four, or five amino groups was also synthesized. These molecules exhibited strong and specific binding to A:T rich regions of DNA. It was found that the iron complexes of these molecules bound and cleaved DNA most efficiently at pH 6.0-6.5, while P5E:Fe bound and cleaved most efficiently at pH 7.5-8.0. Incubating the Bis(Netropsin) Polyamine-EDTA:Fe molecules with K₂PdCl₄ abolished their DNA binding/cleaving activity. It is proposed that the observed negative metalloregulation arises from kinetically inert Bis(Netropsin) Polyamine:Pd(II) complexes or aggregates, which are sterically unsuitable for DNA complexation. Finally, attempts to produce a synthetic metalloregulated DNA binding protein are described. For this study, five derivatives of a synthetic 52 amino acid residue DNA binding/cleaving protein were produced. The synthetic mutant proteins carried a novel pentaether ionophoric amino acid residue at different positions within the primary sequence. The proteins did not exhibit significant DNA binding/cleaving activity, but they served to illustrate the potential for introducing novel amino acid residues within DNA binding protein sequences, and for the development of the tricyclohexyl ester of EDTA as a superior reagent for the introduction of EDT A into synthetic proteins.

Chapter Four describes the discovery and characterization of a new DNA binding/cleaving agent, [SalenMn(III)]OAc. This metal complex produces single- and double-strand cleavage of DNA, with specificity for A:T rich regions, in the presence of oxygen atom donors such as iodosyl benzene, hydrogen peroxide, or peracids. Maximal cleavage by [SalenMn(III)]OAc was produced at pH 6-7. A comparison of DNA singleand double-strand cleavage by [SalenMn(III)]⁺ and other small molecules (Methidiumpropyl-EDTA:Fe, Distamycin-EDTA:Fe, Neocarzinostatin, Bleomycin:Fe) is presented. It was found that DNA cleavage by [SalenMn(III)]⁺ did not require the presence of dioxygen, and that base treatment of DNA subsequent to cleavage by [SalenMn(III)]⁺ afforded greater cleavage and alterations in the cleavage patterns. Analysis of DNA products formed upon DNA cleavage by [SalenMn(III)] indicated that cleavage was due to oxidation of the sugar-phosphate backbone of DNA. Several mechanisms consistent with the observed products and reaction requirements are discussed.

Chapter Five describes progress on some additional studies. In one study, the DNA binding/cleaving specificities of Distamycin-EDTA derivatives bearing pyrrole N-isopropyl substituents were found to be the same as those of derivatives bearing pyrrole N-methyl substituents. In a second study, the design of and synthetic progress towards a series of nucleopeptide activators of transcription are presented. Five synthetic plasmids designed to test for activation of in vitro run-off transcription by DNA triple helix-forming oligonucleotides or nucleopeptides are described.

Chapter Six contains the experimental documentation of the thesis work.

I. The crystal structure of trimesic acid. II. Topics in crystallographic calculations

I. Trimesic acid (1, 3, 5-benzenetricarboxylic acid) crystallizes with a monoclinic unit cell of dimensions a = 26.52 A, b = 16.42 A, c = 26.55 A, and β = 91.53° with 48 molecules /unit cell. Extinctions indicated a space group of Cc or C2/c; a satisfactory structure was obtained in the latter with 6 molecules/asymmetric unit - C₅₄O₃₆H₃₆ with a formula weight of 1261 g. Of approximately 12,000 independent reflections within the CuK_α sphere, intensities of 11,563 were recorded visually from equi-inclination Weissenberg photographs.

The structure was solved by packing considerations aided by molecular transforms and two- and three-dimensional Patterson functions. Hydrogen positions were found on difference maps. A total of 978 parameters were refined by least squares; these included hydrogen parameters and anisotropic temperature factors for the C and O atoms. The final R factor was 0.0675; the final "goodness of fit" was 1.49. All calculations were carried out on the Caltech IBM 7040-7094 computer using the CRYRM Crystallographic Computing System.

The six independent molecules fall into two groups of three nearly parallel molecules. All molecules are connected by carboxylto- carboxyl hydrogen bond pairs to form a continuous array of sixmolecule rings with a chicken-wire appearance. These arrays bend to assume two orientations, forming pleated sheets. Arrays in different orientations interpenetrate - three molecules in one orientation passing through the holes of three parallel arrays in the alternate orientation - to produce a completely interlocking network. One third of the carboxyl hydrogen atoms were found to be disordered.

II. Optical transforms as related to x-ray diffraction patterns are discussed with reference to the theory of Fraunhofer diffraction.

The use of a systems approach in crystallographic computing is discussed with special emphasis on the way in which this has been done at the California Institute of Technology.

An efficient manner of calculating Fourier and Patterson maps on a digital computer is presented. Expressions for the calculation of to-scale maps for standard sections and for general-plane sections are developed; space-group-specific expressions in a form suitable for computers are given for all space groups except the hexagonal ones.

Expressions for the calculation of settings for an Eulerian-cradle diffractometer are developed for both the general triclinic case and the orthogonal case.

Photographic materials on pp. 4, 6, 10, and 20 are essential and will not reproduce clearly on Xerox copies. Photographic copies should be ordered.

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.

Heterometallic Complexes as Models of Enzymatic Active Sites

This dissertation describes studies on two multinucleating ligand architectures: the first scaffold was designed to support tricopper complexes, while the second platform was developed to support tri- and tetrametallic clusters.

In Chapter 2, the synthesis of yttrium (and lanthanide) complexes supported by a tripodal ligand framework designed to bind three copper centers in close proximity is described. Tricopper complexes were shown to react with dioxygen in a 1:1 [Cu₃]/O₂ stoichiometry to form intermediates in which the O–O bond was fully cleaved, as characterized via UV-Vis spectroscopy and determination of the reaction stoichiometry. Pre-arrangement of the three Cu centers was pivotal to cooperative O₂ activation, as mono-copper complexes reacted differently with dioxgyen. The reactivity of the observed intermediates was studied with various substrates (reductants, O-atom acceptors, H-atom donors, Brønsted acids) to determine their properties. In Chapter 3, the reactivity of the same yttrium-tricopper complex with nitric oxide was explored. Reductive coupling to form a trans-hyponitrite complex (characterized by X-ray crystallography) was observed via cooperative reactivity by an yttrium and a copper center on two distinct tetrametallic units. The hyponitrite complex was observed to release nitrous oxide upon treatment with a Brønsted acid, supporting its viability as an intermediate in nitric oxide reduction to nitrous oxide.

In Chapter 4, a different multinucleating ligand scaffold was employed to synthesize heterometallic triiron clusters containing one oxide and one hydroxide bridges. The effects of the redox-inactive, Lewis acidic heterometals on redox potential was studied by cyclic voltammetry, unveiling a linear correlation between redox potential and heterometal Lewis acidity. Further studies on these complexes showed that the Lewis acidity of the redox-inactive metals also affected the oxygen-atom transfer reactivity of these clusters. Comparisons of this reactivity with manganese systems, collaborative efforts to reassign the structures of related manganese oxo-hydroxo clusters, and synthetic attempts to access related dioxo clusters are also described.

In Appendix A, ongoing efforts to synthesize new clusters supported by the same multinucleating ligand platform are described. Studies of novel approaches towards ligand exchange in tetrametallic clusters and incorporation of new supporting and bridging ligand motifs in trinuclear complexes are presented.