{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"institution":[{"name":"bioRxiv"}],"indexed":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T18:45:53Z","timestamp":1768502753976,"version":"3.49.0"},"posted":{"date-parts":[[2022,8,1]]},"group-title":"Biochemistry","reference-count":50,"publisher":"openRxiv","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"accepted":{"date-parts":[[2023,3,6]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                <jats:p>Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of PAXLOVID, a drug approved by FDA for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir\u2019s inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds up in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here we investigated the H172Y Mpro\u2019s conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggests that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weakens the nirma-trelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intel-ligence approaches, and together with biochemical experiments they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the discovery of new antiviral drugs. The presented workflow can be applicable to characterize mutation effects on any protein drug targets.<\/jats:p>","DOI":"10.1101\/2022.07.31.502215","type":"posted-content","created":{"date-parts":[[2022,8,1]],"date-time":"2022-08-01T22:00:16Z","timestamp":1659391216000},"source":"Crossref","is-referenced-by-count":1,"title":["An Integrative Approach to Dissect the Drug Resistance Mechanism of the H172Y Mutation of SARS-CoV-2 Main Protease"],"prefix":"10.64898","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-2652-5994","authenticated-orcid":false,"given":"Joseph","family":"Clayton","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-0927-3825","authenticated-orcid":false,"given":"Vinicius Martins","family":"de Oliveira","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Mohamed Fouad","family":"Ibraham","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2761-2862","authenticated-orcid":false,"given":"Xinyuanyuan","family":"Sun","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Paween","family":"Mahinthichaichan","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7461-3764","authenticated-orcid":false,"given":"Mingzhe","family":"Shen","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8850-2977","authenticated-orcid":false,"given":"Rolf","family":"Hilgenfeld","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3234-0769","authenticated-orcid":false,"given":"Jana","family":"Shen","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"54368","reference":[{"key":"2023030900400845000_2022.07.31.502215v2.1","unstructured":"U.S. Food and Drug Administration, Emergency Use Authorization for Paxlovid [Fact Sheet] https:\/\/www.fda.gov\/media\/155050\/download; 2021."},{"key":"2023030900400845000_2022.07.31.502215v2.2","doi-asserted-by":"publisher","DOI":"10.1007\/s40265-022-01692-5"},{"key":"2023030900400845000_2022.07.31.502215v2.3","doi-asserted-by":"publisher","DOI":"10.1056\/nejmoa2118542"},{"key":"2023030900400845000_2022.07.31.502215v2.4","doi-asserted-by":"crossref","unstructured":"Zhang, L. ; Lin, D. ; Sun, X. ; Curth, U. ; Drosten, C. ; Sauerhering, L. ; Becker, S. ; Rox, K. ; Hilgenfeld, R. Crystal Structure of SARS-CoV-2 Main Protease Provides a Basis for Design of Improved \u03b1-Ketoamide Inhibitors. Science 2020, 368, 409\u2013412.","DOI":"10.1126\/science.abb3405"},{"key":"2023030900400845000_2022.07.31.502215v2.5","doi-asserted-by":"crossref","first-page":"689","DOI":"10.1007\/s13238-021-00883-2","article-title":"Crystal Structure of SARS-CoV-2 Main Protease in Complex with Protease Inhibitor PF-07321332","volume":"13","year":"2022","journal-title":"Protein Cell"},{"key":"2023030900400845000_2022.07.31.502215v2.6","doi-asserted-by":"publisher","DOI":"10.1016\/j.jbc.2022.101972"},{"key":"2023030900400845000_2022.07.31.502215v2.7","doi-asserted-by":"crossref","first-page":"498","DOI":"10.1038\/s41422-022-00640-y","article-title":"The P132H Mutation in the Main Protease of Omicron SARS-CoV-2 Decreases Thermal Stability without Compromising Catalysis or SmallMolecule Drug Inhibition","volume":"32","year":"2022","journal-title":"Cell Res"},{"key":"2023030900400845000_2022.07.31.502215v2.8","doi-asserted-by":"publisher","DOI":"10.1016\/j.bmcl.2022.128629"},{"key":"2023030900400845000_2022.07.31.502215v2.9","doi-asserted-by":"publisher","DOI":"10.1002\/gch2.1018"},{"key":"2023030900400845000_2022.07.31.502215v2.10","doi-asserted-by":"crossref","unstructured":"Hu, Y. ; Lewandowski, E. M. ; Tan, H. ; Morgan, R. T. ; Zhang, X. ; Jacobs, L. M. ; Butler, S. G. ; Mongora, M. V. ; Choy, J. S. ; Chen, Y. , et al. Naturally occurring mutations of SARS-CoV-2 main protease confer drug resistance to nirmatrelvir. bioRxiv 2022,","DOI":"10.1101\/2022.06.28.497978"},{"key":"2023030900400845000_2022.07.31.502215v2.11","first-page":"5.6.1","article-title":"Comparative Protein Structure Modeling Using MODELLER","volume":"54","year":"2016","journal-title":"Current Protocols in Bioinformatics"},{"key":"2023030900400845000_2022.07.31.502215v2.12","doi-asserted-by":"crossref","first-page":"1372","DOI":"10.1021\/acs.jcim.8b00227","article-title":"Generalized Born Based Continuous Constant pH Molecular Dynamics in Amber: Implementation, Benchmarking and Analysis","volume":"58","year":"2018","journal-title":"J. Chem. Inf. Model"},{"key":"2023030900400845000_2022.07.31.502215v2.13","doi-asserted-by":"crossref","first-page":"4821","DOI":"10.1021\/acs.jcim.9b00754","article-title":"GPU-Accelerated Implementation of Continuous Constant pH Molecular Dynamics in Amber: pKa Predictions with Single-pH Simulations","volume":"59","year":"2019","journal-title":"J. Chem. Inf. Model"},{"key":"2023030900400845000_2022.07.31.502215v2.14","doi-asserted-by":"crossref","first-page":"115101","DOI":"10.1063\/5.0020458","article-title":"Assessment of Proton-Coupled Conformational Dynamics of SARS and MERS Coronavirus Papain-like Proteases: Implication for Designing Broad-Spectrum Antiviral Inhibitors","volume":"153","year":"2020","journal-title":"J. Chem. Phys"},{"key":"2023030900400845000_2022.07.31.502215v2.15","doi-asserted-by":"crossref","unstructured":"Verma, N. ; Henderson, J. A. ; Shen, J. Proton-Coupled Conformational Activation of SARS Coronavirus Main Proteases and Opportunity for Designing Small-Molecule Broad-Spectrum Targeted Covalent Inhibitors. J. Am. Chem. Soc. 2020, 21883\u201321890.","DOI":"10.1021\/jacs.0c10770"},{"key":"2023030900400845000_2022.07.31.502215v2.16","doi-asserted-by":"crossref","first-page":"1513","DOI":"10.1039\/D0SC04942E","article-title":"Inhibitor Binding Influences the Protonation States of Histidines in SARS-CoV-2 Main Protease","volume":"12","year":"2021","journal-title":"Chem. Sci"},{"key":"2023030900400845000_2022.07.31.502215v2.17","unstructured":"Case, D. A. et al. AMBER 2020. 2020."},{"key":"2023030900400845000_2022.07.31.502215v2.18","doi-asserted-by":"publisher","DOI":"10.1073\/pnas.1835675100"},{"key":"2023030900400845000_2022.07.31.502215v2.19","doi-asserted-by":"publisher","DOI":"10.1016\/j.jmb.2005.09.012"},{"key":"2023030900400845000_2022.07.31.502215v2.20","doi-asserted-by":"crossref","first-page":"830","DOI":"10.1002\/prot.22921","article-title":"Role of conformational sampling in computing mutation-induced changes in protein structure and stability","volume":"79","year":"2011","journal-title":"Proteins: Structure, Function, and Bioinformatics"},{"key":"2023030900400845000_2022.07.31.502215v2.21","doi-asserted-by":"publisher","DOI":"10.1063\/1.1740409"},{"key":"2023030900400845000_2022.07.31.502215v2.22","doi-asserted-by":"crossref","first-page":"9795","DOI":"10.1021\/jp011878v","article-title":"Can Absolute Free Energies of Association Be Estimated from Molecular Mechanical Simulations? The BiotinStreptavidin System Revisited","volume":"105","year":"2001","journal-title":"The Journal of Physical Chemistry A"},{"key":"2023030900400845000_2022.07.31.502215v2.23","doi-asserted-by":"crossref","first-page":"3794","DOI":"10.1021\/acs.jcim.9b00362","article-title":"Computing Relative Binding Affinity of Ligands to Receptor: An Effective Hybrid Single-Dual-Topology Free-Energy Perturbation Approach in NAMD","volume":"59","year":"2019","journal-title":"J. Chem. Inf. Model"},{"key":"2023030900400845000_2022.07.31.502215v2.24","doi-asserted-by":"publisher","DOI":"10.1002\/jcc.20289"},{"key":"2023030900400845000_2022.07.31.502215v2.25","doi-asserted-by":"publisher","DOI":"10.1021\/acs.jpcb.7b11367"},{"key":"2023030900400845000_2022.07.31.502215v2.26","doi-asserted-by":"crossref","first-page":"3023","DOI":"10.1038\/s41467-021-23246-1","article-title":"Deep Learning the Structural Determinants of Protein Biochemical Properties by Comparing Structural Ensembles with DiffNets","volume":"12","year":"2021","journal-title":"Nature Communications"},{"key":"2023030900400845000_2022.07.31.502215v2.27","doi-asserted-by":"publisher","DOI":"10.1110\/ps.180101"},{"key":"2023030900400845000_2022.07.31.502215v2.28","doi-asserted-by":"publisher","DOI":"10.1093\/nar\/gkp253"},{"key":"2023030900400845000_2022.07.31.502215v2.29"},{"key":"2023030900400845000_2022.07.31.502215v2.30","doi-asserted-by":"publisher","DOI":"10.1063\/1.445869"},{"key":"2023030900400845000_2022.07.31.502215v2.31","doi-asserted-by":"publisher","DOI":"10.1002\/jcc.20035"},{"key":"2023030900400845000_2022.07.31.502215v2.32","doi-asserted-by":"publisher","DOI":"10.1063\/5.0019056"},{"key":"2023030900400845000_2022.07.31.502215v2.33","doi-asserted-by":"publisher","DOI":"10.1063\/1.464397"},{"key":"2023030900400845000_2022.07.31.502215v2.34","unstructured":"Schr\u00f6dinger, LLC,"},{"key":"2023030900400845000_2022.07.31.502215v2.35"},{"key":"2023030900400845000_2022.07.31.502215v2.36","doi-asserted-by":"publisher","DOI":"10.1021\/cr00023a004"},{"key":"2023030900400845000_2022.07.31.502215v2.37","doi-asserted-by":"publisher","DOI":"10.1021\/ct300400x"},{"key":"2023030900400845000_2022.07.31.502215v2.38","doi-asserted-by":"publisher","DOI":"10.1038\/nmeth.4067"},{"key":"2023030900400845000_2022.07.31.502215v2.39","doi-asserted-by":"publisher","DOI":"10.1021\/ci300363c"},{"key":"2023030900400845000_2022.07.31.502215v2.40","doi-asserted-by":"publisher","DOI":"10.1021\/ci3003649"},{"key":"2023030900400845000_2022.07.31.502215v2.41","doi-asserted-by":"publisher","DOI":"10.1016\/0021-9991(77)90098-5"},{"key":"2023030900400845000_2022.07.31.502215v2.42","doi-asserted-by":"publisher","DOI":"10.1063\/1.467468"},{"key":"2023030900400845000_2022.07.31.502215v2.43","doi-asserted-by":"publisher","DOI":"10.1063\/1.470648"},{"key":"2023030900400845000_2022.07.31.502215v2.44","doi-asserted-by":"publisher","DOI":"10.1016\/0263-7855(96)00018-5"},{"key":"2023030900400845000_2022.07.31.502215v2.45","doi-asserted-by":"crossref","first-page":"2606","DOI":"10.1021\/ct300242f","article-title":"A Toolkit for the Analysis of Free-Energy Perturbation Calculations","volume":"8","year":"2012","journal-title":"J. Chem. Theory Comput"},{"key":"2023030900400845000_2022.07.31.502215v2.46","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1016\/0021-9991(76)90078-4","article-title":"Efficient estimation of free energy differences from Monte Carlo data","volume":"22","year":"1976","journal-title":"J. Comput. Phys"},{"key":"2023030900400845000_2022.07.31.502215v2.47","doi-asserted-by":"publisher","DOI":"10.1021\/jp102971x"},{"key":"2023030900400845000_2022.07.31.502215v2.48","doi-asserted-by":"crossref","first-page":"1708","DOI":"10.1021\/acscentsci.8b00717","article-title":"Accurate Estimation of Ligand Binding Affinity Changes upon Protein Mutation","volume":"4","year":"2018","journal-title":"ACS Cent. Sci"},{"key":"2023030900400845000_2022.07.31.502215v2.49","doi-asserted-by":"publisher","DOI":"10.1016\/j.jmb.2008.05.023"},{"key":"2023030900400845000_2022.07.31.502215v2.50","doi-asserted-by":"publisher","DOI":"10.1006\/abio.1998.3014"}],"container-title":[],"original-title":[],"link":[{"URL":"https:\/\/syndication.highwire.org\/content\/doi\/10.1101\/2022.07.31.502215","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T03:26:00Z","timestamp":1768447560000},"score":1,"resource":{"primary":{"URL":"http:\/\/biorxiv.org\/lookup\/doi\/10.1101\/2022.07.31.502215"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,8,1]]},"references-count":50,"URL":"https:\/\/doi.org\/10.1101\/2022.07.31.502215","relation":{"is-preprint-of":[{"id-type":"doi","id":"10.1021\/acs.jcim.3c00344","asserted-by":"subject"}]},"subject":[],"published":{"date-parts":[[2022,8,1]]},"subtype":"preprint"}}