Mechanism of inhibition of SARS-CoV-2 Mpro by N3 peptidyl Michael acceptor explained by QM/MM simulations and design of new derivatives with tunable chemical reactivity

The SARS-CoV-2 main protease (M^pro) is essential for replication of the virus responsible for the COVID-19 pandemic, and one of the main targets for drug design. Here, we simulate the inhibition process of SARS-CoV-2 M^pro with a known Michael acceptor (peptidyl) inhibitor, N3. The free energy landscape for the mechanism of the formation of the covalent enzyme-inhibitor product is computed with QM/MM molecular dynamics methods. The simulations show a two-step mechanism, and give structures and calculated barriers in good agreement with experiment. Using these results and information from our previous investigation on the proteolysis reaction of SARS-CoV-2 M^pro, we design two new, synthetically accessible N3-analogues as potential inhibitors, in which the recognition and warhead motifs are modified. QM/MM modelling of the mechanism of inhibition of M^pro by these novel compounds indicates that both may be promising candidates as drug leads against COVID-19, one as an irreversible inhibitor and one as a potential reversible inhibitor.

QM/MM simulations identify the mechanism of reaction of N3, a covalent peptidyl inhibitor of SARS-CoV-2 main protease. Modelling of two novel proposed compounds, B1 and B2, suggests that reversibility of covalent inhibition could be tailored.     

Selective covalent targeting of SARS-CoV-2 main protease by enantiopure chlorofluoroacetamide

The coronavirus disease 2019 (COVID-19) pandemic has necessitated the development of antiviral agents against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The main protease (M^pro) is a promising target for COVID-19 treatment. Here, we report an irreversible SARS-CoV-2 M^pro inhibitor possessing chlorofluoroacetamide (CFA) as a warhead for the covalent modification of M^pro. Ugi multicomponent reaction using chlorofluoroacetic acid enabled the rapid synthesis of dipeptidic CFA derivatives that identified 18 as a potent inhibitor of SARS-CoV-2 M^pro. Among the four stereoisomers, (R,R)-18 exhibited a markedly higher inhibitory activity against M^pro than the other isomers. Reaction kinetics and computational docking studies suggest that the R configuration of the CFA warhead is crucial for the rapid covalent inhibition of M^pro. Our findings highlight the prominent influence of the CFA chirality on the covalent modification of proteinous cysteines and provide the basis for improving the potency and selectivity of CFA-based covalent inhibitors.

Chlorofluoroacetamide (CFA) was used as the warhead for covalent targeting of SARS-CoV-2 M^pro. The chirality at CFA showed marked influence on inhibitory activity, suggesting stereospecific activation of CFA for cysteine modification in the protein.     

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural M^pro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural M^pro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial M^pro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective ‘peptibitors’ inhibit M^pro competitively. Our combined results provide new insights and highlight opportunities for the development of M^pro inhibitors as anti-COVID-19 drugs.

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target. In silico methods reveal structural aspects of how it binds to its 11 natural cleavage sites, the design of novel peptide inhibitors, and insights into drug design.     

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.     

Synergistic inhibition of SARS-CoV-2 cell entry by otamixaban and covalent protease inhibitors: pre-clinical assessment of pharmacological and molecular properties

SARS-CoV-2, the cause of the COVID-19 pandemic, exploits host cell proteins for viral entry into human lung cells. One of them, the protease TMPRSS2, is required to activate the viral spike protein (S). Even though two inhibitors, camostat and nafamostat, are known to inhibit TMPRSS2 and block cell entry of SARS-CoV-2, finding further potent therapeutic options is still an important task. In this study, we report that a late-stage drug candidate, otamixaban, inhibits SARS-CoV-2 cell entry. We show that otamixaban suppresses TMPRSS2 activity and SARS-CoV-2 infection of a human lung cell line, although with lower potency than camostat or nafamostat. In contrast, otamixaban inhibits SARS-CoV-2 infection of precision cut lung slices with the same potency as camostat. Furthermore, we report that otamixaban''s potency can be significantly enhanced by (sub-) nanomolar nafamostat or camostat supplementation. Dominant molecular TMPRSS2-otamixaban interactions are assessed by extensive 109 μs of atomistic molecular dynamics simulations. Our findings suggest that combinations of otamixaban with supplemental camostat or nafamostat are a promising option for the treatment of COVID-19.

