Conformational adaptation in drug-target interactions and residence time

MicroED Structures of Fluticasone Furoate and Fluticasone Propionate Provide New Insights into Their Function

The detailed understanding of the conformational pathway of fluticasone, a widely prescribed medicine for allergic rhinitis, asthma, and chronic obstructive pulmonary disease (COPD), from formulation to its protein-bound state, has been limited due to a lack of access to its high-resolution structures. The three-dimensional (3D) structure of fluticasone furoate 1 remains unpublished, and the deposited structure of fluticasone propionate 2 could be further refined due to refinement against new data. We applied microcrystal electron diffraction (MicroED) to determine the 3D structures of 1 and 2 in their solid states. The preferred geometries in solution were predicted by using density functional theory (DFT) calculations. A comparative analysis of the structures of 1 and 2 across three states (in solid state, in solution, and protein-bound conformation) revealed the course of the conformational changes during the entire transition. Potential energy plots were calculated for the most dynamic bonds, uncovering their rotational barriers. This study underscores the combined use of MicroED and DFT calculations to provide a comprehensive understanding of conformational and energy changes during drug administration. The quantitative comparison also highlights the subtle structural differences that may lead to significant changes in the pharmaceutical properties.

Programming the Kinetics of Chemical Communication: Induced Fit vs Conformational Selection

Life on Earth depends on chemical communication and the ability of biomolecular switches to integrate various chemical signals that trigger their activation or deactivation over time scales ranging from microseconds to days. The ability to similarly program and control the kinetics of artificial switches would greatly assist the design and optimization of future chemical and nanotechnological systems. Two distinct structure-switching mechanisms are typically employed by biomolecular switches: induced fit (IF) and conformational selection (CS). Despite 60 years of experimental and theoretical investigations, the kinetic and evolutive advantages of these two mechanisms remain unclear. Here, we have created a simple modular DNA switch that can operate through both mechanisms and be easily tuned and adapted to characterize its thermodynamic and kinetic parameters. We show that the fastest activation rate of a switch occurs when the ligand is able to bind its inactive conformation (IF). In contrast, we show that when the ligand can only bind the active conformation of the switch (CS), its activation rate can be easily programmed over many orders of magnitude by a simple tuning of its conformational equilibrium. We demonstrate the programming ability of both these mechanisms by designing a drug delivery vessel that can be programmed to release a drug over different time scales (>1000-fold). Overall, these findings provide a programmable strategy to optimize the kinetics of molecular systems and nanomachines while also illustrating how evolution may have taken advantage of IF and CS mechanisms to optimize the kinetics of biomolecular switches.

Binding kinetics of highly mutated HIV-1 subtype C protease inhibition by Lopinavir and Darunavir in the face of altered conformational dynamics

Highly mutated HIV-1 protease (PR) compromises the efficacy of lopinavir (LPV) and darunavir (DRV) used to formulate salvage regimens in HIV/AIDS management. Here, we report the kinetics of inhibition of lopinavir (LPV) and darunavir (DRV) on highly mutated South African HIV-1 subtype C PR obtained from clinical isolates. The wild-type and mutant South African HIV-1 subtype C PR were cloned and purified. Enzyme inhibition assays and fluorescence spectroscopy were utilized to determine the binding kinetics of LPV and DRV with the wild-type and mutant HIV-1 PR variants. Like DRV, the results of this study show that LPV has a mixed-type inhibition mechanism, which indicates the possibility of a second binding site on HIV-1 PR. Both LPV and DRV poorly inhibited the highly mutated HIV-1 PR variants and had a markedly increased dissociation rate cons bound to the mutant variants compared to the wild type. The fast dissociation of these inhibitors translated into a short residence time of the inhibitor bound to the mutant HIV-1 PR variants. Fluorescent spectroscopy showed that the changes in the tertiary structure of the mutant HIV-1 PR variants were associated with a more open conformation. This open conformation was associated with altered conformational dynamics, which may have resulted in the loss of tight binding of LPV and DRV. This study's findings provide insight into the mechanism of resistance to LPV and DRV by highly mutated HIV-1 PR and provide information supporting the use of binding kinetics measurement in understanding HIV-1 PR inhibitor drug resistance evolution.

The effect of novel β-lactam derivatives synthesized from substituted phenethylamines on resistance genes of MRSA isolates

This study focuses on the activity of previously reported imine and β-lactam derivatives against methicillin-resistant Staphylococcus aureus (MRSA) isolates. The presence of mecA and blaZ genes in the isolates was determined, and the minimum inhibitory concentration (MIC) values were determined based on the antibacterial activity against these isolates. Active compounds were selected and their ability to act against resistant isolates in vitro was determined. Concurrently, biochemical (nitrocefin) and molecular (qRT-PCR) tests were used to investigate the ability of the compounds to induce resistance genes in MRSA isolates. The cytotoxicity of the compounds on human dermal fibroblasts (HDF) was investigated. The MIC values of compounds (10) and (12) against MSSA and MRSA isolates were 7.81 and 15.62 μg ml-1, respectively. The most active compounds were identified as (10) and (12), and it was observed that the isolates did not develop resistance to these compounds in vitro. These compounds were found to inhibit β-lactamase, reduce the expression of resistance genes, and exhibit reduced HDF cell toxicity in a dose-dependent manner. According to the findings of the study, it can be concluded that these compounds show promise as hits with an interesting mechanism of action for further chemical modifications to develop new MRSA inhibitors.

Simplified Method for Kinetic and Thermodynamic Screening of Cardiotonic Steroids through the K+-Dependent Phosphatase Activity of Na+/K+-ATPase with Chromogenic pNPP Substrate

The antitumor effect of cardiotonic steroids (CTS) has stimulated the search for new methods to evaluate both kinetic and thermodynamic aspects of their binding to Na+/K+-ATPase (IUBMB Enzyme Nomenclature). We propose a real-time assay based on a chromogenic substrate for phosphatase activity (pNPPase activity), using only two concentrations with an inhibitory progression curve, to obtain the association rate (kon ), dissociation rate (koff ), and equilibrium (Ki ) constants of CTS for the structure-kinetics relationship in drug screening. We show that changing conditions (from ATPase to pNPPase activity) resulted in an increase of Ki of the cardenolides digitoxigenin, essentially due to a reduction of kon In contrast, the Ki of the structurally related bufadienolide bufalin increased much less due to the reduction of its koff partially compensating the decrease of its kon When evaluating the kinetics of 15 natural and semisynthetic CTS, we observed that both kon and koff correlated with Ki (Spearman test), suggesting that differences in potency depend on variations of both kon and koff A rhamnose in C3 of the steroidal nucleus enhanced the inhibitory potency by a reduction of koff rather than an increase of kon Raising the temperature did not alter the koff of digitoxin, generating a ΔH‡ (koff ) of -10.4 ± 4.3 kJ/mol, suggesting a complex dissociation mechanism. Based on a simple and inexpensive methodology, we determined the values of kon , koff , and Ki of the CTS and provided original kinetics and thermodynamics differences between CTS that could help the design of new compounds. SIGNIFICANCE STATEMENT: This study describes a fast, simple, and cost-effective method for the measurement of phosphatase pNPPase activity enabling structure-kinetics relationships of Na+/K+-ATPase inhibitors, which are important compounds due to their antitumor effect and endogenous role. Using 15 compounds, some of them original, this study was able to delineate the kinetics and/or thermodynamics differences due to the type of sugar and lactone ring present in the steroid structure.

