Étude de la réaction de déamidation dans l'enzyme triosephosphate isomérase au moyen d'outils de calculs en chimie

Deamidation is the posttranslational modification of asparagine (Asn) and glutamine (Glu) residues, which is observed in several proteins and peptides. It has been shown that deamidation limits the lifetime of these macromolecules. In this work, deamidation of asparagine in small peptides and in the enzyme triosephosphate isomerase has been modeled. Deamidation in mammalian triosephosphate isomerase has been observed at two distinct deamidation sites: Asn15 and Asn71. Asn71 deamidates faster than Asn15 and slower than a small peptide. It has been suggested that, deamidation at Asn15 occurs with the influence of deamidated Asn71. In order to explain these experimental findings, microsecond long classical molecular dynamics simulations and free energy calculations using quantum mechanics/molecular mechanics tools combined with umbrella sampling technique have been performed. The sequential deamidation in triosephosphate isomerase has been shown to be related with both global and local effects. These results bring a new perspective to the impact of the high-order structure on deamidation rate. The most plausible route of this reaction was also determined. The pKa shift of backbone amide of the residue adjacent to asparagine has been found to be one of the most crucial factor determining the rate of deamidation. Considering the importance of pKa shifts in protein environment, a computational protocol was suggested in order to obtain accurate and fast pKa predictions. This protocol was applied to small organic molecules, and it has been shown to be applicable to studies concerning aminoacid pKa predictionsLa déamidation est la modification post-traductionnelle de l'asparagine (Asn) et de la glutamine (Glu). Elle est communèment observée dans les peptides et les protéines. Il a été démontré que la déamidation limite la durée de vie de ces macromolécules. Dans ce travail, la déamidation de l'asparagine dans des petits peptides et dans l'enzyme triosephosphate isomérase a été modélisée. La déamidation dans la triosephosphate isomérase de mammifères a été observée sur deux sites distincts: Asn15 et Asn71. Asn71 a une vitesse de déamidation plus élevée que Asn15 et moins grande que pour un petit peptide. Il a été suggéré que la déamidation de Asn15 se produit sous l'influence de la déamidation de Asn71. Pour expliquer ces résultats expérimentaux, des simulations de dynamiques moléculaires classiques à l'échelle de la microseconde et des calculs d'énergie libre, de type umbrella sampling, à l'aide de méthodes combinées mécanique quantique/mécanique moléculaire ont été réalisés. Nous montrons que la déamidation séquentielle dans la triosephosphate isomérase est due à la fois à des effets locaux et globaux. Ces résultats apporte une nouvelle perspective sur l'impact de l'ordre structurel sur la vitesse de déamidation Nous avons également déterminé la voie la plus plausible de cette reaction ainsi que l'influence de la variation du pKa, dans la chaîne principale, de la partie amide du résidu adjacent de l'asparagine sur la vitesse de déamidation. En regard de l'importance des variations de pKa dans l'environnement protéique, nous avons élaboré un protocole informatique permettant d'évaluer de manière rapide et précise des pKa . Ce protocole a été appliqué à des petites molécules organiques et nous avons montré qu'il était également applicable à des études relatives à la prédiction de pKa dans les protéine

Investigation of the deamidation reaction in the enzyme triosephosphate isomerase by means of computational chemistry tools

