Study of pK Values and Effective Dielectric Constants of Ionizable Residues in Pentapeptides and in Staphylococcal Nuclease (SNase) Using a Mean-Field Approach

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Adsorção de polieletrólitos em superfícies opostamente carregadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

Low-dimensional carbon allotropes: an electron microscopy investigation

The research reported in this manuscript concerns the structural characterization of graphene membranes and single-walled carbon nanotubes (SWCNTs). The experimental investigation was performed using a wide range of transmission electron microscopy (TEM) techniques, from conventional imaging and diffraction, to advanced interferometric methods, like electron holography and Geometric Phase Analysis (GPA), using a low-voltage optical set-up, to reduce radiation damage in the samples. Electron holography was used to successfully measure the mean electrostatic potential of an isolated SWCNT and that of a mono-atomically thin graphene crystal. The high accuracy achieved in the phase determination, made it possible to measure, for the first time, the valence-charge redistribution induced by the lattice curvature in an individual SWCNT. A novel methodology for the 3D reconstruction of the waviness of a 2D crystal membrane has been developed. Unlike other available TEM reconstruction techniques, like tomography, this new one requires processing of just a single HREM micrograph. The modulations of the inter-planar distances in the HREM image are measured using Geometric Phase Analysis, and used to recover the waviness of the crystal. The method was applied to the case of a folded FGC, and a height variation of 0.8 nm of the surface was successfully determined with nanometric lateral resolution. The adhesion of SWCNTs to the surface of graphene was studied, mixing shortened SWCNTs of different chiralities and FGC membranes. The spontaneous atomic match of the two lattices was directly imaged using HREM, and we found that graphene membranes act as tangential nano-sieves, preferentially grafting achiral tubes to their surface.

A planar Penning trap

In this thesis I present theoretical and experimental results concern- ing the operation and properties of a new kind of Penning trap, the planar trap. It consists of circular electrodes printed on an isolating surface, with an homogeneous magnetic field pointing perpendicular to that surface. The motivation of such geometry is to be found in the construction of an array of planar traps for quantum informa- tional purposes. The open access to radiation of this geometry, and the long coherence times expected for Penning traps, make the planar trap a good candidate for quantum computation. Several proposals for quantum 2-qubit interactions are studied and estimates for their rates are given. An expression for the electrostatic potential is presented, and its fea- tures exposed. A detailed study of the anharmonicity of the potential is given theoretically and is later demonstrated by experiment and numerical simulations, showing good agreement. Size scalability of this trap has been studied by replacing the original planar trap by a trap twice smaller in the experimental setup. This substitution shows no scale effect apart from those expected for the scaling of the parameters of the trap. A smaller lifetime for trapped electrons is seen for this smaller trap, but is clearly matched to a bigger misalignment of the trap’s surface and the magnetic field, due to its more difficult hand manipulation. I also give a hint that this trap may be of help in studying non-linear dynamics for a sextupolarly perturbed Penning trap.

Studies of degradation of organic solar cell materials using electrical scanning probe microscopy

