De Novo Rational Design of Peptide-Based Protein-Protein Inhibitors (Pep-PPIs) Approach by Mapping the Interaction Motifs of the PP Interface and Physicochemical Filtration: A Case on p25-Cdk5-Mediated Neurodegenerative Diseases.

Protein-protein interactions, or PPIs, are a part of every biological activity and have been linked to a number of diseases, including cancer, infectious diseases, and neurological disorders. As such, targeting PPIs is considered a strategic and vital approach in the development of new medications. Nonetheless, the wide and flat contact interface makes it difficult to find small-molecule PP inhibitors. An alternative strategy would be to use the PPI interaction motifs as building blocks for the design of peptide-based inhibitors. Herein, we designed 12-mer peptide inhibitors to target p25-inducing-cyclin-dependent kinase (Cdk5) hyperregulation, a PPI that has been shown to perpetuate neuroinflammation, which is one of the major causal implications of neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, and frontotemporal dementia. We generated a library of 5 062 500 peptide combination sequences (PCS) derived from the interaction motif of Cdk5/p25 PP interface. The 20 amino acids were differentiated into six groups, namely, hydrophobic (aliphatic), aromatic, basic, acidic, unique, and polar uncharged, on the basis of their physiochemical properties. To preserve the interaction motif necessary for ideal binding, de novo modeling of all possible peptide sequence substitutions was considered. A set of filters, backed by the Support Vector Machine (SVM) algorithm, was then used to create a shortlisted custom peptide library that met specific bioavailability, toxicity, and therapeutic relevance, leading to a refined library of 15 PCS. A greedy algorithm and coarse-grained force field were used to predict peptide structure and folding before subsequent modeling studies. Molecular docking was performed to estimate the relative binding affinities, and out of the top hits, Pep15 was subjected to molecular dynamics simulations and binding free-energy calculations in comparison to a known peptide inhibitor with experimental data (template peptide). Interestingly, the identified peptide through our protocol, Pep15, was found to show a significantly higher binding affinity than the reference template peptide (-48.10 ± 0.23 kcal/mol and -17.53 ± 0.27 kcal/mol, respectively). In comparison to the template peptide, Pep15 was found to possess a more compact and buried surface area, tighter binding landscape, and reduced conformational variability, leading to enhanced structural and kinetic stability of the Cdk5/p25 complex. Notably, both peptide inhibitors were found to have a minimal impact on the architectural integrity of the Cdk5/p25 secondary structure. Herein, we propose Pep15 as a novel and potentially disruptive peptide drug for Cdk5/p25-mediated neurodegenerative phenotypes that require further clinical investigation. The systematic protocol and findings of this report would serve as a valuable tool in the identification of critical PPI interface reactive residues, designing of analogs, and identification of more potent peptide-based PPI inhibitors.

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.

Heat shock protein (Hsp27)-ceramide synthase (Cers1) protein-protein interactions provide a new avenue for unexplored anti-cancer mechanism and therapy.

Hsp27 is a member of the small heat-shock proteins (sHSPs) - the known cellular line of defence against abnormal protein folding behaviors. Nevertheless, its upregulation is linked to a variety of pathological disorders, including several types of cancers. The ceramide synthases (CerS) mediate the synthesis of ceramide, a critical structural and signaling lipid. Functionally, downstream ceramide metabolites are implicated in the apoptosis process and their abnormal functionality has been linked to anticancer resistance. Studies showed that CerS1 are possibly inhibited by Hsp27 leading to biochemical anticancer effects in vitro. Nevertheless, the nature of such protein-protein interaction (PPI) has not been considerably investigated in molecular terms, hence, we present the first description of the dynamics CerS1-Hsp27 interaction landscapes using molecular dynamics simulations. Time-scale molecular dynamics simulation analysis indicated a system-wide conformational events of decreased stability, increased flexibility, reduced compactness, and decreased folding of CerS1. Analysis of binding energy showed a favorable interaction entailing 56 residues at the interface and a total stabilizing energy of -158 KJ/mol. The CerS1 catalytic domain experienced an opposite trend compared to the protein backbone. Yet, these residues adopted a highly compact conformation as per DCCM and DSSP analysis. Furthermore, conserved residues (SER 212, ASP 213, ALA 240, GLY 243, ASP 319) comprising the substrate shuttling machinery showed notable rigidity implying a restrained ceramide precursor access and assembly; hence, a possible inhibitory mechanism. Findings from this report would streamline a better molecular understanding of CerS1-Hsp27 interactions and decipher its potential avenue toward unexplored anti-cancer mechanisms and therapy.

