Molecular mobility, composition and structure analysis in glycerol plasticised chitosan films

This study was developed with the purpose to investigate the effect of polysaccharide/plasticiser concentration on the microstructure and molecular dynamics of polymeric film systems, using transmission electron microscope imaging (TEM) and nuclear magnetic resonance (NMR) techniques. Experiments were carried out in chitosan/glycerol films prepared with solutions of different composition. The films obtained after drying and equilibration were characterised in terms of composition, thickness and water activity. Results show that glycerol quantities used in film forming solutions were responsible for films composition; while polymer/total plasticiser ratio in the solution determined the thickness (and thus structure) of the films. These results were confirmed by TEM. NMR allowed understanding the films molecular rearrangement. Two different behaviours for the two components analysed, water and glycerol were observed: the first is predominantly moving free in the matrix, while glycerol is mainly bounded to the chitosan chain. (C) 2013 Elsevier Ltd. All rights reserved.

Influence of constraining and confinement in the molecular mobility of low molecular weight materials

Dissertation presented to obtain a Ph.D. Degree in Chemical Physics

Molecular details of INH-C-10 binding to wt KatG and Its S315T mutant

Isoniazid (INH) is still one of the two most effective antitubercular drugs and is included in all recommended multitherapeutic regimens. Because of the increasing resistance of Mycobacterium tuberculosis to INH, mainly associated with mutations in the katG gene, new INH-based compounds have been proposed to circumvent this problem. In this work, we present a detailed comparative study of the molecular determinants of the interactions between wt KatG or its S315T mutant form and either INH or INH-C10, a new acylated INH derivative. MD simulations were used to explore the conformational space of both proteins, and results indicate that the S315T mutation did not have a significant impact on the average size of the access tunnel in the vicinity of these residues. Our simulations also indicate that the steric hindrance role assigned to Asp137 is transient and that electrostatic changes can be important in understanding the enzyme activity data of mutations in KatG. Additionally, molecular docking studies were used to determine the preferred modes of binding of the two substrates. Upon mutation, the apparently less favored docking solution for reaction became the most abundant, suggesting that S315T mutation favors less optimal binding modes. Moreover, the aliphatic tail in INH-C10 seems to bring the hydrazine group closer to the heme, thus favoring the apparent most reactive binding mode, regardless of the enzyme form. The ITC data is in agreement with our interpretation of the C10 alkyl chain role and helped to rationalize the significantly lower experimental MIC value observed for INH-C10. This compound seems to be able to counterbalance most of the conformational restrictions introduced by the mutation, which are thought to be responsible for the decrease in INH activity in the mutated strain. Therefore, INH-C10 appears to be a very promising lead compound for drug development.

Structural investigation of the Bacillus subtilis morphogenic factor RodZ

A thesis to obtain a Master degree in Structural and Functional Biochemistry

Fotoionização de moléculas e átomos

Este trabalho, no âmbito da física atómica e molecular, é essencialmente um estudo experimental, por fotoionização, da estrutura electrónica de moléculas e átomos de elevada instabilidade, relevantes na indústria, atmosfera terrestre e astrofísica. Os estudos de física molecular compreendem a caracterização em fase gasosa, à temperatura ambiente e a decomposição térmica controlada até 1100 K, por fotoionização com radiação de ultravioleta de vácuo usando a fonte de He I (21.22 eV), das seguintes moléculas: azidoformato de metilo, N3COOCH3; azidoformato de etilo, N3COOCH2CH3; azidoacetato de etilo, N3CH2COOCH2CH3; azidoetanol, N3CH2CH2OH;2-azidopropionitrilo, CH3CH(N3)CN; e 3-azidopropionitrilo, N3CH2CH2CN. Os resultados de espectroscopia de fotoelectrões, para cada uma destas moléculas, são comparados com os resultados do estudo de pirólise pela técnica de isolamento em matriz de azoto molecular, à temperatura de 12 K, assistida por espectroscopia de infravermelho. Os compostos intermediários e produtos finais observados permitem propor mecanismos de decomposição térmica para cada um dos sistemas. A interpretação e atribuição das bandas espectrais relativas a cada molécula têm por base os cálculos ab initio (HF e MP2) e de funcional da densidade (BLYP e B3LYP) das propriedades moleculares dos seus confórmeros de equilíbrio. A radiação de sincrotrão de ultravioleta de vácuo é usada nos estudos, por fotoionização, do oxigénio e azoto atómicos e dos radicais OH e OD, os quais são produzidos in situ por descarga de microondas e/ou reacção rápida átomo molécula. Para cada um destes sistemas, apresentam-se medidas da secção eficaz parcial de fotoionização, σ i, e do parâmetro de assimetria, β. Nestes estudos é usada a radiação da estação de trabalho 4.2 do Sincrotrão Elettra, em Itália.

Access of extracellular cations to their binding sites in Na,K-ATPase: role of the second extracellular loop of the alpha subunit.

