Development of Potent and Highly Selective Epoxyketone-Based Plasmodium Proteasome Inhibitors.

Here, we present remarkable epoxyketone-based proteasome inhibitors with low nanomolar in vitro potency for blood-stage Plasmodium falciparum and low cytotoxicity for human cells. Our best compound has more than 2,000-fold greater selectivity for erythrocytic-stage P. falciparum over HepG2 and H460 cells, which is largely driven by the accommodation of the parasite proteasome for a D-amino acid in the P3 position and the preference for a difluorobenzyl group in the P1 position. We isolated the proteasome from P. falciparum cell extracts and determined that the best compound is 171-fold more potent at inhibiting the ß5 subunit of P. falciparum proteasome when compared to the same subunit of the human constitutive proteasome. These compounds also significantly reduce parasitemia in a P. berghei mouse infection model and prolong survival of animals by an average of 6 days. The current epoxyketone inhibitors are ideal starting compounds for orally bioavailable anti-malarial drugs.

LIMO: Latent Inceptionism for Targeted Molecule Generation.

Generation of drug-like molecules with high binding affinity to target proteins remains a difficult and resource-intensive task in drug discovery. Existing approaches primarily employ reinforcement learning, Markov sampling, or deep generative models guided by Gaussian processes, which can be prohibitively slow when generating molecules with high binding affinity calculated by computationally-expensive physics-based methods. We present Latent Inceptionism on Molecules (LIMO), which significantly accelerates molecule generation with an inceptionism-like technique. LIMO employs a variational autoencoder-generated latent space and property prediction by two neural networks in sequence to enable faster gradient-based reverse-optimization of molecular properties. Comprehensive experiments show that LIMO performs competitively on benchmark tasks and markedly outperforms state-of-the-art techniques on the novel task of generating drug-like compounds with high binding affinity, reaching nanomolar range against two protein targets. We corroborate these docking-based results with more accurate molecular dynamics-based calculations of absolute binding free energy and show that one of our generated drug-like compounds has a predicted K D (a measure of binding affinity) of 6 · 10-14 M against the human estrogen receptor, well beyond the affinities of typical early-stage drug candidates and most FDA-approved drugs to their respective targets. Code is available at https://github.com/Rose-STL-Lab/LIMO.

Absolute binding free energy calculations improve enrichment of actives in virtual compound screening.

We determined the effectiveness of absolute binding free energy (ABFE) calculations to refine the selection of active compounds in virtual compound screening, a setting where the more commonly used relative binding free energy approach is not readily applicable. To do this, we conducted baseline docking calculations of structurally diverse compounds in the DUD-E database for three targets, BACE1, CDK2 and thrombin, followed by ABFE calculations for compounds with high docking scores. The docking calculations alone achieved solid enrichment of active compounds over decoys. Encouragingly, the ABFE calculations then improved on this baseline. Analysis of the results emphasizes the importance of establishing high quality ligand poses as starting points for ABFE calculations, a nontrivial goal when processing a library of diverse compounds without informative co-crystal structures. Overall, our results suggest that ABFE calculations can play a valuable role in the drug discovery process.

Mechanistic analysis of light-driven overcrowded alkene-based molecular motors by multiscale molecular simulations.

We analyze light-driven overcrowded alkene-based molecular motors, an intriguing class of small molecules that have the potential to generate MHz-scale rotation rates. The full rotation process is simulated at multiple scales by combining quantum surface-hopping molecular dynamics (MD) simulations for the photoisomerization step with classical MD simulations for the thermal helix inversion step. A Markov state analysis resolves conformational substates, their interconversion kinetics, and their roles in the motor's rotation process. Furthermore, motor performance metrics, including rotation rate and maximal power output, are computed to validate computations against experimental measurements and to inform future designs. Lastly, we find that to correctly model these motors, the force field must be optimized by fitting selected parameters to reference quantum mechanical energy surfaces. Overall, our simulations yield encouraging agreement with experimental observables such as rotation rates, and provide mechanistic insights that may help future designs.

Driving torsion scans with wavefront propagation.

The parameterization of torsional/dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields. Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms. To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values. However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses. In this paper, we propose a systematic and versatile workflow called TorsionDrive to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development. The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described. The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.

Enhanced Diffusion and Chemotaxis of Enzymes.

Many enzymes appear to diffuse faster in the presence of substrate and to drift either up or down a concentration gradient of their substrate. Observations of these phenomena, termed enhanced enzyme diffusion (EED) and enzyme chemotaxis, respectively, lead to a novel view of enzymes as active matter. Enzyme chemotaxis and EED may be important in biology and could have practical applications in biotechnology and nanotechnology. They are also of considerable biophysical interest; indeed, their physical mechanisms are still quite uncertain. This review provides an analytic summary of experimental studies of these phenomena and of the mechanisms that have been proposed to explain them and offers a perspective on future directions for the field.

A Thermodynamic Limit on the Role of Self-Propulsion in Enhanced Enzyme Diffusion.

A number of enzymes reportedly exhibit enhanced diffusion in the presence of their substrates, with a Michaelis-Menten-like concentration dependence. Although no definite explanation of this phenomenon has emerged, a physical picture of enzyme self-propulsion using energy from the catalyzed reaction has been widely considered. Here, we present a kinematic and thermodynamic analysis of enzyme self-propulsion that is independent of any specific propulsion mechanism. Using this theory, along with biophysical data compiled for all enzymes so far shown to undergo enhanced diffusion, we show that the propulsion speed required to generate experimental levels of enhanced diffusion exceeds the speeds of well-known active biomolecules, such as myosin, by several orders of magnitude. Furthermore, the minimal power dissipation required to account for enzyme enhanced diffusion by self-propulsion markedly exceeds the chemical power available from enzyme-catalyzed reactions. Alternative explanations for the observation of enhanced enzyme diffusion therefore merit stronger consideration.

Prediction of Drug-Likeness Using Deep Autoencoder Neural Networks.

Due to diverse reasons, most drug candidates cannot eventually become marketed drugs. Developing reliable computational methods for prediction of drug-likeness of candidate compounds is of vital importance to improve the success rate of drug discovery and development. In this study, we used a fully connected neural networks (FNN) to construct drug-likeness classification models with deep autoencoder to initialize model parameters. We collected datasets of drugs (represented by ZINC World Drug), bioactive molecules (represented by MDDR and WDI), and common molecules (represented by ZINC All Purchasable and ACD). Compounds were encoded with MOLD2 two-dimensional structure descriptors. The classification accuracies of drug-like/non-drug-like model are 91.04% on WDI/ACD databases, and 91.20% on MDDR/ZINC, respectively. The performance of the models outperforms previously reported models. In addition, we develop a drug/non-drug-like model (ZINC World Drug vs. ZINC All Purchasable), which distinguishes drugs and common compounds, with a classification accuracy of 96.99%. Our work shows that by using high-latitude molecular descriptors, we can apply deep learning technology to establish state-of-the-art drug-likeness prediction models.