Exploring the composition of protein-ligand binding sites for cancerous inhibitor of PP2A (CIP2A) by inhibitor guided binding analysis: paving a new way for the Discovery of drug candidates against triple negative breast cancer (TNBC).

Triple-negative breast cancer (TNBC) is associated with high-grade invasive carcinoma leading to a 10% to 15% death rate in younger premenopausal women. Targeting cancerous inhibitors of protein phosphatase (CIP2A) has been a highly effective approach for exploring therapeutic drug candidates. Lapatinib, a dual tyrosine kinase inhibitor, has shown promising inhibition properties by inducing apoptosis in TNBC carcinogenesis in vivo. Despite knowledge of the 3D structure of CIP2A, no reports provide insight into CIP2A ligand binding sites. To this effect, we conducted in silico site identification guided by lapatinib binding. Four of the five sites identified were cross-validated, and the stem domain revealed more excellent ligand binding affinity. The binding affinity of lapatinib in these sites was further computed using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) approach. According to MM/PBSA//200 ns MD simulations, lapatinib exhibited a higher binding affinity against CIP2A in site 2 with ΔG critical values of -37.1 kcal/mol. The steadiness and tightness of lapatinib with CIP2A inside the stem domain disclosed glutamic acid-318 as the culprit amino acid with the highest electrostatic energy. These results provide clear information on the CIP2A domain capable of ligand binding and validate lapatinib as a promising CIP2A inhibitor in TNBC carcinogenesis.

The Unusual Architecture of RNA-Dependent RNA Polymerase (RdRp)'s Catalytic Chamber Provides a Potential Strategy for Combination Therapy against COVID-19.

The unusual and interesting architecture of the catalytic chamber of the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) was recently explored using Cryogenic Electron Microscopy (Cryo-EM), which revealed the presence of two distinctive binding cavities within the catalytic chamber. In this report, first, we mapped out and fully characterized the variations between the two binding sites, BS1 and BS2, for significant differences in their amino acid architecture, size, volume, and hydrophobicity. This was followed by investigating the preferential binding of eight antiviral agents to each of the two binding sites, BS1 and BS2, to understand the fundamental factors that govern the preferential binding of each drug to each binding site. Results showed that, in general, hydrophobic drugs, such as remdesivir and sofosbuvir, bind better to both binding sites than relatively less hydrophobic drugs, such as alovudine, molnupiravir, zidovudine, favilavir, and ribavirin. However, suramin, which is a highly hydrophobic drug, unexpectedly showed overall weaker binding affinities in both binding sites when compared to other drugs. This unexpected observation may be attributed to its high binding solvation energy, which disfavors overall binding of suramin in both binding sites. On the other hand, hydrophobic drugs displayed higher binding affinities towards BS1 due to its higher hydrophobic architecture when compared to BS2, while less hydrophobic drugs did not show a significant difference in binding affinities in both binding sites. Analysis of binding energy contributions revealed that the most favorable components are the ΔEele, ΔEvdw, and ΔGgas, whereas ΔGsol was unfavorable. The ΔEele and ΔGgas for hydrophobic drugs were enough to balance the unfavorable ΔGsol, leaving the ΔEvdw to be the most determining factor of the total binding energy. The information presented in this report will provide guidelines for tailoring SARS-CoV-2 inhibitors with enhanced binding profiles.

Unusual chalcogen⋯chalcogen interactions in like⋯like and unlike YCY⋯YCY complexes (Y = O, S, and Se).

