Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

The effect of counterpoise correction and relaxation energy term to the internal rotation barriers: application to the BF3 ...NH3 and C2H4 ...SO2 dimers

The relevance of the fragment relaxation energy term and the effect of the basis set superposition error on the geometry of the BF3⋯NH3 and C2H4⋯SO2 van der Waals dimers have been analyzed. Second-order Møller-Plesset perturbation theory calculations with the d95(d,p) basis set have been used to calculate the counterpoise-corrected barrier height for the internal rotations. These barriers have been obtained by relocating the stationary points on the counterpoise-corrected potential energy surface of the processes involved. The fragment relaxation energy can have a large influence on both the intermolecular parameters and barrier height. The counterpoise correction has proved to be important for these systems

Estabilitat conformacional de variants de la ribonucleasa A: pressió i temperatura com a inductors dels desplegament proteic

A fi d'analitzar la contribució de la regió C-terminal proposada com a iniciadora del plegament (CFIS 106-118) a l'estabilitat de l'RNasa A, els residus alifàtics d'aquesta regió es van substituir, mitjançant mutagènesi dirigida, per altres residus en els quals la cadena lateral alifàtica era rogressivament escurçada. La major part de les substitucions projectades suposaven delecions no disruptives de grups metil(è). A més, es va reemplaçar la Tyr115 per un Trp, de manera que, potencialment, s'introduïa una única sonda fluorescent, no desestabilitzant, per tal de seguir els canvis conformacionals que es poguessin generar en la regió durant el procés de legament/desplegament de la proteïna. Tant els paràmetres cinètics, com els espectres d'FTIR i CD, determinats per cadascuna de les ribonucleases variants, indiquen que els reemplaçaments aminoacídics efectuats presenten, en general, poc o cap efecte en l'estructura nativa i en l'activitat de l'enzim. Es va emprar l'espectroscòpia d'absorció a l'ultraviolat de quarta derivada, la fluorescència (per la variant amb Trp) i l'espectroscòpia d'infraroig per transformada de Fourier, per tal de seguir i caracteritzar, en condicions d'equilibri, les transicions conformacionals de cada variant en funció de la pressió i de la temperatura. Els resultats es van comparar amb els que es van obtenir per la proteïna salvatge. Per determinar més a fons les característiques del procés de desplegament de la variant Y115W, les transicions de desnaturalització induïdes per urea d'aquesta variant i de la proteïna salvatge, van ésser examinades per mitjà d'electroforesi en gradient d'urea i espectroscòpia de fluorescència. Curiosament, els canvis conformacionals que resulten de la desnaturalització per pressió són molt semblants als que s'obtenen per temperatura. Enfront d'un augment gradual tant de pressió com de temperatura, l'estructura terciària i els elements d'estructura secundària de les proteïnes estudiades es perden de manera conjunta i reversible. Aquestes variacions estructurals que es promouen descriuen un procés de desplegament molt cooperatiu i en dos estats. Atès que ambdues tècniques (UV i FTIR) utilitzen cadascuna un règim de concentració proteica molt diferent, els resultats indiquen que el procés de desplegament per pressió i per temperatura és intramolecular. Els resultats obtinguts suggereixen que la hidrofobicitat i el volum de les cadenes laterals del CFIS, juntament amb les interaccions de van der Waals entre elements d'estructura secundària intervenen de manera molt notable en l'estabilització de la proteïna. Entre els diferents aminoàcids alifàtics que pertanyen al CFIS C-terminal, la Val108 és el residu més important per tal de preservar la integritat estructural de l'estat natiu. Els reemplaçaments en aquesta posició causen petites alteracions conformacionals i una gran desestabilització de la proteïna (per exemple, el punt mig de la transició de desnaturalització per pressió i per temperatura de la variant V108G disminueix uns 592 MPa i 25ºC, respectivament, respecte a la proteïna salvatge). D'acord amb els resultats obtinguts, la variant Y115W ofereix una sonda útil per tal de seguir la cinètica de plegament/desplegament de l'RNasa A.

Bis(hydroxy-isoindolinone)s: synthesis, stereochemistry, polymer chemistry, and supramolecular assembly

Pseudoacid chlorides of 2,5-bis(4-fluorobenzoyl) terephthalic acid and 4,6-bis(4-fluorobenzoyl) isophthalic acid condense with primary amines to afford diastereomeric bis(hydroxyindolinone)s in good isolated yields and with diamines to give high molecular weight poly(hydroxyindolinone)s. Bis-N-pyrenemethyl bis(hydroxyindolinone)s assemble, even in dipolar solvents such as DMSO, with macrocyclic diimide-sulfones to give [3]pseudorotaxanes stabilized by electronically complementary aromatic π−π-stacking and shape-complementary van der Waals interactions.

