Synthesis of (±)-7,8-epoxy-4-basmen-6-one by a transannular cyclization strategy

The first synthesis of the cembranoid natural product (±)-7,8-epoxy-4-basmen-6- one (1) is described. Key steps of the synthetic route include the cationic cyclization of the acid chloride from 15 to provide the macrocycle 16, and the photochemical transannular radical cyclization of the ester 41 to form the tricyclic product 50. Product 50 was transformed into 1 in ten steps. Transition-state molecular modeling studies were found to provide accurate predictions of the structural and stereochemical outcomes of cyclization reactions explored experimentally in the development of the synthetic route to 1. These investigations should prove valuable in the development of transannular cyclization as a strategy for synthetic simplification.

Role of feedback and dynamics in a gene regulatory network

Cells exhibit a diverse repertoire of dynamic behaviors. These dynamic functions are implemented by circuits of interacting biomolecules. Although these regulatory networks function deterministically by executing specific programs in response to extracellular signals, molecular interactions are inherently governed by stochastic fluctuations. This molecular noise can manifest as cell-to-cell phenotypic heterogeneity in a well-mixed environment. Single-cell variability may seem like a design flaw but the coexistence of diverse phenotypes in an isogenic population of cells can also serve a biological function by increasing the probability of survival of individual cells upon an abrupt change in environmental conditions. Decades of extensive molecular and biochemical characterization have revealed the connectivity and mechanisms that constitute regulatory networks. We are now confronted with the challenge of integrating this information to link the structure of these circuits to systems-level properties such as cellular decision making. To investigate cellular decision-making, we used the well studied galactose gene-regulatory network in \textit{Saccharomyces cerevisiae}. We analyzed the mechanism and dynamics of the coexistence of two stable on and off states for pathway activity. We demonstrate that this bimodality in the pathway activity originates from two positive feedback loops that trigger bistability in the network. By measuring the dynamics of single-cells in a mixed sugar environment, we observe that the bimodality in gene expression is a transient phenomenon. Our experiments indicate that early pathway activation in a cohort of cells prior to galactose metabolism can accelerate galactose consumption and provide a transient increase in growth rate. Together these results provide important insights into strategies implemented by cells that may have been evolutionary advantageous in competitive environments.

Revealing, illuminating, and modifying proteins in human diseases using noncanonical amino acids

To better understand human diseases, much recent work has focused on proteins to either identify disease targets through proteomics or produce therapeutics via protein engineering. Noncanonical amino acids (ncAAs) are tools for altering the chemical and physical properties of proteins, providing a facile strategy not only to label proteins but also to engineer proteins with novel properties. My thesis research has focused on the development and applications of noncanonical amino acids in identifying, imaging, and engineering proteins for studying human diseases. Chapter 1 introduces the concept of ncAAs and reveals insights to how I chose my thesis projects.

ncAAs have been incorporated to tag and enrich newly synthesized proteins for mass spectrometry through a method termed BONCAT, or bioorthogonal noncanonical amino acid tagging. Chapter 2 describes the investigation of the proteomic response of human breast cancer cells to induced expression of tumor suppressor microRNA miR-126 by combining BONCAT with another proteomic method, SILAC or stable isotope labeling by amino acids in cell culture. This proteomic analysis led to the discovery of a direct target of miR-126, shedding new light on its role in suppressing cancer metastasis.

In addition to mass spectrometry, ncAAs can also be utilized to fluorescently label proteins. Chapter 3 details the synthesis of a set of cell-permeant cyclooctyne probes and demonstration of selective labeling of newly synthesized proteins in live mammalian cells using azidohomoalanine. Similar to live cell imaging, the ability to selectively label a particular cell type within a mixed cell population is important to interrogating many biological systems, such as tumor microenvironments. By taking advantage of the metabolic differences between cancer and normal cells, Chapter 5 discusses efforts to develop selective labeling of cancer cells using a glutamine analogue.

Furthermore, Chapter 4 describes the first demonstration of global replacement at polar amino acid positions and its application in developing an alternative PEGylation strategy for therapeutic proteins. Polar amino acids typically occupy solvent-exposed positions on the protein surface, and incorporation of noncanonical amino acids at these positions should allow easier modification and cause less perturbation compared to replacements at the interior positions of proteins.