Chalcogen⋯chalcogen interactions were investigated within four types of like⋯like and unlike YCY⋯YCY complexes (where Y = O, S, or Se). A plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), surface electrostatic potential extrema, point-of-charge (PoC), quantum theory of atoms in molecules (QTAIM), noncovalent interaction (NCI), and symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) calculations, were executed. The energetic findings revealed a preferential tendency of the studied chalcogen-bearing molecules to engage in type I, II, III, or IV chalcogen⋯chalcogen interactions. Notably, the selenium-bearing molecules exhibited the most potent ability to favorably participate in all the explored chalcogen⋯chalcogen interactions. Among like⋯like complexes, type IV interactions showed the most favorable negative binding energies, whereas type III interactions exhibited the weakest binding energies. Unexpectedly, oxygen-containing complexes within type IV interactions showed an alien pattern of binding energies that decreased along with an increase in the chalcogen atomic size level. QTAIM analysis provided a solo BCP, via chalcogen⋯chalcogen interactions, with no clues as to any secondary ones. SAPT-EDA outlined the domination of the explored interactions by the dispersion forces and indicated the pivotal shares of the electrostatic forces, except type III σ-hole⋯σ-hole and di-σ-hole interactions. These observations demonstrate in better detail all the types of chalcogen⋯chalcogen interactions, providing persuasive reasons for their more intensive use in versatile fields related to materials science and drug design.

Naturally occurring plant-based anticancerous candidates as prospective ABCG2 inhibitors: an in silico drug discovery study.

ATP-binding cassette transporter G2 (ABCG2) is an efflux transporter related to the clinical multidrug resistance (MDR) phenomenon. Identifying ABCG2 inhibitors could help discover extraordinary curative strategies for carcinoma remediation. Hitherto, there is no medication drug inhibiting ABCG2 transporter, notwithstanding that a considerable number of drugs have been submitted to clinical-trial and investigational phases. In the search for unprecedented chemical compounds that could inhibit the ABCG2 transporter, an in silico screening was conducted on the Naturally Occurring Plant-based Anticancer Compound-Activity-Target (NPACT) database containing 1574 compounds. Inhibitor-ABCG2 binding affinities were estimated based on molecular docking and molecular minimization (MM) calculations and compared to a co-crystallized inhibitor (BWQ) acting as a reference inhibitor. Molecular dynamics (MD) simulations pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations were further executed for compounds with MM-GBSA//MM binding energies lower than BWQ (calc. - 60.5 kcal/mol). NPACT00968 and NPACT01545 demonstrated auspicious inhibitory activities according to binding affinities (ΔGbinding) over the 100 ns MD simulations that were nearly one and a half folds compared to BWQ (- 100.4, - 94.7, and - 62.9 kcal/mol, respectively). Throughout the 100 ns MD simulations, structural and energetical analyses unveiled outstanding stability of the ABCG2 transporter when bound with NPACT00968 and NPACT01545. In silico calculations hold a promise for those two inhibitors as drug candidates of ABCG2 transporter and emphasize that further in vitro and in vivo experiments are guaranteed.

Novel Sunifiram-carbamate hybrids as potential dual acetylcholinesterase inhibitor and NMDAR co-agonist: simulation-guided analogue design and pharmacological screening.

An efficient method for synthesising NMDAR co-agonist Sunifiram (DM235), in addition to Sunifram-carbamate and anthranilamide hybrids, has been developed in high yields via protecting group-free stepwise unsymmetric diacylation of piperazine using N-acylbenzotiazole. Compounds 3f, 3d, and 3i exhibited promising nootropic activity by enhancing acetylecholine (ACh) release in A549 cell line. Moreover, the carbamate hybrid 3f was found to exhibit higher in vitro potency than donepezil with IC50 = 18 ± 0.2 nM, 29.9 ± 0.15 nM for 3f and donepezil, respectively. 3f was also found to effectively inhibit AChE activity in rat brain (AChE = 1.266 ng/mL) compared to tacrine (AChE = 1.137 ng/ml). An assessment of the ADMET properties revealed that compounds 3f, 3d, and 3i are drug-like and can penetrate blood-brain barrier. Findings presented here showcase highly potential cholinergic agents, with expected partial agonist activity towards glycine binding pocket of NMDAR which could lead to development and optimisation of novel nootropic drugs.

Co-Binding of JQ1 and Venetoclax Exhibited Synergetic Inhibitory Effect for Cancer Therapy; Potential Line of Treatment for the Waldenström Macroglobulinemia Lymphoma.

