Rationalising the retro-Diels-Alder fragmentation pattern of viscutins using electrospray interface-tandem mass spectrometry coupled to theoretical modelling.

RATIONALE: Although mass spectrometry (MS) is a powerful tool in structural elucidation of unknown flavonoids based on their unique fragmentation patterns, proposing the correct fragmentation mechanism is still a challenge from tandem mass spectrometry data only. In recent years, computational tools such as molecular networking and MS2LDA have played a major role in the identification of structurally related compounds through an in-depth survey of their fragmentation patterns. METHODS: Therefore, in this study, three viscutin molecules in Viscum combreticola Engl. crude extracts were characterised using ultra-high-performance liquid chromatography-quadrupole-time-of-flight tandem mass spectrometry and MS2LDA, a computational tool. Ion-trap mass spectrometry and density functional theoretical modelling were used as confirmatory tools to rationalise the unique fragmentation patterns observed for these molecules. RESULTS: Here, MS2LDA revealed the presence of a unique Mass2Motif in all the three viscutin molecules at m/z 137, which was confirmed to be a 1,3 A- RDA (retro-Diels-Alder) fragmentation product using liquid chromatography-ion-trap mass spectrometry and density functional theoretical modelling. Moreover, MS2LDA proved to be useful in differentiating this spectral feature that was specific to viscutin molecules in the presence of other isobaric ions at m/z 137 occurring in compounds in other molecular families. CONCLUSION: Therefore, the results of the current study showed that computational tools such as MS2LDA are essential in uncovering some gas-phase fragmentation reactions of molecules in MS and that theoretical modelling is a powerful tool in rationalising these reactions in metabolite identification.

New Insights into the (A)Synchronicity of Diels-Alder Reactions: A Theoretical Study Based on the Reaction Force Analysis and Atomic Resolution of Energy Derivatives.

In the present manuscript, we report new insights into the concept of (a)synchronicity in Diels-Alder (DA) reactions in the framework of the reaction force analysis in conjunction with natural population calculations and the atomic resolution of energy derivatives along the intrinsic reaction coordinate (IRC) path. Our findings suggest that the DA reaction transitions from a preferentially concerted mechanism to a stepwise one in a 0.10 Å window of synchronicity indices ranging from 0.90 to 1.00 Å. We have also shown that the relative position of the global minimum of the reaction force constant with respect to the TS is an alternative and quantifiable indicator of the (a)synchronicity in DA reactions. Moreover, the atomic resolution of energy derivatives reveals that the mechanism of the DA reaction involves two inner elementary processes associated with the formation of each of the two C-C bonds. This resolution goes on to indicate that, in asynchronous reactions, the driving and retarding components of the reaction force are mostly due to the fast and slow-forming C-C bonds (elementary processes) respectively, while in synchronous reactions, both elementary processes retard and drive the process concomitantly and equivalently.

Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases.

It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1-PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro-CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.

Force Field Parameters for Fe2+4S2-4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies.

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe2+4S2-4 clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe2+4S2-4 clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe2+ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe2+ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe2+4S2-4 cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

Synthesis and anti-parasitic activity of N-benzylated phosphoramidate Mg2+-chelating ligands.

A series of N-benzylated phosphoramidate esters, containing a 3,4-dihydroxyphenyl Mg2+-chelating group, has been synthesised in five steps as analogues of fosmidomycin, a Plasmodium falciparum 1-deoxy-1-d-xylulose-5-phosphate reductoisomerase (PfDXR) inhibitor. The 3,4-dihydroxyphenyl group effectively replaces the Mg2+-chelating hydroxamic acid group in fosmidomycin. The compounds showed very encouraging anti-parasitic activity with IC50 values of 5.6-16.4 µM against Plasmodium falciparum parasites and IC50 values of 5.2 - 10.2 µM against Trypanosoma brucei brucei (T.b.brucei). Data obtained from in silico docking of the ligands in the PfDXR receptor cavity (3AU9)5 support their potential as PfDXR inhibitors.

Synthesis and anti-parasitic activity of achiral N-benzylated phosphoramidic acid derivatives.

