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.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Spine bone texture assessed by trabecular bone score (TBS) predicts osteoporotic fractures in men: The Manitoba Bone Density Program.

INTRODUCTION: One quarter of osteoporotic fractures occur in men. TBS, a gray-level measurement derived from lumbar spine DXA image texture, is related to microarchitecture and fracture risk independently of BMD. Previous studies reported the ability of spine TBS to predict osteoporotic fractures in women. Our aim was to evaluate the ability of TBS to predict clinical osteoporotic fractures in men. METHODS: 3620 men aged ≥50 (mean 67.6years) at the time of baseline DXA (femoral neck, spine) were identified from a database (Province of Manitoba, Canada). Health service records were assessed for the presence of non-traumatic osteoporotic fracture after BMD testing. Lumbar spine TBS was derived from spine DXA blinded to clinical parameters and outcomes. We used Cox proportional hazard regression to analyze time to first fracture adjusted for clinical risk factors (FRAX without BMD), osteoporosis treatment and BMD (hip or spine). RESULTS: Mean followup was 4.5years. 183 (5.1%) men sustain major osteoporotic fractures (MOF), 91 (2.5%) clinical vertebral fractures (CVF), and 46 (1.3%) hip fractures (HF). Correlation between spine BMD and spine TBS was modest (r=0.31), less than correlation between spine and hip BMD (r=0.63). Significantly lower spine TBS were found in fracture versus non-fracture men for MOF (p<0.001), HF (p<0.001) and CVF (p=0.003). Area under the receiver operating characteristic curve (AUC) for incident fracture discrimination with TBS was significantly better than chance (MOF AUC=0.59, p<0.001; HF AUC=0.67, p<0.001; CVF AUC=0.57, p=0.032). TBS predicted MOF and HF (but not CVF) in models adjusted for FRAX without BMD and osteoporosis treatment. TBS remained a predictor of HF (but not MOF) after further adjustment for hip BMD or spine BMD. CONCLUSION: We observed that spine TBS predicted MOF and HF independently of the clinical FRAX score, HF independently of FRAX and BMD in men. Studies with more incident fractures are needed to confirm these findings.

Melanin-specific life-history strategies.

Abstract The maintenance of genetic variation is a long-standing issue because the adaptive value of life-history strategies associated with each genetic variant is usually unknown. However, evidence for the coexistence of alternative evolutionary fixed strategies at the population level remains scarce. Because in the tawny owl (Strix aluco) heritable melanin-based coloration shows different physiological and behavioral norms of reaction, we investigated whether coloration is associated with investment in maintenance and reproduction. Light melanic owls had lower adult survival compared to dark melanic conspecifics, and color variation was related to the trade-off between offspring number and quality. When we experimentally enlarged brood size, light melanic males produced more fledglings but in poorer condition, and they were less often recruited in the local breeding population than those of darker melanic conspecifics. Our results also suggest that dark melanic males allocate a constant effort to raise their brood independently of environmental conditions, whereas lighter melanic males finely adjust reproductive effort in relation to changes in environmental conditions. Color traits can therefore be associated with life-history strategies, and stochastic environmental perturbation can temporarily favor one phenotype over others. The existence of fixed strategies implies that some phenotypes can sometimes display a "maladapted" strategy. Long-term population monitoring is therefore vital for a full understanding of how different genotypes deal with trade-offs.

Fine structural variations of alphabeta TCRs selected by vaccination with natural versus altered self-antigen in melanoma patients