In recent times, the development of combination therapy has been a focal point in drug discovery. This article explores the potential synergistic effect of co-administration of Bcl2 inhibitor Venetoclax and BET inhibitor JQ1. We envisioned that the 'dual-site'-binding of Bcl2 has significant prospects and paves the way for the next round of rational design of potent Waldenström macroglobulinemia (WM) therapy. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. This study conducted a systematic procedure to explore the preferred binding modes and the structural effects of co-binding at each catalytic active site. Interestingly, a mutual enhanced binding effect was observed - Venetoclax increased the binding affinity of JQ1 by 11.5 %, while JQ1 boosted the binding affinity of Venetoclax by 16.3 % when compared with individual inhibition of each drug. This synergistic binding effect has significantly increased protein stability, with substantial correlated movements and multiple van der Waals interactions. The structural and thermodynamic insights unveiled in this report would assist the future design of improved combined therapy against WM.

An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape.

Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.

Type I-IV Halogen⋯Halogen Interactions: A Comparative Theoretical Study in Halobenzene⋯Halobenzene Homodimers.

In the current study, unexplored type IV halogen⋯halogen interaction was thoroughly elucidated, for the first time, and compared to the well-established types I−III interactions by means of the second-order Møller−Plesset (MP2) method. For this aim, the halobenzene⋯halobenzene homodimers (where halogen = Cl, Br, and I) were designed into four different types, parodying the considered interactions. From the energetic perspective, the preference of scouted homodimers was ascribed to type II interactions (i.e., highest binding energy), whereas the lowest binding energies were discerned in type III interactions. Generally, binding energies of the studied interactions were observed to decline with the decrease in the σ-hole size in the order, C6H5I⋯IC6H5 > C6H5Br⋯BrC6H5 > C6H5Cl⋯ClC6H5 homodimers and the reverse was noticed in the case of type IV interactions. Such peculiar observations were relevant to the ample contributions of negative-belt⋯negative-belt interactions within the C6H5Cl⋯ClC6H5 homodimer. Further, type IV torsional trans → cis interconversion of C6H5X⋯XC6H5 homodimers was investigated to quantify the π⋯π contributions into the total binding energies. Evidently, the energetic features illustrated the amelioration of the considered homodimers (i.e., more negative binding energy) along the prolonged scope of torsional trans → cis interconversion. In turn, these findings outlined the efficiency of the cis configuration over the trans analog. Generally, symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) demonstrated the predominance of all the scouted homodimers by the dispersion forces. The obtained results would be beneficial for the omnipresent studies relevant to the applications of halogen bonds in the fields of materials science and crystal engineering.

Effects of Lewis Basicity and Acidity on σ-Hole Interactions in Carbon-Bearing Complexes: A Comparative Ab Initio Study.

The effects of Lewis basicity and acidity on σ-hole interactions were investigated using two sets of carbon-containing complexes. In Set I, the effect of Lewis basicity was studied by substituting the X3/X atom(s) of the NC-C6H2-X3 and NCX Lewis bases (LB) with F, Cl, Br, or I. In Set II, the W-C-F3 and F-C-X3 (where X and W = F, Cl, Br, and I) molecules were utilized as Lewis acid (LA) centers. Concerning the Lewis basicity effect, higher negative interaction energies (Eint) were observed for the F-C-F3∙∙∙NC-C6H2-X3 complexes compared with the F-C-F3∙∙∙NCX analogs. Moreover, significant Eint was recorded for Set I complexes, along with decreasing the electron-withdrawing power of the X3/X atom(s). Among Set I complexes, the highest negative Eint was ascribed to the F-C-F3∙∙∙NC-C6H2-I3 complex with a value of -1.23 kcal/mol. For Set II complexes, Eint values of F-C-X3 bearing complexes were noted within the -1.05 to -2.08 kcal/mol scope, while they ranged from -0.82 to -1.20 kcal/mol for the W-C-F3 analogs. However, Vs,max quantities exhibited higher values in the case of W-C-F3 molecules compared with F-C-X3; preferable negative Eint were ascribed to the F-C-X3 bearing complexes. These findings were delineated as a consequence of the promoted contributions of the X3 substituents. Dispersion forces (Edisp) were identified as the dominant forces for these interactions. The obtained results provide a foundation for fields such as crystal engineering and supramolecular chemistry studies that focus on understanding the characteristics of carbon-bearing complexes.

Bridging the Gap in Malaria Parasite Resistance, Current Interventions, and the Way Forward from in Silico Perspective: A Review.

