Structurally engineered cytochromes c with novel ligand-binding properties

Semisynthesis of horse heart cytochrome c and site-directed mutagenesis of Saccharomyces cerevisiae (S. c.) iso-1-cytochrome c have been utilized to substitute Ala for the cytochrome c heme axial ligand Met80 to yield ligand-binding proteins (horse heart Ala80cyt c and S.c. Ala80cyt c) with spectroscopic properties remarkably similar to those of myoglobin. Both species of Fe(II)Ala80cyt c form exceptionally stable dioxygen complexes with autoxidation rates 10-30x smaller and O₂ binding constants ~ 3x greater than those of myoglobin. The resistance of O₂-Fe(II)Ala80cyt c to autoxidation is attributed in part to protection of the heme site from solvent as exhibited by the exceptionally slow rate of CO binding to the heme as well as the low quantum yield of CO photodissociation.

UV/vis, EPR, and paramagnetic NMR spectroscopy indicate that at pH 7 the Fe(III)Ala80cyt c heme is low-spin with axial His-OH^- coordination and that below pH ~6.5, Fe(III)Ala80cyt cis high-spin with His-H₂O heme ligation. Significant differences in the pH dependence of the ¹H NMR spectra of S.c. Fe(III)Ala80cyt c compared to wild-type demonstrate that the axial ligands influence the conformational energetics of cytochrome c.

¹H NMR spectroscopy has been utilized to determine the solution structure of the cyanide derivative of S.c. Fe(III)Ala80cyt c. 82% of the resonances in the ¹H NMR spectrum of S.c. CN-Fe(III)Ala80cyt c have been assigned through 1D and 2D experiments. The RMSD values after restrained energy minimization of the family of 17 structures obtained from distance geometry calculations are 0.68 ± 0.11 Å for the backbone and 1.32 ± 0.14 Å for all heavy atoms. The solution structure indicates that a tyrosine in the "distal" pocket of CN-Fe(III)Ala80cyt c forms a hydrogen bond with the Fe(III)-CN unit, suggesting that it may play a role analogous to that of the distal histidine in myoglobin in stabilizing the dioxygen adduct.

Cooperative properties of lipid bilayers : collective director fluctuations and the effects of hydrophobic mismatch in protein/lipid membrane systems

Cooperative director fluctuations in lipid bilayers have been postulated for many years. ^2H-NMR T_1^(-1), T_(1P)^(-1) , and T_2^(-1); measurements have been used identify these motions and to determine the origin of increased slow bilayer motion upon addition of unlike lipids or proteins to a pure lipid bilayer.

The contribution of cooperative director fluctuations to NMR relaxation in lipid bilayers has been expressed mathematically using the approach of Doane et al.^1 and Pace and Chan.^2 The T_2^(-1)’s of pure dimyristoyllecithin (DML) bilayers deuterated at the 2, 9 and 10, and all positions on both lipid hydrocarbon chains have been measured. Several characteristics of these measurements indicate the presence of cooperative director fluctuations. First of all, T_2^(-1) exhibits a linear dependence on S2/CD. Secondly, T_2^(-1) varies across the ^2H-NMR powder pattern as sin^2 (2, β), where , β is the angle between the average bilayer director and the external magnetic field. Furthermore, these fluctuations are restricted near the lecithin head group suggesting that the head group does not participate in these motions but, rather, anchors the hydrocarbon chains in the bilayer.

T_2^(-1)has been measured for selectively deuterated liquid crystalline DML hilayers to which a host of other lipids and proteins have been added. The T_2^(-1) of the DML bilayer is found to increase drastically when chlorophyll a (chl a) and Gramicidin A' (GA') are added to the bilayer. Both these molecules interfere with the lecithin head group spacing in the bilayer. Molecules such as myristic acid, distearoyllecithin (DSL), phytol, and cholesterol, whose hydrocarbon regions are quite different from DML but which have small,neutral polar head groups, leave cooperative fluctuations in the DML bilayer unchanged.

The effect of chl a on cooperative fluctuations in the DML bilayer has been examined in detail using ^2H-NMR T_1^(-1), T_(1P)^(-1) , and T_2^(-1); measurements. Cooperative fluctuations have been modelled using the continuum theory of the nematic state of liquid crystals. Chl a is found to decrease both the correlation length and the elastic constants in the DML bilayer.