SARS-CoV-2, the cause of the COVID-19 pandemic, exploits host proteins for viral entry into human lung cells and is blocked by otamixaban in combination with a covalent protease inhibitor.     

Interaction of the prototypical α-ketoamide inhibitor with the SARS-CoV-2 main protease active site in silico: Molecular dynamic simulations highlight the stability of the ligand-protein complex

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes an illness known as COVID-19, which has been declared a global pandemic with over 2 million confirmed cases and 137,000 deaths in 185 countries and regions at the time of writing (16 April 2020), over a quarter of these cases being in the United States. In the absence of a vaccine, or an approved effective therapeutic, there is an intense interest in repositioning available drugs or designing small molecule antivirals. In this context, in silico modelling has proven to be an invaluable tool. An important target is the SARS-CoV-2 main protease (M^pro), involved in processing translated viral proteins. Peptidomimetic α-ketoamides represent prototypical inhibitors of M^pro. A recent attempt at designing a compound with enhanced pharmacokinetic properties has resulted in the synthesis and evaluation of the α-ketoamide 13b analogue. Here, we performed molecular docking and molecular dynamics simulations to further characterize the interaction of α-ketoamide 13b with the active site of the SARS-CoV-2 M^pro. We included the widely used antibiotic, amoxicillin, for comparison. Our findings indicate that α-ketoamide 13b binds more tightly (predicted GlideScore = -8.7 and -9.2 kcal/mol for protomers A and B, respectively), to the protease active site compared to amoxicillin (-5.0 and -4.8 kcal/mol). Further, molecular dynamics simulations highlight the stability of the interaction of the α-ketoamide 13b ligand with the SARS-CoV-2 M^pro (ΔG = -25.2 and -22.3 kcal/mol for protomers A and B). In contrast, amoxicillin interacts unfavourably with the protease (ΔG = +32.8 kcal/mol for protomer A), with unbinding events observed in several independent simulations. Overall, our findings are consistent with those previously observed, and highlight the need to further explore the α-ketoamides as potential antivirals for this ongoing COVID-19 pandemic.     

Synthesis and Biochemical Evaluation of 8H-Indeno[1,2-d]thiazole Derivatives as Novel SARS-CoV-2 3CL Protease Inhibitors

The COVID-19 pandemic caused by SARS-CoV-2 is a global burden on human health and economy. The 3-Chymotrypsin-like cysteine protease (3CL^pro) becomes an attractive target for SARS-CoV-2 due to its important role in viral replication. We synthesized a series of 8H-indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CL^pro. Among them, the representative compound 7a displayed inhibitory activity with an IC₅₀ of 1.28 ± 0.17 μM against SARS-CoV-2 3CL^pro. Molecular docking of 7a against 3CL^pro was performed and the binding mode was rationalized. These preliminary results provide a unique prototype for the development of novel inhibitors against SARS-CoV-2 3CL^pro.     

Translating daily COVID-19 screening into a simple glucose test: a proof of concept study

Home testing is an attractive emerging strategy to combat the COVID-19 pandemic and prevent overloading of healthcare resources through at-home isolation, screening and monitoring of symptoms. However, current diagnostic technologies of SARS-CoV-2 still suffer from some drawbacks because of the tradeoffs between sensitivity, usability and costs, making the test unaffordable to most users at home. To address these limitations, taking advantage of clustered regularly interspaced short palindromic repeats (CRISPRs) and a portable glucose meter (PGM), we present a proof-of-concept demonstration of a target-responsive CRISPR-PGM system for translating SARS-CoV-2 detection into a glucose test. Using this system, a specific N gene, N protein, and pseudo-viruses of SARS-CoV-2 have been detected quantitatively with a PGM. Given the facile integration of various bioreceptors into the CRISPR-PGM system, the proposed method provides a starting point to provide patients with a single-device solution that can quantitatively monitor multiple COVID-19 biomarkers at home.

COVID-19 glucose test: translating SARS-CoV-2 detection into a glucose test is achieved by incorporating target-responsive rolling circle amplification and a CRISPR-based collateral cleavage module with a portable glucose meter.     

Orally administered bismuth drug together with N-acetyl cysteine as a broad-spectrum anti-coronavirus cocktail therapy