A unified strategy to ent-kauranoid natural products : total syntheses of (–)-maoecrystal Z, (–)-trichorabdal A, and (–)-longikaurin E

The diterpenoid constituents of the Isodon plants have attracted reasearchers interested in both their chemical structures and biological properties for more than a half-century. In recent years, the isolations of new members displaying previously unprecedented ring systems and highly selective biological properties have piqued interest from the synthetic community in this class of natural products.

Reported herein is the first total synthesis of such a recently isolated diterpenoid, (–)-maoecrystal Z. The principal transformations implemented in this synthesis include two highly diastereoselective radical cyclization reactions: a Sm^(II)-mediated reductive cascade cyclization, which forms two rings and establishes four new stereocenters in a single step, and a Ti^(III)-mediated reductive epoxide-acrylate coupling that yields a functionalized spirolactone product, which forms a core bicycle of maoecrystal Z.

The preparation of two additional ent-kauranoid natural products, (–)-trichorabdal A and (–)-longikaurin E, is also described from a derivative of this key spirolactone. These syntheses are additionally enabled by the palladium-mediated oxidative cyclization reaction of a silyl ketene acetal precursor that is used to install the bridgehead all-carbon quaternary stereocenter and bicyclo[3.2.1]octane present in each natural product. These studies have established a synthetic relationship among three architecturally distinct ent-kaurane diterpenoids and have forged a path for the preparation of interesting unnatural ent-kauranoid structural analogs for more thorough biological study.

A cocktail of thermally stable, chemically synthesized capture agents for the efficient detection of anti-gp41 antibodies from human sera and techniques

This thesis reports on a method to improve in vitro diagnostic assays that detect immune response, with specific application to HIV-1. The inherent polyclonal diversity of the humoral immune response was addressed by using sequential in situ click chemistry to develop a cocktail of peptide-based capture agents, the components of which were raised against different, representative anti-HIV antibodies that bind to a conserved epitope of the HIV-1 envelope protein gp41. The cocktail was used to detect anti-HIV-1 antibodies from a panel of sera collected from HIV-positive patients, with improved signal-to-noise ratio relative to the gold standard commercial recombinant protein antigen. The capture agents were stable when stored as a powder for two months at temperatures close to 60°C.

Computationally Guided Monomerization of Red Fluorescent Proteins of the Class Anthozoa

Red fluorescent proteins (RFPs) have attracted significant engineering focus because of the promise of near infrared fluorescent proteins, whose light penetrates biological tissue, and which would allow imaging inside of vertebrate animals. The RFP landscape, which numbers ~200 members, is mostly populated by engineered variants of four native RFPs, leaving the vast majority of native RFP biodiversity untouched. This is largely due to the fact that native RFPs are obligate tetramers, limiting their usefulness as fusion proteins. Monomerization has imposed critical costs on these evolved tetramers, however, as it has invariably led to loss of brightness, and often to many other adverse effects on the fluorescent properties of the derived monomeric variants. Here we have attempted to understand why monomerization has taken such a large toll on Anthozoa class RFPs, and to outline a clear strategy for their monomerization. We begin with a structural study of the far-red fluorescence of AQ143, one of the furthest red emitting RFPs. We then try to separate the problem of stable and bright fluorescence from the design of a soluble monomeric β-barrel surface by engineering a hybrid protein (DsRmCh) with an oligomeric parent that had been previously monomerized, DsRed, and a pre-stabilized monomeric core from mCherry. This allows us to use computational design to successfully design a stable, soluble, fluorescent monomer. Next we took HcRed, which is a previously unmonomerized RFP that has far-red fluorescence (λemission = 633 nm) and attempted to monomerize it making use of lessons learned from DsRmCh. We engineered two monomeric proteins by pre-stabilizing HcRed’s core, then monomerizing in stages, making use of computational design and directed evolution techniques such as error-prone mutagenesis and DNA shuffling. We call these proteins mGinger0.1 (λem = 637 nm / Φ = 0.02) and mGinger0.2 (λem = 631 nm Φ = 0.04). They are the furthest red first generation monomeric RFPs ever developed, are significantly thermostabilized, and add diversity to a small field of far-red monomeric FPs. We anticipate that the techniques we describe will be facilitate future RFP monomerization, and that further core optimization of the mGingers may allow significant improvements in brightness.

