Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

Homogeneous assay technologies in drug screening: Quenching Resonance Energy Transfer (QRET) technique

Information gained from the human genome project and improvements in compound synthesizing have increased the number of both therapeutic targets and potential lead compounds. This has evolved a need for better screening techniques to have a capacity to screen number of compound libraries against increasing amount of targets. Radioactivity based assays have been traditionally used in drug screening but the fluorescence based assays have become more popular in high throughput screening (HTS) as they avoid safety and waste problems confronted with radioactivity. In comparison to conventional fluorescence more sensitive detection is obtained with time-resolved luminescence which has increased the popularity of time-resolved fluorescence resonance energy transfer (TR-FRET) based assays. To simplify the current TR-FRET based assay concept the luminometric homogeneous single-label utilizing assay technique, Quenching Resonance Energy Transfer (QRET), was developed. The technique utilizes soluble quencher to quench non-specifically the signal of unbound fraction of lanthanide labeled ligand. One labeling procedure and fewer manipulation steps in the assay concept are saving resources. The QRET technique is suitable for both biochemical and cell-based assays as indicated in four studies:1) ligand screening study of β2 -adrenergic receptor (cell-based), 2) activation study of Gs-/Gi-protein coupled receptors by measuring intracellular concentration of cyclic adenosine monophosphate (cell-based), 3) activation study of G-protein coupled receptors by observing the binding of guanosine-5’-triphosphate (cell membranes), and 4) activation study of small GTP binding protein Ras (biochemical). Signal-to-background ratios were between 2.4 to 10 and coefficient of variation varied from 0.5 to 17% indicating their suitability to HTS use.

Synechococcus sp. PCC 7002 CpcB lyase null mutations alter phycocyanin chromophore function and, consequently, affect the redistribution of excitation energy via the light state transition /

Phycobilisomes are the major light harvesting complexes for cyanobacteria and phycocyanin is the primary phycobiliprotein of the phycobilisome rod. The phycocyanobilin lyases responsible for chromophorylating the phycocyanin p subunit (CpcB) have been recently identified in the cyanobacterium Synechococcus sp. PCC 7002. Surprisingly, mutants missing the CpcB lyases were nevertheless capable of producing pigmented phycocyanin. 10K absorbance measurements revealed that the energy states of the p phycocyanin chromophores were only subtly shifted; however, 77K steady state fluorescence emission spectroscopy showed excitation energy transfer involving the targeted chromophores to be highly disrupted. Such evidence suggests that phycobilin orientation within the binding domain is specifically modified. We hypothesized that alternate, less specific lyases are able to act on the p binding sites. A phycocyanin linker-polypeptide deficient mutant was similarly characterized. The light state transition, a short term adaptation of the photosynthetic light harvesting apparatus resulting in the redistribution of excitation energy among the photo systems, was shown to be dominated by the reallocation of phycocyanin-absorbed excitation energy. Treatment with a high M phosphate buffer effectively prevented the redistribution of both chlorophyll a- and phycobilisome- absorbed excitation energy, suggesting that the two effects are not strictly independent. The mutant strains required a larger redistribution of excitation energy between light states, perhaps to compensate for their loss in phycobilisome antenna function.

Orientation of Non-Native 2-Monosubstituted and 2,3-Disubstituted 1,4-Naphthoquinones in the A1 binding site of PSI and the effect on the rate of electron transfer : an electron paramagnetic resonance spectroscopy study

The proce-ss ofoxygenic photosynthesis is vital to life on Earth. the central event in photosynthesis is light induced electron transfer that converts light into energy for growth. Ofparticular significance is the membrane bound multisubunit protein known as Photosystem I (PSI). PSI is a reaction centre that is responsible for the transfer of electrons across the membrane to reduce NADP+ to NADPH. The recent publication ofa high resolution X-ray structure of PSI has shown new information about the structure, in particular the electron transfer cofactors, which allows us to study it in more detail. In PSI, the secondary acceptor is crucial for forward electron transfer. In this thesis, the effect of removing the native acceptor phylloquinone and replacing it with a series of structurally related quinones was investigated via transient electron paramagnetic resonance (EPR) experiments. The orientation of non native quinones in the binding site and their ability to function in the electron transfer process was determined. It was found that PSI will readily accept alkyl naphthoquinones and anthraquinone. Q band EPR experiments revealed that the non-native quinones are incorporated into the binding site with the same orientation of the headgroup as in the native system. X band EPR spectra and deuteration experiments indicate that monosubstituted naphthoquinones are bound to the Al site with their side group in the position occupied by the methyl group in native PSI (meta to the hydrogen bonded carbonyl oxygen). X band EPR experiments show that 2, 3- disubstituted methyl naphthoquinones are also incorporated into the Al site in the same orientation as phylloquinone, even with the presence of a halogen- or sulfur-containing side chain in the position normally occupied by the phytyl tail ofphylloquinone. The exception to this is 2-bromo-3-methyl --.- _. -. - -- - - 4 _._ _ _ - _ _ naphthoquinone which has a poorly resolved spectrum, making determination of the orientation difficuh. All of the non-native quinones studied act as efficient electron acceptors. However, forward electron transfer past the quinone could only be demonstrated for anthraquinone, which has a more negative midpoint potential than phylloquinone. In the case of anthraquinone, an increased rate of forward electron transfer compared to native PSI was found. From these results we can conclude that the rate ofelectron transfer from Al to Fx in native PSI lies in the normal region ofthe Marcus Curve.