Intense research is being done in the field of organic photovoltaics in order to synthesize low band-gap organic molecules. These molecules are electron donors which feature in combination with acceptor molecules, typically fullerene derivarntives, forming an active blend. This active blend has phase separated bicontinuous morphology on a nanometer scale. The highest recorded power conversionrnefficiencies for such cells have been 10.6%. Organic semiconductors differ from inorganic ones due to the presence of tightly bonded excitons (electron-hole pairs)resulting from their low dielectric constant (εr ≈2-4). An additional driving force is required to separate such Frenkel excitons since their binding energy (0.3-1 eV) is too large to be dissociated by an electric field alone. This additional driving force arises from the energy difference between the lowest unoccupied molecular orbital (LUMO) of the donor and the acceptor materials. Moreover, the efficiency of the cells also depends on the difference between the highest occupied molecular orbital (HOMO) of the donor and LUMO of the acceptor. Therefore, a precise control and estimation of these energy levels are required. Furthermore any external influences that change the energy levels will cause a degradation of the power conversion efficiency of organic solar cell materials. In particular, the role of photo-induced degradation on the morphology and electrical performance is a major contribution to degradation and needs to be understood on a nanometer scale. Scanning Probe Microscopy (SPM) offers the resolution to image the nanometer scale bicontinuous morphology. In addition SPM can be operated to measure the local contact potential difference (CPD) of materials from which energy levels in the materials can be derived. Thus SPM is an unique method for the characterization of surface morphology, potential changes and conductivity changes under operating conditions. In the present work, I describe investigations of organic photovoltaic materials upon photo-oxidation which is one of the major causes of degradation of these solar cell materials. SPM, Nuclear Magnetic Resonance (NMR) and UV-Vis spectroscopy studies allowed me to identify the chemical reactions occurring inside the active layer upon photo-oxidation. From the measured data, it was possible to deduce the energy levels and explain the various shifts which gave a better understanding of the physics of the device. In addition, I was able to quantify the degradation by correlating the local changes in the CPD and conductivity to the device characteristics, i.e., open circuit voltage and short circuit current. Furthermore, time-resolved electrostatic force microscopy (tr-EFM) allowed us to probe dynamic processes like the charging rate of the individual donor and acceptor domains within the active blend. Upon photo-oxidation, it was observed, that the acceptor molecules got oxidized first preventing the donor polymer from degrading. Work functions of electrodes can be tailored by modifying the interface with monomolecular thin layers of molecules which are made by a chemical reaction in liquids. These modifications in the work function are particularly attractive for opto-electronic devices whose performance depends on the band alignment between the electrodes and the active material. In order to measure the shift in work function on a nanometer scale, I used KPFM in situ, which means in liquids, to follow changes in the work function of Au upon hexadecanethiol adsorption from decane. All the above investigations give us a better understanding of the photo-degradation processes of the active material at the nanoscale. Also, a method to compare various new materials used for organic solar cells for stability is proposed which eliminates the requirement to make fully functional devices saving time and additional engineering efforts.

Modeling and simulation of fibrinogen and its adsorption behavior

In this thesis different approaches for the modeling and simulation of the blood protein fibrinogen are presented. The approaches are meant to systematically connect the multiple time and length scales involved in the dynamics of fibrinogen in solution and at inorganic surfaces. The first part of the thesis will cover simulations of fibrinogen on an all atom level. Simulations of the fibrinogen protomer and dimer are performed in explicit solvent to characterize the dynamics of fibrinogen in solution. These simulations reveal an unexpectedly large and fast bending motion that is facilitated by molecular hinges located in the coiled-coil region of fibrinogen. This behavior is characterized by a bending and a dihedral angle and the distribution of these angles is measured. As a consequence of the atomistic detail of the simulations it is possible to illuminate small scale behavior in the binding pockets of fibrinogen that hints at a previously unknown allosteric effect. In a second step atomistic simulations of the fibrinogen protomer are performed at graphite and mica surfaces to investigate initial adsorption stages. These simulations highlight the different adsorption mechanisms at the hydrophobic graphite surface and the charged, hydrophilic mica surface. It is found that the initial adsorption happens in a preferred orientation on mica. Many effects of practical interest involve aggregates of many fibrinogen molecules. To investigate such systems, time and length scales need to be simulated that are not attainable in atomistic simulations. It is therefore necessary to develop lower resolution models of fibrinogen. This is done in the second part of the thesis. First a systematically coarse grained model is derived and parametrized based on the atomistic simulations of the first part. In this model the fibrinogen molecule is represented by 45 beads instead of nearly 31,000 atoms. The intra-molecular interactions of the beads are modeled as a heterogeneous elastic network while inter-molecular interactions are assumed to be a combination of electrostatic and van der Waals interaction. A method is presented that determines the charges assigned to beads by matching the electrostatic potential in the atomistic simulation. Lastly a phenomenological model is developed that represents fibrinogen by five beads connected by rigid rods with two hinges. This model only captures the large scale dynamics in the atomistic simulations but can shed light on experimental observations of fibrinogen conformations at inorganic surfaces.

Correlation between Accurate Electron Density and Linear Optical Properties in Amino Acid Derivatives: L-Histidinium Hydrogen Oxalate

