Overcoming the Heuristic Nature of k-Means Clustering: Identification and Characterization of Binding Modes from Simulations of Molecular Recognition Complexes

Reconnaissance of Allostery via the Restoration of Native p53 DNA-Binding Domain Dynamics in Y220C Mutant p53 Tumor Suppressor Protein

Allosteric regulation of protein dynamics infers a long-range deliberate propagation of information via micro- and macroscale interactions. The Y220C structural mutant is one of the most frequent cancerous p53 mutants. The mutation is distally located from the DNA-binding site of the p53 DNA-binding domain yet causes changes in DNA recognition. This system presents a unique opportunity to examine the allosteric control of mutated proteins under a drug design paradigm. We focus on the key case study of p53 Y220C mutation restoration by a series of new compounds suggested to have Y220C reactivation properties in comparison to our previous findings on the restorative potential of PK11000, a compound studied extensively for reactivation in vitro and in vivo. Previously, we implemented all-atom molecular dynamics (MD) simulations and our lab's techniques of MD-Sectors and MD-Markov state models on the wild type, the Y220C mutant, and Y220C with PK11000 to characterize the effector's restorative properties in terms of conformational dynamics and hydrogen bonding. In this study, we turn to probing the effects made by docking the battery of a new but less well-tested set of aminobenzothiazole derivative compounds reported by Baud et al., which show promise of Y220C rescue. We find that while complete and precise reconstitution of p53 WT molecular dynamics may not be observed as was the case with PK11000, dispersed local reconstitution of loop dynamics provides evidence of rescuing effects by aminobenzothiazole derivative N,2-dihydroxy-3,5-diiodo-4-(1H-pyrrol-1-yl)benzamide, Effector 22, like what we observed for PK11000. Generalizable insights into the mutation and allosteric reactivation of p53 by various effectors by reconstitution of WT dynamics observed in statistical conformational ensemble analysis and network inference are discussed, considering the development of allosteric drug design rooted in first principles.

k-Means NANI: an improved clustering algorithm for Molecular Dynamics simulations

One of the key challenges of k-means clustering is the seed selection or the initial centroid estimation since the clustering result depends heavily on this choice. Alternatives such as k-means++ have mitigated this limitation by estimating the centroids using an empirical probability distribution. However, with high-dimensional and complex datasets such as those obtained from molecular simulation, k-means++ fails to partition the data in an optimal manner. Furthermore, stochastic elements in all flavors of k-means++ will lead to a lack of reproducibility. K-means N-Ary Natural Initiation (NANI) is presented as an alternative to tackle this challenge by using efficient n-ary comparisons to both identify high-density regions in the data and select a diverse set of initial conformations. Centroids generated from NANI are not only representative of the data and different from one another, helping k-means to partition the data accurately, but also deterministic, providing consistent cluster populations across replicates. From peptide and protein folding molecular simulations, NANI was able to create compact and well-separated clusters as well as accurately find the metastable states that agree with the literature. NANI can cluster diverse datasets and be used as a standalone tool or as part of our MDANCE clustering package.

Insights from incorporating quantum computing into drug design workflows

MOTIVATION

While many quantum computing (QC) methods promise theoretical advantages over classical counterparts, quantum hardware remains limited. Exploiting near-term QC in computer-aided drug design (CADD) thus requires judicious partitioning between classical and quantum calculations.

RESULTS

We present HypaCADD, a hybrid classical-quantum workflow for finding ligands binding to proteins, while accounting for genetic mutations. We explicitly identify modules of our drug-design workflow currently amenable to replacement by QC: non-intuitively, we identify the mutation-impact predictor as the best candidate. HypaCADD thus combines classical docking and molecular dynamics with quantum machine learning (QML) to infer the impact of mutations. We present a case study with the coronavirus (SARS-CoV-2) protease and associated mutants. We map a classical machine-learning module onto QC, using a neural network constructed from qubit-rotation gates. We have implemented this in simulation and on two commercial quantum computers. We find that the QML models can perform on par with, if not better than, classical baselines. In summary, HypaCADD offers a successful strategy for leveraging QC for CADD.

AVAILABILITY AND IMPLEMENTATION

Jupyter Notebooks with Python code are freely available for academic use on GitHub: https://www.github.com/hypahub/hypacadd_notebook.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Agathisflavone, a natural biflavonoid that inhibits SARS-CoV-2 replication by targeting its proteases

Despite the fast development of vaccines, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still circulates through variants of concern (VoC) and escape the humoral immune response. SARS-CoV-2 has provoked over 200,000 deaths/months since its emergence and only a few antiviral drugs showed clinical benefit up to this moment. Thus, chemical structures endowed with anti-SARS-CoV-2 activity are important for continuous antiviral development and natural products represent a fruitful source of substances with biological activity. In the present study, agathisflavone (AGT), a biflavonoid from Anacardium occidentale was investigated as a candidate anti-SARS-CoV-2 compound. In silico and enzymatic analysis indicated that AGT may target mainly the viral main protease (M^pro) and not the papain-like protease (PL^pro) in a non-competitive way. Cell-based assays in type II pneumocytes cell lineage (Calu-3) showed that SARS-CoV-2 is more susceptible to AGT than to apigenin (APG, monomer of AGT), in a dose-dependent manner, with an EC₅₀ of 4.23 ± 0.21 μM and CC₅₀ of 61.3 ± 0.1 μM and with a capacity to inhibit the level of pro-inflammatory mediator tumor necrosis factor-alpha (TNF-α). These results configure AGT as an interesting chemical scaffold for the development of novel semisynthetic antivirals against SARS-CoV-2.

Insights into Rational Design of a New Class of Allosteric Effectors with Molecular Dynamics Markov State Models and Network Theory

The development of drugs to restore protein function has been a major advance facilitated by molecular medicine. Allosteric regulation, a phenomenon widely observed in nature, in which a molecule binds to control a distance active site, holds great promise for regulating proteins, yet how to rationally design such a molecule remains a mystery. Over the past few years, we and others have developed several techniques based on molecular dynamics (MD) simulations: MD-Markov state models to capture global conformational substates, and network theory approach utilizing the interaction energy within the protein to confer local allosteric control. We focus on the key case study of the p53 Y220C mutation restoration by PK11000, a compound experimentally shown to reactivate p53 native function in Y220C mutant present tumors. We gain insights into the mutation and allosteric reactivation of the protein, which we anticipate will be applicable to de novo design to engineer new compounds not only for this mutation, but in other macromolecular systems as well.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.