The Ligand and Membrane-binding Behaviour of the Phosphatidylinositol Transfer Proteins (PITPa & PITPb)

Human Class I phosphatidylinositol transfer proteins (PITPs) exists in two forms: PITPα and PITPβ. PITPs are believed to be lipid transfer proteins based on their capacity to transfer either phosphatidylinositol (PI) or phosphatidylcholine (PC) between membrane compartments in vitro. In Drosophila, the PITP domain is found to be part of a multi-domain protein named retinal degeneration B (RdgBα). The PITP domain of RdgBα shares 40 % sequence identity with PITPα and has been shown to possess PI and PC binding and transfer activity. The detailed molecular mechanism of ligand transfer by the human PITPs and the Drosophila PITP domain remains to be fully established. Here, we investigated the membrane interactions of these proteins using dual polarization interferometry (DPI). DPI is a technique that measures protein binding affinity to a flat immobilized lipid bilayer. In addition, we also measured how quickly these proteins transfer their ligands to lipid vesicles using a fluorescence resonance energy transfer (FRET)-based assay. DPI investigations suggest that PITPβ had a two-fold higher affinity for membranes compared to PITPα. This was reflected by a four-fold faster ligand transfer rate for PITPβ in comparison to PITPα as determined by the FRET assay. Interestingly, DPI analysis also demonstrated that PI-bound human PITPs have lower membrane affinity compared to PC-bound PITPs. In addition, the FRET studies demonstrated the significance of membrane curvature in the ligand transfer rate of PITPs. The ligand transfer rate was higher when the accepting vesicles were highly curved. Furthermore, when the accepting vesicles contained phosphatidic acid (PA) which have smaller head groups, the transfer rate increased. In contrast, when the accepting vesicles contained phosphoinositides which have larger head groups, the transfer rate was diminished. However, PI, the favorite ligand of PITPs, or the presence of anionic lipids did not appear to influence the ligand transfer rate of PITPs. Both DPI and FRET examinations revealed that the PITP domain of RdgBα was able to bind to membranes. However, the RdgBα PITP domain appears to be a poor binder and transporter of PC.

Accumulated metal speciation in earthworm populations with multi-generational exposure to metalliferous soils: cell fractionation and high energy synchrotron analyses

Predicting metal bioaccumulation and toxicity in soil organisms is complicated by site-specific biotic and abiotic parameters. In this study we exploited tissue fractionation and digestion techniques, combined with X-ray absorption spectroscopy (XAS), to investigate the whole-body and subcellular distributions, ligand affinities, and coordination chemistry of accumulated Pb and Zn in field populations of the epigeic earthworm Lumbricus rubellus inhabiting three contrasting metalliferous and two unpolluted soils. Our main findings were (i) earthworms were resident in soils with concentrations of Pb and Zn ranging from 1200 to 27 000 mg kg(-1) and 200 to 34 000 mg kg(-1), respectively; (ii) Pb and Zn primarily accumulated in the posterior alimentary canal in nonsoluble subcellular fractions of earthworms; (iii) site-specific differences in the tissue and subcellular partitioning profiles of populations were observed, with earthworms from a calcareous site partitioning proportionally more Pb to their anterior body segments and Zn to the chloragosome-rich subcellular fraction than their acidic-soil inhabiting counterparts; (iv) XAS indicated that the interpopulation differences in metal partitioning between organs were not accompanied by qualitative differences in ligand-binding speciation, because crystalline phosphate-containing pyromorphite was a predominant chemical species in the whole-worm tissues of all mine soil residents. Differences in metal (Pb, Zn) partitioning at both organ and cellular levels displayed by field populations with protracted histories of metal exposures may reflect their innate ecophysiological responses to essential edaphic variables, such as Ca2+ status. These observations are highly significant in the challenging exercise of interpreting holistic biomarker data delivered by "omic" technologies.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

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.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

Dual role of the beta(2)-adrenergic receptor C terminus for the binding of beta-arrestin and receptor internalization

Homologous desensitization of beta(2)-adrenergic and other G-protein-coupled receptors is a two-step process. After phosphorylation of agonist-occupied receptors by G-protein-coupled receptor kinases, they bind beta-arrestins, which triggers desensitization and internalization of the receptors. Because it is not known which regions of the receptor are recognized by beta-arrestins, we have investigated beta-arrestin interaction and internalization of a set of mutants of the human beta(2)-adrenergic receptor. Mutation of the four serine/threonine residues between residues 355 and 364 led to the loss of agonist-induced receptor-beta-arrestin2 interaction as revealed by fluorescence resonance energy transfer (FRET), translocation of beta-arrestin2 to the plasma membrane, and receptor internalization. Mutation of all seven serine/threonine residues distal to residue 381 did not affect agonist-induced receptor internalization and beta-arrestin2 translocation. A beta(2)-adrenergic receptor truncated distal to residue 381 interacted normally with beta-arrestin2, whereas its ability to internalize in an agonist-dependent manner was compromised. A similar impairment of internalization was observed when only the last eight residues of the C terminus were deleted. Our experiments show that the C terminus distal to residue 381 does not affect the initial interaction between receptor and beta-arrestin, but its last eight amino acids facilitate receptor internalization in concert with beta-arrestin2.

