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.

Excitation energies and decay properties of T = 3/2 states in ^(17)O, ^(17)F and ^(21)Na

The two lowest T = 3/2 levels in ²¹Na have been studied in the ¹⁹F(³He, n), ²⁰Ne (p,p) and ²⁰Ne (p,p’) reactions, and their excitation energies, spins, parities and widths have been determined. In a separate investigation, branching ratios were measured for the isospin-nonconserving particle decays of the lowest T = 3/2 levels in ¹⁷O and ¹⁷F to the ground state and first two excited states of ¹⁶O, by studying the ¹⁵N(³He,n) ¹⁷F*(p) ¹⁶O and ¹⁸O(³He, α)¹⁷O*(n) ¹⁶O reactions.

The ¹⁹F(³He,n) ²¹Na reaction was studied at incident energies between 4.2 and 5.9 MeV using a pulsed-beam neutron-time-of-flight spectrometer. Two T = 3/2 levels were identified at excitation energies of 8.99 ± 0.05 MeV (J > ½) and 9.22 ± 0.015 MeV (J ^π = ½⁺, Γ ˂ 40 keV). The spins and parities were determined by a comparison of the measured angular distributions with the results of DWBA calculations.

These two levels were also obsesrved as isospin-forbidden resonances in the ²⁰Ne(p,p) and ²⁰Ne(p,p’) reactions. Excitation energies were measured and spins, parities, and widths were determined from a single level dispersion theory analysis. The following results were obtained:

E_x = 8.973 ± 0.007 MeV, J ^π = 5/2 ⁺ or 3/2⁺, Γ ≤ 1.2 keV,

^Γp_o = 0.1 ± 0.05 keV; E_x = 9.217 ± 0.007 MeV, J^π = ½ ⁺,

Γ = 2.3 ± 0.5 keV, ^Γp_o = 1.1 ± 0.3 keV.

Isospin assignments were made on the basis of excitation energies, spins, parities, and widths.

Branching ratios for the isospin-nonconserving proton decays of the 11.20 MeV, T = 3/2 level in ¹⁷F were measured by the ¹⁵N(³He,n) ¹⁷ F*(p) ¹⁶O reaction to be 0.088 ± 0.016 to the ground state of ¹⁶O and 0.22 ± 0.04 to the unresolved 6.05 and 6.13 MeV levels of ¹⁶O. Branching ratios for the neutron decays of the analogous T = 3/2 level, at 11.08 MeV in ¹⁷O, were measured by the ¹⁶O(³He, α)¹⁷O*(n)¹⁶O reaction to be 0.91 ± 0.15 to the ground state of ¹⁶O and 0.05 ± 0.02 to the unresolved 6.05 and 6.13 MeV states. By comparing the ratios of reduced widths for the mirror decays, the form of the isospin impurity in the T = 3/2 levels is shown to depend on T_z.

Nuclear Hartree-Fock calculations on parallel computers

Hartree-Fock (HF) calculations have had remarkable success in describing large nuclei at high spin, temperature and deformation. To allow full range of possible deformations, the Skyrme HF equations can be discretized on a three-dimensional mesh. However, such calculations are currently limited by the computational resources provided by traditional supercomputers. To take advantage of recent developments in massively parallel computing technology, we have implemented the LLNL Skyrme-force static and rotational HF codes on Intel's DELTA and GAMMA systems at Caltech.

We decomposed the HF code by assigning a portion of the mesh to each node, with nearest neighbor meshes assigned to nodes connected by communication· channels. This kind of decomposition is well-suited for the DELTA and the GAMMA architecture because the only non-local operations are wave function orthogonalization and the boundary conditions of the Poisson equation for the Coulomb field.

Our first application of the HF code on parallel computers has been the study of identical superdeformed (SD) rotational bands in the Hg region. In the last ten years, many SD rotational bands have been found experimentally. One very surprising feature found in these SD rotational bands is that many pairs of bands in nuclei that differ by one or two mass units have nearly identical deexcitation gamma-ray energies. Our calculations of the five rotational bands in ^(192)Hg and ^(194)Pb show that the filling of specific orbitals can lead to bands with deexcitation gamma-ray energies differing by at most 2 keV in nuclei differing by two mass units and over a range of angular momenta comparable to that observed experimentally. Our calculations of SD rotational bands in the Dy region also show that twinning can be achieved by filling or emptying some specific orbitals.

