Application of per-Residue Energy Decomposition to Design Peptide Inhibitors of PSD95 GK Domain

Identification of CXCR4 inhibitors as a key therapeutic small molecule in renal fibrosis

The final stage of almost all chronic kidney diseases is renal fibrosis. Simple wounds or persistent inflammation can cause tissue inflammation, which, in the case of the kidney, results in scarring. Vascular sclerosis, tubulointerstitial fibrosis and glomerular fibrosis are all types of kidney fibrosis. Renal damage and fibrosis are caused by elevated expression of CXCR4. This study aimed to identify possible pharmacological agents which could bind to and inhibit isoform I of CXCR4 and determine their strength of interactions. The I-TASSER, Galaxyweb and Robetta were used to predict and refine the structure of the CXCR4 protein. ModBase was used to improve the loops, and then the quality was evaluated by using the ERRAT value (92.15) and Ramachandran plot. The improved 3D structure was subjected to small molecule database docking using Maestro (from Schrodinger) and the glide module. GROMACS was used to simulate molecules with the three top low glide scores and the best ADME properties. The best glide score was achieved by ligand ID 4990 (-11.5). Simulations, free energy landscape and residue decomposition analysis revealed that 4990 interacted more consistently with CXCR4 than the other two small molecules.Communicated by Ramaswamy H. Sarma.

PCSK9 inhibitors as safer therapeutics for atherosclerotic cardiovascular disease (ASCVD): Pharmacophore design and molecular dynamics analysis

Cardiovascular disorders are still challenging and are among the deadly diseases. As a major risk factor for atherosclerotic cardiovascular disease, dyslipidemia, and high low-density lipoprotein cholesterol in particular, can be prevented primary and secondary by lipid-lowering medications. Therefore, insights are still needed into designing new drugs with minimal side effects. Proprotein convertase subtilisin/kexin 9 (PCSK9) enzyme catalyses protein-protein interactions with low-density lipoprotein, making it a critical target for designing promising inhibitors compared to statins. Therefore, we screened for potential compounds using a redesigned PCSK9 conformational behaviour to search for a significantly extensive chemical library and investigated the inhibitory mechanisms of the final compounds using integrated computational methods, from ligand essential functional group screening to all-atoms MD simulations and MMGBSA-based binding free energy. The inhibitory mechanisms of the screened compounds compared with the standard inhibitor. K31 and K34 molecules showed stronger interactions for PCSK9, having binding energy (kcal/mol) of -33.39 and -63.51, respectively, against -27.97 of control. The final molecules showed suitable drug-likeness, non-mutagenesis, permeability, and high solubility values. The C-α atoms root mean square deviation and root mean square fluctuation of the bound-PCSK9 complexes showed stable and lower fluctuations compared to apo PCSK9. The findings present a model that unravels the mechanism by which the final molecules proposedly inhibit the PCSK9 function and could further improve the design of novel drugs against cardiovascular diseases.

Revealing the Interaction Mechanism between Mycobacterium tuberculosis GyrB and Novobiocin, SPR719 through Binding Thermodynamics and Dissociation Kinetics Analysis

With the rapid emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), various levels of resistance against existing anti-tuberculosis (TB) drugs have developed. Consequently, the identification of new anti-TB targets and drugs is critically urgent. DNA gyrase subunit B (GyrB) has been identified as a potential anti-TB target, with novobiocin and SPR719 proposed as inhibitors targeting GyrB. Therefore, elucidating the molecular interactions between GyrB and its inhibitors is crucial for the discovery and design of efficient GyrB inhibitors for combating multidrug-resistant TB. In this study, we revealed the detailed binding mechanisms and dissociation processes of the representative inhibitors, novobiocin and SPR719, with GyrB using classical molecular dynamics (MD) simulations, tau-random acceleration molecular dynamics (τ-RAMD) simulations, and steered molecular dynamics (SMD) simulations. Our simulation results demonstrate that both electrostatic and van der Waals interactions contribute favorably to the inhibitors' binding to GyrB, with Asn52, Asp79, Arg82, Lys108, Tyr114, and Arg141 being key residues for the inhibitors' attachment to GyrB. The τ-RAMD simulations indicate that the inhibitors primarily dissociate from the ATP channel. The SMD simulation results reveal that both inhibitors follow a similar dissociation mechanism, requiring the overcoming of hydrophobic interactions and hydrogen bonding interactions formed with the ATP active site. The binding and dissociation mechanisms of GyrB with inhibitors novobiocin and SPR719 obtained in our work will provide new insights for the development of promising GyrB inhibitors.

Investigation of the Binding Interaction of Mfsd2a with NEDD4-2 via Molecular Dynamics Simulations

