An approach to developing a force field for molecular simulation of martensitic phase transitions between phases with subtle differences in energy and structure

Discovery of novel and highly potential inhibitors of glycogen synthase kinase 3-beta (GSK-3β) through structure-based pharmacophore modeling, virtual computational screening, docking and in silico ADMET analysis

The protein Glycogen Synthase Kinase 3-Beta (GSK-3β), is a promising therapeutic target for treating various diseases such as neurodegenerative disorders, diabetes, inflammation and cancer. This study aims to investigate the potential of compounds targeting inflammation or carbohydrate metabolism to selectively inhibit GSK3β by binding to its ATP site. To achieve this goal, we filtered a database of 49367 molecules involved in carbohydrate metabolism or targeting inflammation using various computational analyses, including pharmacophore modeling, molecular docking, dynamic simulation, prime MM-GBSA calculation, and in silico ADME studies. We generated a pharmacophore model (hypo S: AADDHRR) using two different crystallographic complexes of GSK3β and evaluated the model's performance in identifying hits using various parameters, including EF, GH, ROC, AUC and BEDROC. Subsequently, we performed various dockings (HTVS, SP, XP and IFD) for the retrieved hits and found that, 5 out of the top 10 ranked compounds had the scaffold of pyrazolidine 3,5-dione, which has never been reported to inhibit kinases. We also conducted ADMET studies to and concluded that compound N6 exhibited the best pharmacokinetic profile passing the blood-brain barrier, possessing high lipophilicity and a high coefficient of skin permeability in the intestines, along with good bioavailability and low toxicity risk assessment. Dynamic simulation were also performed indicating that compounds N6 derived from pyrazolidine 3,5-dione demonstrated better binding potential for GSK3β during the simulation period. Therefore, we propose that compounds derived from pyrazolidine-3,5-dione, which modulate the activity of lysosomal alpha-glucosidase could serve as a novel scaffold for the selective inhibition of GSK-3β.Communicated by Ramaswamy H. Sarma.

Employing comparative QSAR techniques for the recognition of dibenzofuran and dibenzothiophene derivatives toward MMP-12 inhibition

Among various matrix metalloproteinases (MMPs), MMP-12 is one of the potential targets for cancer and other diseases. However, none of the MMP-12 inhibitors has passed the clinical trials to date. Therefore, designing potential MMP-12 inhibitors as new drug molecules can provide effective therapeutic strategies for several diseases. In this study, a series of dibenzofuran and dibenzothiophene derivatives were subjected to different 2D and 3D-QSAR techniques to point out the crucial structural contributions highly influential toward the MMP-12 inhibitory activity. These techniques identified some structural attributes of these compounds that are responsible for influencing their MMP-12 inhibition. The carboxylic group may enhance proper binding with catalytic Zn2+ ion at the MMP-12 active site. Again, the i-propyl sulfonamido carboxylic acid function contributed positively toward MMP-12 inhibition. Moreover, the dibenzofuran moiety conferred stable binding at the S1' pocket for higher MMP-12 inhibition. The steric and hydrophobic groups were found favourable near the furan ring substituted at the dibenzofuran moiety. Besides these ligand-based approaches, molecular docking and molecular dynamic (MD) simulation studies not only elucidated the importance of several aspects of these MMP-12 inhibitors while disclosing the significance of the finding of these QSAR studies and their influences toward MMP-12 inhibition. The MD simulation study also revealed stable and compact binding between such compounds at the MMP-12 active site. Therefore, the findings of these validated ligand-based and structure-based molecular modeling studies can aid the development of selective and potent lead molecules that can be used for the treatment of MMP-12-associated diseases.Communicated by Ramaswamy H. Sarma.

Data Mining and in Silico Analysis of Ethiopian Traditional Medicine: Unveiling the Therapeutic Potential of Rumex abyssinicus Jacq

Multicomponent traditional medicine prescriptions are widely used in Ethiopia for disease treatment. However, inconsistencies across practitioners, cultures, and locations have hindered the development of reliable therapeutic medicines. Systematic analysis of traditional medicine data is crucial for identifying consistent and reliable medicinal materials. In this study, we compiled and analyzed a dataset of 505 prescriptions, encompassing 567 medicinal materials used for treating 106 diseases. Using association rule mining, we identified significant associations between diseases and medicinal materials. Notably, wound healing-the most frequently treated condition-was strongly associated with Rumex abyssinicus Jacq., showing a high support value. This association led to further in silico and network analysis of R. abyssinicus Jacq. compounds, revealing 756 therapeutic targets enriched in various KEGG pathways and biological processes. The Random-Walk with Restart (RWR) algorithm applied to the CODA PPI network identified these targets as linked to diseases such as cancer, inflammation, and metabolic, immune, respiratory, and neurological disorders. Many hub target genes from the PPI network were also directly associated with wound healing, supporting the traditional use of R. abyssinicus Jacq. for treating wounds. In conclusion, this study uncovers significant associations between diseases and medicinal materials in Ethiopian traditional medicine, emphasizing the therapeutic potential of R. abyssinicus Jacq. These findings provide a foundation for further research, including in vitro and in vivo studies, to explore and validate the efficacy of traditional and natural product-derived medicines.

Network pharmacology, molecular docking, and dynamics analyses to predict the antiviral activity of ginger constituents against coronavirus infection

COVID-19 is a global pandemic that caused a dramatic loss of human life worldwide, leading to accelerated research for antiviral drug discovery. Herbal medicine is one of the most commonly used alternative medicine for the prevention and treatment of many conditions including respiratory system diseases. In this study, a computational pipeline was employed, including network pharmacology, molecular docking simulations, and molecular dynamics simulations, to analyze the common phytochemicals of ginger rhizomes and identify candidate constituents as viral inhibitors. Furthermore, experimental assays were performed to analyze the volatile and non-volatile compounds of ginger and to assess the antiviral activity of ginger oil and hydroalcoholic extract. Network pharmacology analysis showed that ginger compounds target human genes that are involved in related cellular processes to the viral infection. Docking analysis highlighted five pungent compounds and zingiberenol as potential inhibitors for the main protease (Mpro), spike receptor-binding domain (RBD), and human angiotensin-converting enzyme 2 (ACE2). Then, (6)-gingerdiacetate was selected for molecular dynamics (MD) simulations as it exhibited the best binding interactions and free energies over the three target proteins. Trajectories analysis of the three complexes showed that RBD and ACE2 complexes with the ligand preserved similar patterns of root mean square deviation (RMSD) and radius of gyration (Rg) values to their respective native structures. Finally, experimental validation of the ginger hydroalcoholic extract confirmed the existence of (6)-gingerdiacetate and revealed the strong antiviral activity of the hydroalcoholic extract with IC50 of 2.727 μ g / ml . Our study provides insights into the potential antiviral activity of (6)-gingerdiacetate that may enhance the host immune response and block RBD binding to ACE2, thereby, inhibiting SARS-CoV-2 infection.

A combined in vitro and molecular dynamics simulation studies unveil the molecular basis of the anticancer potential of piperine targeting AKT1 against prostate cancer

