Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

Crystallization and preliminary crystallographic characterization of an SH3 domain from the IB1 scaffold protein.

IB1 is a mammalian scaffold protein that interacts with components of the c-Jun N-terminal kinase (JNK) signal-transduction pathway mainly via its protein-protein interaction domains. Crystallization of the key Src homology 3 (SH3) domain of IB1 has been achieved. Crystallization experiments with unmodified protein and deliberately oxidized protein have led to different crystal forms. X-ray data have been collected to 3.0 A resolution from a crystal form with rectangular prism morphology. These crystals are orthorhombic (P2(1)2(1)2(1)), with unit-cell parameters a = 45.9, b = 57.0, c = 145.5 A. These are the first crystallographic data on a scaffold molecule such as IB1 to be reported.

Neutron structure of type-III antifreeze protein allows the reconstruction of AFP-ice interface.

Antifreeze proteins (AFPs) inhibit ice growth at sub-zero temperatures. The prototypical type-III AFPs have been extensively studied, notably by X-ray crystallography, solid-state and solution NMR, and mutagenesis, leading to the identification of a compound ice-binding surface (IBS) composed of two adjacent ice-binding sections, each which binds to particular lattice planes of ice crystals, poisoning their growth. This surface, including many hydrophobic and some hydrophilic residues, has been extensively used to model the interaction of AFP with ice. Experimentally observed water molecules facing the IBS have been used in an attempt to validate these models. However, these trials have been hindered by the limited capability of X-ray crystallography to reliably identify all water molecules of the hydration layer. Due to the strong diffraction signal from both the oxygen and deuterium atoms, neutron diffraction provides a more effective way to determine the water molecule positions (as D(2) O). Here we report the successful structure determination at 293 K of fully perdeuterated type-III AFP by joint X-ray and neutron diffraction providing a very detailed description of the protein and its solvent structure. X-ray data were collected to a resolution of 1.05 Å, and neutron Laue data to a resolution of 1.85 Å with a "radically small" crystal volume of 0.13 mm(3). The identification of a tetrahedral water cluster in nuclear scattering density maps has allowed the reconstruction of the IBS-bound ice crystal primary prismatic face. Analysis of the interactions between the IBS and the bound ice crystal primary prismatic face indicates the role of the hydrophobic residues, which are found to bind inside the holes of the ice surface, thus explaining the specificity of AFPs for ice versus water.

EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

Molecular mechanism of a green-shifted, pH-dependent red fluorescent protein mKate variant.

Fluorescent proteins that can switch between distinct colors have contributed significantly to modern biomedical imaging technologies and molecular cell biology. Here we report the identification and biochemical analysis of a green-shifted red fluorescent protein variant GmKate, produced by the introduction of two mutations into mKate. Although the mutations decrease the overall brightness of the protein, GmKate is subject to pH-dependent, reversible green-to-red color conversion. At physiological pH, GmKate absorbs blue light (445 nm) and emits green fluorescence (525 nm). At pH above 9.0, GmKate absorbs 598 nm light and emits 646 nm, far-red fluorescence, similar to its sequence homolog mNeptune. Based on optical spectra and crystal structures of GmKate in its green and red states, the reversible color transition is attributed to the different protonation states of the cis-chromophore, an interpretation that was confirmed by quantum chemical calculations. Crystal structures reveal potential hydrogen bond networks around the chromophore that may facilitate the protonation switch, and indicate a molecular basis for the unusual bathochromic shift observed at high pH. This study provides mechanistic insights into the color tuning of mKate variants, which may aid the development of green-to-red color-convertible fluorescent sensors, and suggests GmKate as a prototype of genetically encoded pH sensors for biological studies.

Protein pocket and ligand shape comparison and its application in virtual screening.

Understanding molecular recognition is one major requirement for drug discovery and design. Physicochemical and shape complementarity between two binding partners is the driving force during complex formation. In this study, the impact of shape within this process is analyzed. Protein binding pockets and co-crystallized ligands are represented by normalized principal moments of inertia ratios (NPRs). The corresponding descriptor space is triangular, with its corners occupied by spherical, discoid, and elongated shapes. An analysis of a selected set of sc-PDB complexes suggests that pockets and bound ligands avoid spherical shapes, which are, however, prevalent in small unoccupied pockets. Furthermore, a direct shape comparison confirms previous studies that on average only one third of a pocket is filled by its bound ligand, supplemented by a 50 % subpocket coverage. In this study, we found that shape complementary is expressed by low pairwise shape distances in NPR space, short distances between the centers-of-mass, and small deviations in the angle between the first principal ellipsoid axes. Furthermore, it is assessed how different binding pocket parameters are related to bioactivity and binding efficiency of the co-crystallized ligand. In addition, the performance of different shape and size parameters of pockets and ligands is evaluated in a virtual screening scenario performed on four representative targets.

