Thermodynamics of the Hydroxyl Radical Addition to Isoprene

Oxidation of isoprene by the hydroxyl radical leads to tropospheric ozone formation. Consequently, a more complete understanding of this reaction could lead to better models of regional air quality, a better understanding of aerosol formation, and a better understanding of reaction kinetics and dynamics. The most common first step in the oxidation of isoprene is the formation of an adduct, with the hydroxyl radical adding to one of four unsaturated carbon atoms in isoprene. In this paper, we discuss how the initial conformations of isoprene, s-trans and s-gauche, influences the pathways to adduct formation. We explore the formation of pre-reactive complexes at low and high temperatures, which are often invoked to explain the negative temperature dependence of this reaction’s kinetics. We show that at higher temperatures the free energy surface indicates that a pre-reactive complex is unlikely, while at low temperatures the complex exists on two reaction pathways. The theoretical results show that at low temperatures all eight pathways possess negative reaction barriers, and reaction energies that range from −36.7 to −23.0 kcal·mol−1. At temperatures in the lower atmosphere, all eight pathways possess positive reaction barriers that range from 3.8 to 6.0 kcal·mol−1 and reaction energies that range from −28.8 to −14.4 kcal·mol−1.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1 clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Comparison of CBS-QB3, CBS-APNO, G2, and G3 thermochemical predictions with experiment for formation of ionic clusters of hydronium and hydroxide ions complexed with water

The GAUSSIAN 2, GAUSSIAN 3, complete basis set-QB3, and complete basis set-APNO methods have been used to calculate ΔH∘ and ΔG∘ values for ionic clusters of hydronium and hydroxide ions complexed with water. Results for the clusters H3O+(H2O)n andOH−(H2O)n, where n=1–4 are reported in this paper, and compared against experimental values contained in the National Institutes of Standards and Technology (NIST) database. Agreement with experiment is excellent for the three ab initio methods for formation of these clusters. The high accuracy of these methods makes them reliable for calculating energetics for the formation of ionic clusters containing water. In addition this allows them to serve as a valuable check on the accuracy of experimental data reported in the NIST database, and makes them useful tools for addressing unresolved issues in atmospheric chemistry.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Annulation of Tetrathiafulvalene to the Bay Region of Perylenediimide: Fast Electron-Transfer Processes in Polar and Nonpolar Solvents

A tetrathiafulvalene donor has been annulated to the bay region of perylenediimide through a 1H-benzo-[d]pyrrolo[1,2-a]imidazol-1-one spacer affording an extended pi-conjugated molecular dyad (TTF-PDI). To gain insight into its ground- and excited-state electronic properties, the reference compound Ph-PDI has been prepared via a direct Schiff-base condensation of N,N'-bis(1-octylnonyl) benzoperylene-1',2':3,4:9,10-hexacarboxylic-1',2'-anhydride-3,4:9,10-bis (imide) with benzene-1,2-diamine. Both the experimental and the computational (DFT) results indicate that TTF-PDI exhibits significant intramolecular electronic interactions giving rise to an efficient photoinduced charge-separation process. Free-energy calculations verify that the process from TTF to the singlet-excited state of PDI is exothermic in both polar and nonpolar solvents. Fast adiabatic electron-transfer processes of a compactly fused, pi-conjugated TTF-PDI dyad in benzonitrile, 2-methyltetrahydrofuran, anisole and toluene were observed by femtosecond transient absorption spectral measurements. The lifetimes of radical-ion pairs slightly increase with decreasing the solvent polarities, suggesting that the charge-recombination occurs in the Marcus inverted region. By utilizing the nanosecond transient absorption technique, the intermolecular electron-transfer process in a mixture of has been observed via the triplet excited PDI for the first time.

Engineering Tocopherol Selectivity in alpha-TTP: A Combined In Vitro/In Silico Study

