MLP Tools: a PyMOL plugin for using the molecular lipophilicity potential in computer-aided drug design

CageCavityCalc (C3): A Computational Tool for Calculating and Visualizing Cavities in Molecular Cages

Organic(porous) and metal-organic cages are promising biomimetic platforms with diverse applications spanning recognition, sensing, and catalysis. The key to the emergence of these functions is the presence of well-defined inner cavities capable of binding a wide range of guest molecules and modulating their properties. However, despite the myriad cage architectures currently available, the rational design of structurally diverse and functional cages with specific host-guest properties remains challenging. Efficiently predicting such properties is critical for accelerating the discovery of novel functional cages. Herein, we introduce CageCavityCalc (C3), a Python-based tool for calculating the cavity size of molecular cages. The code is available on GitHub at https://github.com/VicenteMartiCentelles/CageCavityCalc. C3 utilizes a novel algorithm that enables the rapid calculation of cavity sizes for a wide range of molecular structures and porous systems. Moreover, C3 facilitates easy visualization of the computed cavity size alongside hydrophobic and electrostatic potentials, providing insights into host-guest interactions within the cage. Furthermore, the calculated cavity can be visualized using widely available visualization software, such as PyMol, VMD, or ChimeraX. To enhance user accessibility, a PyMol plugin has been created, allowing nonspecialists to use this tool without requiring computer programming expertise. We anticipate that the deployment of this computational tool will significantly streamline cage cavity calculations, thereby accelerating the discovery of functional cages.

Structure and Ligand Based Virtual Screening and MPO Topological Analysis of Triazolo Thiadiazepine-fused Coumarin Derivatives as Anti-Parkinson Drug Candidates

Parkinson's disease (PD) is a debilitating condition that can cause locomotor problems in affected patients, such as tremors and body rigidity. PD therapy often includes the use of monoamine oxidase B (MAOB) inhibitors, particularly phenylhalogen compounds and coumarin-based semi-synthetic compounds. The objective of this study was to analyze the structural, pharmacokinetic, and pharmacodynamic profile of a series of Triazolo Thiadiazepine-fused Coumarin Derivatives (TDCDs) against MAOB, in comparison with the inhibitor safinamide. To achieve this goal, we utilized structure-based virtual screening techniques, including target prediction and absorption, distribution, metabolism, and excretion (ADME) prediction based on multi-parameter optimization (MPO) topological analysis, as well as ligand-based virtual screening techniques, such as docking and molecular dynamics. The findings indicate that the TDCDs exhibit structural similarity to other bioactive compounds containing coumarin and MAOB-binding azoles, which are present in the ChEMBL database. The topological analyses suggest that TDCD3 has the best ADME profile, particularly due to the alignment between low lipophilicity and high polarity. The coumarin and triazole portions make a strong contribution to this profile, resulting in a permeability with Papp estimated at 2.15 × 10-5 cm/s, indicating high cell viability. The substance is predicted to be metabolically stable. It is important to note that this is an objective evaluation based on the available data. Molecular docking simulations showed that the ligand has an affinity energy of - 8.075 kcal/mol with MAOB and interacts with biological substrate residues such as Pro102 and Phe103. The results suggest that the compound has a safe profile in relation to the MAOB model, making it a promising active ingredient for the treatment of PD.

Melatonin Inhibits AGS Cell Proliferation by Binding to the ATP Binding Site of CDK2 Under Hyperglycemic Conditions

Cancer cells utilize glucose as their primary energy source. The aggressive nature of cancer cells is therefore enhanced in hyperglycemic conditions. This study has been adopted to investigate the therapeutic potential of melatonin against such aggressive proliferation of AGS cells-a human gastric cancer cell line, under hyperglycemic conditions. AGS cells were incubated with high glucose-containing media, and the effects of melatonin have been evaluated, therein. Cell proliferation, ROS generation, flow-cytometric analysis for cell cycle and apoptosis, wound healing, immunoblotting, zymography, reverse zymography assays, in-silico analysis, and kinase activity assays were performed to evaluate the effects of melatonin. We observed that melatonin inhibited the hyperglycemia-induced cell proliferation in a dose-dependent manner. It further altered the expression and activity of MMP-9 and TIMP-1. Moreover, melatonin inhibited AGS cell proliferation by arresting AGS cells in the G0/G1 phase after binding in the ATP binding site of CDK-2, thereby inhibiting its kinase activity. In association, a significant decrease in the expression of cyclin D1, cyclin E, CDK-4, and CDK-2 were observed. In conclusion, these findings suggest that melatonin has anti-gastric cancer potential. Melatonin could therefore be included in future drug designs for gastric cancer-hyperglycemia co-morbidity treatment.

Ramaswamy V, Vargiu AV, Malloci G, Bosin A, Ruggerone P. Common recognition topology of mex transporters of Pseudomonas aeruginosa revealed by molecular modelling