Development of metabolic labeling strategies to study changes in protein expression

Résumé : Les progrès techniques de la spectrométrie de masse (MS) ont contribué au récent développement de la protéomique. Cette technique peut actuellement détecter, identifier et quantifier des milliers de protéines. Toutefois, elle n'est pas encore assez puissante pour fournir une analyse complète des modifications du protéome corrélées à des phénomènes biologiques. Notre objectif était le développement d'une nouvelle stratégie pour la détection spécifique et la quantification des variations du protéome, basée sur la mesure de la synthèse des protéines plutôt que sur celle de la quantité de protéines totale. Pour cela, nous volions associer le marquage pulsé des protéines par des isotopes stables avec une méthode d'acquisition MS basée sur le balayage des ions précurseurs (precursor ion scan, ou PIS), afin de détecter spécifiquement les protéines ayant intégré les isotopes et d'estimer leur abondance par rapport aux protéines non marquées. Une telle approche peut identifier les protéines avec les plus hauts taux de synthèse dans une période de temps donnée, y compris les protéines dont l'expression augmente spécifiquement suite à un événement précis. Nous avons tout d'abord testé différents acides aminés marqués en combinaison avec des méthodes PIS spécifiques. Ces essais ont permis la détection spécifique des protéines marquées. Cependant, en raison des limitations instrumentales du spectromètre de masse utilisé pour les méthodes PIS, la sensibilité de cette approche s'est révélée être inférieure à une analyse non ciblée réalisée sur un instrument plus récent (Chapitre 2.1). Toutefois, pour l'analyse différentielle de deux milieux de culture conditionnés par des cellules cancéreuses humaines, nous avons utilisé le marquage métabolique pour distinguer les protéines d'origine cellulaire des protéines non marquées du sérum présentes dans les milieux de culture (Chapitre 2.2). Parallèlement, nous avons développé une nouvelle méthode de quantification nommée IBIS, qui utilise des paires d'isotopes stables d'acides aminés capables de produire des ions spécifiques qui peuvent être utilisés pour la quantification relative. La méthode IBIS a été appliquée à l'analyse de deux lignées cellulaires cancéreuses complètement marquées, mais de manière différenciée, par des paires d'acides aminés (Chapitre 2.3). Ensuite, conformément à l'objectif initial de cette thèse, nous avons utilisé une variante pulsée de l'IBIS pour détecter des modifications du protéome dans des cellules HeLa infectée par le virus humain Herpes Simplex-1 (Chapitre 2.4). Ce virus réprime la synthèse des protéines des cellules hôtes afin d'exploiter leur mécanisme de traduction pour la production massive de virions. Comme prévu, de hauts taux de synthèse ont été mesurés pour les protéines virales détectées, attestant de leur haut niveau d'expression. Nous avons de plus identifié un certain nombre de protéines humaines dont le rapport de synthèse et de dégradation (S/D) a été modifié par l'infection virale, ce qui peut donner des indications sur les stratégies utilisées par les virus pour détourner la machinerie cellulaire. En conclusion, nous avons montré dans ce travail que le marquage métabolique peut être employé de façon non conventionnelle pour étudier des dimensions peu explorées en protéomique. Summary : In recent years major technical advancements greatly supported the development of mass spectrometry (MS)-based proteomics. Currently, this technique can efficiently detect, identify and quantify thousands of proteins. However, it is not yet sufficiently powerful to provide a comprehensive analysis of the proteome changes correlated with biological phenomena. The aim of our project was the development of ~a new strategy for the specific detection and quantification of proteomé variations based on measurements of protein synthesis rather than total protein amounts. The rationale for this approach was that changes in protein synthesis more closely reflect dynamic cellular responses than changes in total protein concentrations. Our starting idea was to couple "pulsed" stable-isotope labeling of proteins with a specific MS acquisition method based on precursor ion scan (PIS), to specifically detect proteins that incorporated the label and to simultaneously estimate their abundance, relative to the unlabeled protein isoform. Such approach could highlight proteins with the highest synthesis rate in a given time frame, including proteins specifically up-regulated by a given biological stimulus. As a first step, we tested different isotope-labeled amino acids in combination with dedicated PIS methods and showed that this leads to specific detection of labeled proteins. Sensitivity, however, turned out to be lower than an untargeted analysis run on a more recent instrument, due to MS hardware limitations (Chapter 2.1). We next used metabolic labeling to distinguish the proteins of cellular origin from a high background of unlabeled (serum) proteins, for the differential analysis of two serum-containing culture media conditioned by labeled human cancer cells (Chapter 2.2). As a parallel project we developed a new quantification method (named ISIS), which uses pairs of stable-isotope labeled amino acids able to produce specific reporter ions, which can be used for relative quantification. The ISIS method was applied to the analysis of two fully, yet differentially labeled cancer cell lines, as described in Chapter 2.3. Next, in line with the original purpose of this thesis, we used a "pulsed" variant of ISIS to detect proteome changes in HeLa cells after the infection with human Herpes Simplex Virus-1 (Chapter 2.4). This virus is known to repress the synthesis of host cell proteins to exploit the translation machinery for the massive production of virions. As expected, high synthesis rates were measured for the detected viral proteins, confirming their up-regulation. Moreover, we identified a number of human proteins whose synthesis/degradation ratio (S/D) was affected by the viral infection and which could provide clues on the strategies used by the virus to hijack the cellular machinery. Overall, in this work, we showed that metabolic labeling can be employed in alternative ways to investigate poorly explored dimensions in proteomics.