MicroED Structures of Fluticasone Furoate and Fluticasone Propionate Provide New Insights to Their Function

The detailed understanding of fluticasone, a widely prescribed medicine for allergic rhinitis, asthma, and chronic obstructive pulmonary disease (COPD), has not been complete due to challenges in structural elucidation. The three-dimensional (3D) structure of fluticasone furoate 1 remained undetermined for decades, while the existing structures of fluticasone propionate 2 required refinement against improved data. In this study, we applied microcrystal electron diffraction (MicroED) to determine the 3D structures of 1 and 2 in their drug formulation state. Density functional theory (DFT) calculations were utilized to model solvent effects to determine the preferred geometries in solution. A comparative analysis of structures of 1 and 2 across three states (drug formulation state, in solution, and biologically active state) revealed major conformational changes during the entire transition. Potential energy plots were calculated for the most dynamic bonds, uncovering their rotational barriers. This study underscores the combined use of MicroED and DFT calculations to provide a comprehensive understanding of conformational and energy changes during drug functioning in humans. The quantitative comparison highlights the subtle structural differences that can lead to significant functional changes in pharmaceutical properties.

Efficient molecular conformation generation with quantum-inspired algorithm

CONTEXT

Conformation generation, also known as molecular unfolding (MU), is a crucial step in structure-based drug design, remaining a challenging combinatorial optimization problem. Quantum annealing (QA) has shown great potential for solving certain combinatorial optimization problems over traditional classical methods such as simulated annealing (SA). However, a recent study showed that a 2000-qubit QA hardware was still unable to outperform SA for the MU problem. Here, we propose the use of quantum-inspired algorithm to solve the MU problem, in order to go beyond traditional SA. We introduce a highly compact phase encoding method which can exponentially reduce the representation space, compared with the previous one-hot encoding method. For benchmarking, we tested this new approach on the public QM9 dataset generated by density functional theory (DFT). The root-mean-square deviation between the conformation determined by our approach and DFT is negligible (less than about 0.5Å), which underpins the validity of our approach. Furthermore, the median time-to-target metric can be reduced by a factor of five compared to SA. Additionally, we demonstrate a simulation experiment by MindQuantum using quantum approximate optimization algorithm (QAOA) to reach optimal results. These results indicate that quantum-inspired algorithms can be applied to solve practical problems even before quantum hardware becomes mature.

METHODS

The objective function of MU is defined as the sum of all internal distances between atoms in the molecule, which is a high-order unconstrained binary optimization (HUBO) problem. The degree of freedom of variables is discretized and encoded with binary variables by the phase encoding method. We employ the quantum-inspired simulated bifurcation algorithm for optimization. The public QM9 dataset is generated by DFT. The simulation experiment of quantum computation is implemented by MindQuantum using QAOA.

Illuminating Neuropeptide Y Y4 Receptor Binding: Fluorescent Cyclic Peptides with Subnanomolar Binding Affinity as Novel Molecular Tools

The neuropeptide Y (NPY) Y4 receptor (Y4R), a member of the family of NPY receptors, is physiologically activated by the linear 36-amino acid peptide pancreatic polypeptide (PP). The Y4R is involved in the regulation of various biological processes, most importantly pancreatic secretion, gastrointestinal motility, and regulation of food intake. So far, Y4R binding affinities have been mostly studied in radiochemical binding assays. Except for a few fluorescently labeled PP derivatives, fluorescence-tagged Y4R ligands with high affinity have not been reported. Here, we introduce differently fluorescence-labeled (Sulfo-Cy5, Cy3B, Py-1, Py-5) Y4R ligands derived from recently reported cyclic hexapeptides showing picomolar Y4R binding affinity. With pKi values of 9.22-9.71 (radioligand competition binding assay), all fluorescent ligands (16-19) showed excellent Y4R affinity. Y4R saturation binding, binding kinetics, and competition binding with reference ligands were studied using different fluorescence-based methods: flow cytometry (Sulfo-Cy5, Cy3B, and Py-1 label), fluorescence anisotropy (Cy3B label), and NanoBRET (Cy3B label) binding assays. These experiments confirmed the high binding affinity to Y4R (equilibrium pKd: 9.02-9.9) and proved the applicability of the probes for fluorescence-based Y4R competition binding studies and imaging techniques such as single-receptor molecule tracking.

De Novo Discovery of Cysteine Frameworks for Developing Multicyclic Peptide Libraries for Ligand Discovery

Conserved cysteine frameworks are essential components of disulfide-rich peptides (DRPs), which dominantly define the structural diversity of both naturally occurring and de novo-designed DRPs. However, there are only very limited numbers of conserved cysteine frameworks, and general methods enabling de novo discovery of cysteine frameworks with robust foldability are still not available. Here, we devised a "touchstone"-based strategy that relies on chasing oxidative foldability between two individual disulfide-rich folds on the phage surface to discover new cysteine frameworks from random sequences. Unique cysteine frameworks with a high degree of compatibility with phage display systems and broad sequence tolerance were successfully identified, which were subsequently exploited for the development of multicyclic DRP libraries, enabling the rapid discovery of new peptide ligands with low-nanomolar and picomolar binding affinity. This study provides an unprecedented method for exploring and exploiting the sequence and structure space of DRPs that is not readily accessible by existing strategies, holding the potential to revolutionize the study of DRPs and significantly advance the design and discovery of multicyclic peptide ligands and drugs.

Entrectinib a Plausible Inhibitor for Osteopontin (SPP1) in Cervical Cancer-Integrated Bioinformatic Approach

Cervical cancer is one of the major causes of death in women, especially in developing countries bearing more than a quarter of the global burden. Secreted phosphoprotein-1, also known as OPN (osteopontin), is an integrin-binding glycophosphoprotein that is overexpressed in a variety of tumors. OPN is a chemokine-like calcified ECM-associated protein that plays a crucial role in evaluating the metastatic potential of various cancers. However, the role of SPP1 in the tumor microenvironment and associated signaling pathways in CC is still unclear. In our study, three CC microarray datasets (GSE9750, GSE46857, and GSE67522) were obtained from the GEO database to identify the differentially expressed genes. Enrichment analysis was carried out by Enrichr and ShinyGO and the PPI interaction network was created by using String and Cytoscape. GEPIA datasets were used to validate the top 10 hub genes, and virtual screening, docking, and dynamic simulation studies were used to identify a suitable inhibitor against the OPN protein using MVD, PyRx, and GROMACS respectively. Our results show that a total of 11 DEGs were common for three datasets and gene ontology pathway enrichment analysis revealed that 2 biological processes i.e. programmed cell death and animal organ development commonly affected mechanisms in all three datasets. Docking and dynamic studies revealed that Entrectinib showed excellent binding affinity against OPN protein. Based on the results, we conclude that OPN is one of the most upregulated genes in cervical cancer and Entrectinib emerges to be a promising potential OPN inhibitor to curtail cervical cancer progression. Schematic representation: The schematic representation of methodology steps is illustrated in the graphical abstract. Schematic representation of methodology.

