Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Revisiting G3BP1 as a RasGAP binding protein: sensitization of tumor cells to chemotherapy by the RasGAP 317-326 sequence does not involve G3BP1.

RasGAP is a multifunctional protein that controls Ras activity and that is found in chromosomal passenger complexes. It also negatively or positively regulates apoptosis depending on the extent of its cleavage by caspase-3. RasGAP has been reported to bind to G3BP1 (RasGAP SH3-domain-binding protein 1), a protein regulating mRNA stability and stress granule formation. The region of RasGAP (amino acids 317-326) thought to bind to G3BP1 corresponds exactly to the sequence within fragment N2, a caspase-3-generated fragment of RasGAP, that mediates sensitization of tumor cells to genotoxins. While assessing the contribution of G3BP1 in the anti-cancer function of a cell-permeable peptide containing the 317-326 sequence of RasGAP (TAT-RasGAP₃₁₇₋₃₂₆), we found that, in conditions where G3BP1 and RasGAP bind to known partners, no interaction between G3BP1 and RasGAP could be detected. TAT-RasGAP₃₁₇₋₃₂₆ did not modulate binding of G3BP1 to USP10, stress granule formation or c-myc mRNA levels. Finally, TAT-RasGAP₃₁₇₋₃₂₆ was able to sensitize G3BP1 knock-out cells to cisplatin-induced apoptosis. Collectively these results indicate that G3BP1 and its putative RasGAP binding region have no functional influence on each other. Importantly, our data provide arguments against G3BP1 being a genuine RasGAP-binding partner. Hence, G3BP1-mediated signaling may not involve RasGAP.

A role for septins in the interaction between the Listeria monocytogenes INVASION PROTEIN InlB and the Met receptor.

Septins are conserved GTPases that form filaments and are required for cell division. During interphase, septin filaments associate with cellular membrane and cytoskeleton networks, yet the functional significance of these associations have, to our knowledge, remained unknown. We recently discovered that different septins, SEPT2 and SEPT11, regulate the InlB-mediated entry of Listeria monocytogenes into host cells. Here we address the role of SEPT2 and SEPT11 in the InlB-Met interactions underlying Listeria invasion to explore how septins modulate surface receptor function. We observed that differences in InlB-mediated Listeria entry correlated with differences in Met surface expression caused by septin depletion. Using atomic force microscopy on living cells, we show that septin depletion significantly reduced the unbinding force of InlB-Met interaction and the viscosity of membrane tethers at locations where the InlB-Met interaction occurs. Strikingly, the same order of difference was observed for cells in which the actin cytoskeleton was disrupted. Consistent with a proposed role of septins in association with the actin cytoskeleton, we show that cell elasticity is decreased upon septin or actin inactivation. Septins are therefore likely to participate in anchorage of the Met receptor to the actin cytoskeleton, and represent a critical determinant in surface receptor function.

Clic4, a novel protein that sensitizes ß-cells to apoptosis.

