Molecular diversity of bla genes in Klebsiella pneumoniae and Escherichia coli isolates

Dissertation presented to obtain a Ph.D. degree in Biology, speciality Microbiology, by Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

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

Riesz potential versus fractional Laplacian

This paper starts by introducing the Grünwald–Letnikov derivative, the Riesz potential and the problem of generalizing the Laplacian. Based on these ideas, the generalizations of the Laplacian for 1D and 2D cases are studied. It is presented as a fractional version of the Cauchy–Riemann conditions and, finally, it is discussed with the n-dimensional Laplacian.

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.

Riesz potential versus fractional Laplacian

This paper starts by introducing the Grünwald–Letnikov derivative, the Riesz potential and the problem of generalizing the Laplacian. Based on these ideas, the generalizations of the Laplacian for 1D and 2D cases are studied. It is presented as a fractional version of the Cauchy–Riemann conditions and, finally, it is discussed with the n-dimensional Laplacian.

A breakthrough on Amanita phalloides poisoning: an effective antidotal effect by polymyxin B

Amanita phalloides is responsible for more than 90 % of mushroom-related fatalities, and no effective antidote is available. a-Amanitin, the main toxin of A. phalloides, inhibits RNA polymerase II (RNAP II), causing hepatic and kidney failure. In silico studies included docking and molecular dynamics simulation coupled to molecular mechanics with generalized Born and surface area method energy decomposition on RNAP II. They were performed with a clinical drug that shares chemical similarities to a-amanitin, polymyxin B. The results show that polymyxin B potentially binds to RNAP II in the same interface of a-amanitin, preventing the toxin from binding to RNAP II. In vivo, the inhibition of the mRNA transcripts elicited by a-amanitin was efficiently reverted by polymyxin B in the kidneys. Moreover, polymyxin B significantly decreased the hepatic and renal a-amanitin-induced injury as seen by the histology and hepatic aminotransferases plasma data. In the survival assay, all animals exposed to a-amanitin died within 5 days, whereas 50 % survived up to 30 days when polymyxin B was administered 4, 8, and 12 h post-a-amanitin. Moreover, a single dose of polymyxin B administered concomitantly with a-amanitin was able to guarantee 100 % survival. Polymyxin B protects RNAP II from inactivation leading to an effective prevention of organ damage and increasing survival in a-amanitin-treated animals. The present use of clinically relevant concentrations of an already human-use-approved drug prompts the use of polymyxin B as an antidote for A. phalloides poisoning in humans.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

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

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.

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.

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.

Random-field Potts model with dipolarlike interactions: Hysteresis avalanches and microstructure

A model for the study of hysteresis and avalanches in a first-order phase transition from a single variant phase to a multivariant phase is presented. The model is based on a modification of the random-field Potts model with metastable dynamics by adding a dipolar interaction term truncated at nearest neighbors. We focus our study on hysteresis loop properties, on the three-dimensional microstructure formation, and on avalanche statistics.

Landau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate

We have studied the interaction between the low-lying transverse collective oscillations and the thermal excitations of an elongated Bose-Einstein condensate by means of perturbation theory. We consider a cylindrical trapped condensate and calculate the transverse elementary excitations at zero temperature by solving the linearized Gross-Pitaevskii equations in two dimensions (2D). We use them to calculate the matrix elements between the thermal excited states and the quasi-2D collective modes. The Landau damping of transverse collective modes is studied as a function of temperature. At low temperatures, the corresponding damping rate is in agreement with the experimental data for the decay of the transverse quadrupole mode, but it is too small to explain the observed slow decay of the transverse breathing mode. The reason for this discrepancy is discussed.