Role of PheE15 gate in ligand entry and nitric oxide detoxification function of Mycobacterium tuberculosis truncated hemoglobin N

The truncated hemoglobin N, HbN, of Mycobacterium tuberculosis is endowed with a potent nitric oxide dioxygenase (NOD) activity that allows it to relieve nitrosative stress and enhance in vivo survival of its host. Despite its small size, the protein matrix of HbN hosts a two-branched tunnel, consisting of orthogonal short and long channels, that connects the heme active site to the protein surface. A novel dual-path mechanism has been suggested to drive migration of O(2) and NO to the distal heme cavity. While oxygen migrates mainly by the short path, a ligand-induced conformational change regulates opening of the long tunnel branch for NO, via a phenylalanine (PheE15) residue that acts as a gate. Site-directed mutagenesis and molecular simulations have been used to examine the gating role played by PheE15 in modulating the NOD function of HbN. Mutants carrying replacement of PheE15 with alanine, isoleucine, tyrosine and tryptophan have similar O(2)/CO association kinetics, but display significant reduction in their NOD function. Molecular simulations substantiated that mutation at the PheE15 gate confers significant changes in the long tunnel, and therefore may affect the migration of ligands. These results support the pivotal role of PheE15 gate in modulating the diffusion of NO via the long tunnel branch in the oxygenated protein, and hence the NOD function of HbN

Molecular simulations of globins: Exploring the relationship between structure, dynamics and function

Podeu consultar el llibre complet a: http://hdl.handle.net/2445/32392The discovery in the last two decades of novel members of the globin superfamily has challenged the conventional view about the structure and function of globins. Thus, peculiar structural differences are expected to have direct influence on properties related to ligand migration, binding affinity and heme reactivity. Molecular simulations are a valuable tool to gain insigth into the molecular mechanisms that underlie those structural differences, and their relationship with the diversity of functional roles. In this work, the impact of molecular simulations in exploring the linkage between structure, dynamics and function is highlighted for three representative cases: the migration of ligands through the protein matrix of truncated hemoglobins, the modulation of binding affinity by heme distortion in protoglobin, and finally the functional implications due to the equilibrium between penta- and hexacoordination of the heme with distal histidine in neuroglobin

CO rebinding kinetics and molecular dynamics simulations highlight dynamic regulation of internal cavities in human cytoglobin

Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein

Flexibilidad conformacional de moléculas bioactivas: implicaciones en diseño de fármacos y función de globinas