"We present a combined in vitro/in silico study to determine the molecular origin of the selectivity of a-tocopherol transfer" "protein (a-TTP) towards a-tocopherol. Molecular dynamics simulations combined to free energy perturbation calculations predict a binding free energy for a-tocopherol to a-TTP 8.26+2.13 kcal mol{1 lower than that of c-tocopherol. Our calculations show that c-tocopherol binds to a-TTP in a significantly distorted geometry as compared to that of the natural ligand. Variations in the hydration of the binding pocket and in the protein structure are found as well. We propose a mutation, A156L, which significantly modifies the selectivity properties of a-TTP towards the two tocopherols. In particular, our simulations predict that A156L binds preferentially to c-tocopherol, with striking structural similarities to the wild-type- a-tocopherol complex. The affinity properties are confirmed by differential scanning fluorimetry as well as in vitro competitive binding assays. Our data indicate that residue A156 is at a critical position for determination of the selectivity of a-TTP. The engineering of TTP mutants with modulating binding properties can have potential impact at industrial level for easier purification of single tocopherols from vitamin E mixtures coming from natural oils or synthetic processes. Moreover," "the identification of a c-tocopherol selective TTP offers the possibility to challenge the hypotheses for the evolutionary development of a mechanism for a-tocopherol selection in omnivorous animals."

Wetting and spreading of a surfactant film on solid particles: influence of sharp edges and surface irregularities

In addition to particle size and surface chemistry, the shape of particles plays an important role in their wetting and displacement by the surfactant film in the lung. The role of particle shape was the subject of our investigations using a model system consisting of a modified Langmuir-Wilhelmy surface balance. We measured the influence of sharp edges (lines) and other highly curved surfaces, including sharp corners or spikes, of different particles on the spreading of a dipalmitoylphosphatidyl (DPPC) film. The edges of cylindrical sapphire plates (circular curved edges, 1.65 mm radius) were wetted at a surface tension of 10.7 mJ/m2 (standard error (SE) = 0.45, n = 20) compared with that of 13.8 mJ/m2 (SE = 0.20, n = 20) for cubic sapphire plates (straight linear edges, edge length 3 mm) (p < 0.05). The top surfaces of the sapphire plates (cubic and cylindrical) were wetted at 8.4 mJ/m2 (SE = 0.54, n = 20) and 9.1 mJ/m2 (SE = 0.50, n = 20), respectively, but the difference was not significant (p > 0.05). The surfaces of the plates showed significantly higher resistance to spreading compared to that of the edges, as substantially lower surface tensions were required to initiate wetting (p < 0.05). Similar results were found for talc particles, were the edges of macro- and microcrystalline particles were wetted at 7.2 mJ/m2 (SE = 0.52, n = 20) and 8.2 mJ/m2 (SE = 0.30, n = 20) (p > 0.05), respectively, whereas the surfaces were wetted at 3.8 mJ/m2 (SE = 0.89, n = 20) and 5.8 mJ/m2 (SE = 0.52, n = 20) (p < 0.05), respectively. Further experiments with pollen of malvaceae and maize (spiky and fine knobbly surfaces) were wetted at 10.0 mJ/m2 (SE = 0.52, n = 10) and 22.75 mJ/m2 (SE = 0.81, n = 10), respectively (p < 0.05). These results show that resistance to spreading of a DPPC film on various surfaces is dependent on the extent these surfaces are curved. This is seen with cubic sapphire plates which have at their corners a radius of curvature of about 0.75 microm, spiky malvaceae pollen with an even smaller radius on top of their spikes, or talc with various highly curved surfaces. These highly curved surfaces resisted wetting by the DPPC film to a higher degree than more moderately curved surfaces such as those of cylindrical sapphire plates, maize pollens, or polystyrene spheres, which have a surface free energy similar to that of talc but a smooth surface. The macroscopic plane surfaces of the particles demonstrated the greatest resistance to spreading. This was explained by the extremely fine grooves in the nanometer range, as revealed by electron microscopy. In summary, to understand the effects of airborne particles retained on the surfaces of the respiratory tract, and ultimately their pathological potential, not only the particle size and surface chemistry but also the particle shape should be taken in consideration.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

Optical hysteresis of ridge waveguide magnetic garnet films and cavities in two-dimensional magneto-photonic crystal slabs