Introduction: La préservation et/ou l'expansion de la masse des cellules ß pourraient constituer des approches prometteuses dans le traitement du diabète. L'une des stratégies clés serait de réduire l'apoptose des cellules ß. Le chloride intracellular channel protein 4 (Clic4) est une protéine exprimée de manière ubiquitaire et supposée agir dans de nombreux processus cellulaires tels que le contrôle du cycle cellulaire, la différenciation cellulaire et l'apoptose. Ici, nous avons étudié le rôle de Clic4 dans l'apoptose des cellules ß pancréatiques en utilisant des cellules ßTC-tet et des îlots de Langerhans issus de souris knockout pour Clic4 (ßClic4KO). Résultats: L'expression de l'ARNm et de la protéine Clic4 était augmentée par un traitement aux cytokines dans les cellules ßTC-tet et encore plus fortement dans des îlots isolés de souris. De plus, la sous-expression de Clic4 dans les cellules ßTC-tet diminuait leur sensibilité à l'apoptose induite par les cytokines. La sous-expression de Clic4 dans les cellules ßTC-tet n'affectait pas l'expression des ARNm de Bcl-2 et Bad, mais augmentait leur expression protéique ainsi que la forme phosphorylée de Bad. Les mêmes résultats ont été obtenus sur des îlots isolés de souris contrôles et ßClic4KO. De plus, les îlots issus de souris ßClic4KO présentaient une augmentation de l'expression de la protéine Bcl-xL. Dans le but de déterminer si Clic4 augmentait l'expression de Bcl-2 et Bad via une interaction protéique directe, nous avons immunoprécipité Clic4 à partir de cellules ßTC-tet à l'aide d'anticorps dirigés contre la partie C- ou N-terminale de la protéine, puis nous avons soumis les immunoprécipités à une analyse de spectrométrie de masse. Aucune co- immunoprécipitation avec Bcl-2 ou d'autres protéines de la famille Bcl-2 n'a été détectée. Cependant, de manière intéressante, Clic4 était co-purifié avec plusieurs protéines du protéasome suggérant un rôle de Clic4 dans la dégradation des protéines. Par conséquent, nous avons étudié la demi-vie de Bcl-2 et Bad, et avons observé que la sous-expression de Clic4 dans les cellules ßTC-tet augmentait la demi-vie de ces protéines. De plus, l'expression de l'ARNm et de la protéine Clic4 était également augmentée lors d'un stress du réticulum endoplasmique induit par la thapsigargine dans les cellules ßTC-tet. La sous-expression de Clic4 dans les cellules ßTC-tet ou chez les KO diminuait la sensibilité des cellules ß à l'apoptose induite par la thapsigargine ou l'acide palmitique, respectivement. Conclusion: Ces résultats suggèrent que Clic4 sensibilise les cellules ß à l'apoptose induite par les cytokines ou l'acide palmitique/thapsigargine (stress du réticulum endoplasmique). De plus, la sous-expression de Clic4 améliore la survie des cellules ß en diminuant la dégradation de Bcl-2 et Bcl-xL, et en augmentant le niveau total de Bad phosphorylé, peut-être suite à une interaction de Clic4 avec le protéasome.

KIAA1985, a Protein Mutant in Charcot-Marie-Tooth Neuropathy, Links Peripheral Nerve Myelination to Endosomal Recycling Pathways

Charcot-Marie-Tooth neuropathy (CMT) represents a heterogenous group of inherited disorders of the peripheral nervous system. One form of autosomal recessive demyelinating CMT (CMT4C, 5q32) is caused by mutations in the gene encoding KIAA1985, a protein of so far unknown function. Here we show that KIAA1985 is exclusively expressed in Schwann cells. KIAA1985 is tethered to cellular membranes through an N-terminal myristic acid anchor and localizes to the perinuclear recycling compartment. A search for proteins that interact with KIAA1985 identified the small GTPase Rab11, a key regulator of recycling endosome functions. CMT4C-related missense mutations disrupt the KIAA1985/Rab11 interaction. Protein binding studies indicate that KIAA1985 functions as a Rab11 effector, as it interacts only with active forms of Rab11 (WT and Q70L) and does not interact with the GDP locked mutant (S25N). Consistent with a function of Rab11 in Schwann cell myelination, myelin formation was strongly impaired when dorsal root ganglion neurons were co-cultured with Schwann cells infected with Rab11 S25N. Our data indicate that the KIAA1985/Rab11 interaction is relevant for peripheral nerve pathophysiology and place endosomal recycling on the list of cellular mechanisms involved in Schwann cell myelination.

Overexpression of guanylate cyclase activating protein 2 in rod photoreceptors in vivo leads to morphological changes at the synaptic ribbon.

Guanylate cyclase activating proteins are EF-hand containing proteins that confer calcium sensitivity to retinal guanylate cyclase at the outer segment discs of photoreceptor cells. By making the rate of cGMP synthesis dependent on the free intracellular calcium levels set by illumination, GCAPs play a fundamental role in the recovery of the light response and light adaptation. The main isoforms GCAP1 and GCAP2 also localize to the synaptic terminal, where their function is not known. Based on the reported interaction of GCAP2 with Ribeye, the major component of synaptic ribbons, it was proposed that GCAP2 could mediate the synaptic ribbon dynamic changes that happen in response to light. We here present a thorough ultrastructural analysis of rod synaptic terminals in loss-of-function (GCAP1/GCAP2 double knockout) and gain-of-function (transgenic overexpression) mouse models of GCAP2. Rod synaptic ribbons in GCAPs−/− mice did not differ from wildtype ribbons when mice were raised in constant darkness, indicating that GCAPs are not required for ribbon early assembly or maturation. Transgenic overexpression of GCAP2 in rods led to a shortening of synaptic ribbons, and to a higher than normal percentage of club-shaped and spherical ribbon morphologies. Restoration of GCAP2 expression in the GCAPs−/− background (GCAP2 expression in the absence of endogenous GCAP1) had the striking result of shortening ribbon length to a much higher degree than overexpression of GCAP2 in the wildtype background, as well as reducing the thickness of the outer plexiform layer without affecting the number of rod photoreceptor cells. These results indicate that preservation of the GCAP1 to GCAP2 relative levels is relevant for maintaining the integrity of the synaptic terminal. Our demonstration of GCAP2 immunolocalization at synaptic ribbons at the ultrastructural level would support a role of GCAPs at mediating the effect of light on morphological remodeling changes of synaptic ribbons.