Applications of a Concise and General Strategy for the Syntheses of Transtaganolide and Basiliolide Natural Products

Herein are described the total syntheses of all members of the transtaganolide and basiliolide natural product family. Utilitzation of an Ireland–Claisen rearrangement/Diels–Alder cycloaddition cascade (ICR/DA) allowed for rapid assembly of the transtaganolide and basiliolide oxabicyclo[2.2.2]octane core. This methodology is general and was applicable to all members of the natural product family.

A brief introduction outlines all the synthetic progress previously disclosed by Lee, Dudley, and Johansson. This also includes the initial syntheses of transtaganolides C and D, as well as basiliolide B and epi-basiliolide B accomplished by Stoltz in 2011. Lastly, we discuss our racemic synthesis of basililide C and epi-basiliolide C, which utilized an ICR/DA cascade to constuct the oxabicyclo[2.2.2]octane core and formal [5+2] annulation to form the ketene-acetal containing 7-membered C-ring.

Next, we describe a strategy for an asymmetric ICR/DA cascade, by incorporation of a chiral silane directing group. This allowed for enantioselective construction of the C8 all-carbon quaternary center formed in the Ireland–Claisen rearrangement. Furthermore, a single hydride reduction and subsequent translactonization of a C4 methylester bearing oxabicyclo[2.2.2]octane core demonstrated a viable strategy for the desired skeletal rearrangement to obtain pentacyclic transtaganolides A and B. Application of the asymmetric strategy culminated in the total syntheses of (–)-transtaganolide A, (+)-transtaganolide B, (+)-transtaganolide C, and (–)-transtaganolide D. Comparison of the optical rotation data of the synthetically derived transtaganolides to that from the isolated counterparts has overarching biosynthetic implications which are discussed.

Lastly, improvement to the formal [5+2] annulation strategy is described. Negishi cross-coupling of methoxyethynyl zinc chloride using a palladium Xantphos catalyst is optimized for iodo-cyclohexene. Application of this technology to an iodo-pyrone geranyl ester allowed for formation and isolation of the eneyne product. Hydration of the enenye product forms natural metabolite basiliopyrone. Furthermore, the eneyne product can undergo an ICR/DA cascade and form transtaganolides C and D in a single step from an achiral monocyclic precursor.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Quantum Mechanics Studies of Fuel Cell Catalysts and Proton Conducting Ceramics with Validation by Experiment

We carried out quantum mechanics (QM) studies aimed at improving the performance of hydrogen fuel cells. This led to predictions of improved materials, some of which were subsequently validated with experiments by our collaborators.

In part I, the challenge was to find a replacement for the Pt cathode that would lead to improved performance for the Oxygen Reduction Reaction (ORR) while remaining stable under operational conditions and decreasing cost. Our design strategy was to find an alloy with composition Pt3M that would lead to surface segregation such that the top layer would be pure Pt, with the second and subsequent layers richer in M. Under operating conditions we expect the surface to have significant O and/or OH chemisorbed on the surface, and hence we searched for M that would remain segregated under these conditions. Using QM we examined surface segregation for 28 Pt₃M alloys, where M is a transition metal. We found that only Pt₃Os and Pt₃Ir showed significant surface segregation when O and OH are chemisorbed on the catalyst surfaces. This result indicates that Pt₃Os and Pt3Ir favor formation of a Pt-skin surface layer structure that would resist the acidic electrolyte corrosion during fuel cell operation environments. We chose to focus on Os because the phase diagram for Pt-Ir indicated that Pt-Ir could not form a homogeneous alloy at lower temperature. To determine the performance for ORR, we used QM to examine all intermediates, reaction pathways, and reaction barriers involved in the processes for which protons from the anode reactions react with O₂ to form H₂O. These QM calculations used our Poisson-Boltzmann implicit solvation model include the effects of the solvent (water with dielectric constant 78 with pH 7 at 298K). We found that the rate determination step (RDS) was the O_ad hydration reaction (O_ad + H₂O_ad -> OH_ad + OH_ad) in both cases, but that the barrier for pure Pt of 0.50 eV is reduced to 0.48 eV for Pt₃Os, which at 80 degrees C would increase the rate by 218%. We collaborated with the Pu-Wei Wu’s group to carry out experiments, where we found that the dealloying process-treated Pt2Os catalyst showed two-fold higher activity at 25 degrees C than pure Pt and that the alloy had 272% improved stability, validating our theoretical predictions.

