Whole-body protein turnover and resting energy expenditure in obese, prepubertal children.

BACKGROUND: Obesity is becoming more frequent in children; understanding the extent to which this condition affects not only carbohydrate and lipid metabolism but also protein metabolism is of paramount importance. OBJECTIVE: We evaluated the kinetics of protein metabolism in obese, prepubertal children in the static phase of obesity. DESIGN: In this cross-sectional study, 9 obese children (x +/- SE: 44+/-4 kg, 30.9+/-1.5% body fat) were compared with 8 lean (28+/-2 kg ,16.8+/-1.2% body fat), age-matched (8.5+/-0.2 y) control children. Whole-body nitrogen flux, protein synthesis, and protein breakdown were calculated postprandially over 9 h from 15N abundance in urinary ammonia by using a single oral dose of [15N]glycine; resting energy expenditure (REE) was assessed by indirect calorimetry (canopy) and body composition by multiple skinfold-thickness measurements. RESULTS: Absolute rates of protein synthesis and breakdown were significantly greater in obese children than in control children (x +/- SE: 208+/-24 compared with 137+/-14 g/d, P < 0.05, and 149+/-20 compared with 89+/-13 g/d, P < 0.05, respectively). When these variables were adjusted for fat-free mass by analysis of covariance, however, the differences between groups disappeared. There was a significant relation between protein synthesis and fat-free mass (r = 0.83, P < 0.001) as well as between protein synthesis and REE (r = 0.79, P < 0.005). CONCLUSIONS: Obesity in prepubertal children is associated with an absolute increase in whole-body protein turnover that is consistent with an absolute increase in fat-free mass, both of which contribute to explaining the greater absolute REE in obese children than in control children.

Precision of a new tool to measure visceral adipose tissue (VAT) using dual-energy X-Ray absorptiometry (DXA).

OBJECTIVE: A new tool to quantify visceral adipose tissue (VAT) over the android region of a total body dual-energy x-ray absorptiometry (DXA) scan has recently been reported. The measurement, CoreScan, is currently available on Lunar iDXA densitometers. The purpose of the study was to determine the precision of the CoreScan VAT measurement, which is critical for understanding the utility of this measure in longitudinal trials. DESIGN AND METHODS: VAT precision was characterized in both an anthropomorphic imaging phantom (measured on 10 Lunar iDXA systems) and a clinical population consisting of obese women (n = 32). RESULTS: The intrascanner precision for the VAT phantom across 9 quantities of VAT mass (0-1,800 g) ranged from 28.4 to 38.0 g. The interscanner precision ranged from 24.7 to 38.4 g. There was no statistical dependence on the quantity of VAT for either the inter- or intrascanner precision result (p = 0.670). Combining inter- and intrascanner precision yielded a total phantom precision estimate of 47.6 g for VAT mass, which corresponds to a 4.8% coefficient of variance (CV) for a 1 kg VAT mass. Our clinical population, who completed replicate total body scans with repositioning between scans, showed a precision of 56.8 g on an average VAT mass of 1110.4 g. This corresponds to a 5.1% CV. Hence, the in vivo precision result was similar to the phantom precision result. CONCLUSIONS: The study suggests that CoreScan has a relatively low precision error in both phantoms and obese women and therefore may be a useful addition to clinical trials where interventions are targeted towards changes in visceral adiposity.

Effect of a novel beta-adrenoceptor agonist (Ro 40-2148) on resting energy expenditure in obese women.

The aim of this single-blind, placebo-controlled study was to investigate the effects of the new beta-adrenergic compound Ro 40-2148 on resting energy expenditure (REE) at rest and after an oral glucose load in non-diabetic obese women before and after two weeks of treatment. After one week of placebo administration and after an overnight fast and one hour rest, REE and glucose and lipid oxidation rates were measured by indirect calorimetry (hood system) before and for 6 h after a single dose of placebo solution. A 75 g oral glucose tolerance test (OGTT) was performed during this period starting 90 min after the placebo administration. During the following two weeks, using a randomization design, six patients received Ro 40-2148 at a dose of 400 mg diluted in 100 ml water twice a day (i.e. 800 mg per day), while six others continued with the placebo administration. The same tests and measurements were repeated after two weeks, except for the treatment group which received the drug instead of the placebo. The 14-day period of drug administration did not increase REE measured in post-absorptive conditions. Similarly, there was no acute effect on REE of a 400 mg dose of Ro 40-2148. In contrast, glucose-induced thermogenesis was significantly increased after two weeks in the treatment group (means +/- s.e.m.: 3.7 +/- 1.3%, P = 0.047), while no change was observed in the placebo group (-0.8 +/- 0.7%, not significant). Since there was no significant change in the respiratory quotient, the increase in energy expenditure observed in the treatment group was due to stimulation of both lipid and glucose oxidation. The drug induced no variations in heart rate, blood pressure, axillary temperature or in plasma glucose, insulin and free fatty acid levels. In conclusion, this study shows that Ro 40-2148 activates glucose-induced thermogenesis in obese non-diabetic patients.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Correlations between trabecular bone score, measured using anteroposterior dual-energy x-ray absorptiometry acquisition, and 3-dimensional parameters of bone microarchitecture: an experimental study on human cadaver vertebrae.