QMrebind: incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations

Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparameterization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.

How Robust Is the Ligand Binding Transition State?

For many drug targets, it has been shown that the kinetics of drug binding (e.g., on rate and off rate) is more predictive of drug efficacy than thermodynamic quantities alone. This motivates the development of predictive computational models that can be used to optimize compounds on the basis of their kinetics. The structural details underpinning these computational models are found not only in the bound state but also in the short-lived ligand binding transition states. Although transition states cannot be directly observed experimentally due to their extremely short lifetimes, recent successes have demonstrated that modeling the ligand binding transition state is possible with the help of enhanced sampling molecular dynamics methods. Previously, we generated unbinding paths for an inhibitor of soluble epoxide hydrolase (sEH) with a residence time of 11 min. Here, we computationally modeled unbinding events with the weighted ensemble method REVO (resampling of ensembles by variation optimization) for five additional inhibitors of sEH with residence times ranging from 14.25 to 31.75 min, with average prediction accuracy within an order of magnitude. The unbinding ensembles are analyzed in detail, focusing on features of the ligand binding transition state ensembles (TSEs). We find that ligands with similar bound poses can show significant differences in their ligand binding TSEs, in terms of their spatial distribution and protein-ligand interactions. However, we also find similarities across the TSEs when examining more general features such as ligand degrees of freedom. Together these findings show significant challenges for rational, kinetics-based drug design.

Simulation of the ligand-leaving process of the human heat shock protein

The human heat shock protein plays a critical role in various diseases and is an important target for pharmacological modulation. Simulation of conformational changes and free energy profiles of the human heat shock protein derived by the ligand-leaving process is a challenging issue. In this work, steered molecular dynamics simulation was adopted to simulate the ligand-leaving process. Two composite systems of heat shock protein NHSP90 and small molecules 6FJ and 6G7 are selected as research objects. The free energy during the leaving of ligand small molecules is calculated using conventional molecular dynamics simulation, steered molecular dynamics simulation (SMD), and the umbrella sampling method. We found that the a slower pulling velocity (0.001 nm ns-1) will result in 2.19 kcal mol-1, and the umbrella sampling method gives a value of 3.26 kcal mol-1 for the free energy difference for the two systems, which reasonably agrees with experimental results. A faster-pulling velocity (0.01 nm ns-1) leads to a large overestimation of free energy. At the same time, the conformational analysis indicated that the faster pulling velocity may lead to the conformational change of NHSP90, which was proved to be false by the slower pulling velocity and the umbrella sampling method.

[3H]UR-JG102-A Radiolabeled Cyclic Peptide with High Affinity and Excellent Selectivity for the Neuropeptide Y Y4 Receptor

The family of human neuropeptide Y receptors (YRs) comprises four subtypes (Y1R, Y2R, Y4R, and Y5R) that are involved in the regulation of numerous physiological processes. Until now, Y4R binding studies have been predominantly performed in hypotonic sodium-free buffers using 125I-labeled derivatives of the endogenous YR agonists pancreatic polypeptide or peptide YY. A few tritium-labeled Y4R ligands have been reported; however, when used in buffers containing sodium at a physiological concentration, their Y4R affinities are insufficient. Based on the cyclic hexapeptide UR-AK86C, we developed a new tritium-labeled Y4R radioligand ([3H]UR-JG102, [3H]20). In sodium-free buffer, [3H]20 exhibits a very low Y4R dissociation constant (Kd 0.012 nM). In sodium-containing buffer (137 mM Na+), the Y4R affinity is lower (Kd 0.11 nM) but still considerably higher compared to previously reported tritiated Y4R ligands. Therefore, [3H]20 represents a useful tool compound for the determination of Y4R binding affinities under physiological-like conditions.

Evolution of Slow-Binding Inhibitors Targeting Histone Deacetylase Isoforms

Because the overexpression of histone deacetylase enzymes (HDACs) has been linked to numerous diseases, including various cancers and neurodegenerative disorders, HDAC inhibitors have emerged as promising therapeutic agents. However, most HDAC inhibitors lack both subclass and isoform selectivity, which leads to potential toxicity. Unlike classical hydroxamate HDAC inhibitors, slow-binding HDAC inhibitors form tight and prolonged bonds with HDAC enzymes. This distinct mechanism of action improves both selectivity and toxicity profiles, which makes slow-binding HDAC inhibitors a promising class of therapeutic agents for various diseases. Therefore, the development of slow-binding HDAC inhibitors that can effectively target a wide range of HDAC isoforms is crucial. This Perspective provides valuable insights into the potential and progress of slow-binding HDAC inhibitors as promising drug candidates for the treatment of various diseases.

Estimation of Drug-Target Residence Time by Targeted Molecular Dynamics Simulations

Drug-target residence time has emerged as a key selection factor in drug discovery since the binding duration of a drug molecule to its protein target can significantly impact its in vivo efficacy. The challenge in studying the residence time, in early drug discovery stages, lies in how to cost-effectively determine the residence time for the systematic assessment of compounds. Currently, there is still a lack of computational protocols to quickly estimate such a measure, particularly for large and flexible protein targets and drugs. Here, we report an efficient computational protocol, based on targeted molecular dynamics, to rank drug candidates by their residence time and to obtain insights into ligand-target dissociation mechanisms. The method was assessed on a dataset of 10 arylpyrazole inhibitors of CDK8, a large, flexible, and clinically important target, for which the experimental residence time of the inhibitors ranges from minutes to hours. The compounds were correctly ranked according to their estimated residence time scores compared to their experimental values. The analysis of protein-ligand interactions along the dissociation trajectories highlighted the favorable contribution of hydrophobic contacts to residence time and revealed key residues that strongly affect compound residence time.

Click Chemistry-Enabled Conjugation Strategy for Producing Dibenzodiazepinone-Type Fluorescent Probes To Target M2 Acetylcholine Receptors

The development of fluorescently labeled receptor-targeting compounds represents a powerful pharmacological tool to study and characterize ligand-receptor interactions. Despite significant advances in developing sub-type-specific antagonists for muscarinic acetylcholine receptors (mAChRs), reports on antagonists feasible for click chemistry are less common. Here, we designed and synthesized an antagonist suitable for probe attachment through click chemistry, namely, dibenzodiazepinone (DIBA)-alkyne, based on a previously reported DIBA scaffold with a high binding affinity to type-2 mAChR (M2R). To demonstrate the versatility of DIBA-alkyne as a building block for bioconjugates, we assembled DIBA-alkyne with Cyanine5 fluorophores (Cy5) and polyethylene glycol (PEG) biomolecules to obtain fluorescent DIBA antagonist (DIBA-Cy5) and fluorescent DIBA PEG derivatives. Flow cytometric analysis showed that DIBA-Cy5 possessed a high binding affinity to M2R (Kd = 1.80 nM), a two-order magnitude higher binding affinity than M1R. Fluorescent DIBA PEG derivatives maintained a potent binding to the M2R (Kd ≤ 4 nM), confirmed by confocal microscopic imaging. Additionally, DIBA-Cy5 can serve as a fluorescent ligand in the receptor-ligand competitive binding assay for other mAChR ligands, an attractive alternative to the traditional radioligand-based assay. The competitive binding mode between DIBA-Cy5 and orthosteric antagonist atropine/allosteric modulator LY2119620 indicated a dualsteric binding mode of the DIBA-type antagonist to M2R. Lastly, we demonstrated the direct staining of DIBA-Cy5 to M2R receptors in the sinoatrial node of a mouse heart. The adaptability of the clickable DIBA antagonist to a wide range of fluorophores and biomolecules can facilitate its use in various biomedical applications such as binding assays that screen compounds for M2R as the receptor target.

