Insight into the Binding Energetics of Targeted Reversible Covalent Inhibitors of the SARS-CoV-2 Main Protease

The main protease (Mpro) of the SARS-CoV-2 virus is an attractive therapeutic target for developing antivirals to combat COVID-19. Mpro is essential for the replication cycle of the SARS-CoV-2 virus, so inhibiting Mpro blocks a vital piece of the cell replication machinery of the virus. A promising strategy to disrupt the viral replication cycle is to design inhibitors that bind to the active site cysteine (Cys145) of the Mpro. Cysteine targeted covalent inhibitors are gaining traction in drug discovery owing to the benefits of improved potency and extended drug-target engagement. An interesting aspect of these inhibitors is that they can be chemically tuned to form a covalent, but reversible bond, with their targets of interest. Several small-molecule cysteine-targeting covalent inhibitors of the Mpro have been discovered—some of which are currently undergoing evaluation in early phase human clinical trials. Understanding the binding energetics of these inhibitors could provide new insights to facilitate the design of potential drug candidates against COVID-19. Motivated by this, we employed rigorous absolute binding free energy calculations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to estimate the energetics of binding of some promising reversible covalent inhibitors of the Mpro. We find that the inclusion of enhanced sampling techniques such as replica-exchange algorithm in binding free energy calculations can improve the convergence of predicted non-covalent binding free energy estimates of inhibitors binding to the Mpro target. In addition, our results indicate that binding free energy calculations coupled with multiscale simulations can be a useful approach to employ in ranking covalent inhibitors to their targets. This approach may be valuable in prioritizing and refining covalent inhibitor compounds for lead discovery efforts against COVID-19 and future coronavirus infections.

Insight into the Binding Energetics of Targeted Reversible Covalent Inhibitors of the SARS-CoV-2 Main Protease

The main protease (Mpro) of the SARS-CoV-2 virus is an attractive therapeutic target for developing antivirals to combat COVID-19. Mpro is essential for the replication cycle of the SARS-CoV-2 virus, so inhibiting Mpro blocks a vital piece of the cell replication machinery of the virus. A promising strategy to disrupt the viral replication cycle is to design inhibitors that bind to the active site cysteine (Cys145) of the Mpro. Cysteine targeted covalent inhibitors are gaining traction in drug discovery owing to the benefits of improved potency and extended drug-target engagement. An interesting aspect of these inhibitors is that they can be chemically tuned to form a covalent, but reversible bond, with their targets of interest. Several small-molecule cysteine-targeting covalent inhibitors of the Mpro have been discovered—some of which are currently undergoing evaluation in early phase human clinical trials. Understanding the binding energetics of these inhibitors could provide new insights to facilitate the design of potential drug candidates against COVID-19. Motivated by this, we employed rigorous absolute binding free energy calculations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to estimate the energetics of binding of some promising reversible covalent inhibitors of the Mpro. We find that the inclusion of enhanced sampling techniques such as replica-exchange algorithm in binding free energy calculations can improve the convergence of predicted non-covalent binding free energy estimates of inhibitors binding to the Mpro target. In addition, our results indicate that binding free energy calculations coupled with multiscale simulations can be a useful approach to employ in ranking covalent inhibitors to their targets. This approach may be valuable in prioritizing and refining covalent inhibitor compounds for lead discovery efforts against COVID-19 and future coronavirus infections.

Molecular insights into receptor binding energetics and neutralization of SARS-CoV-2 variants

Despite an unprecedented global gain in knowledge since the emergence of SARS-CoV-2, almost all mechanistic knowledge related to the molecular and cellular details of viral replication, pathology and virulence has been generated using early prototypic isolates of SARS-CoV-2. Here, using atomic force microscopy and molecular dynamics, we investigated how these mutations quantitatively affected the kinetic, thermodynamic and structural properties of RBD—ACE2 complex formation. We observed for several variants of concern a significant increase in the RBD—ACE2 complex stability. While the N501Y and E484Q mutations are particularly important for the greater stability, the N501Y mutation is unlikely to significantly affect antibody neutralization. This work provides unprecedented atomistic detail on the binding of SARS-CoV-2 variants and provides insight into the impact of viral mutations on infection-induced immunity.

Molecular insights into receptor binding energetics and neutralization of SARS-CoV-2 variants