We also carried out similar QM studies followed by experimental validation for the Os/Pt core-shell catalyst fabricated by the underpotential deposition (UPD) method. The QM results indicated that the RDS for ORR is a compromise between the OOH formation step (0.37 eV for Pt, 0.23 eV for Pt_2ML/Os core-shell) and H₂O formation steps (0.32 eV for Pt, 0.22 eV for Pt_2ML/Os core-shell). We found that Pt_2ML/Os has the highest activity (compared to pure Pt and to the Pt₃Os alloy) because the 0.37 eV barrier decreases to 0.23 eV. To understand what aspects of the core shell structure lead to this improved performance, we considered the effect on ORR of compressing the alloy slab to the dimensions of pure Pt. However this had little effect, with the same RDS barrier 0.37 eV. This shows that the ligand effect (the electronic structure modification resulting from the Os substrate) plays a more important role than the strain effect, and is responsible for the improved activity of the core- shell catalyst. Experimental materials characterization proves the core-shell feature of our catalyst. The electrochemical experiment for Pt_2ML/Os/C showed 3.5 to 5 times better ORR activity at 0.9V (vs. NHE) in 0.1M HClO₄ solution at 25 degrees C as compared to those of commercially available Pt/C. The excellent correlation between experimental half potential and the OH binding energies and RDS barriers validate the feasibility of predicting catalyst activity using QM calculation and a simple Langmuir–Hinshelwood model.

In part II, we used QM calculations to study methane stream reforming on a Ni-alloy catalyst surfaces for solid oxide fuel cell (SOFC) application. SOFC has wide fuel adaptability but the coking and sulfur poisoning will reduce its stability. Experimental results suggested that the Ni4Fe alloy improves both its activity and stability compared to pure Ni. To understand the atomistic origin of this, we carried out QM calculations on surface segregation and found that the most stable configuration for Ni₄Fe has a Fe atom distribution of (0%, 50%, 25%, 25%, 0%) starting at the bottom layer. We calculated that the binding of C atoms on the Ni4Fe surface is 142.9 Kcal/mol, which is about 10 Kcal/mol weaker compared to the pure Ni surface. This weaker C binding energy is expected to make coke formation less favorable, explaining why Ni₄Fe has better coking resistance. This result confirms the experimental observation. The reaction energy barriers for CHx decomposition and C binding on various alloy surface, Ni₄X (X=Fe, Co, Mn, and Mo), showed Ni₄Fe, Ni₄Co, and Fe₄Mn all have better coking resistance than pure Ni, but that only Ni₄Fe and Fe₄Mn have (slightly) improved activity compared to pure Ni.

In part III, we used QM to examine the proton transport in doped perovskite-ceramics. Here we used a 2x2x2 supercell of perovskite with composition Ba₈X₇M₁(OH)₁O₂₃ where X=Ce or Zr and M=Y, Gd, or Dy. Thus in each case a 4⁺ X is replace by a 3⁺ M plus a proton on one O. Here we predicted the barriers for proton diffusion allowing both includes intra-octahedron and inter-octahedra proton transfer. Without any restriction, we only observed the inter-octahedra proton transfer with similar energy barrier as previous computational work but 0.2 eV higher than experimental result for Y doped zirconate. For one restriction in our calculations is that the O_donor-O_acceptor atoms were kept at fixed distances, we found that the barrier difference between cerates/zirconates with various dopants are only 0.02~0.03 eV. To fully address performance one would need to examine proton transfer at grain boundaries, which will require larger scale ReaxFF reactive dynamics for systems with millions of atoms. The QM calculations used here will be used to train the ReaxFF force field.

Robust Control of Evolutionary Dynamics