Natural Products as Mcl-1 Inhibitors: A Comparative Study of Experimental and Computational Modelling Data

The human myeloid leukemia cell differentiation protein (hMcl-1) is an anti-apoptotic multi-partner protein, belonging to the B-cell lymphoma-2 (Bcl-2) family of proteins. Studies have linked hMcl-1 alleviated expression with resistance to hemopoietic chemotherapeutics, which makes it a key drug target in blood cancers. However, most of the developed small- to medium-sized hMcl-1 inhibitors have typical off-target activity towards other members of the Bcl-2 family. To improve the hMcl-1 inhibitor design, especially exploring a suitable scaffold with pharmacophoric features, we focused on natural hMcl-1 inhibitors. To date, seven classes of natural compounds have been isolated, which display a low micromolar affinity for hMcl-1 and have limited biophysical studies. We screened hMcl-1 co-crystal structures, and identified nine co-crystal structures of hMcl-1 protein, which were later evaluated by multiple receptor conformations (which indicates that the differences between hMcl-1 in crystal structures are low (RMSD values between 0.52 and 1.13 Å, average RMSD of 0.638–0.888 Å, with a standard deviation of 0.102–0.185Å)), and multiple ligand conformations (which led to the selection of the PDB structure, 3WIX (RMSD value = 0.879 Å, standard deviation 0.116 Å), to accommodate various Mcl-1 ligands from a range of co-crystal PDB files) methods. Later, the three adopted docking methods were assessed for their ability to reproduce the conformation bound to the crystal as well as predict trends in Ki values based on calculated RMSD and docking energies. Iterative docking and clustering of the docked pose within ≤1.0 Å was used to evaluate the reproducibility of the adopted docking methods and compared with their experimentally determined hMcl-1 affinity data.

Large-Scale Ligand Perturbations of the Protein Conformational Landscape Reveal State-Specific Interaction Hotspots

Protein flexibility is important for ligand binding but often ignored in drug design. Considering proteins as ensembles rather than static snapshots creates opportunities to target dynamic proteins that lack FDA-approved drugs, such as the human chaperone, heat shock protein 90 (Hsp90). Hsp90α accommodates ligands with a dynamic lid domain, yet no comprehensive analysis relating lid conformations to ligand properties is available. To date, ∼300 ligand-bound Hsp90α crystal structures are deposited in the Protein Data Bank, which enables us to consider ligand binding as a perturbation of the protein conformational landscape. By estimating binding site volumes, we classified structures into distinct major and minor lid conformations. Supported by retrospective docking, each conformation creates unique hotspots that bind chemically distinguishable ligands. Clustering revealed insightful exceptions and the impact of crystal packing. Overall, Hsp90α’s plasticity provides a cautionary tale of overinterpreting individual crystal structures and motivates an ensemble-based view of drug design.

Targeting Glutaminase by Natural Compounds: Structure-Based Virtual Screening and Molecular Dynamics Simulation Approach to Suppress Cancer Progression

Cancer cells change their glucose and glutamine (GLU) metabolism to obtain the energy required to continue growing. Glutaminase (GLS) plays a crucial role in promoting cell metabolism for cancer cell growth; targeting GLU metabolism by inhibiting GLS has attracted interest as a potential cancer management strategy. Herein, we employed a sequential screening of traditional Chinese medicine (TCM) database followed by drug-likeness and molecular dynamics simulations against the active site of GLS. We report 12 potent compounds after screening the TCM database against GLS, followed by a drug-likeness filter with Lipinski and Veber rule criteria. Among them, ZINC03978829 and ZINC32296657 were found to have higher binding energy (BE) values than the control compound 6-Diazo-5-Oxo-L-Norleucine, with BEs of -9.3 and -9.7 kcal/mol, respectively, compared to the BE of 6-Diazo-5-Oxo-L-Norleucine (-4.7 kcal/mol) with GLS. Molecular dynamics simulations were used to evaluate the results further, and a 100 ns MD simulation revealed that the hits form stable complexes with GLS and formed 2-5 hydrogen bond interactions. This study indicates that these hits might be employed as GLS inhibitors in the battle against cancer. However, more laboratory tests are a prerequisite to optimize them as GLS inhibitors.

A Target Engagement Assay to Assess Uptake, Potency, and Retention of Antibiotics in Living Bacteria

New antibiotics are urgently needed to counter the emergence of antimicrobial-resistant pathogenic bacteria. A major challenge in antibiotic drug discovery is to turn potent biochemical inhibitors of essential bacterial components into effective antimicrobials. This difficulty is underpinned by a lack of methods to investigate the physicochemical properties needed for candidate antibiotics to permeate the bacterial cell envelope and avoid clearance by the action of bacterial efflux pumps. To address these issues, here we used a target engagement assay to measure the equilibrium and kinetic binding parameters of antibiotics targeting dihydrofolate reductase (DHFR) in live bacteria. We also used this assay to identify novel DHFR ligands having antimicrobial activity. We validated this approach using the Gram-negative bacteria Escherichia coli and the emerging human pathogen Mycobacterium abscessus. We expect the use of target engagement assays in bacteria to expedite the discovery and progression of novel, cell-permeable antibiotics with on-target activity.

Development of an intracellular quantitative assay to measure compound binding kinetics

Contemporary drug discovery typically quantifies the effect of a molecule on a biological target using the equilibrium-derived measurements of IC50, EC50, or KD. Kinetic descriptors of drug binding are frequently linked with the effectiveness of a molecule in modulating a disease phenotype; however, these parameters are yet to be fully adopted in early drug discovery. Nanoluciferase bioluminescence resonance energy transfer (NanoBRET) can be used to measure interactions between fluorophore-conjugated probes and luciferase fused target proteins. Here, we describe an intracellular NanoBRET competition assay that can be used to quantify cellular kinetic rates of compound binding to nanoluciferase-fused bromodomain and extra-terminal (BET) proteins. Comparative rates are generated using a cell-free NanoBRET assay and by utilizing orthogonal recombinant protein-based methodologies. A screen of known pan-BET inhibitors is used to demonstrate the value of this approach in the investigation of kinetic selectivity between closely related proteins.

Enhanced Sampling Approach to the Induced-Fit Docking Problem in Protein-Ligand Binding: The Case of Mono-ADP-Ribosylation Hydrolase Inhibitors