The fourth extracellular loop of the alpha subunit of Na,K-ATPase. Functional evidence for close proximity with the second extracellular loop.

Na,K-ATPase is the main active transport system that maintains the large gradients of Na(+) and K(+) across the plasma membrane of animal cells. The crystal structure of a K(+)-occluding conformation of this protein has been recently published, but the movements of its different domains allowing for the cation pumping mechanism are not yet known. The structure of many more conformations is known for the related calcium ATPase SERCA, but the reliability of homology modeling is poor for several domains with low sequence identity, in particular the extracellular loops. To better define the structure of the large fourth extracellular loop between the seventh and eighth transmembrane segments of the alpha subunit, we have studied the formation of a disulfide bond between pairs of cysteine residues introduced by site-directed mutagenesis in the second and the fourth extracellular loop. We found a specific pair of cysteine positions (Y308C and D884C) for which extracellular treatment with an oxidizing agent inhibited the Na,K pump function, which could be rapidly restored by a reducing agent. The formation of the disulfide bond occurred preferentially under the E2-P conformation of Na,K-ATPase, in the absence of extracellular cations. Using recently published crystal structure and a distance constraint reproducing the existence of disulfide bond, we performed an extensive conformational space search using simulated annealing and showed that the Tyr(308) and Asp(884) residues can be in close proximity, and simultaneously, the SYGQ motif of the fourth extracellular loop, known to interact with the extracellular domain of the beta subunit, can be exposed to the exterior of the protein and can easily interact with the beta subunit.

Proteomic and Metabolomic Analyses of Mitochondrial Complex I-deficient Mouse Model Generated by Spontaneous B2 Short Interspersed Nuclear Element (SINE) Insertion into NADH Dehydrogenase (Ubiquinone) Fe-S Protein 4 (Ndufs4) Gene.

Eukaryotic cells generate energy in the form of ATP, through a network of mitochondrial complexes and electron carriers known as the oxidative phosphorylation system. In mammals, mitochondrial complex I (CI) is the largest component of this system, comprising 45 different subunits encoded by mitochondrial and nuclear DNA. Humans diagnosed with mutations in the gene NDUFS4, encoding a nuclear DNA-encoded subunit of CI (NADH dehydrogenase ubiquinone Fe-S protein 4), typically suffer from Leigh syndrome, a neurodegenerative disease with onset in infancy or early childhood. Mitochondria from NDUFS4 patients usually lack detectable NDUFS4 protein and show a CI stability/assembly defect. Here, we describe a recessive mouse phenotype caused by the insertion of a transposable element into Ndufs4, identified by a novel combined linkage and expression analysis. Designated Ndufs4(fky), the mutation leads to aberrant transcript splicing and absence of NDUFS4 protein in all tissues tested of homozygous mice. Physical and behavioral symptoms displayed by Ndufs4(fky/fky) mice include temporary fur loss, growth retardation, unsteady gait, and abnormal body posture when suspended by the tail. Analysis of CI in Ndufs4(fky/fky) mice using blue native PAGE revealed the presence of a faster migrating crippled complex. This crippled CI was shown to lack subunits of the "N assembly module", which contains the NADH binding site, but contained two assembly factors not present in intact CI. Metabolomic analysis of the blood by tandem mass spectrometry showed increased hydroxyacylcarnitine species, implying that the CI defect leads to an imbalanced NADH/NAD(+) ratio that inhibits mitochondrial fatty acid β-oxidation.