Magnetic iron garnets as well as magnetic photonic crystals are of great interests in magneto-optic applications such as isolators, current captors, circulators, TE-TM mode conversion, wavelength accordable filters, optical sensors and switches, all of which provide a promising platform for future integrated optical circuits. In the present work, two topics are studied based on magnetic iron garnet films. In the first part, the characteristics of the magnetization are investigated for ridge waveguides fabricated on (100) oriented iron garnet thin films. The magnetic response in magneto-optic waveguides patterned on epitaxial magnetic garnet films depends on the crystallographic orientation of the waveguides and the magnetic anisotropy of the material. These can be studied by polarization rotation hysteresis loops, which are related to the component of magnetization parallel to the light propagation direction and the linear birefringence. Polarization rotation hysteresis loops for low birefringence waveguides with different orientations are experimentally investigated. Asymmetric stepped curves are obtained from waveguides along, due to the large magnetocrystalline anisotropy in the plane. A model based on the free energy density is developed to demonstrate the motion of the magnetization and can be used in the design of magneto-optic devices. The second part of this thesis focuses on the design and fabrication of high-Q cavities in two-dimensional magneto-photonic crystal slabs. The device consists of a layer of silicon and a layer of iron garnet thin film. Triangular lattice elliptical air holes are patterned in the slab. The fundamental TM band gap overlaps with the first-order TE band gap from 0374~0.431(a/λ) showing that both TE and TM polarization light can be confined in the photonic crystals. A nanocavity is designed to obtain both TE and TM defect modes in the band gaps. Additional work is needed to overlap the TE and TM defect modes and obtain a high-Q cavity so as to develop miniaturized Faraday rotators.

How fluids bend: the elastic expansion for higher-dimensional black holes

Hydrodynamics can be consistently formulated on surfaces of arbitrary co-dimension in a background space-time, providing the effective theory describing long-wavelength perturbations of black branes. When the co-dimension is non-zero, the system acquires fluid-elastic properties and constitutes what is called a fluid brane. Applying an effective action approach, the most general form of the free energy quadratic in the extrinsic curvature and extrinsic twist potential of stationary fluid brane configurations is constructed to second order in a derivative expansion. This construction generalizes the Helfrich-Canham bending energy for fluid membranes studied in theoretical biology to the case in which the fluid is rotating. It is found that stationary fluid brane configurations are characterized by a set of 3 elastic response coefficients, 3 hydrodynamic response coefficients and 1 spin response coefficient for co-dimension greater than one. Moreover, the elastic degrees of freedom present in the system are coupled to the hydrodynamic degrees of freedom. For co-dimension-1 surfaces we find a 8 independent parameter family of stationary fluid branes. It is further shown that elastic and spin corrections to (non)-extremal brane effective actions can be accounted for by a multipole expansion of the stress-energy tensor, therefore establishing a relation between the different formalisms of Carter, Capovilla-Guven and Vasilic-Vojinovic and between gravity and the effective description of stationary fluid branes. Finally, it is shown that the Young modulus found in the literature for black branes falls into the class predicted by this approach - a relation which is then used to make a proposal for the second order effective action of stationary blackfolds and to find the corrected horizon angular velocity of thin black rings.

RNA duplexes with biphenyl substituents as base replacements are less stable than DNA duplexes

Introduction of a hydrophobic biphenyl-C-nucleotide pair into a 11-mer RNA duplex is associated with a net penalty in the free energy of duplex formation of 2.0 kcal mol(-1) or 10 degrees C in T(m), relative to DNA. These differential stabilities are of relevance with respect to the transcriptional and translational aspects of hydrophobic base-pairs

Tracking individual membrane proteins and their biochemistry: The power of direct observation

The advent of single molecule fluorescence microscopy has allowed experimental molecular biophysics and biochemistry to transcend traditional ensemble measurements, where the behavior of individual proteins could not be precisely sampled. The recent explosion in popularity of new super-resolution and super-localization techniques coupled with technical advances in optical designs and fast highly sensitive cameras with single photon sensitivity and millisecond time resolution have made it possible to track key motions, reactions, and interactions of individual proteins with high temporal resolution and spatial resolution well beyond the diffraction limit. Within the purview of membrane proteins and ligand gated ion channels (LGICs), these outstanding advances in single molecule microscopy allow for the direct observation of discrete biochemical states and their fluctuation dynamics. Such observations are fundamentally important for understanding molecular-level mechanisms governing these systems. Examples reviewed here include the effects of allostery on the stoichiometry of ligand binding in the presence of fluorescent ligands; the observation of subdomain partitioning of membrane proteins due to microenvironment effects; and the use of single particle tracking experiments to elucidate characteristics of membrane protein diffusion and the direct measurement of thermodynamic properties, which govern the free energy landscape of protein dimerization. The review of such characteristic topics represents a snapshot of efforts to push the boundaries of fluorescence microscopy of membrane proteins to the absolute limit.

Probing the moduli dependence of refined topological amplitudes