Enhanced sampling methods can predict free-energy landscapes associated with protein/ligand binding, characterizing the involved intermolecular interactions in a precise way. However, these in silico approaches can be challenged by induced-fit effects. Here, we present a variant of volume-based metadynamics tailored to tackle this problem in a general and efficient way. The validity of the approach is established by applying it to substrate/enzyme complexes of pharmacological relevance: mono-ADP-ribose (ADPr) in complex with mono-ADP-ribosylation hydrolases (MacroD1 and MacroD2), where induced-fit phenomena are known to be significant. The calculated binding free energies are consistent with experiments, with an absolute error smaller than 0.5 kcal/mol. Our simulations reveal that in all circumstances, the active loops, delimiting the boundaries of the binding site, undergo significant conformation rearrangements upon ligand binding. The calculations further provide, for the first time, the molecular basis of ADPr specificity and the relative changes in its experimental binding affinity on passing from MacroD1 to MacroD2 and all its mutants. Our study paves the way to the quantitative description of induced-fit events in molecular recognition.

Screening and Identification of Potential iNOS Inhibitors to Curtail Cervical Cancer Progression: an In Silico Drug Repurposing Approach

Cervical cancer is the second most common cause of cancer deaths in women worldwide and remains the main reason of mortality among women of reproductive age in developing countries. Nitric oxide is involved in several physiological functions inclusive of inflammatory and immune responses. However, the function of NO in tumor biology is debatable. The inducible NOS (iNOS/NOS2) isoform is the one responsible to maintain the levels of NO, and it exhibits pleotropic effects in various cancers with concentration-dependent pro- and anti-tumor effects. iNOS triggers angiogenesis and endothelial cell migration in tumors by regulating the levels of vascular endothelial growth factor (VEGF). In drug discovery, drug repurposing involves investigations of approved drug candidates to treat various other diseases. In this study, we used anti-cancer drugs and small molecules to target iNOS and identify a potential selective iNOS inhibitor. The structures of ligands were geometrically optimized and energy minimized using Hyperchem software. Molecular docking was performed using Molegro virtual docker, and ligands were selected based on MolDock score, Rerank score, and H-bonding energy. In the study shown, venetoclax compound demonstrated excellent binding affinity to iNOS protein. This compound exhibited the lowest MolDock score and Rerank score with better H-bonding energy to iNOS. The binding efficacy of venetoclax was analyzed by performing molecular docking and molecular dynamic simulations. Multiple parameters were used to analyze the simulation trajectory, like root mean square deviation (RMSD), radius of gyration (Rg), and hydrogen bond interactions. Based on the results, venetoclax emerges to be a promising potential iNOS inhibitor to curtail cervical cancer progression.

Induced fit versus conformational selection: From rate constants to fluxes… and back to rate constants

Induced fit- (IF) and conformational selection (CS) binding mechanisms have long been regarded to be mutually exclusive. Yet, they are now increasingly considered to produce the final ligand-target complex alongside within a thermodynamic cycle. This viewpoint benefited from the introduction of binding fluxes as a tool for analyzing the overall behavior of such cycle. This study aims to provide more vivid and applicable insights into this emerging field. In this respect, combining differential equation- based simulations and hitherto little explored alternative modes of calculation provide concordant information about the intricate workings of such cycle. In line with previous reports, we observe that the relative contribution of IF increases with the ligand concentration at equilibrium. Yet the baseline contribution may vary from one case to another and simulations as well as calculations show that this parameter is essentially regulated by the dissociation rate of both pathways. Closer attention should be paid to how the contributions of IF and CS compare at physiologically relevant drug/ligand concentrations. To this end, a simple equation discloses how changing a limited set of "microscopic" rate constants can extend the concentration range at which CS contributes most effectively. Finally, it could also be beneficial to extend the utilization of flux- based approaches to more physiologically relevant time scales and alternative binding models.

Computational Screening of Natural Compounds for Identification of Potential Anti-Cancer Agents Targeting MCM7 Protein

Minichromosome maintenance complex component 7 (MCM7) is involved in replicative licensing and the synthesis of DNA, and its overexpression is a fascinating biomarker for various cancer types. There is currently no effective agent that can prevent the development of cancer caused by the MCM7 protein. However, on the molecular level, inhibiting MCM7 lowers cancer-related cellular growth. With this purpose, this study screened 452 biogenic compounds extracted from the UEFS Natural Products dataset against MCM protein by using the in silico art of technique. The hit compounds UEFS99, UEFS137, and UEFS428 showed good binding with the MCM7 protein with binding energy values of −9.95, −8.92, and −8.71 kcal/mol, which was comparatively higher than that of the control compound ciprofloxacin (−6.50). The hit (UEFS99) with the minimum binding energy was picked for molecular dynamics (MD) simulation investigation, and it demonstrated stability at 30 ns. Computational prediction of physicochemical property evaluation revealed that these hits are non-toxic and have good drug-likeness features. It is suggested that hit compounds UEFS99, UEFS137, and UEFS428 pave the way for further bench work validation in novel inhibitor development against MCM7 to fight the cancers.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gα_q protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gα_q protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gα_q proteins in cells in which the native Gα_q proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the "lipophilic anchors" of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

Potent Phenylpyridine and Oxodihydrofuran Inhibitors of Cyclooxygenase-2: Optimization toward a Long Residence Time with Balanced Internal Energetics

Long residence time enzyme inhibitors with a two-step binding mechanism are characterized by a high internal energy barrier for target association. This raises the question of whether optimizing residence time via further increasing this internal energy barrier would inevitably lead to insufficient target occupancy in vivo due to slow, time-dependent binding. We attempted to address this question during optimization of cyclooxygenase-2 (COX-2) inhibitors. Defining long residence time drugs with acceptable association and dissociation rate constants required for sufficient target occupancy and sustained efficacy, which we termed "balanced internal energetics", provides an important criterion for successful progression during lead optimization. Despite the advancement of several COX-2 inhibitors to marketed drugs, their detailed inhibition kinetics have been surprisingly limiting especially during the structure-activity relationship process mainly due to the lack of robust kinetic assays. Herein, we describe a reoptimized COX enzymatic assay and a novel MS-based assay enabling detailed mechanistic studies for identifying long residence time COX-2 inhibitors with balanced internal energetics. These efforts led to the discovery of promising leads possessing dissociation half-lives of ≤40 h, much greater than the values of 6 and 0.71 h for two marketed drugs, etoricoxib and celecoxib, respectively. Importantly, the inhibition rate constants remain comparable to those of the marketed drugs and above the lower limits set by the criteria of balanced internal energetics, predicting sufficient target occupancy required for efficacy. Taken together, this study demonstrates the feasibility of increasing the internal energy barrier as a viable approach for lead optimization toward discovering long residence time drug candidates.

High Throughput Virtual Screening and Molecular Dynamics Simulation for Identifying a Putative Inhibitor of Bacterial CTX-M-15