The accurate electron density and linear optical properties of L-histidinium hydrogen oxalate are discussed. Two high-resolution single crystal X-ray diffraction experiments were performed and compared with density functional calculations in the solid state as well as in the gas phase. The crystal packing and the hydrogen bond network are accurately investigated using topological analysis based on quantum theory of atoms in molecules, Hirshfeld surface analysis, and electrostatic potential mapping. The refractive indices are computed from couple perturbed Kohn-Sham calculations and measured experimentally. Moreover, distributed atomic polarizabilities are used to analyze the origin of the linear susceptibility in the crystal, in order to separate molecular and intermolecular causes. The optical properties are also correlated with the electron density distribution. This compound also offers the possibility to test the electron density building block approach for material science and different refinement schemes for accurate positions and displacement parameters of hydrogen atoms, in the absence of neutron diffraction data.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Electrostatic interaction between helical macromolecules in dense aggregates: An impetus for DNA poly- and meso-morphism

DNA exhibits a surprising multiplicity of structures when it is packed into dense aggregates. It undergoes various polymorphous transitions (e.g., from the B to A form) and mesomorphous transformations (from hexagonal to orthorhombic or monoclinic packing, changes in the mutual alignment of nearest neighbors, etc). In this report we show that such phenomena may have their origin in the specific helical symmetry of the charge distribution on DNA surface. Electrostatic interaction between neighboring DNA molecules exhibits strong dependence on the patterns of molecular surface groups and adsorbed counter-ions. As a result, it is affected by such structural parameters as the helical pitch, groove width, the number of base pairs per helical turn, etc. We derive expressions which relate the energy of electrostatic interaction with these parameters and with the packing variables characterizing the axial and azimuthal alignment between neighboring macromolecules. We show, in particular, that the structural changes upon the B-to-A transition reduce the electrostatic energy by ≈kcal/mol per base pair, at a random adsorption of counter ions. Ion binding into the narrow groove weakens or inverts this effect, stabilizing B-DNA, as it is presumably the case in Li+-DNA assemblies. The packing symmetry and molecular alignment in DNA aggregates are shown to be affected by the patterns of ion binding.

Solution structure of the antiapoptotic protein bcl-2

The structures of two isoforms of Bcl-2 that differ by two amino acids have been determined by NMR spectroscopy. Because wild-type Bcl-2 behaved poorly in solution, the structures were determined by using Bcl-2/Bcl-xL chimeras in which part of the putative unstructured loop of Bcl-2 was replaced with a shortened loop from Bcl-xL. These chimeric proteins have a low pI compared with the wild-type protein and are soluble. The structures of the two Bcl-2 isoforms consist of 6 α-helices with a hydrophobic groove on the surface similar to that observed for the homologous protein, Bcl-xL. Comparison of the Bcl-2 structures to that of Bcl-xL shows that although the overall fold is the same, there are differences in the structural topology and electrostatic potential of the binding groove. Although the structures of the two isoforms of Bcl-2 are virtually identical, differences were observed in the ability of the proteins to bind to a 25-residue peptide from the proapoptotic Bad protein and a 16-residue peptide from the proapoptotic Bak protein. These results suggest that there are subtle differences in the hydrophobic binding groove in Bcl-2 that may translate into differences in antiapoptotic activity for the two isoforms.

Computation of Surface Electrical Potentials of Plant Cell Membranes : Correspondence to Published Zeta Potentials from Diverse Plant Sources

A Gouy-Chapman-Stern model has been developed for the computation of surface electrical potential (ψ0) of plant cell membranes in response to ionic solutes. The present model is a modification of an earlier version developed to compute the sorption of ions by wheat (Triticum aestivum L. cv Scout 66) root plasma membranes. A single set of model parameters generates values for ψ0 that correlate highly with published ζ potentials of protoplasts and plasma membrane vesicles from diverse plant sources. The model assumes ion binding to a negatively charged site (R− = 0.3074 μmol m−2) and to a neutral site (P0 = 2.4 μmol m−2) according to the reactions R− + IΖ ⇌ RIΖ−1 and P0 + IΖ ⇌ PIΖ, where IΖ represents an ion of charge Ζ. Binding constants for the negative site are 21,500 m−1 for H+, 20,000 m−1 for Al3+, 2,200 m−1 for La3+, 30 m−1 for Ca2+ and Mg2+, and 1 m−1 for Na+ and K+. Binding constants for the neutral site are 1/180 the value for binding to the negative site. Ion activities at the membrane surface, computed on the basis of ψ0, appear to determine many aspects of plant-mineral interactions, including mineral nutrition and the induction and alleviation of mineral toxicities, according to previous and ongoing studies. A computer program with instructions for the computation of ψ0, ion binding, ion concentrations, and ion activities at membrane surfaces may be requested from the authors.