The secondary transporters of the resistance-nodulation-cell division (RND) superfamily mediate multidrug resistance in Gram-negative bacteria like Pseudomonas aeruginosa. Among these RND transporters, MexB, MexF, and MexY, with partly overlapping specificities, have been implicated in pathogenicity. Only the structure of the former has been resolved experimentally, which together with the lack of data about the functional dynamics of the full set of transporters, limited a systematic investigation of the molecular determinants defining their peculiar and shared features. In a previous work (Ramaswamy et al., Front. Microbiol., 2018, 9, 1144), we compared at an atomistic level the two main putative recognition sites (named access and deep binding pockets) of MexB and MexY. In this work, we expand the comparison by performing extended molecular dynamics (MD) simulations of these transporters and the pathologically relevant transporter MexF. We employed a more realistic model of the inner phospholipid membrane of P. aeruginosa and more accurate force-fields. To elucidate structure/dynamics-activity relationships we performed physico-chemical analyses and mapped the binding propensities of several organic probes on all transporters. Our data revealed the presence, also in MexF, of a few multifunctional sites at locations equivalent to the access and deep binding pockets detected in MexB. Furthermore, we report for the first time about the multidrug binding abilities of two out of five gates of the channels deputed to peripheral (early) recognition of substrates. Overall, our findings help to define a common “recognition topology” characterizing Mex transporters, which can be exploited to optimize transport and inhibition propensities of antimicrobial compounds.

In silico proteolysis and molecular interaction of tilapia (Oreochromis niloticus) skin collagen-derived peptides for environmental remediation

Fish skin collagen hydrolyzate has demonstrated the potent inhibition of dipeptidyl peptidase-IV (DPP-IV), one of the treatments for type-2 diabetes mellitus (type-2 DM), but the precise mechanism is still unclear. This study used in silico method to evaluate the potential of the active peptides from tilapia skin collagen (Oreochromis niloticus) for DPP-IV inhibitor. The methodology includes collagen hydrolysis using BIOPEP, which is the database of bioactive peptides; active peptide selection; toxicity, allergenicity, sensory analysis of active peptides; and binding of active peptides to DPP-IV compared with linagliptin. The result indicated that in silico enzymatic hydrolysis of collagen produced active peptides with better prediction of biological activity than intact collagen. There are 13 active peptides were predicted as non-toxic and non-allergenic, some of which have a bitter, salty, and undetectable taste. Docking simulations showed all active peptides interacted with DPP-IV through hydrogen bonds, van der Waals force, hydrophobic interaction, electrostatic force, π-sulfur, and unfavorable interaction, where WF (Trp-Phe), VW (Val-Trp), WY (Trp-Tyr), and WG (Trp-Gly) displayed higher binding affinities of 0.8; 0.5; 0.4; and 0.3 kcal/mol compared with linagliptin. In this study, we successfully demonstrated antidiabetic type-2 DM potential of the active peptides from tilapia skin collagen. The obtained data provided preliminary data for further research in the utilization of fish skin waste as a functional compound to treat the type-2 DM patients. Alternatively, this treatment can be synergistically combined with the available antidiabetic drugs to improve the insulin secretion of the type-2 DM patients.

PyVibMS: a PyMOL plugin for visualizing vibrations in molecules and solids

Visualizing vibrational motions calculated with different ab initio packages requires dedicated post-processing tools. Here, we present a PyMOL plugin called PyVibMS for visualizing the vibrational motions for both molecular and solid systems calculated by mainstream quantum chemical computer programs including Gaussian, Q-Chem, VASP, and CRYSTAL. Benefiting from the continuing development of the PyMOL platform, PyVibMS provides powerful functionalities and user-friendly interface. PyVibMS was written in Python and its open-source nature makes it flexible and sustainable. As an example, the motions of the Konkoli-Cremer local vibrational modes are shown in this work for the first time. PyVibMS is freely available at https://github.com/smutao/PyVibMS . Graphical abstract In this work, a PyMOL plugin named PyVibMS is developed to visualize molecular and lattice vibrations.

Molecular Determinants of the Promiscuity of MexB and MexY Multidrug Transporters of Pseudomonas aeruginosa

Secondary multidrug transporters of the resistance-nodulation-cell division (RND) superfamily contribute crucially to antibiotic resistance in Gram-negative bacteria. Compared to the most studied transporter AcrB of Escherichia coli, little is known about the molecular determinants of distinct polyspecificities of the most important RND transporters MexB and MexY of Pseudomonas aeruginosa. In an effort to add knowledge on this topic, we performed an exhaustive atomic-level comparison of the main putative recognition sites (access and deep binding pockets) in these two Mex transporters. We identified an underlying link between some structural, chemical and dynamical features of the binding pockets and the physicochemical nature of the corresponding substrates recognized by either one or both pumps. In particular, mosaic-like lipophilic and electrostatic surfaces of the binding pockets provide for both proteins several multifunctional sites for diffuse binding of diverse substrates. Specific lipophilicity signatures of the weakly conserved deep pocket suggest a key role of this site as a selectivity filter as in Acr transporters. Finally, the different dynamics of the bottom-loop in MexB and MexY support its possible role in binding of large substrates. Our work represents the first comparative study of the major RND transporters in P. aeruginosa and also the first structure-based study of MexY, for which no experimental structure is available yet.