The emergence of SARS-CoV-2 variants of concern compromises vaccine efficacy and emphasizes the need for further development of anti-SARS-CoV-2 therapeutics, in particular orally administered take-home therapies. Cocktail therapy has shown great promise in the treatment of viral infection. Herein, we reported the potent preclinical anti-SARS-CoV-2 efficacy of a cocktail therapy consisting of clinically used drugs, e.g. colloidal bismuth subcitrate (CBS) or bismuth subsalicylate (BSS), and N-acetyl-l-cysteine (NAC). Oral administration of the cocktail reduced viral loads in the lung and ameliorated virus-induced pneumonia in a hamster infection model. The mechanistic studies showed that NAC prevented the hydrolysis of bismuth drugs at gastric pH via the formation of the stable component [Bi(NAC)₃], and optimized the pharmacokinetics profile of CBS in vivo. Combination of bismuth drugs with NAC suppressed the replication of a panel of medically important coronaviruses including Middle East respiratory syndrome-related coronavirus (MERS-CoV), Human coronavirus 229E (HCoV-229E) and SARS-CoV-2 Alpha variant (B.1.1.7) with broad-spectrum inhibitory activities towards key viral cysteine enzymes/proteases including papain-like protease (PL^pro), main protease (M^pro), helicase (Hel) and angiotensin-converting enzyme 2 (ACE2). Importantly, our study offered a potential at-home treatment for combating SARS-CoV-2 and future coronavirus infections.

A cocktail therapy comprising bismuth drugs and N-acetyl-l-cysteine is reported to suppress the replication of SARS-CoV-2 via the oral route. The broad-spectrum inhibitory activities of the combination upon key viral cysteine enzymes are verified.     

Interaction of small molecules with the SARS-CoV-2 main protease in silico and in vitro validation of potential lead compounds using an enzyme-linked immunosorbent assay

Caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the COVID-19 pandemic is ongoing, with no proven safe and effective vaccine to date. Further, effective therapeutic agents for COVID-19 are limited, and as a result, the identification of potential small molecule antiviral drugs is of particular importance. A critical antiviral target is the SARS-CoV-2 main protease (M^pro), and our aim was to identify lead compounds with potential inhibitory effects. We performed an initial molecular docking screen of 300 small molecules, which included phenolic compounds and fatty acids from our OliveNet™ library (224), and an additional group of curated pharmacological and dietary compounds. The prototypical α-ketoamide 13b inhibitor was used as a control to guide selection of the top 30 compounds with respect to binding affinity to the M^pro active site. Further studies and analyses including blind docking were performed to identify hypericin, cyanidin-3-O-glucoside and SRT2104 as potential leads. Molecular dynamics simulations demonstrated that hypericin (ΔG = -18.6 and -19.3 kcal/mol), cyanidin-3-O-glucoside (ΔG = -50.8 and -42.1 kcal/mol), and SRT2104 (ΔG = -8.7 and -20.6 kcal/mol), formed stable interactions with the M^pro active site. An enzyme-linked immunosorbent assay indicated that, albeit, not as potent as the covalent positive control (GC376), our leads inhibited the M^pro with activity in the micromolar range, and an order of effectiveness of hypericin and cyanidin-3-O-glucoside > SRT2104 > SRT1720. Overall, our findings, and those highlighted by others indicate that hypericin and cyanidin-3-O-glucoside are suitable candidates for progress to in vitro and in vivo antiviral studies.     

In-silico investigation of phenolic compounds from leaves of Phillyrea angustifolia L. as a potential inhibitor against the SARS-CoV-2 main protease (Mpro PDB ID:5R83) using a virtual screening method

There is currently a global COVID-19 pandemic caused by the severe acute respiratory syndrome Coronavirus-2 (SARS-CoV-2) and its variants. This highly contagious viral disease continues to pose a major health threat global. The discovery of vaccinations is not enough to prevent their spread and dire consequences. To take advantage of the current drugs and isolated compounds, and immediately qualifying approach is required. The aim of our research is evaluation the potency for natural antiviral compounds against the SARS CoV-2 M^pro. Molecular docking of four phenolic compounds from Phillyrea angustifolia leaves with SARS-CoV-2 M^pro has been conducted. Similarly, the stability of selected ligand–protein interactions has been determined using MD simulations. Moreover, the quantitative structure–activity relationship (QSAR), MMGBSA binding energies, pharmacokinetics, and drug-likeness predictions for selected phenolic have been reported. The selected phenolic compounds (Luteolin-7-O-glucoside, Apigenin-7-O-glucoside, Demethyl-oleuropein, and Oleuropein aglycone) revealed strong binding contacts in the two active pockets of a target protein of SARS-CoV-2 M^pro with the docking scores and highest binding energies with a binding energy of ?8.2 kcal/mol; ?7.8 kcal/mol; ?7.2 kcal/mol and ?7.0 kcal/mol respectively. Both Demethyloleoeuropein and Oleuropein aglycone can interact with residues His41 and Cys145 (catalytic dyad) and other amino acids of the binding pocket of M^pro. According to QSAR, studies on pharmacokinetics and drug-like properties suggested that oleuropein aglycone could be the best inhibitor of SARS-CoV-2 for new drug design and development. Further in vivo, in vitro, and clinical studies are highly needed to examine the potential of these phenolic compounds in the fight against COVID-19.     

