Olefin cyclopropanation and carbon-hydrogen amination via carbene and nitrene transfers catalyzed by engineered cytochrome P450 enzymes

Synthetic biology promises to transform organic synthesis by enabling artificial catalysis in living cells. I start by reviewing the state of the art in this young field and recognizing that new approaches are required for designing enzymes that catalyze nonnatural reactions, in order to expand the scope of biocatalytic transformations. Carbene and nitrene transfers to C=C and C-H bonds are reactions of tremendous synthetic utility that lack biological counterparts. I show that various heme proteins, including cytochrome P450BM3, will catalyze promiscuous levels of olefin cyclopropanation when provided with the appropriate synthetic reagents (e.g., diazoesters and styrene). Only a few amino acid substitutions are required to install synthetically useful levels of stereoselective cyclopropanation activity in P450BM3. Understanding that the ferrous-heme is the active species for catalysis and that the artificial reagents are unable to induce a spin-shift-dependent increase in the redox potential of the ferric P450, I design a high-potential serine-heme ligated P450 (P411) that can efficiently catalyze cyclopropanation using NAD(P)H. Intact E. coli whole-cells expressing P411 are highly efficient asymmetric catalysts for olefin cyclopropanation. I also show that engineered P450s can catalyze intramolecular amination of benzylic C-H bonds from arylsulfonyl azides. Finally, I review other examples of where synthetic reagents have been used to drive the evolution of novel enzymatic activity in the environment and in the laboratory. I invoke preadaptation to explain these observations and propose that other man-invented reactions may also be transferrable to natural enzymes by using a mechanism-based approach for choosing the enzymes and the reagents. Overall, this work shows that existing enzymes can be readily adapted for catalysis of synthetically important reactions not previously observed in nature.

Proton transfers at the air-water interface

Proton transfer reactions at the interface of water with hydrophobic media, such as air or lipids, are ubiquitous on our planet. These reactions orchestrate a host of vital phenomena in the environment including, for example, acidification of clouds, enzymatic catalysis, chemistries of aerosol and atmospheric gases, and bioenergetic transduction. Despite their importance, however, quantitative details underlying these interactions have remained unclear. Deeper insight into these interfacial reactions is also required in addressing challenges in green chemistry, improved water quality, self-assembly of materials, the next generation of micro-nanofluidics, adhesives, coatings, catalysts, and electrodes. This thesis describes experimental and theoretical investigation of proton transfer reactions at the air-water interface as a function of hydration gradients, electrochemical potential, and electrostatics. Since emerging insights hold at the lipid-water interface as well, this work is also expected to aid understanding of complex biological phenomena associated with proton migration across membranes.

Based on our current understanding, it is known that the physicochemical properties of the gas-phase water are drastically different from those of bulk water. For example, the gas-phase hydronium ion, H₃O⁺(g), can protonate most (non-alkane) organic species, whereas H₃O⁺(aq) can neutralize only relatively strong bases. Thus, to be able to understand and engineer water-hydrophobe interfaces, it is imperative to investigate this fluctuating region of molecular thickness wherein the ‘function’ of chemical species transitions from one phase to another via steep gradients in hydration, dielectric constant, and density. Aqueous interfaces are difficult to approach by current experimental techniques because designing experiments to specifically sample interfacial layers (< 1 nm thick) is an arduous task. While recent advances in surface-specific spectroscopies have provided valuable information regarding the structure of aqueous interfaces, but structure alone is inadequate to decipher the function. By similar analogy, theoretical predictions based on classical molecular dynamics have remained limited in their scope.

Recently, we have adapted an analytical electrospray ionization mass spectrometer (ESIMS) for probing reactions at the gas-liquid interface in real time. This technique is direct, surface-specific,and provides unambiguous mass-to-charge ratios of interfacial species. With this innovation, we have been able to investigate the following:

1. How do anions mediate proton transfers at the air-water interface?

2. What is the basis for the negative surface potential at the air-water interface?

3. What is the mechanism for catalysis ‘on-water’?

In addition to our experiments with the ESIMS, we applied quantum mechanics and molecular dynamics to simulate our experiments toward gaining insight at the molecular scale. Our results unambiguously demonstrated the role of electrostatic-reorganization of interfacial water during proton transfer events. With our experimental and theoretical results on the ‘superacidity’ of the surface of mildly acidic water, we also explored implications on atmospheric chemistry and green chemistry. Our most recent results explained the basis for the negative charge of the air-water interface and showed that the water-hydrophobe interface could serve as a site for enhanced autodissociation of water compared to the condensed phase.

Macaque lateral intraparietal area and oculomotor behaviors