The present study investigates the activity of the natural compound piperine on prostate cancer cell line (PC-3), followed by exploring its mechanistic inhibition on the RAC-alpha serine/threonine-protein kinase (AKT1) protein. The 3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide (MTT) assay showed that after 24 hrs of exposure to piperine (15 µmol/ml), cell viability fell to 50% compared to the standard drug flutamide (SDF) (51 µmol/ml) with a lower IC50 concentration. However, the Dual acridine orange/ethidium bromide (AO/EtBr) staining demonstrated that, as compared to the SDF, piperine caused substantial cellular death in PC-3 cells, presumably by triggering DNA fragmentation. In addition, compared to untreated cells, the proportion of the sub-G0/G1 and G2/M stages population increased considerably in piperine-treated cells. The cell cycle's sub-G0/G1 and G2/M phases were also arrested in piperine-treated cells compared to the SDF in cell cycle analysis. Based on our systems pharmacology and molecular docking studies, AKT1 is predicted as a potential target against piperine. The complementary charge between AKT1 and piperine was emphasized in the transient ligand-protein binding interaction in molecular dynamic modeling over 100 ns, and stable hydrogen bond interaction between Lys268 and Ser205 amino acid residues of the active pocket was hypothesized. Overall, the findings from our in vitro and MD simulations provide insights into the mechanism of piperine targeting AKT1 and offer a possible candidate for future prostate cancer therapeutic development.Communicated by Ramaswamy H. Sarma.

Antiurease Activity of Antibiotics: In Vitro, In Silico, Structure Activity Relationship, and MD Simulations of Cephalosporins and Fluoroquinolones

Helicobacter pylori infection is widespread in 50% of the world's population and is associated with gastric ulcers and related disorders that ultimately culminate in gastric cancer. Levofloxacin-based, or clarithromycin-based, triple therapy is frequently used to inhibit the bacterial urease enzyme for the eradication of H. pylori. A comprehensive investigation based on the urease inhibitory profiles of antibiotics and their computational implications is lacking in the scientific literature. The present study was aimed specifically to determine the antiurease activities within the realms of cephalosporins and fluoroquinolones by in vitro methods supported with in silico investigations. The results demonstrate the jack bean urease inhibitory activity of cephalosporins, wherein cefadroxil, cefpodoxime, cefotaxime, and cefaclor displayed inhibitions (IC50 21.35 ± 0.64 to 62.86 ± 0.78 μM) compared with the standard thiourea (IC50 21.25 ± 0.15 μM). Among fluoroquinolones, levofloxacin, ofloxacin, and gemifloxacin (IC50 7.24 ± 0.29 to 16.53 ± 0.85 μM) unveiled remarkable inhibitory profiles. Levofloxacin and ofloxacin exhibited competitive inhibition against the said enzyme. Ciprofloxacin and moxifloxacin displayed weak urease inhibitions. During molecular docking studies, Asp362, Gly279, Arg338, Asn168, Asp223, Gln364, and Met366 were involved in hydrogen bonding in fluoroquinolones, and hydrogen bonding was established with Arg338, His248, Asn168 residues, and metal Ni601 and Ni602 of the enzyme. MD simulations and MMPBSA results demonstrated the existence of significant protein-ligand binding. Overall, these results warrant further investigations into the significance of these active molecules in relation to their inhibitory potential against the targeted urease enzyme.

Anti-enzymatic and DNA docking studies of montelukast: A multifaceted molecular scaffold with in vitro investigations, molecular expression analysis and molecular dynamics simulations

Montelukast, an approved leukotriene receptor 1 (Cys-LT 1) antagonist with anti-inflammatory properties is used for the treatment of asthma and allergic rhinitis. In the present studies, montelukast was subjected to in vitro inhibitory assays followed by kinetic and in silico investigations. Montelukast demonstrated inhibitory activity against yeast α-glucosidase (IC50 44.31 ± 1.21 μM), jack bean urease (JB urease, IC50 8.72 ± 0.23 μM), human placental alkaline phosphatase (hPAP, IC50 17.53 ± 0.19 μM), bovine intestinal alkaline phosphatase (bIAP, IC50 15.18 ± 0.23 μM) and soybean 15-lipoxygenase (15-LOX, IC50 2.41 ± 0.13 μM). Kinetic studies against α-glucosidase and urease enzymes revealed its competitive mode of inhibition. Molecular expression analysis of montelukast in breast cancer cell line MCF-7 down-regulated AP by a factor of 0.27 (5 μM) compared with the 0.26 value for standard inhibitor levamisole (10 μM). Molecular docking estimated a binding affinity ranging -8.82 to -15.65 kcal/mol for the enzymes. Docking against the DNA dodecamer (ID: 1BNA) observed -9.13 kcal/mol via minor groove binding. MD simulations suggested stable binding between montelukast and the target proteins predicting strong inhibitory potential of the ligand. Montelukast features a chloroquinoline, phenyl ring, a cyclopropane group, a carboxylic group and a sulfur atom all of which collectively enhance its inhibitory potential against the said enzymes. These in vitro and computational investigations demonstrate that it is possible and suggested that the interactions of montelukast with more than one targets presented herein may be linked with the side effects presented by this drug and necessitate additional work. The results altogether suggest montelukast as an important structural scaffold possessing multitargeted features and warrant further investigations in repurposing beyond its traditional pharmacological use.

Finite temperature string by K-means clustering sampling with order parameters as collective variables for molecular crystals: application to polymorphic transformation between β-CL-20 and ε-CL-20

Polymorphic transformation of molecular crystals is a fundamental phase transition process, and it is important practically in the chemical, material, biopharmaceutical, and energy storage industries. However, understanding of the transformation mechanism at the molecular level is poor due to the extreme simulating challenges in enhanced sampling and formulating order parameters (OPs) as the collective variables that can distinguish polymorphs with quite similar and complicated structures so as to describe the reaction coordinate. In this work, two kinds of OPs for CL-20 were constructed by the bond distances, bond orientations and relative orientations. A K-means clustering algorithm based on the Euclidean distance and sample weight was used to smooth the initial finite temperature string (FTS), and the minimum free energy path connecting β-CL-20 and ε-CL-20 was sketched by the string method in collective variables, and the free energy profile along the path and the nucleation kinetics were obtained by Markovian milestoning with Voronoi tessellations. In comparison with the average-based sampling, the K-means clustering algorithm provided an improved convergence rate of FTS. The simulation of transformation was independent of OP types but was affected greatly by finite-size effects. A surface-mediated local nucleation mechanism was confirmed and the configuration located at the shoulder of potential of mean force, rather than overall maximum, was confirmed to be the critical nucleus formed by the cooperative effect of the intermolecular interactions. This work provides an effective way to explore the polymorphic transformation of caged molecular crystals at the molecular level.

In Silico and In Vitro Analyses to Repurpose Quercetin as a Human Pancreatic α-Amylase Inhibitor

Human pancreatic α-amylase (HPA), situated at the apex of the starch digestion hierarchy, is an attractive therapeutic approach to precisely regulate blood glucose levels, thereby efficiently managing diabetes. Polyphenols offer a natural and multifaceted approach to moderate postprandial sugar spikes, with their slight modulation in carbohydrate digestion and potential secondary benefits, such as antioxidant and anti-inflammatory effects. Taking into consideration the unfavorable side effects of currently available commercial medications, we aimed to study a library of polyphenols attributed to their remarkable antidiabetic properties and screened the most potent HPA inhibitor via a comprehensive in silico study encompassing molecular docking, molecular mechanics with generalized Born and surface area solvation (MM/GBSA) calculation, molecular dynamics (MD) simulation, density functional theory (DFT) study, and pharmacokinetic properties followed by an in vitro assay. Significant hydrogen bonding with the catalytic triad residues of HPA, prominent MM/GBSA binding energy of -27.03 kcal/mol, and the stable nature of the protein-ligand complex with regard to 100 ns MD simulation screened quercetin as the best HPA inhibitor. Additionally, quercetin showed strong reactivity in the substrate-binding pocket of HPA and exhibited favorable pharmacokinetic properties with a considerable inhibitory concentration (IC50) of 57.37 ± 0.9 μg/mL against α-amylase. This study holds prospects for HPA inhibition and suggests quercetin as an approach to therapy for diabetes; however, it is imperative to conduct further research.

Exploring STAT3 stimulatory potential of novel wound healing molecules by virtual screening and molecular dynamics simulations

STAT3 signaling is a major intrinsic pathway for cell proliferation owing to its frequent activation in injured tissues. Various STAT3-regulated genes encode cytokines and growth factors, the receptors of which in turn activate the same STAT3 pathways, thereby regulating cell proliferation. In present study, we aimed to analyze several compounds for their wound healing and tissue repair potential by computer-aided virtual screening and Molecular dynamics (MD) simulation. Based on literature studies, a total of 36 drug molecules were selected having critical functions in wound healing and tissue repair. The pharmacological features (ADME and toxicity) of these molecules were predicted to find lead molecules among them. Further, a comparative study was performed to screen binding efficiency of STAT3 with many conventional wound healers by molecular docking. Among all, W6S, Strychnin, Prednisone and N-(6-(4-(3-(4-((4-Methylpiperazin-1-yl) methyl)-3- (trifluoromethyl)phenyl)ureido)phenoxy)pyrimidin-4-yl)cyclopropanecarboxamide showed best docking with STAT3 protein. The calculated binding energy of these molecules with STAT3 was found to be -8.9 Kca/mol for N-(6-(4-(3-(4-((4-Methylpiperazin-1-yl) methyl)-3-(trifluoromethyl) phenyl)ureido)phenoxy)pyrimidin-4-yl)cyclopropanecarboxamide, -8.7 Kcal/mol for W6S, -8.5 Kcal/mol for Strychnine and -8.4 Kcal/mol for Prednisone . The result was reconsidered for MD simulation. The simulation result showed stable binding of the ligand with STAT3 protein for 100 ns. These compounds showed better interaction potential with STAT3 was compared to known tissue repair molecules. Our data paves way for further exploration of these molecules as novel cell proliferators to be tested in various types of wound and tissue injuries.Communicated by Ramaswamy H. Sarma.

Deciphering Molecular Aspects of Potential α-Glucosidase Inhibitors within Aspergillus terreus: A Computational Odyssey of Molecular Docking-Coupled Dynamics Simulations and Pharmacokinetic Profiling