Immunotherapy of cancer is often performed with altered "analog" peptide Ags optimized for HLA class I binding, resulting in enhanced immunogenicity, but the induced T cell responses require further evaluation. Recently, we demonstrated fine specificity differences and enhanced recognition of naturally presented Ag by T cells after vaccination with natural Melan-A/MART-1 peptide, as compared with analog peptide. In this study, we compared the TCR primary structures of 1489 HLA-A*0201/Melan-A26-35-specific CD8 T cells derived from both cohorts of patients. Although a strong preference for TRAV12-2 segment usage was present in nearly all patients, usage of particular TRAJ gene segments and CDR3 composition differed slightly after vaccination with natural vs analog peptide. Moreover, TCR β-chain repertoires were broader after natural than analog peptide vaccination. In all patients, we observed a marked conservation of the CDR3β amino acid composition with recurrent sequences centered on a glycyl-leucyl/valyl/alanyl-glycyl motif. In contrast to viral-specific TCR repertoires, such "public" motifs were primarily expressed by nondominant T cell clonotypes, which contrasted with "private" CDR3β signatures frequently found in T cell clonotypes that dominated repertoires of individual patients. Interestingly, no differences in functional avidity were observed between public and private T cell clonotypes. Collectively, our data indicate that T cell repertoires generated against natural or analog Melan-A peptide exhibited slightly distinct but otherwise overlapping and structurally conserved TCR features, suggesting that the differences in binding affinity/avidity of TCRs toward pMHC observed in the two cohorts of patients are caused by subtle structural TCR variations.

Mineral zoning and geochemistry of epithermal polymetallic Zn-Pb-Ag-Cu-Bi mineralization at Cerro de Pasco, Peru