The interpretation of future precise experiments on atomic parity nonconservation (PNC) in terms of parameters of the Standard Model could be hampered by uncertainties in the atomic and nuclear structure. As a further application of the massively parallel HF calculations, we calculated the proton and neutron densities of the Cesium isotopes from A = 125 to A = 139. Based on our good agreement with experimental charge radii, binding energies, and ground state spins, we conclude that the uncertainties in the ratios of weak charges are less than 10^(-3), comfortably smaller than the anticipated experimental error.

Sequence-function relationships in E. coli transcriptional regulation

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

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

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

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Measurement of the band bending and surface dipole at chemically functionalized Si(111)/vacuum interfaces

The core-level energy shifts observed using X-ray photoelectron spectroscopy (XPS) have been used to determine the band bending at Si(111) surfaces terminated with Si-Br, Si-H, and Si-CH₃ groups, respectively. The surface termination influenced the band bending, with the Si 2p_3/2 binding energy affected more by the surface chemistry than by the dopant type. The highest binding energies were measured on Si(111)-Br (whose Fermi level was positioned near the conduction band at the surface), followed by Si(111)-H, followed by Si(111)-CH₃ (whose Fermi level was positioned near mid-gap at the surface). Si(111)-CH₃ surfaces exposed to Br₂(g) yielded the lowest binding energies, with the Fermi level positioned between mid-gap and the valence band. The Fermi level position of Br₂(g)-exposed Si(111)-CH₃ was consistent with the presence of negatively charged bromine-containing ions on such surfaces. The binding energies of all of the species detected on the surface (C, O, Br) shifted with the band bending, illustrating the importance of isolating the effects of band bending when measuring chemical shifts on semiconductor surfaces. The influence of band bending was confirmed by surface photovoltage (SPV) measurements, which showed that the core levels shifted toward their flat-band values upon illumination. Where applicable, the contribution from the X-ray source to the SPV was isolated and quantified. Work functions were measured by ultraviolet photoelectron spectroscopy (UPS), allowing for calculation of the sign and magnitude of the surface dipole in such systems. The values of the surface dipoles were in good agreement with previous measurements as well as with electronegativity considerations. The binding energies of the adventitious carbon signals were affected by band bending as well as by the surface dipole. A model of band bending in which charged surface states are located exterior to the surface dipole is consistent with the XPS and UPS behavior of the chemically functionalized Si(111) surfaces investigated herein.

Biophysics of V(D)J recombination and Genome Packaging: In singulo studies on RAG, HMGB1, and TFAM

The recombination-activating gene products, RAG1 and RAG2, initiate V(D)J recombination during lymphocyte development by cleaving DNA adjacent to conserved recombination signal sequences (RSSs). The reaction involves DNA binding, synapsis, and cleavage at two RSSs located on the same DNA molecule and results in the assembly of antigen receptor genes. Since their discovery full-length, RAG1 and RAG2 have been difficult to purify, and core derivatives are shown to be most active when purified from adherent 293-T cells. However, the protein yield from adherent 293-T cells is limited. Here we develop a human suspension cell purification and change the expression vector to boost RAG production 6-fold. We use these purified RAG proteins to investigate V(D)J recombination on a mechanistic single molecule level. As a result, we are able to measure the binding statistics (dwell times and binding energies) of the initial RAG binding events with or without its co-factor high mobility group box protein 1 (HMGB1), and to characterize synapse formation at the single-molecule level yielding insights into the distribution of dwell times in the paired complex and the propensity for cleavage upon forming the synapse. We then go on to investigate HMGB1 further by measuring it compact single DNA molecules. We observed concentration dependent DNA compaction, differential DNA compaction depending on the divalent cation type, and found that at a particular HMGB1 concentration the percentage of DNA compacted is conserved across DNA lengths. Lastly, we investigate another HMGB protein called TFAM, which is essential for packaging the mitochondrial genome. We present crystal structures of TFAM bound to the heavy strand promoter 1 (HSP1) and to nonspecific DNA. We show TFAM dimerization is dispensable for DNA bending and transcriptional activation, but is required for mtDNA compaction. We propose that TFAM dimerization enhances mtDNA compaction by promoting looping of mtDNA.

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

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

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

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

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

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

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