The past decade has seen most antimalarial drugs lose their clinical potency stemming from parasite resistance. Despite immense efforts by researchers to mitigate this global scourge, a breakthrough is yet to be achieved, as most current malaria chemotherapies suffer the same fate. Though the etiology of parasite resistance is not well understood, the parasite's complex life has been implicated. A drug-combination therapy with artemisinin as the central drug, artemisinin-based combination therapy (ACT), is currently the preferred malaria chemotherapy in most endemic zones. The emerging concern of parasite resistance to artemisinin, however, has compromised this treatment paradigm. Membrane-bound Ca2+-transporting ATPase and endocytosis pathway protein, Kelch13, among others, are identified as drivers in plasmodium parasite resistance to artemisinin. To mitigate parasite resistance to current chemotherapy, computer-aided drug design (CADD) techniques have been employed in the discovery of novel drug targets and the development of small molecule inhibitors to provide an intriguing alternative for malaria treatment. The evolution of plasmepsins, a class of aspartyl acid proteases, has gained tremendous attention in drug discovery, especially the non-food vacuole. They are expressed at multi-stage of the parasite's life cycle and involve in hepatocytes' egress, invasion, and dissemination of the parasite within the human host, further highlighting their essentiality. In silico exploration of non-food vacuole plasmepsin, PMIX and PMX unearthed the dual enzymatic inhibitory mechanism of the WM382 and 49c, novel plasmepsin inhibitors presently spearheading the search for potent antimalarial. These inhibitors impose structural compactness on the protease, distorting the characteristic twist motion. Pharmacophore modeling and structure activity of these compounds led to the generation of hits with better affinity and inhibitory prowess towards PMIX and PMX. Despite these headways, the major obstacle in targeting PM is the structural homogeneity among its members and to human Cathepsin D. The incorporation of CADD techniques described in the study at early stages of drug discovery could help in selective inhibition to augment malaria chemotherapy.

External Electric Field Effect on the Strength of σ-Hole Interactions: A Theoretical Perspective in Like⋯Like Carbon-Containing Complexes.

For the first time, σ-hole interactions within like⋯like carbon-containing complexes were investigated, in both the absence and presence of the external electric field (EEF). The effects of the directionality and strength of the utilized EEF were thoroughly unveiled in the (F-C-F3)2, (F-C-H3)2, and (H-C-F3)2 complexes. In the absence of the EEF, favorable interaction energies, with negative values, are denoted for the (F-C-F3)2 and (H-C-F3)2 complexes, whereas the (F-C-H3)2 complex exhibits unfavorable interactions. Remarkably, the strength of the applied EEF exhibits a prominent role in turning the repulsive forces within the latter complex into attractive ones. The symmetrical nature of the considered like⋯like carbon-containing complexes eradicated the effect of directionality of the EEF. The quantum theory of atoms in molecules (QTAIM), and the noncovalent interaction (NCI) index, ensured the occurrence of the attractive forces, and also outlined the substantial contributions of the three coplanar atoms to the total strength of the studied complexes. Symmetry-adapted perturbation theory (SAPT) results show the dispersion-driven nature of the interactions.

Exploring Natural Product Activity and Species Source Candidates for Hunting ABCB1 Transporter Inhibitors: An In Silico Drug Discovery Study.

The P-glycoprotein (P-gp/ABCB1) is responsible for a xenobiotic efflux pump that shackles intracellular drug accumulation. Additionally, it is included in the dud of considerable antiviral and anticancer chemotherapies because of the multidrug resistance (MDR) phenomenon. In the search for prospective anticancer drugs that inhibit the ABCB1 transporter, the Natural Product Activity and Species Source (NPASS) database, containing >35,000 molecules, was explored for identifying ABCB1 inhibitors. The performance of AutoDock4.2.6 software to anticipate ABCB1 docking score and pose was first assessed according to available experimental data. The docking scores of the NPASS molecules were predicted against the ABCB1 transporter. Molecular dynamics (MD) simulations were conducted for molecules with docking scores lower than taxol, a reference inhibitor, pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. On the basis of MM-GBSA calculations, five compounds revealed promising binding affinities as ABCB1 inhibitors with ΔGbinding < −105.0 kcal/mol. The binding affinity and stability of the identified inhibitors were compared to the chemotherapeutic agent. Structural and energetical analyses unveiled great steadiness of the investigated inhibitors within the ABCB1 active site throughout 100 ns MD simulations. Conclusively, these findings point out that NPC104372, NPC475164, NPC2313, NPC197736, and NPC477344 hold guarantees as potential ABCB1 drug candidates and warrant further in vitro/in vivo tests.

