Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug–transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.     

The Impact of P-Glycoprotein on Opioid Analgesics: What’s the Real Meaning in Pain Management and Palliative Care?

Opioids are widely used in cancer and non-cancer pain management. However, many transporters at the blood–brain barrier (BBB), such as P-glycoprotein (P-gp, ABCB1/MDR1), may impair their delivery to the brain, thus leading to opioid tolerance. Nonetheless, opioids may regulate P-gp expression, thus altering the transport of other compounds, namely chemotherapeutic agents, resulting in pharmacoresistance. Other kinds of painkillers (e.g., acetaminophen, dexamethasone) and adjuvant drugs used for neuropathic pain may act as P-gp substrates and modulate its expression, thus making pain management challenging. Inflammatory conditions are also believed to upregulate P-gp. The role of P-gp in drug–drug interactions is currently under investigation, since many P-gp substrates may also act as substrates for the cytochrome P450 enzymes, which metabolize a wide range of xenobiotics and endobiotics. Genetic variability of the ABCB1/MDR1 gene may be accountable for inter-individual variation in opioid-induced analgesia. P-gp also plays a role in the management of opioid-induced adverse effects, such as constipation. Peripherally acting mu-opioid receptors antagonists (PAMORAs), such as naloxegol and naldemedine, are substrates of P-gp, which prevent their penetration in the central nervous system. In our review, we explore the interactions between P-gp and opioidergic drugs, with their implications in clinical practice.     

CYP 2D6 Binding Affinity Predictions Using Multiple Ligand and Protein Conformations

Because of the large flexibility and malleability of Cytochrome P450 enzymes (CYPs), in silico prediction of CYP binding affinities to drugs and other xenobiotic compounds is a true challenge. In the current work, we use an iterative linear interaction energy (LIE) approach to compute CYP binding affinities from molecular dynamics (MD) simulation. In order to improve sampling of conformational space, we combine results from simulations starting with different relevant protein-ligand geometries. For calculated binding free energies of a set of thiourea compounds binding to the flexible CYP 2D6 isoform, improved correlation with experiment was obtained by combining results ofMDruns starting from distinct protein conformations and ligand-binding orientations. This accuracy was obtained from relatively short MD simulations, which makes our approach computationally attractive for automated calculations of ligand-binding affinities to flexible proteins such as CYPs.     

Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations

The novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CL^pro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CL^pro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CL^pro–PF-07321332 and 3CL^pro–α-ketoamide complexes remained stable compared with 3CL^pro–ritonavir and 3CL^pro–lopinavir. Investigating the dynamic behavior of ligand–protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41–Cys145 of 3CL^pro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CL^pro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.     

In Silico Screening of Novel α1-GABAA Receptor PAMs towards Schizophrenia Based on Combined Modeling Studies of Imidazo [1,2-a]-Pyridines

The ionotropic GABA_A receptor (GABA_AR) has been proven to be an important target of atypical antipsychotics. A novel series of imidazo [1,2-a]-pyridine derivatives, as selective positive allosteric modulators (PAMs) of α1-containing GABA_ARs with potent antipsychotic activities, have been reported recently. To better clarify the pharmacological essentiality of these PAMs and explore novel antipsychotics hits, three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) were performed on 33 imidazo [1,2-a]-pyridines. The constructed 3D-QSAR models exhibited good predictive abilities. The dockings results and MD simulations demonstrated that hydrogen bonds, π–π stackings, and hydrophobic interactions play essential roles in the binding of these novel PAMs in the GABA_AR binding pocket. Four hit compounds (DS01–04) were then screened out by the combination of the constructed models and computations, including the pharmacophore model, Topomer Search, molecular dockings, ADME/T predictions, and MD simulations. The compounds DS03 and DS04, with higher docking scores and better predicted activities, were also found to be relatively stable in the binding pocket by MD simulations. These results might provide a significant theoretical direction or information for the rational design and development of novel α1-GABA_AR PAMs with antipsychotic activities.     

Discovery of a Flavonoid FM04 as a Potent Inhibitor to Reverse P-Glycoprotein-Mediated Drug Resistance in Xenografts and Improve Oral Bioavailability of Paclitaxel

Biotransformation of flavonoid dimer FD18 resulted in an active metabolite FM04. It was more druggable because of its improved physicochemical properties. FM04 (EC₅₀ = 83 nM) was 1.8-fold more potent than FD18 in reversing P-glycoprotein (P-gp)-mediated paclitaxel (PTX) resistance in vitro. Similar to FD18, FM04 chemosensitized LCC6MDR cells towards multiple anticancer drugs by inhibiting the transport activity of P-gp and restoring intracellular drug levels. It stimulated the P-gp ATPase by 3.3-fold at 100 μM. Different from FD18, FM04 itself was not a transport substrate of P-gp and presumably, it cannot work as a competitive inhibitor. In the human melanoma MDA435/LCC6MDR xenograft, the co-administration of FM04 (28 mg/kg, I.P.) with PTX (12 mg/kg, I.V.) directly modulated P-gp-mediated PTX resistance and caused a 56% (*, p < 0.05) reduction in tumor volume without toxicity or animal death. When FM04 was administered orally at 45 mg/kg as a dual inhibitor of P-gp/CYP2C8 or 3A4 enzymes in the intestine, it increased the intestinal absorption of PTX from 0.2% to 14% in mice and caused about 57- to 66-fold improvement of AUC as compared to a single oral dose of PTX. Oral co-administration of FM04 (45 mg/kg) with PTX (40, 60 or 70 mg/kg) suppressed the human melanoma MDA435/LCC6 tumor growth with at least a 73% (***, p < 0.001) reduction in tumor volume without serious toxicity. Therefore, FM04 can be developed into a novel combination chemotherapy to treat cancer by directly targeting the P-gp overexpressed tumors or potentiating the oral bioavailability of P-gp substrate drugs.     

Conformational Landscape of Cytochrome P450 Reductase Interactions

Oxidative reactions catalyzed by Cytochrome P450 enzymes (CYPs), which constitute the most relevant group of drug-metabolizing enzymes, are enabled by their redox partner Cytochrome P450 reductase (CPR). Both proteins are anchored to the membrane of the endoplasmic reticulum and the CPR undergoes a conformational change in order to interact with the respective CYP and transfer electrons. Here, we conducted over 22 microseconds of molecular dynamics (MD) simulations in combination with protein–protein docking to investigate the conformational changes necessary for the formation of the CPR–CYP complex. While some structural features of the CPR and the CPR–CYP2D6 complex that we highlighted confirmed previous observations, our simulations revealed additional mechanisms for the conformational transition of the CPR. Unbiased simulations exposed a movement of the whole protein relative to the membrane, potentially to facilitate interactions with its diverse set of redox partners. Further, we present a structural mechanism for the susceptibility of the CPR to different redox states based on the flip of a glycine residue disrupting the local interaction network that maintains inter-domain proximity. Simulations of the CPR–CYP2D6 complex pointed toward an additional interaction surface of the FAD domain and the proximal side of CYP2D6. Altogether, this study provides novel structural insight into the mechanism of CPR–CYP interactions and underlying conformational changes, improving our understanding of this complex machinery relevant for drug metabolism.     

5-Arylidenerhodanines as P-gp Modulators: An Interesting Effect of the Carboxyl Group on ABCB1 Function in Multidrug-Resistant Cancer Cells

Multidrug resistance (MDR) is considered one of the major mechanisms responsible for the failure of numerous anticancer and antiviral chemotherapies. Various strategies to overcome the MDR phenomenon have been developed, and one of the most attractive research directions is focused on the inhibition of MDR transporters, membrane proteins that extrude cytotoxic drugs from living cells. Here, we report the results of our studies on a series newly synthesized of 5-arylidenerhodanines and their ability to inhibit the ABCB1 efflux pump in mouse T-lymphoma cancer cells. In the series, compounds possessing a triphenylamine moiety and the carboxyl group in their structure were of particular interest. These amphiphilic compounds showed over 17-fold stronger efflux pump inhibitory effects than verapamil. The cytotoxic and antiproliferative effects of target rhodanines on T-lymphoma cells were also investigated. A putative binding mode for 11, one of the most potent P-gp inhibitors tested here, was predicted by molecular docking studies and discussed with regard to the binding mode of verapamil.     

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.     

The Two-Step Clustering Approach for Metastable States Learning

Understanding the energy landscape and the conformational dynamics is crucial for studying many biological or chemical processes, such as protein–protein interaction and RNA folding. Molecular Dynamics (MD) simulations have been a major source of dynamic structure. Although many methods were proposed for learning metastable states from MD data, some key problems are still in need of further investigation. Here, we give a brief review on recent progresses in this field, with an emphasis on some popular methods belonging to a two-step clustering framework, and hope to draw more researchers to contribute to this area.     

