The thermodynamics and statistical mechanics of protein folding

The physics of self-organization and complexity is manifested on a variety of biological scales, from large ecosystems to the molecular level. Protein molecules exhibit characteristics of complex systems in terms of their structure, dynamics, and function. Proteins have the extraordinary ability to fold to a specific functional three-dimensional shape, starting from a random coil, in a biologically relevant time. How they accomplish this is one of the secrets of life. In this work, theoretical research into understanding this remarkable behavior is discussed. Thermodynamic and statistical mechanical tools are used in order to investigate the protein folding dynamics and stability. Theoretical analyses of the results from computer simulation of the dynamics of a four-helix bundle show that the excluded volume entropic effects are very important in protein dynamics and crucial for protein stability. The dramatic effects of changing the size of sidechains imply that a strategic placement of amino acid residues with a particular size may be an important consideration in protein engineering. Another investigation deals with modeling protein structural transitions as a phase transition. Using finite size scaling theory, the nature of unfolding transition of a four-helix bundle protein was investigated and critical exponents for the transition were calculated for various hydrophobic strengths in the core. It is found that the order of the transition changes from first to higher order as the strength of the hydrophobic interaction in the core region is significantly increased. Finally, a detailed kinetic and thermodynamic analysis was carried out in a model two-helix bundle. The connection between the structural free-energy landscape and folding kinetics was quantified. I show how simple protein engineering, by changing the hydropathy of a small number of amino acids, can enhance protein folding by significantly changing the free energy landscape so that kinetic traps are removed. The results have general applicability in protein engineering as well as understanding the underlying physical mechanisms of protein folding. ^

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

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.

Estudo teórico de intermediários tetraédricos acidez / basicidade e estereosseletividade enzimáticos

The present work aimed first, the theoretical study of tetrahedral intermediate stability formed from carbonyl addition reactions using the second (MP2) and third (MP3) order Møller–Plesset perturbation theory. Linear correlations between electronic energy difference of reactions with Wiberg Indexes and C-O bond lengths were obtained, and was observed that the stability of adducts formed depends directly of electronic density involved between these atoms. The knowing of electronic parameters of these structures has an important hole due to the large use on reactions that in his course forms this tetrahedral intermediate. Employing the ONIOM (B3LYP:AMBER) methodology, was evaluated the stereoselectivity of a enzymatic reaction between CAL B enzyme and a long chain ester. In this study, were obtained the electronic energies of ground state and intermediate state of transesterification rate-determing step from two possible proquirals faces Re and Si. The objective was study the enantioselectivity of CAL B and rationalizes it using quantum theory of atoms in molecules (QTAIM). A theoretical study employing inorganic compounds was performed using ab initio CBS-QB3 method aiming to find a link between thermodynamic and equilibrium involving acids and bases. The results observed showed an excellent relationship between difference in Gibbs free energy, ΔG of acid dissociation reaction and ΔG of hydrolysis reaction of the corresponding conjugate base. It was also observed, a relationship between ΔG of hydrolysis reaction of conjugate acids and their corresponding atomic radius showing that stability plays an important role in hydrolysis reactions. The importance of solvation in acid/base behavior when compared to theoretical and experimental ΔG´s also was evaluated.

Síntese e caracterização de materiais mesoporosos para a captura de CO2: influência do óxido de níquel