A mismatch between the hydrophobic length of a lipid bilayer and that of an added protein has also been found to change the cooperative properties of the lecithin bilayer. Hydrophobic mismatch has been studied in a series GA' / lecithin bilayers. The dependence of 2H-NMR order parameters and relaxation rates on GA' concentration has been measured in selectively deuterated DML, dipalmitoyllecithin (DPL), and DSL systems. Order parameters, cooperative lengths, and elastic constants of the DML bilayer are most disrupted by GA', while the DSL bilayer is the least perturbed by GA'. Thus, it is concluded that the hydrophobic length of GA' best matches that of the DSL bilayer. Preliminary Raman spectroscopy and Differential Scanning Calorimetry experiments of GA' /lecithin systems support this conclusion. Accommodation of hydrophobic mismatch is used to rationalize the absence of H_(II) phase formation in GA' /DML systems and the observation of H_(II) phase in GA' /DPL and GA' /DSL systems.

1. J. W. Doane and D. L. Johnson, Chem. Phy3. Lett., 6, 291-295 (1970). 2. R. J. Pace and S. I. Chan, J. Chem. Phy3., 16, 4217-4227 (1982).

Synthetic, kinetic and mechanistic studies in early transition metal systems. I. α- and β-hydrogen elimination reactions in derivatives of permethyltitanocene, -zirconocene, and -hafnocene. II. The reactions of aluminum alkyls with derivatives of permethylniobocene

The thermal decomposition of Cp*Ti(CH_3)_2 (Cp*≡ ƞ^5-C_5Me_5) toluene solution follows cleanly first-order kinetics and produces a single titanium product Cp*(C_5Me_4CH_2)Ti(CH_3) concurrent with the evolution of one equivalent of methane. Labeling studies using Cp*_2Ti- (CD_3)_2 and (Cp*-d_(15))_2Ti(CH_3)_2 show the decomposition to be intramolecular and the methane to be produced by the coupling of a methyl group with a hydrogen from the other TiCH_3 group. Activation parameters, ΔH^‡ and ΔS^‡, and kinetic deuterium isotope effects have been measured. The alternative decomposition pathways of α-hydrogen abstraction and a-hydrogen elimination, both leading to a titanium-methylidene intermediate, are discussed.

The insertion of unactivated acetylenes into the metal-hydride bonds of Cp*_2MH_2 (M = Zr, Hf) proceeds rapidly at low temperature to form monoand/ or bisinsertion products, dependent upon the steric bulk of the acetylene substituents. Cp*_2M(H)(C(Me)=CHMe), Cp*_2M(H)(CH=CHCMe_3), Cp*_2M(H)-(CH=CHPh), Cp*_2M(CH=CHPh)_2, Cp*_2M(CH=CHCH_3)_2 and Cp*_2Zr- (CH=CHCH_2CH_3)_2 have been isolated and characterized. To extend the study of unsaturated-carbon ligands, Cp*_2M(C≡CCH_3)_2 have been prepared by treating Cp*_2MCl_2 with LiC≡CCH_3. The reactivity of many of these complexes with carbon monoxide and dihydrogen is surveyed. The mono(2- butenyl) complexes Cp*_2M(H)(C(Me)=CHMe) rearrange at room temperature, forming the crotyl-hydride species Cp*_2M(H)(ƞ^3-C_4H_7). The bis(propenyl) and bis(l-butenyl) zirconium complexes Cp*_2Zr(CH=CHR)_2 (R = CH_3, CH_2CH_3) also rearrange, forming zirconacyclopentenes. Labeling studies, reaction chemistry, and kinetic measurements, including deuterium isotope effects, demonstrate that the unusual 6-hydrogen elimination from an sp^2-hybridized carbon is the first step in these latter rearrangements but is not observed in the former. Details of these mechanisms and the differences in reactivity of the zirconium and hafnium complexes are discussed.