Neurons in the primate lateral intraparietal area (area LIP) carry visual, saccade-related and eye position activities. The visual and saccade activities are anchored in a retinotopic framework and the overall response magnitude is modulated by eye position. It was proposed that the modulation by eye position might be the basis of a distributed coding of target locations in a head-centered space. Other recording studies demonstrated that area LIP is involved in oculomotor planning. These results overall suggest that area LIP transforms sensory information for motor functions. In this thesis I further explore the role of area LIP in processing saccadic eye movements by observing the effects of reversible inactivation of this area. Macaque monkeys were trained to do visually guided and memory saccades and a double saccade task to examine the use of eye position signal. Finally, by intermixing visual saccades with trials in which two targets were presented at opposite sides of the fixation point, I examined the behavior of visual extinction.

In chapter 2, I will show that lesion of area LIP results in increased latency of contralesional visual and memory saccades. Contralesional memory saccades are also hypometric and slower in velocity. Moreover, the impairment of memory saccades does not vary with the duration of the delay period. This suggests that the oculomotor deficits observed after inactivation of area LIP is not due to the disruption of spatial memory.

In chapter 3, I will show that lesion of area LIP does not severely affect the processing of spontaneous eye movement. However, the monkeys made fewer contralesional saccades and tended to confine their gaze to the ipsilesional field after inactivation of area LIP. On the other hand, lesion of area LIP results in extinction of the contralesional stimulus. When the initial fixation position was varied so that the retinal and spatial locations of the targets could be dissociated, it was found that the extinction behavior could best be described in a head-centered coordinate.

In chapter 4, I will show that inactivation of area LIP disrupts the use of eye position signal to compute the second movement correctly in the double saccade task. If the first saccade steps into the contralesional field, the error rate and latency of the second saccade are both increased. Furthermore, the direction of the first eye movement largely does not have any effect on the impairment of the second saccade. I will argue that this study provides important evidence that the extraretinal signal used for saccadic localization is eye position rather than a displacement vector.

In chapter 5, I will demonstrate that in parietal monkeys the eye drifts toward the lesion side at the end of the memory saccade in darkness. This result suggests that the eye position activity in the posterior parietal cortex is active in nature and subserves gaze holding.

Overall, these results further support the view that area LIP neurons encode spatial locations in a craniotopic framework and is involved in processing voluntary eye movements.

The strategic behavior of rational novices

There is a growing amount of experimental evidence that suggests people often deviate from the predictions of game theory. Some scholars attempt to explain the observations by introducing errors into behavioral models. However, most of these modifications are situation dependent and do not generalize. A new theory, called the rational novice model, is introduced as an attempt to provide a general theory that takes account of erroneous behavior. The rational novice model is based on two central principals. The first is that people systematically make inaccurate guesses when they are evaluating their options in a game-like situation. The second is that people treat their decisions similar to a portfolio problem. As a result, non optimal actions in a game theoretic sense may be included in the rational novice strategy profile with positive weights.

The rational novice model can be divided into two parts: the behavioral model and the equilibrium concept. In a theoretical chapter, the mathematics of the behavioral model and the equilibrium concept are introduced. The existence of the equilibrium is established. In addition, the Nash equilibrium is shown to be a special case of the rational novice equilibrium. In another chapter, the rational novice model is applied to a voluntary contribution game. Numerical methods were used to obtain the solution. The model is estimated with data obtained from the Palfrey and Prisbrey experimental study of the voluntary contribution game. It is found that the rational novice model explains the data better than the Nash model. Although a formal statistical test was not used, pseudo R^2 analysis indicates that the rational novice model is better than a Probit model similar to the one used in the Palfrey and Prisbrey study.

The rational novice model is also applied to a first price sealed bid auction. Again, computing techniques were used to obtain a numerical solution. The data obtained from the Chen and Plott study were used to estimate the model. The rational novice model outperforms the CRRAM, the primary Nash model studied in the Chen and Plott study. However, the rational novice model is not the best amongst all models. A sophisticated rule-of-thumb, called the SOPAM, offers the best explanation of the data.