Synthetic pathways have been developed to access a series of N-benzylated phosphoramidic acid derivatives as novel, achiral analogues of the established Plasmodium falciparum 1-deoxy-d-xylulose-5-phosphate reductase (PfDXR) enzyme inhibitor, FR900098. Bioassays of the targeted compounds and their synthetic precursors have revealed minimal antimalarial activity but encouraging anti-trypanosomal activity - in one case with an IC50 value of 5.4 µM against Trypanosoma brucei, the parasite responsible for Nagana (African cattle sleeping sickness). The results of relevant in silico modelling and docking studies undertaken in the design and evaluation of these compounds are discussed.

Establishing Computational Approaches Towards Identifying Malarial Allosteric Modulators: A Case Study of Plasmodium falciparum Hsp70s.

Combating malaria is almost a never-ending battle, as Plasmodium parasites develop resistance to the drugs used against them, as observed recently in artemisinin-based combination therapies. The main concern now is if the resistant parasite strains spread from Southeast Asia to Africa, the continent hosting most malaria cases. To prevent catastrophic results, we need to find non-conventional approaches. Allosteric drug targeting sites and modulators might be a new hope for malarial treatments. Heat shock proteins (HSPs) are potential malarial drug targets and have complex allosteric control mechanisms. Yet, studies on designing allosteric modulators against them are limited. Here, we identified allosteric modulators (SANC190 and SANC651) against P. falciparum Hsp70-1 and Hsp70-x, affecting the conformational dynamics of the proteins, delicately balanced by the endogenous ligands. Previously, we established a pipeline to identify allosteric sites and modulators. This study also further investigated alternative approaches to speed up the process by comparing all atom molecular dynamics simulations and dynamic residue network analysis with the coarse-grained (CG) versions of the calculations. Betweenness centrality (BC) profiles for PfHsp70-1 and PfHsp70-x derived from CG simulations not only revealed similar trends but also pointed to the same functional regions and specific residues corresponding to BC profile peaks.

Synthetic analogues of the marine bisindole deoxytopsentin: potent selective inhibitors of MRSA pyruvate kinase.

As part of an ongoing study to elucidate the SAR of bisindole alkaloid inhibitors against the evolutionary conserved MRSA pyruvate kinase (PK), we present here the synthesis and biological activity of six dihalogenated analogues of the naturally occurring sponge metabolite deoxytopsentin, including the naturally occurring dibromodeoxytopsentin. The most active compounds displayed potent low nanomolar inhibitory activity against MRSA PK with concomitant significant selectivity for MRSA PK over human PK orthologues. Computational studies suggest that these potent MRSA PK inhibitors occupy a region of the small interface of the enzyme tetramer where amino acid sequence divergence from common human PK orthologues may contribute to the observed selectivity.

A Computational Study of Green Tea Extracts and their Derivatives as Potential Inhibitors for Squalene Monooxygenase.

BACKGROUND: According to the World Health Organisation, cardiovascular complications have been recognized as the leading course of death between 2000 and 2019. Cardiovascular complications are caused by excess LDL cholesterol in the body or arteries that can build up to form a plaque. There are drugs currently in clinical use called statins that target HMGCoA reductase. However, these drugs result in several side effects. This work investigated using computational approaches to lower cholesterol by investigating green tea extracts as an inhibitors for squalene monooxygenase (the second-rate-controlling step in cholesterol synthesis). METHODS: Pharmacophore modeling was done to identify possible pharmacophoric sites based on the pIC50 values. The best hypothesis generated by pharmacophore modeling was further validated by atom-based 3D QSAR, where 70% of the data set was treated as the training set. Prior molecular docking ADMET studies were done to investigate the physiochemical properties of these molecules. Glide docking was performed, followed by molecular dynamics to evaluate the protein conformational changes. RESULTS: Pharmacophore results suggest that the best molecules to interact with the biological target should have at least one hydrogen acceptor (A5), two hydrogen donors (D9 and D10), and two benzene rings (R14 and R15) for green tea polyphenols and theasinensin A. ADMET result shows that all molecules in this class have low oral adsorption. Molecular docking results showed that some green tea polyphenols have good binding affinities, with most of these structures having a docking score of less than -10 kcal/mol. Molecular dynamics further illustrated that the best-docked ligands perfectly stay within the active site over a 100 ns simulation. CONCLUSION: The results obtained from this study suggest that green tea polyphenols have the potential for inhibition of squalene monooxygenase, except for theasinensin A.