The reactions of hydride- and alkyl-carbonyl derivatives of permethylniobocene with equimolar amounts of trialkylaluminum reagents occur rapidly producing the carbonyl adducts Cp*_2Nb(R)(COAlR'_3) (R = H, CH_3, CH_2CH_3, CH_2CH_2Ph, C(Me)=CHMe; R' = Me, Et). The hydride adduct Cp*_2NbH_3•AlEt_3 has also been formed. In solution, each of these compounds exists in equilibrium with the uncomplexed species. The formation constants for Cp*_2Nb(H)(COA1R'_R) have been measured. They indicate the steric bulk of the Cp* ligands plays a deciding factor in the isolation of the first example of an aluminum Lewis acid bound to a carbonyl-oxygen in preference to a metalhydride. Reactions of Cp*_2Nb(H)CO with other Lewis acids and of the one:one adducts with H_2, CO and C_2H_4 are also discussed.

Cp*_2Nb(H)(C_2H_4) also reacts with equimolar amounts of trialkylaluminum reagents, forming a one:one complex that ^1H NMR spectroscopy indicates contains a Nb-CH_2CH_2-Al bridge. This adduct also exists in equilibrium with the uncomplexed species in solution. The formation constant for Cp*_2N+/b(H)(CH_2CH_2ĀlEt_3) has been measured. Reactions of Cp*_2Nb(H)(C_2H_4) with other Lewis acids and the reactions of Cp*_2N+b(H)- (CH_2CH_2ĀlEt_3) with CO and C_2H_4 are described, as are the reactions of Cp_*2Nb(H)(CH_2=CHR) (R = Me, Ph), Cp*_2Nb(H)(CH_3C≡CCH_3) and Cp*_2Ti-(C_2H_4) with AlEt_3.

Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

Edge function method applied to thin plates resting on Winkler foundation

The Edge Function method formerly developed by Quinlan⁽²⁵⁾ is applied to solve the problem of thin elastic plates resting on spring supported foundations subjected to lateral loads the method can be applied to plates of any convex polygonal shapes, however, since most plates are rectangular in shape, this specific class is investigated in this thesis. The method discussed can also be applied easily to other kinds of foundation models (e.g. springs connected to each other by a membrane) as long as the resulting differential equation is linear. In chapter VII, solution of a specific problem is compared with a known solution from literature. In chapter VIII, further comparisons are given. The problems of concentrated load on an edge and later on a corner of a plate as long as they are far away from other boundaries are also given in the chapter and generalized to other loading intensities and/or plates springs constants for Poisson's ratio equal to 0.2

Studies of phonon anharmonicity in solids

Today our understanding of the vibrational thermodynamics of materials at low temperatures is emerging nicely, based on the harmonic model in which phonons are independent. At high temperatures, however, this understanding must accommodate how phonons interact with other phonons or with other excitations. We shall see that the phonon-phonon interactions give rise to interesting coupling problems, and essentially modify the equilibrium and non-equilibrium properties of materials, e.g., thermodynamic stability, heat capacity, optical properties and thermal transport of materials. Despite its great importance, to date the anharmonic lattice dynamics is poorly understood and most studies on lattice dynamics still rely on the harmonic or quasiharmonic models. There have been very few studies on the pure phonon anharmonicity and phonon-phonon interactions. The work presented in this thesis is devoted to the development of experimental and computational methods on this subject.

Modern inelastic scattering techniques with neutrons or photons are ideal for sorting out the anharmonic contribution. Analysis of the experimental data can generate vibrational spectra of the materials, i.e., their phonon densities of states or phonon dispersion relations. We obtained high quality data from laser Raman spectrometer, Fourier transform infrared spectrometer and inelastic neutron spectrometer. With accurate phonon spectra data, we obtained the energy shifts and lifetime broadenings of the interacting phonons, and the vibrational entropies of different materials. The understanding of them then relies on the development of the fundamental theories and the computational methods.

We developed an efficient post-processor for analyzing the anharmonic vibrations from the molecular dynamics (MD) calculations. Currently, most first principles methods are not capable of dealing with strong anharmonicity, because the interactions of phonons are ignored at finite temperatures. Our method adopts the Fourier transformed velocity autocorrelation method to handle the big data of time-dependent atomic velocities from MD calculations, and efficiently reconstructs the phonon DOS and phonon dispersion relations. Our calculations can reproduce the phonon frequency shifts and lifetime broadenings very well at various temperatures.

To understand non-harmonic interactions in a microscopic way, we have developed a numerical fitting method to analyze the decay channels of phonon-phonon interactions. Based on the quantum perturbation theory of many-body interactions, this method is used to calculate the three-phonon and four-phonon kinematics subject to the conservation of energy and momentum, taking into account the weight of phonon couplings. We can assess the strengths of phonon-phonon interactions of different channels and anharmonic orders with the calculated two-phonon DOS. This method, with high computational efficiency, is a promising direction to advance our understandings of non-harmonic lattice dynamics and thermal transport properties.

These experimental techniques and theoretical methods have been successfully performed in the study of anharmonic behaviors of metal oxides, including rutile and cuprite stuctures, and will be discussed in detail in Chapters 4 to 6. For example, for rutile titanium dioxide (TiO₂), we found that the anomalous anharmonic behavior of the B_1g mode can be explained by the volume effects on quasiharmonic force constants, and by the explicit cubic and quartic anharmonicity. For rutile tin dioxide (SnO₂), the broadening of the B_2g mode with temperature showed an unusual concave downwards curvature. This curvature was caused by a change with temperature in the number of down-conversion decay channels, originating with the wide band gap in the phonon dispersions. For silver oxide (Ag₂O), strong anharmonic effects were found for both phonons and for the negative thermal expansion.

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.

Measurement of the neutron beta decay asymmetry using ultracold neutrons

The free neutron beta decay correlation A₀ between neutron polarization and electron emission direction provides the strongest constraint on the ratio λ = g_A/g_V of the Axial-vector to Vector coupling constants in Weak decay. In conjunction with the CKM Matrix element V_ud and the neutron lifetime τ_n, λ provides a test of Standard Model assumptions for the Weak interaction. Leading high-precision measurements of A₀ and τ_n in the 1995-2005 time period showed discrepancies with prior measurements and Standard Model predictions for the relationship between λ, τ_n, and V_ud. The UCNA experiment was developed to measure A₀ from decay of polarized ultracold neutrons (UCN), providing a complementary determination of λ with different systematic uncertainties from prior cold neutron beam experiments. This dissertation describes analysis of the dataset collected by UCNA in 2010, with emphasis on detector response calibrations and systematics. The UCNA measurement is placed in the context of the most recent τ_n results and cold neutron A₀ experiments.

Pulsed neutron experiments in graphite

Detailed pulsed neutron measurements have been performed in graphite assemblies ranging in size from 30.48 cm x 38.10 cm x 38.10 cm to 91.44 cm x 66.67 cm x 66.67 cm. Results of the measurement have been compared to a modeled theoretical computation.

In the first set of experiments, we measured the effective decay constant of the neutron population in ten graphite stacks as a function of time after the source burst. We found the decay to be non-exponential in the six smallest assemblies, while in three larger assemblies the decay was exponential over a significant portion of the total measuring interval. The decay in the largest stack was exponential over the entire ten millisecond measuring interval. The non-exponential decay mode occurred when the effective decay constant exceeded 1600 sec^( -1).

In a second set of experiments, we measured the spatial dependence of the neutron population in four graphite stacks as a function of time after the source pulse. By doing an harmonic analysis of the spatial shape of the neutron distribution, we were able to compute the effective decay constants of the first two spatial modes. In addition, we were able to compute the time dependent effective wave number of neutron distribution in the stacks.

Finally, we used a Laplace transform technique and a simple modeled scattering kernel to solve a diffusion equation for the time and energy dependence of the neutron distribution in the graphite stacks. Comparison of these theoretical results with the results of the first set of experiments indicated that more exact theoretical analysis would be required to adequately describe the experiments.

The implications of our experimental results for the theory of pulsed neutron experiments in polycrystalline media are discussed in the last chapter.

Electrocatalytic activity of transition metal substituted heteropolytungstates

The electrochemical and electrocatalytic behavior of a series of heteropolytungstate anions in which a tungsten atom in the well known Keggin structure has been replaced by an iron atom is described. All of the iron substituted ions exhibit a one electron reversible couple associated with the Fe³⁺ center and a pair of two electron waves attributed to electron addition and removal from the tungsten oxo framework. The pH and ionic strength effects upon the various electrochemical processes are discussed and interpreted in terms of a competition between protonation and ion pairing of the highly negatively charged ions.

The anions are efficient catalysts for the electroreduction of hydrogen peroxide. A catalytic mechanism involving a formally Fe(IV) intermediate is proposed. Pulse radiolysis experiments were employed to detect the intermediate and evaluate the rate constants for the reactions in which it is formed and decomposed. A chain mechanism for the catalytic decomposition of hydrogen peroxide in which the Fe center shuttles between the +2, +3, and +4 oxidation states is proposed to explain the non-integral stoichiometry observed for the iron substituted polytungstate catalyzed electroreduction of hydrogen peroxide.

The anions are also efficient electrocatalyst for the electrochemical conversion of nitric oxide to ammonia. The catalyzed reduction does not produce hydroxylamine as an intermediate and appears to depend upon the ability of the multiply reduced heteropolytungstates to deliver several electrons to the bound NO group in a concerted step. A valuable feature of the heteropolytungstates is the ease at which the formal potentials of the several redox couples they exhibit may be shifted by changing the identity of the central heteroatom. Exploitation of this feature provided diagnostic information that was decisive in establishing the mechanism of electrocatalytic reduction.

The iron substituted heteropolytungstates are not degraded by repeated cycling between their oxidized and reduced states. They also show superior activity compared to their unsubstituted analogues, indicating that the Fe center acts as a binding site that facilitates inner-sphere electron transfer processes. The basic electrochemistry of several other transition metal substituted Keggin ions is also described.

Wastewater Electrolysis Cell for Environmental Pollutants Degradation and Molecular Hydrogen Generation

This study proposes a wastewater electrolysis cell (WEC) for on-site treatment of human waste coupled with decentralized molecular H₂ production. The core of the WEC includes mixed metal oxides anodes functionalized with bismuth doped TiO₂ (BiO_x/TiO₂). The BiO_x/TiO₂ anode shows reliable electro-catalytic activity to oxidize Cl- to reactive chlorine species (RCS), which degrades environmental pollutants including chemical oxygen demand (COD), protein, NH4⁺, urea, and total coliforms. The WEC experiments for treatment of various kinds of synthetic and real wastewater demonstrate sufficient water quality of effluent for reuse for toilet flushing and environmental purposes. Cathodic reduction of water and proton on stainless steel cathodes produced molecular H2 with moderate levels of current and energy efficiency. This thesis presents a comprehensive environmental analysis together with kinetic models to provide an in-depth understanding of reaction pathways mediated by the RCS and the effects of key operating parameters. The latter part of this thesis is dedicated to bilayer hetero-junction anodes which show enhanced generation efficiency of RCS and long-term stability.

Chapter 2 describes the reaction pathway and kinetics of urea degradation mediated by electrochemically generated RCS. The urea oxidation involves chloramines and chlorinated urea as reaction intermediates, for which the mass/charge balance analysis reveals that N₂ and CO₂ are the primary products. Chapter 3 investigates direct-current and photovoltaic powered WEC for domestic wastewater treatment, while Chapter 4 demonstrates the feasibility of the WEC to treat model septic tank effluents. The results in Chapter 2 and 3 corroborate the active roles of chlorine radicals (Cl•/Cl₂^-•) based on iR-compensated anodic potential (thermodynamic basis) and enhanced pseudo-first-order rate constants (kinetic basis). The effects of operating parameters (anodic potential and [Cl^-] in Chapter 3; influent dilution and anaerobic pretreatment in Chapter 4) on the rate and current/energy efficiency of pollutants degradation and H₂ production are thoroughly discussed based on robust kinetic models. Chapter 5 reports the generation of RCS on Ir_0.7Ta_0.3O_y/Bi_xTi_1-xO_z hetero-junction anodes with enhanced rate, current efficiency, and long-term stability compared to the Ir_0.7Ta_0.3O_y anode. The effects of surficial Bi concentration are interrogated, focusing on relative distributions between surface-bound hydroxyl radical and higher oxide.

Convex Relaxation for Low-Dimensional Representation: Phase Transitions and Limitations

There is a growing interest in taking advantage of possible patterns and structures in data so as to extract the desired information and overcome the curse of dimensionality. In a wide range of applications, including computer vision, machine learning, medical imaging, and social networks, the signal that gives rise to the observations can be modeled to be approximately sparse and exploiting this fact can be very beneficial. This has led to an immense interest in the problem of efficiently reconstructing a sparse signal from limited linear observations. More recently, low-rank approximation techniques have become prominent tools to approach problems arising in machine learning, system identification and quantum tomography.

In sparse and low-rank estimation problems, the challenge is the inherent intractability of the objective function, and one needs efficient methods to capture the low-dimensionality of these models. Convex optimization is often a promising tool to attack such problems. An intractable problem with a combinatorial objective can often be "relaxed" to obtain a tractable but almost as powerful convex optimization problem. This dissertation studies convex optimization techniques that can take advantage of low-dimensional representations of the underlying high-dimensional data. We provide provable guarantees that ensure that the proposed algorithms will succeed under reasonable conditions, and answer questions of the following flavor:

More specifically, i) Focusing on linear inverse problems, we generalize the classical error bounds known for the least-squares technique to the lasso formulation, which incorporates the signal model. ii) We show that intuitive convex approaches do not perform as well as expected when it comes to signals that have multiple low-dimensional structures simultaneously. iii) Finally, we propose convex relaxations for the graph clustering problem and give sharp performance guarantees for a family of graphs arising from the so-called stochastic block model. We pay particular attention to the following aspects. For i) and ii), we aim to provide a general geometric framework, in which the results on sparse and low-rank estimation can be obtained as special cases. For i) and iii), we investigate the precise performance characterization, which yields the right constants in our bounds and the true dependence between the problem parameters.

On the Role of Delays in Biological Systems : Analysis and Design

This work quantifies the nature of delays in genetic regulatory networks and their effect on system dynamics. It is known that a time lag can emerge from a sequence of biochemical reactions. Applying this modeling framework to the protein production processes, delay distributions are derived in a stochastic (probability density function) and deterministic setting (impulse function), whilst being shown to be equivalent under different assumptions. The dependence of the distribution properties on rate constants, gene length, and time-varying temperatures is investigated. Overall, the distribution of the delay in the context of protein production processes is shown to be highly dependent on the size of the genes and mRNA strands as well as the reaction rates. Results suggest longer genes have delay distributions with a smaller relative variance, and hence, less uncertainty in the completion times, however, they lead to larger delays. On the other hand large uncertainties may actually play a positive role, as broader distributions can lead to larger stability regions when this formalization of the protein production delays is incorporated into a feedback system.

Furthermore, evidence suggests that delays may play a role as an explicit design into existing controlling mechanisms. Accordingly, the reccurring dual-feedback motif is also investigated with delays incorporated into the feedback channels. The dual-delayed feedback is shown to have stabilizing effects through a control theoretic approach. Lastly, a distributed delay based controller design method is proposed as a potential design tool. In a preliminary study, the dual-delayed feedback system re-emerges as an effective controller design.

Effects of self energy of the ions on the double layer structure and properties at the dielectric interface

Although numerous theoretical efforts have been put forth, a systematic, unified and predictive theoretical framework that is able to capture all the essential physics of the interfacial behaviors of ions, such as the Hofmeister series effect, Jones-Ray effect and the salt effect on the bubble coalescence remain an outstanding challenge. The most common approach to treating electrostatic interactions in the presence of salt ions is the Poisson-Boltzmann (PB) theory. However, there are many systems for which the PB theory fails to offer even a qualitative explanation of the behavior, especially for ions distributed in the vicinity of an interface with dielectric contrast between the two media (like the water-vapor/oil interface). A key factor missing in the PB theory is the self energy of the ion.

In this thesis, we develop a self-consistent theory that treats the electrostatic self energy (including both the short-range Born solvation energy and the long-range image charge interactions), the nonelectrostatic contribution of the self energy, the ion-ion correlation and the screening effect systematically in a single framework. By assuming a finite charge spread of the ion instead of using the point-charge model, the self energy obtained by our theory is free of the divergence problems and gives a continuous self energy across the interface. This continuous feature allows ions on the water side and the vapor/oil side of the interface to be treated in a unified framework. The theory involves a minimum set of parameters of the ion, such as the valency, radius, polarizability of the ions, and the dielectric constants of the medium, that are both intrinsic and readily available. The general theory is first applied to study the thermodynamic property of the bulk electrolyte solution, which shows good agreement with the experiment result for predicting the activity coefficient and osmotic coefficient.