Several materials are currently under study for the CO2 capture process, like the metal oxides and mixed metal oxides, zeolites, carbonaceous materials, metal-organic frameworks (MOF's) organosilica and modified silica surfaces. In this work, evaluated the adsorption capacity of CO2 in mesoporous materials of different structures, such as MCM-48 and SBA- 15 without impregnating and impregnated with nickel in the proportions 5 %, 10 % and 20 % (m/m), known as 5Ni-MCM-48, 10Ni-MCM-48, 20Ni-MCM-48 and 5Ni-SBA-15, 10NiSBA-15, 20Ni-SBA-15. The materials were characterized by means of X-ray diffraction (XRD), thermal analysis (TG and DTG), Fourier transform infrared spectroscopy (FT-IR), N2 adsorption and desorption (BET) and scanning electron microscopy (SEM) with EDS. The adsorption process was performed varying the pressure of 100 - 4000 kPa and keeping the temperature constant and equal to 298 K. At a pressure of 100 kPa, higher concentrations of adsorption occurred for the materials 5Ni-MCM-48 (0.795 mmol g-1 ) and SBA-15 (0.914 mmol g-1 ) is not impregnated, and at a pressure of 4000 kPa for MCM-48 materials (14.89 mmol g-1) and SBA-15 (9.97 mmol g-1) not impregnated. The results showed that the adsorption capacity varies positively with the specific area, however, has a direct dependency on the type and geometry of the porous structure of channels. The data were fitted using the Langmuir and Freundlich models and were evaluated thermodynamic parameters Gibbs free energy and entropy of the adsorption system

Contribuição ao estudo de treliças não lineares considerando efeitos térmicos

This work presents the numerical analysis of nonlinear trusses summited to thermomechanical actions with Finite Element Method (FEM). The proposed formulation is so-called positional FEM and it is based on the minimum potential energy theorem written according to nodal positions, instead of displacements. The study herein presented considers the effects of geometric and material nonlinearities. Related to dynamic problems, a comparison between different time integration algorithms is performed. The formulation is extended to impact problems between trusses and rigid wall, where the nodal positions are constrained considering nullpenetration condition. In addition, it is presented a thermodynamically consistent formulation, based on the first and second law of thermodynamics and the Helmholtz free-energy for analyzing dynamic problems of truss structures with thermoelastic and thermoplastic behavior. The numerical results of the proposed formulation are compared with examples found in the literature.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Phase coexistence of cluster crystals: beyond the Gibbs phase rule

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants.

Numerical detection of the Gardner transition in a mean-field glass former.

Recent theoretical advances predict the existence, deep into the glass phase, of a novel phase transition, the so-called Gardner transition. This transition is associated with the emergence of a complex free energy landscape composed of many marginally stable sub-basins within a glass metabasin. In this study, we explore several methods to detect numerically the Gardner transition in a simple structural glass former, the infinite-range Mari-Kurchan model. The transition point is robustly located from three independent approaches: (i) the divergence of the characteristic relaxation time, (ii) the divergence of the caging susceptibility, and (iii) the abnormal tail in the probability distribution function of cage order parameters. We show that the numerical results are fully consistent with the theoretical expectation. The methods we propose may also be generalized to more realistic numerical models as well as to experimental systems.

Phonon anharmonicity and negative thermal expansion in SnSe

The anharmonic phonon properties of SnSe in the Pnma phase were investigated with a combination of experiments and first-principles simulations. Using inelastic neutron scattering (INS) and nuclear resonant inelastic X-ray scattering (NRIXS), we have measured the phonon dispersions and density of states (DOS) and their temperature dependence, which revealed a strong, inhomogeneous shift and broadening of the spectrum on warming. First-principles simulations were performed to rationalize these measurements, and to explain the previously reported anisotropic thermal expansion, in particular the negative thermal expansion within the Sn-Se bilayers. Including the anisotropic strain dependence of the phonon free energy, in addition to the electronic ground state energy, is essential to reproduce the negative thermal expansion. From the phonon DOS obtained with INS and additional calorimetry measurements, we quantify the harmonic, dilational, and anharmonic components of the phonon entropy, heat capacity, and free energy. The origin of the anharmonic phonon thermodynamics is linked to the electronic structure.

Étude des mécanismes d’extraction lipidique par le peptide mélittine et la protéine BSP1

Les peptides et protéines extracteurs de lipides (PEL) se lient aux membranes lipidiques puis en extraient des lipides en formant de plus petits auto-assemblages, un phénomène qui peut aller jusqu'à la fragmentation des membranes. Dans la nature, cette extraction se produit sur une gamme de cellules et entraîne des conséquences variées, comme la modification de la composition de la membrane et la mort de la cellule. Cette thèse se penche sur l’extraction lipidique, ou fragmentation, induite par le peptide mélittine et la protéine Binder-of-SPerm 1 (BSP1) sur des membranes lipidiques modèles. Pour ce faire, des liposomes de différentes compositions sont préparés et incubés avec la mélittine ou la BSP1. L'association aux membranes est déterminée par la fluorescence intrinsèque des PEL, tandis que l'extraction est caractérisée par une plateforme analytique combinant des tests colorimétriques et des analyses en chromatographie en phase liquide et spectrométrie de masse (LCMS). La mélittine fait partie des peptides antimicrobiens cationiques, un groupe de PEL très répandu chez les organismes vivants. Ces peptides sont intéressants du point du vue médical étant donné leur mode d’action qui vise directement les lipides des membranes. Plusieurs de ceux-ci agissent sur les membranes des bactéries selon le mécanisme dit « en tapis », par lequel ils s’adsorbent à leur surface, forment des pores et ultimement causent leur fragmentation. Dans cette thèse, la mélittine est utilisée comme peptide modèle afin d’étudier le mécanisme par lequel les peptides antimicrobiens cationiques fragmentent les membranes. Les résultats montrent que la fragmentation des membranes de phosphatidylcholines (PC) est réduite par une déméthylation graduelle de leur groupement ammonium. L'analyse du matériel fragmenté révèle que les PC sont préférentiellement extraites des membranes, dû à un enrichissement local en PC autour de la mélittine à l'intérieur de la membrane. De plus, un analogue de la mélittine, dont la majorité des résidus cationiques sont neutralisés, est utilisé pour évaluer le rôle du caractère cationique de la mélittine native. La neutralisation augmente l'affinité du peptide pour les membranes neutres et anioniques, réduit la fragmentation des membranes neutres et augmente la fragmentation des membranes anioniques. Malgré les interactions électrostatiques entre le peptide cationique et les lipides anioniques, aucune spécificité lipidique n'est observée dans l'extraction. La BSP1 est la protéine la plus abondante du liquide séminal bovin et constitue un autre exemple de PEL naturel important. Elle se mélange aux spermatozoïdes lors de l’éjaculation et extrait des lipides de leur membrane, notamment le cholestérol et les phosphatidylcholines. Cette étape cruciale modifie la composition lipidique de la membrane du spermatozoïde, ce qui faciliterait par la suite la fécondation de l’ovule. Cependant, le contact prolongé de la protéine avec les spermatozoïdes endommagerait la semence. Cette thèse cherche donc à approfondir notre compréhension de ce délicat phénomène en étudiant le mécanisme moléculaire par lequel la protéine fragmente les membranes lipidiques. Les résultats des présents travaux permettent de proposer un mécanisme d’extraction lipidique en 3 étapes : 1) L'association à l’interface des membranes; 2) La relocalisation de l’interface vers le cœur lipidique; 3) La fragmentation des membranes. La BSP1 se lie directement à deux PC à l'interface; une quantité suffisante de PC dans les membranes est nécessaire pour permettre l'association et la fragmentation. Cette liaison spécifique ne mène généralement pas à une extraction lipidique sélective. L'impact des insaturations des chaînes lipidiques, de la présence de lysophosphatidylcholines, de phosphatidyléthanolamine, de cholestérol et de lipides anioniques est également évalué. Les présentes observations soulignent la complexe relation entre l'affinité d'un PEL pour une membrane et le niveau de fragmentation qu'il induit. L'importance de la relocalisation des PEL de l'interface vers le cœur hydrophobe des membranes pour permettre leur fragmentation est réitérée. Cette fragmentation semble s'accompagner d'une extraction lipidique préférentielle seulement lorsqu'une séparation de phase est induite au niveau de la membrane, nonobstant les interactions spécifiques PEL-lipide. Les prévalences des structures amphiphiles chez certains PEL, ainsi que de la fragmentation en auto-assemblages discoïdaux sont discutées. Finalement, le rôle des interactions électrostatiques entre les peptides antimicrobiens cationiques et les membranes bactériennes anioniques est nuancé : les résidus chargés diminueraient l'association des peptides aux membranes neutres suite à l'augmentation de leur énergie de solvatation.