Kinetic and crystallographic studies of glucopyranose spirohydantoin and glucopyranosylamine analogs inhibitors of glycogen phosphorylase

Glycogen phosphorylase (GP) is currently exploited as a target for inhibition of hepatic glycogenolysis under high glucose conditions. Spirohydantoin of glucopyranose and N-acetyl-beta-D-glucopyranosylamine have been identified as the most potent inhibitors of GP that bind at the catalytic site. Four spirohydantoin and three beta-D-glucopyranosylamine analogs have been designed, synthesized and tested for inhibition of GP in kinetic experiments. Depending on the functional group introduced, the K(i) values varied from 16.5 microM to 1200 microM. In order to rationalize the kinetic results, we determined the crystal structures of the analogs in complex with GP. All the inhibitors bound at the catalytic site of the enzyme, by making direct and water-mediated hydrogen bonds with the protein and by inducing minor movements of the side chains of Asp283 and Asn284, of the 280s loop that blocks access of the substrate glycogen to the catalytic site, and changes in the water structure in the vicinity of the site. The differences observed in the Ki values of the analogs can be interpreted in terms of variations in hydrogen bonding and van der Waals interactions, desolvation effects, ligand conformational entropy, and displacement of water molecules on ligand binding to the catalytic site.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background: MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results: A large dataset comprising MHC-peptide structural complexes was created by remodelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion: The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Design of a turn-linker-turn foldamer by incorporating meta-amino benzoic acid in the middle of a helix forming hexapeptide sequence: A helix breaking approach

Single crystal X-ray diffraction studies reveal that the incorporation of meta-amino benzoic acid in the middle of a helix forming hexapeptide sequence such as in peptide I Boc-Ile(1)-Aib(2)-Val(3)-m-ABA(4)-Ile(5)-Aib(6)-Leu(7)-OMe (Aib: alpha-amino isobutyric acid: m-ABA: meta-amino benzoic acid) breaks the helix propagation to produce a turn-linker-turn (T-L-T) foldamer in the solid state. In the crystalline state two conformational isomers of peptide I self-assemble in antiparallel fashion through intermolecular hydrogen bonds and aromatic pi-pi interactions to form a molecular duplex. The duplexes are further interconnected through intermolecular hydrogen bonds to form a layer of peptides. The layers are stacked one on top of the other through van der Waals interactions to form hydrophilic channels filled with solvent methanol. (C) 2009 Elsevier B.V. All rights reserved.

Formation of fibrillar structures through self-assembly of designed peptide turns

Three tripeptides Boc-Phe-Aib-Val-OMe (1), Boc-Leu-Aib-p-NA-NO2 (2) and Boc-Pro-Aib-m-NA-NO2 (3) (Aib: alpha-aminoisobutyric acid; p- and m-NA: para- and meta-nitroaniline) have been designed by incorporating aromatic rings to study the self-assembly and fibril formation. Single crystal X-ray diffraction studies show that all the peptides adopt turn-like structures that are self-assembled through intermolecular hydrogen bonds and van der Waals interactions to create layers of beta-sheets. Solvent dependent NMR titration and CD studies show that the turn structures of the peptides also exist in the solution phase. The field emission scanning electron microscopic (FE-SEM) images of the peptides in the solid state reveal fibrillar structures of flat morphology that are formed through beta-sheet mediated self-assembly of the preorganized turn building blocks.

Design of supramolecular beta-sheet forming beta-turns containing aromatic rings and non-coded alpha-aminoisobutyric acid and the formation of flat fibrillar structures through self-assembly

Single crystal X-ray diffraction studies show that the three designed tripeptides Boc-Leu-Aib-m-NA-NO2 (I), Boc-Phe-Aib-m-NA-NO2 (II) and Boc-Pro-Aib-m-ABA-OMe (III) (Aib, -aminoisobutyric acid; m-NA, m-nitroaniline; m-ABA, m-aminobenzoic acid; Boc, t-butyloxycarbonyl) containing aromatic rings in the backbones adopt -turn structures that are self-assembled through intermolecular hydrogen bonds and van der Waals interactions to create layers of -sheets. Solvent-dependent NMR titration and CD studies show that the -turn structures of the peptides also exist in the solution phase. The field emission scanning electron microscopic and transmission electron microscopic images of the peptides in the solid state reveal fibrillar structures of flat morphology that are formed through -sheet mediated self-assembly of the preorganised -turn building blocks.