Despite an unprecedented global gain in knowledge since the emergence of SARS-CoV-2, almost all mechanistic knowledge related to the molecular and cellular details of viral replication, pathology and virulence has been generated using early prototypic isolates of SARS-CoV-2. Here, using atomic force microscopy and molecular dynamics, we investigated how these mutations quantitatively affected the kinetic, thermodynamic and structural properties of RBD—ACE2 complex formation. We observe a direct link between increased RBD—ACE2 complex stability and the greater transmissibility observed for the variants of concern. While the N501Y and E484Q mutations are particularly important for the greater stability, the N501Y mutation is unlikely to significantly affect antibody neutralization. This work provides unprecedented atomistic detail on the binding of SARS-CoV-2 variants and provides insight into the impact of viral mutations on infection-induced immunity.

Energetics of a protein disorder-order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered lid region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Synaptotagmin-1 membrane binding is driven by the C2B domain and assisted cooperatively by the C2A domain

Synaptotagmin interaction with anionic lipid (phosphatidylserine/phosphatidylinositol) containing membranes, both in the absence and presence of calcium ions (Ca2+), is critical to its central role in orchestrating neurotransmitter release. The molecular surfaces involved, namely the conserved polylysine motif in the C2B domain and Ca2+-binding aliphatic loops on both C2A and C2B domains, are known. Here we use surface force apparatus combined with systematic mutational analysis of the functional surfaces to directly measure Syt1-membrane interaction and fully map the site-binding energetics of Syt1 both in the absence and presence of Ca2+. By correlating energetics data with the molecular rearrangements measured during confinement, we find that both C2 domains cooperate in membrane binding, with the C2B domain functioning as the main energetic driver, and the C2A domain acting as a facilitator.

Optimization of a bacterial three-hybrid assay through in vivo titration of an RNA-DNA adapter-protein

ABSTRACTNon-coding RNAs regulate gene expression in every domain of life. In bacteria, small RNAs (sRNAs) regulate gene expression in response to stress and are often assisted by RNA-chaperone proteins, such as Hfq. We have recently developed a bacterial three-hybrid (B3H) assay that detects the strong binding interactions of certain E. coli sRNAs with proteins Hfq and ProQ. Despite the promise of this system, the signal-to-noise has made it challenging to detect weaker interactions. In this work, we use Hfq-sRNA interactions as a model system to optimize the B3H assay, so that weaker RNA-protein interactions can be more reliably detected. We find that the concentration of the RNA-DNA adapter is an important parameter in determining the signal in the system, and have modified the plasmid expressing this component to tune its concentration to optimal levels. In addition, we have systematically perturbed the binding affinity of Hfq-RNA interactions to define, for the first time, the relationship between B3H signal and in vitro binding energetics. The new pAdapter construct presented here substantially expands the range of detectable interactions in the B3H assay, broadening its utility. This improved assay will increase the likelihood of identifying novel protein-RNA interactions with the B3H system, and will facilitate exploration of the binding mechanisms of these interactions.

Polyamine blockade and binding energetics in the MthK potassium channel

Polyamines such as spermidine and spermine are found in nearly all cells, at concentrations ranging up to 0.5 mM. These cations are endogenous regulators of cellular K+ efflux, binding tightly in the pores of inwardly rectifying K+ (Kir) channels in a voltage-dependent manner. Although the voltage dependence of Kir channel polyamine blockade is thought to arise at least partially from the energetically coupled movements of polyamine and K+ ions through the pore, the nature of physical interactions between these molecules is unclear. Here we analyze the polyamine-blocking mechanism in the model K+ channel MthK, using a combination of electrophysiology and computation. Spermidine (SPD3+) and spermine (SPM4+) each blocked current through MthK channels in a voltage-dependent manner, and blockade by these polyamines was described by a three-state kinetic scheme over a wide range of polyamine concentrations. In the context of the scheme, both SPD3+ and SPM4+ access a blocking site with similar effective gating valences (0.84 ± 0.03 e0 for SPD3+ and 0.99 ± 0.04 e0 for SPM4+), whereas SPM4+ binds in the blocked state with an ∼20-fold higher affinity than SPD3+ (Kd = 28.1 ± 3.1 µM for SPD3+ and 1.28 ± 0.20 µM for SPM4+), consistent with a free energy difference of 1.8 kcal/mol. Molecular simulations of the MthK pore in complex with either SPD3+ or SPM4+ are consistent with the leading amine interacting with the hydroxyl groups of T59, at the selectivity filter threshold, with access to this site governed by outward movement of K+ ions. These coupled movements can account for a large fraction of the voltage dependence of blockade. In contrast, differences in binding energetics between SPD3+ and SPM4+ may arise from distinct electrostatic interactions between the polyamines and carboxylate oxygens on the side chains of E92 and E96, located in the pore-lining helix.