Identification of potential RapJ hits as sporulation pathway inducer candidates in Bacillus coagulans via structure-based virtual screening and molecular dynamics simulation studies.

Phytocompounds of Senecio candicans as potential acetylcholinesterase inhibitors targeting Alzheimer's disease: A structure-based virtual screening and molecular dynamics simulation study.

Identification of novel serotonin reuptake inhibitors targeting central and allosteric binding sites: A virtual screening and molecular dynamics simulations study

Correction to: Identification of potential RapJ hits as sporulation pathway inducer candidates in Bacillus coagulans via structure-based virtual screening and molecular dynamics simulation studies.

Alzheimer is a severe memory and cognitive impairment neurodegenerative disease that is the most common cause of dementia worldwide and characterized by the pathological accumulation of tau protein and amyloid-beta peptides. In this study, we have developed E-pharmacophore modeling to screen the eMolecules database with the help of a reported co-crystal structure bound with Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE-1). Flumemetamol, florbetaben, and florbetapir are currently approved drugs for use in the clinical diagnosis of Alzheimer’s disease. Despite the benefits of commercially approved drugs, there is still a need for novel diagnostic agents with enhanced physicochemical and pharmacokinetic properties compared to those currently used in clinical practice and research. In the E-pharmacophore modeling results, it is revealed that two aromatic rings (R19, R20), one donor (D12), and one acceptor (A8) are obtained, and also that similar pharmacophoric features of compounds are identified from pharmacophore-based virtual screening. The identified screened hits were filtered for further analyses using structure-based virtual screening and MM/GBSA. From the analyses, top hits such as ZINC39592220 and en1003sfl.46293 are selected based on their top docking scores (−8.182 and −7.184 Kcal/mol, respectively) and binding free energy (−58.803 and −56.951 Kcal/mol, respectively). Furthermore, a molecular dynamics simulation and MMPBSA study were performed, which revealed admirable stability and good binding free energy throughout the simulation period. Moreover, Qikprop results revealed that the selected, screened hits have good drug-likeness and pharmacokinetic properties. The screened hits ZINC39592220 and en1003sfl.46293 could be used to develop drug molecules against Alzheimer’s disease.

A mechanistic approach to explore novel HDAC1 inhibitor using pharmacophore modeling, 3D- QSAR analysis, molecular docking, density functional and molecular dynamics simulation study

Dengue fever has been a global health concern. Mitigation is a challenging problem due to non-availability of workable treatments. The most difficult objective is to design a perfect anti-dengue agent capable of inhibiting infections caused by all four serotypes. Various tactics have been employed in the past to discover dengue antivirals, including screening of chemical compounds against dengue virus enzymes. The objective of the current study is to investigate phytocompounds as anti-dengue remedies that target the non-structural 2B and non-structural 3 protease (NS2B-NS3pro), a possible therapeutic target for dengue fever. Initially, 300 + antiviral phytocompounds were collected from Duke’s phytochemical and ethnobotanical database and 30 phytocompounds with anti-dengue properties were identified from previously reported studies, which were virtually screened against NS2B-NS3pro using molecular docking and toxicity evaluation. The top five most screened ligands were naringin, hesperidin, gossypol, maslinic acid and rhodiolin with binding affinities of − 8.7 kcal/mol, − 8.5 kcal/mol, − 8.5 kcal/mol, − 8.5 kcal/mol and − 8.1 kcal/mol, respectively. The finest docked compounds complexed with NS2B-NS3pro were subjected for molecular dynamics (MD) simulations and binding free energy estimations through molecular mechanics generalized born surface area–based calculations. The results of the study are intriguing in the context of computer-aided screening and the binding affinities of the phytocompounds, proposing maslinic acid (MAS) as a potent bioactive antiviral for the development of phytocompound-based anti-dengue agent.Graphical abstract Supplementary InformationThe online version contains supplementary material available at 10.1007/s00894-022-05355-w.