Hyperglycemia, as a hallmark of the metabolic malady diabetes mellitus, has been an overwhelming healthcare burden owing to its high rates of comorbidity and mortality, as well as prospective complications affecting different body organs. Available therapeutic agents, with α-glucosidase inhibitors as one of their cornerstone arsenal, control stages of broad glycemia while showing definitive characteristics related to their low clinical efficiency and off-target complications. This has propelled the academia and industrial section into discovering novel and safer candidates. Herein, we provided a thorough computational exploration of identifying candidates from the marine-derived Aspergillus terreus isolates. Combined structural- and ligand-based approaches using a chemical library of 275 metabolites were adopted for pinpointing promising α-glucosidase inhibitors, as well as providing guiding insights for further lead optimization and development. Structure-based virtual screening through escalating precision molecular docking protocol at the α-glucosidase canonical pocket identified 11 promising top-docked hits, with several being superior to the market drug reference, acarbose. Comprehensive ligand-based investigations of these hits' pharmacokinetics ADME profiles, physiochemical characterizations, and obedience to the gold standard Lipinski's rule of five, as well as toxicity and mutagenicity profiling, proceeded. Under explicit conditions, a molecular dynamics simulation identified the top-stable metabolites: butyrolactone VI (SK-44), aspulvinone E (SK-55), butyrolactone I 4''''-sulfate (SK-72), and terrelumamide B (SK-173). They depicted the highest free binding energies and steadiest thermodynamic behavior. Moreover, great structural insights have been revealed, including the advent of an aromatic scaffold-based interaction for ligand-target complex stability. The significance of introducing balanced hydrophobic/polar moieties, like triazole and other bioisosteres of carboxylic acid, has been highlighted across docking, ADME/Tox profiling, and molecular dynamics studies for maximizing binding interactions while assuring safety and optimal pharmacokinetics for targeting the intestinal-localized α-glucosidase enzyme. Overall, this study provided valuable starting points for developing new α-glucosidase inhibitors based on nature-derived unique scaffolds, as well as guidance for prospective lead optimization and development within future pre-clinical and clinical investigations.

Evaluating chemotherapeutic potential of soya-isoflavonoids against high penetrance genes in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive molecular subtype of breast cancer (BC) associated with a poor prognosis. Owing to the structural similarity with 17-β-estradiol, consumption of soya-isoflavonoids are associated with a reduced rate of hormone-receptive BC incidence, but their role in TNBC is not deciphered in detail. This present study thus aims to investigate the therapeutic binding dynamics of dietary soya-flavonoids with the six high penetrance (HP) receptors in TNBC, viz. BRCA1, BRCA2, PALB2, PTEN, STK11 and TP53. Out of the 14 soya-flavonoids screened based on ADMET descriptors and several other physicochemical, bioavailability, drug and lead-likeness properties, four hits were shortlisted (Daidzein, Genistein, Glycitein and Biochanin A). Docking and molecular dynamics (MD) simulation revealed Genistein as the most potential multi-target inhibitor of the six TNBC HP genes. Additionally, Genistein exhibited excellent binding specificity with PTEN, a potent mediator of the PI3K signaling pathway in TNBC. The binding interaction of PTEN and Genistein was further compared against a standardized FDA-approved chemotherapeutic inhibitor, Olaparib, computed through various MD trajectory analysis, principal component analysis and computation of free energy landscape. This study reveals a comparatively better binding dynamics of PTEN-Genistein than PTEN-Olaparib. With a significant global surge in biomarker-based precision therapeutics in oncology, the results of this exhaustive in-silico study thus encourage the prospect of validating PTEN as a druggable target of Genistein, a unique drug-receptor combination in the future.Communicated by Ramaswamy H. Sarma.

Picloram binds to the h1 and h4 helices of HSA domain IIIA at drug binding site 2

Picloram (PC) is a systemic herbicide that controls herbaceous weeds and woody plants. HSA, the most abundant protein in human physiology, binds to all exogenic and endogenic ligands. PC is a stable molecule (t_1/2~157-513 days) and a potential threat to human health via the food chain. HSA and PC binding study has been done to decipher the location and thermodynamics of binding. It has been studied with prediction tools like autodocking and MD simulation and then confirmed with fluorescence spectroscopy. HSA fluorescence was quenched by PC at pH 7.4 (N state), pH 3.5 (F state), and pH 7.4 with 4.5 M urea (I state) at temperatures 283 K, 297 K, and 303 K. The location of binding was found to be interdomain between II and III which overlaps with drug binding site 2. The binding was spontaneous, and entropy-driven that show a noticeable increase in binding with the increase in temperature. No secondary structure change at the native state has been observed due to binding. The binding results are important to understand the physiological assimilation of PC. In silico predictions and the results of spectroscopic studies unambiguously indicate the locus and nature of the binding.

Integration of Chemoinformatics and Multi-Omics Analysis Defines ECT2 as a Potential Target for Cancer Drug Therapy

Epithelial cell transforming 2 (ECT2) is a potential oncogene and a number of recent studies have correlated it with the progression of several human cancers. Despite this elevated attention for ECT2 in oncology-related reports, there is no collective study to combine and integrate the expression and oncogenic behavior of ECT2 in a panel of human cancers. The current study started with a differential expression analysis of ECT2 in cancerous versus normal tissue. Following that, the study asked for the correlation between ECT2 upregulation and tumor stage, grade, and metastasis, along with its effect on patient survival. Moreover, the methylation and phosphorylation status of ECT2 in tumor versus normal tissue was assessed, in addition to the investigation of the ECT2 effect on the immune cell infiltration in the tumor microenvironment. The current study revealed that ECT2 was upregulated as mRNA and protein levels in a list of human tumors, a feature that allowed for the increased filtration of myeloid-derived suppressor cells (MDSC) and decreased the level of natural killer T (NKT) cells, which ultimately led to a poor prognosis survival. Lastly, we screened for several drugs that could inhibit ECT2 and act as antitumor agents. Collectively, this study nominated ECT2 as a prognostic and immunological biomarker, with reported inhibitors that represent potential antitumor drugs.

In silico investigation of binding propensity of hematoxylin derivative and damnacanthal for their potential inhibitory effect on HIV-1 Vpr from different subtypes

HIV-1, the causative agent of AIDS leads to many deaths worldwide though few options are available as therapeutics. To deal with the continuous mutation in the virus genome, requirement of new drugs is always there. Subtype variation plays a crucial role in case of HIV-1 therapeutics development. In this study, we want to investigate some pre examined molecules that can be effective for HIV-1 VPR. Inhibition of several protein-protein interactions with the small molecules will lead to identify some molecules as therapeutics other than the conventional drugs. We retrieved the sequences of different subtypes from the database and representative sequences were identified. Representative structures were modelled and validated using MD simulations. Forty molecules, showing anti Vpr activity in vitro were identified from literature survey and those were docked with each subtype representative structures. Two molecules a stable Hematoxylin Derivative (SHD) and Damnacanthal (D3), these were shown to be bind more effectively for all the subtypes. The stability of the protein and those two small molecule complexes were identified again with MD simulation followed by the binding energy calculation. Thus, these molecules can be thought as any option other than the conventional drug targeting HIV-1 Vpr.Communicated by Ramaswamy H. Sarma.

Molecular and Biological Investigation of Isolated Marine Fungal Metabolites as Anticancer Agents: A Multi-Target Approach

Cancer is the leading cause of death globally, with an increasing number of cases being annually reported. Nature-derived metabolites have been widely studied for their potential programmed necrosis, cytotoxicity, and anti-proliferation leading to enrichment for the modern medicine, particularly within the last couple of decades. At a more rapid pace, the concept of multi-target agents has evolved from being an innovative approach into a regular drug development procedure for hampering the multi-fashioned pathophysiology and high-resistance nature of cancer cells. With the advent of the Red Sea Penicillium chrysogenum strain S003-isolated indole-based alkaloids, we thoroughly investigated the molecular aspects for three major metabolites: meleagrin (MEL), roquefortine C (ROC), and isoroquefortine C (ISO) against three cancer-associated biological targets Cdc-25A, PTP-1B, and c-Met kinase. The study presented, for the first time, the detailed molecular insights and near-physiological affinity for these marine indole alkaloids against the assign targets through molecular docking-coupled all-atom dynamic simulation analysis. Findings highlighted the superiority of MEL's binding affinity/stability being quite in concordance with the in vitro anticancer activity profile conducted via sulforhodamine B bioassay on different cancerous cell lines reaching down to low micromolar or even nanomolar potencies. The advent of lengthy structural topologies via the metabolites' extended tetracyclic cores and aromatic imidazole arm permitted multi-pocket accommodation addressing the selectivity concerns. Additionally, the presence decorating polar functionalities on the core hydrophobic tetracyclic ring contributed compound's pharmacodynamic preferentiality. Introducing ionizable functionality with more lipophilic characters was highlighted to improve binding affinities which was also in concordance with the conducted drug-likeness/pharmacokinetic profiling for obtaining a balanced pharmacokinetic/dynamic profile. Our study adds to the knowledge regarding drug development and optimization of marine-isolated indole-based alkaloids for future iterative synthesis and pre-clinical investigations as multi-target anticancer agents.

In silico activity and ADMET profiling of phytochemicals from Ethiopian indigenous aloes using pharmacophore models

In silico profiling is used in identification of active compounds and guide rational use of traditional medicines. Previous studies on Ethiopian indigenous aloes focused on documentation of phytochemical compositions and traditional uses. In this study, ADMET and drug-likeness properties of phytochemicals from Ethiopian indigenous aloes were evaluated, and pharmacophore-based profiling was done using Discovery Studio to predict therapeutic targets. The targets were examined using KEGG pathway, gene ontology and network analysis. Using random-walk with restart algorithm, network propagation was performed in CODA network to find diseases associated with the targets. As a result, 82 human targets were predicted and found to be involved in several molecular functions and biological processes. The targets also were linked to various cancers and diseases of immune system, metabolism, neurological system, musculoskeletal system, digestive system, hematologic, infectious, mouth and dental, and congenital disorder of metabolism. 207 KEGG pathways were enriched with the targets, and the main pathways were metabolism of steroid hormone biosynthesis, lipid and atherosclerosis, chemical carcinogenesis, and pathways in cancer. In conclusion, in silico target fishing and network analysis revealed therapeutic activities of the phytochemicals, demonstrating that Ethiopian indigenous aloes exhibit polypharmacology effects on numerous genes and signaling pathways linked to many diseases.

Bio-Guided Isolation of SARS-CoV-2 Main Protease Inhibitors from Medicinal Plants: In Vitro Assay and Molecular Dynamics

Since the emergence of the pandemic of the coronavirus disease (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the discovery of antiviral phytoconstituents from medicinal plants against SARS-CoV-2 has been comprehensively researched. In this study, thirty-three plants belonging to seventeen different families used traditionally in Saudi Arabia were tested in vitro for their ability to inhibit the SARS-CoV-2 main protease (M^PRO). Major constituents of the bio-active extracts were isolated and tested for their inhibition potential against this enzyme; in addition, their antiviral activity against the SARS-CoV-2 Egyptian strain was assessed. Further, the thermodynamic stability of the best active compounds was studied through focused comparative insights for the active metabolites regarding ligand–target binding characteristics at the molecular level. Additionally, the obtained computational findings provided useful directions for future drug optimization and development. The results revealed that Psiadia punctulata, Aframomum melegueta, and Nigella sativa extracts showed a high percentage of inhibition of 66.4, 58.7, and 31.5%, against SARS-CoV-2 M^PRO, respectively. The major isolated constituents of these plants were identified as gardenins A and B (from P. punctulata), 6-gingerol and 6-paradol (from A. melegueta), and thymoquinone (from N. sativa). These compounds are the first to be tested invitro against SARS-CoV-2 M^PRO. Among the isolated compounds, only thymoquinone (THY), gardenin A (GDA), 6-gingerol (GNG), and 6-paradol (PAD) inhibited the SARS-CoV-2 M^PRO enzyme with inhibition percentages of 63.21, 73.80, 65.2, and 71.8%, respectively. In vitro assessment of SARS-CoV-2 (hCoV-19/Egypt/NRC-03/2020 (accession number on GSAID: EPI_ISL_430820) revealed a strong-to-low antiviral activity of the isolated compounds. THY showed relatively high cytotoxicity and was anti-SARS-CoV-2, while PAD demonstrated a cytotoxic effect on the tested VERO cells with a selectivity index of CC₅₀/IC₅₀ = 1.33 and CC₅₀/IC₅₀ = 0.6, respectively. Moreover, GNG had moderate activity at non-cytotoxic concentrations in vitro with a selectivity index of CC₅₀/IC₅₀ = 101.3/43.45 = 2.3. Meanwhile, GDA showed weak activity with a selectivity index of CC₅₀/IC₅₀ = 246.5/83.77 = 2.9. The thermodynamic stability of top-active compounds revealed preferential stability and SARS-CoV-2 M^PRO binding affinity for PAD through molecular-docking-coupled molecular dynamics simulation. The obtained results suggest the treating potential of these plants and/or their active metabolites for COVID-19. However, further in-vivo and clinical investigations are required to establish the potential preventive and treatment effectiveness of these plants and/or their bio-active compounds in COVID-19.

Discovery of New 1,3,4-Oxadiazoles with Dual Activity Targeting the Cholinergic Pathway as Effective Anti-Alzheimer Agents

Finding an effective anti-Alzheimer agent is quite challenging due to its multifactorial nature. As such, multitarget directed ligands (MTDLs) could be a promising paradigm for finding potential therapeutically effective new small-molecule bioactive agents against Alzheimer's disease (AD). We herein present the design, synthesis, and biological evaluation of a new series of compounds based on a 5-pyrid-3-yl-1,3,4-oxadiazole scaffold. Our synthesized compounds displayed excellent in vitro enzyme inhibitory activity at nanomolar (nM) concentrations against two major AD disease-modifying targets, i.e., acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). Among our compounds, 5e was considered the best dual inhibitor of both AChE (IC₅₀ = 50.87 nM) and BuChE (IC₅₀ = 4.77 nM), where these values surpassed those of rivastagmine (the only FDA-approved dual AChE and BuChE inhibitor) in our study. Furthermore, in vivo and ex vivo testing of the hit compound 5e highlighted its significant AD-biotargeting effects including reducing the elevated levels of lipid peroxidation and glutathione (GSH), normalizing levels of 8-OHdG, and, most importantly, decreasing the levels of the well-known AD hallmark β-amyloid protein. Finally, the binding ability of 5e to each of our targets, AChE and BuChE, was confirmed through additional molecular docking and molecular dynamics (MD) simulations that reflected good interactions of 5e to the active site of both targets. Hence, we herein present a series of new 1,3,4-oxadiazoles that are promising leads for the development of dual-acting AChE and BuChE inhibitors for the management of AD.

Molecular basis of the potential interaction of SARS-CoV-2 spike protein to CD147 in COVID-19 associated-lymphopenia

Lymphopenia is considered one of the most characteristic clinical features of the coronavirus disease 2019 (COVID-19). SARS-CoV-2 infects host cells via the interaction of its spike protein with the human angiotensin-converting enzyme 2 (hACE2) receptor. Since T lymphocytes display a very low expression level of hACE2, a novel receptor might be involved in the entry of SARS-CoV-2 into T cells. The transmembrane glycoprotein CD147 is highly expressed by activated T lymphocytes, and was recently proposed as a probable route for SARS-CoV-2 invasion. To understand the molecular basis of the potential interaction of SARS-CoV-2 to CD147, we have investigated the binding of the viral spike protein to this receptor in-silico. The results showed that this binding is dominated by electrostatic interactions involving residues Arg403, Asn481, and the backbone of Gly502. The overall binding arrangement shows the CD147 C-terminal domain interacting with the spike external subdomain in the grove between the short antiparallel β strands, β1' and β2', and the small helix α1'. This proposed interaction was further confirmed using MD simulation and binding free energy calculation. These data contribute to a better understanding of the mechanism of infection of SARS-CoV-2 to T lymphocytes and could provide valuable insights for the rational design of adjuvant treatment for COVID-19. Communicated by Ramaswamy H. Sarma.

In Silico and In Vitro Studies for Benzimidazole Anthelmintics Repurposing as VEGFR-2 Antagonists: Novel Mebendazole-Loaded Mixed Micelles with Enhanced Dissolution and Anticancer Activity

Cancer is a leading cause of death worldwide and its incidence is unfortunately anticipated to rise in the next years. On the other hand, vascular endothelial growth factor receptor 2 (VEGFR-2) is highly expressed in tumor-associated endothelial cells, where it affects tumor-promoting angiogenesis. Therefore, VEGFR-2 is considered one of the most promising therapeutic targets for cancer treatment. Furthermore, some FDA-approved benzimidazole anthelmintics have already shown potential anticancer activities. Therefore, repurposing them against VEGFR-2 can provide a rapid and effective alternative that can be implicated safely for cancer treatment. Hence, 13 benzimidazole anthelmintic drugs were subjected to molecular docking against the VEGFR-2 receptor. Among the tested compounds, fenbendazole (FBZ, 1), mebendazole (MBZ, 2), and albendazole (ABZ, 3) were proposed as potential VEGFR-2 antagonists. Furthermore, molecular dynamics simulations were carried out at 200 ns, giving more information on their thermodynamic and dynamic properties. Besides, the anticancer activity of the aforementioned drugs was tested in vitro against three different cancer cell lines, including liver cancer (HUH7), lung cancer (A549), and breast cancer (MCF7) cell lines. The results depicted potential cytotoxic activity especially against both HUH7 and A549 cell lines. Furthermore, to improve the aqueous solubility of MBZ, it was formulated in the form of mixed micelles (MMs) which showed an enhanced drug release with better promising cytotoxicity results compared to the crude MBZ. Finally, an in vitro quantification for VEGFR-2 concentration in treated HUH7 cells has been conducted based on the enzyme-linked immunosorbent assay. The results disclosed that FBZ, MBZ, and ABZ significantly (p < 0.001) reduced the concentration of VEGFR-2, while the lowest inhibition was achieved in MBZ-loaded MMs, which was even much better than the reference drug sorafenib. Collectively, the investigated benzimidazole anthelmintics could be encountered as lead compounds for further structural modifications and thus better anticancer activity, and that was accomplished through studying their structure-activity relationships.

Recent advances in drug polymorphs: Aspects of pharmaceutical properties and selective crystallization

Drug polymorphism, an established term used to describe the phenomenon that a drug can exist in different crystalline phases, has attracted great interests in pharmaceutical field in consideration of its important role in affecting the pharmaceutical performance of oral formulations. This paper presents an overview of recent advances in the research on polymorphic drug systems including understandings on nucleation, crystal growth, dissolution, mechanical properties, polymorphic transformation, etc. Moreover, new strategies and mechanisms in the control of polymorphic forms are also highlighted in this review. Furthermore, challenges and trends in the development of polymorphic drugs are briefly discussed, aiming at developing effective and efficient pharmaceutical formulations containing the polymorphic drugs.

Calendulaglycoside A showing potential activity against SARS-CoV-2 main protease: Molecular docking, molecular dynamics, and SAR studies

BACKGROUND AND AIM

The discovery of drugs capable of inhibiting SARS-CoV-2 is a priority for human beings due to the severity of the global health pandemic caused by COVID-19. To this end, natural products can provide therapeutic alternatives that could be employed as an effective safe treatment for COVID-19.

EXPERIMENTAL PROCEDURE

Twelve compounds were isolated from the aerial parts of C. officinalis L. and investigated for their inhibitory activities against SARS-CoV-2 Mpro compared to its co-crystallized N3 inhibitor using molecular docking studies. Furthermore, a 100 ns MD simulation was performed for the most active two promising compounds, Calendulaglycoside A (SAP5) and Osteosaponin-I (SAP8).

RESULTS AND CONCLUSION

At first, molecular docking studies showed interesting binding scores as compared to the N3 inhibitor. Calendulaglycoside A (SAP5) achieved a superior binding than the co-crystallized inhibitor indicating promising affinity and intrinsic activity towards the Mpro of SARS-CoV-2 as well. Moreover, findings illustrated preferential stability for SAP5 within the Mpro pocket over that of N3 beyond the 40 ns MD simulation course. Structural preferentiality for triterpene-Mpro binding highlights the significant role of 17β-glucosyl and carboxylic 3α-galactosyl I moieties through high electrostatic interactions across the MD simulation trajectories. Furthermore, this study clarified a promising SAR responsible for the antiviral activity against the SARS-CoV-2 Mpro and the design of new drug candidates targeting it as well. The above findings could be promising for fast examining the previously isolated triterpenes both pre-clinically and clinically for the treatment of COVID-19.

In Silico Designing of a Multitope Vaccine against Rhizopus microsporus with Potential Activity against Other Mucormycosis Causing Fungi

During the current era of the COVID-19 pandemic, the dissemination of Mucorales has been reported globally, with elevated rates of infection in India, and because of the high rate of mortality and morbidity, designing an effective vaccine against mucormycosis is a major health priority, especially for immunocompromised patients. In the current study, we studied shared Mucorales proteins, which have been reported as virulence factors, and after analysis of several virulent proteins for their antigenicity and subcellular localization, we selected spore coat (CotH) and serine protease (SP) proteins as the targets of epitope mapping. The current study proposes a vaccine constructed based on top-ranking cytotoxic T lymphocyte (CTL), helper T lymphocyte (HTL), and B cell lymphocyte (BCL) epitopes from filtered proteins. In addition to the selected epitopes, β-defensins adjuvant and PADRE peptide were included in the constructed vaccine to improve the stimulated immune response. Computational tools were used to estimate the physicochemical and immunological features of the proposed vaccine and validate its binding with TLR-2, where the output data of these assessments potentiate the probability of the constructed vaccine to stimulate a specific immune response against mucormycosis. Here, we demonstrate the approach of potential vaccine construction and assessment through computational tools, and to the best of our knowledge, this is the first study of a proposed vaccine against mucormycosis based on the immunoinformatics approach.

β-Blockers bearing hydroxyethylamine and hydroxyethylene as potential SARS-CoV-2 Mpro inhibitors: rational based design, in silico, in vitro, and SAR studies for lead optimization

The global COVID-19 pandemic became more threatening especially after the introduction of the second and third waves with the current large expectations for a fourth one as well. This urged scientists to rapidly develop a new effective therapy to combat SARS-CoV-2. Based on the structures of β-adrenergic blockers having the same hydroxyethylamine and hydroxyethylene moieties present in the HIV-1 protease inhibitors which were found previously to inhibit the replication of SARS-CoV, we suggested that they may decrease the SARS-CoV-2 entry into the host cell through their ability to decrease the activity of RAAS and ACE2 as well. Herein, molecular docking of twenty FDA-approved β-blockers was performed targeting SARS-CoV-2 Mpro. Results showed promising inhibitory activities especially for Carvedilol (CAR) and Nebivolol (NEB) members. Moreover, these two drugs together with Bisoprolol (BIS) as an example from the lower active ones were subjected to molecular dynamics simulations at 100 ns. Great stability across the whole 100 ns timeframe was observed for the top docked ligands, CAR and NEB, over BIS. Conformational analysis of the examined drugs and hydrogen bond investigation with the pocket's crucial residues confirm the great affinity and confinement of CAR and NEB within the Mpro binding site. Moreover, the binding-free energy analysis and residue-wise contribution analysis highlight the nature of ligand–protein interaction and provide guidance for lead development and optimization. Furthermore, the examined three drugs were tested for their in vitro inhibitory activities towards SARS-CoV-2. It is worth mentioning that NEB achieved the most potential anti-SARS-CoV-2 activity with an IC₅₀ value of 0.030 μg ml⁻¹. Besides, CAR was found to have a promising inhibitory activity with an IC₅₀ of 0.350 μg ml⁻¹. Also, the IC₅₀ value of BIS was found to be as low as 15.917 μg ml⁻¹. Finally, the SARS-CoV-2 Mpro assay was performed to evaluate and confirm the inhibitory effects of the tested compounds (BIS, CAR, and NEB) towards the SARS-CoV-2 Mpro enzyme. The obtained results showed very promising SARS-CoV-2 Mpro inhibitory activities of BIS, CAR, and NEB (IC₅₀ = 118.50, 204.60, and 60.20 μg ml⁻¹, respectively) compared to lopinavir (IC₅₀ = 73.68 μg ml⁻¹) as a reference standard.

Hydroxyethylamine and hydroxyethylene moieties of β-blockers exert potential SARS-CoV-2 inhibitory effects: rational-based design and in silico, in vitro, and SAR Studies.

Marine Brominated Tyrosine Alkaloids as Promising Inhibitors of SARS-CoV-2

There have been more than 150 million confirmed cases of SARS-CoV-2 since the beginning of the pandemic in 2019. By June 2021, the mortality from such infections approached 3.9 million people. Despite the availability of a number of vaccines which provide protection against this virus, the evolution of new viral variants, inconsistent availability of the vaccine around the world, and vaccine hesitancy, in some countries, makes it unreasonable to rely on mass vaccination alone to combat this pandemic. Consequently, much effort is directed to identifying potential antiviral treatments. Marine brominated tyrosine alkaloids are recognized to have antiviral potential. We test here the antiviral capacity of fourteen marine brominated tyrosine alkaloids against five different target proteins from SARS-CoV-2, including main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H). These marine alkaloids, particularly the hexabrominated compound, fistularin-3, shows promising docking interactions with predicted binding affinities (S-score = −7.78, −7.65, −6.39, −6.28, −8.84 Kcal/mol) for the main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H), respectively, where it forms better interactions with the protein pockets than the native interaction. It also shows promising molecular dynamics, pharmacokinetics, and toxicity profiles. As such, further exploration of the antiviral properties of fistularin-3 against SARS-CoV-2 is merited.

Investigating the structure-activity relationship of marine natural polyketides as promising SARS-CoV-2 main protease inhibitors

Since its first report in December 2019, the novel coronavirus virus, SARS-CoV-2, has caused an unprecedented global health crisis and economic loss imposing a tremendous burden on the worldwide finance, healthcare system, and even daily life. Even with the introduction of different preventive vaccines, there is still a dire need for effective antiviral therapeutics. Nature has been considered as the historical trove of drug discovery and development, particularly in cases of worldwide crises. Herein, a comprehensive in silico investigation of a highly focused chemical library of 34 pederin-structurally related marine compounds, belonging to four polyketides families, was initiated against the SARS-CoV-2 main protease, Mpro, being the key replicating element of the virus and main target in many drugs development programs. Two of the most potent SARS-CoV-2 Mpro co-crystallized inhibitors, O6K and N3, were added to the tested database as reference standards. Through molecular docking simulation, promising compounds including Pederin (1), Dihydro-onnamide A (11), Onnamide C (14), Pseudo-onnamide A (17), and Theopederin G (29) have been identified from different families based on their superior ligand–protein energies and relevant binding profiles with the key Mpro pocket residues. Thermodynamic behaviors of the identified compounds were investigated through 200 ns all-atom molecular dynamics simulation illustrating their significant stability and pocket accommodation. Furthermore, structural activity preferentiality was identified for the pederin-based marine compounds highlighting the importance of the terminal guanidine and cyclic hemiacetal linker, and the length of the sidechain. Our findings highlight the challenges of targeting SARS-CoV-2 Mpro as well as recommending further in vitro and in vivo studies regarding the examined marine products either alone or in combination paving the way for promising lead molecules.

Marine natural polyketides showed promising SARS-CoV-2 main protease inhibitory activities.

Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of ACEIs Against SARS-CoV-2 Targeting the hACE2 Receptor

The rapid and global spread of a new human coronavirus, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has produced an immediate urgency to discover promising targets for the treatment of COVID-19. Here, we consider drug repurposing as an attractive approach that can facilitate the drug discovery process by repurposing existing pharmaceuticals to treat illnesses other than their primary indications. We review current information concerning the global health issue of COVID-19 including promising approved drugs, e.g., human angiotensin-converting enzyme inhibitors (hACEIs). Besides, we describe computational approaches to be used in drug repurposing and highlight examples of in-silico studies of drug development efforts against SARS-CoV-2. Alacepril and lisinopril were found to interact with human angiotensin-converting enzyme 2 (hACE2), the host entranceway for SARS-CoV-2 spike protein, through exhibiting the most acceptable rmsd_refine values and the best binding affinity through forming a strong hydrogen bond with Asn90, which is assumed to be essential for the activity, as well as significant extra interactions with other receptor-binding residues. Furthermore, molecular dynamics (MD) simulations followed by calculation of the binding free energy were also carried out for the most promising two ligand-pocket complexes from docking studies (alacepril and lisinopril) to clarify some information on their thermodynamic and dynamic properties and confirm the docking results as well. These results we obtained probably provided an excellent lead candidate for the development of therapeutic drugs against COVID-19. Eventually, animal experiments and accurate clinical trials are needed to confirm the potential preventive and treatment effect of these compounds.

Molecular Docking and Dynamics Simulation Study of Hyrtios erectus Isolated Scalarane Sesterterpenes as Potential SARS-CoV-2 Dual Target Inhibitors

Presently, the world is under the toll of pandemic coronavirus disease-2019 (COVID-19) outbreak caused by SARS-CoV-2. Lack of effective and safe therapeutics has stressed the scientific community for developing novel therapeutics capable of alleviating and stopping this pandemic. Within the presented study, molecular docking, ADME properties and all-atom molecular dynamic (MD) simulation, along with two standard antiviral agents (lopinavir and benzopurpurin-4B), were applied to investigate 15 scalaranes sesterterpenes natural compounds, purified from the Red Sea marine sponge Hyrtios erectus, as potential COVID-19 dual-target inhibitors. Following multi-step docking within COVID-19 main protease and Nsp15 endoribonuclease cavities, nine promising drug-like compounds exhibited higher docking scores as well as better interactions with the target's crucial residues than those of reference ligands. Compounds 2, 6, 11, and 15, were predicted to simultaneously subdue the activity of the two COVID-19 targets. Dynamics behavior of the best-docked molecules, compounds 15 and 6, within COVID-19 target pockets showed substantial stability of ligand-protein complexes as presented via several MD simulation parameters. Furthermore, calculated free-binding energies from MD simulation illustrated significant ligand's binding affinity towards respective target pockets. All provided findings supported the utility of scalarane-based sesterterpenes, particularly compounds 15 and 6, as promising lead candidates guiding the development of effective therapeutics against SARS-CoV-2.

In a search for potential drug candidates for combating COVID-19: computational study revealed salvianolic acid B as a potential therapeutic targeting 3CLpro and spike proteins

The global prevalence of COVID-19 disease and the overwhelming increase in death toll urge scientists to discover new effective drugs. Although the drug discovery process is a challenging and time-consuming, fortunately, the plant kingdom was found to have many active therapeutics possessing broad-spectrum antiviral activity including those candidates active against severe acute respiratory syndrome coronaviruses (SARS-CoV). Herein, nine traditional Chinese medicinal plant constituents from different origins (Glycyrrhizin 1, Lycorine 2, Puerarin 3, Daidzein 4, Daidzin 5, Salvianolic acid B 6, Dihydrotanshinone I 7, Tanshinone I 8, Tanshinone IIa 9) previously reported to exhibit antiviral activity against SARS-CoV were virtually screened in silico (molecular docking) as potential inhibitors of SARS-CoV-2 target proteins. The tested medicinal plant compounds were in silico screened for their activity against two key SARS-CoV-2 target proteins; 3CLpro, and Spike binding-domain proteins. Among the tested medicinal plant compounds, Salvianolic acid B 6 (Sal-B) showed promising binding affinities against the two specified SARS-CoV-2 target proteins compared to the reference standards used. Hence molecular dynamics simulations followed by calculating the free-binding energy were carried out for Sal-B providing information on its affinity, stability, and thermodynamic behavior within the two SARS-CoV-2 target proteins as well as key ligand-protein binding aspects. Besides, the quantum mechanical calculations showed that Sal-B can adopt different conformations due to the existence of various rotatable bonds. Therefore, the enhanced antiviral activity of Sal-B among other studied compounds can be also attributed to the structural flexibility of Sal-B. Our study gives an explanation of the structure activity relationship required for targeting SARS-CoV-2 3CLpro and Spike proteins and also facilitates the future design and synthesis of new potential drugs exhibiting better affinity and specificity. Besides, an ADME study was carried out on screened compounds and reference controls revealing their pharmacokinetics properties.Communicated by Ramaswamy H. Sarma.

On the search for COVID-19 therapeutics: identification of potential SARS-CoV-2 main protease inhibitors by virtual screening, pharmacophore modeling and molecular dynamics

COVID-19 also known as SARS-CoV-2 outbreak in late 2019 and its worldwide pandemic spread has taken the world by surprise. The minute-to-minute increasing coronavirus cases (>85 M) and progressive deaths (≈1.8 M) calls for finding a cure to this devastating pandemic. While there have been many attempts to find biologically active molecules targeting SARS-CoV-2 for treatment of this viral infection, none has found a way to the clinic yet. In this study, a 3-feature structure-based pharmacophore model was designed for SARS-CoV-2 main protease (M^Pro) that plays a vital role in the viral cellular penetration. High throughput virtual screening of the lead-like ZINC library was then performed to find a potent inhibitor employing the predesigned pharmacophore. In-silico pharmacokinetics/toxicity prediction study was subsequently applied towards the best hits. Finally, a 50 ns molecular dynamics simulation was carried out for the best hit and compared to the co-crystallized ligand where the hit compound displayed high binding and comparable interactions. The results identified new hits for SARS-CoV-2 M^Pro inhibition showing good docking score, pharmacokinetics and toxicity profile, drug-likeness, high binding energy in addition to a promising synthetic accessibility. Identifying new small compounds as potential leads for inhibiting SARS-CoV-2 is a very important step towards designing a synthesizing of promising drug candidates.Communicated by Ramaswamy H. Sarma.

Exploring the solid state phase transition in dl-norvaline with terahertz spectroscopy

Experimental and theoretical demonstration of the power of terahertz spectroscopy to provide novel insights into solid-state phase-transformations in organic materials.

Polymorphic phase transitions: Macroscopic theory and molecular simulation

Transformations in the solid state are of considerable interest, both for fundamental reasons and because they underpin important technological applications. The interest spans a wide spectrum of disciplines and application domains. For pharmaceuticals, a common issue is unexpected polymorphic transformation of the drug or excipient during processing or on storage, which can result in product failure. A more ambitious goal is that of exploiting the advantages of metastable polymorphs (e.g. higher solubility and dissolution rate) while ensuring their stability with respect to solid state transformation. To address these issues and to advance technology, there is an urgent need for significant insights that can only come from a detailed molecular level understanding of the involved processes. Whilst experimental approaches at best yield time- and space-averaged structural information, molecular simulation offers unprecedented, time-resolved molecular-level resolution of the processes taking place. This review aims to provide a comprehensive and critical account of state-of-the-art methods for modelling polymorph stability and transitions between solid phases. This is flanked by revisiting the associated macroscopic theoretical framework for phase transitions, including their classification, proposed molecular mechanisms, and kinetics. The simulation methods are presented in tutorial form, focusing on their application to phase transition phenomena. We describe molecular simulation studies for crystal structure prediction and polymorph screening, phase coexistence and phase diagrams, simulations of crystal-crystal transitions of various types (displacive/martensitic, reconstructive and diffusive), effects of defects, and phase stability and transitions at the nanoscale. Our selection of literature is intended to illustrate significant insights, concepts and understanding, as well as the current scope of using molecular simulations for understanding polymorphic transitions in an accessible way, rather than claiming completeness. With exciting prospects in both simulation methods development and enhancements in computer hardware, we are on the verge of accessing an unprecedented capability for designing and developing dosage forms and drug delivery systems in silico, including tackling challenges in polymorph control on a rational basis.

Do solid-to-solid polymorphic transitions in DL-norleucine proceed through nucleation?

DL-Norleucine is a molecular crystal exhibiting two enantiotropic phase transitions. The high temperature α ↔ γ transition has been shown to proceed through nucleation and growth [Mnyukh et al., J. Phys. Chem. Solids, 1975, 36, 127]. We focus on the low temperature β ↔ α transition in a combined computational and experimental study. The temperature dependence of the structural and energetic properties of both polymorphic forms is nearly identical. Molecular dynamics simulations and nudged elastic band calculations of the transition process itself, suggest that the transition is governed by cooperative movements of bilayers over relatively large energy barriers.

Water mediated reactions: TiO2and ZnO nanoparticle catalyzed multi component domino reaction in the synthesis of tetrahydroacridinediones, acridindiones, xanthenones and xanthenes

Tetrahydroacridine-1,8-(2H,5H,9H,10H)-diones4from 1,3-cyclohexanedione and/or dimedone 1,2-chloro-3-formylquinoline2and anilines3in water at 90 °C were obtained by domino reaction approach.

Size-dependent phase stability of a molecular nanocrystal: a proxy for investigating the early stages of crystallization

We make the link between the size-dependent phase stability of a nanocrystal and the phase-transition behavior of emerging crystallites during the earliest stages of crystallization, by using the former as a proxy for the latter. We outline an extension of the classical nucleation theory to describe crystal nucleation and subsequent transformations of competing polymorphic phases that characterize Ostwald's rule of stages. The theoretical framework reveals that the relative stability of the competing phases is a function of cluster size, which in turn varies with time, and therefore explains the complex transformation behavior observed for some systems. We investigated the stability of a nanocrystal of dl-norleucine by means of molecular simulation as a proxy for post-nucleation phase-transformation behavior in emerging crystallites. The simulations reveal that, for nanocrystals, the surface energy of the transition state of a transformation can dominate the barrier to phase change, thus causing metastable phases to be stabilized, not because they are thermodynamically stable, but rather due to kinetic hindering. Therefore, in the context of the earliest stages of crystal growth, not only does phase stability vary as a function of cluster size, and hence time, but thermodynamically feasible transformations are also prone to kinetic hindering.

Prediction of the mechanical behaviour of crystalline solids

PURPOSE

To explore the use of crystal inter-planar d-spacings and slip-plane interaction energies for predicting and characterising mechanical properties of crystalline solids.

METHODS

Potential relationships were evaluated between mechanical properties and inter-planar d-spacing, inter-planar interaction energy, and dispersive surface energy as determined using inverse gas chromatography (IGC) for a set of pharmaceutical materials. Inter-planar interaction energies were determined by molecular modelling.

RESULTS

General trends were observed between mechanical properties and the largest inter-planar d-spacing, inter-planar interaction energies, and IGC dispersive surface energy. A number of materials showed significant deviations from general trends. Weak correlations and outliers were rationalised.

CONCLUSIONS

Results suggest that the highest d-spacing of a material could serve as a first-order indicator for ranking mechanical behaviour of pharmaceutical powders, but with some reservation. Inter-planar interaction energy normalised for surface area shows only a weak link with mechanical properties and does not appear to capture essential physics of deformation. A novel framework linking mechanical properties of crystals to the distinct quantities, slip-plane energy barrier and inter-planar interaction (detachment) energy is proposed.

Uncovering molecular processes in crystal nucleation and growth by using molecular simulation

Exploring nucleation processes by molecular simulation provides a mechanistic understanding at the atomic level and also enables kinetic and thermodynamic quantities to be estimated. However, whilst the potential for modeling crystal nucleation and growth processes is immense, there are specific technical challenges to modeling. In general, rare events, such as nucleation cannot be simulated using a direct "brute force" molecular dynamics approach. The limited time and length scales that are accessible by conventional molecular dynamics simulations have inspired a number of advances to tackle problems that were considered outside the scope of molecular simulation. While general insights and features could be explored from efficient generic models, new methods paved the way to realistic crystal nucleation scenarios. The association of single ions in solvent environments, the mechanisms of motif formation, ripening reactions, and the self-organization of nanocrystals can now be investigated at the molecular level. The analysis of interactions with growth-controlling additives gives a new understanding of functionalized nanocrystals and the precipitation of composite materials.

Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials

Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.