Background: Multidrug resistant bacteria are a major therapeutic challenge. CTX-M-type enzymes are an important group of class A extended-spectrum β-lactamases (ESBLs). ESBLs are the enzymes that arm bacterial pathogens with drug resistance to an array of antibiotics, notably the advanced-generation cephalosporins. The current need for an effective CTX-M-inhibitor is high. Objective: The aim of the current study was to identify a promising anti-CTX-M-15 ligand whose chemical skeleton could be used as a ‘seed-molecule’ for future drug design against resistant bacteria. Methods: Virtual screening of 5,000,000 test molecules was performed by ‘MCULE Drug Discovery Platform’. ‘ADME analyses’ was performed by ‘SWISS ADME’. TOXICITY CHECKER of MCULE was employed to predict the safety profile of the test molecules. The complex of the ‘Top inhibitor’ with the ‘bacterial CTX-M-15 enzyme’ was subjected to 102.25 ns molecular dynamics simulation. This simulation was run for 3 days on a HP ZR30w workstation. Trajectory analyses were performed by employing the macro ‘md_analyze.mcr’ of YASARA STRUCTURE version 20.12.24.W.64 using AMBER14 force field. YANACONDA macro language was used for complex tasks. Figures, including RMSD and RMSF plots, were generated. Snapshots were acquired after every 250 ps. Finally, two short videos of ‘41 s’ and ‘1 min and 22 s’ duration were recorded. Results: 5-Amino-1-(2H-[1,2,4]triazino[5,6-b]indol-3-yl)-1H-pyrazole-4-carbonitrile, denoted by the MCULE-1352214421-0-56, displayed the most efficient binding with bacterial CTX-M-15 enzyme. This screened molecule significantly interacted with CTX-M-15 via 13 amino acid residues. Notably, nine amino acid residues were found common to avibactam binding (the reference ligand). Trajectory analysis yielded 410 snapshots. The RMSD plot revealed that around 26 ns, equilibrium was achieved and, thereafter, the complex remained reasonably stable. After a duration of 26 ns and onwards until 102.25 ns, the backbone RMSD fluctuations were found to be confined within a range of 0.8–1.4 Å. Conclusion: 5-Amino-1-(2H-[1,2,4]triazino[5,6-b]indol-3-yl)-1H-pyrazole-4-carbonitrile could emerge as a promising seed molecule for CTX-M-15-inhibitor design. It satisfied ADMET features and displayed encouraging ‘simulation results’. Advanced plots obtained by trajectory analyses predicted the stability of the proposed protein-ligand complex. ‘Hands on’ wet laboratory validation is warranted.

Identification of New Mycobacterium tuberculosis Proteasome Inhibitors Using a Knowledge-Based Computational Screening Approach

Mycobacterium tuberculosis (Mtb) is a deadly tuberculosis (TB)-causing pathogen. The proteasome is vital to the survival of Mtb and is therefore validated as a potential target for anti-TB therapy. Mtb resistance to existing antibacterial agents has enhanced drastically, becoming a worldwide health issue. Therefore, new potential therapeutic agents need to be developed that can overcome the complications of TB. With this purpose, in the present study, 224,205 natural compounds from the ZINC database have been screened against the catalytic site of Mtb proteasome by the computational approach. The best scoring hits, ZINC3875469, ZINC4076131, and ZINC1883067, demonstrated robust interaction with Mtb proteasome with binding energy values of −7.19, −7.95, and −7.21 kcal/mol for the monomer (K-chain) and −8.05, −9.10, and −7.07 kcal/mol for the dimer (both K and L chains) of the beta subunit, which is relatively higher than that of reference compound HT1171 (−5.83 kcal/mol (monomer) and −5.97 kcal/mol (dimer)). In-depth molecular docking of top-scoring compounds with Mtb proteasome reveals that amino acid residues Thr1, Arg19, Ser20, Thr21, Gln22, Gly23, Asn24, Lys33, Gly47, Asp124, Ala126, Trp129, and Ala180 are crucial in binding. Furthermore, a molecular dynamics study showed steady-state interaction of hit compounds with Mtb proteasome. Computational prediction of physicochemical property assessment showed that these hits are non-toxic and possess good drug-likeness properties. This study proposed that these compounds could be utilized as potential inhibitors of Mtb proteasome to combat TB infection. However, there is a need for further bench work experiments for their validation as inhibitors of Mtb proteasome.

Red-Edge Excitation Shift Spectroscopy (REES): Application to Hidden Bound States of Ligands in Protein-Ligand Complexes

Ligand-protein binding is responsible for the vast majority of bio-molecular functions. Most experimental techniques examine the most populated ligand-bound state. The determination of less populated, intermediate, and transient bound states is experimentally challenging. However, hidden bound states are also important because these can strongly influence ligand binding and unbinding processes. Here, we explored the use of a classical optical spectroscopic technique, red-edge excitation shift spectroscopy (REES) to determine the number, population, and energetics associated with ligand-bound states in protein–ligand complexes. We describe a statistical mechanical model of a two-level fluorescent ligand located amongst a finite number of discrete protein microstates. We relate the progressive emission red shift with red-edge excitation to thermodynamic parameters underlying the protein–ligand free energy landscape and to photo-physical parameters relating to the fluorescent ligand. We applied the theoretical model to published red-edge excitation shift data from small molecule inhibitor–kinase complexes. The derived thermodynamic parameters allowed dissection of the energetic contribution of intermediate bound states to inhibitor–kinase interactions.

Medicinal chemistry strategies for discovering antivirals effective against drug-resistant viruses

During the last forty years we have witnessed impressive advances in the field of antiviral drug discovery culminating with the introduction of therapies able to stop human immunodeficiency virus (HIV) replication, or cure hepatitis C virus infections in people suffering from liver disease. However, there are important viral diseases without effective treatments, and the emergence of drug resistance threatens the efficacy of successful therapies used today. In this review, we discuss strategies to discover antiviral compounds specifically designed to combat drug resistance. Currently, efforts in this field are focused on targeted proteins (e.g. multi-target drug design strategies), but also on drug conformation (either improving drug positioning in the binding pocket or introducing conformational constraints), in the introduction or exploitation of new binding sites, or in strengthening interaction forces through the introduction of multiple hydrogen bonds, covalent binding, halogen bonds, additional van der Waals forces or multivalent binding. Among the new developments, proteolysis targeting chimeras (PROTACs) have emerged as a valid approach taking advantage of intracellular mechanisms involving protein degradation by the ubiquitin-proteasome system. Finally, several molecules targeting host factors (e.g. human dihydroorotate dehydrogenase and DEAD-box polypeptide 3) have been identified as broad-spectrum antiviral compounds. Implementation of herein described medicinal chemistry strategies are expected to contribute to the discovery of new drugs effective against current and future threats due to emerging and re-emerging viral pandemics.

High-Throughput Quantitative Assay Technologies for Accelerating the Discovery and Optimization of Targeted Protein Degradation Therapeutics