Molecular Dynamics Simulations of Acylpeptide Hydrolase Bound to Chlorpyrifosmethyl Oxon and Dichlorvos

Acylpeptide hydrolases (APHs) catalyze the removal of N-acylated amino acids from blocked peptides. Like other prolyloligopeptidase (POP) family members, APHs are believed to be important targets for drug design. To date, the binding pose of organophosphorus (OP) compounds of APH, as well as the different OP compounds binding and inducing conformational changes in two domains, namely, α/β hydrolase and β-propeller, remain poorly understood. We report a computational study of APH bound to chlorpyrifosmethyl oxon and dichlorvos. In our docking study, Val471 and Gly368 are important residues for chlorpyrifosmethyl oxon and dichlorvos binding. Molecular dynamics simulations were also performed to explore the conformational changes between the chlorpyrifosmethyl oxon and dichlorvos bound to APH, which indicated that the structural feature of chlorpyrifosmethyl oxon binding in APH permitted partial opening of the β-propeller fold and allowed the chlorpyrifosmethyl oxon to easily enter the catalytic site. These results may facilitate the design of APH-targeting drugs with improved efficacy.     

Nanobody Paratope Ensembles in Solution Characterized by MD Simulations and NMR

Variable domains of camelid antibodies (so-called nanobodies or V_HH) are the smallest antibody fragments that retain complete functionality and therapeutic potential. Understanding of the nanobody-binding interface has become a pre-requisite for rational antibody design and engineering. The nanobody-binding interface consists of up to three hypervariable loops, known as the CDR loops. Here, we structurally and dynamically characterize the conformational diversity of an anti-GFP-binding nanobody by using molecular dynamics simulations in combination with experimentally derived data from nuclear magnetic resonance (NMR) spectroscopy. The NMR data contain both structural and dynamic information resolved at various timescales, which allows an assessment of the quality of protein MD simulations. Thus, in this study, we compared the ensembles for the anti-GFP-binding nanobody obtained from MD simulations with results from NMR. We find excellent agreement of the NOE-derived distance maps obtained from NMR and MD simulations and observe similar conformational spaces for the simulations with and without NOE time-averaged restraints. We also compare the measured and calculated order parameters and find generally good agreement for the motions observed in the ps–ns timescale, in particular for the CDR3 loop. Understanding of the CDR3 loop dynamics is especially critical for nanobodies, as this loop is typically critical for antigen recognition.     

Discovery of Novel HIV Protease Inhibitors Using Modern Computational Techniques

The human immunodeficiency virus type 1 (HIV-1) has continued to be a global concern. With the new HIV incidence, the emergence of multi-drug resistance and the untoward side effects of currently used anti-HIV drugs, there is an urgent need to discover more efficient anti-HIV drugs. Modern computational tools have played vital roles in facilitating the drug discovery process. This research focuses on a pharmacophore-based similarity search to screen 111,566,735 unique compounds in the PubChem database to discover novel HIV-1 protease inhibitors (PIs). We used an in silico approach involving a 3D-similarity search, physicochemical and ADMET evaluations, HIV protease-inhibitor prediction (IC₅₀/percent inhibition), rigid receptor–molecular docking studies, binding free energy calculations and molecular dynamics (MD) simulations. The 10 FDA-approved HIV PIs (saquinavir, lopinavir, ritonavir, amprenavir, fosamprenavir, atazanavir, nelfinavir, darunavir, tipranavir and indinavir) were used as reference. The in silico analysis revealed that fourteen out of the twenty-eight selected optimized hit molecules were within the acceptable range of all the parameters investigated. The hit molecules demonstrated significant binding affinity to the HIV protease (PR) when compared to the reference drugs. The important amino acid residues involved in hydrogen bonding and п-п stacked interactions include ASP25, GLY27, ASP29, ASP30 and ILE50. These interactions help to stabilize the optimized hit molecules in the active binding site of the HIV-1 PR (PDB ID: 2Q5K). HPS/002 and HPS/004 have been found to be most promising in terms of IC₅₀/percent inhibition (90.15%) of HIV-1 PR, in addition to their drug metabolism and safety profile. These hit candidates should be investigated further as possible HIV-1 PIs with improved efficacy and low toxicity through in vitro experiments and clinical trial investigations.     