Activity and functional interaction of alternative oxidase and uncoupling protein in mitochondria from tomato fruit

Cyanide-resistant alternative oxidase (AOX) is not limited to plant mitochondria and is widespread among several types of protists. The uncoupling protein (UCP) is much more widespread than previously believed, not only in tissues of higher animals but also in plants and in an amoeboid protozoan. The redox energy-dissipating pathway (AOX) and the proton electrochemical gradient energy-dissipating pathway (UCP) lead to the same final effect, i.e., a decrease in ATP synthesis and an increase in heat production. Studies with green tomato fruit mitochondria show that both proteins are present simultaneously in the membrane. This raises the question of a specific physiological role for each energy-dissipating system and of a possible functional connection between them (shared regulation). Linoleic acid, an abundant free fatty acid in plants which activates UCP, strongly inhibits cyanide-resistant respiration mediated by AOX. Moreover, studies of the evolution of AOX and UCP protein expression and of their activities during post-harvest ripening of tomato fruit show that AOX and plant UCP work sequentially: AOX activity decreases in early post-growing stages and UCP activity is decreased in late ripening stages. Electron partitioning between the alternative oxidase and the cytochrome pathway as well as H+ gradient partitioning between ATP synthase and UCP can be evaluated by the ADP/O method. This method facilitates description of the kinetics of energy-dissipating pathways and of ATP synthase when state 3 respiration is decreased by limitation of oxidizable substrate.

The interaction of housing condition and acute immobilization stress on the elevated plus-maze behaviors of protein-malnourished rats

Protein malnutrition induces structural, neurochemical and functional alterations in the central nervous system, leading to behavioral alterations. In the present study, we used the elevated plus-maze (EPM) as a measure of anxiety to evaluate the interaction between acute immobilization and housing conditions on the behavior of malnourished rats. Pups (6 males and 2 females) were fed by Wistar lactating dams receiving a 6% (undernourished) or 16% (well-nourished) protein diet. After weaning, the animals continued to receive the same diets ad libitum until 49 days of age when they started to receive a regular lab chow diet. From weaning to the end of the tests on day 70, the animals were housed under two different conditions, i.e., individual or in groups of three. On the 69th day, half of the animals were submitted to immobilization for 2 h, while the other half were undisturbed, and both groups were tested 24 h later for 5 min in the EPM. Independent of other factors, protein malnutrition increased, while immobilization and social isolation per se decreased, EPM exploration. Analysis of the interaction of diet vs immobilization vs housing conditions showed that the increased EPM exploration presented by the malnourished group was reversed by acute immobilization in animals reared in groups but not in animals reared individually. The interaction between immobilization and housing conditions suggests that living for a long time in social isolation is sufficiently stressful to reduce the responses to another anxiogenic procedure (immobilization), while living in groups prompts the animals to react to acute stress. Thus, it is suggested that housing condition can modulate the effects of an anxiogenic procedure on behavioral responses of malnourished rats in the EPM.

Interaction between human cytomegalovirus UL136 protein and ATP1B1 protein

Interplay between the host and human cytomegalovirus (HCMV) has a pivotal role in the outcome of infection. A region (referred to as UL/b’) present in the Toledo strain of HCMV and low passage clinical isolates contains 19 additional genes, which are absent in the highly passaged laboratory strain AD169. Products of the UL/b’ genes may determine the manifestations of HCMV infection in vivo. However, little is known about the host factors, which interact with UL/b’ proteins. This study was conducted to investigate the function of the HCMV UL136 protein. By yeast two-hybrid screening, the β1 subunit of the host Na+/K+-ATPase (ATP1B1) was identified to be a candidate protein, which interacts with the HCMV UL136 protein. The interaction was further evaluated both in vitro by pull-down assay and in vivo by immunofluorescent co-localization. The results showed that the UL136 protein can interact with ATP1B1 in vitro. Co-localization of UL136-EGFP and ATP1B1-DsRed in cell membranes suggests that ATP1B1 was a partner of the UL136 protein. It can be proposed that the HCMV UL136 protein may have important roles in processes such as cell-to-cell spread, and in maintaining cell osmotic pressure and intracellular ion homeostasis during HCMV infection.