Multiplex metal-detection based assay (MMDA) for COVID-19 diagnosis and identification of disease severity biomarkers

The ongoing COVID-19 pandemic caused by SARS-CoV-2 highlights the urgent need to develop sensitive methods for diagnosis and prognosis. To achieve this, multidimensional detection of SARS-CoV-2 related parameters including virus loads, immune response, and inflammation factors is crucial. Herein, by using metal-tagged antibodies as reporting probes, we developed a multiplex metal-detection based assay (MMDA) method as a general multiplex assay strategy for biofluids. This strategy provides extremely high multiplexing capability (theoretically over 100) compared with other reported biofluid assay methods. As a proof-of-concept, MMDA was used for serologic profiling of anti-SARS-CoV-2 antibodies. The MMDA exhibits significantly higher sensitivity and specificity than ELISA for the detection of anti-SARS-CoV-2 antibodies. By integrating the high dimensional data exploration/visualization tool (tSNE) and machine learning algorithms with in-depth analysis of multiplex data, we classified COVID-19 patients into different subgroups based on their distinct antibody landscape. We unbiasedly identified anti-SARS-CoV-2-nucleocapsid IgG and IgA as the most potently induced types of antibodies for COVID-19 diagnosis, and anti-SARS-CoV-2-spike IgA as a biomarker for disease severity stratification. MMDA represents a more accurate method for the diagnosis and disease severity stratification of the ongoing COVID-19 pandemic, as well as for biomarker discovery of other diseases.

A MMDA platform is developed by using metal-tagged antibodies as reporting probes combined with machine learning algorithms, as a general strategy for highly multiplexed biofluid assay.     

Mechanism of covalent binding of ibrutinib to Bruton's tyrosine kinase revealed by QM/MM calculations

Ibrutinib is the first covalent inhibitor of Bruton''s tyrosine kinase (BTK) to be used in the treatment of B-cell cancers. Understanding the mechanism of covalent inhibition will aid in the design of safer and more selective covalent inhibitors that target BTK. The mechanism of covalent inhibition in BTK has been uncertain because there is no appropriate residue nearby that can act as a base to deprotonate the cysteine thiol prior to covalent bond formation. We investigate several mechanisms of covalent modification of C481 in BTK by ibrutinib using combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics reaction simulations. The lowest energy pathway involves direct proton transfer from C481 to the acrylamide warhead in ibrutinib, followed by covalent bond formation to form an enol intermediate. There is a subsequent rate-limiting keto–enol tautomerisation step (ΔG^‡ = 10.5 kcal mol⁻¹) to reach the inactivated BTK/ibrutinib complex. Our results represent the first mechanistic study of BTK inactivation by ibrutinib to consider multiple mechanistic pathways. These findings should aid in the design of covalent drugs that target BTK and other similar targets.

QM/MM simulations show that covalent modification of BTK by ibrutinib proceeds via an intramolecular proton transfer from C481 to the acrylamide warhead of ibrutinib, followed by covalent bond formation and subsequent keto–enol tautomerisation.      

Design of Novel Enantiopure Dispirooxindolopyrrolidine-Piperidones as Promising Candidates toward COVID-19: Asymmetric Synthesis,Crystal Structure and In Silico Studies