La déamidation est la modification post-traductionnelle de l'asparagine (Asn) et de la glutamine (Glu). Elle est communèment observée dans les peptides et les protéines. Il a été démontré que la déamidation limite la durée de vie de ces macromolécules. Dans ce travail, la déamidation de l'asparagine dans des petits peptides et dans l'enzyme triosephosphate isomérase a été modélisée. La déamidation dans la triosephosphate isomérase de mammifères a été observée sur deux sites distincts: Asn15 et Asn71. Asn71 a une vitesse de déamidation plus élevée que Asn15 et moins grande que pour un petit peptide. Il a été suggéré que la déamidation de Asn15 se produit sous l'influence de la déamidation de Asn71. Pour expliquer ces résultats expérimentaux, des simulations de dynamiques moléculaires classiques à l'échelle de la microseconde et des calculs d'énergie libre, de type umbrella sampling, à l'aide de méthodes combinées mécanique quantique/mécanique moléculaire ont été réalisés. Nous montrons que la déamidation séquentielle dans la triosephosphate isomérase est due à la fois à des effets locaux et globaux. Ces résultats apporte une nouvelle perspective sur l'impact de l'ordre structurel sur la vitesse de déamidation Nous avons également déterminé la voie la plus plausible de cette reaction ainsi que l'influence de la variation du pKa, dans la chaîne principale, de la partie amide du résidu adjacent de l'asparagine sur la vitesse de déamidation. En regard de l'importance des variations de pKa dans l'environnement protéique, nous avons élaboré un protocole informatique permettant d'évaluer de manière rapide et précise des pKa . Ce protocole a été appliqué à des petites molécules organiques et nous avons montré qu'il était également applicable à des études relatives à la prédiction de pKa dans les protéinesDeamidation is the posttranslational modification of asparagine (Asn) and glutamine (Glu) residues, which is observed in several proteins and peptides. It has been shown that deamidation limits the lifetime of these macromolecules. In this work, deamidation of asparagine in small peptides and in the enzyme triosephosphate isomerase has been modeled. Deamidation in mammalian triosephosphate isomerase has been observed at two distinct deamidation sites: Asn15 and Asn71. Asn71 deamidates faster than Asn15 and slower than a small peptide. It has been suggested that, deamidation at Asn15 occurs with the influence of deamidated Asn71. In order to explain these experimental findings, microsecond long classical molecular dynamics simulations and free energy calculations using quantum mechanics/molecular mechanics tools combined with umbrella sampling technique have been performed. The sequential deamidation in triosephosphate isomerase has been shown to be related with both global and local effects. These results bring a new perspective to the impact of the high-order structure on deamidation rate. The most plausible route of this reaction was also determined. The pKa shift of backbone amide of the residue adjacent to asparagine has been found to be one of the most crucial factor determining the rate of deamidation. Considering the importance of pKa shifts in protein environment, a computational protocol was suggested in order to obtain accurate and fast pKa predictions. This protocol was applied to small organic molecules, and it has been shown to be applicable to studies concerning aminoacid pKa prediction

Proton relay network in P450cam formed upon docking of putidaredoxin

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key OO cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular mechanics computations, both on the well-studied FeOO state and on Cpd(0) (compound 0). Our simulations include (all camphor-bound): (a) WT (wild type), FeOO state. (b) WT, Cpd(0). (c) Pdx (Putidaredoxin, redox partner of P450)-docked-WT, FeOO state. (d) Pdx-docked WT, Cpd(0). (e) Pdx-docked T252A mutant, Cpd(0). Among our key findings: (a) Effect of Pdx docking appears to go far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (b) Specific proton relay networks we identify are: FeOO(H)MIDLINE HORIZONTAL ELLIPSIST252MIDLINE HORIZONTAL ELLIPSISnH2OMIDLINE HORIZONTAL ELLIPSISD251 in WT; FeOO(H)MIDLINE HORIZONTAL ELLIPSISnH2OMIDLINE HORIZONTAL ELLIPSISD251 in T252A mutant; both occur with Pdx docking. (c) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (d) In WT, T252 is in the proton relay network. (e) Positioning of camphor appears significant: when camphor is part of H-bonding network, second protonation appears to be facilitated

Proton Relay Network in P450cam is Formed Upon Docking of its Redox Partner Putidaredoxin (Pdx)

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O-O cleavage step remain elusive to date; effects observed in various enzyme mutants remain unexplained. We have carried out extended (up to 1000 ns) MM-MD and follow-on QM/MM computations, both on the well-studied FeOO state and, for the first time, on Cpd(0). Our simulations include (all camphor-bound) : (1) WT (wild type), FeOO state. (2) WT, Cpd(0). (3) Pdx-docked-WT, FeOO state. (4) Pdx-docked WT, Cpd(0). (5) Pdx-docked T252A mutant, Cpd(0). Among our key findings, for the first time to our knowledge: (1) Effect of Pdx docking goes far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (2) The specific proton relay networks we identify are FeOO(H)---T252---nH2O---D251 in WT and FeOO(H)---nH2O---D251 in T252A mutant; both occur with Pdx docking. (3) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (4) T252 is in the proton relay network. (5) Positioning of camphor is crucial: when camphor is part of H-bonding network, coupling is facilitated.<br /

Molecular Modelling Reveals Eight Novel Druggable Binding Sites in SARS-CoV-2’s Spike Protein

Spike glycoprotein (S), one of the signature proteins of the SARS-CoV-2, initiates the membrane fusion and virus entry to the host cell. The S protein’s key role in virus viability makes it an attractive candidate for drug design studies. Besides the recent structural characterization of the S protein, information fundamental to drug design such as possible binding sites or molecular fragments with high affinity towards the protein is unknown. We explored the druggability of this protein, focusing on its S1 and S2 domains. We performed virtual screening studies on both closed and open forms of the protein, using both cryo-EM structures and geometries obtained from molecular dynamic simulations. We targeted 20 distinct ligand binding centres with a set of about 9,000 molecules. Our docking calculations followed by molecular mechanics-based refinement of ligand/protein complexes led us to detect eight binding sites that were so far undocumented. By further focusing on a subset of approximately 1,000 approved and marketed drugs, we aimed at suggesting a new direction for drug repurposing strategies that were not considered so far. Within this approved set, our best hits include a number of antibacterial and antiviral drugs (e.g., Streptomycin, Nelfinavir), which were not yet investigated clinically in treating COVID-19. We also identified some molecules (e.g., folic acid, Famotidine) that were already suggested to be effective towards SARS-CoV-2, yet without molecular explanation. Our results also indicate a great affinity of SARS-CoV-2’s S protein towards nucleoside analogues, either approved or experimental.</p

Investigation of the deamidation reaction in the enzyme triosephosphate isomerase by means of computational chemistry tools

La déamidation est la modification post-traductionnelle de l'asparagine (Asn) et de la glutamine (Glu). Elle est communèment observée dans les peptides et les protéines. Il a été démontré que la déamidation limite la durée de vie de ces macromolécules. Dans ce travail, la déamidation de l'asparagine dans des petits peptides et dans l'enzyme triosephosphate isomérase a été modélisée. La déamidation dans la triosephosphate isomérase de mammifères a été observée sur deux sites distincts: Asn15 et Asn71. Asn71 a une vitesse de déamidation plus élevée que Asn15 et moins grande que pour un petit peptide. Il a été suggéré que la déamidation de Asn15 se produit sous l'influence de la déamidation de Asn71. Pour expliquer ces résultats expérimentaux, des simulations de dynamiques moléculaires classiques à l'échelle de la microseconde et des calculs d'énergie libre, de type umbrella sampling, à l'aide de méthodes combinées mécanique quantique/mécanique moléculaire ont été réalisés. Nous montrons que la déamidation séquentielle dans la triosephosphate isomérase est due à la fois à des effets locaux et globaux. Ces résultats apporte une nouvelle perspective sur l'impact de l'ordre structurel sur la vitesse de déamidation Nous avons également déterminé la voie la plus plausible de cette reaction ainsi que l'influence de la variation du pKa, dans la chaîne principale, de la partie amide du résidu adjacent de l'asparagine sur la vitesse de déamidation. En regard de l'importance des variations de pKa dans l'environnement protéique, nous avons élaboré un protocole informatique permettant d'évaluer de manière rapide et précise des pKa . Ce protocole a été appliqué à des petites molécules organiques et nous avons montré qu'il était également applicable à des études relatives à la prédiction de pKa dans les protéinesDeamidation is the posttranslational modification of asparagine (Asn) and glutamine (Glu) residues, which is observed in several proteins and peptides. It has been shown that deamidation limits the lifetime of these macromolecules. In this work, deamidation of asparagine in small peptides and in the enzyme triosephosphate isomerase has been modeled. Deamidation in mammalian triosephosphate isomerase has been observed at two distinct deamidation sites: Asn15 and Asn71. Asn71 deamidates faster than Asn15 and slower than a small peptide. It has been suggested that, deamidation at Asn15 occurs with the influence of deamidated Asn71. In order to explain these experimental findings, microsecond long classical molecular dynamics simulations and free energy calculations using quantum mechanics/molecular mechanics tools combined with umbrella sampling technique have been performed. The sequential deamidation in triosephosphate isomerase has been shown to be related with both global and local effects. These results bring a new perspective to the impact of the high-order structure on deamidation rate. The most plausible route of this reaction was also determined. The pKa shift of backbone amide of the residue adjacent to asparagine has been found to be one of the most crucial factor determining the rate of deamidation. Considering the importance of pKa shifts in protein environment, a computational protocol was suggested in order to obtain accurate and fast pKa predictions. This protocol was applied to small organic molecules, and it has been shown to be applicable to studies concerning aminoacid pKa predictionsNANCY-INPL-Bib. électronique (545479901) / SudocSudocFranceF