Major facilitator superfamily domain-containing 2a (Mfsd2a) is a sodium-dependent lysophosphatidylcholine cotransporter that plays an important role in maintaining the integrity of the blood-brain barrier and neurological function. Abnormal degradation of Mfsd2a often leads to dysfunction of the blood-brain barrier, while upregulation of Mfsd2a can retrieve neurological damage. It has been reported that Mfsd2a can be specifically recognized and ubiquitinated by neural precursor cell-expressed developmentally downregulated gene 4 type 2 (NEDD4-2) ubiquitin ligase and finally degraded through the proteasome pathway. However, the structural basis for the specific binding of Mfsd2a to NEDD4-2 is unclear. In this work, we combined deep learning and molecular dynamics simulations to obtain a Mfsd2a structure with high quality and a stable Mfsd2a/NEDD4-2-WW3 interaction model. Moreover, molecular mechanics generalized Born surface area (MM-GBSA) methods coupled with per-residue energy decomposition studies were carried out to analyze the key residues that dominate the binding interaction. Based on these results, we designed three peptides containing the key residues by truncating the Mfsd2a sequences. One of them was found to significantly inhibit Mfsd2a ubiquitination, which was further validated in an oxygen-glucose deprivation (OGD) model in a human microvascular endothelial cell line. This work provides some new insights into the understanding of Mfsd2a and NEDD4-2 interaction and might promote further development of drugs targeting Mfsd2a ubiquitination.

Design of Inhibitors That Target the Menin-Mixed-Lineage Leukemia Interaction

The prognosis of mixed-lineage leukemia (MLL) has remained a significant health concern, especially for infants. The minimal treatments available for this aggressive type of leukemia has been an ongoing problem. Chromosomal translocations of the KMT2A gene are known as MLL, which expresses MLL fusion proteins. A protein called menin is an important oncogenic cofactor for these MLL fusion proteins, thus providing a new avenue for treatments against this subset of acute leukemias. In this study, we report results using the structure-based drug design (SBDD) approach to discover potential novel MLL-mediated leukemia inhibitors from natural products against menin. The three-dimensional (3D) protein model was derived from Protein Databank (Protein ID: 4GQ4), and EasyModeller 4.0 and I-TASSER were used to fix missing residues during rebuilding. Out of the ten protein models generated (five from EasyModeller and I-TASSER each), one model was selected. The selected model demonstrated the most reasonable quality and had 75.5% of residues in the most favored regions, 18.3% of residues in additionally allowed regions, 3.3% of residues in generously allowed regions, and 2.9% of residues in disallowed regions. A ligand library containing 25,131 ligands from a Chinese database was virtually screened using AutoDock Vina, in addition to three known menin inhibitors. The top 10 compounds including ZINC000103526876, ZINC000095913861, ZINC000095912705, ZINC000085530497, ZINC000095912718, ZINC000070451048, ZINC000085530488, ZINC000095912706, ZINC000103580868, and ZINC000103584057 had binding energies of -11.0, -10.7, -10.6, -10.2, -10.2, -9.9, -9.9, -9.9, -9.9, and -9.9 kcal/mol, respectively. To confirm the stability of the menin-ligand complexes and the binding mechanisms, molecular dynamics simulations including molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) computations were performed. The amino acid residues that were found to be potentially crucial in ligand binding included Phe243, Met283, Cys246, Tyr281, Ala247, Ser160, Asn287, Asp185, Ser183, Tyr328, Asn249, His186, Leu182, Ile248, and Pro250. MI-2-2 and PubChem CIDs 71777742 and 36294 were shown to possess anti-menin properties; thus, this justifies a need to experimentally determine the activity of the identified compounds. The compounds identified herein were found to have good pharmacological profiles and had negligible toxicity. Additionally, these compounds were predicted as antileukemic, antineoplastic, chemopreventive, and apoptotic agents. The 10 natural compounds can be further explored as potential novel agents for the effective treatment of MLL-mediated leukemia.

Exploring the potential of biologically active phenolic acids from marine natural products as anticancer agents targeting the epidermal growth factor receptor

The epidermal growth factor receptor (EGFR) dimerizes upon ligand bindings to the extracellular domain that initiates the downstream signaling cascades and activates intracellular kinase domain. Thus, activation of autophosphorylation through kinase domain results in metastasis, cell proliferation, and angiogenesis. The main objective of this research is to discover more promising anti-cancer lead compound against EGRF from the phenolic acids of marine natural products using in-silico approaches. Phenolic compounds reported from marine sources are reviewed from previous literatures. Furthermore, molecular docking was carried out using the online tool CB-Dock. The molecules with good docking and binding energies scores were subjected to ADME, toxicity and drug-likeness analysis. Subsequently, molecules from the docking experiments were also evaluated using the acute toxicity and MD simulation studies. Fourteen phenolic compounds from the reported literatures were reviewed based on the findings, isolation, characterized and applications. Molecular docking studies proved that the phenolic acids have good binding fitting by forming hydrogen bonds with amino acid residues at the binding site of EGFR. Chlorogenic acid, Chicoric acid and Rosmarinic acid showed the best binding energies score and forming hydrogen bonds with amino acid residues compare to the reference drug Erlotinib. Among these compounds, Rosmarinic acid showed the good pharmacokinetics profiles as well as acute toxicity profile. The MD simulation study further revealed that the lead complex is stable and could be future drug to treat the cancer disease. Furthermore, in a wet lab environment, both in-vitro and in-vivo testing will be employed to validate the existing computational results.Communicated by Ramaswamy H. Sarma.