Despite the effectiveness of COVID-19 vaccines, there is still an urgent need for discovering new anti-viral drugs to address the awful spread and transmission of the rapidly modifiable virus. In this study, the ability of a small library of enantiomerically pure spirooxindolopyrrolidine-grafted piperidones to inhibit the main protease of SARS-CoV-2 (M^pro) is evaluated. These spiroheterocycles were synthesized by 1,3-dipolar cycloaddition of various stabilized azomethine ylides with chiral dipolarophiles derived from N-[(S)-(-)-methylbenzyl]-4-piperidone. The absolute configuration of contiguous carbons was confirmed by a single crystal X-ray diffraction analysis. The binding of these compounds to SARS-CoV-2 M^pro was investigated using molecular docking and molecular dynamics simulation. Three compounds 4a, 4b and 4e exhibited stable binding modes interacting with the key subsites of the substrate-binding pocket of SARS-CoV-2 M^pro. The synthesized compounds represent potential leads for the development of novel inhibitors of SARS-CoV-2 main protease protein for COVID-19 treatment.     

The role of the active site tyrosine in the mechanism of lytic polysaccharide monooxygenase

Catalytic breakdown of polysaccharides can be achieved more efficiently by means of the enzymes lytic polysaccharide monooxygenases (LPMOs). However, the LPMO mechanism has remained controversial, preventing full exploitation of their potential. One of the controversies has centered around an active site tyrosine, present in most LPMO classes. Recent investigations have for the first time obtained direct (spectroscopic) evidence for the possibility of chemical modification of this tyrosine. However, the spectroscopic features obtained in the different investigations are remarkably different, with absorption maximum at 420 and 490 nm, respectively. In this paper we use density functional theory (DFT) in a QM/MM formulation to reconcile these (apparently) conflicting results. By modeling the spectroscopy as well as the underlying reaction mechanism we can show how formation of two isomers (both involving deprotonation of tyrosine) explains the difference in the observed spectroscopic features. Both isomers have a [TyrO–Cu–OH]⁺ moiety with the OH in either the cis- or trans-position to a deprotonated tyrosine. Although the cis-[TyrO–Cu–OH]⁺ moiety is well positioned for oxidation of the substrate, preliminary calculations with the substrate reveal that the reactivity is at best moderate, making a protective role of tyrosine more likely.

With QM/MM, we investigate the mechanism of tyrosine deprotonation in lytic polysaccharide monooxygenases. Our results support deprotonation and our calculated UV-vis spectra show that two isomers must be formed to match recent experiments.     

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein,Mpro and PLpro of SARS-CoV-2: Protein–Protein Interactions,Dynamics Simulations and Free Energy Calculations

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (−16.8 ± 0.02 kcal/mol, −12.3 ± 0.03 kcal/mol and −13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (M^pro) and the papainlike protease (PL^pro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.     

High-resolution mining of the SARS-CoV-2 main protease conformational space: supercomputer-driven unsupervised adaptive sampling

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs). The global exploration problem is decomposed into a set of separate MD trajectories that can be restarted within a selective process to achieve sufficient phase-space sampling. Accurate statistical properties can be obtained through reweighting. Within this highly parallel setup, the Tinker-HP package can be powered by an arbitrary large number of GPUs on supercomputers, reducing exploration time from years to days. This approach is used to tackle the urgent modeling problem of the SARS-CoV-2 Main Protease (M^pro) producing more than 38 μs of all-atom simulations of its apo (ligand-free) dimer using the high-resolution AMOEBA PFF. The first 15.14 μs simulation (physiological pH) is compared to available non-PFF long-timescale simulation data. A detailed clustering analysis exhibits striking differences between FFs, with AMOEBA showing a richer conformational space. Focusing on key structural markers related to the oxyanion hole stability, we observe an asymmetry between protomers. One of them appears less structured resembling the experimentally inactive monomer for which a 6 μs simulation was performed as a basis for comparison. Results highlight the plasticity of the M^pro active site. The C-terminal end of its less structured protomer is shown to oscillate between several states, being able to interact with the other protomer, potentially modulating its activity. Active and distal site volumes are found to be larger in the most active protomer within our AMOEBA simulations compared to non-PFFs as additional cryptic pockets are uncovered. A second 17 μs AMOEBA simulation is performed with protonated His172 residues mimicking lower pH. Data show the protonation impact on the destructuring of the oxyanion loop. We finally analyze the solvation patterns around key histidine residues. The confined AMOEBA polarizable water molecules are able to explore a wide range of dipole moments, going beyond bulk values, leading to a water molecule count consistent with experimental data. Results suggest that the use of PFFs could be critical in drug discovery to accurately model the complexity of the molecular interactions structuring M^pro.

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs).