Modelling the super-equilibria in thermal biomass conversion : applications and limitations of the constrained free energy method

Bioenergi ses som en viktig del av det nu- och framtida sortimentet av inhemsk energi. Svartlut, bark och skogsavfall täcker mer än en femtedel av den inhemska energianvändningen. Produktionsanläggningar kan fungera ofullständigt och en mängd gas-, partikelutsläpp och tjära produceras samtidigt och kan leda till beläggningsbildning och korrosion. Orsaken till dessa problem är ofta obalans i processen: vissa föreningar anrikas i processen och superjämviktstillstånd är bildas. I denna doktorsavhandling presenteras en ny beräkningsmetod, med vilken man kan beskriva superjämviktstillståndet, de viktigaste kemiska reaktionerna, processens värmeproduktion och tillståndsstorheter samtidigt. Beräkningsmetoden grundar sig på en unik frienergimetod med bivillkor som har utvecklats vid VTT. Den här så kallade CFE-metoden har tidigare utnyttjats i pappers-, metall- och kemiindustrin. Applikationer för bioenergi, vilka är demonstrerade i doktorsavhandlingen, är ett nytt användingsområde för metoden. Studien visade att beräkningsmetoden är väl lämpad för högtemperaturenergiprocesser. Superjämviktstillstånden kan uppstå i dessa processer och det kemiska systemet kan definieras med några bivillkor. Typiska tillämpningar är förbränning av biomassa och svartlut, förgasning av biomassa och uppkomsten av kväveoxider. Också olika sätt att definiera superjämviktstillstånd presenterades i doktorsavhandlingen: empiriska konstanter, empiriska hastighetsuttryck eller reaktionsmekanismer kan användas. Resultaten av doktorsavhandlingen kan utnyttjas i framtiden i processplaneringen och i undersökning av nya tekniska lösningar för förgasning, förbränningsteknik och biobränslen. Den presenterade metoden är ett bra alternativ till de traditionella mekanistiska och fenomenmodeller och kombinerar de bästa delarna av både. --------------------------------------------------------------- Bioenergia on tärkeä osa nykyistä ja tulevaa kotimaista energiapalettia. Mustalipeä, kuori ja metsätähteet kattavat yli viidenneksen kotimaisesta energian kulutuksesta. Tuotantolaitokset eivät kuitenkaan aina toimi täydellisesti ja niiden prosesseissa syntyy erilaisia kaasu- ja hiukkaspäästöjä, tervoja sekä prosessilaitteita kuluttavia saostumia ja ruostumista. Usein syy näihin ongelmiin on prosessissa esiintyvä epätasapainotila: tietyt yhdisteet rikastuvat prosessissa ja muodostavat supertasapainotiloja. Väitöstyössä kehitettiin uusi laskentamenetelmä, jolla voidaan kuvata nämä supertasapainotilat, tärkeimmät niihin liittyvät kemialliset reaktiot, prosessin lämmöntuotanto ja tilansuureet yhtä aikaa. Laskentamenetelmä perustuu VTT:llä kehitettyyn ainutlaatuiseen rajoitettuun vapaaenergiamenetelmään. Tätä niin kutsuttua CFE-menetelmää on aiemmin sovelluttu onnistuneesti muun muassa paperi-, metalli- ja kemianteollisuudessa. Väitöstyössä esitetyt bioenergiasovellukset ovat uusi sovellusalue menetelmälle. Työ osoitti laskentatavan soveltuvan hyvin korkealämpöisiin energiatekniikan prosesseihin, joissa kemiallista systeemiä rajoittavia tekijöitä oli rajallinen määrä ja siten super-tasapainotila saattoi muodostua prosessin aikana. Tyypillisiä sovelluskohteita ovat biomassan ja mustalipeän poltto, biomassan kaasutus ja typpioksidipäästöt. Työn aikana arvioitiin myös erilaisia tapoja määritellä super-tasapainojen muodostumista rajoittavat tekijät. Rajoitukset voitiin tehdä teollisiin mittauksiin pohjautuen, kokeellisia malleja hyödyntäen tai mekanistiseen reaktiokinetiikkaan perustuen. Tulevaisuudessa väitöstyön tuloksia voidaan hyödyntää prosessisuunnittelussa ja tutkittaessa uusia teknisiä ratkaisuja kaasutus- ja polttotekniikoissa sekä biopolttoaineiden tutkimuksessa. Kehitetty menetelmä tarjoaa hyvän vaihtoehdon perinteisille mekanistisille ja ilmiömalleille yhdistäen näiden parhaita puolia.

Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Strand Analysis, a free online program for the computational identification of the best RNA interference (RNAi) targets based on Gibbs free energy

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

Relativistic deformed mean-field calculation of binding energy differences of mirror nuclei

Binding energy differences of mirror nuclei for A = 15, 17, 27, 29, 31, 33, 39 and 41 are calculated in the framework of relativistic deformed mean-field theory. To fully include the effects of the polarization of the nuclear core due to the extra particle or hole, the spatial components of the vector meson fields and the photon are taken into account in a self-consistent manner. The calculated binding energy differences are systematically smaller than the experimental values and lend support to the existency of the Okamoto-Nolen-Schiffer anomaly found decades ago in nonrelativistic calculations, For the majority of the nuclei studied, however, the results are such that the anomaly is significantly smaller than the one obtained within state-of-the-art nonrelativistic calculations.

Nonanalyticity of the free energy in thermal field theory

We study, in a d-dimensional space-time, the nonanalyticity of the thermal free energy in the scalar phi(4) theory as well as in QED. We find that the infrared divergent contributions induce, when d is even, a nonanalyticity in the coupling alpha of the form (alpha)((d-1)/2) whereas when d is odd the nonanalyticity is only logarithmic.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

Microemulsion: prediction of the phase diagram with a modified Helfrich free energy

Microemulsions are thermodynamically stable, macroscopically homogeneous but microscopically heterogeneous, mixtures of water and oil stabilised by surfactant molecules. They have unique properties like ultralow interfacial tension, large interfacial area and the ability to solubilise other immiscible liquids. Depending on the temperature and concentration, non-ionic surfactants self assemble to micelles, flat lamellar, hexagonal and sponge like bicontinuous morphologies. Microemulsions have three different macroscopic phases (a) 1phase- microemulsion (isotropic), (b) 2phase-microemulsion coexisting with either expelled water or oil and (c) 3phase- microemulsion coexisting with expelled water and oil.rnrnOne of the most important fundamental questions in this field is the relation between the properties of the surfactant monolayer at water-oil interface and those of microemulsion. This monolayer forms an extended interface whose local curvature determines the structure of the microemulsion. The main part of my thesis deals with the quantitative measurements of the temperature induced phase transitions of water-oil-nonionic microemulsions and their interpretation using the temperature dependent spontaneous curvature [c0(T)] of the surfactant monolayer. In a 1phase- region, conservation of the components determines the droplet (domain) size (R) whereas in 2phase-region, it is determined by the temperature dependence of c0(T). The Helfrich bending free energy density includes the dependence of the droplet size on c0(T) as

Arterio-venous gradients of free energy change for assessment of systemic and splanchnic perfusion in cardiac surgery patients

OBJECTIVE: Adequacy of organ perfusion depends on sufficient oxygen supply in relation to the metabolic needs. The aim of this study was to evaluate the relationship between gradients of free energy change, and the more commonly used parameter for the evaluation of the adequacy of organ perfusion, such as oxygen-extraction in patients undergoing valve replacement surgery using normothermic cardiopulmonary bypass (CPB). METHODS: In 43 cardiac patients, arterial, mixed venous, and hepato-venous blood samples were taken synchronously after induction of anaesthesia (preCPB), during CPB, and 2 and 7 h after admission to the intensive care unit (ICU+2, ICU+7). Blood gas analysis, cardiac output, and hepato-splanchnic blood flow were measured. Free energy change gradients between mixed venous and arterial (-deltadeltaG(v - a)) and hepato-venous and arterial (-deltadeltaG(hv - a)) compartments were calculated. MEASUREMENTS AND RESULTS: Cardiac index (CI) increased from 1.9 (0.7) to 2.8 (1.3) L/min/m (median, inter-quartile range) (p = 0.001), and hepato-splanchnic blood flow index (HBFI) from 0.6 (0.22) to 0.8 (0.53) L/min/m (p = 0.001). Despite increasing flow, systemic oxygen extraction increased after CPB from 24 (10)% to 35 (10)% at ICU+2 (p = 0.002), and splanchnic oxygen extraction increased during CPB from 37 (19)% to 52 (14)% (p = 0.001), and remained high thereafter. After CPB, high splanchnic and systemic gradients of free energy change gradients were associated with high splanchnic and systemic oxygen extraction, respectively (p = 0.001, 0.033, respectively). CONCLUSION: Gradients of free energy change may be helpful in characterising adequacy of perfusion in cardiac surgery patients independently from measurements or calculations of data from oxygen transport.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

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.

Magnetohydrodynamics (MHD). Hearing, Ninety-first Congress, first [and second] sessions, on pollution-free production of electrical energy from low-grade coal: magnetohydrodynamics (MHD).

Hearings held Dec. 18, 1969-