Na,K-ATPase, the main active transport system for monovalent cations in animal cells, is responsible for maintaining Na(+) and K(+) gradients across the plasma membrane. During its transport cycle it binds three cytoplasmic Na(+) ions and releases them on the extracellular side of the membrane, and then binds two extracellular K(+) ions and releases them into the cytoplasm. The fourth, fifth, and sixth transmembrane helices of the alpha subunit of Na,K-ATPase are known to be involved in Na(+) and K(+) binding sites, but the gating mechanisms that control the access of these ions to their binding sites are not yet fully understood. We have focused on the second extracellular loop linking transmembrane segments 3 and 4 and attempted to determine its role in gating. We replaced 13 residues of this loop in the rat alpha1 subunit, from E314 to G326, by cysteine, and then studied the function of these mutants using electrophysiological techniques. We analyzed the results using a structural model obtained by homology with SERCA, and ab initio calculations for the second extracellular loop. Four mutants were markedly modified by the sulfhydryl reagent MTSET, and we investigated them in detail. The substituted cysteines were more readily accessible to MTSET in the E1 conformation for the Y315C, W317C, and I322C mutants. Mutations or derivatization of the substituted cysteines in the second extracellular loop resulted in major increases in the apparent affinity for extracellular K(+), and this was associated with a reduction in the maximum activity. The changes produced by the E314C mutation were reversed by MTSET treatment. In the W317C and I322C mutants, MTSET also induced a moderate shift of the E1/E2 equilibrium towards the E1(Na) conformation under Na/Na exchange conditions. These findings indicate that the second extracellular loop must be functionally linked to the gating mechanism that controls the access of K(+) to its binding site.

The Molecular Basis for Antimicrobial Activity of Pore-Forming Cyclic Peptides

The mechanism of action of antimicrobial peptides is, to our knowledge, still poorly understood. To probe the biophysical characteristics that confer activity, we present here a molecular-dynamics and biophysical study of a cyclic antimicrobial peptide and its inactive linear analog. In the simulations, the cyclic peptide caused large perturbations in the bilayer and cooperatively opened a disordered toroidal pore, 1–2 nm in diameter. Electrophysiology measurements confirm discrete poration events of comparable size. We also show that lysine residues aligning parallel to each other in the cyclic but not linear peptide are crucial for function. By employing dual-color fluorescence burst analysis, we show that both peptides are able to fuse/aggregate liposomes but only the cyclic peptide is able to porate them. The results provide detailed insight on the molecular basis of activity of cyclic antimicrobial peptides

Exploring chromium (VI) dioxodihalides chemistry: is density functional theory the most suitable tool?

A comparative systematic study of the CrO2F2 compound has been performed using different conventional ab initio methodologies and density functional procedures. Two points have been analyzed: first, the accuracy of results yielded by each method under study, and second, the computational cost required to reach such results. Weighing up both aspects, density functional theory has been found to be more appropriate than the Hartree-Fock (HF) and the analyzed post-HF methods. Hence, the structural characterization and spectroscopic elucidation of the full CrO2X2 series (X=F,Cl,Br,I) has been done at this level of theory. Emphasis has been given to the unknown CrO2I2 species, and specially to the UV/visible spectra of all four compounds. Furthermore, a topological analysis in terms of charge density distributions has revealed why the valence shell electron pair repulsion model fails in predicting the molecular shape of such CrO2X2 complexes

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Molecular mechanisms involved in the activation and regulation of the alpha 1-adrenergic receptor subtypes.

The adrenergic receptors (ARs) belong to the superfamily of membrane-bound G protein coupled receptors (GPCRs). Our investigation has focused on the structure-function relationship of the alpha 1b-AR subtype used as the model system for other GPCRs. Site-directed mutagenesis studies have elucidated the structural domains of the alpha 1b-AR involved in ligand binding, G protein coupling or desensitization. In addition, a combined approach using site-directed mutagenesis and molecular dynamics analysis of the alpha 1b-AR has provided information about the potential mechanisms underlying the activation process of the receptor, i.e. its transition from the 'inactive' to the 'active' conformation.

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Molecular determinants of desensitization in an ENaC/degenerin channel.

Epithelial Na(+) channel (ENaC)/degenerin family members are involved in mechanosensation, blood pressure control, pain sensation, and the expression of fear. Several of these channel types display a form of desensitization that allows the channel to limit Na(+) influx during prolonged stimulation. We used site-directed mutagenesis and chemical modification, functional analysis, and molecular dynamics simulations to investigate the role of the lower palm domain of the acid-sensing ion channel 1, a member of the ENaC/degenerin family. The lower palm domains of this trimeric channel are arranged around a central vestibule, at ∼20 Å above the plasma membrane and are covalently linked to the transmembrane channel parts. We show that the lower palm domains approach one another during desensitization. Residues in the palm co-determine the pH dependence of desensitization, its kinetics, and the stability of the desensitized state. Mutations of palm residues impair desensitization by preventing the closing movement of the palm. Overexpression of desensitization-impaired channel mutants in central neurons allowed--in contrast to overexpression of wild type--a sustained signaling response to rapid pH fluctuations. We identify and describe here the function of an important regulatory domain that most likely has a conserved role in ENaC/degenerin channels.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Static structure and dynamics of the liquid Li-Na and Li-Mg alloys

We present calculations for the static structure and ordering properties of two lithium-based s-p bonded liquid alloys, Li-Na and Li-Mg. Our theoretical approach is based on the neutral pseudoatom method to derive the interatomic pair potentials, and on the modified-hypernetted-chain theory of liquids to obtain the liquid static structure, leading to a whole combination that is free of adjustable parameters. The study is complemented by performing molecular dynamics simulations which, besides checking the theoretical static structural results, also allow a calculation of some dynamical properties. The obtained results are compared with the available experimental data.