Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.

GHMP Kinases - Structures, Mechanisms and Potential for Therapeutically Relevant Inhibition

The GHMP kinases are a structurally related family of small molecule kinases named after four of its members - galactokinase, homoserine kinase, mevalonate kinase and phosphomevalonate kinase. The group also includes the enzymes N-acetylgalactosamine kinase, arabinose kinase, mevalonate 5-diphosphate decarboxylase, archeal shikimate kinase and 4-(cytidine 5'-diphospho)-2-c-methyl-D-erythritol kinase. In addition the group includes two members not known to be catalytically active, the Caenorhabditis elegans sex-fate determining protein XOL-1 and the Saccharomyces cerevisiae transcriptional activator Gal3p. Two catalytic mechanisms have been proposed for GHMP kinases. The structure of mevalonate kinase suggests that an aspartate residue acts as an active site base, removing a proton from the substrate to facilitate attack on the ? phosphate of MgATP. In contrast, in homoserine kinase there is no potential catalytic base and it is proposed that catalysis is driven by transition state stabilisation. Potential chemotherapeutic interventions against GHMP kinases fall into three main categories: inhibition of galactokinase to assist suffers of galactosemia, inhibition of mevalonate kinase or mevalonate 5-diphosphate decarboxylase to reduce flux through the cholesterol biosynthesis pathway and inhibition of bacterial GHMP kinases for novel anti-microbial therapies. These are in the early stages of development, but the accumulation of structural and mechanistic data will assist future progress.

The structure of an unconventional HD-GYP protein from Bdellovibrio reveals the roles of conserved residues in this class of cyclic-di-GMP phosphodiesterases

UNLABELLED: Cyclic-di-GMP is a near-ubiquitous bacterial second messenger that is important in localized signal transmission during the control of various processes, including virulence and switching between planktonic and biofilm-based lifestyles. Cyclic-di-GMP is synthesized by GGDEF diguanylate cyclases and hydrolyzed by EAL or HD-GYP phosphodiesterases, with each functional domain often appended to distinct sensory modules. HD-GYP domain proteins have resisted structural analysis, but here we present the first structural representative of this family (1.28 Å), obtained using the unusual Bd1817 HD-GYP protein from the predatory bacterium Bdellovibrio bacteriovorus. Bd1817 lacks the active-site tyrosine present in most HD-GYP family members yet remains an excellent model of their features, sharing 48% sequence similarity with the archetype RpfG. The protein structure is highly modular and thus provides a basis for delineating domain boundaries in other stimulus-dependent homologues. Conserved residues in the HD-GYP family cluster around a binuclear metal center, which is observed complexed to a molecule of phosphate, providing information on the mode of hydroxide ion attack on substrate. The fold and active site of the HD-GYP domain are different from those of EAL proteins, and restricted access to the active-site cleft is indicative of a different mode of activity regulation. The region encompassing the GYP motif has a novel conformation and is surface exposed and available for complexation with binding partners, including GGDEF proteins.

IMPORTANCE: It is becoming apparent that many bacteria use the signaling molecule cyclic-di-GMP to regulate a variety of processes, most notably, transitions between motility and sessility. Importantly, this regulation is central to several traits implicated in chronic disease (adhesion, biofilm formation, and virulence gene expression). The mechanisms of cyclic-di-GMP synthesis via GGDEF enzymes and hydrolysis via EAL enzymes have been suggested by the analysis of several crystal structures, but no information has been available to date for the unrelated HD-GYP class of hydrolases. Here we present the multidomain structure of an unusual member of the HD-GYP family from the predatory bacterium Bdellovibrio bacteriovorus and detail the features that distinguish it from the wider structural family of general HD fold hydrolases. The structure reveals how a binuclear iron center is formed from several conserved residues and provides a basis for understanding HD-GYP family sequence requirements for c-di-GMP hydrolysis.

A predatory patchwork:membrane and surface structures of Bdellovibrio bacteriovorus

Predatory Bdellovibrio bacteriovorus bacteria are remarkable in that they attach to, penetrate and digest other Gram-negative bacteria, living and replicating within them until all resources are exhausted, when they escape the prey ghost to invade fresh prey. Remarkable remodeling of both predator and prey cell occurs during this process to allow the Bdellovibrio to exploit the intracellular niche they have worked so hard to enter, keeping the prey "bdelloplast" intact until the end of predatory growth. If one views motile non-predatory bacteria in a light microscope, one is immediately struck by how rare it is for bacteria to collide. This highlights how the cell surface of Bdellovibrio must be specialized and adapted to allow productive collisions and further to allow entry into the prey periplasm and subsequent secretion of hydrolytic enzymes to digest it. Bdellovibrio can, however, also be made to grow artificially without prey; thus, they have a large genome containing both predatory genes and genes for saprophytic heterotrophic growth. Thus, the membrane and outer surface layers are a patchwork of proteins encompassing not only those that have a sole purpose in heterotrophic growth but also many more that are specialized or employed to attach to, enter, remodel, kill and ultimately digest prey cells. There is much that is as yet not understood, but molecular genetic and post-genomic approaches to microbial physiology have enhanced the pioneering biochemical work of four decades ago in characterizing some of the key events and surface protein requirements for prey attack.

General, PDB-based collective variables for protein folding

New, automated forms of data-analysis are required in order to understand the high-dimensional trajectories that are obtained from molecular dynamics simulations on proteins. Dimensionality reduction algorithms are particularly appealing in this regard as they allow one to construct unbiased, low-dimensional representations of the trajectory using only the information encoded in the trajectory. The downside of this approach is that different sets of coordinates are required for each different chemical systems under study precisely because the coordinates are constructed using information from the trajectory. In this paper we show how one can resolve this problem by using the sketch-map algorithm that we recently proposed to construct a low-dimensional representation of the structures contained in the protein data bank (PDB). We show that the resulting coordinates are as useful for analysing trajectory data as coordinates constructed using landmark configurations taken from the trajectory and that these coordinates can thus be used for understanding protein folding across a range of systems.

Synthesis, characterization, electrochemical behavior and in vitro protein tyrosine kinase inhibitory activity of the cymene-halogenobenzohydroxamato [Ru(eta(6)-cymene)(bha)Cl] complexes

The ruthenium(II)-cymene complexes [Ru(eta(6)-cymene)(bha)Cl] with substituted halogenobenzohydroxamato (bha) ligands (substituents = 4-F, 4-Cl, 4-Br, 2,4-F-2, 3,4-F-2, 2,5-F-2, 2,6-F-2) have been synthesized and characterized by elemental analysis, IR, H-1 NMR, C-13 NMR, cyclic voltammetry and controlled-potential electrolysis, and density functional theory (DFT) studies. The compositions of their frontier molecular orbitals (MOs) were established by DFT calculations, and the oxidation and reduction potentials are shown to follow the orders of the estimated vertical ionization potential and electron affinity, respectively. The electrochemical E-L Lever parameter is estimated for the first time for the various bha ligands, which can thus be ordered according to their electron-donor character. All complexes exhibit very strong protein tyrosine kinase (PTK) inhibitory activity, even much higher than that of genistein, the clinically used PTK inhibitory drug. The complex containing the 2,4-difluorobenzohydroxamato ligand is the most active one, and the dependences of the PTK activity of the complexes and of their redox potentials on the ring substituents are discussed. (C) 2012 Elsevier B.V. All rights reserved.

A minimal version of a Protein Tyrosine Phosphatase

Phosphatase and tensin homologue (PTEN) protein belongs to the family of protein tyrosine phos-phatase. Mutations on the phosphatase and tensin homologue (PTEN) protein are highly observed in diverse types of human tumors, being mostly identified on the phosphatase domain of the protein. Although PTEN is a modular protein composed by a phosphatase domain and a C2 domain for mem-brane anchoring, this work aimed at developing a minimal version of PTEN´s phosphatase domain. The minimal version (Small Domain) comprises a 28 residue peptide, with the PTEN 8-mer catalytic peptide accommodated between a α-helix and β-turn as observed in PTEN native structure. Firstly, a de novo prediction of the Small Domain´s secondary structure was carried out by molecular modeling tools. The stability of the predicted structures were then evaluated by Molecular Dynamics. Automated molecular docking of PTEN natural substrate PIP3, its analogue (Inositol) and a PTEN inhibitor (L-tar-tare) were performed with the modeled structure, and PTEN used as a positive control. The gene en-coding for Small Domain was designed and cloned into an expression vector at N-terminal of Green Fluorescence Protein (GFP) encoding gene. The fusion protein was then expressed in Escherichia coli cells. Different expression conditions have been explored for the production of the fusion protein to minimize the formation of inclusion bodies.

Peroxisome proliferator-activated receptor structures: ligand specificity, molecular switch and interactions with regulators.

Peroxisome proliferator-activated receptors (PPARs) compose a family of nuclear receptors that mediate the effects of lipidic ligands at the transcriptional level. In this review, we highlight advances in the understanding of the PPAR ligand binding domain (LBD) structure at the atomic level. The overall structure of PPARs LBD is described, and important protein ligand interactions are presented. Structure-activity relationships between isotypes structures and ligand specificity are addressed. It is shown that the numerous experimental three-dimensional structures available, together with in silico simulations, help understanding the role played by the activating function-2 (AF-2) in PPARs activation and its underlying molecular mechanism. The relation between the PPARs constitutive activity and the intrinsic stability of the active conformation is discussed. Finally, the interactions of PPARs LBD with co-activators or co-repressors, as well as with the retinoid X receptor (RXR) are described and considered in relation to PPARs activation.

Profiling of T-cell receptor signaling complex assembly in human CD4 T-lymphocytes using RP protein arrays.

The RP protein (RPP) array approach immobilizes minute amounts of cell lysates or tissue protein extracts as distinct microspots on NC-coated slide. Subsequent detection with specific antibodies allows multiplexed quantification of proteins and their modifications at a scale that is beyond what traditional techniques can achieve. Cellular functions are the result of the coordinated action of signaling proteins assembled in macromolecular complexes. These signaling complexes are highly dynamic structures that change their composition with time and space to adapt to cell environment. Their comprehensive analysis requires until now relatively large amounts of cells (>5 x 10(7)) due to their low abundance and breakdown during isolation procedure. In this study, we combined small scale affinity capture of the T-cell receptor (TCR) and RPP arrays to follow TCR signaling complex assembly in human ex vivo isolated CD4 T-cells. Using this strategy, we report specific recruitment of signaling components to the TCR complex upon T-cell activation in as few as 0.5 million of cells. Second- to fourth-order TCR interacting proteins were accurately quantified, making this strategy specially well-suited to the analysis of membrane-associated signaling complexes in limited amounts of cells or tissues, e.g., ex vivo isolated cells or clinical specimens.

Visualisation of human rad52 protein and its complexes with hRad51 and DNA.

The human Rad52 protein stimulates joint molecule formation by hRad51, a homologue of Escherichia coli RecA protein. Electron microscopic analysis of hRad52 shows that it self-associates to form ring structures with a diameter of approximately 10 nm. Each ring contains a hole at its centre. hRad52 binds to single and double-stranded DNA. In the ssDNA-hRad52 complexes, hRad52 was distributed along the length of the DNA, which exhibited a characteristic "beads on a string" appearance. At higher concentrations of hRad52, "super-rings" (approximately 30 nm) were observed and the ssDNA was collapsed upon itself. In contrast, in dsDNA-hRad52 complexes, some regions of the DNA remained protein-free while others, containing hRad52, interacted to form large protein-DNA networks. Saturating concentrations of hRad51 displaced hRad52 from ssDNA, whereas dsDNA-Rad52 complexes (networks) were more resistant to hRad51 invasion and nucleoprotein filament formation. When Rad52-Rad51-DNA complexes were probed with gold-conjugated hRad52 antibodies, the presence of globular hRad52 structures within the Rad51 nucleoprotein filament was observed. These data provide the first direct visualisation of protein-DNA complexes formed by the human Rad51 and Rad52 recombination/repair proteins.

Myocardial angiogenesis induction with bone protein derived growth factors (animal experiment).

Myocardial angiogenesis induction with vascular growth factors constitutes a potential strategy for patients whose coronary artery disease is refractory to conventional treatment. The importance of angiogenesis in bone formation has led to the development of growth factors derived from bovine bone protein. Twelve pigs (mean weight, 73 +/- 3 kg) were chosen for the study. In the first group (n = 6, growth factor group) five 100 micrograms boluses of growth factors derived from bovine bone protein, diluted in Povidone 5%, were injected in the lateral wall of the left ventricle. In the second group (n = 6, control group), the same operation was performed but only the diluting agent was injected. All the animals were sacrificed after 28 days and the vascular density of the left lateral wall (expressed as the number of vascular structures per mm2) as well as the area of blood vessel profiles per myocardial area analysed were determined histologically with a computerised system. The growth factor group had a capillary density which was significantly higher than that of the control group: 12.6 +/- 0.9/mm2 vs 4.8 +/- 0.5/mm2 (p < 0.01). The same holds true for the arteriolar density: 1 +/- 0.2/mm2 vs 0.3 +/- 0.1/mm2 (p < 0.01). The surface ratios of blood vessel profiles per myocardial area were 4900 +/- 800 micron 2/mm2 and 1550 +/- 400 micron 2/mm2 (p < 0.01) respectively. In this experimental model, bovine bone protein derived growth factors induce a significant neovascularisation in healthy myocardium, and appear therefore as promising candidates for therapeutic angiogenesis.

Modeling Metal Protein Complexes from Experimental Extended X-ray Absorption Fine Structure using Computational Intelligence

Experimental Extended X-ray Absorption Fine Structure (EXAFS) spectra carry information about the chemical structure of metal protein complexes. However, pre- dicting the structure of such complexes from EXAFS spectra is not a simple task. Currently methods such as Monte Carlo optimization or simulated annealing are used in structure refinement of EXAFS. These methods have proven somewhat successful in structure refinement but have not been successful in finding the global minima. Multiple population based algorithms, including a genetic algorithm, a restarting ge- netic algorithm, differential evolution, and particle swarm optimization, are studied for their effectiveness in structure refinement of EXAFS. The oxygen-evolving com- plex in S1 is used as a benchmark for comparing the algorithms. These algorithms were successful in finding new atomic structures that produced improved calculated EXAFS spectra over atomic structures previously found.

Représentation et recherche de motifs cycliques et structuraux d’ARN connus dans les structures secondaires

L'acide désoxyribonucléique (ADN) et l'acide ribonucléique (ARN) sont des polymères de nucléotides essentiels à la cellule. À l'inverse de l'ADN qui sert principalement à stocker l'information génétique, les ARN sont impliqués dans plusieurs processus métaboliques. Par exemple, ils transmettent l’information génétique codée dans l’ADN. Ils sont essentiels pour la maturation des autres ARN, la régulation de l’expression génétique, la prévention de la dégradation des chromosomes et le ciblage des protéines dans la cellule. La polyvalence fonctionnelle de l'ARN résulte de sa plus grande diversité structurale. Notre laboratoire a développé MC-Fold, un algorithme pour prédire la structure des ARN qu'on représente avec des graphes d'interactions inter-nucléotidiques. Les sommets de ces graphes représentent les nucléotides et les arêtes leurs interactions. Notre laboratoire a aussi observé qu'un petit ensemble de cycles d'interactions à lui seul définit la structure de n'importe quel motif d'ARN. La formation de ces cycles dépend de la séquence de nucléotides et MC-Fold détermine les cycles les plus probables étant donnée cette séquence. Mon projet de maîtrise a été, dans un premier temps, de définir une base de données des motifs structuraux et fonctionnels d'ARN, bdMotifs, en terme de ces cycles. Par la suite, j’ai implanté un algorithme, MC-Motifs, qui recherche ces motifs dans des graphes d'interactions et, entre autres, ceux générés par MC-Fold. Finalement, j’ai validé mon algorithme sur des ARN dont la structure est connue, tels que les ARN ribosomaux (ARNr) 5S, 16S et 23S, et l'ARN utilisé pour prédire la structure des riborégulateurs. Le mémoire est divisé en cinq chapitres. Le premier chapitre présente la structure chimique, les fonctions cellulaires de l'ARN et le repliement structural du polymère. Dans le deuxième chapitre, je décris la base de données bdMotifs. Dans le troisième chapitre, l’algorithme de recherche MC-Motifs est introduit. Le quatrième chapitre présente les résultats de la validation et des prédictions. Finalement, le dernier chapitre porte sur la discussion des résultats suivis d’une conclusion sur le travail.

Structural insights into the allosteric activation mechanism of cGMP-dependent protein kinase I

Cyclic GMP-dependent protein kinase (PKG) is a key transducer in the NO-cGMP signaling pathway. In this line, PKG has been considered an important drug target for treating hypertensive cardiovascular and pulmonary diseases. However, the investigation of PKG’s allosteric activation mechanism has been hampered by a lack of structural information. One of the fundamental questions on the cGMP-dependent activation of PKG is how the enzyme can distinguish cGMP over cAMP and selectively respond to cGMP. To ensure proper signaling, PKG must have developed unique features to ensure its activation upon the right activation signal. In this thesis, the cGMP-selective activation mechanism of PKG was studied through determining crystal structures of three truncated constructs of the regulatory domain [CNB-A (92-227), CNB-B (271-369), and CNB-A/B (92-351)] of PKG Iβ in the absence or presence of cyclic nucleotides. Herein, two individual CNB domain structures with biochemical data revealed that the C-terminal CNB domain (CNB-B) is responsible for cGMP selectivity, while the N-terminal CNB-domain (CNB-A) has a higher binding affinity for both cGMP and cAMP without showing any selectivity. Based on these crystal structures, mutagenesis studies were performed in which the critical residues for cyclic nucleotide selectivity and activation were identified. Furthermore, we discovered that the conformational changes of the C-terminal helix of the CNB-B that bridges between the regulatory and catalytic domains including the hydrophobic capping interaction are crucial for PKG activation. In addition, to observe the global conformation of the activated R-domain, I solved a co-crystal structure of the CNB-A/B with cGMP. Although a monomeric construct was crystallized, the structure displays a dimer. Strikingly, the CNB-A domain and its bound cGMP provide a key interface for this dimeric interaction. Using small angle X-ray scattering (SAXS), the existence of the cGMP-mediated dimeric interface within the CNB domains was confirmed. Furthermore, measuring cGMP-binding affinities (EC50) of the dimeric interface mutants as well as determining activation constants (Ka) revealed that the interface formation is important for PKG activation. To conclude, this thesis study provides a new mechanistic insight in PKG activation along with a newly found interface that can be targeted for designing PKG-specific activity modulators.

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.