The application of principles from evolutionary biology has long been used to gain new insights into the progression and clinical control of both infectious diseases and neoplasms. This iterative evolutionary process consists of expansion, diversification and selection within an adaptive landscape - species are subject to random genetic or epigenetic alterations that result in variations; genetic information is inherited through asexual reproduction and strong selective pressures such as therapeutic intervention can lead to the adaptation and expansion of resistant variants. These principles lie at the center of modern evolutionary synthesis and constitute the primary reasons for the development of resistance and therapeutic failure, but also provide a framework that allows for more effective control.

A model system for studying the evolution of resistance and control of therapeutic failure is the treatment of chronic HIV-1 infection by broadly neutralizing antibody (bNAb) therapy. A relatively recent discovery is that a minority of HIV-infected individuals can produce broadly neutralizing antibodies, that is, antibodies that inhibit infection by many strains of HIV. Passive transfer of human antibodies for the prevention and treatment of HIV-1 infection is increasingly being considered as an alternative to a conventional vaccine. However, recent evolution studies have uncovered that antibody treatment can exert selective pressure on virus that results in the rapid evolution of resistance. In certain cases, complete resistance to an antibody is conferred with a single amino acid substitution on the viral envelope of HIV.

The challenges in uncovering resistance mechanisms and designing effective combination strategies to control evolutionary processes and prevent therapeutic failure apply more broadly. We are motivated by two questions: Can we predict the evolution to resistance by characterizing genetic alterations that contribute to modified phenotypic fitness? Given an evolutionary landscape and a set of candidate therapies, can we computationally synthesize treatment strategies that control evolution to resistance?

To address the first question, we propose a mathematical framework to reason about evolutionary dynamics of HIV from computationally derived Gibbs energy fitness landscapes -- expanding the theoretical concept of an evolutionary landscape originally conceived by Sewall Wright to a computable, quantifiable, multidimensional, structurally defined fitness surface upon which to study complex HIV evolutionary outcomes.

To design combination treatment strategies that control evolution to resistance, we propose a methodology that solves for optimal combinations and concentrations of candidate therapies, and allows for the ability to quantifiably explore tradeoffs in treatment design, such as limiting the number of candidate therapies in the combination, dosage constraints and robustness to error. Our algorithm is based on the application of recent results in optimal control to an HIV evolutionary dynamics model and is constructed from experimentally derived antibody resistant phenotypes and their single antibody pharmacodynamics. This method represents a first step towards integrating principled engineering techniques with an experimentally based mathematical model in the rational design of combination treatment strategies and offers predictive understanding of the effects of combination therapies of evolutionary dynamics and resistance of HIV. Preliminary in vitro studies suggest that the combination antibody therapies predicted by our algorithm can neutralize heterogeneous viral populations despite containing resistant mutations.

On the Lyapunov transformation for stable matrices

The matrices studied here are positive stable (or briefly stable). These are matrices, real or complex, whose eigenvalues have positive real parts. A theorem of Lyapunov states that A is stable if and only if there exists H ˃ 0 such that AH + HA* = I. Let A be a stable matrix. Three aspects of the Lyapunov transformation L_A :H → AH + HA* are discussed.

1. Let C₁ (A) = {AH + HA* :H ≥ 0} and C₂ (A) = {H: AH+HA* ≥ 0}. The problems of determining the cones C₁(A) and C₂(A) are still unsolved. Using solvability theory for linear equations over cones it is proved that C₁(A) is the polar of C₂(A*), and it is also shown that C₁ (A) = C₁(A^-1). The inertia assumed by matrices in C₁(A) is characterized.

2. The index of dissipation of A was defined to be the maximum number of equal eigenvalues of H, where H runs through all matrices in the interior of C₂(A). Upper and lower bounds, as well as some properties of this index, are given.

3. We consider the minimal eigenvalue of the Lyapunov transform AH+HA*, where H varies over the set of all positive semi-definite matrices whose largest eigenvalue is less than or equal to one. Denote it by ψ(A). It is proved that if A is Hermitian and has eigenvalues μ₁ ≥ μ₂…≥ μ_n ˃ 0, then ψ(A) = -(μ₁-μ_n)²/(4(μ₁ + μ_n)). The value of ψ(A) is also determined in case A is a normal, stable matrix. Then ψ(A) can be expressed in terms of at most three of the eigenvalues of A. If A is an arbitrary stable matrix, then upper and lower bounds for ψ(A) are obtained.