[spa] El estudio de la unión de un ligando con un receptor es de particular relevancia en el desarrollo de fármacos, donde se persigue encontrar un ligando con elevada afinidad hacia un determinado receptor. El ligando se une al receptor en el denominado “sitio de unión”. El proceso de unión puede dividirse en 2 etapas: i) difusión del ligando a través del receptor hasta el sitio de unión, y ii) la selección de aquella conformación (bioactiva) que cumple con las características geométricas y fisicoquímicas que le permiten unirse al receptor. En la presente memoria se ha estudiado la difusión molecular del ligando a través de la matriz proteica hasta el sitio de unión, así como la selección de la conformación bioactiva implicada en el reconocimiento ligando-receptor. En la primera parte se ha desarrollado una metodología Multinivel que proporciona un método eficiente para la exploración conformacional de moléculas tipo fármaco, permitiendo identificar estructuras estables y su población relativa. Hace uso de una exploración conformacional a nivel RM1 y posterior refinamiento a nivel más robusto de teoría: B3LYP y MP2. Se consigue así un buen balance costo/calidad para la estimación precisa de la población relativa de los confórmeros relevantes. Para compuestos en solución se utiliza la metodología de solvente continuo MST. Esto hace necesario derivar una nueva versión parametrizada del modelo MST para el hamiltoniano RM1. El modelo parametrizado ha probado proveer estimaciones precisas de la energía de solvatación para compuestos neutros. El tratamiento de compuestos iónicos es más delicado, dado que existe una dependencia fuerte entre el factor de escalado utilizado para modular la frontera electrostática entre el solvente y el soluto y el grupo ionizable. El uso de un factor de escalado adaptado al entorno químico provee una estrategia computacional viable para obtener estimaciones precisas de la energía de solvatación para compuestos cargados. Los resultados obtenidos para la histamina neutra son alentadores considerando la complejidad conformacional y tautomérica que presenta y coinciden con la mayoría de los estudios previos. En la segunda parte, se ha estudiado el impacto de las distorsiones de la planaridad del grupo hemo en la afinidad por O(2). Se halló que éstas pueden modular la afinidad de una globina por ligandos bimoleculares. Las distorsiones resultan en una disminución de la afinidad por O(2).con excepción del modo breathing, que implica la compresión-expansión del grupo hemo en el plano. El “breathing” positivo del hemo de la protoglobina “Methanosarcina Acetivorans” (MaPgb) podría contribuir a la alta afinidad por O(2) reportada dado que no presenta efecto distal, habitualmente responsable de la alta afinidad hallada en otras globinas. Por otra parte, la accesibilidad del ligando a esta proteína se logra a través de un sistema de túneles novedoso definido entre las hélices G y B (túnel 1) y B y E (túnel 2). El túnel 2, a diferencia del 1, está siempre abierto. Phe(145)G8 se halla en 2 conformaciones: abierta y cerrada, que regula la migración de ligandos a través del 1. La dimerización y la presencia de ligando unido al hemo facilitan la apertura del 1. La presencia del ligando unido al hemo es detectada por Phe(93)E11, y transmitida a Phe(145)G8 a través del cambio de las hélices B y E. Esto proporciona una explicación a los estudios cinéticos de Cristiano Viappiani. Cabe pensar que las conformaciones de ligación lenta y rápida reflejan la apertura del túnel 1 cuando se ha fijado un primer ligando en el grupo hemo. MaPgb podría participar en un proceso bimolecular, donde la entrada de un ligando facilitara de un segundo.[eng] The study of ligand-receptor interaction is relevant to drug design. In this work we have studied: ligand diffusion through the receptor to the active site and selection of the bioactive conformation implied in ligand-receptor recognition In the first part we have developed a Multilevel methodology to efficiently explore the conformational space of drug-like molecules. It is based on a conformational search at a RM1 level and a “a posteriori” refinement with B3LYP and MP2. For molecules in solution a parametrization of the MST solvation model for RM1 was needed. Precise free energies of solvation have been obtained with the new set of parameters for neutral compounds. For ionic compounds there is a strong dependence with the scale factor used to modulate the electrostatic frontier between solvent and solute. The use of environment-adapted scale factor provides a suitable strategy. The results obtained for neutral histamine are encouraging considering its conformational and tautomeric complexity and agree with most of the previous studies. In the second part, we have studied the impact of heme distortions on its O2 affinity. Distortions reduce affinity with the exception of the breathing mode (compressionexpansion) Positive breathing found in “Methanosarcina Acetivorans” protoglobin (MaPgb) could contribute to the reported high O(2) affinity since there is no distal effect, generally responsible of high affinity in globins. A new tunnel system defined between G and B helices (tunnel 1) and B and E helices (tunnel 2) allows ligand access. Contrary to tunnel 1, tunnel 2 is always open. Phe(145)G8 is found in open and closed conformations that regulate migration through tunnel 1. Dimerization and the presence of a heme-bound ligand facilitate opening of tunnel 1. The bound ligand is detected by Phe(93)E11 and transferred to Phe(145)G8 through changes in helices B and E. This could shed light on the kinetic studies by Cristiano Viappiani. Slow and fast conformations could reflect the opening of tunnel 1 when another ligand is bound. MaPgb could be involved in a bimolecular process where the entrance of a ligand would facilitate binding of a second one

Flexibilidad conformacional de moléculas bioactivas: implicaciones en diseño de fármacos y función de globinas