Exploring the structural basis and atomistic binding mechanistic of the selective antagonist blockade at D3 dopamine receptor over D2 dopamine receptor.

More recently, there has been a paradigm shift toward selective drug targeting in the treatment of neurological disorders, including drug addiction, schizophrenia, and Parkinson's disease mediated by the different dopamine receptor subtypes. Antagonists with higher selectivity for D3 dopamine receptor (D3DR) over D2 dopamine receptor (D2DR) have been shown to attenuate drug-seeking behavior and associated side effects compared to non-subtype selective antagonists. However, high conservations among constituent residues of both proteins, particularly at the ligand-binding pockets, remain a challenge to therapeutic drug design. Recent studies have reported the discovery of two small-molecules R-VK4-40 and Y-QA31 which substantially inhibited D3DR with >180-fold selectivity over D2DR. Therefore, in this study, we seek to provide molecular and structural insights into these differential binding mechanistic using meta-analytic computational simulation methods. Findings revealed that R-VK4-40 and Y-QA31 adopted shallow binding modes and were more surface-exposed at D3DR while on the contrary, they exhibited deep hydrophobic pocket binding at D2DR. Also, two non-conserved residues; Tyr361.39 and Ser18245.51 were identified in D3DR, based on their crucial roles and contributions to the selective binding of R-VK4-40 and Y-QA31. Importantly, both antagonists exhibited high affinities in complex with D3DR compared to D2DR, while van der Waals energies contributed majorly to their binding and stability. Structural analyses also revealed the distinct stabilizing effects of both compounds on D3DR secondary architecture relative to D2DR. Therefore, findings herein pinpointed the origin and mechanistic of selectivity of the compounds, which may assist in the rational design of potential small molecules of the D2 -like dopamine family receptor subtype with improved potency and selectivity.

Dual-Knockout of Mutant Isocitrate Dehydrogenase 1 and 2 Subtypes Towards Glioma Therapy: Structural Mechanistic Insights on the Role of Vorasidenib.

Recently, Vorasidenib (AG-881) has been reported as a therapeutic alternative that exerts potent dual inhibitory activity against mIDH1/2 towards the treatment of low-grade glioma. However, structural and dynamic events associated with its dual inhibition mechanism remain unclear. As such, we employ integrative computer-assisted atomistic techniques to provide thorough structural and dynamic insights. Our analysis proved that the dual-targeting ability of AG-881 is mediated by Val255/Val294 within the binding pockets of both mIDH1 and mIDH2 which are shown to elicit a strong intermolecular interaction, thus favoring binding affinity. The structural orientations of AG-881 within the respective hydrophobic pockets allowed favorable interactions with binding site residues which accounted for its high binding free energy of -28.69 kcal/mol and -19.89 kcal/mol towards mIDH1 and mIDH2, respectively. Interestingly, upon binding, AG-881 was found to trigger systemic alterations of mIDH1 and mIDH2 characterized by restricted residue flexibility and a reduction in exposure of residues to the solvent surface area. As a result of these structural alterations, crucial interactions of the mutant enzymes were inhibited, a phenomenon that results in a suppression of the production of oncogenic stimulator 2-HG. Findings therefore provide thorough structural and dynamic insights associated with the dual inhibitory activity of AG-881 towards glioma therapy.

Mechanistic Insights into the Selective Dual BET and PLK1 Inhibitory Activity of a Novel Benzamide Compound in Castration-Resistant Prostrate Cancer.