Inhibitory effects of selected cannabinoids against dipeptidyl peptidase IV, an enzyme linked to type 2 diabetes.

Ethnopharmacological relevance: In recent times the decriminalisation of cannabis globally has increased its use as an alternative medication. Where it has been used in modern medicinal practises since the 1800s, there is limited scientific investigation to understand the biological activities of this plant. Aim of the study: Dipeptidyl peptidase IV (DPP-IV) plays a key role in regulating glucose homeostasis, and inhibition of this enzyme has been used as a therapeutic approach to treat type 2 diabetes. However, some of the synthetic inhibitors for this enzyme available on the market may cause undesirable side effects. Therefore, it is important to identify new inhibitors of DPP-IV and to understand their interaction with this enzyme. Methods: In this study, four cannabinoids (cannabidiol, cannabigerol, cannabinol and Δ9-tetrahydrocannabinol) were evaluated for their inhibitory effects against recombinant human DPP-IV and their potential inhibition mechanism was explored using both in vitro and in silico approaches. Results: All four cannabinoids resulted in a dose-dependent response with IC50 values of between 4.0 and 6.9 µg/mL. Kinetic analysis revealed a mixed mode of inhibition. CD spectra indicated that binding of cannabinoids results in structural and conformational changes in the secondary structure of the enzyme. These findings were supported by molecular docking studies which revealed best docking scores at both active and allosteric sites for all tested inhibitors. Furthermore, molecular dynamics simulations showed that cannabinoids formed a stable complex with DPP-IV protein via hydrogen bonds at an allosteric site, suggesting that cannabinoids act by either inducing conformational changes or blocking the active site of the enzyme. Conclusion: These results demonstrated that cannabinoids may modulate DPP-IV activity and thereby potentially assist in improving glycaemic regulation in type 2 diabetes.

Exploring DOXP-reductoisomerase binding limits using phosphonated N-aryl and N-heteroarylcarboxamides as DXR inhibitors.

DOXP-reductoisomerase (DXR) is a validated target for the development of antimalarial drugs to address the increase in resistant strains of Plasmodium falciparum. Series of aryl- and heteroarylcarbamoylphosphonic acids, their diethyl esters and disodium salts have been prepared as analogues of the potent DXR inhibitor fosmidomycin. The effects of the carboxamide N-substituents and the length of the methylene linker have been explored using in silico docking studies, saturation transfer difference NMR spectroscopy and enzyme inhibition assays using both EcDXR and PfDXR. These studies indicate an optimal linker length of two methylene units and have confirmed the importance of an additional binding pocket in the PfDXR active site. Insights into the constraints of the PfDXR binding site provide additional scope for the rational design of DXR inhibitors with increased ligand-receptor interactions.

The formation of 2,2,4-trimethyl-2,3-dihydro-1H-1,5-benzodiazepine from 1,2-diaminobenzene in the presence of acetone.

In an attempt to synthesize a 2-substituted benzimidazole from the reaction of o-phenylenediamine and isophthalic acid in the presence of acetone and ethanol under microwave irradiation, a salt of the isophthalate ion and 2,2,4-trimethyl-2,3-dihydro-1H-1,5-benzodiazepin-5-ium ion was obtained. The condensation of two moles of acetone with the amine groups resulted in the formation of the benzodiazepine which crystallized as an iminium cation forming a salt with the isophthalate anion. The formation of benzodiazepine was also confirmed by performing the reaction of o-phenylenediamine with excess acetone in ethanol under conventional heating conditions. The compounds were characterized by NMR, FTIR, HRMS and microanalysis as well as X-ray crystallography. The reaction mechanism leading to the formation of benzodiazepine is also discussed.

In Silico Substrate-Binding Profiling for SARS-CoV-2 Main Protease (Mpro) Using Hexapeptide Substrates.