Geometrical isomers of [TEAH][Co(L-Se)(2)]center dot xH(2)O: synthesis, structural, spectroscopic and computational studies

Trans-1, [HNEt3][Co-III(L-Se)(2)]center dot H2O and cis-1, [HNEt3][Co-III(L-Se)(2)]center dot 3H(2)O have been synthesized and characterized by single-crystal X-ray studies. The counter ion Et3NH+ plays a crucial role in the crystal packing leading to the formation of two distinctly different supramolecular assemblies in the two complexes. In trans-1, Co-bisphenolate units and triethylamine molecules are arranged in a linear fashion leading to a supramolecular columnar assembly along the crystallographic a-axis. In this assembly, triethylammonium ions are sandwiched between successive Co-bisphenolate units and act as gluing agents joining Co-bisphenolate units on either side through C-H center dot center dot center dot pi interactions. In sharp contrast to trans-1, Co-bisphenolate units and triethylammonium ions in cis-1 are arranged in a helical supramolecular assembly through similar C-H center dot center dot center dot pi interactions along the crystallographic b-axis. The Se center dot center dot center dot Se van der Waals interactions may be responsible for the predominant occurrence of the cis-isomer. The cyclic voltammetric studies showed quasi-reversible waves for the cobalt(III) -> cobalt(II) reductions with E-1/2 = 0.635 and 0.628 V vs. Ag/AgCl for cis-1 (at similar to 5 degrees C) and trans-1 (at similar to 25 degrees C), respectively. DFT calculations show that the trans-form is the thermodynamic product with higher stability than the cis-one, which is consistent with the variable temperature H-1 NMR studies

Development of novel brush-type chiral stationary phases based on terpenoid selectors: HPLC evaluation and theoretical investigation of enantioselective binding interactions

The terpenoid chiral selectors dehydroabietic acid, 12,14-dinitrodehydroabietic acid and friedelin have been covalently linked to silica gel yielding three chiral stationary phases CSP 1, CSP 2 and CSP 3, respectively. The enantiodiscriminating capability of each one of these phases was evaluated by HPLC with four families of chiral aromatic compounds composed of alcohols, amines, phenylalanine and tryptophan amino acid derivatives and beta-lactams. The CSP 3 phase, containing a selector with a large friedelane backbone is particularly suitable for resolving free alcohols and their derivatives bearing fluorine substituents, while CSP 2 with a dehydroabietic architecture is the only phase that efficiently discriminates 1, 1'-binaphthol atropisomers. CSP 3 also gives efficient resolution of the free amines. All three phases resolve well the racemates of N-trifluoracetyl and N-3,5-dinitrobenzoyl phenylalanine amino acid ester derivatives. Good enantioseparation of beta-lactams and N-benzoyl tryptophan amino acid derivatives was achieved on CSP 1. In order to understand the structural factors that govern the chiral molecular recognition ability of these phases, molecular dynamics simulations were carried out in the gas phase with binary diastereomeric complexes formed by the selectors of CSP 1 and CSP 2 and several amino acid derivatives. Decomposition of molecular mechanics energies shows that van der Waals interactions dominate the formation of the diastereomeric transient complexes while the electrostatic binding interactions are primarily responsible for the enantioselective binding of the (R)- and (S)-analytes. Analysis of the hydrogen bonds shows that electrostatic interactions are mainly associated with the formation of N-(HO)-O-...=C enantio selective hydrogen bonds between the amide binding sites from the selectors and the carbonyl groups of the analytes. The role of mobile phase polarity, a mixture of n-hexane and propan-2-ol in different ratios, was also evaluated through molecular dynamics simulations in explicit solvent. (c) 2006 Elsevier Ltd. All rights reserved.

The impact of charge transfer and structural disorder on the thermoelectric properties of cobalt intercalated TiS2