Molecular Rationale behind the Differential Substrate Specificity of Bacterial RND Multi-Drug Transporters

Resistance-Nodulation-cell Division (RND) transporters AcrB and AcrD of Escherichia coli expel a wide range of substrates out of the cell in conjunction with AcrA and TolC, contributing to the onset of bacterial multidrug resistance. Despite sharing an overall sequence identity of ~66% (similarity ~80%), these RND transporters feature distinct substrate specificity patterns whose underlying basis remains elusive. We performed exhaustive comparative analyses of the putative substrate binding pockets considering crystal structures, homology models and conformations extracted from multi-copy μs-long molecular dynamics simulations of both AcrB and AcrD. The impact of physicochemical and topographical properties (volume, shape, lipophilicity, electrostatic potential, hydration and distribution of multi-functional sites) within the pockets on their substrate specificities was quantitatively assessed. Differences in the lipophilic and electrostatic potentials among the pockets were identified. In particular, the deep pocket of AcrB showed the largest lipophilicity convincingly pointing out its possible role as a lipophilicity-based selectivity filter. Furthermore, we identified dynamic features (not inferable from sequence analysis or static structures) such as different flexibilities of specific protein loops that could potentially influence the substrate recognition and transport profile. Our findings can be valuable for drawing structure (dynamics)-activity relationship to be employed in drug design.

IPP51, a chalcone acting as a microtubule inhibitor with in vivo antitumor activity against bladder carcinoma

We previously identified 1-(2,4-dimethoxyphenyl)-3-(1-methylindolyl) propenone (IPP51), a new chalcone derivative that is capable of inducing prometaphase arrest and subsequent apoptosis of bladder cancer cells. Here, we demonstrate that IPP51 selectively inhibits proliferation of tumor-derived cells versus normal non-tumor cells. IPP51 interfered with spindle formation and mitotic chromosome alignment. Accumulation of cyclin B1 and mitotic checkpoint proteins Bub1 and BubR1 on chromosomes in IPP51 treated cells indicated the activation of spindle-assembly checkpoint, which is consistent with the mitotic arrest. The antimitotic actions of other chalcones are often associated with microtubule disruption. Indeed, IPP51 inhibited tubulin polymerization in an in vitro assay with purified tubulin. In cells, IPP51 induced an increase in soluble tubulin. Furthermore, IPP51 inhibited in vitro capillary-like tube formation by endothelial cells, indicating that it has anti-angiogenic activity. Molecular docking showed that the indol group of IPP51 can be accommodated in the colchicine binding site of tubulin. This characteristic was confirmed by an in vitro competition assay demonstrating that IPP51 can compete for colchicine binding to soluble tubulin. Finally, in a human bladder xenograft mouse model, IPP51 inhibited tumor growth without signs of toxicity. Altogether, these findings suggest that IPP51 is an attractive new microtubule-targeting agent with potential chemotherapeutic value.

New pseudodimeric aurones as palm pocket inhibitors of Hepatitis C virus RNA-dependent RNA polymerase

The NS5B RNA-dependent RNA polymerase (RdRp) is a key enzyme for Hepatitis C Virus (HCV) replication. In addition to the catalytic site, this enzyme is characterized by the presence of at least four allosteric pockets making it an interesting target for development of inhibitors as potential anti-HCV drugs. Based on a previous study showing the potential of the naturally occurring aurones as inhibitors of NS5B, we pursued our efforts to focus on pseudodimeric aurones that have never been investigated so far. Hence, 14 original compounds characterized by the presence of a spacer between the benzofuranone moieties were synthesized and investigated as HCV RdRp inhibitors by means of an in vitro assay. The most active inhibitor, pseudodimeric aurone 4, induced high inhibition activity (IC50 = 1.3 μM). Mutagenic and molecular modeling studies reveal that the binding site for the most active derivatives probably is the palm pocket I instead of the thumb pocket I as for the monomeric derivatives.

Sevoflurane-Sulfobutylether-β-Cyclodextrin Complex: Preparation, Characterization, Cellular Toxicity, Molecular Modeling and Blood-Brain Barrier Transport Studies

The objective of the present investigation was to study the ability of sulfobutylether-β-cyclodextrin (SBEβCD) to form an inclusion complex with sevoflurane (SEV), a volatile anesthetic with poor water solubility. The inclusion complex was prepared, characterized and its cellular toxicity and blood-brain barrier (BBB) permeation potential of the formulated SEV have also been examined for the purpose of controlled drug delivery. The SEV-SBEβCD complex was nontoxic to the primary brain microvascular endothelial (pEND) cells at a clinically relevant concentration of sevoflurane. The inclusion complex exhibited significantly higher BBB permeation profiles as compared with the reference substance (propranolol) concerning calculated apparent permeability values (P_app). In addition, SEV binding affinity to SBEβCD was confirmed by a minimal Gibbs free energy of binding (ΔG_bind) value of −1.727 ± 0.042 kcal·mol⁻¹ and an average binding constant (K_b) of 53.66 ± 9.24 mM indicating rapid drug liberation from the cyclodextrin amphiphilic cavity.