Developing a novel technique for the efficient, noninvasive clinical evaluation of bone microarchitecture remains both crucial and challenging. The trabecular bone score (TBS) is a new gray-level texture measurement that is applicable to dual-energy X-ray absorptiometry (DXA) images. Significant correlations between TBS and standard 3-dimensional (3D) parameters of bone microarchitecture have been obtained using a numerical simulation approach. The main objective of this study was to empirically evaluate such correlations in anteroposterior spine DXA images. Thirty dried human cadaver vertebrae were evaluated. Micro-computed tomography acquisitions of the bone pieces were obtained at an isotropic resolution of 93μm. Standard parameters of bone microarchitecture were evaluated in a defined region within the vertebral body, excluding cortical bone. The bone pieces were measured on a Prodigy DXA system (GE Medical-Lunar, Madison, WI), using a custom-made positioning device and experimental setup. Significant correlations were detected between TBS and 3D parameters of bone microarchitecture, mostly independent of any correlation between TBS and bone mineral density (BMD). The greatest correlation was between TBS and connectivity density, with TBS explaining roughly 67.2% of the variance. Based on multivariate linear regression modeling, we have established a model to allow for the interpretation of the relationship between TBS and 3D bone microarchitecture parameters. This model indicates that TBS adds greater value and power of differentiation between samples with similar BMDs but different bone microarchitectures. It has been shown that it is possible to estimate bone microarchitecture status derived from DXA imaging using TBS.

Inhibition of Na(+)-K(+)-ATPase by digoxin and its relation with energy expenditure and nutrient oxidation rate.

This study investigates the effects of digoxin, an inhibitor of the Na+ pump (Na(+)-K(+)-ATPase), on resting metabolic rate (RMR), respiratory quotient (RQ), and nutrient oxidation rate. Twelve healthy male subjects followed a double-blind protocol design and received either 1 mg/day digoxin or a placebo 2 days before indirect calorimetry measurements. Digoxin induced a 0.22 +/- 0.07 kJ/min or 3.8 +/- 1.5% (mean +/- SE, P = 0.01) decrease in RMR and a 0.40 +/- 0.13 kJ/min (P = 0.01) decrease in fat oxidation rate, whereas carbohydrate and protein oxidation rates did not change significantly. A dose-response relationship between serum digoxin and RQ was observed. These results suggest that digoxin reduces not only RMR but also fat oxidation rate by mechanisms that remain to be elucidated. Because a linkage and an association between genes coding the Na(+)-K(+)-ATPase and the RQ have been previously observed, the present demonstration of an effect of Na(+)-K(+)-ATPase inhibition on fat oxidation rate strengthens the concept that the activity of this enzyme may play a role in body weight regulation.

Brain energy consumption induced by electrical stimulation promotes systemic glucose uptake.

BACKGROUND: Controlled transcranial stimulation of the brain is part of clinical treatment strategies in neuropsychiatric diseases such as depression, stroke, or Parkinson's disease. Manipulating brain activity by transcranial stimulation, however, inevitably influences other control centers of various neuronal and neurohormonal feedback loops and therefore may concomitantly affect systemic metabolic regulation. Because hypothalamic adenosine triphosphate-sensitive potassium channels, which function as local energy sensors, are centrally involved in the regulation of glucose homeostasis, we tested whether transcranial direct current stimulation (tDCS) causes an excitation-induced transient neuronal energy depletion and thus influences systemic glucose homeostasis and related neuroendocrine mediators.METHODS: In a crossover design testing 15 healthy male volunteers, we increased neuronal excitation by anodal tDCS versus sham and examined cerebral energy consumption with (31)phosphorus magnetic resonance spectroscopy. Systemic glucose uptake was determined by euglycemic-hyperinsulinemic glucose clamp, and neurohormonal measurements comprised the parameters of the stress systems.RESULTS: We found that anodic tDCS-induced neuronal excitation causes an energetic depletion, as quantified by (31)phosphorus magnetic resonance spectroscopy. Moreover, tDCS-induced cerebral energy consumption promotes systemic glucose tolerance in a standardized euglycemic-hyperinsulinemic glucose clamp procedure and reduces neurohormonal stress axes activity.CONCLUSIONS: Our data demonstrate that transcranial brain stimulation not only evokes alterations in local neuronal processes but also clearly influences downstream metabolic systems regulated by the brain. The beneficial effects of tDCS on metabolic features may thus qualify brain stimulation as a promising nonpharmacologic therapy option for drug-induced or comorbid metabolic disturbances in various neuropsychiatric diseases.

Rational design of 4-aryl-1,2,3-triazoles for indoleamine 2,3-dioxygenase 1 inhibition.

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. Starting from the scaffold of our previously discovered IDO1 inhibitor 4-phenyl-1,2,3-triazole, we used computational structure-based methods to design more potent ligands. This approach yielded highly efficient low molecular weight inhibitors, the most active being of nanomolar potency both in an enzymatic and in a cellular assay, while showing no cellular toxicity and a high selectivity for IDO1 over tryptophan 2,3-dioxygenase (TDO). A quantitative structure-activity relationship based on the electrostatic ligand-protein interactions in the docked binding modes and on the quantum chemically derived charges of the triazole ring demonstrated a good explanatory power for the observed activities.