Though multifactorial, BET and PLK1 proteins have been found to be key players in the oncogenic process leading to castration-resistant prostate cancer through regulation of AR and MYC-mediated transcription. Hence, dual inhibition of these proteins appears to be an auspicious approach for CRPC therapy. WNY0824 has been reported to exhibit nanomolar range inhibition as well as significant anti-proliferative activity on AR-positive CRPC cells in vitro. However, structural, and mechanistic events associated with its dual inhibitory and anti-proliferative mechanisms remain unclear. Utilizing integrative computer-assisted atomistic techniques, analyses revealed that the dual-inhibitory activity of WNY0824 against BRD4 and PLK1 proteins is mediated by conserved residues present in the binding cavities of both proteins which are shown to elicit various strong intermolecular interactions and thus favour binding affinity. Also, binding orientation of the ligand at the protein binding cavities allowed for important hydrophobic interactions which resulted in high binding free energy of -42.50 kcal/mol and -51.64 kcal/mol towards BRD4 and PLK1, respectively. While van der Waals interactions are very important to ligand binding in BRD4-WNY complex, electrostatic interactions are pertinent to PLK1-WNY complex. Intriguingly, WNY0824 triggered conformational alterations in both proteins through increased structural instability, decreased structural compactness and mitigation in exposure of residues to solvent surface area. Consequently, critical interactions peculiar to the oncogenic activities of BRD4 and PLK1 were inhibited, a phenomenon that results in an antagonism of CRPC progression. The mechanistic insights presented in this report would further assist in the structure-based design of improved inhibitors useful in CRPC therapy.

Insight into the Therapeutic Potential of a Bicyclic Hydroxypyridone Compound 2-[(2,4-Dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one as COMT Inhibitor in the Treatment of Parkinson's Disease: A Molecular Dynamic Simulation Approach.

Parkinson's disease (PD) is one of the most targeted neurodegenerative diseases in clinical research. Awareness of research is due to its increasing number of affected people worldwide. The pathology of PD has been linked to several key proteins upregulation such as the catechol O-Methyltransferase (COMT). Hence, the synthesis of compounds possessing inhibitory capacity has been the frontline of research in recent years. Several compounds have been synthesized among which is the nitrocatechol. However, major limitations associated with the nitrocatechol scaffold include the inability to possess adequate CNS penetration properties and hepatic toxicity associated with the compounds. However, a series of bicyclic hydroxypyridones compounds were synthesized to evaluate their inhibitory potentials on COMT protein with compound 38 (c38) 2-[(2,4-dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one shown to have a 40 fold increase level coverage in its IC50 over brain exposure when compared to the other synthesized compound. The molecular dynamics method was employed to understand the nature of interaction exhibited by c38. Molecular mechanics of c38 revealed a disruptive effect on the secondary structure of COMT protein. Per residue decomposition analysis revealed similar crucial residues involved in the favorable binding of c38 and tolcapone implicated its increased inhibitory capacity on COMT in preventing PD. Free binding energy (ΔGbind ) of c38 further revealed the inhibitory capacity towards COMT protein in comparison to the FDA approved tolcapone. Ligand mobility analysis of both compounds showed a timewise different mobility pattern across the simulation time frame at the active site pocket of the protein connoting the different inhibitory potency exhibited by c38 and tolcapone. Findings from this study revealed optimization of c38 could facilitate the discovery of new compounds with enhanced inhibitory properties towards COMT in treating PD.

Investigating the Mechanistic Inhibitory Discrepancies of Novel Halogen and Alkyl Di-Substituted Oxadiazole-Based Dibenzo-Azepine-Dione Derivatives on Poly (ADP-Ribose) Polymerase-1.

Numerous studies have established the involvement of Poly (ADP-ribose) Polymerase-1 (PARP-1) in cancer presenting it as an important therapeutic target over recent years. Although homology among the PARP protein family makes selective targeting difficult, two compounds [d11 (0.939 µM) and d21 (0.047 µM)] with disparate inhibitory potencies against PARP-1 were recently identified. In this study, free energy calculations and molecular simulations were used to decipher underlying mechanisms of differential PARP-1 inhibition exhibited by the two compounds. The thermodynamics calculation revealed that compound d21 had a relatively higher ΔGbind than d11. High involvement of van der Waal and electrostatic effects potentiated the affinity of d21 at PARP-1 active site. More so, incorporated methyl moiety in d11 accounted for steric hindrance which, in turn, prevented complementary interactions of key site residues such as TYR889, MET890, TYR896, TYR907. Conformational studies also revealed that d21 is more stabilized for interactions in the active site compared to d11. We believe that findings from this study would provide an important avenue for the development of selective PARP-1 inhibitors.