Tracing the Chromospheric and Coronal Magnetic Field with AIA, IRIS, IBIS, and ROSA Data

The aim of this study is to explore the suitability of chromospheric images for magnetic modeling of active regions. We use high-resolutionimages (≈0.2"-0.3"), from the Interferometric Bidimensional Spectrometer in the Ca II 8542 Å line, the Rapid Oscillations in the Solar Atmosphere instrument in the Hα 6563Å line, the Interface Region Imaging Spectrograph in the 2796Å line, and compare non-potential magnetic field models obtainedfrom those chromospheric images with those obtained from images of the Atmospheric Imaging Assembly in coronal (171 Å, etc.) and inchromospheric (304 Å) wavelengths. Curvi-linear structures are automatically traced in those images with the OCCULT-2 code, to which we forward-fitted magnetic field lines computed with the Vertical-current Approximation Nonlinear Force Free Field code. We find that the chromospheric images: (1) reveal crisp curvi-linear structures (fibrils, loop segments, spicules) that are extremely well-suited for constraining magnetic modeling; (2) that these curvi-linear structures arefield-aligned with the best-fit solution by a median misalignment angle of μ2 ≈ 4°–7° (3) the free energy computed from coronal data may underestimate that obtained from chromospheric data by a factor of ≈2–4, (4) the height range of chromospheric features is confined to h≲4000 km, while coronal features are detected up to h = 35,000 km; and (5) the plasma-β parameter is β ≈ 10^-5 - 10^-1 for all traced features. We conclude that chromospheric images reveal important magnetic structures that are complementary to coronal images and need to be included in comprehensive magnetic field models, something that is currently not accomodated in standard NLFFF codes.