Chapter Six contains the experimental documentation of the thesis work.

Unnatural amino acids in the synthesis and semisynthesis of metalloprotein motifs

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

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

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

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

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

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

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

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

Insights into the mechanism of human erythrocyte hexose transport : a transferred NOE study of glucose binding to GLUT1

This study examines binding of α- and β-D-glucose in their equilibrium mixture to the glucose transporter (GLUT1) in human erythrocyte membrane preparations by an ^1H NMR method, the transferred NOE (TRNOE). This method is shown theoretically and experimentally to be a sensitive probe of weak ligand-macromolecule interactions. The TRNOEs observed are shown to arise solely from glucose binding to GLUT1. Sites at both membrane faces contribute to the TRNOEs. Binding curves obtained are consistent with a homogeneous class of sugar sites, with an apparent K_D which varies (from ~30 mM to ~70 mM for both anomers) depending on the membrane preparation examined. Preparations with a higher proportion of the cytoplasmic membrane face exposed to bulk solution yield higher apparent KK_Ds. The glucose transport inhibitor cytochalasin B essentially eliminates the TRNOE. Nonlinearity was found in the dependence on sugar concentration of the apparent inhibition constant for cytochalasin B reversal of the TRNOE observed in the α anomer (and probably the β anomer); such nonlinearity implies the existence of ternary complexes of sugar, inhibitor and transporter. The inhibition results furthermore imply the presence of a class of relatively high-affinity (K_D < 2mM) sugar sites specific for the α anomer which do not contribute to NMR-observable binding. The presence of two classes of sugar-sensitive cytochalasin B sites is also indicated. These results are compared with predictions of the alternating conformer model of glucose transport. Variation of apparent K_D in the NMR-observable sites, the formation of ternary complexes and the presence of an anomer-specific site are shown to be inconsistent with this model. An alternate model is developed which reconciles these results with the known transport behavior of GLUT1. In this model, the transporter possesses (at minimum) three classes of sugar sites: (i) transport sites, which are alternately exposed to the cytoplasmic or the extracellular compartment, but never to both simultaneously, (ii) a class of sites (probably relatively low-affinity) which are confined to one compartment, and (iii) the high-affinity α anomer-specific sites, which are confined to the cytoplasmic compartment.

The relationship between the radio and gamma-ray emission of blazars

Blazars are active galaxies with a jet closely oriented to our line of sight. They are powerful, variable emitters from radio to gamma-ray wavelengths. Although the general picture of synchrotron emission at low energies and inverse Compton at high energies is well established, important aspects of blazars are not well understood. In particular, the location of the gamma-ray emission region is not clearly established, with some theories favoring a location close to the central engine, while others place it at parsec scales in the radio jet.

We developed a program to locate the gamma-ray emission site in blazars, through the study of correlated variations between their gamma-ray and radio-wave emission. Correlated variations are expected when there is a relation between emission processes at both bands, while delays tell us about the relative location of their energy generation zones. Monitoring at 15 GHz using the Owens Valley Radio Observatory 40 meter telescope started in mid-2007. The program monitors 1593 blazars twice per week, including all blazars detected by the Fermi Gamma-ray Space Telescope (Fermi) north of -20 degrees declination. This program complements the continuous monitoring of gamma-rays by Fermi.

Three year long gamma-ray light curves for bright Fermi blazars are cross-correlated with four years of radio monitoring. The significance of cross-correlation peaks is investigated using simulations that account for the uneven sampling and noise properties of the light curves, which are modeled as red-noise processes with a simple power-law power spectral density. We found that out of 86 sources with high quality data, only three show significant correlations (AO 0235+164, B2 2308+34 and PKS 1502+106). Additionally, we find a significant correlation for Mrk 421 when including the strong gamma-ray/radio flare of late 2012. In all four cases radio variations lag gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. For PKS 1502+106 we locate the gamma-ray emission site parsecs away from the central engine, thus disfavoring the model of Blandford and Levinson (1995), while other cases are inconclusive. These findings show that continuous monitoring over long time periods is required to understand the cross-correlation between gamma-ray and radio-wave variability in most blazars.

Binding site size limitations of imidazole-pyrrole polyamides for recognition in the minor groove of DNA

The discovery that the three ring polyamide Im-Py-Py-Dp containing imidazole (Im) and pyrrole (Py) carboxamides binds the DNA sequence 5'-(A,T)G(A,T)C(A,T)-3' as an antiparallel dimer offers a new model for the design of ligands for specific recognition of sequences in the minor groove containing both G,C and A,T base pairs. In Chapter 2, experiments are described in which the sequential addition of five N- methylpyrrolecarboxamides to the imidazole-pyrrole polyamide Im-Py-Py-Dp affords a series of six homologous polyamides, Im-(Py)_2-7-Dp, that differ in the size of their binding site, apparent first order binding affinity, and sequence specificity. These results demonstrate that DNA sequences up to nine base pairs in length can be specifically recognized by imidazole-pyrrole polyamides containing three to seven rings by 2:1 polyamide-DNA complex formation in the minor groove. Recognition of a nine base pair site defines the new lower limit of the binding site size that can be recognized by polyamides containing exclusively imidazole and pyrrolecarboxamides. The results of this study should provide useful guidelines for the design of new polyamides that bind longer DNA sites with enhanced affinity and specificity.

In Chapter 3 the design and synthesis of the hairpin polyamide Im-Py-Im-Py-γ-Im- Py-Im-Py-Dp is described. Quantitative DNase I footprint titration experiments reveal that Im-Py-Im-Py-γ-Im-Py-Im-Py-Dp binds six base pair 5'-(A,T)GCGC(A,T)-3' sequences with 30-fold higher affinity than the unlinked polyamide Im-Py-Im-Py-Dp. The hairpin polyamide does not discriminate between A•T and T•A at the first and sixth positions of the binding site as three sites 5'-TGCGCT-3', 5'-TGCGCA-3', and 5 'AGCGCT- 3' are bound with similar affinity. However, Im-Py-Im-Py-γ-Im-Py-Im-PyDp is specific for and discriminates between G•C and C•G base pairs in the 5'-GCGC-3' core as evidenced by lower affinities for the mismatched sites 5'-AACGCA-3', 5'- TGCGTT-3', 5'-TGCGGT-3', and 5'-ACCGCT-3'.

In Chapter 4, experiments are described in which a kinetically stable hexa-aza Schiff base La³⁺ complex is covalently attached to a Tat(49-72) peptide which has been shown to bind the HIV-1 TAR RNA sequence. Although these metallo-peptides cleave TAR site-specifically in the hexanucleotide loop to afford products consistent with hydrolysis, a series of control experiments suggests that the observed cleavage is not caused by a sequence-specifically bound Tat(49-72)-La(L)³⁺ peptide.

Functional evaluation of noncovalent interactions in neuroreceptors and progress toward the expansion of unnatural amino acid methodology

This dissertation primarily describes chemical-scale studies of G protein-coupled receptors and Cys-loop ligand-gated ion channels to better understand ligand binding interactions and the mechanism of channel activation using recently published crystal structures as a guide. These studies employ the use of unnatural amino acid mutagenesis and electrophysiology to measure subtle changes in receptor function.

In chapter 2, the role of a conserved aromatic microdomain predicted in the D3 dopamine receptor is probed in the closely related D2 and D4 dopamine receptors. This domain was found to act as a structural unit near the ligand binding site that is important for receptor function. The domain consists of several functionally important noncovalent interactions including hydrogen bond, aromatic-aromatic, and sulfur-π interactions that show strong couplings by mutant cycle analysis. We also assign an alternate interpretation for the linear fluorination plot observed at W6.48, a residue previously thought to participate in a cation-π interaction with dopamine.

Chapter 3 outlines attempts to incorporate chemically synthesized and in vitro acylated unnatural amino acids into mammalian cells. While our attempts were not successful, method optimizations and data for nonsense suppression with an in vivo acylated tRNA are included. This chapter is aimed to aid future researchers attempting unnatural amino acid mutagenesis in mammalian cells.

Chapter 4 identifies a cation-π interaction between glutamate and a tyrosine residue on loop C in the GluClβ receptor. Using the recently published crystal structure of the homologous GluClα receptor, other ligand-binding and protein-protein interactions are probed to determine the similarity between this invertebrate receptor and other more distantly related vertebrate Cys-loop receptors. We find that many of the interactions previously observed are conserved in the GluCl receptors, however care must be taken when extrapolating structural data.

Chapter 5 examines inherent properties of the GluClα receptor that are responsible for the observed glutamate insensitivity of the receptor. Chimera synthesis and mutagenesis reveal the C-terminal portion of the M4 helix and the C-terminus as contributing to formation of the decoupled state, where ligand binding is incapable of triggering channel gating. Receptor mutagenesis was unable to identify single residue mismatches or impaired protein-protein interactions within this domain. We conclude that M4 helix structure and/or membrane dynamics are likely the cause of ligand insensitivity in this receptor and that the M4 helix has an role important in the activation process.

Design, synthesis, and biological activity of rhodium metalloinsertors

Deficiencies in the mismatch repair (MMR) pathway are associated with several types of cancers, as well as resistance to commonly used chemotherapeutics. Rhodium metalloinsertors have been found to bind DNA mismatches with high affinity and specificity in vitro, and also exhibit cell-selective cytotoxicity, targeting MMR-deficient cells over MMR-proficient cells.

Here we examine the biological fate of rhodium metalloinsertors bearing dipyridylamine ancillary ligands. These complexes are shown to exhibit accelerated cellular uptake which permits the observation of various cellular responses, including disruption of the cell cycle and induction of necrosis, which occur preferentially in the MMR-deficient cell line. These cellular responses provide insight into the mechanisms underlying the selective activity of this novel class of targeted anti-cancer agents.

In addition, ten distinct metalloinsertors with varying lipophilicities are synthesized and their mismatch binding affinities and biological activities studied. While they are found to have similar binding affinities, their cell-selective antiproliferative and cytotoxic activities vary significantly. Inductively coupled plasma mass spectrometry (ICP-MS) experiments show that all of these metalloinsertors localize in the nucleus at sufficient concentrations for binding to DNA mismatches. Furthermore, metalloinsertors with high rhodium localization in the mitochondria show toxicity that is not selective for MMR-deficient cells. This work supports the notion that specific targeting of the metalloinsertors to nuclear DNA gives rise to their cytotoxic and antiproliferative activities that are selective for cells deficient in MMR.

To explore further the basis of the unique selectivity of the metlloinsertors in targeting MMR-deficient cells, experiments were conducted using engineered NCI-H23 lung adenocarcinoma cells that contain a doxycycline-inducible shRNA which suppresses the expression of the MMR gene MLH1. Here we use this new cell line to further validate rhodium metalloinsertors as compounds capable of differentially inhibiting the proliferation of MMR-deficient cancer cells over isogenic MMR-proficient cells. General DNA damaging agents, such as cisplatin and etoposide, in contrast, are less effective in the induced cell line defective in MMR.

Finally, we describe a new subclass of metalloinsertors with enhanced potency and selectivity, in which the complexes show Rh-O coordination. In particular, it has been found that both Δ and Λ enantiomers of [Rh(chrysi)(phen)(DPE)]²⁺ bind to DNA with similar affinities, suggesting a possible different binding conformation than previous metalloinsertors. Remarkably, all members of this new family of compounds have significantly increased potency in a range of cellular assays; indeed, all are more potent than the FDA-approved anticancer drugs cisplatin and MNNG. Moreover, these activities are coupled with high levels of selectivity for MMR-deficient cells.