The aberrant regulation of protein expression and function can drastically alter cellular physiology and lead to numerous pathophysiological conditions such as cancer, inflammatory diseases, and neurodegeneration. The steady-state expression levels of endogenous proteins are controlled by a balance of de novo synthesis rates and degradation rates. Moreover, the levels of activated proteins in signaling cascades can be further modulated by a variety of posttranslational modifications and protein-protein interactions. The field of targeted protein degradation is an emerging area for drug discovery in which small molecules are used to recruit E3 ubiquitin ligases to catalyze the ubiquitination and subsequent degradation of disease-causing target proteins by the proteasome in both a dose- and time-dependent manner. Traditional approaches for quantifying protein level changes in cells, such as Western blots, are typically low throughput with limited quantification, making it hard to drive the rapid development of therapeutics that induce selective, rapid, and sustained protein degradation. In the last decade, a number of techniques and technologies have emerged that have helped to accelerate targeted protein degradation drug discovery efforts, including the use of fluorescent protein fusions and reporter tags, flow cytometry, time-resolved fluorescence energy transfer (TR-FRET), and split luciferase systems. Here we discuss the advantages and disadvantages associated with these technologies and their application to the development and optimization of degraders as therapeutics.

An update into the medicinal chemistry of translocator protein (TSPO) ligands

The Translocator Protein 18 kDa (TSPO) has been discovered in 1977 as an alternative binding site for the benzodiazepine diazepam. It is an evolutionary well-conserved and tryptophan-rich 169-amino acids protein with five alpha helical transmembrane domains stretching the outer mitochondrial membrane, with the carboxyl-terminus in the cytosol and a short amino-terminus in the intermembrane space of mitochondrion. At this level, together with the voltage-dependent anion channel (VDAC) and the adenine nucleotide translocase (ANT), it forms the mitochondrial permeability transition pore (MPTP). TSPO expression is ubiquitary, with higher levels in steroid producing tissues; in the central nervous system, it is mainly expressed in glial cells and in neurons. TSPO is implicated in a variety of fundamental cellular processes including steroidogenesis, heme biosynthesis, mitochondrial respiration, mitochondrial membrane potential, cell proliferation and differentiation, cell life/death balance, oxidative stress. Altered TSPO expression has been found in some pathological conditions. In particular, high TSPO expression levels have been documented in cancer, neuroinflammation, and brain injury. Conversely, low TSPO expression levels have been evidenced in anxiety disorders. Therefore, TSPO is not only an interesting drug target for therapeutic purpose (anticonvulsant, anxiolytic, etc.), but also a valid diagnostic marker of related-diseases detectable by fluorescent or radiolabeled ligands. The aim of this report is to present an update of previous reviews dealing with the medicinal chemistry of TSPO and to highlight the most outstanding advances in the development of TSPO ligands as potential therapeutic or diagnostic tools, especially referring to the last five years.

Slow-binding reversible inhibitor of acetylcholinesterase with long-lasting action for prophylaxis of organophosphate poisoning

Organophosphorus (OP) compounds represent a serious health hazard worldwide. The dominant mechanism of their action results from covalent inhibition of acetylcholinesterase (AChE). Standard therapy of acute OP poisoning is partially effective. However, prophylactic administration of reversible or pseudo-irreversible AChE inhibitors before OP exposure increases the efficiency of standard therapy. The purpose of the study was to test the duration of the protective effect of a slow-binding reversible AChE inhibitor (C547) in a mouse model against acute exposure to paraoxon (POX). It was shown that the rate of inhibition of AChE by POX in vitro after pre-inhibition with C547 was several times lower than without C547. Ex vivo pre-incubation of mouse diaphragm with C547 significantly prevented the POX-induced muscle weakness. Then it was shown that pre-treatment of mice with C547 at the dose of 0.01 mg/kg significantly increased survival after poisoning by 2xLD50 POX. The duration of the pre-treatment was effective up to 96 h, whereas currently used drug for pre-exposure treatment, pyridostigmine at a dose of 0.15 mg/kg was effective less than 24 h. Thus, long-lasting slow-binding reversible AChE inhibitors can be considered as new potential drugs to increase the duration of pre-exposure treatment of OP poisoning.

Targeting the SARS-CoV-2 main protease using FDA-approved Isavuconazonium, a P2-P3 α-ketoamide derivative and Pentagastrin: An in-silico drug discovery approach

The SARS-CoV-2 main protease (M^pro) is an attractive target towards discovery of drugs to treat COVID-19 because of its key role in virus replication. The atomic structure of M^pro in complex with an α-ketoamide inhibitor (Lig13b) is available (PDB ID:6Y2G). Using 6Y2G and the prior knowledge that protease inhibitors could eradicate COVID-19, we designed a computational study aimed at identifying FDA-approved drugs that could interact with M^pro. We searched the DrugBank and PubChem for analogs and built a virtual library containing ∼33,000 conformers. Using high-throughput virtual screening and ligand docking, we identified Isavuconazonium, a ketoamide inhibitor (α-KI) and Pentagastrin as the top three molecules (Lig13b as the benchmark) based on docking energy. The ΔG_bind of Lig13b, Isavuconazonium, α-KI, Pentagastrin was −28.1, −45.7, −44.7, −34.8 kcal/mol, respectively. Molecular dynamics simulation revealed that these ligands are stable within the M^pro active site. Binding of these ligands is driven by a variety of non-bonded interaction, including polar bonds, H-bonds, van der Waals and salt bridges. The overall conformational dynamics of the complexed-M^pro was slightly altered relative to apo-M^pro. This study demonstrates that three distinct classes molecules, Isavuconazonium (triazole), α-KI (ketoamide) and Pentagastrin (peptide) could serve as potential drugs to treat patients with COVID-19.

•

Using computational modelling, we show that SARS-CoV-2 main protease (M^pro) interacts with peptidomimetic drugs.

•

The interaction between the M^pro and the peptidomimetics is energetically favourable.

•

Isavuconazonium, a P2–P3 α-ketoamide derivative and Pentagastrin showed the tightest and most favourable binding to M^pro.

Slow-binding inhibitors of acetylcholinesterase of medical interest

Certain ligands slowly bind to acetylcholinesterase. As a result, there is a slow establishment of enzyme-inhibitor equilibrium characterized by a slow onset of inhibition prior reaching steady state. Three mechanisms account for slow-binding inhibition: a) slow binding rate constant k_on, b) slow ligand induced-fit following a fast binding step, c) slow conformational selection of an enzyme form. The slow equilibrium may be followed by a chemical step. This later that can be irreversible has been observed with certain alkylating agents and substrate transition state analogs. Slow-binding inhibitors present long residence times on target. This results in prolonged pharmacological or toxicological action. Through several well-known molecules (e.g. huperzine) and new examples (tocopherol, trifluoroacetophenone and a 6-methyluracil alkylammonium derivative), we show that slow-binding inhibitors of acetylcholinesterase are promising drugs for treatment of neurological diseases such as Alzheimer disease and myasthenia gravis. Moreover, they may be of interest for neuroprotection (prophylaxis) against organophosphorus poisoning. This article is part of the special issue entitled 'Acetylcholinesterase Inhibitors: From Bench to Bedside to Battlefield'.

Red-Emitting Dibenzodiazepinone Derivatives as Fluorescent Dualsteric Probes for the Muscarinic Acetylcholine M2 Receptor