Next, we address the effect of local Born solvation energy on the bulk thermodynamics and interfacial properties of electrolyte solution mixtures. We show that difference in the solvation energy between the cations and anions naturally gives rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The miscibility of the mixture can either increases or decreases depending on the competition between the solvation energy and translation entropy of the ions. The interfacial tension shows a non-monotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations, and decreases approximately as the square root of the salt concentration for dilute solutions, which is in agreement with the Jones-Ray effect observed in experiment.

Next, we investigate the image effects on the double layer structure and interfacial properties near a single charged plate. We show that the image charge repulsion creates a depletion boundary layer that cannot be captured by a regular perturbation approach. The correct weak-coupling theory must include the self-energy of the ion due to the image charge interaction. The image force qualitatively alters the double layer structure and properties, and gives rise to many non-PB effects, such as nonmonotonic dependence of the surface energy on concentration and charge inversion. The image charge effect is then studied for electrolyte solutions between two plates. For two neutral plates, we show that depletion of the salt ions by the image charge repulsion results in short-range attractive and long-range repulsive forces. If cations and anions are of different valency, the asymmetric depletion leads to the formation of an induced electrical double layer. For two charged plates, the competition between the surface charge and the image charge effect can give rise to like- charge attraction.

Then, we study the inhomogeneous screening effect near the dielectric interface due to the anisotropic and nonuniform ion distribution. We show that the double layer structure and interfacial properties is drastically affected by the inhomogeneous screening if the bulk Debye screening length is comparable or smaller than the Bjerrum length. The width of the depletion layer is characterized by the Bjerrum length, independent of the salt concentration. We predict that the negative adsorption of ions at the interface increases linearly with the salt concentration, which cannot be captured by either the bulk screening approximation or the WKB approximation. For asymmetric salt, the inhomogeneous screening enhances the charge separation in the induced double layer and significantly increases the value of the surface potential.