Search for color transparency in A(e,e'p) at high momentum transfer

Rates for A(e, e'p) on the nuclei ^2H, C, Fe, and Au have been measured at momentum transfers Q^2 = 1, 3, 5, and 6.8 (GeV fc)^2 . We extract the nuclear transparency T, a measure of the importance of final state interactions (FSI) between the outgoing proton and the recoil nucleus. Some calculations based on perturbative QCD predict an increase in T with momentum transfer, a phenomenon known as Color Transparency. No statistically significant rise is seen in the present experiment.

Limits and tradeoffs in the control of autocatalytic systems

Despite the complexity of biological networks, we find that certain common architectures govern network structures. These architectures impose fundamental constraints on system performance and create tradeoffs that the system must balance in the face of uncertainty in the environment. This means that while a system may be optimized for a specific function through evolution, the optimal achievable state must follow these constraints. One such constraining architecture is autocatalysis, as seen in many biological networks including glycolysis and ribosomal protein synthesis. Using a minimal model, we show that ATP autocatalysis in glycolysis imposes stability and performance constraints and that the experimentally well-studied glycolytic oscillations are in fact a consequence of a tradeoff between error minimization and stability. We also show that additional complexity in the network results in increased robustness. Ribosome synthesis is also autocatalytic where ribosomes must be used to make more ribosomal proteins. When ribosomes have higher protein content, the autocatalysis is increased. We show that this autocatalysis destabilizes the system, slows down response, and also constrains the system’s performance. On a larger scale, transcriptional regulation of whole organisms also follows architectural constraints and this can be seen in the differences between bacterial and yeast transcription networks. We show that the degree distributions of bacterial transcription network follow a power law distribution while the yeast network follows an exponential distribution. We then explored the evolutionary models that have previously been proposed and show that neither the preferential linking model nor the duplication-divergence model of network evolution generates the power-law, hierarchical structure found in bacteria. However, in real biological systems, the generation of new nodes occurs through both duplication and horizontal gene transfers, and we show that a biologically reasonable combination of the two mechanisms generates the desired network.

A measurement of inclusive quasi-elastic electron scattering cross sections at high momentum transfer

We have measured inclusive electron-scattering cross sections for targets of ^(4)He, C, Al, Fe, and Au, for kinematics spanning the quasi-elastic peak, with squared, four­ momentum transfers (q^2) between 0.23 and 2.89 (GeV/c)^2. Additional data were measured for Fe with q^2's up to 3.69 (GeV/c)^2 These cross sections were analyzed for the y-scaling behavior expected from a simple, impulse-approximation model, and are found to approach a scaling limit at the highest q^2's. The q^2 approach to scaling is compared with a calculation for infinite nuclear matter, and relationships between the scaling function and nucleon momentum distributions are discussed. Deviations from perfect scaling are used to set limits on possible changes in the size of nucleons inside the nucleus.

Frictional properties of faults : from observation on the Longitudinal Valley Fault, Taiwan, to dynamic simulations

Faults can slip either aseismically or through episodic seismic ruptures, but we still do not understand the factors which determine the partitioning between these two modes of slip. This challenge can now be addressed thanks to the dense set of geodetic and seismological networks that have been deployed in various areas with active tectonics. The data from such networks, as well as modern remote sensing techniques, indeed allow documenting of the spatial and temporal variability of slip mode and give some insight. This is the approach taken in this study, which is focused on the Longitudinal Valley Fault (LVF) in Eastern Taiwan. This fault is particularly appropriate since the very fast slip rate (about 5 cm/yr) is accommodated by both seismic and aseismic slip. Deformation of anthropogenic features shows that aseismic creep accounts for a significant fraction of fault slip near the surface, but this fault also released energy seismically, since it has produced five M_w>6.8 earthquakes in 1951 and 2003. Moreover, owing to the thrust component of slip, the fault zone is exhumed which allows investigation of deformation mechanisms. In order to put constraint on the factors that control the mode of slip, we apply a multidisciplinary approach that combines modeling of geodetic observations, structural analysis and numerical simulation of the "seismic cycle". Analyzing a dense set of geodetic and seismological data across the Longitudinal Valley, including campaign-mode GPS, continuous GPS (cGPS), leveling, accelerometric, and InSAR data, we document the partitioning between seismic and aseismic slip on the fault. For the time period 1992 to 2011, we found that about 80-90% of slip on the LVF in the 0-26 km seismogenic depth range is actually aseismic. The clay-rich Lichi M\'elange is identified as the key factor promoting creep at shallow depth. Microstructural investigations show that deformation within the fault zone must have resulted from a combination of frictional sliding at grain boundaries, cataclasis and pressure solution creep. Numerical modeling of earthquake sequences have been performed to investigate the possibility of reproducing the results from the kinematic inversion of geodetic and seismological data on the LVF. We first investigate the different modeling strategy that was developed to explore the role and relative importance of different factors on the manner in which slip accumulates on faults. We compare the results of quasi dynamic simulations and fully dynamic ones, and we conclude that ignoring the transient wave-mediated stress transfers would be inappropriate. We therefore carry on fully dynamic simulations and succeed in qualitatively reproducing the wide range of observations for the southern segment of the LVF. We conclude that the spatio-temporal evolution of fault slip on the Longitudinal Valley Fault over 1997-2011 is consistent to first order with prediction from a simple model in which a velocity-weakening patch is embedded in a velocity-strengthening area.

High energy photoproduction of neutral mesons

Cross sections for the photoproduction of neutral pi, eta, rho and phi mesons on hydrogen have been measured at the Stanford Linear Accelerator Center using a missing mass spectrometer technique. The data cover photon energies between 5.0 and 17.8 GeV and four momentum transfer squared t between -.12 and -1.38 (GeV/c)².

Pion differential cross sections at lower energies show a peak at low momentum transfers, a distinctive dip and secondary maximum for t in the region -.4 to -.9 (GeV /c)², and a smooth decrease at higher momentum transfers. As photon energy increases, the dip becomes less pronounced, in contradiction to the expectations of simple Regge theories based on the exchange of omega and B trajectories only.

Eta photoproduction was measured only below 10 GeV. The cross section has about the same magnitude as the pion production cross section, but decreases exponentially with t, showing no dip.

Rho mesons appear to be diffractively produced. The differential cross section varies approximately as exp(8.5t + 2t²). It falls slowly with energy, decreasing about 35 percent from 6 GeV to 17.8 GeV. A simple quark model relation appears to describe the data well.

Phi meson cross sections are also consistent with diffraction production. The differential cross section varies approximately as exp(4t). The cross section tends to decrease slightly with photon energy.

Chelation-enforced metal-arene interactions : insights into substrate binding and catalysis by late transition metal complexes

Understanding and catalyzing chemical reactions requiring multiple electron transfers is an endeavor relevant to many outstanding challenges in the field of chemistry. To study multi-electron reactions, a terphenyl diphosphine framework was designed to support one or more metals in multiple redox states via stabilizing interactions with the central arene of the terphenyl backbone. A variety of unusual compounds and reactions and their relevance toward prominent research efforts in chemistry are the subject of this dissertation.

Chapter 2 introduces the para-terphenyl diphosphine framework and its coordination chemistry with group 10 transition metal centers. Both mononuclear and dinuclear compounds are characterized. In many cases, the metal center(s) are stabilized by the terphenyl central arene. These metal–arene interactions are characterized both statically, in the solid state, and fluxionally, in solution. As a proof-of-principle, a dinickel framework is shown to span multiple redox states, showing that multielectron chemistry can be supported by the coordinatively flexible terphenyl diphosphine.

Chapter 3 presents reactivity of the terphenyl diphosphine when bound to a metal center. Because of the dearomatizing effect of the metal center, the central arene of the ligand is susceptible to reactions that do not normally affect arenes. In particular, Ni-to-arene H-transfer and arene dihydrogenation reactions are presented. Additionally, evidence for reversibility of the Ni-to-arene H-transfer is discussed.

Chapter 4 expands beyond the chelated metal-arene interactions of the previous chapters. A dipalladium(I) terphenyl diphosphine framework is used to bind a variety of exogenous organic ligands including arenes, dienes, heteroarenes, thioethers, and anionic ligands. The compounds are structurally characterized, and many ligands exhibit unprecedented bindng modes across two metal centers. The relative binding affinities are evaluated spectroscopically, and equilibrium binding constants for the examined ligands are determined to span over 13 orders of magnitude. As an application of this framework, mild hydrogenation conditions of bound thiophene are presented.

Chapter 5 studies nickel-mediated C–O bond cleavage of aryl alkyl ethers, a transformation with emerging applications in fields such as lignin biofuels and organic methodology. Other group members have shown the mechanism of C–O bond cleavage of an aryl methyl ether incorporated into a meta-terphenyl diphosphine framework to proceed through β-H elimination of an alkoxide. First, the electronic selectivity of the model system is examined computationally and compared with catalytic systems. The lessons learned from the model system are then applied to isotopic labeling studies for catalytic aryl alkyl ether cleavage under dihydrogen. Results from selective deuteration experiments and mass spectrometry draw a clear analogy between the mechanisms of the model and catalytic systems that does not require dihydrogen for C–O bond cleavage, although dihydrogen is proposed to play a role in catalyst activation and catalytic turnover.

Appendix A presents initial efforts toward heterodinuclear complexes as models for CO dehydrogenase and Fischer Tropsch chemistry. A catechol-incorporating terphenyl diphosphine is reported, and metal complexes thereof are discussed.

Appendix B highlights some structurally characterized terphenyl diphosphine complexes that either do not thematically belong in the research chapters or proved to be difficult to reproduce. These compounds show unusual coordination modes of the terphenyl diphosphine from which other researchers may glean insights.