New insights into the binding mode of melanin concentrating hormone receptor-1 antagonists: homology modeling and explicit membrane molecular dynamics simulation study

Identification of novel inhibitors of Neisseria gonorrhoeae MurI using homology modeling, structure-based pharmacophore, molecular docking, and molecular dynamics simulation-based approach

MurI is one of the most significant role players in the biosynthesis of the peptidoglycan layer in Neisseria gonorrhoeae (Ng). We attempted to highlight the structural and functional relationship between Ng-MurI and D-glutamate to design novel molecules targeting this interaction. The three-dimensional (3D) model of the protein was constructed by homology modeling and the quality and consistency of generated model were assessed. The binding site of the protein was identified by molecular docking studies and a pharmacophore was identified using the interactions of the control ligand. The structure-based pharmacophore model was validated and employed for high-throughput virtual screening and molecular docking to identify novel Ng-MurI inhibitors. Finally, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with the substrate glutamate and novel molecules facilitated us to confirm the stability of the protein-ligand docked complexes. The 100 ns MD simulations of the potential lead compounds with protein confirmed that the modeled complexes were stable. This study identifies novel potential compounds with good fitness and docking scores, which made the interactions of biological significance within the protein active site. Hence, the identified compounds may act as new leads to design and develop Ng-MurI inhibitors.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.

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.

Identification of a new small molecule chemotype of Melanin Concentrating Hormone Receptor-1 antagonists using pharmacophore-based virtual screening

MCH receptor is a G protein-coupled receptor with two subtypes R1 and R2. Many studies have demonstrated the role of MCH-R1 in feeding and energy homeostasis. It has been proven that oral administration of small molecule MCH-R1 antagonists significantly reduces food intake and causes a dose-dependent weight loss. In this study, two ligand-based pharmacophores were developed and validated based on recently published MCH-R1 antagonists with diverse structures. Successful pharmacophores had one hydrogen bond acceptor, one positive ionizable, one ring aromatic and two or three hydrophobic groups. These 3D-QSAR models were used for virtual screening of the ZINC chemical database resulting in the identification of nine compounds with more than 50% displacement of radiolabeled MCH at a 20 μM concentration. Moreover, four of these compounds showed antagonistic activities in Aequorin functional assay, including MH-3 which is the first MCH-R1 antagonist based on a diazaspiro[4.5]decane scaffold. The most active compounds were also docked into our previously published MCH-R1 homology model to gain insights into their binding determinants. These compounds could represent a viable starting scaffold for the design of potent MCH-R1 antagonists with improved pharmacokinetic properties as an effective treatment for obesity.

Profiling the interaction mechanism of quinoline/quinazoline derivatives as MCHR1 antagonists: an in silico method

Melanin concentrating hormone receptor 1 (MCHR1), a crucial regulator of energy homeostasis involved in the control of feeding and energy metabolism, is a promising target for treatment of obesity. In the present work, the up-to-date largest set of 181 quinoline/quinazoline derivatives as MCHR1 antagonists was subjected to both ligand- and receptor-based three-dimensional quantitative structure–activity (3D-QSAR) analysis applying comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The optimal predictable CoMSIA model exhibited significant validity with the cross-validated correlation coefficient (Q²) = 0.509, non-cross-validated correlation coefficient (R²_ncv) = 0.841 and the predicted correlation coefficient (R²_pred) = 0.745. In addition, docking studies and molecular dynamics (MD) simulations were carried out for further elucidation of the binding modes of MCHR1 antagonists. MD simulations in both water and lipid bilayer systems were performed. We hope that the obtained models and information may help to provide an insight into the interaction mechanism of MCHR1 antagonists and facilitate the design and optimization of novel antagonists as anti-obesity agents.

Structurally conserved binding sites of hemagglutinin as targets for influenza drug and vaccine development

ProBiS is a new method to identify the binding site of protein through local structural alignment against the nonredundant Protein Data Bank (PDB), which may result in unique findings compared to the energy-based, geometry-based, and sequence-based predictors. In this work, binding sites of Hemagglutinin (HA), which is an important target for drugs and vaccines in influenza treatment, have been revisited by ProBiS. For the first time, the identification of conserved binding sites by local structural alignment across all subtypes and strains of HA available in PDB is presented. ProBiS finds three distinctive conserved sites on HA's structure (named Site 1, Site 2, and Site 3). Compared to other predictors, ProBiS is the only one that accurately defines the receptor binding site (Site 1). Apart from that, Site 2, which is located slightly above the TBHQ binding site, is proposed as a potential novel conserved target for membrane fusion inhibitor. Lastly, Site 3, located around Helix A at the stem domain and recently targeted by cross-reactive antibodies, is predicted to be conserved in the latest H7N9 China 2013 strain as well. The further exploration of these three sites provides valuable insight in optimizing the influenza drug and vaccine development.

3D QSAR pharmacophore modeling, in silico screening, and density functional theory (DFT) approaches for identification of human chymase inhibitors

Human chymase is a very important target for the treatment of cardiovascular diseases. Using a series of theoretical methods like pharmacophore modeling, database screening, molecular docking and Density Functional Theory (DFT) calculations, an investigation for identification of novel chymase inhibitors, and to specify the key factors crucial for the binding and interaction between chymase and inhibitors is performed. A highly correlating (r = 0.942) pharmacophore model (Hypo1) with two hydrogen bond acceptors, and three hydrophobic aromatic features is generated. After successfully validating "Hypo1", it is further applied in database screening. Hit compounds are subjected to various drug-like filtrations and molecular docking studies. Finally, three structurally diverse compounds with high GOLD fitness scores and interactions with key active site amino acids are identified as potent chymase hits. Moreover, DFT study is performed which confirms very clear trends between electronic properties and inhibitory activity (IC(50)) data thus successfully validating "Hypo1" by DFT method. Therefore, this research exertion can be helpful in the development of new potent hits for chymase. In addition, the combinational use of docking, orbital energies and molecular electrostatic potential analysis is also demonstrated as a good endeavor to gain an insight into the interaction between chymase and inhibitors.