Molecular dynamics simulations as a guide for modulating small molecule aggregation.

Small colloidally aggregating molecules (SCAMs) can be problematic for biological assays in drug discovery campaigns. However, the self-associating properties of SCAMs have potential applications in drug delivery and analytical biochemistry. Consequently, the ability to predict the aggregation propensity of a small organic molecule is of considerable interest. Chemoinformatics-based filters such as ChemAGG and Aggregator Advisor offer rapid assessment but are limited by the assay quality and structural diversity of their training set data. Complementary to these tools, we explore here the ability of molecular dynamics (MD) simulations as a physics-based method capable of predicting the aggregation propensity of diverse chemical structures. For a set of 32 molecules, using simulations of 100 ns in explicit solvent, we find a success rate of 97% (one molecule misclassified) as opposed to 75% by Aggregator Advisor and 72% by ChemAGG. These short timescale MD simulations are representative of longer microsecond trajectories and yield an informative spectrum of aggregation propensities across the set of solutes, capturing the dynamic behaviour of weakly aggregating compounds. Implicit solvent simulations using the generalized Born model were less successful in predicting aggregation propensity. MD simulations were also performed to explore structure-aggregation relationships for selected molecules, identifying chemical modifications that reversed the predicted behaviour of a given aggregator/non-aggregator compound. While lower throughput than rapid cheminformatics-based SCAM filters, MD-based prediction of aggregation has potential to be deployed on the scale of focused subsets of moderate size, and, depending on the target application, provide guidance on removing or optimizing a compound's aggregation propensity.

Brazilin is a natural product inhibitor of the NLRP3 inflammasome.

Excessive or aberrant NLRP3 inflammasome activation has been implicated in the progression and initiation of many inflammatory conditions; however, currently no NLRP3 inflammasome inhibitors have been approved for therapeutic use in the clinic. Here we have identified that the natural product brazilin effectively inhibits both priming and activation of the NLRP3 inflammasome in cultured murine macrophages, a human iPSC microglial cell line and in a mouse model of acute peritoneal inflammation. Through computational modeling, we predict that brazilin can adopt a favorable binding pose within a site of the NLRP3 protein which is essential for its conformational activation. Our results not only encourage further evaluation of brazilin as a therapeutic agent for NLRP3-related inflammatory diseases, but also introduce this small-molecule as a promising scaffold structure for the development of derivative NLRP3 inhibitor compounds.

A neural network potential based on pairwise resolved atomic forces and energies.