Novel ketone diet enhances physical and cognitive performance

Ketone bodies are the most energy-efficient fuel and yield more ATP per mole of substrate than pyruvate and increase the free energy released from ATP hydrolysis. Elevation of circulating ketones via high-fat, low-carbohydrate diets has been used for the treatment of drug-refractory epilepsy and for neurodegenerative diseases, such as Parkinson's disease. Ketones may also be beneficial for muscle and brain in times of stress, such as endurance exercise. The challenge has been to raise circulating ketone levels by using a palatable diet without altering lipid levels. We found that blood ketone levels can be increased and cholesterol and triglycerides decreased by feeding rats a novel ketone ester diet: chow that is supplemented with (R)-3-hydroxybutyl (R)-3-hydroxybutyrate as 30% of calories. For 5 d, rats on the ketone diet ran 32% further on a treadmill than did control rats that ate an isocaloric diet that was supplemented with either corn starch or palm oil (P < 0.05). Ketone-fed rats completed an 8-arm radial maze test 38% faster than did those on the other diets, making more correct decisions before making a mistake (P < 0.05). Isolated, perfused hearts from rats that were fed the ketone diet had greater free energy available from ATP hydrolysis during increased work than did hearts from rats on the other diets as shown by using [(31)P]-NMR spectroscopy. The novel ketone diet, therefore, improved physical performance and cognitive function in rats, and its energy-sparing properties suggest that it may help to treat a range of human conditions with metabolic abnormalities.

X-Ray Diffraction of Highway Materials, HR-128, 1969

X-ray diffraction data during adsorption of water vapor on Na- and Ca-montmorillonites show that interlayer expansion is continuous but nonuniform. X-ray and adsorption isotherm data indicate an ice-like configuration of water molecules is completed with the fourth layer of interlayer water for the Ca-clay; a fifth layer intrudes to give a less ordered structure. Data for the Na-clay indicate a laminar stacking arrangement for up to three layers of interlayer water. The Na-clay adsorbs more than twice as much water and undergoes four times as large a volume change than the Ca-clay. The free energy change during adsorption of water vapor on the Ca-clay is nearly twice that for the Na-clay. Free energy changes with increasing relative pressure reflect interlayer expansion increments.