|
Chargé de Recherche au CNRS (Section 16) ENS – Département de chimie Email: nicolas.cheron@ens.psl.eu Website: https://nicolascheron.wordpress.com/ |
Education and professional experience
- 2004-2008 : Education at ENS Lyon
- 2008-2012 : PhD with Paul Fleurat-Lessard, Laurence Grimaud and Laurent El Kaïm (ENS Lyon and ENSTA ParisTech)
- 2012-2015 : Postdoc with Eugene Shakhnovich (Harvard University)
- 2015-2017 : Postdoc with Damien Laage (ENS)
- 2017-now : CNRS Researcher in UMR 8640 at ENS
Research interests
- Chemical and biochemical reactivity
- CO2 reactivity
- Photosynthesis and Calvin-Benson cycle
- DNA Three-Way Junction
- Drug design
Teaching
- L1 Sciences pour un monde durable (SMD)
- L2 Cycle pluridisciplinaire d’études supérieures (CPES)
- Préparation à l’agrégation de chimie
Publications
2024 |
Structural Optimization of Azacryptands for Targeting Three-Way DNA Junctions Article de journal Angélique Pipier; Titouan Chetot; Apollonia Kalamatianou; Nicolas Martin; Maëlle Caroff; Sébastien Britton; Nicolas Chéron; Lukáš Trantírek; Anton Granzhan; David Monchaud Angewandte Chemie International Edition, n/a (n/a), p. e202409780, 2024. @article{https://doi.org/10.1002/anie.202409780, title = {Structural Optimization of Azacryptands for Targeting Three-Way DNA Junctions}, author = {Ang\'{e}lique Pipier and Titouan Chetot and Apollonia Kalamatianou and Nicolas Martin and Ma\"{e}lle Caroff and S\'{e}bastien Britton and Nicolas Ch\'{e}ron and Luk\'{a}\v{s} Trant\'{i}rek and Anton Granzhan and David Monchaud}, url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/anie.202409780}, doi = {https://doi.org/10.1002/anie.202409780}, year = {2024}, date = {2024-06-14}, journal = {Angewandte Chemie International Edition}, volume = {n/a}, number = {n/a}, pages = {e202409780}, abstract = {Transient melting of the duplex-DNA (B-DNA) during DNA transactions allows repeated sequences to fold into non-B DNA structures, including DNA junctions and G-quadruplexes. These noncanonical structures can act as impediments to DNA polymerase progression along the duplex, thereby triggering DNA damage and ultimately jeopardizing genomic stability. Their stabilization by ad hoc ligands is currently being explored as a putative anticancer strategy since it might represent an efficient way to inflict toxic DNA damage specifically to rapidly dividing cancer cells. The relevance of this strategy is only emerging for three-way DNA junctions (TWJs) and, to date, no molecule has been recognized as a reference TWJ ligand, featuring both high affinity and selectivity. Herein, we characterize such reference ligands through a combination of in vitro techniques comprising affinity and selectivity assays (competitive FRET-melting and TWJ Screen assays), functional tests (qPCR and Taq stop assays), and structural analyses (molecular dynamics and NMR investigations). We identify novel azacryptands TrisNP-amphi and TrisNP-ana as the most promising ligands, interacting with TWJs with high affinity and selectivity. These ligands represent new molecular tools to investigate the cellular roles of TWJs and explore how they can be exploited in innovative anticancer therapies.}, keywords = {}, pubstate = {published}, tppubtype = {article} } Transient melting of the duplex-DNA (B-DNA) during DNA transactions allows repeated sequences to fold into non-B DNA structures, including DNA junctions and G-quadruplexes. These noncanonical structures can act as impediments to DNA polymerase progression along the duplex, thereby triggering DNA damage and ultimately jeopardizing genomic stability. Their stabilization by ad hoc ligands is currently being explored as a putative anticancer strategy since it might represent an efficient way to inflict toxic DNA damage specifically to rapidly dividing cancer cells. The relevance of this strategy is only emerging for three-way DNA junctions (TWJs) and, to date, no molecule has been recognized as a reference TWJ ligand, featuring both high affinity and selectivity. Herein, we characterize such reference ligands through a combination of in vitro techniques comprising affinity and selectivity assays (competitive FRET-melting and TWJ Screen assays), functional tests (qPCR and Taq stop assays), and structural analyses (molecular dynamics and NMR investigations). We identify novel azacryptands TrisNP-amphi and TrisNP-ana as the most promising ligands, interacting with TWJs with high affinity and selectivity. These ligands represent new molecular tools to investigate the cellular roles of TWJs and explore how they can be exploited in innovative anticancer therapies. |
Binding Sites of Bicarbonate in Phosphoenolpyruvate Carboxylase Article de journal Nicolas Chéron Journal of Chemical Information and Modeling, 64 (8), p. 3375-3385, 2024, (PMID: 38533570). @article{doi:10.1021/acs.jcim.3c01830, title = {Binding Sites of Bicarbonate in Phosphoenolpyruvate Carboxylase}, author = {Nicolas Ch\'{e}ron}, url = {https://doi.org/10.1021/acs.jcim.3c01830}, doi = {10.1021/acs.jcim.3c01830}, year = {2024}, date = {2024-01-01}, journal = {Journal of Chemical Information and Modeling}, volume = {64}, number = {8}, pages = {3375-3385}, note = {PMID: 38533570}, keywords = {}, pubstate = {published}, tppubtype = {article} } |
Are changes in antibiotic prophylaxis recommendations responsible for an increased risk of cefazolin allergy? Article de journal Nicolas Chéron; Luc de Chaisemartin; Simon Aubert; Felix Laborier; Philippe Montravers; Catherine Neukirch; Aurélie Gouel-Chéron Anaesthesia Critical Care & Pain Medicine, 43 (2), p. 101349, 2024, ISSN: 2352-5568. @article{CHERON2024101349, title = {Are changes in antibiotic prophylaxis recommendations responsible for an increased risk of cefazolin allergy?}, author = {Nicolas Ch\'{e}ron and Luc de Chaisemartin and Simon Aubert and Felix Laborier and Philippe Montravers and Catherine Neukirch and Aur\'{e}lie Gouel-Ch\'{e}ron}, url = {https://www.sciencedirect.com/science/article/pii/S2352556824000079}, doi = {https://doi.org/10.1016/j.accpm.2024.101349}, issn = {2352-5568}, year = {2024}, date = {2024-01-01}, journal = {Anaesthesia Critical Care & Pain Medicine}, volume = {43}, number = {2}, pages = {101349}, abstract = {Background The first line of prevention of surgical site infection relies on the timely administration of antibiotic prophylaxis. First- and second-generation cephalosporins are the most recommended antibiotics in elective surgery. The incidence of cefazolin allergy has increased worldwide over the years. The sensitization mechanism of cefazolin is currently unknown, and data supporting cross-reactivity between penicillins and cephalosporins are lacking. Sensitization could occur through previous exposure either to cefazolin or to structurally related chemical agents. The objective of this study was to evaluate sensitization agents towards cefazolin. Methods The OpenBabel chemoinformatics toolbox was used to search for similarities between cefazolin and other molecules in an extensive drug database. Using the pholcodine-rocuronium similarity score as a threshold, we selected drugs with the most similar structure to that of cefazolin. Exposure to those drugs and cefazolin was assessed in a cohort of patients with skin test-proven cefazolin allergy at a specialized allergy centre via a self-administered anonymous questionnaire. Results Using the pholcodine-rocuronium similarity score as a threshold (score≥0.7), 42 molecules were found to be similar to cefazolin (all cephalosporins). Only 8 were marketed in France. None of the 14 cefazolin-allergic patients who answered the questionnaire (65% female, median age 56 years) reported exposure to any identified antibiotics. In contrast, 11 (78%) had at least one previous surgery requiring cefazolin before the index case. Conclusion Direct previous cefazolin exposure was identified in 78% of cefazolin-allergic patients. Cefazolin started to take a central place in antibiotic prophylaxis after 2010, when cefamandole usage decreased drastically. Changes in antibiotic prophylaxis over the past 14 years in France could have been the turning point for the increased incidence of cefazolin allergy.}, keywords = {}, pubstate = {published}, tppubtype = {article} } Background The first line of prevention of surgical site infection relies on the timely administration of antibiotic prophylaxis. First- and second-generation cephalosporins are the most recommended antibiotics in elective surgery. The incidence of cefazolin allergy has increased worldwide over the years. The sensitization mechanism of cefazolin is currently unknown, and data supporting cross-reactivity between penicillins and cephalosporins are lacking. Sensitization could occur through previous exposure either to cefazolin or to structurally related chemical agents. The objective of this study was to evaluate sensitization agents towards cefazolin. Methods The OpenBabel chemoinformatics toolbox was used to search for similarities between cefazolin and other molecules in an extensive drug database. Using the pholcodine-rocuronium similarity score as a threshold, we selected drugs with the most similar structure to that of cefazolin. Exposure to those drugs and cefazolin was assessed in a cohort of patients with skin test-proven cefazolin allergy at a specialized allergy centre via a self-administered anonymous questionnaire. Results Using the pholcodine-rocuronium similarity score as a threshold (score≥0.7), 42 molecules were found to be similar to cefazolin (all cephalosporins). Only 8 were marketed in France. None of the 14 cefazolin-allergic patients who answered the questionnaire (65% female, median age 56 years) reported exposure to any identified antibiotics. In contrast, 11 (78%) had at least one previous surgery requiring cefazolin before the index case. Conclusion Direct previous cefazolin exposure was identified in 78% of cefazolin-allergic patients. Cefazolin started to take a central place in antibiotic prophylaxis after 2010, when cefamandole usage decreased drastically. Changes in antibiotic prophylaxis over the past 14 years in France could have been the turning point for the increased incidence of cefazolin allergy. |
2022 |
Structure of the photosynthetic Calvin-Benson-Bassham sedoheptulose-1,7-bisphosphatase SBPase from the model microalgaChlamydomonas reinhardtii Article de journal Théo Le Moigne; Martina Santoni; Lucile Jomat; Stéphane D Lemaire; Mirko Zaffagnini; Nicolas Chéron; Julien Henri 2022. @article{Moigne_2022, title = {Structure of the photosynthetic Calvin-Benson-Bassham sedoheptulose-1,7-bisphosphatase SBPase from the model microalgaChlamydomonas reinhardtii}, author = {Th\'{e}o Le Moigne and Martina Santoni and Lucile Jomat and St\'{e}phane D Lemaire and Mirko Zaffagnini and Nicolas Ch\'{e}ron and Julien Henri}, url = {http://dx.doi.org/10.1101/2022.10.28.514230}, doi = {10.1101/2022.10.28.514230}, year = {2022}, date = {2022-01-01}, publisher = {Cold Spring Harbor Laboratory}, keywords = {}, pubstate = {published}, tppubtype = {article} } |
Click-Chemistry-Based Biomimetic Ligands Efficiently Capture G-Quadruplexes In Vitro and Help Localize Them at DNA Damage Sites in Human Cells Article de journal Francesco Rota Sperti; Baptiste Dupouy; Jérémie Mitteaux; Angélique Pipier; Marc Pirrotta; Nicolas Chéron; Ibai E Valverde; David Monchaud JACS Au, 2 (7), p. 1588-1595, 2022. @article{doi:10.1021/jacsau.2c00082, title = {Click-Chemistry-Based Biomimetic Ligands Efficiently Capture G-Quadruplexes In Vitro and Help Localize Them at DNA Damage Sites in Human Cells}, author = {Francesco Rota Sperti and Baptiste Dupouy and J\'{e}r\'{e}mie Mitteaux and Ang\'{e}lique Pipier and Marc Pirrotta and Nicolas Ch\'{e}ron and Ibai E Valverde and David Monchaud}, url = {https://doi.org/10.1021/jacsau.2c00082}, doi = {10.1021/jacsau.2c00082}, year = {2022}, date = {2022-01-01}, journal = {JACS Au}, volume = {2}, number = {7}, pages = {1588-1595}, keywords = {}, pubstate = {published}, tppubtype = {article} } |
2021 |
Dual targeting of higher-order DNA structures by azacryptands induces DNA junction-mediated DNA damage in cancer cells Article de journal Joanna Zell; Katerina Duskova; Leïla Chouh; Madeleine Bossaert; Nicolas Chéron; Anton Granzhan; Sébastien Britton; David Monchaud Nucleic Acids Research, 49 (18), p. 10275-10288, 2021, ISSN: 0305-1048. @article{10.1093/nar/gkab796, title = {Dual targeting of higher-order DNA structures by azacryptands induces DNA junction-mediated DNA damage in cancer cells}, author = {Joanna Zell and Katerina Duskova and Le\"{i}la Chouh and Madeleine Bossaert and Nicolas Ch\'{e}ron and Anton Granzhan and S\'{e}bastien Britton and David Monchaud}, url = {https://doi.org/10.1093/nar/gkab796}, doi = {10.1093/nar/gkab796}, issn = {0305-1048}, year = {2021}, date = {2021-01-01}, journal = {Nucleic Acids Research}, volume = {49}, number = {18}, pages = {10275-10288}, abstract = {DNA is intrinsically dynamic and folds transiently into alternative higher-order structures such as G-quadruplexes (G4s) and three-way DNA junctions (TWJs). G4s and TWJs can be stabilised by small molecules (ligands) that have high chemotherapeutic potential, either as standalone DNA damaging agents or combined in synthetic lethality strategies. While previous approaches have claimed to use ligands that specifically target either G4s or TWJs, we report here on a new approach in which ligands targeting both TWJs and G4s in vitro demonstrate cellular effects distinct from that of G4 ligands, and attributable to TWJ targeting. The DNA binding modes of these new, dual TWJ-/G4-ligands were studied by a panel of in vitro methods and theoretical simulations, and their cellular properties by extensive cell-based assays. We show here that cytotoxic activity of TWJ-/G4-ligands is mitigated by the DNA damage response (DDR) and DNA topoisomerase 2 (TOP2), making them different from typical G4-ligands, and implying a pivotal role of TWJs in cells. We designed and used a clickable ligand, TrisNP-α, to provide unique insights into the TWJ landscape in cells and its modulation upon co-treatments. This wealth of data was exploited to design an efficient synthetic lethality strategy combining dual ligands with clinically relevant DDR inhibitors.}, keywords = {}, pubstate = {published}, tppubtype = {article} } DNA is intrinsically dynamic and folds transiently into alternative higher-order structures such as G-quadruplexes (G4s) and three-way DNA junctions (TWJs). G4s and TWJs can be stabilised by small molecules (ligands) that have high chemotherapeutic potential, either as standalone DNA damaging agents or combined in synthetic lethality strategies. While previous approaches have claimed to use ligands that specifically target either G4s or TWJs, we report here on a new approach in which ligands targeting both TWJs and G4s in vitro demonstrate cellular effects distinct from that of G4 ligands, and attributable to TWJ targeting. The DNA binding modes of these new, dual TWJ-/G4-ligands were studied by a panel of in vitro methods and theoretical simulations, and their cellular properties by extensive cell-based assays. We show here that cytotoxic activity of TWJ-/G4-ligands is mitigated by the DNA damage response (DDR) and DNA topoisomerase 2 (TOP2), making them different from typical G4-ligands, and implying a pivotal role of TWJs in cells. We designed and used a clickable ligand, TrisNP-α, to provide unique insights into the TWJ landscape in cells and its modulation upon co-treatments. This wealth of data was exploited to design an efficient synthetic lethality strategy combining dual ligands with clinically relevant DDR inhibitors. |
2020 |
Protein Preferential Solvation in Water:Glycerol Mixtures Article de journal Nicolas Chéron; Margaux Naepels; Eva Pluhařová; Damien Laage The Journal of Physical Chemistry B, 124 (8), p. 1424-1437, 2020. @article{doi:10.1021/acs.jpcb.9b11190, title = {Protein Preferential Solvation in Water:Glycerol Mixtures}, author = {Nicolas Ch\'{e}ron and Margaux Naepels and Eva Pluha\v{r}ov\'{a} and Damien Laage}, url = {https://doi.org/10.1021/acs.jpcb.9b11190}, doi = {10.1021/acs.jpcb.9b11190}, year = {2020}, date = {2020-01-01}, journal = {The Journal of Physical Chemistry B}, volume = {124}, number = {8}, pages = {1424-1437}, keywords = {}, pubstate = {published}, tppubtype = {article} } |
2017 |
A Hybrid Knowledge-Based and Empirical Scoring Function for Protein-Ligand Interaction: SMoG2016 Article de journal T Debroise; E I Shakhnovich; N Chéron Journal of Chemical Information and Modeling, 57 (3), p. 584–593, 2017. @article{Debroise:2017, title = {A Hybrid Knowledge-Based and Empirical Scoring Function for Protein-Ligand Interaction: SMoG2016}, author = {T Debroise and E I Shakhnovich and N Ch\'{e}ron}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85025107228&doi=10.1021%2facs.jcim.6b00610&partnerID=40&md5=6b16dec8e3efb91720fd713ee75e44b5}, doi = {10.1021/acs.jcim.6b00610}, year = {2017}, date = {2017-01-01}, journal = {Journal of Chemical Information and Modeling}, volume = {57}, number = {3}, pages = {584--593}, abstract = {We present the third generation of our scoring function for the prediction of protein-ligand binding free energy. This function is now a hybrid between a knowledge-based potential and an empirical function. We constructed a diversified set of ∼1000 complexes from the PDBBinding-CN database for the training of the function, and we show that this number of complexes generates enough data to build the potential. The occurrence of 420 different types of atomic pairwise interactions is computed in up to five different ranges of distances to derive the knowledge-based part. All of the parameters were optimized, and we were able to considerably improve the accuracy of the scoring function with a Pearson correlation coefficient against experimental binding free energies of up to 0.57, which ranks our new scoring function as one of the best currently available and the second-best in terms of standard deviation (SD = 1.68 kcal/mol). The function was then further improved by inclusion of different terms taking into account repulsion and loss of entropy upon binding, and we show that it is capable of recovering native binding poses up to 80% of the time. All of the programs, tools, and protein sets are released in the Supporting Information or as open-source programs. © 2017 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } We present the third generation of our scoring function for the prediction of protein-ligand binding free energy. This function is now a hybrid between a knowledge-based potential and an empirical function. We constructed a diversified set of ∼1000 complexes from the PDBBinding-CN database for the training of the function, and we show that this number of complexes generates enough data to build the potential. The occurrence of 420 different types of atomic pairwise interactions is computed in up to five different ranges of distances to derive the knowledge-based part. All of the parameters were optimized, and we were able to considerably improve the accuracy of the scoring function with a Pearson correlation coefficient against experimental binding free energies of up to 0.57, which ranks our new scoring function as one of the best currently available and the second-best in terms of standard deviation (SD = 1.68 kcal/mol). The function was then further improved by inclusion of different terms taking into account repulsion and loss of entropy upon binding, and we show that it is capable of recovering native binding poses up to 80% of the time. All of the programs, tools, and protein sets are released in the Supporting Information or as open-source programs. © 2017 American Chemical Society. |
Chapter 13: Effects of Water and Non-Aqueous Solvents on Enzyme Activity Livre E Pluhav rová; N Chéron; D Laage 2017. @book{Pluharova:2017a, title = {Chapter 13: Effects of Water and Non-Aqueous Solvents on Enzyme Activity}, author = {E Pluha{v r}ov\'{a} and N Ch\'{e}ron and D Laage}, doi = {10.1039/9781782626831-00436}, year = {2017}, date = {2017-01-01}, volume = {2017-January}, series = {RSC Theoretical and Computational Chemistry Series}, abstract = {In this chapter, we review the available experimental data and molecular models describing the effect of different solvents, including water, on the catalytic activity of enzymes. While a popular picture suggests that water acts as a lubricant of the protein conformational motions required for catalysis, we show that this dynamical picture is not supported by recent results and we discuss alternative models. textcopyright 2017 The Royal Society of Chemistry.}, keywords = {}, pubstate = {published}, tppubtype = {book} } In this chapter, we review the available experimental data and molecular models describing the effect of different solvents, including water, on the catalytic activity of enzymes. While a popular picture suggests that water acts as a lubricant of the protein conformational motions required for catalysis, we show that this dynamical picture is not supported by recent results and we discuss alternative models. textcopyright 2017 The Royal Society of Chemistry. |
Effect of sampling on BACE-1 ligands binding free energy predictions via MM-PBSA calculations Article de journal N Chéron; E I Shakhnovich Journal of Computational Chemistry, 38 (22), p. 1941–1951, 2017. @article{Cheron:2017, title = {Effect of sampling on BACE-1 ligands binding free energy predictions via MM-PBSA calculations}, author = {N Ch\'{e}ron and E I Shakhnovich}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85020100583&doi=10.1002%2fjcc.24839&partnerID=40&md5=516fbbbbe873be6ebb343d12fe5ead0b}, doi = {10.1002/jcc.24839}, year = {2017}, date = {2017-01-01}, journal = {Journal of Computational Chemistry}, volume = {38}, number = {22}, pages = {1941--1951}, abstract = {The BACE-1 enzyme is a prime target to find a cure to Alzheimer's disease. In this article, we used the MM-PBSA approach to compute the binding free energies of 46 reported ligands to this enzyme. After showing that the most probable protonation state of the catalytic dyad is mono-protonated (on ASP32), we performed a thorough analysis of the parameters influencing the sampling of the conformational space (in total, more than 35 μs of simulations were performed). We show that ten simulations of 2 ns gives better results than one of 50 ns. We also investigated the influence of the protein force field, the water model, the periodic boundary conditions artifacts (box size), as well as the ionic strength. Amber03 with TIP3P, a minimal distance of 1.0 nm between the protein and the box edges and a ionic strength of I = 0.2 M provides the optimal correlation with experiments. Overall, when using these parameters, a Pearson correlation coefficient of R = 0.84 (R2 = 0.71) is obtained for the 46 ligands, spanning eight orders of magnitude of Kd (from 0.017 nm to 2000 μM, i.e., from −14.7 to −3.7 kcal/mol), with a ligand size from 22 to 136 atoms (from 138 to 937 g/mol). After a two-parameter fit of the binding affinities for 12 of the ligands, an error of RMSD = 1.7 kcal/mol was obtained for the remaining ligands. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.}, keywords = {}, pubstate = {published}, tppubtype = {article} } The BACE-1 enzyme is a prime target to find a cure to Alzheimer's disease. In this article, we used the MM-PBSA approach to compute the binding free energies of 46 reported ligands to this enzyme. After showing that the most probable protonation state of the catalytic dyad is mono-protonated (on ASP32), we performed a thorough analysis of the parameters influencing the sampling of the conformational space (in total, more than 35 μs of simulations were performed). We show that ten simulations of 2 ns gives better results than one of 50 ns. We also investigated the influence of the protein force field, the water model, the periodic boundary conditions artifacts (box size), as well as the ionic strength. Amber03 with TIP3P, a minimal distance of 1.0 nm between the protein and the box edges and a ionic strength of I = 0.2 M provides the optimal correlation with experiments. Overall, when using these parameters, a Pearson correlation coefficient of R = 0.84 (R2 = 0.71) is obtained for the 46 ligands, spanning eight orders of magnitude of Kd (from 0.017 nm to 2000 μM, i.e., from −14.7 to −3.7 kcal/mol), with a ligand size from 22 to 136 atoms (from 138 to 937 g/mol). After a two-parameter fit of the binding affinities for 12 of the ligands, an error of RMSD = 1.7 kcal/mol was obtained for the remaining ligands. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc. |
Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent Article de journal E Duboué-Dijon; E Pluhařová; D Domin; K Sen; A C Fogarty; N Chéron; D Laage Journal of Physical Chemistry B, 121 (29), p. 7027–7041, 2017. @article{Duboue-Dijon:2017, title = {Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent}, author = {E Dubou\'{e}-Dijon and E Pluha\v{r}ov\'{a} and D Domin and K Sen and A C Fogarty and N Ch\'{e}ron and D Laage}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85026553265&doi=10.1021%2facs.jpcb.7b03102&partnerID=40&md5=b660afe385e614dde371ae0c9cdcfd33}, doi = {10.1021/acs.jpcb.7b03102}, year = {2017}, date = {2017-01-01}, journal = {Journal of Physical Chemistry B}, volume = {121}, number = {29}, pages = {7027--7041}, abstract = {Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site. © 2017 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site. © 2017 American Chemical Society. |
2016 |
OpenGrowth: An Automated and Rational Algorithm for Finding New Protein Ligands Article de journal N Chéron; N Jasty; E I Shakhnovich Journal of Medicinal Chemistry, 59 (9), p. 4171-4188, 2016. @article{Cheron:2016a, title = {OpenGrowth: An Automated and Rational Algorithm for Finding New Protein Ligands}, author = {N Ch\'{e}ron and N Jasty and E I Shakhnovich}, doi = {10.1021/acs.jmedchem.5b00886}, year = {2016}, date = {2016-01-01}, journal = {Journal of Medicinal Chemistry}, volume = {59}, number = {9}, pages = {4171-4188}, abstract = {We present a new open-source software, called OpenGrowth, which aims to create de novo ligands by connecting small organic fragments in the active site of proteins. Molecule growth is biased to produce structures that statistically resemble drugs in an input training database. Consequently, the produced molecules have superior synthetic accessibility and pharmacokinetic properties compared with randomly grown molecules. The growth process can take into account the flexibility of the target protein and can be started from a seed to mimic R-group strategy or fragment-based drug discovery. Primary applications of the software on the HIV-1 protease allowed us to quickly identify new inhibitors with a predicted Kd as low as 18 nM. We also present a graphical user interface that allows a user to select easily the fragments to include in the growth process. OpenGrowth is released under the GNU GPL license and is available free of charge on the authors website and at http://opengrowth.sourceforge.net/. textcopyright 2015 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } We present a new open-source software, called OpenGrowth, which aims to create de novo ligands by connecting small organic fragments in the active site of proteins. Molecule growth is biased to produce structures that statistically resemble drugs in an input training database. Consequently, the produced molecules have superior synthetic accessibility and pharmacokinetic properties compared with randomly grown molecules. The growth process can take into account the flexibility of the target protein and can be started from a seed to mimic R-group strategy or fragment-based drug discovery. Primary applications of the software on the HIV-1 protease allowed us to quickly identify new inhibitors with a predicted Kd as low as 18 nM. We also present a graphical user interface that allows a user to select easily the fragments to include in the growth process. OpenGrowth is released under the GNU GPL license and is available free of charge on the authors website and at http://opengrowth.sourceforge.net/. textcopyright 2015 American Chemical Society. |
Evolutionary Dynamics of Viral Escape under Antibodies Stress: A Biophysical Model Article de journal N Chéron; A W R Serohijos; J -M Choi; E I Shakhnovich Protein Science, p. 1332-1340, 2016. @article{Cheron:2016, title = {Evolutionary Dynamics of Viral Escape under Antibodies Stress: A Biophysical Model}, author = {N Ch\'{e}ron and A W R Serohijos and J -M Choi and E I Shakhnovich}, doi = {10.1002/pro.2915}, year = {2016}, date = {2016-01-01}, journal = {Protein Science}, pages = {1332-1340}, abstract = {Viruses constantly face the selection pressure of antibodies, either from innate immune response of the host or from administered antibodies for treatment. We explore the interplay between the biophysical properties of viral proteins and the population and demographic variables in the viral escape. The demographic and population genetics aspect of the viral escape have been explored before; however one important assumption was the a priori distribution of fitness effects (DFE). Here, we relax this assumption by instead considering a realistic biophysics-based genotype-phenotype relationship for RNA viruses escaping antibodies stress. In this model the DFE is itself an evolvable property that depends on the genetic background (epistasis) and the distribution of biophysical effects of mutations, which is informed by biochemical experiments and theoretical calculations in protein engineering. We quantitatively explore in silico the viability of viral populations under antibodies pressure and derive the phase diagram that defines the fate of the virus population (extinction or escape from stress) in a range of viral mutation rates and antibodies concentrations. We find that viruses are most resistant to stress at an optimal mutation rate (OMR) determined by the competition between supply of beneficial mutation to facilitate escape from stressors and lethal mutagenesis caused by excess of destabilizing mutations. We then show the quantitative dependence of the OMR on genome length and viral burst size. We also recapitulate the experimental observation that viruses with longer genomes have smaller mutation rate per nucleotide. textcopyright 2016 The Protein Society}, keywords = {}, pubstate = {published}, tppubtype = {article} } Viruses constantly face the selection pressure of antibodies, either from innate immune response of the host or from administered antibodies for treatment. We explore the interplay between the biophysical properties of viral proteins and the population and demographic variables in the viral escape. The demographic and population genetics aspect of the viral escape have been explored before; however one important assumption was the a priori distribution of fitness effects (DFE). Here, we relax this assumption by instead considering a realistic biophysics-based genotype-phenotype relationship for RNA viruses escaping antibodies stress. In this model the DFE is itself an evolvable property that depends on the genetic background (epistasis) and the distribution of biophysical effects of mutations, which is informed by biochemical experiments and theoretical calculations in protein engineering. We quantitatively explore in silico the viability of viral populations under antibodies pressure and derive the phase diagram that defines the fate of the virus population (extinction or escape from stress) in a range of viral mutation rates and antibodies concentrations. We find that viruses are most resistant to stress at an optimal mutation rate (OMR) determined by the competition between supply of beneficial mutation to facilitate escape from stressors and lethal mutagenesis caused by excess of destabilizing mutations. We then show the quantitative dependence of the OMR on genome length and viral burst size. We also recapitulate the experimental observation that viruses with longer genomes have smaller mutation rate per nucleotide. textcopyright 2016 The Protein Society |
2015 |
Repurposing of Rutin for the Inhibition of Norovirus Replication Article de journal N Chéron; C Yu; A O Kolawole; E I Shakhnovich; C E Wobus Archives of Virology, 160 (9), p. 2353-2358, 2015. @article{Cheron:2015, title = {Repurposing of Rutin for the Inhibition of Norovirus Replication}, author = {N Ch\'{e}ron and C Yu and A O Kolawole and E I Shakhnovich and C E Wobus}, doi = {10.1007/s00705-015-2495-y}, year = {2015}, date = {2015-01-01}, journal = {Archives of Virology}, volume = {160}, number = {9}, pages = {2353-2358}, abstract = {Drug repurposing is a strategy employed to circumvent some of the bottlenecks involved in drug development, such as the cost and time needed for developing new molecular entities. Noroviruses cause recurrent epidemics and sporadic outbreaks of gastroenteritis associated with significant mortality and economic costs, but no treatment has been approved to date. Herein, a library of molecules previously used in humans was screened to find compounds with anti-noroviral activity. Antiviral testing for four selected compounds against murine norovirus infection revealed that rutin has anti-murine norovirus activity in cell-based assays. textcopyright 2015, Springer-Verlag Wien.}, keywords = {}, pubstate = {published}, tppubtype = {article} } Drug repurposing is a strategy employed to circumvent some of the bottlenecks involved in drug development, such as the cost and time needed for developing new molecular entities. Noroviruses cause recurrent epidemics and sporadic outbreaks of gastroenteritis associated with significant mortality and economic costs, but no treatment has been approved to date. Herein, a library of molecules previously used in humans was screened to find compounds with anti-noroviral activity. Antiviral testing for four selected compounds against murine norovirus infection revealed that rutin has anti-murine norovirus activity in cell-based assays. textcopyright 2015, Springer-Verlag Wien. |
2013 |
Substituent effects in ugi-smiles reactions Article de journal N Chéron; R Ramozzi; L El Kaïm; L Grimaud; P Fleurat-Lessard Journal of Physical Chemistry A, 117 (33), p. 8035–8042, 2013. @article{Cheron:2013, title = {Substituent effects in ugi-smiles reactions}, author = {N Ch\'{e}ron and R Ramozzi and L El Ka\"{i}m and L Grimaud and P Fleurat-Lessard}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-84883175431&doi=10.1021%2fjp4052227&partnerID=40&md5=ebecae65f895e7263d4627bc0de0a528}, doi = {10.1021/jp4052227}, year = {2013}, date = {2013-01-01}, journal = {Journal of Physical Chemistry A}, volume = {117}, number = {33}, pages = {8035--8042}, abstract = {In a recent communication, we described the mechanism of the well-known Ugi-type reactions with a model system (J. Org. Chem. 2012, 77, 1361-1366). Herein, focusing on the Ugi-Smiles coupling, we study the effects of each of the four reactants on the energy profile to further explain the experimental results. The variations observed with different carbonyl compounds rely on their influence on the formation of the aryl-imidate, whereas the variations on the amine preferentially affect the Smiles rearrangement. The effect of substituents on the phenol derivative is seen upon both aryl-imidate formation and the rearrangement. The effect of the isocyanide substituents is less pronounced. © 2013 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } In a recent communication, we described the mechanism of the well-known Ugi-type reactions with a model system (J. Org. Chem. 2012, 77, 1361-1366). Herein, focusing on the Ugi-Smiles coupling, we study the effects of each of the four reactants on the energy profile to further explain the experimental results. The variations observed with different carbonyl compounds rely on their influence on the formation of the aryl-imidate, whereas the variations on the amine preferentially affect the Smiles rearrangement. The effect of substituents on the phenol derivative is seen upon both aryl-imidate formation and the rearrangement. The effect of the isocyanide substituents is less pronounced. © 2013 American Chemical Society. |
2012 |
A Qualitative Failure of B3LYP for Textbook Organic Reactions Article de journal N Chéron; D Jacquemin; P Fleurat-Lessard Physical Chemistry Chemical Physics, 14 (19), p. 7170-7175, 2012. @article{Cheron:2012, title = {A Qualitative Failure of B3LYP for Textbook Organic Reactions}, author = {N Ch\'{e}ron and D Jacquemin and P {Fleurat-Lessard}}, doi = {10.1039/c2cp40438a}, year = {2012}, date = {2012-01-01}, journal = {Physical Chemistry Chemical Physics}, volume = {14}, number = {19}, pages = {7170-7175}, abstract = {Depending on the selected DFT functional, two different mechanisms are found for two organic reactions (an intramolecular nucleophilic aromatic substitution and a nucleophilic addition on a carbonyl moiety). Indeed, B3LYP predicts a concerted mechanism whereas M06-2X foresees a multistep one. Calculations at the MP4(SDQ) level proved the mechanisms to be stepwise. We studied these reactions with a large panel of exchange-correlation functionals and demonstrated that the amount of exact exchange is of first importance. For some borderline cases, the form of the functional has also an impact, e.g. the Meisenheimer $sigma$-adduct of the intramolecular nucleophilic aromatic substitution can be located with B3PW91 but not with B3LYP. These results stress the need to use recently proposed functionals to investigate chemical reactivity. textcopyright 2012 the Owner Societies.}, keywords = {}, pubstate = {published}, tppubtype = {article} } Depending on the selected DFT functional, two different mechanisms are found for two organic reactions (an intramolecular nucleophilic aromatic substitution and a nucleophilic addition on a carbonyl moiety). Indeed, B3LYP predicts a concerted mechanism whereas M06-2X foresees a multistep one. Calculations at the MP4(SDQ) level proved the mechanisms to be stepwise. We studied these reactions with a large panel of exchange-correlation functionals and demonstrated that the amount of exact exchange is of first importance. For some borderline cases, the form of the functional has also an impact, e.g. the Meisenheimer $sigma$-adduct of the intramolecular nucleophilic aromatic substitution can be located with B3PW91 but not with B3LYP. These results stress the need to use recently proposed functionals to investigate chemical reactivity. textcopyright 2012 the Owner Societies. |
A Valence Bond View of Isocyanides' Electronic Structure Article de journal R Ramozzi; N Chéron; B Bra"ida; P C Hiberty; P Fleurat-Lessard New Journal of Chemistry, 36 (5), p. 1137-1140, 2012. @article{Ramozzi:2012, title = {A Valence Bond View of Isocyanides' Electronic Structure}, author = {R Ramozzi and N Ch\'{e}ron and B Bra{"i}da and P C Hiberty and P {Fleurat-Lessard}}, doi = {10.1039/c2nj40050b}, year = {2012}, date = {2012-01-01}, journal = {New Journal of Chemistry}, volume = {36}, number = {5}, pages = {1137-1140}, abstract = {High level Valence Bond calculations support a predominantly carbenic electronic structure for isocyanides, with a secondary zwitterionic character, despite their linear geometry. This geometry results from the significant energetic stabilization due to nitrogen $pi$ lone pair donation. Results are not changed by substitution or solvation effects. textcopyright The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2012.}, keywords = {}, pubstate = {published}, tppubtype = {article} } High level Valence Bond calculations support a predominantly carbenic electronic structure for isocyanides, with a secondary zwitterionic character, despite their linear geometry. This geometry results from the significant energetic stabilization due to nitrogen $pi$ lone pair donation. Results are not changed by substitution or solvation effects. textcopyright The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2012. |
Challenging 50 Years of Established Views on Ugi Reaction: A Theoretical Approach Article de journal N Chéron; R Ramozzi; L E Ka"im; L Grimaud; P Fleurat-Lessard Journal of Organic Chemistry, 77 (3), p. 1361-1366, 2012. @article{Cheron:2012a, title = {Challenging 50 Years of Established Views on Ugi Reaction: A Theoretical Approach}, author = {N Ch\'{e}ron and R Ramozzi and L E Ka{"i}m and L Grimaud and P {Fleurat-Lessard}}, doi = {10.1021/jo2021554}, year = {2012}, date = {2012-01-01}, journal = {Journal of Organic Chemistry}, volume = {77}, number = {3}, pages = {1361-1366}, abstract = {The Ugi reaction is one of the most famous multicomponent couplings, and its efficiency is still explained by the original mechanism suggested by Ugi in the 60s. This article aims to present a thorough theoretical study of this reaction. It describes how the imine is activated and how the new stereogenic center is formed. Our calculations strongly suggest alternatives to some commonly accepted features, such as the reversibility of the intermediate steps, and temper the nature of the driving force of the reaction. textcopyright 2012 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } The Ugi reaction is one of the most famous multicomponent couplings, and its efficiency is still explained by the original mechanism suggested by Ugi in the 60s. This article aims to present a thorough theoretical study of this reaction. It describes how the imine is activated and how the new stereogenic center is formed. Our calculations strongly suggest alternatives to some commonly accepted features, such as the reversibility of the intermediate steps, and temper the nature of the driving force of the reaction. textcopyright 2012 American Chemical Society. |
2011 |
Evidences for the Key Role of Hydrogen Bonds in Nucleophilic Aromatic Substitution Reactions Article de journal N Chéron; L El Ka"im; L Grimaud; P Fleurat-Lessard Chemistry - A European Journal, 17 (52), p. 14929-14934, 2011. @article{Cheron:2011, title = {Evidences for the Key Role of Hydrogen Bonds in Nucleophilic Aromatic Substitution Reactions}, author = {N Ch\'{e}ron and L El Ka{"i}m and L Grimaud and P {Fleurat-Lessard}}, doi = {10.1002/chem.201102463}, year = {2011}, date = {2011-01-01}, journal = {Chemistry - A European Journal}, volume = {17}, number = {52}, pages = {14929-14934}, abstract = {The effect of hydrogen bonds on the fate of nucleophilic aromatic substitutions (S NAr) has been studied in silico using a density functional theory approach in the condensed phase. The importance of these hydrogen bonds can explain the "built-in solvation" model of Bunnett concerning intermolecular processes between halogenonitrobenzenes and amines. It is also demonstrated that it can explain experimental results for a multicomponent reaction (the Ugi-Smiles coupling), involving an intramolecular S NAr (the Smiles rearrangement) as the key step of the process. Modeling reveals that when an intramolecular hydrogen bond is present, it lowers the activation barrier of this step and enables the multicomponent reaction to proceed. textcopyright 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.}, keywords = {}, pubstate = {published}, tppubtype = {article} } The effect of hydrogen bonds on the fate of nucleophilic aromatic substitutions (S NAr) has been studied in silico using a density functional theory approach in the condensed phase. The importance of these hydrogen bonds can explain the "built-in solvation" model of Bunnett concerning intermolecular processes between halogenonitrobenzenes and amines. It is also demonstrated that it can explain experimental results for a multicomponent reaction (the Ugi-Smiles coupling), involving an intramolecular S NAr (the Smiles rearrangement) as the key step of the process. Modeling reveals that when an intramolecular hydrogen bond is present, it lowers the activation barrier of this step and enables the multicomponent reaction to proceed. textcopyright 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim. |
A Density Functional Theory Study of the Nef-Isocyanide Reaction: Mechanism, Influence of Parameters and Scope Article de journal N Chéron; L El Ka"im; L Grimaud; P Fleurat-Lessard Journal of Physical Chemistry A, 115 (35), p. 10106-10112, 2011. @article{Cheron:2011a, title = {A Density Functional Theory Study of the Nef-Isocyanide Reaction: Mechanism, Influence of Parameters and Scope}, author = {N Ch\'{e}ron and L El Ka{"i}m and L Grimaud and P {Fleurat-Lessard}}, doi = {10.1021/jp205909d}, year = {2011}, date = {2011-01-01}, journal = {Journal of Physical Chemistry A}, volume = {115}, number = {35}, pages = {10106-10112}, abstract = {The Nef reaction between isocyanides and acyl chlorides is studied at the M06-2X/6-311+G(d,p) level of theory in toluene. After proving that the reaction follows a concerted mechanism instead of an addition-elimination path, we study the influences of the solvent, the isocyanide, the acyl moiety and the leaving group on the energy profile of the reaction. The calculated data can be rationalized with the pKa of the leaving group, or more generally with the population of the oxygen lone pairs of the acyl moiety. textcopyright 2011 American Chemical Society.}, keywords = {}, pubstate = {published}, tppubtype = {article} } The Nef reaction between isocyanides and acyl chlorides is studied at the M06-2X/6-311+G(d,p) level of theory in toluene. After proving that the reaction follows a concerted mechanism instead of an addition-elimination path, we study the influences of the solvent, the isocyanide, the acyl moiety and the leaving group on the energy profile of the reaction. The calculated data can be rationalized with the pKa of the leaving group, or more generally with the population of the oxygen lone pairs of the acyl moiety. textcopyright 2011 American Chemical Society. |
2009 |
How Does Microhydration Impact on Structure, Spectroscopy and Formation of Disulfide Radical Anions? An Ab Initio Investigation on Dimethyldisulfide Article de journal I Geronimo; N Chéron; P Fleurat-Lessard; E Dumont Chemical Physics Letters, 481 (4-6), p. 173-179, 2009. @article{Geronimo:2009, title = {How Does Microhydration Impact on Structure, Spectroscopy and Formation of Disulfide Radical Anions? An Ab Initio Investigation on Dimethyldisulfide}, author = {I Geronimo and N Ch\'{e}ron and P {Fleurat-Lessard} and E Dumont}, doi = {10.1016/j.cplett.2009.09.079}, year = {2009}, date = {2009-01-01}, journal = {Chemical Physics Letters}, volume = {481}, number = {4-6}, pages = {173-179}, abstract = {The effect of microhydration on the one-electron uptake of a prototypical disulfide is examined at the MP2/DZP++ level of theory. After systematic generation of low-lying conformers of dimethyldisulfide-water clusters (DMDS:nH2O}, keywords = {}, pubstate = {published}, tppubtype = {article} } The effect of microhydration on the one-electron uptake of a prototypical disulfide is examined at the MP2/DZP++ level of theory. After systematic generation of low-lying conformers of dimethyldisulfide-water clusters (DMDS:nH2O |