The large Cerro de Pasco Cordilleran base metal deposit in central Peru is located on the eastern margin of a middle Miocene diatreme-dome complex and comprises two mineralization stages. The first stage consists of a large pyrite-quartz body replacing Lower Mesozoic Pucara carbonate rocks and, to a lesser extent, diatreme breccia. This body is composed of pyrite with pyrrhotite inclusions, quartz, and black and red chalcedony (containing hypogene hematite). At the contact with the pyrite-quartz body, the diatreme breccia is altered to pyrite-quartz-sericite-pyrite. This body was, in part, replaced by pipelike pyrrhotite bodies zoned outward to carbonate-replacement Zn-Pb ores hearing Fe-rich sphalerite (up to 24 mol % Fes). The second mineralization stage is partly superimposed on the first and consists of zoned east-west-trending Cu-Ag-(Au-Zn-Pb) enargite-pyrite veins hosted in the diatreme breccia in the western part of the deposit and well-zoned Zn-Pb-(Bi-Ag-Cu) carbonate-replacement orebodies; in both cases, sphalerite is Fe poor and the inner parts of the orebodies show typically advanced argillic alteration assemblages, including aluminum phosphate Sulfate (APS) minerals. The zoned enargite-pyrite veins display mineral zoning, from a core of enargite-pyrite +/- alunite with traces of Au, through an intermediate zone of tennantite, chalcopyrite, and Bi minerals to a poorly developed Outer zone hearing sphalerite-galena +/- kaolinite. The carbonate-hosted replacement ores are controlled along N 35 degrees E, N 90 degrees E, N 120 degrees E, and N 170 degrees E faults. They form well-zoned upward-flaring pipelike orebodies with a core of famatinite-pyrite and alunite, an intermediate zone with tetrahedrite-pyrite, chalcopyrite, matildite, cuprobismutite, emplectite, and other Bi minerals accompanied by APS minerals, kaolinite, and dickite, and an outer zone composed of Fe-poor sphalerite (in the range of 0.05-3.5 mol % Fes) and galena. The outermost zone consists of hematite, magnetite, and Fe-Mn-Zn-Ca-Mg carbonates. Most of the second-stage carbonate-replacement orebodies plunge between 25 degrees and 60 degrees to the west, suggesting that the hydrothermal fluids ascended from deeper levels and that no lateral feeding from the veins to the carbonate-replacement orebodies took place. In the Venencocha and Santa Rosa areas, located 2.5 km northwest of the Cerro de Pasco open pit and in the southern part of the deposit, respectively, advanced argillic altered dacitic domes and oxidized veins with advanced argillic alteration halos occur. The latter veins are possibly the oxidized equivalent of the second-stage enargite-pyrite veins located in the western part of the deposit. The alteration assemblage quartz-muscovite-pyrite associated with the pyrite-quartz body suggests that the first stage precipitated at slightly, acidic fin. The sulfide mineral assemblages define an evolutionary path close to the pyrite-pyrrhotite boundary and are characteristic of low-sulfidation states; they suggest that the oxidizing slightly acidic hydrothermal fluid was buffered by phyllite, shale, and carbonate host rock. However, the presence in the pyrite-quartz body of hematite within quartz suggests that, locally, the fluids were less buffered by the host rock. The mineral assemblages of the second mineralization stage are characteristic of high- to intermediate-sulfidation states. High-sulfidation states and oxidizing conditions were achieved and maintained in the cores of the second-stage orebodies, even in those replacing carbonate rocks. The observation that, in places, second-stage mineral assemblages are found in the inner and outer zones is explained in terms of the hydrothermal fluid advancing and waning. Microthermometric data from fluid inclusions in quartz indicate that the different ores of the first mineralization stage formed at similar temperatures and moderate salinities (200 degrees-275 degrees C and 0.2-6.8 wt % NaCl equiv in the pyrite-quartz body; 192 degrees-250 degrees C and 1.1-4.3 wt % NaCl equiv in the pyrrhotite bodies; and 183 degrees-212 degrees C and 3.2-4.0 wt % NaCl equiv in the Zn-Pb ores). These values are similar to those obtained for fluid inclusions in quartz and sphalerite from the second-stage ores (187 degrees-293 degrees C and 0.2-5.2 wt % NaCl equiv in the enargite-pyrite veins: 178 degrees-265 degrees C and 0.2-7.5 wt % NaCl equiv in quartz of carbonate-replacement orebodies; 168 degrees-999 degrees C and 3-11.8 wt % NaCl equiv in sphalerite of carbonate-replacement orebodies; and 245 degrees-261 degrees C and 3.2-7.7 wt % NaCl equiv in quartz from Venencocha). Oxygen and hydrogen isotope compositions oil kaolinite from carbonate-replacement orebodies (delta(18)O = 5.3-11.5%o, delta D = -82 to -114%o) and on alunite from the Venencocha and Santa Rosa areas (delta(18)O = 1.9-6.9%o, delta D = -56 to -73%o). Oxygen isotope compositions of quartz from the first and second stages have 6180 values from 9.1 to 1.7.8 per mil. Calculated fluids in equilibrium with kaolinite have delta(18)O values of 2.0 to 8.2 and delta D values of -69 to -97 per mil; values in equilibrium with alunite are -1.4 to -6.4 and -62 to -79 per mil. Sulfur isotope compositions of sulfides from both stages have a narrow range of delta(34)S values, between -3.7 and +4.2 per mil; values for sulfates from the second stage are between 4.2 and 31.2 per mil. These results define two mixing trends for the ore-forming fluids. The first trend reflects mixing between a moderately saline (similar to 10 wt % NaCl equiv) magmatic end member that had degassed (as indicated by the low delta D values) and meteoric water. The second mixing indicates condensation of magmatic vapor with HCl and SO(2) into meteoric water, which formed alunite. The hydrothermal system at Cerro de Pasco was emplaced at a shallow depth (similar to 500 m) in the epithermal and upper part of a porphyry environment. The similar temperatures and salinities obtained for the first stage and second stages, together with the stable isotope data, indicate that both stages are linked and represent successive stages of epithermal polymetallic mineralization in the upper part of a porphyry system.

Prepregnancy body mass index and resting metabolic rate during pregnancy.