A Lumenal Loop Associated with Catalytic Asymmetry in Plant Vacuolar H+-Translocating Pyrophosphatase

Membrane-integral inorganic pyrophosphatases (mPPases) couple pyrophosphate hydrolysis with H⁺ and Na⁺ pumping in plants and microbes. mPPases are homodimeric transporters with two catalytic sites facing the cytoplasm and demonstrating highly different substrate-binding affinities and activities. The structural aspects of the functional asymmetry are still poorly understood because the structure of the physiologically relevant dimer form with only one active site occupied by the substrate is unknown. We addressed this issue by molecular dynamics (MD) simulations of the H⁺-transporting mPPase of Vigna radiata, starting from its crystal structure containing a close substrate analog (imidodiphosphate, IDP) in both active sites. The MD simulations revealed pre-existing subunit asymmetry, which increased upon IDP binding to one subunit and persisted in the fully occupied dimer. The most significant asymmetrical change caused by IDP binding is a ‘rigid body’-like displacement of the lumenal loop connecting α-helices 2 and 3 in the partner subunit and opening its exit channel for water. This highly conserved 14–19-residue loop is found only in plant vacuolar mPPases and may have a regulatory function, such as pH sensing in the vacuole. Our data define the structural link between the loop and active sites and are consistent with the published structural and functional data.     

Dynophore-Based Approach in Virtual Screening: A Case of Human DNA Topoisomerase IIα

In this study, we utilized human DNA topoisomerase IIα as a model target to outline a dynophore-based approach to catalytic inhibitor design. Based on MD simulations of a known catalytic inhibitor and the native ATP ligand analog, AMP-PNP, we derived a joint dynophore model that supplements the static structure-based-pharmacophore information with a dynamic component. Subsequently, derived pharmacophore models were employed in a virtual screening campaign of a library of natural compounds. Experimental evaluation identified flavonoid compounds with promising topoisomerase IIα catalytic inhibition and binding studies confirmed interaction with the ATPase domain. We constructed a binding model through docking and extensively investigated it with molecular dynamics MD simulations, essential dynamics, and MM-GBSA free energy calculations, thus reconnecting the new results to the initial dynophore-based screening model. We not only demonstrate a new design strategy that incorporates a dynamic component of molecular recognition, but also highlight new derivates in the established flavonoid class of topoisomerase II inhibitors.     

Melatonin Inhibits hIAPP Oligomerization by Preventing β-Sheet and Hydrogen Bond Formation of the Amyloidogenic Region Revealed by Replica-Exchange Molecular Dynamics Simulation

The pathogenesis of type 2 diabetes (T2D) is highly related to the abnormal self-assembly of the human islet amyloid polypeptide (hIAPP) into amyloid aggregates. To inhibit hIAPP aggregation is considered a promising therapeutic strategy for T2D treatment. Melatonin (Mel) was reported to effectively impede the accumulation of hIAPP aggregates and dissolve preformed fibrils. However, the underlying mechanism at the atomic level remains elusive. Here, we performed replica-exchange molecular dynamics (REMD) simulations to investigate the inhibitory effect of Mel on hIAPP oligomerization by using hIAPP_20–29 octamer as templates. The conformational ensemble shows that Mel molecules can significantly prevent the β-sheet and backbone hydrogen bond formation of hIAPP_20–29 octamer and remodel hIAPP oligomers and transform them into less compact conformations with more disordered contents. The interaction analysis shows that the binding behavior of Mel is dominated by hydrogen bonding with a peptide backbone and strengthened by aromatic stacking and CH–π interactions with peptide sidechains. The strong hIAPP–Mel interaction disrupts the hIAPP_20–29 association, which is supposed to inhibit amyloid aggregation and cytotoxicity. We also performed conventional MD simulations to investigate the influence and binding affinity of Mel on the preformed hIAPP_1–37 fibrillar octamer. Mel was found to preferentially bind to the amyloidogenic region hIAPP_20–29, whereas it has a slight influence on the structural stability of the preformed fibrils. Our findings illustrate a possible pathway by which Mel alleviates diabetes symptoms from the perspective of Mel inhibiting amyloid deposits. This work reveals the inhibitory mechanism of Mel against hIAPP_20–29 oligomerization, which provides useful clues for the development of efficient anti-amyloid agents.     

Toward the Discovery of a Novel Class of Leads for High Altitude Disorders by Virtual Screening and Molecular Dynamics Approaches Targeting Carbonic Anhydrase