beta-arrestin binding to the beta(2)-adrenergic receptor requires both receptor phosphorylation and receptor activation

Homologous desensitization of beta(2)-adrenergic receptors has been shown to be mediated by phosphorylation of the agonist-stimulated receptor by G-protein-coupled receptor kinase 2 (GRK2) followed by binding of beta-arrestins to the phosphorylated receptor. Binding of beta-arrestin to the receptor is a prerequisite for subsequent receptor desensitization, internalization via clathrin-coated pits, and the initiation of alternative signaling pathways. In this study we have investigated the interactions between receptors and beta-arrestin2 in living cells using fluorescence resonance energy transfer. We show that (a) the initial kinetics of beta-arrestin2 binding to the receptor is limited by the kinetics of GRK2-mediated receptor phosphorylation; (b) repeated stimulation leads to the accumulation of GRK2-phosphorylated receptor, which can bind beta-arrestin2 very rapidly; and (c) the interaction of beta-arrestin2 with the receptor depends on the activation of the receptor by agonist because agonist withdrawal leads to swift dissociation of the receptor-beta-arrestin2 complex. This fast agonist-controlled association and dissociation of beta-arrestins from prephosphorylated receptors should permit rapid control of receptor sensitivity in repeatedly stimulated cells such as neurons.

Supranutritional selenium induces alterations in molecular targets related to energy metabolism in skeletal muscle and visceral adipose tissue of pigs

While selenium (Se) is an essential micronutrient for humans, epidemiological studies have raised concern that supranutritional Se intake may increase the risk to develop Type 2 diabetes mellitus (T2DM). We aimed to determine the impact of Se at a dose and source frequently ingested by humans on markers of insulin sensitivity and signalling. Male pigs were fed either a Se-adequate (0.17 mg Se/kg) or a Se-supranutritional (0.50 mg Se/kg; high-Se) diet. After 16 weeks of intervention, fasting plasma insulin and cholesterol levels were non-significantly increased in the high-Se pigs, whereas fasting glucose concentrations did not differ between the two groups. In skeletal muscle of high-Se pigs, glutathione peroxidase activity was increased, gene expression of forkhead box O1 transcription factor and peroxisomal proliferator-activated receptor- coactivator 1 were increased and gene expression of the glycolytic enzyme pyruvate kinase was decreased. In visceral adipose tissue of high-Se pigs, mRNA levels of sterol regulatory element-binding transcription factor 1 were increased, and the phosphorylation of Akt, AMP-activated kinase and mitogen-activated protein kinases was affected. In conclusion, dietary Se oversupply may affect expression and activity of proteins involved in energy metabolism in major insulin target tissues, though this is probably not sufficient to induce diabetes.

Proteins and their interacting partners: an introduction to protein–ligand binding site prediction methods

Elucidating the biological and biochemical roles of proteins, and subsequently determining their interacting partners, can be difficult and time consuming using in vitro and/or in vivo methods, and consequently the majority of newly sequenced proteins will have unknown structures and functions. However, in silico methods for predicting protein–ligand binding sites and protein biochemical functions offer an alternative practical solution. The characterisation of protein–ligand binding sites is essential for investigating new functional roles, which can impact the major biological research spheres of health, food, and energy security. In this review we discuss the role in silico methods play in 3D modelling of protein–ligand binding sites, along with their role in predicting biochemical functionality. In addition, we describe in detail some of the key alternative in silico prediction approaches that are available, as well as discussing the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated Model EvaluatiOn (CAMEO) projects, and their impact on developments in the field. Furthermore, we discuss the importance of protein function prediction methods for tackling 21st century problems.

Mapping the Intramolecular Vibrational Energy Flow in Proteins Reveals Functionally Important Residues

Unveiling the mechanisms of energy relaxation in biomolecules is key to our understanding of protein stability, allostery, intramolecular signaling, and long-lasting quantum coherence phenomena at ambient temperatures. Yet, the relationship between the pathways of energy transfer and the functional role of the residues involved remains largely unknown. Here, we develop a simulation method of mapping out residues that are highly efficient in relaxing an initially localized excess vibrational energy and perform site-directed mutagenesis functional assays to assess the relevance of these residues to protein function. We use the ligand binding domains of thyroid hormone receptor (TR) subtypes as a test case and find that conserved arginines, which are critical to TR transactivation function, are the most effective heat diffusers across the protein structure. These results suggest a hitherto unsuspected connection between a residue`s ability to mediate intramolecular vibrational energy redistribution and its functional relevance.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.