Fluorescently labeled dibenzodiazepinone-type muscarinic acetylcholine receptor (MR) antagonists, including dimeric ligands, were prepared using red-emitting cyanine dyes. Probes containing a fluorophore with negative charge showed high M2R affinities (pKi (radioligand competition binding): 9.10-9.59). Binding studies at M1 and M3-M5 receptors indicated a M2R preference. Flow cytometric and high-content imaging saturation and competition binding (M1R, M2R, and M4R) confirmed occupation of the orthosteric site. Confocal microscopy revealed that fluorescence was located mainly at the cell membrane (CHO-hM2R cells). Results from dissociation and saturation binding experiments (M2R) in the presence of allosteric M2R modulators (dissociation: W84, LY2119620, and alcuronium; saturation binding: W84) were consistent with a competitive mode of action between the fluorescent probes and the allosteric ligands. Taken together, these lines of evidence indicate that these ligands are useful fluorescent molecular tools to label the M2R in imaging and binding studies and suggest that they have a dualsteric mode of action.

KOFFI and Anabel 2.0-a new binding kinetics database and its integration in an open-source binding analysis software

The kinetics of featured interactions (KOFFI) database is a novel tool and resource for binding kinetics data from biomolecular interactions. While binding kinetics data are abundant in literature, finding valuable information is a laborious task. We used text extraction methods to store binding rates (association, dissociation) as well as corresponding meta-information (e.g. methods, devices) in a novel database. To date, over 270 articles were manually curated and binding data on over 1705 interactions was collected and stored in the (KOFFI) database. Moreover, the KOFFI database application programming interface was implemented in Anabel (open-source software for the analysis of binding interactions), enabling users to directly compare their own binding data analyses with related experiments described in the database.

Perspective: Implications of Ligand-Receptor Binding Kinetics for Therapeutic Targeting of G Protein-Coupled Receptors

The concept of ligand-receptor binding kinetics has been broadly applied in drug development pipelines focusing on G protein-coupled receptors (GPCRs). The ligand residence time (RT) for a receptor describes how long a ligand-receptor complex exists, and is defined as the reciprocal of the dissociation rate constant (k _off). RT has turned out to be a valuable parameter for GPCR researchers focusing on drug development as a good predictor of in vivo efficacy. The positive correlation between RT and in vivo efficacy has been established for several drugs targeting class A GPCRs (e.g., the neurokinin-1 receptor (NK₁R), the β₂ adrenergic receptor (β₂AR), and the muscarinic 3 receptor (M₃R)) and for drugs targeting class B1 (e.g., the glucagon-like peptide 1 receptor (GLP-1R)). Recently, the association rate constant (k _on) has gained similar attention as another parameter affecting in vivo efficacy. In the current perspective, we address the importance of studying ligand-receptor binding kinetics for therapeutic targeting of GPCRs, with an emphasis on how binding kinetics can be altered by subtle molecular changes in the ligands and/or the receptors and how such changes affect treatment outcome. Moreover, we speculate on the impact of binding kinetic parameters for functional selectivity and sustained receptor signaling from endosomal compartments; phenomena that have gained increasing interest in attempts to improve therapeutic targeting of GPCRs.

Validation of the Slow Off-Kinetics of Sirtuin-Rearranging Ligands (SirReals) by Means of Label-Free Electrically Switchable Nanolever Technology

We have discovered the sirtuin-rearranging ligands (SirReals) to be highly potent and selective inhibitors of the NAD⁺ -dependent lysine deacetylase Sirt2. Using a biotinylated SirReal in combination with biolayer interferometry, we previously observed a slow dissociation rate of the inhibitor-enzyme complex; this had been postulated to be the key to the high affinity and selectivity of SirReals. However, to attach biotin to the SirReal core, we introduced a triazole as a linking moiety; this was shown by X-ray co-crystallography to interact with Arg97 of the cofactor binding loop. Herein, we aim to elucidate whether the observed long residence time of the SirReals is induced mainly by triazole incorporation or is an inherent characteristic of the SirReal inhibitor core. We used the novel label-free switchSENSE® technology, which is based on electrically switchable DNA nanolevers, to prove that the long residence time of the SirReals is indeed caused by the core scaffold.

[Ligand-receptor binding kinetics in drug design]

Traditionally, the thermodynamic values of affinity are considered as the main criterion for the development of new drugs. Usually, these values for drugs are measured <i>in vitro</i> at steady concentrations of the receptor and ligand, which are differed from <i>in vivo</i> environment. Recent studies have shown that the kinetics of the process of drug binding to its receptor make significant contribution in the drug effectiveness. This has increased attention in characterizing and predicting the rate constants of association and dissociation of the receptor ligand at the stage of preclinical studies of drug candidates. A drug with a long residence time can determine ligand-receptor selectivity (kinetic selectivity), maintain pharmacological activity of the drug at its low concentration in vivo. The paper discusses the theoretical basis of protein-ligand binding, molecular determinants that control the kinetics of the drug-receptor binding. Understanding the molecular features underlying the kinetics of receptor-ligand binding will contribute to the rational design of drugs with desired properties.

A Variable Neighbourhood Descent Heuristic for Conformational Search Using a Quantum Annealer

Discovering the low-energy conformations of a molecule is of great interest to computational chemists, with applications in in silico materials design and drug discovery. In this paper, we propose a variable neighbourhood search heuristic for the conformational search problem. Using the structure of a molecule, neighbourhoods are chosen to allow for the efficient use of a binary quadratic optimizer for conformational search. The method is flexible with respect to the choice of molecular force field and the number of discretization levels in the search space, and can be further generalized to take advantage of higher-order binary polynomial optimizers. It is well-suited for the use of devices such as quantum annealers. After carefully defining neighbourhoods, the method easily adapts to the size and topology of these devices, allowing for seamless scaling alongside their future improvements.

Molecular interaction of investigational ligands with human brain acetylcholinesterase

Alzheimer's disease, a neurodegenerative disorder continues to be an area of investigation by the international researchers' fraternity. Despite all the ongoing efforts, the effective set of promising cholinesterase inhibitors available in the market for patients' use is limited. Furthermore, the currently available drugs could provide only a palliative type of treatment instead of providing a complete cure or foolproof prevention. Hence, design/discovery of fresh drug molecules as acetylcholinesterase (AChE) inhibitors still remains an urgent requirement. The drug discovery platform, MCULE in the "structure-based virtual screening" (SBVS) mode was used for high throughput ligand screening of over five million structures targeted against the AChE catalytic site. A stepwise query was made for the SBVS input. The number of hits was narrowed down in consecutive succession via varied filtration criteria as AutoDock-Vina rankings, MCULE toxicity filtration, exclusion of ligands having less than four H-bond acceptors, filtration by ΔG cutoff, rule-of-five violation and SWISS ADME profiling. This was followed by holistic analysis of all the results, thereby leading to one promising ligand. The screened out drug molecule, MCULE-5872671137-0-1 exhibited a robust interaction with the AChE catalytic site involving 20 amino acid residues, an acceptable binding free energy of -10.2 kcal/mol in addition to a favorable SWISS ADME-profie showing no harmful effects on the human body. It can be carefully stated that the molecule, MCULE-5872671137-0-1, which is chemically (3S)-N-{4-[(4-chlorophenyl)sulfanyl]phenyl}-3-hydroxypyrrolidine-1-carboxamide could function as a significant "seed" ligand for future design of potent AChE inhibitors and/or novel neuro drugs built upon the seed-scaffold.