Biochemical characterization and DNA binding properties of the MocR family member GabR, a PLP-dependent transcription factor

GabR è un fattore di trascrizione chimerico appartenente alla famiglia dei MocR/GabR, costituito da un dominio N-terminale elica-giro-elica di legame al DNA e un dominio effettore e/o di oligomerizzazione al C-terminale. I due domini sono connessi da un linker flessibile di 29 aminoacidi. Il dominio C-terminale è strutturalmente omologo agli enzimi aminotransferasici fold-type I, i quali, utilizzando il piridossal-5’-fosfato (PLP) come cofattore, sono direttamente coinvolti nel metabolismo degli aminoacidi. L’interazione contemporanea di PLP e acido γ-aminobutirrico (GABA) a GabR fa sì che questa promuova la trascrizione di due geni, gabT e gabD, implicati nel metabolismo del GABA. GabR cristallizza come un omodimero con una configurazione testa-coda. Il legame con la regione promotrice gabTD avviene attraverso il riconoscimento specifico di due sequenze dirette e ripetute (ATACCA), separate da uno spacer di 34 bp. In questo studio sono state indagate le proprietà biochimiche, strutturali e di legame al DNA della proteina GabR di Bacillus subtilis. L’analisi spettroscopica dimostra che GabR interagisce con il PLP formando l’aldimina interna, mentre in presenza di GABA si ottiene l’aldimina esterna. L’interazione fra il promotore gabTD e le forme holo e apo di GabR è stata monitorata mediante Microscopia a Forza atomica (AFM). In queste due condizioni di legame è stata stimata una Kd di circa 40 ηM. La presenza di GABA invece, determinava un incremento di circa due volte della Kd, variazioni strutturali nei complessi GabR-DNA e una riduzione del compattamento del DNA alla proteina, indipendentemente dalla sequenza del promotore in esame. Al fine di valutare il ruolo delle caratteristiche topologiche del promotore, sono state inserite cinque e dieci bp all’interno della regione spacer che separa le due sequenze ripetute dirette riconosciute da GabR. I significativi cambiamenti topologici riscontrati nel frammento aggiunto di cinque bp si riflettono anche sulla forte riduzione dell’affinità di legame verso la proteina. Al contrario, l’inserzione di 10 bp provoca solamente l’allontanamento delle sequenze ripetute dirette. L’assenza quindi di cambiamenti significativi nella topologia di questo promotore fa sì che l’affinità di legame per GabR rimanga pressoché inalterata rispetto al promotore non mutato. L’analisi del potenziale elettrostatico superficiale di GabR mostra la presenza di una fascia carica positivamente che si estende lungo un’intera faccia della proteina. Per verificare l’importanza di questa caratteristica di GabR nel meccanismo di interazione al DNA, sono stati preparati ed indagati i mutanti R129Q e K362-366Q, in cui la carica positiva superficiale risultava indebolita. L’affinità di legame dei mutanti di GabR per il DNA era inferiore rispetto alla proteina non mutata, in particolar modo nel mutante K362-366Q. Le evidenze acquisite suggeriscono che la curvatura intrinseca del promotore ed il corretto orientamento delle sequenze sulla doppia elica, più della distanza che le separa, siano critici per sostenere l’interazione con GabR. Oltre a questo, la superficie positiva di GabR è richiesta per accomodare la curvatura del DNA sul corpo della proteina. Alla luce di questo, l’interazione GabR-gabTD è un esempio di come il riconoscimento specifico di sequenze, la topologia del DNA e le caratteristiche strutturali della proteina siano contemporaneamente necessarie per sostenere un’interazione proteina-DNA stabile.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.

Zastosowanie modelu hantli w badaniach właściwości fizykochemicznych i strukturalnych układów homo- i heterogenicznych metodą Monte Carlo

Wydział Chemii