MM-GBSA binding free energy decomposition and T cell receptor engineering.

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. The major determinant of T cell activation is the affinity of the TCR for the peptide-MHC complex, though kinetic parameters are also important. A study of the 2C TCR/SIYR/H-2Kb system using a binding free energy decomposition (BFED) based on the MM-GBSA approach had been performed to assess the performance of the approach on this system. The results showed that the TCR-p-MHC BFED including entropic terms provides a detailed and reliable description of the energetics of the interaction (Zoete and Michielin, 2007). Based on these results, we have developed a new approach to design sequence modifications for a TCR recognizing the human leukocyte antigen (HLA)-A2 restricted tumor epitope NY-ESO-1. NY-ESO-1 is a cancer testis antigen expressed not only in melanoma, but also on several other types of cancers. It has been observed at high frequencies in melanoma patients with unusually positive clinical outcome and, therefore, represents an interesting target for adoptive transfer with modified TCR. Sequence modifications of TCR potentially increasing the affinity for this epitope have been proposed and tested in vitro. T cells expressing some of the proposed TCR mutants showed better T cell functionality, with improved killing of peptide-loaded T2 cells and better proliferative capacity compared to the wild type TCR expressing cells. These results open the door of rational TCR design for adoptive transfer cancer therapy.

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.

Efficient total variation algorithm for fetal brain MRI reconstruction.

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy, Total Variation (TV)- based energies and more recently non-local means. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm or fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n2) and O(1/√ε), while existing techniques are in O(1/n2) and O(1/√ε). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Development of a dose-controlled multiculture cell exposure chamber for efficient delivery of airborne and engineered nanoparticles

In order to study the various health influencing parameters related to engineered nanoparticles as well as to soot emitted b diesel engines, there is an urgent need for appropriate sampling devices and methods for cell exposure studies that simulate the respiratory system and facilitate associated biological and toxicological tests. The objective of the present work was the further advancement of a Multiculture Exposure Chamber (MEC) into a dose-controlled system for efficient delivery of nanoparticles to cells. It was validated with various types of nanoparticles (diesel engine soot aggregates, engineered nanoparticles for various applications) and with state-of-the-art nanoparticle measurement instrumentation to assess the local deposition of nanoparticles on the cell cultures. The dose of nanoparticles to which cell cultures are being exposed was evaluated in the normal operation of the in vitro cell culture exposure chamber based on measurements of the size specific nanoparticle collection efficiency of a cell free device. The average efficiency in delivering nanoparticles in the MEC was approximately 82%. The nanoparticle deposition was demonstrated by Transmission Electron Microscopy (TEM). Analysis and design of the MEC employs Computational Fluid Dynamics (CFD) and true to geometry representations of nanoparticles with the aim to assess the uniformity of nanoparticle deposition among the culture wells. Final testing of the dose-controlled cell exposure system was performed by exposing A549 lung cell cultures to fluorescently labeled nanoparticles. Delivery of aerosolized nanoparticles was demonstrated by visualization of the nanoparticle fluorescence in the cell cultures following exposure. Also monitored was the potential of the aerosolized nanoparticles to generate reactive oxygen species (ROS) (e.g. free radicals and peroxides generation), thus expressing the oxidative stress of the cells which can cause extensive cellular damage or damage on DNA.

Energy gap in the aetiology of body weight gain and obesity: a challenging concept with a complex evaluation and pitfalls.

The concept of energy gap(s) is useful for understanding the consequence of a small daily, weekly, or monthly positive energy balance and the inconspicuous shift in weight gain ultimately leading to overweight and obesity. Energy gap is a dynamic concept: an initial positive energy gap incurred via an increase in energy intake (or a decrease in physical activity) is not constant, may fade out with time if the initial conditions are maintained, and depends on the 'efficiency' with which the readjustment of the energy imbalance gap occurs with time. The metabolic response to an energy imbalance gap and the magnitude of the energy gap(s) can be estimated by at least two methods, i.e. i) assessment by longitudinal overfeeding studies, imposing (by design) an initial positive energy imbalance gap; ii) retrospective assessment based on epidemiological surveys, whereby the accumulated endogenous energy storage per unit of time is calculated from the change in body weight and body composition. In order to illustrate the difficulty of accurately assessing an energy gap we have used, as an illustrative example, a recent epidemiological study which tracked changes in total energy intake (estimated by gross food availability) and body weight over 3 decades in the US, combined with total energy expenditure prediction from body weight using doubly labelled water data. At the population level, the study attempted to assess the cause of the energy gap purported to be entirely due to increased food intake. Based on an estimate of change in energy intake judged to be more reliable (i.e. in the same study population) and together with calculations of simple energetic indices, our analysis suggests that conclusions about the fundamental causes of obesity development in a population (excess intake vs. low physical activity or both) is clouded by a high level of uncertainty.

EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a 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 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, 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 in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.