The SARS-CoV-2 main protease (Mpro) is essential for the life cycle of the COVID-19 virus. It cleaves the two polyproteins at 11 positions to generate mature proteins for virion formation. The cleavage site on these polyproteins is known to be Leu-Gln↓(Ser/Ala/Gly). A range of hexapeptides that follow the known sequence for recognition and cleavage was constructed using RDKit libraries and complexed with the crystal structure of Mpro (PDB ID 6XHM) through extensive molecular docking calculations. A subset of 131 of these complexes underwent 20 ns molecular dynamics simulations. The analyses of the trajectories from molecular dynamics included principal component analysis (PCA), and a method to compare PCA plots from separate trajectories was developed in terms of encoding PCA progression during the simulations. The hexapeptides formed stable complexes as expected, with reproducible molecular docking of the substrates given the extensiveness of the procedure. Only Lys-Leu-Gln*** (KLQ***) sequence complexes were studied for molecular dynamics. In this subset of complexes, the PCA analysis identified four classifications of protein motions across these sequences. KLQ*** complexes illustrated the effect of changes in substrate on the active site, with implications for understanding the substrate recognition of Mpro and informing the development of small molecule inhibitors.

Drug Resistance in the HIV-1 Subtype C Protease Enzyme: A High Throughput Virtual Screening Approach in Search of New Ligands with Activity.

BACKGROUND: HIV-1 subtype C protease is a strategic target for antiretroviral treatment. However, resistance to protease inhibitors appears after months of treatment. Chromones and 2- biscoumarin derivatives show potential for inhibition of the HIV- subtype C protease. OBJECTIVE: Different heterocyclic structures from the ZINC database were docked against Human Immunodeficiency Virus-1 (HIV) subtype C protease crystal structure 2R5Q and 2R5P. The 5 best molecules were selected to be docked against 62 homology models based on HIV-protease sequences from infants failing antiretroviral protease treatment. This experimentation was performed with two molecular docking programs: Autodock and Autodock Vina. These molecules were modified by substituting protons with different moieties, and the derivatives were docked against the same targets. Ligand-protein interactions, physical/chemical proprieties of the molecules, and dynamics simulations were analyzed. METHODS: Docking of all of the molecules was performed to find out the binding sites of HIV-1 subtype C proteases. An in-house script was made to substitute protons of molecules with different moieties. According to the Lipinski rule of five, physical and chemical properties were determined. Complexes of certain ligands-protease were compared to the protein alone in molecular dynamics simulations. RESULTS: From the first docking results, the 5 best (lowest energy) ligands (dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[a]acridine) were selected. The binding energy of the modified ligands increased, including the poorest-performing molecules. A correlation between nature, the position, and the resulting binding energy was observed. According to the Lipinski rules, the physico-chemical characteristics of the five best-modified ligands are ideal for oral bioavailability. Molecular dynamics simulations show that some lead-protease complexes were stable. CONCLUSION: Dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[ a]acridine and their derivatives might be considered as promising HIV-1 subtype C protease inhibitors. This could be confirmed through synthesis and subsequent in vitro assays.

The Molecular Basis of the Effect of Temperature on the Structure and Function of SARS-CoV-2 Spike Protein.

The remarkable rise of the current COVID-19 pandemic to every part of the globe has raised key concerns for the current public healthcare system. The spike (S) protein of SARS-CoV-2 shows an important part in the cell membrane fusion and receptor recognition. It is a key target for vaccine production. Several researchers studied the nature of this protein under various environmental conditions. In this work, we applied molecular modeling and extensive molecular dynamics simulation approaches at 0°C (273.15 K), 20°C (293.15 K), 40°C (313.15 K), and 60°C (333.15 K) to study the detailed conformational alterations in the SARS-CoV-2 S protein. Our aim is to understand the influence of temperatures on the structure, function, and dynamics of the S protein of SARS-CoV-2. The structural deviations, and atomic and residual fluctuations were least at low (0°C) and high (60°C) temperature. Even the internal residues of the SARS-CoV-2 S protein are not accessible to solvent at high temperature. Furthermore, there was no unfolding of SARS-CoV-2 spike S reported at higher temperature. The most stable conformations of the SARS-CoV-2 S protein were reported at 20°C, but the free energy minimum region of the SARS-CoV-2 S protein was sharper at 40°C than other temperatures. Our findings revealed that higher temperatures have little or no influence on the stability and folding of the SARS-CoV-2 S protein.

AMADAR: a python-based package for large scale prediction of Diels-Alder transition state geometries and IRC path analysis.