Effect of ATP and 2-oxoglutarate on the in vitro interaction between the NifA GAF domain and the GlnB protein of Azospirillum brasilense

Azospirillum brasilense is a diazotroph that associates with important agricultural crops and thus has potential to be a nitrogen biofertilizer. The A. brasilense transcription regulator NifA, which seems to be constitutively expressed, activates the transcription of nitrogen fixation genes. It has been suggested that the nitrogen status-signaling protein GlnB regulates NifA activity by direct interaction with the NifA N-terminal GAF domain, preventing the inhibitory effect of this domain under conditions of nitrogen fixation. In the present study, we show that an N-terminal truncated form of NifA no longer required GlnB for activity and lost regulation by ammonium. On the other hand, in trans co-expression of the N-terminal GAF domain inhibited the N-truncated protein in response to fixed nitrogen levels. We also used pull-down assays to show in vitro interaction between the purified N-terminal GAF domain of NifA and the GlnB protein. The results showed that A. brasilense GlnB interacts directly with the NifA N-terminal domain and this interaction is dependent on the presence of ATP and 2-oxoglutarate.

Stromal interaction molecule 1 and its role in the regulation, proliferation and protein expression in follicular ML-1 thyroid cancer cells

Calcium (Ca2+) is involved in the regulation of variety of cellular functions including hallmarks of cancer development such as cellular migration and cellular proliferation. Store-operated calcium entry (SOCE) is a central mechanism in cellular calcium signaling and in maintaining the cellular calcium balance. Stromal interaction molecule 1(STIM1) has been identified as an important constituent of SOCE. In this thesis , the STIM1 proteins are studied for their importance in cellular processes and their effects on the expression of S1P1, S1P2, S1P3, VEGFR-2, and TRPC-1 in follicular ML-1 thyroid cancer cells. The results show the importance of STIM1 proteins in SOCE in these cells. The SOCE is significantly reduced in the STIM1 knockdown cells. The results also show the importance of STIM1 proteins in the expression of S1P2 and VEGFR-2 in these cells, as knockdown of STIM1 was shown to upregulate the expression of S1P2 and VEGFR-2. The migration and proliferation is also considerably reduced in the cells in which STIM1 has been knocked down showing the significance of STIM1 in the migration and proliferation in these cells.

From protein structure to function with bioinformatics

It has long been known that amino acids are the building blocks for proteins and govern their folding into specific three-dimensional structures. However, the details of this process are still unknown and represent one of the main problems in structural bioinformatics, which is a highly active research area with the focus on the prediction of three-dimensional structure and its relationship to protein function. The protein structure prediction procedure encompasses several different steps from searches and analyses of sequences and structures, through sequence alignment to the creation of the structural model. Careful evaluation and analysis ultimately results in a hypothetical structure, which can be used to study biological phenomena in, for example, research at the molecular level, biotechnology and especially in drug discovery and development. In this thesis, the structures of five proteins were modeled with templatebased methods, which use proteins with known structures (templates) to model related or structurally similar proteins. The resulting models were an important asset for the interpretation and explanation of biological phenomena, such as amino acids and interaction networks that are essential for the function and/or ligand specificity of the studied proteins. The five proteins represent different case studies with their own challenges like varying template availability, which resulted in a different structure prediction process. This thesis presents the techniques and considerations, which should be taken into account in the modeling procedure to overcome limitations and produce a hypothetical and reliable three-dimensional structure. As each project shows, the reliability is highly dependent on the extensive incorporation of experimental data or known literature and, although experimental verification of in silico results is always desirable to increase the reliability, the presented projects show that also the experimental studies can greatly benefit from structural models. With the help of in silico studies, the experiments can be targeted and precisely designed, thereby saving both money and time. As the programs used in structural bioinformatics are constantly improved and the range of templates increases through structural genomics efforts, the mutual benefits between in silico and experimental studies become even more prominent. Hence, reliable models for protein three-dimensional structures achieved through careful planning and thoughtful executions are, and will continue to be, valuable and indispensable sources for structural information to be combined with functional data.