Why Does Asn71 Deamidate Faster Than Asn15 in the Enzyme Triosephosphate Isomerase? Answers from Microsecond Molecular Dynamics Simulation and QM/MM Free Energy Calculations

Deamidation is the uncatalyzed process by which asparagine or glutamine can be transformed into aspartic acid or glutamic acid, respectively. In its active homodimeric form, mammalian triosephosphate isomerase (TPI) contains two deamidation sites per monomer. Experimental evidence shows that the primary deamidation site (Asn71-Gly72) deamidates faster than the secondary deamidation site (Asn15-Gly16). To evaluate the factors controlling the rates of these two deamidation sites in TPI, we have performed graphics processing unit-enabled microsecond long molecular dynamics simulations of rabbit TPI. The kinetics of asparagine dipeptide and two deamidation sites in mammalian TPI are also investigated using quantum mechanical/molecular mechanical tools with the umbrella sampling technique. Analysis of the simulations has been performed using independent global and local descriptors that can influence the deamidation rates: desolvation effects, backbone acidity, and side chain conformations. Our findings show that all the descriptors add up to favor the primary deamidation site over the secondary one in mammalian TPI: Asn71 deamidates faster because it is more solvent accessible, the adjacent glycine NH backbone acidity is enhanced, and the Asn side chain has a preferential near attack conformation. The crucial impact of the backbone amide acidity of the adjacent glycine on the deamidation rate is shown by kinetic analysis. Our findings also shed light on the effect of high-order structure on deamidation: the deamidation in a small peptide is favored first because of the higher reactivity of the asparagine residue and then because of the stronger stability of the tetrahedral intermediate

Predicting the bioactive conformations of macrocycles: a molecular dynamics-based docking procedure with DynaDock

Macrocyclic compounds are of growing interest as a new class of therapeutics, especially as inhibitors binding to protein-protein interfaces. As molecular modeling is a well-established complimentary tool in modern drug design, the number of attempts to develop reliable docking strategies and algorithms to accurately predict the binding mode of macrocycles is rising continuously. Standard molecular docking approaches need to be adapted to this application, as a comprehensive yet efficient sampling of all ring conformations of the macrocycle is necessary. To overcome this issue, we designed a molecular dynamics-based docking protocol for macrocycles, in which the challenging sampling step is addressed by conventional molecular dynamics (750ns) simulations performed at moderately high temperature (370K). Consecutive flexible docking with the DynaDock approach based on multiple, pre-sampled ring conformations yields highly accurate poses with ligand RMSD values lower than 1.8 angstrom. We further investigated the value of molecular dynamics-based complex stability estimations for pose selection and discuss its applicability in combination with standard binding free energy estimations for assessing the quality of poses in future blind docking studies