DNA gyrase brings negative supercoils into DNA and loosens up certain positive supercoils that collect during replication and transcription and is a notable antibacterial target. To fight against the menace of antibiotic-resistant bacterial infections, we have employed various computational tools like high throughput virtual screening (HTVS), standard precision (SP) docking, extra precision (XP) docking, and molecular dynamics (MD) simulation studies to identify some potential DNA gyrase inhibitors. A focused library of 5968 anti-bacterial compounds was screened using the HTVS docking protocol of the glide module of Maestro. The top 200 docked compounds were further filtered using SP and XP docking protocols, and their free binding energies were calculated using MM-GBSA studies. The binding and stability of the top two compounds which showed better docking scores than the co-crystallized ligand (Clorobiocin) of DNA gyrase (PDB ID: 1KZN) were further probed by MD simulation of 100 ns using GROMACS. MD simulation study suggested that the compounds AM1 and AM5 form a stable complex with DNA gyrase with a good number of hydrogen bonds. XP docking study showed that interaction with the crucial amino acids for compounds AM1 and AM5 was like the co-crystallized ligand. These compounds were also predicted to be drug-like molecules with good water solubility and excellent absorption profiles. Based on the above studies, herein we report compounds AM1 (1R,3S)-1-(2-((3-(ammoniomethyl)phenyl)amino)-2-oxoethyl)-3-carbamoylpiperidin-1-ium and AM5 (1'S,2 s,4R)-4-ammonio-6-ethyl-1'-methylspiro[chromane-2,4'-piperidin]-1'-ium as potential DNA gyrase inhibitors which can be further developed as a potential lead against the menace of antibiotic resistance.

Antiviral drug therapy against SARS-CoV-2 is not yet established and posing a serious global health issue. Remdesivir is the first antiviral compound approved by the US FDA for the SARS-CoV-2 treatment for emergency use, targeting RNA-dependent RNA polymerase (RdRp) enzyme. In this work, we have examined the action of remdesivir and other two ligands screened from the library of nucleotide analogues using docking and molecular dynamics (MD) simulation studies. The MD simulations have been performed for all the ligand-bound RdRp complexes for the 30 ns time scale. This is one of the earlier reports to perform the MD simulations studies using the SARS-CoV-2 RdRp crystal structure (PDB ID 7BTF). The MD trajectories were analyzed and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were performed to calculate the binding free energy. The binding energy data reveal that compound-17 (-59.6 kcal/mol) binds more strongly as compared to compound-8 (-46.3 kcal/mol) and remdesivir (-29.7 kcal/mol) with RdRp. The detailed analysis of trajectories shows that the remdesivir binds in the catalytic site and forms a hydrogen bond with the catalytic residues from 0 to 0.46 ns. Compound-8 binds in the catalytic site but does not form direct hydrogen bonds with catalytic residues. Compound-17 showed the formation of hydrogen bonds with catalytic residues throughout the simulation process. The MD simulation results such as hydrogen bonding, the center of mass distance analysis, snapshots at a different time interval, and binding energy suggest that compound-17 binds strongly with RdRp of SARS-CoV-2 and has the potential to develop as a new antiviral against COVID-19. Further, the frontier molecular orbital analysis and molecular electrostatic potential (MESP) iso-surface analysis using DFT calculations shed light on the superior binding of compound-17 with RdRp compared to remdesivir and compound-8. The computed as well as the experimentally reported pharmacokinetics and toxicity parameters of compound-17 is encouraging and therefore can be one of the potential candidates for the treatment of COVID-19.

Background: Multiple myeloma (MM) is a hematological malignancy of plasma cells that produce a monoclonal immunoglobulin protein. Despite significant advances in the treatment of MM, currently available therapies are associated with toxicity and resistance. As a result, there is an increasing demand for novel, effective therapeutics. Inhibition of histone deacetylases (HDACs) is emerging as a potential method for treating cancer. HDAC6 is one of 18 different HDAC isoforms that regulate tubulin lysine 40 and function in the microtubule network. HDAC6 participates in tumorigenesis and metastasis through protein ubiquitination, tubulin, and Hsp90. Several studies have found that inhibiting HDAC6 causes AKT and ERK dephosphorylation, which leads to decreased cell proliferation and promotes cancer cell death via the PI3K/AKT and MAPK/ERK signaling pathways. Objective: The objective of this study is to target HDAC6 and identify potent inhibitors for the treatment of multiple myeloma by employing computer-aided drug design. Materials and Methods: A total of 199,611,439 molecules from five different chemical databases, such as CHEMBL25, ChemSpace, Mcule, MolPort, and ZINC, have been screened against HDAC6 by structure- based virtual screening, followed by filtering for various drug-likeness, ADME, toxicity, consensus molecular docking, and 100 ns MD simulation. Results: Our research work resulted in three molecules that have shown strong binding affinity (CHEMBL2425964 -9.99 kcal/mol, CHEMBL2425966 -9.89 kcal/mol, and CSC067477144 -9.86 kcal/mol) at the active site HDAC6, along with effective ADME properties, low toxicity, and high stability. Inhibiting HDAC6 with these identified molecules will induce AKT and ERK dephosphorylation linked to reduced cell proliferation and promote cancer cell death. Conclusion: CHEMBL2425964, CHEMBL2425966, and CSC067477144 could be effective against multiple myeloma.

Both receptor-binding domain in spike protein (S-RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and human neuropilin-1 (NRP1) are important in the virus entry, and their concomitant inhibition may become a potential strategy against the SARS-CoV-2 infection. Herein, five novel dual S-RBD/NRP1-targeting peptides with nanomolar binding affinities were identified by structure-based virtual screening. Particularly, RN-4 was found to be the most promising peptide targeting S-RBD (K d = 7.4 ± 0.5 nM) and NRP1-BD (the b1 domain of NRP1) (K d = 16.1 ± 1.1 nM) proteins. Further evidence in the pseudovirus infection assay showed that RN-4 can significantly inhibit the SARS-CoV-2 pseudovirus entry into 293 T cells (EC50 = 0.39 ± 0.09 μM) without detectable side effects. These results suggest that RN-4, a novel dual S-RBD/NRP1-targeting agent, holds potential as an effective therapeutic to combat the SARS-CoV-2 infection.

The use of US FDA-approved drugs is preferred due to the need for lower costs and less time. In in silico medicine, repurposing is a quick and accurate way to screen US FDA-approved medications to find a therapeutic option for COVID-19 infection. Dual inhibitors possess dual inhibitory activity, which may be due to the inhibition of two different enzymes, and are considered better than combination therapy from the developmental and clinical perspectives. In this study, a molecular docking simulation was performed to identify the interactions of antiviral drugs with the critical residues in the binding site of the main SARS-CoV-2 protease, spike glycoprotein, and papain-like protease receptors compared to the angiotensin-converting enzyme-related carboxypeptidase (ACE2) receptor of host cells. Each of the receptors was docked with 70 US FDA-approved antiviral drugs using AutoDock Vina. A molecular dynamics (MD) simulation study was also used for 100 ns to confirm the stability behaviour of the ligand receptor complexes. Among the drugs that had the strongest interaction with the SARS-CoV-2 main protease, spike glycoprotein and papain-like protease receptors, and host cell ACE2 receptors, Simeprevir, Maraviroc and Saquinavir had dual inhibitory effects. The MD simulation study confirmed the stability of the strongest interactions between the antiviral drugs and the main protease, ACE2, spike glycoprotein, and papain-like protease receptors to 100 ns. However the results of MMPBSA analysis showed that the bond between Saquinavir and the ACE2 receptor was weak. Simeprevir and Maraviroc drugs had acceptable binding energies with dual receptors, especially the Simeprevir. Communicated by Ramaswamy H. Sarma

Recent outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to a pandemic of COVID-19. The absence of a therapeutic drug and vaccine is causing severe loss of life and economy worldwide. SARS-CoV and SARS-CoV-2 employ the host cellular serine protease TMPRSS2 for spike (S) protein priming for viral entry into host cells. A potential way to reduce the initial site of SARS-CoV-2 infection may be to inhibit the activity of TMPRSS2. In the current study, the three-dimensional structure of TMPRSS2 was generated by homology modelling and subsequently validated with a number of parameters. The structure-based virtual screening of Selleckchem database was performed through ‘Virtual Work Flow’ (VSW) to find out potential lead-like TMPRSS2 inhibitors. Camostat and bromhexine are known TMPRSS2 inhibitor drugs, hence these were used as control molecules throughout the study. Based on better dock score, binding-free energy and binding interactions compared to the control molecules, six molecules (Neohesperidin, Myricitrin, Quercitrin, Naringin, Icariin, and Ambroxol) were found to be promising against the TMPRSS2. Binding interactions analysis revealed a number of significant binding interactions with binding site amino residues of TMPRSS2. The all-atoms molecular dynamics (MD) simulation study indicated that all proposed molecules retain inside the receptor in dynamic states. The binding energy calculated from the MD simulation trajectories also favour the strong affinity of the molecules towards the TMPRSS2. Proposed molecules belong to the bioflavonoid class of phytochemicals and are reported to possess antiviral activity, our study indicates their possible potential for application in COVID-19.

Cytochrome P450 3A5 (CYP3A5) is one of the crucial CYP family members and has already proven to be an important drug target for cardiovascular diseases. In the current study, the PubChem database was screened through molecular docking and high-affinity molecules were adopted for further assessment. A negative image-based (NIB) model was used for a similarity search by considering the complementary shape and electrostatics of the target and small molecules. Further, the molecules were segregated into active and inactive groups through six machine learning (ML) matrices. The active molecules found in each ML model were used for in silico pharmacokinetics and toxicity assessments. A total of five molecules followed the acceptable pharmacokinetics and toxicity profiles. Several potential binding interactions between the proposed molecules and CYP3A5 were observed. The dynamic behavior of the selected molecules in the CYP3A5 was explored through a molecular dynamics (MD) simulation study. Several parameters obtained from the MD simulation trajectory explained the stability of the protein–ligand complexes in dynamic states. The high binding affinity of each molecule was revealed by the binding free energy calculation through the MM-GBSA methods. Therefore, it can be concluded that the proposed molecules might be potential CYP3A5 molecules for therapeutic application in cardiovascular diseases subjected to in vitro/in vivo validations.

The global rise of antibiotic-resistant infections has been driven in part by the spread of bacteria producing metallo-β-lactamase (MBL), particularly New Delhi metallo-β-lactamase-1 (NDM-1). Currently, there are no clinically approved inhibitors targeting NDM-1 or other MBLs, highlighting the urgent need for novel therapeutic agents. This study addresses this gap by identifying potential NDM-1 inhibitors through a comprehensive in silico workflow. A total of 4,561 natural product compounds were screened using a machine learning (ML)-based quantitative structure-activity relationship (QSAR) model. Molecular docking was then performed to prioritize hits, followed by Tanimoto similarity-based clustering to identify representative compounds. The three most promising compounds identified were S721-1034, S904-0022, and N118-0137. 300 ns molecular dynamics (MD) simulation was used to examine binding interactions and stability of a control molecule (meropenem (0RV)) and the three selected compounds (S721-1034, S904-0022, and N118-0137) with the target protein. Among the three compounds evaluated, S904-0022 demonstrated consistent root mean square deviation (RMSD) values throughout the molecular dynamics (MD) simulation compared to the other two ligands. Additionally, S904-0022 exhibited considerable affinity with key residues, including Gln123, His250, Trp93, and Val73, indicating robust interactions with NDM-1. The strength of this interaction was further validated by a significantly favorable binding free energy of -35.77 kcal/mol, markedly better than the control compound (-18.90 kcal/mol). The strength of this interaction was further validated by a significantly favorable binding free energy of -35.77 kcal/mol, markedly better than the control compound (-18.90 kcal/mol). The findings of this study provide valuable insights into the molecular interactions and stability of these compounds, which can be used to improve drug development and explore the interactions between proteins and ligands. The study concluded that S904-0022 exhibited substantial therapeutic potential and requires additional experimental exploration as a potential NDM-1 inhibitor.

The cause of Alzheimer's disease is related to aggregates such as oligomers and amyloid fibrils consisting of amyloid-β (Aβ) peptides. Molecular dynamics (MD) simulation studies have been conducted to understand the molecular mechanism of the formation and disruption of Aβ aggregates. In this Perspective, the MD simulation studies are classified into four categories, focusing on the target systems: aggregation of Aβ peptides in bulk solution, Aβ aggregation at the interface, aggregation inhibitor against Aβ peptides, and nonequilibrium MD simulation of Aβ aggregates. MD simulation studies in these categories are first reviewed. Future perspectives in each category are then presented. Finally, the overall perspective is presented on how MD simulations of Aβ aggregates can be utilized for developing Alzheimer's disease treatment.