El estudio de la unión de un ligando con un receptor es de particular relevancia en el desarrollo de fármacos, donde se persigue encontrar un ligando con elevada afinidad hacia un determinado receptor. El ligando se une al receptor en el denominado “sitio de unión”. El proceso de unión puede dividirse en 2 etapas: i) difusión del ligando a través del receptor hasta el sitio de unión, y ii) la selección de aquella conformación (bioactiva) que cumple con las características geométricas y fisicoquímicas que le permiten unirse al receptor. En la presente memoria se ha estudiado la difusión molecular del ligando a través de la matriz proteica hasta el sitio de unión, así como la selección de la conformación bioactiva implicada en el reconocimiento ligando-receptor. En la primera parte se ha desarrollado una metodología Multinivel que proporciona un método eficiente para la exploración conformacional de moléculas tipo fármaco, permitiendo identificar estructuras estables y su población relativa. Hace uso de una exploración conformacional a nivel RM1 y posterior refinamiento a nivel más robusto de teoría: B3LYP y MP2. Se consigue así un buen balance costo/calidad para la estimación precisa de la población relativa de los confórmeros relevantes. Para compuestos en solución se utiliza la metodología de solvente continuo MST. Esto hace necesario derivar una nueva versión parametrizada del modelo MST para el hamiltoniano RM1. El modelo parametrizado ha probado proveer estimaciones precisas de la energía de solvatación para compuestos neutros. El tratamiento de compuestos iónicos es más delicado, dado que existe una dependencia fuerte entre el factor de escalado utilizado para modular la frontera electrostática entre el solvente y el soluto y el grupo ionizable. El uso de un factor de escalado adaptado al entorno químico provee una estrategia computacional viable para obtener estimaciones precisas de la energía de solvatación para compuestos cargados. Los resultados obtenidos para la histamina neutra son alentadores considerando la complejidad conformacional y tautomérica que presenta y coinciden con la mayoría de los estudios previos. En la segunda parte, se ha estudiado el impacto de las distorsiones de la planaridad del grupo hemo en la afinidad por O(2). Se halló que éstas pueden modular la afinidad de una globina por ligandos bimoleculares. Las distorsiones resultan en una disminución de la afinidad por O(2).con excepción del modo breathing, que implica la compresión-expansión del grupo hemo en el plano. El “breathing” positivo del hemo de la protoglobina “Methanosarcina Acetivorans” (MaPgb) podría contribuir a la alta afinidad por O(2) reportada dado que no presenta efecto distal, habitualmente responsable de la alta afinidad hallada en otras globinas. Por otra parte, la accesibilidad del ligando a esta proteína se logra a través de un sistema de túneles novedoso definido entre las hélices G y B (túnel 1) y B y E (túnel 2). El túnel 2, a diferencia del 1, está siempre abierto. Phe(145)G8 se halla en 2 conformaciones: abierta y cerrada, que regula la migración de ligandos a través del 1. La dimerización y la presencia de ligando unido al hemo facilitan la apertura del 1. La presencia del ligando unido al hemo es detectada por Phe(93)E11, y transmitida a Phe(145)G8 a través del cambio de las hélices B y E. Esto proporciona una explicación a los estudios cinéticos de Cristiano Viappiani. Cabe pensar que las conformaciones de ligación lenta y rápida reflejan la apertura del túnel 1 cuando se ha fijado un primer ligando en el grupo hemo. MaPgb podría participar en un proceso bimolecular, donde la entrada de un ligando facilitara de un segundo.The study of ligand-receptor interaction is relevant to drug design. In this work we have studied: ligand diffusion through the receptor to the active site and selection of the bioactive conformation implied in ligand-receptor recognition In the first part we have developed a Multilevel methodology to efficiently explore the conformational space of drug-like molecules. It is based on a conformational search at a RM1 level and a “a posteriori” refinement with B3LYP and MP2. For molecules in solution a parametrization of the MST solvation model for RM1 was needed. Precise free energies of solvation have been obtained with the new set of parameters for neutral compounds. For ionic compounds there is a strong dependence with the scale factor used to modulate the electrostatic frontier between solvent and solute. The use of environment-adapted scale factor provides a suitable strategy. The results obtained for neutral histamine are encouraging considering its conformational and tautomeric complexity and agree with most of the previous studies. In the second part, we have studied the impact of heme distortions on its O2 affinity. Distortions reduce affinity with the exception of the breathing mode (compressionexpansion) Positive breathing found in “Methanosarcina Acetivorans” protoglobin (MaPgb) could contribute to the reported high O(2) affinity since there is no distal effect, generally responsible of high affinity in globins. A new tunnel system defined between G and B helices (tunnel 1) and B and E helices (tunnel 2) allows ligand access. Contrary to tunnel 1, tunnel 2 is always open. Phe(145)G8 is found in open and closed conformations that regulate migration through tunnel 1. Dimerization and the presence of a heme-bound ligand facilitate opening of tunnel 1. The bound ligand is detected by Phe(93)E11 and transferred to Phe(145)G8 through changes in helices B and E. This could shed light on the kinetic studies by Cristiano Viappiani. Slow and fast conformations could reflect the opening of tunnel 1 when another ligand is bound. MaPgb could be involved in a bimolecular process where the entrance of a ligand would facilitate binding of a second one

Experience in Designing and Evaluating a Teleconsultation System Supporting Shared Care of Oncological Patients

This poster presents our experience in designing, developing and deploying a Web-based Teleconsultation System based on a Patient Centred Oncological Electronic Medical Record (PEMR) specifically designed to provide clinicians a cooperative work tool supporting the oncological patient management. An evaluation phase in a clinical setting was performed when the system was deployed in the hospitals. A second evaluation phase after two year of utilization has been carrying on