Neue Methoden zur Computersimulation von Polymersystemen auf verschiedenen Längenskalen und ihre Anwendungen

In dieser Dissertation stellen wir einen neuen Ansatz zurModellierungvon Polymersystemen vor. Es werden (von methodischer Seiteher) zweiautomatisierte Iterationschemata dazu eingeführt,Kraftfeldparametermesoskopischer Polymersysteme systematisch zu optimieren:DasSimplex-Verfahren und das Struktur-Differenzen-Verfahren. Sowerdendiejenigen Freiheitsgrade aus Polymersystemen eliminiert,die einehohe Auflösung erfordern, was die Modellierung größerersystemeermöglicht. Nach Tests an einfachen Flüssigkeiten werdenvergröberteModelle von drei prototypischen Polymeren (Polyacrylsäure,Polyvinylalkohol und Polyisopren) in unterschiedlichenUmgebungen(gutes Lösungsmittel und Schmelze) entwickelt und ihrVerhalten aufder Mesoskala ausgiebig geprüft. Die zugehörige Abbildung(vonphysikalischer Seite her) so zu gestalten, daß sie dieunverwechselbaren Charakteristiken jedes systems auf diemesoskopischeLängenskala überträgt, stellt eine entscheidende Anforderungan dieautomatisierten Verfahren dar. Unsere Studien belegen, daß mesoskopische Kraftfeldertemperatur- unddichtespezifisch sind und daher bei geändernden Bedingungennachoptimiert werden müssen. Gleichzeitig läßt sichabschätzen, beiwelchen Umgebungsbedingungen dies noch nicht notwendig wird.In allenFällen reichen effektive Paarpotentiale aus, einrealistischesmesoskopisches Modell zu konstruieren. VergröberteSimulationenwerden im Falle der Polyacrylsäure erfolgreich gegenexperimentelleLichtstreudaten getestet. Wir erzielen für Molmassen bis zu300000g/mol eine hervorragende Übereinstimmung für denhydrodynamischenRadius. Unsere Ergebnisse erklären auch Korrekturen zudessenVerhalten als Funktion der Kettenlänge ('Skalenverhalten'). Im Fallevon Polyisopren untersuchen wir sowohl statische als auchdynamischeGrößen und stellen klare Unterschiede unserer Ergebnisse zudeneneines einfachen semi-flexiblen Mesoskalenmodells fest. InderProteinforschung werden aus Datenbanken gewonnene effektivePaarwechselwirkungen dazu verwendet, die freie Energie einesneuensystems vorherzusagen. Wir belegen in einem Exkurs mittelsGittersimulationen, daß es selbst in einfachsten Fällennicht gelingt,dies auch nur qualitativ korrekt zu bewerkstelligen.

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using Complete Basis Set and Gaussian-n Models combined with CPCM Continuum Solvation Methods

Complete basis set and Gaussian-n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pKa values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol, to calculate pKa values with cycle 1. The complete basis set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001

Accurate pKa Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pKa values for six carboxylic acids. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol to calculate pKa values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit.

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

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.

Novel poly(3-hydroxybutyrate) nanocomposites containing WS2 inorganic nanotubes with improved thermal, mechanical and tribological properties

Poly(3-hydroxybutyrate) (PHB) nanocomposites containing environmentally-friendly tungsten disulphide inorganic nanotubes (INTeWS2) have been successfully prepared by a simple solution blending method. The dynamic and isothermal crystallization studies by differential scanning calorimetry (DSC) demonstrated that the INTeWS2 exhibits much more prominent nucleation activity on the crystallization of PHB than specific nucleating agents or other nanoscale fillers. Both crystallization rate and crystallinity significantly increase in the nanocomposites compared to neat PHB. These changes occur without modifying the crystalline structure of PHB in the nanocomposites, as shown by wide-angle X-ray diffraction (WAXS) and infrared/Raman spectroscopy. Other parameters such as the Avrami exponent, the equilibrium melting temperature, global rate constant and the fold surface free energy of PHB chains in the nanocomposites were obtained from the calorimetric data in order to determine the influence of the INTeWS2 filler. The addition of INTeWS2 remarkably influences the energetics and kinetics of nucleation and growth of PHB, reducing the fold surface free energy by up to 20%. Furthermore, these nanocomposites also show an improvement in both tribological and mechanical (hardness and modulus) properties with respect to pure PHB evidenced by friction and nanoindentation tests, which is of important potential interest for industrial and medical applications.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Polarity and permeation profiles in lipid membranes