A family of phases, CoxTiS2 (0 ≤ x ≤ 0.75) has been prepared and characterised by powder X-ray and neutron diffraction, electrical and thermal transport property measurements, thermal analysis and SQUID magnetometry. With increasing cobalt content, the structure evolves from a disordered arrangement of cobalt ions in octahedral sites located in the van der Waals’ gap (x ≤ 0.2), through three different ordered vacancy phases, to a second disordered phase at x ≥ 0.67. Powder neutron diffraction reveals that both octahedral and tetrahedral inter-layer sites are occupied in Co0.67TiS2. Charge transfer from the cobalt guest to the TiS2 host affords a systematic tuning of the electrical and thermal transport properties. At low levels of cobalt intercalation (x < 0.1), the charge transfer increases the electrical conductivity sufficiently to offset the concomitant reduction in |S|. This, together with a reduction in the overall thermal conductivity leads to thermoelectric figures of merit that are 25 % higher than that of TiS2, ZT reaching 0.30 at 573 K for CoxTiS2 with 0.04 ≤ x ≤ 0.08. Whilst the electrical conductivity is further increased at higher cobalt contents, the reduction in |S| is more marked due to the higher charge carrier concentration. Furthermore both the charge carrier and lattice contributions to the thermal conductivity are increased in the electrically conductive ordered-vacancy phases, with the result that the thermoelectric performance is significantly degraded. These results illustrate the competition between the effects of charge transfer from guest to host and the disorder generated when cobalt cations are incorporated in the inter-layer space.

MODELLING OF HIGH-PRESSURE PHASE EQUILIBRIUM IN SYSTEMS OF INTEREST IN THE FOOD ENGINEERING FIELD USING THE PENG-ROBINSON EQUATION OF STATE WITH TWO DIFFERENT MIXING RULES

In this work, thermodynamic models for fitting the phase equilibrium of binary systems were applied, aiming to predict the high pressure phase equilibrium of multicomponent systems of interest in the food engineering field, comparing the results generated by the models with new experimental data and with those from the literature. Two mixing rules were used with the Peng-Robinson equation of state, one with the mixing rule of van der Waals and the other with the composition-dependent mixing rule of Mathias et al. The systems chosen are of fundamental importance in food industries, such as the binary systems CO(2)-limonene, CO(2)-citral and CO(2)-linalool, and the ternary systems CO(2)-Limonene-Citral and CO(2)-Limonene-Linalool, where high pressure phase equilibrium knowledge is important to extract and fractionate citrus fruit essential oils. For the CO(2)-limonene system, some experimental data were also measured in this work. The results showed the high capability of the model using the composition-dependent mixing rule to model the phase equilibrium behavior of these systems.

A model for the stabilization of atomic hydrogen centers in borate glasses

A previously proposed model describing the trapping site of the interstitial atomic hydrogen in borate glasses is analyzed. In this model the atomic hydrogen is stabilized at the centers of oxygen polygons belonging to B-O ring structures in the glass network by van der Waals forces. The previously reported atomic hydrogen isothermal decay experimental data are discussed in the light of this microscopic model. A coupled differential equation system of the observed decay kinetics was solved numerically using the Runge Kutta method. The experimental untrapping activation energy of 0.7 x 10(-19) J is in good agreement with the calculated results of dispersion interaction between the stabilized atomic hydrogen and the neighboring oxygen atoms at the vertices of hexagonal ring structures. (C) 2009 Elsevier B.V. All rights reserved.

Spectroscopic Investigation on the Formation of a Dinuclear Me(2)SO-Ru(2)-Mercaptopyridine Bridged Complex: [Ru(2)Cl(3)(mu-pyS)(mu-dmso)(dmso)(4)]center dot 2H(2)O

A dinuclear ruthenium(II) complex double-bridged by an N-aromatic ligand 2-mercaptopyridine (2-pyridinethiol or 2-pyridyl mercaptan) and a methyl sulfoxide (dmso) have been characterized by X-ray crystallography. The reported compound with formula [Ru(2)Cl(3) (mu-pyS)(mu-dmso)(dmso)(4)] center dot 2H(2)O, [C(15)H(36)Cl(3)NO(7)S(6)Ru(2)] (P2/c, a = 13.8175(2) angstrom, b = 10.5608(2) angstrom, c = 21.3544 (3) angstrom, beta = 106.090(1)degrees, V = 2,994.05(8) angstrom(3), Z = 4) represents a seven-membered ring system with both rutheniums in an octahedral geometry. All the hydrogen bonds (C-H-Cl) and the van der Waals contacts give rise to three-dimensional network in the structure and add stability to the dinuclear compound. To our knowledge, this is the first time that the formation of a dinuclear ruthenium(II) complex double-bridged by an N-aromatic ligand 2-mercaptopyridine and dmso have been reported. The study also provided valuable insight into bioinorganic chemistry as continuing efforts are being made to develop metal-based cancer chemotherapeutics. A major feature of this paper is the resolution of a double bridged ruthenium structure which contributes to a better understanding of ruthenium reactivity.