Resting metabolic rate and protein turnover in apparently healthy elderly Gambian men.

Body composition, resting energy expenditure (REE), and whole body protein metabolism were studied in 26 young and 28 elderly Gambian men matched for body mass index during the dry season in a rural village in The Gambia. REE was measured by indirect calorimetry (hood system) in the fasting state and after five successive meals. Rates of whole body nitrogen flux, protein synthesis, and protein breakdown were determined in the fed state from the level of isotopic enrichment of urinary ammonia over a period of 12 h after a single oral dose of [15N]glycine. Expressed in absolute value, REE was significantly lower in the elderly compared with the young group (3.21 +/- 0.07 vs. 4.04 +/- 0.07 kJ/min, P < 0.001) and when adjusted to body weight (3.29 +/- 0.05 vs. 3.96 +/- 0.05 kJ/min, P < 0.0001) and fat-free mass (FFM; 3.38 +/- 0.01 vs. 3.87 +/- 0.01 kJ/min, P < 0.0001). The rate of protein synthesis averaged 207 +/- 13 g protein/day in the elderly and 230 +/- 13 g protein/day in the young group, whereas protein breakdown averaged 184 +/- 13 g protein/day in the elderly and 203 +/- 13 g protein/day in the young group (nonsignificant). When values were adjusted for body weight or FFM, they did not reveal any difference between the two groups. It is concluded that the reduced REE adjusted for body composition observed in elderly Gambian men is not explained by a decrease in protein turnover.

Expression of the alpha 1b-adrenergic receptor and G protein subunits in mammalian cell lines using the Semliki Forest virus expression system.

Semliki Forest virus (SFV) vectors have been efficiently used for rapid high level expression of several G protein-coupled receptors. Here we describe the use of SFV vectors to express the alpha 1b-adrenergic receptor (AR) alone or in the presence of the G protein alpha q and/or beta 2 and gamma 2 subunits. Infection of baby hamster kidney (BHK) cells with recombinant SFV-alpha 1b-AR particles resulted in high specific binding activity of the alpha 1b-AR (24 pmol receptor/mg protein). Time-course studies indicated that the highest level of receptor expression was obtained 30 hours post-infection. The stimulation of BHK cells, with epinephrine led to a 5-fold increase in inositol phosphate (IP) accumulation, confirming the functional coupling of the receptor to G protein-mediated activation of phospholipase C. The SFV expression system represents a rapid and reproducible system to study the pharmacological properties and interactions of G protein coupled receptors and of G protein subunits.

The protein phosphatase 7 regulates phytochrome signaling in Arabidopsis.

The psi2 mutant of Arabidopsis displays amplification of the responses controlled by the red/far red light photoreceptors phytochrome A (phyA) and phytochrome B (phyB) but no apparent defect in blue light perception. We found that loss-of-function alleles of the protein phosphatase 7 (AtPP7) are responsible for the light hypersensitivity in psi2 demonstrating that AtPP7 controls the levels of phytochrome signaling. Plants expressing reduced levels of AtPP7 mRNA display reduced blue-light induced cryptochrome signaling but no noticeable deficiency in phytochrome signaling. Our genetic analysis suggests that phytochrome signaling is enhanced in the AtPP7 loss of function alleles, including in blue light, which masks the reduced cryptochrome signaling. AtPP7 has been found to interact both in yeast and in planta assays with nucleotide-diphosphate kinase 2 (NDPK2), a positive regulator of phytochrome signals. Analysis of ndpk2-psi2 double mutants suggests that NDPK2 plays a critical role in the AtPP7 regulation of the phytochrome pathway and identifies NDPK2 as an upstream element involved in the modulation of the salicylic acid (SA)-dependent defense pathway by light. Thus, cryptochrome- and phytochrome-specific light signals synchronously control their relative contribution to the regulation of plant development. Interestingly, PP7 and NDPK are also components of animal light signaling systems.