Finally, to account for the ion specificity, we study the self energy of a single ion across the dielectric interface. The ion is considered to be polarizable: its charge distribution can be self-adjusted to the local dielectric environment to minimize the self energy. Using intrinsic parameters of the ions, such as the valency, radius, and polarizability, we predict the specific ion effect on the interfacial affinity of halogen anions at the water/air interface, and the strong adsorption of hydrophobic ions at the water/oil interface, in agreement with experiments and atomistic simulations.

The theory developed in this work represents the most systematic theoretical technique for weak-coupling electrolytes. We expect the theory to be more useful for studying a wide range of structural and dynamic properties in physicochemical, colloidal, soft-matter and biophysical systems.

Investigation of competitive antagonist binding to the nicotinic acetylcholine receptor using voltage-jump and light-flash techniques

1. The effect of 2,2’-bis-[α-(trimethylammonium)methyl]azobenzene (2BQ), a photoisomerizable competitive antagonist, was studied at the nicotinic acetycholine receptor of Electrophorus electroplaques using voltage-jump and light-flash techniques.

2. 2BQ, at concentrations below 3 μΜ, reduced the amplitude of voltage-jump relaxations but had little effect on the voltage-jump relaxation time constants under all experimental conditions. At higher concentrations and voltages more negative than -150 mV, 2BQ caused significant open channel blockade.

3. Dose-ratio studies showed that the cis and trans isomers of 2BQ have equilibrium binding constants (K _ᵢ) of .33 and 1.0 μΜ, respectively. The binding constants determined for both isomers are independent of temperature, voltage, agonist concentration, and the nature of the agonist.

4. In a solution of predominantly cis-2BQ, visible-light flashes led to a net cis→trans isomerization and caused an increase in the agonist-induced current. This increase had at least two exponential components; the larger amplitude component had the same time constant as a subsequent voltage-jump relaxation; the smaller amplitude component was investigated using ultraviolet light flashes.

5. In a solution of predominantly trans-2BQ, UV-light flashes led to a net trans→cis isomerization and caused a net decrease in the agonist-induced current. This effect had at least two exponential components. The smaller and faster component was an increase in agonist-induced current and had a similar time constant to the voltage-jump relaxation. The larger component was a slow decrease in the agonist-induced current with rate constant approximately an order of magnitude less than that of the voltage-jump relaxation. This slow component provided a measure of the rate constant for dissociation of cis-2BQ (k_ = 60/s at 20°C). Simple modelling of the slope of the dose-rate curves yields an association rate constant of 1.6 x 10⁸/M/s. This agrees with the association rate constant of 1.8 x 10⁸/M/s estimated from the binding constant (K_i). The Q₁₀ of the dissociation rate constant of cis-2BQ was 3.3 between 6° and 20°C. The rate constants for association and dissociation of cis-28Q at receptors are independent of voltage, agonist concentration, and the nature of the agonist.

6. We have measured the molecular rate constants of a competitive antagonist which has roughly the same K _ᵢ as d-tubocurarine but interacts more slowly with the receptor. This leads to the conclusion that curare itself has an association rate constant of 4 x 10⁹/M/s or roughly as fast as possible for an encounter-limited reaction.