The isotropic 14N-hyperfine coupling constant, a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{o}^{N}}}\end{equation*}\end{document}, of nitroxide spin labels is dependent on the local environmental polarity. The dependence of a\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{o}^{N}}}\end{equation*}\end{document} in fluid phospholipid bilayer membranes on the C-atom position, n, of the nitroxide in the sn-2 chain of a spin-labeled diacyl glycerophospholipid therefore determines the transmembrane polarity profile. The polarity variation in phospholipid membranes, with and without equimolar cholesterol, is characterized by a sigmoidal, trough-like profile of the form {1 + exp [(n − no)/λ]}−1, where n = no is the point of maximum gradient, or polarity midpoint, beyond which the free energy of permeation decreases linearly with n, on a characteristic length-scale, λ. Integration over this profile yields a corresponding expression for the permeability barrier to polar solutes. For fluid membranes without cholesterol, no ≈ 8 and λ ≈ 0.5–1 CH2 units, and the permeability barrier introduces an additional diffusive resistance that is equivalent to increasing the effective membrane thickness by 35–80%, depending on the lipid. For membranes containing equimolar cholesterol, no ≈ 9–10, and the total change in polarity is greater than for membranes without cholesterol, increasing the permeability barrier by a factor of 2, whereas the decay length remains similar. The permeation of oxygen into fluid lipid membranes (determined by spin-label relaxation enhancements) displays a profile similar to that of the transmembrane polarity but of opposite sense. For fluid membranes without cholesterol no ≈ 8 and λ ≈ 1 CH2 units, also for oxygen. The permeation profile for polar paramagnetic ion complexes is closer to a single exponential decay, i.e., no lies outside the acyl-chain region of the membrane. These results are relevant not only to the permeation of water and polar solutes into membranes and their permeabilities, but also to depth determinations of site-specifically spin-labeled protein residues by using paramagnetic relaxation agents.

The equilibrium between carbon oxysulfide, carbon monoxide and sulfur ...

Thesis t.-p. pasted over the cover of the issue reprinted from the Journal of the American chemical society, vol. XXXVII, no. 9, September, 1915, which has title: The equilibrium between carbon oxysulfide, carbon monoxide and sulfur, by Gilbert N. Lewis and William N. Lacey.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Structural flexibility and oxygen diffusion pathways in monomeric fluorescent proteins

Fluorescent proteins are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. We used molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variants, and its GFP homolog citrine. We employed implicit ligand sampling and locally enhanced sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. The pathway contains several oxygen hosting pockets, which were identified by the amino acid residues that form the pocket. We calculated the free-energy of an oxygen molecule at points along the path. We also investigated an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. We showed that oxygen pathways can be blocked or altered, and barrel fluctuations can be reduced by strategic amino acid substitutions. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

VanB-type Enterococcus faecium clinical isolate successively inducibly resistant to, dependent on, and constitutively resistant to vancomycin

Three Enterococcus faecium strains isolated successively from the same patient, vancomycin-resistant strain BM4659, vancomycin-dependent strain BM4660, and vancomycin-revertant strain BM4661, were indistinguishable by pulsed-field gel electrophoresis and harbored plasmid pIP846, which confers VanB-type resistance. The vancomycin dependence of strain BM4660 was due to mutation P(175)L, which suppressed the activity of the host Ddl D-Ala:D-Ala ligase. Reversion to resistance in strain BM4661 was due to a G-to-C transversion in the transcription terminator of the vanRS(B) operon that lowered the free energy of pairing from -13.08 to -6.65 kcal/mol, leading to low-level constitutive expression of the resistance genes from the P(RB) promoter, as indicated by analysis of peptidoglycan precursors and of VanX(B) D,D-dipeptidase activity. Transcription of the resistance genes, studied by Northern hybridization and reverse transcription, initiated from the P(YB) resistance promoter, was inducible in strains BM4659 and BM4660, whereas it started from the P(RB) regulatory promoter in strain BM4661, where it was superinducible. Strain BM4661 provides the first example of reversion to vancomycin resistance of a VanB-type dependent strain not due to a compensatory mutation in the ddl or vanS(B) gene. Instead, a mutation in the transcription terminator of the regulatory genes resulted in transcriptional readthrough of the resistance genes from the P(RB) promoter in the absence of vancomycin.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

Use of binding energy by an RNA enzyme for catalysis by positioning and substrate destabilization.

A fundamental catalytic principle for protein enzymes in the use of binding interactions away from the site of chemical transformation for catalysis. We have compared the binding and reactivity of a series of oligonucleotide substrates and products of the Tetrahymena ribozyme, which catalyzes a site-specific phosphodiester cleavage reaction: CCCUCUpA+G<-->CCCUCU-OH+GpA. The results suggest that this RNA enzyme, like protein enzymes, can utilize binding interactions to achieve substantial catalysis via entropic fixation and substrate destabilization. The stronger binding of the all-ribose oligonucleotide product compared to an analog with a terminal 3' deoxyribose residue gives an effective concentration of 2200 M for the 3' hydroxyl group, a value approaching those obtained with protein enzymes and suggesting the presence of a structurally well defined active site capable of precise positioning. The stabilization from tertiary binding interactions is 40-fold less for the oligonucleotide substrate than the oligonucleotide product, despite the presence of the reactive phosphoryl group in the substrate. This destabilization is accounted for by a model in which tertiary interactions away from the site of bond cleavage position the electron-deficient 3' bridging phosphoryl oxygen of the oligonucleotide substrate next to an electropositive Mg ion. As the phosphodiester bond breaks and this 3' oxygen atom develops a negative charge in the transition state, the weak interaction of the substrate with Mg2+ becomes strong. These strategies of "substrate destabilization" and "transition state stabilization" provide estimated rate enhancements of approximately 280- and approximately 60-fold, respectively. Analogous substrate destabilization by a metal ion or hydrogen bond donor may be used more generally by RNA and protein enzymes catalyzing reactions of phosphate esters.