Molecular simulations have become a key tool in molecular and materials design. Machine learning (ML)-based potential energy functions offer the prospect of simulating complex molecular systems efficiently at quantum chemical accuracy. In previous work, we have introduced the ML-based PairF-Net approach to neural network potentials, that adopts a pairwise interatomic scheme to predicting forces within a molecular system. Here, we further develop the PairF-Net model to intrinsically incorporate energy conservation and couple the model to a molecular mechanical (MM) environment within the OpenMM package. The updated PairF-Net model yields energy and force predictions and dynamical distributions in good agreement with the rMD17 dataset of ten small organic molecules in the gas-phase. We further show that these in vacuo ML models of small molecules can be applied to force predictions in aqueous solution via hybrid ML/MM simulations. We present a new benchmark dataset for these ten molecules in solution, obtained from QM/MM simulations, which we denote as rMD17-aq (https://zenodo.org/records/10048644); and assess the ability of PairF-Net to reproduce the molecular energy, atomic forces and dynamical distributions of these solution conformations via ML/MM simulations.

Cellulose Iß microfibril interaction with pristine graphene in water: Effects of amphiphilicity by molecular simulation.

Graphene-cellulose interactions have considerable potential in the development of new materials. In previous computational work (Biomacromolecules2016, 16, 1771), we predicted that the model 100 hydrophobic surface of cellulose interacted favourably with pristine graphene in aqueous solution molecular dynamics simulations; conversely, a model of the hydrophilic 010 surface of cellulose exhibited progressive rearrangement to present a more hydrophobic face with the graphene, with weakened hydrogen bonds between cellulose chains and partial permeation of water. Here, we extend this work by simulating the interaction in aqueous solution of the amphiphilic 110 surface of a cellulose Iß microfibril model, comprising 36 chains of 40 glucosyl residues, with an infinite sheet of pristine graphene. This face of the microfibril is of intermediate hydrophilicity and progressively associates with graphene over replicate simulations. As cellulose chains adhere to the graphene surface, forming interactions via its CH and OH groups, we observe a degree of local and global untwisting of the microfibril. Complementary rippling of the graphene surface is also observed, as it adapts to interaction with the microfibril. This adsorption process is accompanied by increased exclusion of water between cellulose and graphene although some water localises between chains at the immediate interface. The predicted propensity of a cellulose microfibril to adsorb spontaneously on the graphene surface, with mutual structural accommodation, highlights the amphiphilic nature of cellulose and the types of interactions that can be harnessed to design new graphene-carbohydrate biopolymer materials.

A Selective Review and Virtual Screening Analysis of Natural Product Inhibitors of the NLRP3 Inflammasome.

The NLRP3 inflammasome is currently an exciting target for drug discovery due to its role in various inflammatory diseases; however, to date, no NLRP3 inhibitors have reached the clinic. Several studies have used natural products as hit compounds to facilitate the design of novel selective NLRP3 inhibitors. Here, we review selected natural products reported in the literature as NLRP3 inhibitors, with a particular focus on those targeting gout. To complement this survey, we also report a virtual screen of the ZINC20 natural product database, predicting favored chemical features that can aid in the design of novel small molecule NLRP3 inhibitors.

Modeling pyranose ring pucker in carbohydrates using machine learning and semi-empirical quantum chemical methods.

Pyranose ring pucker is a key coordinate governing the structure, interactions and reactivity of carbohydrates. We assess the ability of the machine learning potentials, ANI-1ccx and ANI-2x, and the GFN2-xTB semiempirical quantum chemical method, to model ring pucker conformers of five monosaccharides and oxane in the gas phase. Relative to coupled-cluster quantum mechanical calculations, we find that ANI-1ccx most accurately reproduces the ring pucker energy landscape for these molecules, with a correlation coefficient r2 of 0.83. This correlation in relative energies lowers to values of 0.70 for ANI-2x and 0.60 for GFN2-xTB. The ANI-1ccx also provides the most accurate estimate of the energetics of the 4 C1 -to-1 C4 minimum energy pathway for the six molecules. All three models reproduce chair more accurately than non-chair geometries. Analysis of small model molecules suggests that the ANI-1ccx model favors puckers with equatorial hydrogen bonding substituents; that ANI-2x and GFN2-xTB models overstabilize conformers with axially oriented groups; and that the endo-anomeric effect is overestimated by the machine learning models and underestimated via the GFN2-xTB method. While the pucker conformers considered in this study correspond to a gas phase environment, the accuracy and computational efficiency of the ANI-1ccx approach in modeling ring pucker in vacuo provides a promising basis for future evaluation and application to condensed phase environments.

Probing the effect of NEK7 and cofactor interactions on dynamics of NLRP3 monomer using molecular simulation.

The NLRP3 inflammasome is a cytoplasmic complex that regulates the activation of inflammatory cytokines and, given its implication in a range of diseases, is an important therapeutic target. The cofactor ATP and the centrosomal kinase NEK7 are important for NLRP3 activation. Here we have constructed and simulated computational models of full-length monomeric NLRP3 to shed light on the importance of NEK7 and cofactor interactions for its conformation and dynamics in aqueous solution. We find that molecular dynamics simulation reproduces well the features of the recently published cryo-EM structure of the ADP-bound NLRP3-NEK7 complex; on the removal of NEK7, the NLRP3 molecule adopts a more compact closed form during simulations. Replacement of ADP by ATP promotes a rearrangement of hydrogen-bonding interactions, domain interfaces, and a degree of opening of the NLRP3 conformation. We also examine the dynamics of an acidic loop of the LRR domain of NLRP3, which samples in a region observed in the NEK7-bound cryo-EM structure but not in an oligomeric form of inactive NLRP3. During the molecular dynamics simulations of NLRP3, we find some plasticity in its topology that suggests access routes for ATP to the cofactor pocket not immediately evident from the existing NEK7-bound cryo-EM structure. These computed dynamical trajectories of NLRP3 provide insight into coordinates of deformation that may be key for cofactor binding and inflammasome activation.

Structure-based assessment and druggability classification of protein-protein interaction sites.

The featureless interface formed by protein-protein interactions (PPIs) is notorious for being considered a difficult and poorly druggable target. However, recent advances have shown PPIs to be druggable, with the discovery of potent inhibitors and stabilizers, some of which are currently being clinically tested and approved for medical use. In this study, we assess the druggability of 12 commonly targeted PPIs using the computational tool, SiteMap. After evaluating 320 crystal structures, we find that the PPI binding sites have a wide range of druggability scores. This can be attributed to the unique structural and physiochemical features that influence their ligand binding and concomitantly, their druggability predictions. We then use these features to propose a specific classification system suitable for assessing PPI targets based on their druggability scores and measured binding-affinity. Interestingly, this system was able to distinguish between different PPIs and correctly categorize them into four classes (i.e. very druggable, druggable, moderately druggable, and difficult). We also studied the effects of protein flexibility on the computed druggability scores and found that protein conformational changes accompanying ligand binding in ligand-bound structures result in higher protein druggability scores due to more favorable structural features. Finally, the drug-likeness of many published PPI inhibitors was studied where it was found that the vast majority of the 221 ligands considered here, including orally tested/marketed drugs, violate the currently acceptable limits of compound size and hydrophobicity parameters. This outcome, combined with the lack of correlation observed between druggability and drug-likeness, reinforces the need to redefine drug-likeness for PPI drugs. This work proposes a PPI-specific classification scheme that will assist researchers in assessing the druggability and identifying inhibitors of the PPI interface.

SAR of Novel 3-Arylisoquinolinones: meta-Substitution on the Aryl Ring Dramatically Enhances Antiproliferative Activity through Binding to Microtubules.

A set of meta-substituted 3-arylisoquinolinones have been identified that show substantial cytotoxicity in breast, liver, lung and colon cancer cell lines; these are up to 700-fold more active than the corresponding para analogues. These compounds were initially proposed as inhibitors of N-ribosyl dihydronicotinamide (NRH): quinone oxidoreductase 2 (NQO2) but were found to be inactive against the enzyme. Instead, COMPARE analysis suggested that 6-fluoro-3-(meta-fluorophenyl)isoquinolin-1(2H)-one (4) could mimic colchicine and interact with microtubules, a recognized target for cancer therapy. Subsequent docking, molecular dynamics simulations, and free energy analysis further suggested that compound 4 bound well into the colchicine-binding pocket of tubulin. Indeed, 4 suppressed tubulin polymerization, caused G2/M cell cycle arrest, and induced apoptosis. Also, 4 inhibited the formation of endothelial cell capillary-like tubes and further disrupted the structure of preestablished tubes; the effects were not observed with para analogue 5. In accordance with this, the computed free energy of binding of 5 to tubulin was lower in magnitude than that for 4 and appeared to arise in part from the inability of the para substituent to occupy a tubulin subpocket, which is possible in the meta orientation. In conclusion, the antiproliferative potential of the novel 3-arylisoquinolinones is markedly influenced by a subtle change in the structure (meta versus para). The meta-substituted isoquinolinone 4 is a microtubule-destabilizing agent with potential tumor-selectivity and antiangiogenic and vascular disrupting features.

Direct Deoxydehydration of Cyclic trans-Diol Substrates: An Experimental and Computational Study of the Reaction Mechanism of Vanadium(V)-based Catalysis*.

The deoxydehydration of carbohydrates represents a key target to leverage renewable biomass resources chemically. Using a vanadium(V)-based catalyst, it was possible to directly deoxydehydrate cyclic trans-diol substrates. Accompanying mechanistic characterisation of this process by density functional calculations pointed to an energetically tractable route for deoxydehydration of cyclic trans-diol substrates involving stepwise cleavage of the diol C-O bonds via the triplet state; experimentally, this was supported by light dependence of the reaction. Calculations also indicated that cyclic cis-diols and a linear diol substrate could additionally proceed by a concerted singlet DODH mechanism. This work potentially opens a new and cost-effective way to efficiently convert carbohydrates of trans-diol stereochemistry into alkenes.

Generation of Phenothiazine with Potent Anti-TLK1 Activity for Prostate Cancer Therapy.

Through in vitro kinase assays and docking studies, we report the synthesis and biological evaluation of a phenothiazine analog J54 with potent TLK1 inhibitory activity for prostate cancer (PCa) therapy. Most PCa deaths result from progressive failure in standard androgen deprivation therapy (ADT), leading to metastatic castration-resistant PCa. Treatments that can suppress the conversion to mCRPC have high potential to be rapidly implemented in the clinics. ADT results in increased expression of TLK1B, a key kinase upstream of NEK1 and ATR and mediating the DNA damage response that typically results in temporary cell-cycle arrest of androgen-responsive PCa cells, whereas its abrogation leads to apoptosis. We studied J54 as a potent inhibitor of this axis and as a mediator of apoptosis in vitro and in LNCaP xenografts, which has potential for clinical investigation in combination with ADT. J54 has low affinity for the dopamine receptor in modeling and competition studies and weak detrimental behavioral effects in mice and C. elegans.

How Do Small Molecule Aggregates Inhibit Enzyme Activity? A Molecular Dynamics Study.

Small molecule compounds which form colloidal aggregates in solution are problematic in early drug discovery; adsorption of the target protein by these aggregates can lead to false positives in inhibition assays. In this work, we probe the molecular basis of this inhibitory mechanism using molecular dynamics simulations. Specifically, we examine in aqueous solution the adsorption of the enzymes ß-lactamase and PTP1B onto aggregates of the drug miconazole. In accordance with experiment, molecular dynamics simulations observe formation of miconazole aggregates as well as subsequent association of these aggregates with ß-lactamase and PTP1B. When complexed with aggregate, the proteins do not exhibit significant alteration in protein tertiary structure or dynamics on the microsecond time scale of the simulations, but they do indicate persistent occlusion of the protein active site by miconazole molecules. MD simulations further suggest this occlusion can occur via surficial interactions of protein with miconazole but also potentially by envelopment of the protein by miconazole. The heterogeneous polarity of the miconazole aggregate surface seems to underpin its activity as an invasive and nonspecific inhibitory agent. A deeper understanding of these protein/aggregate systems has implications not only for drug design but also for their exploitation as tools in drug delivery and analytical biochemistry.

What Next for Quantum Mechanics in Structure-Based Drug Discovery?

There is significant potential for electronic structure methods to improve the quality of the predictions furnished by the tools of computer-aided drug design, which typically rely on empirically derived functions. In this perspective, we consider some recent examples of how quantum mechanics has been applied in predicting protein-ligand geometries, protein-ligand binding affinities and ligand strain on binding. We then outline several significant developments in quantum mechanics methodology likely to influence these approaches: in particular, we note the advent of more computationally expedient ab initio quantum mechanical methods that can provide chemical accuracy for larger molecular systems than hitherto possible. We highlight the emergence of increasingly accurate semiempirical quantum mechanical methods and the associated role of machine learning and molecular databases in their development. Indeed, the convergence of improved algorithms for solving and analyzing electronic structure, modern machine learning methods, and increasingly comprehensive benchmark data sets of molecular geometries and energies provides a context in which the potential of quantum mechanics will be increasingly realized in driving future developments and applications in structure-based drug discovery.

Lead Optimization of Dehydroemetine for Repositioned Use in Malaria.

Drug repositioning offers an effective alternative to de novo drug design to tackle the urgent need for novel antimalarial treatments. The antiamoebic compound emetine dihydrochloride has been identified as a potent in vitro inhibitor of the multidrug-resistant strain K1 of Plasmodium falciparum (50% inhibitory concentration [IC50], 47 nM ± 2.1 nM [mean ± standard deviation]). Dehydroemetine, a synthetic analogue of emetine dihydrochloride, has been reported to have less-cardiotoxic effects than emetine. The structures of two diastereomers of dehydroemetine were modeled on the published emetine binding site on the cryo-electron microscopy (cryo-EM) structure with PDB code 3J7A (P. falciparum 80S ribosome in complex with emetine), and it was found that (-)-R,S-dehydroemetine mimicked the bound pose of emetine more closely than did (-)-S,S-dehydroisoemetine. (-)-R,S-dehydroemetine (IC50 71.03 ± 6.1 nM) was also found to be highly potent against the multidrug-resistant K1 strain of P. falciparum compared with (-)-S,S-dehydroisoemetine (IC50, 2.07 ± 0.26 µM), which loses its potency due to the change of configuration at C-1'. In addition to its effect on the asexual erythrocytic stages of P. falciparum, the compound exhibited gametocidal properties with no cross-resistance against any of the multidrug-resistant strains tested. Drug interaction studies showed (-)-R,S-dehydroemetine to have synergistic antimalarial activity with atovaquone and proguanil. Emetine dihydrochloride and (-)-R,S-dehydroemetine failed to show any inhibition of the hERG potassium channel and displayed activity affecting the mitochondrial membrane potential, indicating a possible multimodal mechanism of action.

Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations.

The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 µs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on ß-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of ß-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the 1C4 ring conformation; this observation may have implications for understanding the structural basis of the biological function of ß-GalNAc-containing glycosaminoglycans such as dermatan sulfate.

Discovery of potent 4-aminoquinoline hydrazone inhibitors of NRH:quinoneoxidoreductase-2 (NQO2).

(NRH):quinone oxidoreductase 2 (NQO2) is associated with various processes involved in cancer initiation and progression probably via the production of ROS during quinone metabolism. Thus, there is a need to develop inhibitors of NQO2 that are active in vitro and in vivo. As part of a strategy to achieve this we have used the 4-aminoquinoline backbone as a starting point and synthesized 21 novel analogues. The syntheses utilised p-anisidine with Meldrum's acid and trimethyl orthoacetate or trimethyl orthobenzoate to give the 4-hydrazin-quinoline scaffold, which was derivatised with aldehydes or acid chlorides to give hydrazone or hydrazide analogues, respectively. The hydrazones were the most potent inhibitors of NQO2 in cell free systems, some with low nano-molar IC50 values. Structure-activity analysis highlighted the importance of a small substituent at the 2-position of the 4-aminoquinoline ring, to reduce steric hindrance and improve engagement of the scaffold within the NQO2 active site. Cytotoxicity and NQO2-inhibitory activity in vitro was evaluated using ovarian cancer SKOV-3 and TOV-112 cells (expressing high and low levels of NQO2, respectively). Generally, the hydrazones were more toxic than hydrazide analogues and further, toxicity is unrelated to cellular NQO2 activity. Pharmacological inhibition of NQO2 in cells was measured using the toxicity of CB1954 as a surrogate end-point. Both the hydrazone and hydrazide derivatives are functionally active as inhibitors of NQO2 in the cells, but at different inhibitory potency levels. In particular, 4-((2-(6-methoxy-2-methylquinolin-4-yl)hydrazono)methyl)phenol has the greatest potency of any compound yet evaluated (53 nM), which is 50-fold lower than its toxicity IC50. This compound and some of its analogues could serve as useful pharmacological probes to determine the functional role of NQO2 in cancer development and response to therapy.

Tracing the GSAP-APP C-99 Interaction Site in the ß-Amyloid Pathway Leading to Alzheimer's Disease.

Gamma secretase activating protein (GSAP) present in ß-amyloid pathway orchestrates the formation of ß-amyloid plaques by γ-secretase activation and is an emerging therapeutic target for the treatment of Alzheimer's disease. It forms a ternary complex with γ-secretase and APP C-99. However, there are limited reports for the interaction of APP C-99 with GSAP. Here, we report the characterization of purified maltose binding protein (MBP) tagged human GSAP and its interaction with synthetic APP C-99 peptide fragments (712IATVIVITLVMLKKQ727 (712IQ727), 719TLVMLKKKQYTSIHHGVVEVDAAVT743 (719TT743) 734GVVEVDAAVTPEERHLSKMQQNGY757 (734GY757), and 746ERHLSKMQQNGYENPTYKFFEQMQN770 (746EN770)). The results emphasize the selective interaction of peptide (719TT743) with MBP-GSAP with a dissociation constant of 0.136 µM. Further, computational modeling of the GSAP-719TT743 complex finds an optimal bound pose of 719TT743 within an extended groove on the surface of GSAP. The preliminary results highlight the interaction between the two major proteins in the plausible ternary complex: APP C-99-GSAP-γ-secretase. It paves a futuristic path to investigate the GSAP-APP C-99 binding in detail and accentuates the role of GSAP in the ß-amyloid pathway.

Negative Cooperativity in NAD(P)H Quinone Oxidoreductase 1 (NQO1).

NAD(P)H quinone oxidoreductase-1 (NQO1) is a homodimeric protein that acts as a detoxifying enzyme or as a chaperone protein. Dicourmarol interacts with NQO1 at the NAD(P)H binding site and can both inhibit enzyme activity and modulate the interaction of NQO1 with other proteins. We show that the binding of dicoumarol and related compounds to NQO1 generates negative cooperativity between the monomers. This does not occur in the presence of the reducing cofactor, NAD(P)H, alone. Alteration of Gly150 (but not Gly149 or Gly174) abolished the dicoumarol-induced negative cooperativity. Analysis of the dynamics of NQO1 with the Gaussian network model indicates a high degree of collective motion by monomers and domains within NQO1. Ligand binding is predicted to alter NQO1 dynamics both proximal to the ligand binding site and remotely, close to the second binding site. Thus, drug-induced modulation of protein motion might contribute to the biological effects of putative inhibitors of NQO1.

Evaluation of analogues of furan-amidines as inhibitors of NQO2.

Inhibitors of the enzyme NQO2 (NRH: quinone oxidoreductase 2) are of potential use in cancer chemotherapy and malaria. We have previously reported that non-symmetrical furan amidines are potent inhibitors of NQO2 and here novel analogues are evaluated. The furan ring has been changed to other heterocycles (imidazole, N-methylimidazole, oxazole, thiophene) and the amidine group has been replaced with imidate, reversed amidine, N-arylamide and amidoxime to probe NQO2 activity, improve solubility and decrease basicity of the lead furan amidine. All compounds were fully characterised spectroscopically and the structure of the unexpected product N-hydroxy-4-(5-methyl-4-phenylfuran-2-yl)benzamidine was established by X-ray crystallography. The analogues were evaluated for inhibition of NQO2, which showed lower activity than the lead furan amidine. The observed structure-activity relationship for the furan-amidine series with NQO2 was rationalized by preliminary molecular docking and binding mode analysis. In addition, the oxazole-amidine analogue inhibited the growth of Plasmodium falciparum with an IC50 value of 0.3 µM.