AIM: The resting metabolic rate (RMR) varies among pregnant women. The factors responsible for this variability are unknown. This study aimed to assess the influence of the prepregnancy body mass index (BMI) on the RMR during late pregnancy. METHODS: RMR, height, weight, and total (TEE) and activity (AEE) energy expenditures were measured in 46 healthy women aged 31 ± 5 years (mean ± SD) with low (<19.8), normal (19.8-26.0), and high (>26.0) prepregnancy BMI at 38.2 ± 1.5 weeks of gestation (t(gest)) and 40 ± 7 weeks postpartum (t(post)) (n = 27). RESULTS: The mean t(gest) RMR for the low-, normal-, and high-BMI groups was 1,373, 1,807, and 2,191 kcal/day, respectively (p = 0.001). The overall mean t(gest) RMR was 316 ± 183 kcal/day (21%), higher than the overall mean t(post) value and this difference was correlated with gestational weight gain (r = 0.78, p < 0.001). The scaled metabolic rate by allometry (RMR/kilograms⁰·⁷³) was similar in the low-, normal-, and high-BMI groups, respectively (p = 0.45). Changes in t(gest) TEE closely paralleled changes in t(gest) RMR (r = 0.84, p < 0.001). AEE was similar among the BMI groups. CONCLUSION: The RMR is significantly increased in the third trimester of pregnancy. The absolute gestational RMR is higher in women with high prepregnancy BMI due to increased body weight. The scaled metabolic rate (RMR/kilograms⁰·⁷³) is similar among the BMI groups of pregnant women.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Intensified chemotherapy with ACVBP plus rituximab versus standard CHOP plus rituximab for the treatment of diffuse large B-cell lymphoma (LNH03-2B): an open-label randomised phase 3 trial.

BACKGROUND: The outcome of diffuse large B-cell lymphoma has been substantially improved by the addition of the anti-CD20 monoclonal antibody rituximab to chemotherapy regimens. We aimed to assess, in patients aged 18-59 years, the potential survival benefit provided by a dose-intensive immunochemotherapy regimen plus rituximab compared with standard treatment plus rituximab. METHODS: We did an open-label randomised trial comparing dose-intensive rituximab, doxorubicin, cyclophosphamide, vindesine, bleomycin, and prednisone (R-ACVBP) with subsequent consolidation versus standard rituximab, doxorubicin, cyclophosphamide, vincristine, and prednisone (R-CHOP). Random assignment was done with a computer-assisted randomisation-allocation sequence with a block size of four. Patients were aged 18-59 years with untreated diffuse large B-cell lymphoma and an age-adjusted international prognostic index equal to 1. Our primary endpoint was event-free survival. Our analyses of efficacy and safety were of the intention-to-treat population. This study is registered with ClinicalTrials.gov, number NCT00140595. FINDINGS: One patient withdrew consent before treatment and 54 did not complete treatment. After a median follow-up of 44 months, our 3-year estimate of event-free survival was 81% (95% CI 75-86) in the R-ACVBP group and 67% (59-73) in the R-CHOP group (hazard ratio [HR] 0·56, 95% CI 0·38-0·83; p=0·0035). 3-year estimates of progression-free survival (87% [95% CI, 81-91] vs 73% [66-79]; HR 0·48 [0·30-0·76]; p=0·0015) and overall survival (92% [87-95] vs 84% [77-89]; HR 0·44 [0·28-0·81]; p=0·0071) were also increased in the R-ACVBP group. 82 (42%) of 196 patients in the R-ACVBP group experienced a serious adverse event compared with 28 (15%) of 183 in the R-CHOP group. Grade 3-4 haematological toxic effects were more common in the R-ACVBP group, with a higher proportion of patients experiencing a febrile neutropenic episode (38% [75 of 196] vs 9% [16 of 183]). INTERPRETATION: Compared with standard R-CHOP, intensified immunochemotherapy with R-ACVBP significantly improves survival of patients aged 18-59 years with diffuse large B-cell lymphoma with low-intermediate risk according to the International Prognostic Index. Haematological toxic effects of the intensive regimen were raised but manageable. FUNDING: Groupe d'Etudes des Lymphomes de l'Adulte and Amgen.

Traité sur l'art de chasser avec le chien courant / par le Cte Boisrot de Lacour ; illustrations de Carle Vernet

Collection : Les maîtres de la vénerie ; VII

Développer les habiletés sociales par le jeu

L'entraînement des habiletés sociales est une intervention utile pour améliorer le fonctionnement interpersonnel de personnes qui ont un diagnostic de schizophrénie. Les nouvelles connaissances sur la cognition sociale doivent être intégrées dans l'entraînement. Cet article présente des ajouts en ce sens apportés au jeu Compétence et les résultats préliminaires d'une étude pilote d'évaluation de cette nouvelle version.