Predicting transition state geometries is one of the most challenging tasks in computational chemistry, which often requires expert-based knowledge and permanent human intervention. This short communication reports technical details and preliminary results of a python-based tool (AMADAR) designed to generate any Diels-Alder (DA) transition state geometry (TS) and analyze determined IRC paths in a (quasi-)automated fashion, given the product SMILES. Two modules of the package are devoted to performing, from IRC paths, reaction force analyses (RFA) and atomic (fragment) decompositions of the reaction force F and reaction force constant [Formula: see text]. The performance of the protocol has been assessed using a dataset of 2000 DA cycloadducts retrieved from the ZINC database. The sequential location of the corresponding TSs was achieved with a success rate of 95%. RFA plots confirmed the reaction force constant [Formula: see text] to be a good indicator of the (non)synchronicity of the associated DA reactions. Moreover, the atomic decomposition of [Formula: see text] allows for the rationalization of the (a)synchronicity of each DA reaction in terms of contributions stemming from pairs of interacting atoms. The source code of the AMADAR tool is available on GitHub [ CMCDD/AMADAR(github.com) ] and can be used directly with minor customizations, mostly regarding the local working environment of the user.

Solvent promoted tautomerism in thione-containing tetraazatricyclics: evidence from 1H NMR spectroscopy and transition state studies.

Tautomerism in the nitro substituted thione-containing traazatricyclics has been investigated. Evidence from 1H NMR indicating the existence of the tautomers has been augmented with computational studies providing evidence of the stability or otherwise of these tautomers. The role of water and DMSO in the formation of the tautomers has been explained. The role of the nitro group in assisting in the formation of the tautomers has been discussed.

A multiscale ONIOM study of the buckminsterfullerene (C60) Diels-Alder reaction: from model design to reaction path analysis.

The hybrid ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) formalism is employed to investigate the Diels-Alder reaction of the buckminsterfullerene C60. Our computations suggest that the ONIOM2(M06-2X/6-31G(d): SVWN/STO-3G) model, enclosing both the diene and the pyracyclene fragment of C60 in the higher-layer, provides a reasonable trade-off between accuracy and computational cost as it comes to predicting reaction energetics. Moreover, the frontier molecular orbital (FMO) theory and activation strain model (ASM) are jointly relied on to rationalize the effect of -OH and -CN substituents on the activation barrier of this reaction. Finally, reaction paths are scrutinized to get insight into the various forces underpinning the process of cycloadduct formation.

Synthesis and evaluation of phosphonated N-heteroarylcarboxamides as DOXP-reductoisomerase (DXR) inhibitors.

The diethyl esters and disodium salts of a range of heteroarylcarbamoylphosphonic acids have been prepared and evaluated as analogues of the highly active DOXP-reductoisomerase (DXR) inhibitor, fosmidomycin. Computer-simulated docking studies, Saturation Transfer Difference (STD) NMR analysis and enzyme inhibition assays have been used to explore enzyme-binding and -inhibition potential, while in silico analysis of the DXR active site has highlighted the importance of including a well-parameterised metal co-factor in docking studies and has revealed the availability of an additional binding pocket to guide future drug design.

Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through proteome scale computational drug discovery and in vitro assay.

Malaria elimination can benefit from time and cost-efficient approaches for antimalarials such as drug repurposing. In this work, 796 DrugBank compounds were screened against 36 Plasmodium falciparum targets using QuickVina-W. Hits were selected after rescoring using GRaph Interaction Matching (GRIM) and ligand efficiency metrics: surface efficiency index (SEI), binding efficiency index (BEI) and lipophilic efficiency (LipE). They were further evaluated in Molecular dynamics (MD). Twenty-five protein-ligand complexes were finally retained from the 28,656 (36 × 796) dockings. Hit GRIM scores (0.58 to 0.78) showed their molecular interaction similarity to co-crystallized ligands. Minimum LipE (3), SEI (23) and BEI (7) were in at least acceptable thresholds for hits. Binding energies ranged from -6 to -11 kcal/mol. Ligands showed stability in MD simulation with good hydrogen bonding and favorable protein-ligand interactions energy (the poorest being -140.12 kcal/mol). In vitro testing showed 4 active compounds with two having IC50 values in the single-digit µM range.