With the aim of providing a worldsheet description of the refined topological string, we continue the study of a particular class of higher derivative couplings Fg,n in the type II string effective action compactified on a Calabi–Yau threefold. We analyse first order differential equations in the anti-holomorphic moduli of the theory, which relate the Fg,n to other component couplings. From the point of view of the topological theory, these equations describe the contribution of non-physical states to twisted correlation functions and encode an obstruction for interpreting the Fg,n as the free energy of the refined topological string theory. We investigate possibilities of lifting this obstruction by formulating conditions on the moduli dependence under which the differential equations simplify and take the form of generalised holomorphic anomaly equations. We further test this approach against explicit calculations in the dual heterotic theory.

Solvation-Driven Charge Transfer and Localization in Metal Complexes

In any physicochemical process in liquids, the dynamical response of the solvent to the solutes out of equilibrium plays a crucial role in the rates and products: the solvent molecules react to the changes in volume and electron density of the solutes to minimize the free energy of the solution, thus modulating the activation barriers and stabilizing (or destabilizing) intermediate states. In charge transfer (CT) processes in polar solvents, the response of the solvent always assists the formation of charge separation states by stabilizing the energy of the localized charges. A deep understanding of the solvation mechanisms and time scales is therefore essential for a correct description of any photochemical process in dense phase and for designing molecular devices based on photosensitizers with CT excited states. In the last two decades, with the advent of ultrafast time-resolved spectroscopies, microscopic models describing the relevant case of polar solvation (where both the solvent and the solute molecules have a permanent electric dipole and the mutual interaction is mainly dipole−dipole) have dramatically progressed. Regardless of the details of each model, they all assume that the effect of the electrostatic fields of the solvent molecules on the internal electronic dynamics of the solute are perturbative and that the solvent−solute coupling is mainly an electrostatic interaction between the constant permanent dipoles of the solute and the solvent molecules. This well-established picture has proven to quantitatively rationalize spectroscopic effects of environmental and electric dynamics (time-resolved Stokes shifts, inhomogeneous broadening, etc.). However, recent computational and experimental studies, including ours, have shown that further improvement is required. Indeed, in the last years we investigated several molecular complexes exhibiting photoexcited CT states, and we found that the current description of the formation and stabilization of CT states in an important group of molecules such as transition metal complexes is inaccurate. In particular, we proved that the solvent molecules are not just spectators of intramolecular electron density redistribution but significantly modulate it. Our results solicit further development of quantum mechanics computational methods to treat the solute and (at least) the closest solvent molecules including the nonperturbative treatment of the effects of local electrostatics and direct solvent−solute interactions to describe the dynamical changes of the solute excited states during the solvent response.

Role of phosphatidylethanolamine in the function and assembly of phenylalanine-specific permease of Escherichia coli

Transmembrane segments of polytopic membrane proteins once inserted are generally considered stably oriented due to the large free energy barrier for topological reorientation of adjacent extra-membrane domains. However, proper topology and function of the polytopic membrane protein lactose permease (LacY) of Escherichia coli is dependent on the membrane phospholipid composition revealing topological dynamics of transmembrane domains (Bogdanov, M., Heacock, P. N., and Dowhan, W. (2002) EMBO J. 21, 2107–2116). The high affinity phenylalanine permease PheP shares many topological similarities with LacY. In this study, mutant E. coli cells lacking phosphatidylethanolamine (PE) as a membrane component were used to evaluate the role of PE in the function and assembly of PheP. Active transport of phenylalanine by cells lacking PE was severely inhibited (both Vmax and Km were altered), whereas the PheP protein level in membranes was unaffected. Cysteine residues were introduced into predicted periplasmic or cytoplasmic segments of cysteine-less PheP, and the topology of the protein was explored using a membrane-impermeable thiol-specific biotinylated probe. Based on the biotinylation patterns of PheP in whole cells, the N-terminus and adjoining transmembrane hairpin of PheP adopted an inverted topological orientation in PE-lacking cells. Introduction of PE following the assembly of PheP triggered a reorientation of the N-terminus and adjacent hairpin to their native orientation associated with regain of wild type transport function. These results coupled with the results for LacY support a specific role for membrane lipid composition in determining topological organization and function of membrane proteins. Several other secondary symporters are compromised for activity in PE-lacking cells suggesting that lipid-assisted topogenesis is a general property of such transporters. The reversible orientation of these secondary transport proteins in response to a change of phospholipid composition might be a result of inherent conformational flexibility necessary for transport function or during protein assembly. ^