For decades, carbonic anhydrase (CA) inhibitors, most notably the acetazolamide-bearing 1,3,4-thiadiazole moiety, have been exploited at high altitudes to alleviate acute mountain sickness, a syndrome of symptomatic sensitivity to the altitude characterized by nausea, lethargy, headache, anorexia, and inadequate sleep. Therefore, inhibition of CA may be a promising therapeutic strategy for high-altitude disorders. In this study, co-crystallized inhibitors with 1,3,4-thiadiazole, 1,3-benzothiazole, and 1,2,5-oxadiazole scaffolds were employed for pharmacophore-based virtual screening of the ZINC database, followed by molecular docking and molecular dynamics simulation studies against CA to find possible ligands that may emerge as promising inhibitors. Compared to the co-crystal ligands of PDB-1YDB, 6BCC, and 6IC2, ZINC12336992, ZINC24751284, and ZINC58324738 had the highest docking scores of −9.0, −9.0, and −8.9 kcal/mol, respectively. A molecular dynamics (MD) simulation analysis of 100 ns was conducted to verify the interactions of the top-scoring molecules with CA. The system’s backbone revealed minor fluctuations, indicating that the CA–ligand complex was stable during the simulation period. Simulated trajectories were used for the MM-GBSA analysis, showing free binding energies of −16.00 ± 0.19, −21.04 ± 0.17, and −19.70 ± 0.18 kcal/mol, respectively. In addition, study of the frontier molecular orbitals of these compounds by DFT-based optimization at the level of B3LYP and the 6-311G(d,p) basis set showed negative values of the HOMO and LUMO, indicating that the ligands are energetically stable, which is essential for forming a stable ligand–protein complex. These molecules may prove to be a promising therapy for high-altitude disorders, necessitating further investigations.     

Essential dynamics of the cellular retinol-binding protein evidence for ligand-induced conformational changes

The cellular retinol-binding protein (CRBP) is an intracellularretinol carrier protein belonging to a family of hydrophobicligand–binding proteins. It transports retinol to specificlocations in the cell where, for instance, it is esterifiedfor storage. Recently solved crystallographic structures ofCRBP homologues with and without bound ligand do not provideevidence for a ligand–induced conformational change. However,it has been shown that there is a difference in binding of holo–CRBPand apo–CRBP to lecithin–retinol acyltransferase.Moreover, proteolysis of holo–CRBP and apo–CRBPyields different products, indicating a difference in structureor dynamics between the two forms. Here, we present the resultsof molecular dynamics simulations of holo–CRBP and apo–CRBP.The simulations show a significant difference in conformation,in agreement with experimental results. The essential dynamicsmethod was used to study differences in dynamics between theapo and holo forms of CRBP, and showed inhibition of essentialmotions upon ligand binding. It also revealed large correlatedmotions of retinol with regions of the protein, pointing toa possible retinol entry/exit site.     

On the Role of Temperature in the Depolymerization of PET by FAST-PETase: An Atomistic Point of View on Possible Active Site Pre-Organization and Substrate-Destabilization Effects

Enzyme FAST-PETase, recently obtained by a machine learning approach, can depolymerize poly(ethylene terephthalate) (PET), a synthetic resin employed in plastics and in clothing fibers. Therefore it represents a promising solution for the recycling of PET-based materials. In this study, a model of PET was adopted to describe the substrate, and all-atoms classical molecular dynamics (MD) simulations on apo- and substrate-bound FAST-PETase were carried out at 30 and 50 °C to provide atomistic details on the binding step of the catalytic cycle. Comparative analysis shed light on the interactions occurring between the FAST-PETase and 4PET at 50 °C, the optimal working conditions of the enzyme. Pre-organization of the enzyme active and binding sites has been highlighted, while MD simulations of FAST-PETase:4PET pointed out the occurrence of solvent-inaccessible conformations of the substrate promoted by the enzyme. Indeed, neither of these conformations was observed during MD simulations of the substrate alone in solution performed at 30, 50 and 150 °C. The analysis led us to propose that, at 50 °C, the FAST-PETase is pre-organized to bind the PET and that the interactions occurring in the binding site can promote a more reactive conformation of PET substrate, thus enhancing the catalytic activity of the enzyme.     

Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3

Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3_D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein–ligand binding affinity by employing molecular mechanics Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure–activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC₅₀ values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.