Iron-regulated surface determinant protein a mediates adhesion of staphylococcus aureus to human corneocyte envelope proteins

The ability of Staphylococcus aureus to colonize the human nares is a crucial prerequisite for disease. IsdA is a major S. aureus surface protein that is expressed during human infection and required for nasal colonization and survival on human skin. In this work, we show that IsdA binds to involucrin, loricrin, and cytokeratin K10, proteins that are present in the cornified envelope of human desquamated epithelial cells. To measure the forces and dynamics of the interaction between IsdA and loricrin (the most abundant protein of the cornified envelope), single-molecule force spectroscopy was used, demonstrating high-specificity binding. IsdA acts as a cellular adhesin to the human ligands, promoting whole-cell binding to immobilized proteins, even in the absence of other S. aureus components (as shown by heterologous expression in Lactococcus lactis). Inhibition experiments revealed the binding of the human ligands to the same IsdA region. This region was mapped to the NEAT domain of IsdA. The NEAT domain also was found to be required for S. aureus whole-cell binding to the ligands as well as to human nasal cells. Thus, IsdA is an important adhesin to human ligands, which predominate in its primary ecological niche.

The structure of echovirus type 12 bound to a two-domain fragment of its cellular attachment protein decay-accelerating factor (CD 55)

Echovirus type 12 (EV12), an enterovirus of the Picornaviridae family, uses the complement regulator, decay-accelerating factor (DAF, CD55) as a cellular receptor. We have calculated a three-dimensional reconstruction of EV12 bound to a fragment of DAF, consisting of short consensus repeat domains 3 and 4, from cryo-negative stain electron microscopy data (EMD #1057). This shows that, as for an earlier reconstruction of the related echovirus type 7 bound to DAF, attachment is not within the viral canyon but occurs close to the two-fold symmetry axes. Despite this general similarity, our reconstruction reveals a receptor interaction that is quite different from that observed for EV7. Fitting of the crystallographic co-ordinates for DAF34 and EV11 into the reconstruction shows a close agreement between the crystal structure of the receptor fragment and the density for the virus-bound receptor, allowing unambiguous positioning of the receptor with respect to the virion (PDB #1UPN). Our finding that the mode of virus-receptor interaction in EV12 is distinct from that seen for EV7 raises interesting questions regarding the evolution and biological significance of the DAF-binding phenotype in these viruses.

Glucokinase Regulatory Protein genetic variant interacts with omega-3 PUFA to influence insulin resistance and inflammation in metabolic syndrome

Glucokinase Regulatory Protein (GCKR) plays a central role regulating both hepatic triglyceride and glucose metabolism. Fatty acids are key metabolic regulators, which interact with genetic factors and influence glucose metabolism and other metabolic traits. Omega-3 polyunsaturated fatty acids (n-3 PUFA) have been of considerable interest, due to their potential to reduce metabolic syndrome (MetS) risk. Objective To examine whether genetic variability at the GCKR gene locus was associated with the degree of insulin resistance, plasma concentrations of C-reactive protein (CRP) and n-3 PUFA in MetS subjects. Design Homeostasis model assessment of insulin resistance (HOMA-IR), HOMA-B, plasma concentrations of C-peptide, CRP, fatty acid composition and the GCKR rs1260326-P446L polymorphism, were determined in a cross-sectional analysis of 379 subjects with MetS participating in the LIPGENE dietary cohort. Results Among subjects with n-3 PUFA levels below the population median, carriers of the common C/C genotype had higher plasma concentrations of fasting insulin (P = 0.019), C-peptide (P = 0.004), HOMA-IR (P = 0.008) and CRP (P = 0.032) as compared with subjects carrying the minor T-allele (Leu446). In contrast, homozygous C/C carriers with n-3 PUFA levels above the median showed lower plasma concentrations of fasting insulin, peptide C, HOMA-IR and CRP, as compared with individuals with the T-allele. Conclusions We have demonstrated a significant interaction between the GCKR rs1260326-P446L polymorphism and plasma n-3 PUFA levels modulating insulin resistance and inflammatory markers in MetS subjects. Further studies are needed to confirm this gene-diet interaction in the general population and whether targeted dietary recommendations can prevent MetS in genetically susceptible individuals.