In silico approaches for the identification of novel ULK1 inhibitors: pharmacophore model, molecular docking and molecular dynamics simulations

The serine/threonine kinase unc-51-like autophagy activating kinase 1 (ULK1) has been regarded as an attractive target for tumor therapy. In this study, in silico approaches, such as the pharmacophore-based virtual screening strategy, molecular docking and molecular dynamics (MD) simulations, were applied to develop novel potential ULK1 inhibitors. The pharmacophore models based on known aminopyrimidine ULK1 inhibitors were constructed to screen the dataset of 1.68 million compounds, which were obtained via screening the 2.30 million compounds in ChEMBL database by Lipinski's rule of five. Seven novel compounds and 1 known ULK1 inhibitor stand out for the strong virtual biological activity by molecular docking, cluster analysis, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation and Absorption Distribution Metabolism Excretion Toxicity (ADMET) prediction. Their results of MD included principal component analysis (PCA) and Free Energy Landscapes surface (FELs) indicated that the protein-ligand complex was stable in simulated trajectories of 100 ns. The binding free energy (BFE) calculations showed that a total of 6 novel compounds (CL130, CL834, CL961, CL966, CL163 and CL329) with the stable binding state and stronger BFE (-61.17 to -37.01 kcal/mol) than that of original ligand 3RF (-36.66 kcal/mol). With reference to the ULK1 inhibition of 3RF (IC50 = 160 nM), it can be inferred that these compounds could be used as a new type of potential ULK1 inhibitors and be worthy of further investigation for tumor treatments.Communicated by Ramaswamy H. Sarma.

Docking and Molecular Dynamics Simulations Clarify Binding Sites for Interactions of Novel Marine Sulfated Glycans with SARS-CoV-2 Spike Glycoprotein

The entry of SARS-CoV-2 into the host cell is mediated by its S-glycoprotein (SGP). Sulfated glycans bind to the SGP receptor-binding domain (RBD), which forms a ternary complex with its receptor angiotensin converting enzyme 2. Here, we have conducted a thorough and systematic computational study of the binding of four oligosaccharide building blocks from novel marine sulfated glycans (isolated from Pentacta pygmaea and Isostichopus badionotus) to the non-glycosylated and glycosylated RBD. Blind docking studies using three docking programs identified five potential cryptic binding sites. Extensive site-targeted docking and molecular dynamics simulations using two force fields confirmed only two binding sites (Sites 1 and 5) for these novel, highly charged sulfated glycans, which were also confirmed by previously published reports. This work showed the structural features and key interactions driving ligand binding. A previous study predicted Site 2 to be a potential binding site, which was not observed here. The use of several molecular modeling approaches gave a comprehensive assessment. The detailed comparative study utilizing multiple modeling approaches is the first of its kind for novel glycan-SGP interaction characterization. This study provided insights into the key structural features of these novel glycans as they are considered for development as potential therapeutics.

Design of Tracers in Fluorescence Polarization Assay for Extensive Application in Small Molecule Drug Discovery

Development of fluorescence polarization (FP) assays, especially in a competitive manner, is a potent and mature tool for measuring the binding affinities of small molecules. This approach is suitable for high-throughput screening (HTS) for initial ligands and is also applicable for further study of the structure-activity relationships (SARs) of candidate compounds for drug discovery. Buffer and tracer, especially rational design of the tracer, play a vital role in an FP assay system. In this perspective, we provided different kinds of approaches for tracer design based on successful cases in recent years. We classified these tracers by different types of ligands in tracers, including peptide, nucleic acid, natural product, and small molecule. To make this technology accessible for more targets, we briefly described the basic theory and workflow, followed by highlighting the design and application of typical FP tracers from a perspective of medicinal chemistry.

Simplified Derivatives of Tetrandrine as Potent and Specific P-gp Inhibitors to Reverse Multidrug Resistance in Cancer Chemotherapy

Targeted inhibition of a drug efflux transporter P-glycoprotein (P-gp) is an important strategy to reverse multidrug resistance in cancer chemotherapy. In this study, a rationally structural simplification to natural tetrandrine was performed based on molecular dynamics simulation and fragment growth, leading to an easily prepared, novel, and simplified compound OY-101 with high reversal activity and low cytotoxicity. Its excellent synergistic anti-cancer effect with vincristine (VCR) against drug-resistant cells Eca109/VCR was confirmed by reversal activity assay, flow cytometry, plate clone formation assay, and drug synergism analysis (IC50 = 9.9 nM, RF = 690). Further mechanism study confirmed that the OY-101 was a specific and efficient P-gp inhibitor. Importantly, OY-101 increased VCR sensitization in vivo without obvious toxicity. Overall, our findings may provide an alternative strategy for the design of novel specific P-gp inhibitor as an anti-tumor chemotherapy sensitizer.