Ensemble Docking in Drug Discovery: How Many Protein Configurations from Molecular Dynamics Simulations are Needed To Reproduce Known Ligand Binding?

Exploring the Potential of Oleanolic Acid Dimers-Cytostatic and Antioxidant Activities, Molecular Docking, and ADMETox Profile

The presented work aimed to explore the potential of oleanolic acid dimers (OADs): their cytostatic and antioxidant activities, molecular docking, pharmacokinetics, and ADMETox profile. The cytostatic properties of oleanolic acid (1) and its 14 synthesised dimers (2a-2n) were evaluated against 10 tumour types and expressed as IC50 values. Molecular docking was performed with the CB-Dock2 server. Antioxidant properties were evaluated with the CUPRAC method. ADMETox properties were evaluated with the ADMETlab Manual (2.0) database. The results indicate that the obtained OADs can be effective cytostatic agents, for which the IC50 not exceeded 10.00 for many tested cancer cell lines. All OADs were much more active against all cell lines than the mother compound (1). All dimers can inhibit the interaction between the 1MP8 protein and cellular proteins with the best results for compounds 2f and 2g with unsaturated bonds within the linker. An additional advantage of the tested OADs was a high level of antioxidant activity, with Trolox equivalent for OADs 2c, 2d, 2g-2j, 2l, and 2m of approximately 0.04 mg/mL, and beneficial pharmacokinetics and ADMETox properties. The differences in the DPPH and CUPRAC assay results obtained for OADs may indicate that these compounds may be effective antioxidants against different radicals.

A Robust Induced Fit Docking Approach with the Combination of the Hybrid All-Atom/United-Atom/Coarse-Grained Model and Simulated Annealing

Molecular docking remains an indispensable tool in computational biology and structure-based drug discovery. However, the correct prediction of binding poses remains a major challenge for molecular docking, especially for target proteins where a substrate binding induces significant reorganization of the active site. Here, we introduce an Induced Fit Docking (IFD) approach named AA/UA/CG-SA-IFD, which combines a hybrid All-Atom/United-Atom/Coarse-Grained model with Simulated Annealing. In this approach, the core region is represented by the All-Atom(AA) model, while the protein environment beyond the core region and the solvent are treated with either the United-Atom (UA) or the Coarse-Grained (CG) model. By combining the Elastic Network Model (ENM) for the CG region, the hybrid model ensures a reasonable description of ligand binding and the environmental effects of the protein, facilitating highly efficient and reliable sampling of ligand binding through Simulated Annealing (SA) at a high temperature. Upon validation with two testing sets, the AA/UA/CG-SA-IFD approach demonstrates remarkable accuracy and efficiency in induced fit docking, even for challenging cases where the docked poses significantly deviate from crystal structures.

Exploring the interaction between a fluorescent Ag(I)-biscarbene complex and non-canonical DNA structures: a multi-technique investigation

Silver compounds are mainly studied as antimicrobial agents, but they also have anticancer properties, with the latter, in some cases, being better than their gold counterparts. Herein, we analyse the first example of a new Ag(I)-biscarbene that can bind non-canonical structures of DNA, more precisely G-quadruplexes (G4), with different binding signatures depending on the type of G4. Moreover, we show that this Ag-based carbene binds the i-motif DNA structure. Alternatively, its Au(I) counterpart, which was investigated for comparison, stabilises mitochondrial G4. Theoretical in silico studies elucidated the details of different binding modes depending on the geometry of G4. The two complexes showed increased cytotoxic activity compared to cisplatin, overcoming its resistance in ovarian cancer. The binding of these new drug candidates with other relevant biosubstrates was studied to afford a more complete picture of their possible targets. In particular, the Ag(I) complex preferentially binds DNA structures over RNA structures, with higher binding constants for the non-canonical nucleic acids with respect to natural calf thymus DNA. Regarding possible protein targets, its interaction with the albumin model protein BSA was also tested.

Repurposing of drug molecules from FDA database against Hepatitis C virus E2 protein using ensemble docking approach

Hepatitis C virus, a member of the Flaviviridae family and genus Hepacivirus, is an enveloped, positively single stranded RNA virus. Its surface consists of a heterodimer of E1 and E2 proteins which play a crucial role in receptor binding and membrane fusion. In this study we have used in silico virtual screening by utilizing ensemble docking on the approved drugs. These drugs can bind with high efficiency to the 36 prominent conformations of the CD81 binding site clustered from a total of 3 µs MD simulation data on the E2 protein. We started with 9213 compounds from the FDA list of drugs and progressively came down to 5 compounds which have been seen to bind with very high efficiency to not only all the conformations but also the two predicted druggable pockets that encompass the CD81 binding site. MM/PBSA binding energy calculations also point to the highly stable interaction of the compounds to the E2 protein. This study may in future broaden the arsenal of therapeutics for use against HCV infection and lead to more effective care against the virus.

DiffBindFR: an SE(3) equivariant network for flexible protein-ligand docking

Molecular docking, a key technique in structure-based drug design, plays pivotal roles in protein-ligand interaction modeling, hit identification and optimization, in which accurate prediction of protein-ligand binding mode is essential. Conventional docking approaches perform well in redocking tasks with known protein binding pocket conformation in the complex state. However, in real-world docking scenario without knowing the protein binding conformation for a new ligand, accurately modeling the binding complex structure remains challenging as flexible docking is computationally expensive and inaccurate. Typical deep learning-based docking methods do not explicitly consider protein side chain conformations and fail to ensure the physical plausibility and detailed atomic interactions. In this study, we present DiffBindFR, a full-atom diffusion-based flexible docking model that operates over the product space of ligand overall movements and flexibility and pocket side chain torsion changes. We show that DiffBindFR has higher accuracy in producing native-like binding structures with physically plausible and detailed interactions than available docking methods. Furthermore, in the Apo and AlphaFold2 modeled structures, DiffBindFR demonstrates superior advantages in accurate ligand binding pose and protein binding conformation prediction, making it suitable for Apo and AlphaFold2 structure-based drug design. DiffBindFR provides a powerful flexible docking tool for modeling accurate protein-ligand binding structures.

Efficient Refinement of Complex Structures of Flexible Histone Peptides Using Post-Docking Molecular Dynamics Protocols

Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To address this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In this present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32% improvement over the docked structures in terms of the change in root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements are also discussed to help with their implementation in future methods and applications.

Tepotinib and tivantinib as potential inhibitors for the serine/threonine kinase of the mpox virus: insights from structural bioinformatics analysis

The serine/threonine kinase (STK) plays a central role as the primary kinase in poxviruses, directing phosphoryl transfer reactions. Such reactions are pivotal for the activation of certain proteins during viral replication, assembly, and maturation. Therefore, targeting this key protein is anticipated to impede virus replication. In this work, a structural bioinformatics approach was employed to evaluate the potential of drug-like kinase inhibitors in binding to the ATP-binding pocket on the STK of the Mpox virus. Virtual screening of known kinase inhibitors revealed that the top 10 inhibitors exhibited binding affinities ranging from -8.59 to -12.05 kcal/mol. The rescoring of compounds using the deep-learning default model in GNINA was performed to predict accurate binding poses. Subsequently, the top three inhibitors underwent unbiased molecular dynamics (MD) simulations for 100 ns. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) analysis and Principal Component Analysis (PCA) suggested tepotinib as a competitive inhibitor for Mpox virus STK as evidenced by its binding free energy and the induction of similar conformational behavior of the enzyme. Nevertheless, it is sensible to experimentally test all top 10 compounds, as scoring functions and energy calculations may not consistently align with experimental findings. These insights are poised to provide an attempt to identify an effective inhibitor for the Mpox virus.Communicated by Ramaswamy H. Sarma.

Spectroscopy and molecular simulation on the interaction of Nano-Kaempferol prepared by oil-in-water with two carrier proteins: An investigation of protein-protein interaction

In this work, the interaction of human serum albumin (HSA) and human holo-transferrin (HTF) with the prepared Nano-Kaempferol (Nano-KMP) through oil-in-water procedure was investigated in the form of binary and ternary systems by the utilization of different spectroscopy techniques along with molecular simulation and cancer cell experiments. According to fluorescence spectroscopy outcomes, Nano-KMP is capable of quenching both proteins as binary systems by a static mechanism, while in the form of (HSA-HTF) Nano-KMP as the ternary system, an unlinear Stern-Volmer plot was elucidated with the occurrence of both dynamic and static fluorescence quenching mechanisms in the binding interaction. In addition, the two acquired Ksv values in the ternary system signified the existence of two sets of binding sites with two different interaction behaviors. The binding constant values of HSA-Nano KMP, HTF-Nano-KMP, and (HSA-HTF) Nano-KMP complexes formation were (2.54 ± 0.03) × 104, (2.15 ± 0.02) × 104 and (1.43 ± 0.04) × 104M-1at the first set of binding sites and (4.68 ± 0.05) × 104 M-1 at the second set of binding sites, respectively. The data of thermodynamic parameters confirmed the major roles of hydrogen binding and van der Waals forces in the formation of HSA-Nano KMP and HTF-Nano KMP complexes. The thermodynamic parameter values of (HSA-HTF) Nano KMP revealed the dominance of hydrogen binding and van der Waals forces in the first set of binding sites and hydrophobic forces for the second set of binding sites. Resonance light scattering (RLS) analysis displayed the existence of a different interaction behavior for HSA-HTF complex in the presence of Nano-KMP as the ternary system. Moreover, circular dichroism (CD) technique affirmed the conformational changes of the secondary structure of proteins as binary and ternary systems. Molecular docking and molecular dynamics simulations (for 100 ns) were performed to investigate the mechanism of KMP binding to HSA, HTF, and HSA-HTF. Next to observing a concentration and time-dependent cytotoxicity, the down regulation of PI3K/AkT/mTOR pathway resulted in cell cycle arrest in SW480 cells.

An ensemble docking-based virtual screening and molecular dynamics simulation of phytochemical compounds from Malaysian Kelulut Honey (KH) against SARS-CoV-2 target enzyme, human angiotensin-converting enzyme 2 (ACE-2)

The human angiotensin-converting enzyme 2 (ACE-2) receptor is a metalloenzyme that plays an important role in regulating blood pressure by modulating angiotensin II. This receptor facilitates SARS-CoV-2 entry into human cells via receptor-mediated endocytosis, causing the global COVID-19 pandemic and a major health crisis. Kelulut honey (KH), one of Malaysian honey recently gained attention for its distinct flavour and taste while having many nutritional and medicinal properties. Recent study demonstrates the antiviral potential of KH against SARS-CoV-2 by inhibiting ACE-2 in vitro, but the bioactive compound pertaining to the ACE-2 inhibition is yet unknown. An ensemble docking-based virtual screening was employed to screen the phytochemical compounds from KH with high binding affinity against the 10 best representative structures of ACE-2 that mostly formed from MD simulation. From 110 phytochemicals previously identified in KH, 27 compounds passed the ADMET analysis and proceeded to docking. Among the docked compound, SDC and FMN consistently exhibited strong binding to ACE-2's active site (-9.719 and -9.473 kcal/mol) and allosteric site (-7.305 and -7.464 kcal/mol) as compared to potent ACE-2 inhibitor, MLN 4760. Detailed trajectory analysis of MD simulation showed stable binding interaction towards active and allosteric sites of ACE-2. KH's compounds show promise in inhibiting SARS-CoV-2 binding to ACE-2 receptors, indicating potential for preventive use or as a supplement to other COVID-19 treatments. Additional research is needed to confirm KH's antiviral effects and its role in SARS-CoV-2 therapy, including prophylaxis and adjuvant treatment with vaccination.Communicated by Ramaswamy H. Sarma.

Biochemical and molecular docking-based strategies of Acalypha indica and Boerhavia diffusa extract by targeting bacterial strains and cancer proteins

Antibiotic-resistant microbes have emerged around the world, presenting a risk to health. Plant-derived drugs have become a potential source for the production of antibiotic-resistant drugs and cancer therapies. In this study, we investigated the antibacterial, cytotoxic and antioxidant properties of Acalypha indica and Boerhavia diffusa, and conducted in silico molecular docking experiments against EGFR and VEGFR-2 proteins. The metabolic extract of A. indica inhibited Streptococcus iniae and Staphylococcus sciuri with inhibition zones of 21.66 ± 0.57 mm and 20.33 ± 0.57 mm, respectively. The B. diffusa leaf extract produced inhibition zones of 20.3333 ± 0.5773 mm and 20.33 ± 0.57 mm against Streptococcus iniae and Edwardsiella anguillarum, respectively. A. indica and B. diffusa extracts had toxicities of 162.01 μg/ml and 175.6 μg/ml, respectively. Moreover, B. diffusa (IC50 =154.42 µg/ml) leaf extract exhibited moderately higher antioxidant activity compared with the A. indica (IC50 = 218.97 µg/ml) leaf extract. Multiple interactions were observed at Leu694, Met769 and Leu820 sites for EGFR and at Asp1046 and Cys1045 sites for VEGFR during the molecular docking study. CID-235030, CID-70825 and CID-156619353 had binding energies of -7.6 kJ/mol, -7.5 kJ/mol and -7.6 kJ/mol, respectively, with EGFR protein. VEGFR-2 protein had docking energies of -7.5 kJ/mol, -7.6 kJ/mol and -7.3 kJ/mol, respectively, for CID-6420353, CID-156619353 and CID-70825 compounds. The MD simulation trajectories revealed the hit compound; CID-235030 and EGFR complex, CID-6420353 and VEGFR-2 exhibit stable profile in the root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), hydrogen bond and root mean square fluctuation (RMSF) and the binding free energy by MM-PBSA method. This study indicates that methanol extracts of A. indica and B. diffusa may play a crucial role in developing antibiotic-resistant and cancer drugs.Communicated by Ramaswamy H. Sarma.

Acridones as promising drug candidates against Oropouche virus

Oropouche virus (OROV) is an emerging vector-borne arbovirus found in South America that causes Oropouche fever, a febrile infection similar to dengue fever. It has a high epidemic potential, causing illness in over 500,000 cases diagnosed since the virus was first discovered in 1955. Currently, the prevention of human viral infection depends on vaccination, but availability for many viruses is limited, and they are classified as neglected viruses. At present, there are no vaccines or antiviral treatments available. An alternative approach to limiting the spread of the virus is to selectively disrupt viral replication mechanisms. Here, we demonstrate the inhibitory effect of acridones, which efficiently inhibited viral replication by 99.9 % in vitro. To evaluate possible mechanisms of action, we conducted tests with dsRNA, an intermediate in virus replication, as well as MD simulations, docking, and binding free energy analysis. The results showed a strong interaction between FAC21 and the OROV endonuclease, which possibly limits the interaction of viral RNA with other proteins. Therefore, our results suggest a dual mechanism of antiviral action, possibly caused by ds-RNA intercalation. In summary, our findings demonstrate that a new generation of antiviral drugs could be developed based on the selective optimization of molecules.

The indole motif is essential for the antitrypanosomal activity of N5-substituted paullones

Severe infections with potentially fatal outcomes are caused by parasites from the genera Trypanosoma and Leishmania (class Kinetoplastea). The diseases affect people of remote areas in the tropics and subtropics with limited access to adequate health care. Besides insufficient diagnostics, treatment options are limited, with tenuous developments in recent years. Therefore, new antitrypanosomal antiinfectives are required to fight these maladies. In the presented approach, new compounds were developed and tested on the target trypanothione synthetase (TryS). This enzyme is crucial to the kinetoplastids' unique trypanothione-based thiol redox metabolism and thus for pathogen survival. Preceding studies have shown that N5-substituted paullones display antitrypanosomal activity as well as TryS inhibition. Herein, this compound class was further examined regarding the structure-activity relationships (SAR). Diverse benzazepinone derivatives were designed and tested in cell-based assays on bloodstream Trypanosoma brucei brucei (T. b. brucei) and intracellular amastigotes of Leishmania infantum (L. infantum) as well as in enzyme-based assays on L. infantum TryS (LiTryS) and T. b. brucei TryS (TbTryS). While an exchange of just the substituent in the 9-position of paullones led to potent inhibitors on LiTryS and T. b. brucei parasites, new compounds lacking the indole moiety showed a total loss of activity in both assays. Conclusively, the indole as part of the paullone structure is pivotal for keeping the TryS inhibitory and antitrypanosomal activity of this substance class.

Enhancement by pyrazolones of colistin efficacy against mcr-1-expressing E. coli: an in silico and in vitro investigation

Owing to the emergence of antibiotic resistance, the polymyxin colistin has been recently revived to treat acute, multidrug-resistant Gram-negative bacterial infections. Positively charged colistin binds to negatively charged lipids and damages the outer membrane of Gram-negative bacteria. However, the MCR-1 protein, encoded by the mobile colistin resistance (mcr) gene, is involved in bacterial colistin resistance by catalysing phosphoethanolamine (PEA) transfer onto lipid A, neutralising its negative charge, and thereby reducing its interaction with colistin. Our preliminary results showed that treatment with a reference pyrazolone compound significantly reduced colistin minimal inhibitory concentrations in Escherichia coli expressing mcr-1 mediated colistin resistance (Hanpaibool et al. in ACS Omega, 2023). A docking-MD combination was used in an ensemble-based docking approach to identify further pyrazolone compounds as candidate MCR-1 inhibitors. Docking simulations revealed that 13/28 of the pyrazolone compounds tested are predicted to have lower binding free energies than the reference compound. Four of these were chosen for in vitro testing, with the results demonstrating that all the compounds tested could lower colistin MICs in an E. coli strain carrying the mcr-1 gene. Docking of pyrazolones into the MCR-1 active site reveals residues that are implicated in ligand-protein interactions, particularly E246, T285, H395, H466, and H478, which are located in the MCR-1 active site and which participate in interactions with MCR-1 in ≥ 8/10 of the lowest energy complexes. This study establishes pyrazolone-induced colistin susceptibility in E. coli carrying the mcr-1 gene, providing a method for the development of novel treatments against colistin-resistant bacteria.

Computational Exploration of Potential CFTR Binding Sites for Type I Corrector Drugs

Cystic fibrosis (CF) is a recessive genetic disease that is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) protein. The recent development of a class of drugs called "correctors", which repair the structure and function of mutant CFTR, has greatly enhanced the life expectancy of CF patients. These correctors target the most common disease causing CFTR mutant F508del and are exemplified by the FDA-approved VX-809. While one binding site of VX-809 to CFTR was recently elucidated by cryo-electron microscopy, four additional binding sites have been proposed in the literature and it has been theorized that VX-809 and structurally similar correctors may engage multiple CFTR binding sites. To explore these five binding sites, ensemble docking was performed on wild-type CFTR and the F508del mutant using a large library of structurally similar corrector drugs, including VX-809 (lumacaftor), VX-661 (tezacaftor), ABBV-2222 (galicaftor), and a host of other structurally related molecules. For wild-type CFTR, we find that only one site, located in membrane spanning domain 1 (MSD1), binds favorably to our ligand library. While this MSD1 site also binds our ligand library for F508del-CFTR, the F508del mutation also opens a binding site in nucleotide binding domain 1 (NBD1), which enables strong binding of our ligand library to this site. This NBD1 site in F508del-CFTR exhibits the strongest overall binding affinity for our library of corrector drugs. This data may serve to better understand the structural changes induced by mutation of CFTR and how correctors bind to the protein. Additionally, it may aid in the design of new, more effective CFTR corrector drugs.

CHARMM-GUI-Based Induced Fit Docking Workflow to Generate Reliable Protein-Ligand Binding Modes

Molecular docking is a preferred method to predict ligand binding modes and their binding energy to target protein receptors, which is critical in early phase structure-based drug discovery. However, there is a persistent challenge in docking that can be attributed to the induced fit effect, as receptor binding sites undergo induced fit conformational changes upon ligand binding to achieve better binding modes. In this work, based on CHARMM-GUI LBS Finder& Refiner and High-Throughput Simulator, we present a straightforward CHARMM-GUI induced fit docking (CGUI-IFD) workflow to generate reliable protein-ligand binding modes. The CGUI-IFD workflow generates an ensemble of receptor binding site conformations through ligand-binding site (LBS) refinement, runs rigid receptor docking, and performs high-throughput molecular dynamics (MD) simulations of protein-ligand complex structures in explicit solvents. The results are evaluated based on the ligand root-mean-square deviation (RMSD)-based binding stability and the molecular mechanics generalized Born surface area binding energy. For a benchmark test, we used 258 cross-docking protein-ligand pairs across 41 target proteins from the Schrodinger IFD-MD data set. The application of CGUI-IFD on this data set shows 80% success rate (within 2.5 Å RMSD from the experimental structures). We expect that the CGUI-IFD workflow can be useful to generate reliable ligand binding modes for cross-docking cases.

Structure-Based Design of Novel MAO-B Inhibitors: A Review

With the significant growth of patients suffering from neurodegenerative diseases (NDs), novel classes of compounds targeting monoamine oxidase type B (MAO-B) are promptly emerging as distinguished structures for the treatment of the latter. As a promising function of computer-aided drug design (CADD), structure-based virtual screening (SBVS) is being heavily applied in processes of drug discovery and development. The utilization of molecular docking, as a helping tool for SBVS, is providing essential data about the poses and the occurring interactions between ligands and target molecules. The current work presents a brief discussion of the role of MAOs in the treatment of NDs, insight into the advantages and drawbacks of docking simulations and docking software, and a look into the active sites of MAO-A and MAO-B and their main characteristics. Thereafter, we report new chemical classes of MAO-B inhibitors and the essential fragments required for stable interactions focusing mainly on papers published in the last five years. The reviewed cases are separated into several chemically distinct groups. Moreover, a modest table for rapid revision of the revised works including the structures of the reported inhibitors together with the utilized docking software and the PDB codes of the crystal targets applied in each study is provided. Our work could be beneficial for further investigations in the search for novel, effective, and selective MAO-B inhibitors.

Antiviral Activity of Quercetin Hydrate against Zika Virus

Zika virus (ZIKV) has re-emerged in recent decades, leading to outbreaks of Zika fever in Africa, Asia, and Central and South America. Despite its drastic re-emergence and clinical impact, no vaccines or antiviral compounds are available to prevent or control ZIKV infection. This study evaluated the potential antiviral activity of quercetin hydrate against ZIKV infection and demonstrated that this substance inhibits virus particle production in A549 and Vero cells under different treatment conditions. In vitro antiviral activity was long-lasting (still observed 72 h post-infection), suggesting that quercetin hydrate affects multiple rounds of ZIKV replication. Molecular docking indicates that quercetin hydrate can efficiently interact with the specific allosteric binding site cavity of the NS2B-NS3 proteases and NS1-dimer. These results identify quercetin as a potential compound to combat ZIKV infection in vitro.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

In silico investigation of the interactions of certain drugs proposed for the treatment of Covid-19 with the paraoxonase-1

Coronavirus disease 2019 (Covid-19) has caused one of the biggest pandemics of modern times, infected over 240 million people and killed over 4.9 million people, and continues to do so. Although many drugs are widely recommended in the treatment of this disease, the interactions of these drugs with an anti-atherosclerotic enzyme, paraoxonase-1 (PON1), are not well known. In our study, we investigated the interactions of 18 different drugs, which are claimed to be effective against covid-19, with the PON1 enzyme and its genetics variants L55M and Q192R with molecular docking, molecular dynamics simulation and free energy calculation method MM/PBSA. We found that ruxolitinib, dexamethasone, colchicine; dexamethasone, sitagliptin, baricitinib and galidesivir, ruxolitinib, hydroxychloroquine were the most effective compounds in binding PON1-w, PON1L55M and PON1Q192R respectively. Mainly, sitagliptin, galidesivir and hydroxychloroquine have attracted attention by showing very high affinity (<-300 kJ/mol) according to the MM/PBSA method. We concluded that the drug interactions should be considered and more attention should be paid in the use of these drugs.Communicated by Ramaswamy H. Sarma.

An ensemble docking-based virtual screening according to different TRPV1 pore states toward identifying phytochemical activators

Identifying phytochemical activators for TRPV1 using ensemble-based virtual screening, machine learning, and MD simulation.

Targeting human thymidylate synthase: Ensemble-based virtual screening for drug repositioning and the role of water

A drug repositioning computational approach was carried to search inhibitors for human thymidylate synthase. An ensemble-based virtual screening of FDA-approved drugs showed the drugs Imatinib, Lumacaftor and Naldemedine to be likely candidates for repurposing. The role of water in the drug-receptor interactions was revealed by the application of an extended AutoDock scoring function that included the water forcefield. The binding affinity scores when hydrated ligands were docked were improved in the drugs considered. Further binding free energy calculations based on the Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that Imatinib, Lumacaftor and Naldemedine scored -130.7 ± 28.1, -210.6 ± 29.9 and -238.0 ± 25.4 kJ/mol, respectively, showing good binding affinity for the candidates considered. Overall, the analysis of the molecular dynamics trajectory of the receptor-drug complexes revealed stable structures for Imatinib, Lumacaftor and Naldemedine, for the entire simulation time.

In silico design of HDAC6 inhibitors with neuroprotective effects

HDAC6 has emerged as a molecular target to treat neurodegenerative disorders, due to its participation in protein aggregate degradation, oxidative stress process, mitochondrial transport, and axonal transport. Thus, in this work we have designed a set of 485 compounds with hydroxamic and bulky-hydrophobic moieties that may function as HDAC6 inhibitors with a neuroprotective effect. These compounds were filtered by their predicted ADMET properties and their affinity to HDAC6 demonstrated by molecular docking and molecular dynamics simulations. The combination of in silico with in vitro neuroprotective results allowed the identification of a lead compound (FH-27) which shows neuroprotective effect that could be due to HDAC6 inhibition. Further, FH-27 chemical moiety was used to design a second series of compounds improving the neuroprotective effect from 2- to 10-fold higher (YSL-99, YSL-109, YSL-112, YSL-116 and YSL-121; 1.25 ± 0.67, 1.82 ± 1.06, 7.52 ± 1.78, 5.59 and 5.62 ± 0.31 µM, respectively). In addition, the R enantiomer of FH-27 (YSL-106) was synthesized, showing a better neuroprotective effect (1.27 ± 0.60 µM). In conclusion, we accomplish the in silico design, synthesis, and biological evaluation of hydroxamic acid derivatives with neuroprotective effect as suggested by an in vitro model. Communicated by Ramaswamy H. Sarma.

Beyond sequence: Structure-based machine learning

Recent breakthroughs in protein structure prediction demarcate the start of a new era in structural bioinformatics. Combined with various advances in experimental structure determination and the uninterrupted pace at which new structures are published, this promises an age in which protein structure information is as prevalent and ubiquitous as sequence. Machine learning in protein bioinformatics has been dominated by sequence-based methods, but this is now changing to make use of the deluge of rich structural information as input. Machine learning methods making use of structures are scattered across literature and cover a number of different applications and scopes; while some try to address questions and tasks within a single protein family, others aim to capture characteristics across all available proteins. In this review, we look at the variety of structure-based machine learning approaches, how structures can be used as input, and typical applications of these approaches in protein biology. We also discuss current challenges and opportunities in this all-important and increasingly popular field.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Modeling receptor flexibility in the structure-based design of KRASG12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS^G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS^G12C—sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS^G12C. In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS^G12C, we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS^G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

Emerging frontiers in virtual drug discovery: From quantum mechanical methods to deep learning approaches

Virtual screening-based approaches to discover initial hit and lead compounds have the potential to reduce both the cost and time of early drug discovery stages, as well as to find inhibitors for even challenging target sites such as protein-protein interfaces. Here in this review, we provide an overview of the progress that has been made in virtual screening methodology and technology on multiple fronts in recent years. The advent of ultra-large virtual screens, in which hundreds of millions to billions of compounds are screened, has proven to be a powerful approach to discover highly potent hit compounds. However, these developments are just the tip of the iceberg, with new technologies and methods emerging to propel the field forward. Examples include novel machine-learning approaches, which can reduce the computational costs of virtual screening dramatically, while progress in quantum-mechanical approaches can increase the accuracy of predictions of various small molecule properties.

Extensive Sampling of Molecular Dynamics Simulations to Identify Reliable Protein Structures for Optimized Virtual Screening Studies: The Case of the hTRPM8 Channel

(1) Background: Virtual screening campaigns require target structures in which the pockets are properly arranged for binding. Without these, MD simulations can be used to relax the available target structures, optimizing the fine architecture of their binding sites. Among the generated frames, the best structures can be selected based on available experimental data. Without experimental templates, the MD trajectories can be filtered by energy-based criteria or sampled by systematic analyses. (2) Methods: A blind and methodical analysis was performed on the already reported MD run of the hTRPM8 tetrameric structures; a total of 50 frames underwent docking simulations by using a set of 1000 ligands including 20 known hTRPM8 modulators. Docking runs were performed by LiGen program and involved the frames as they are and after optimization by SCRWL4.0. For each frame, all four monomers were considered. Predictive models were developed by the EFO algorithm based on the sole primary LiGen scores. (3) Results: On average, the MD simulation progressively enhances the performance of the extracted frames, and the optimized structures perform better than the non-optimized frames (EF1% mean: 21.38 vs. 23.29). There is an overall correlation between performances and volumes of the explored pockets and the combination of the best performing frames allows to develop highly performing consensus models (EF1% = 49.83). (4) Conclusions: The systematic sampling of the entire MD run provides performances roughly comparable with those previously reached by using rationally selected frames. The proposed strategy appears to be helpful when the lack of experimental data does not allow an easy selection of the optimal structures for docking simulations. Overall, the reported docking results confirm the relevance of simulating all the monomers of an oligomer structure and emphasize the efficacy of the SCRWL4.0 method to optimize the protein structures for docking calculations.

Dynamic Docking of Macrocycles in Bound and Unbound Protein Structures with DynaDock

Macrocycles are interesting molecules with unique features due to their conformationally constrained yet flexible ring structure. This characteristic poses a difficult challenge for computational modeling studies since they rely on accurate structural descriptions. In particular, molecular docking calculations suffer from the lack of ring flexibility during pose generation, which is often compensated by using pregenerated ligand conformer ensembles. Moreover, receptor structures are mainly treated rigidly, which limits the use of many docking tools. In this study, we optimized our previous molecular dynamics-based sampling and docking pipeline specifically designed for the accurate prediction of macrocyclic compounds. We developed a dihedral classification procedure for in-depth conformational analysis of the macrocyclic rings and extracted structural ensembles that were subsequently docked in both bound and unbound protein structures employing a fully flexible approach. Our results suggest that including a ring conformer close to the bound state in the starting ensemble increases the chance of successful docking. The bioactive conformations of a diverse set of ligands could be predicted with high and decent accuracy in bound and unbound protein structures, respectively, due to the incorporation of full molecular flexibility in our approach. The remaining unsuccessful docking calculations were mainly caused by large flexible substituents that bind to surface-exposed binding sites, rather than the macrocyclic ring per se and could be further improved by explicit molecular dynamics simulations of the docked complex.

Functional inhibition of c-Myc using novel inhibitors identified through “hot spot” targeting

Protein–protein interactions drive various biological processes in healthy as well as disease states. The transcription factor c-Myc plays a crucial role in maintaining cellular homeostasis, and its deregulated expression is linked to various human cancers; therefore, it can be considered a viable target for cancer therapeutics. However, the structural heterogeneity of c-Myc due to its disordered nature poses a major challenge to drug discovery. In the present study, we used an in silico alanine scanning mutagenesis approach to identify “hot spot” residues within the c-Myc/Myc-associated factor X interface, which is highly disordered and has not yet been systematically analyzed for potential small molecule binding sites. We then used the information gained from this analysis to screen potential inhibitors using a conformation ensemble approach. The fluorescence-based biophysical experiments showed that the identified hit molecules displayed noncovalent interactions with these hot spot residues, and further cell-based experiments showed substantial in vitro potency against diverse c-Myc-expressing cancer/stem cells by deregulating c-Myc activity. These biophysical and computational studies demonstrated stable binding of the hit compounds with the disordered c-Myc protein. Collectively, our data indicated effective drug targeting of the disordered c-Myc protein via the determination of hot spot residues in the c-Myc/Myc-associated factor X heterodimer.

From Data to Knowledge: Systematic Review of Tools for Automatic Analysis of Molecular Dynamics Output

An increasing number of crystal structures available on one side, and the boost of computational power available for computer-aided drug design tasks on the other, have caused that the structure-based drug design tools are intensively used in the drug development pipelines. Docking and molecular dynamics simulations, key representatives of the structure-based approaches, provide detailed information about the potential interaction of a ligand with a target receptor. However, at the same time, they require a three-dimensional structure of a protein and a relatively high amount of computational resources. Nowadays, as both docking and molecular dynamics are much more extensively used, the amount of data output from these procedures is also growing. Therefore, there are also more and more approaches that facilitate the analysis and interpretation of the results of structure-based tools. In this review, we will comprehensively summarize approaches for handling molecular dynamics simulations output. It will cover both statistical and machine-learning-based tools, as well as various forms of depiction of molecular dynamics output.

Benchmarking ensemble docking methods in D3R Grand Challenge 4

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.

Ensemble learning from ensemble docking: revisiting the optimum ensemble size problem

Despite considerable advances obtained by applying machine learning approaches in protein–ligand affinity predictions, the incorporation of receptor flexibility has remained an important bottleneck. While ensemble docking has been used widely as a solution to this problem, the optimum choice of receptor conformations is still an open question considering the issues related to the computational cost and false positive pose predictions. Here, a combination of ensemble learning and ensemble docking is suggested to rank different conformations of the target protein in light of their importance for the final accuracy of the model. Available X-ray structures of cyclin-dependent kinase 2 (CDK2) in complex with different ligands are used as an initial receptor ensemble, and its redundancy is removed through a graph-based redundancy removal, which is shown to be more efficient and less subjective than clustering-based representative selection methods. A set of ligands with available experimental affinity are docked to this nonredundant receptor ensemble, and the energetic features of the best scored poses are used in an ensemble learning procedure based on the random forest method. The importance of receptors is obtained through feature selection measures, and it is shown that a few of the most important conformations are sufficient to reach 1 kcal/mol accuracy in affinity prediction with considerable improvement of the early enrichment power of the models compared to the different ensemble docking without learning strategies. A clear strategy has been provided in which machine learning selects the most important experimental conformers of the receptor among a large set of protein–ligand complexes while simultaneously maintaining the final accuracy of affinity predictions at the highest level possible for available data. Our results could be informative for future attempts to design receptor-specific docking-rescoring strategies.

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via ¹³C-metabolic flux analysis (¹³C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated ¹³C labeling data were used with ¹³C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

DINC-COVID: A webserver for ensemble docking with flexible SARS-CoV-2 proteins

An unprecedented research effort has been undertaken in response to the ongoing COVID-19 pandemic. This has included the determination of hundreds of crystallographic structures of SARS-CoV-2 proteins, and numerous virtual screening projects searching large compound libraries for potential drug inhibitors. Unfortunately, these initiatives have had very limited success in producing effective inhibitors against SARS-CoV-2 proteins. A reason might be an often overlooked factor in these computational efforts: receptor flexibility. To address this issue we have implemented a computational tool for ensemble docking with SARS-CoV-2 proteins. We have extracted representative ensembles of protein conformations from the Protein Data Bank and from in silico molecular dynamics simulations. Twelve pre-computed ensembles of SARS-CoV-2 protein conformations have now been made available for ensemble docking via a user-friendly webserver called DINC-COVID (dinc-covid.kavrakilab.org). We have validated DINC-COVID using data on tested inhibitors of two SARS-CoV-2 proteins, obtaining good correlations between docking-derived binding energies and experimentally-determined binding affinities. Some of the best results have been obtained on a dataset of large ligands resolved via room temperature crystallography, and therefore capturing alternative receptor conformations. In addition, we have shown that the ensembles available in DINC-COVID capture different ranges of receptor flexibility, and that this diversity is useful in finding alternative binding modes of ligands. Overall, our work highlights the importance of accounting for receptor flexibility in docking studies, and provides a platform for the identification of new inhibitors against SARS-CoV-2 proteins.

Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

Identification of ligand binding sites in intrinsically disordered proteins with a differential binding score

Screening ligands directly binding to an ensemble of intrinsically disordered proteins (IDP) to discover potential hits or leads for new drugs is an emerging but challenging area as IDPs lack well-defined and ordered 3D-protein structures. To explore a new IDP-based rational drug discovery strategy, a differential binding score (DIBS) is defined. The basis of DIBS is to quantitatively determine the binding preference of a ligand to an ensemble of conformations specified by IDP versus such preferences to an ensemble of random coil conformations of the same protein. Ensemble docking procedures performed on repeated sampling of conformations, and the results tested for statistical significance determine the preferential ligand binding sites of the IDP. The results of this approach closely reproduce the experimental data from recent literature on the binding of the ligand epigallocatechin gallate (EGCG) to the intrinsically disordered N-terminal domain of the tumor suppressor p53. Combining established approaches in developing a new method to screen ligands against IDPs could be valuable as a screening tool for IDP-based drug discovery.

In silico investigation to identify potential small molecule inhibitors of the RNA-dependent RNA polymerase (RdRp) nidovirus RdRp-associated nucleotidyltransferase domain

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp) is a promising target for antiviral drugs. In this study, a chemical library (n = 300) was screened against the nidovirus RdRp-associated nucleotidyltransferase (NiRAN) domain. Blind docking was performed using a selection of 30 compounds and nine ligands were chosen based on their docking scores, safety profile, and availability. Using cluster analysis on a 10 microsecond molecular dynamics simulation trajectory (from D.E. Shaw Research), the compounds were docked to the different conformations. On the basis of our modelling studies, oleuropein was identified as a potential lead compound.

A computational study of cooperative binding to multiple SARS-CoV-2 proteins

Structure-based drug design targeting the SARS-CoV-2 virus has been greatly facilitated by available virus-related protein structures. However, there is an urgent need for effective, safe small-molecule drugs to control the spread of the virus and variants. While many efforts are devoted to searching for compounds that selectively target individual proteins, we investigated the potential interactions between eight proteins related to SARS-CoV-2 and more than 600 compounds from a traditional Chinese medicine which has proven effective at treating the viral infection. Our original ensemble docking and cooperative docking approaches, followed by a total of over 16-micorsecond molecular simulations, have identified at least 9 compounds that may generally bind to key SARS-CoV-2 proteins. Further, we found evidence that some of these compounds can simultaneously bind to the same target, potentially leading to cooperative inhibition to SARS-CoV-2 proteins like the Spike protein and the RNA-dependent RNA polymerase. These results not only present a useful computational methodology to systematically assess the anti-viral potential of small molecules, but also point out a new avenue to seek cooperative compounds toward cocktail therapeutics to target more SARS-CoV-2-related proteins.

Exploring the stability of inhibitor binding to SIK2 using molecular dynamics simulation and binding free energy calculation

Salt inducible kinase 2 (SIK2) is a calcium/calmodulin-dependent protein kinase-like kinase that is implicated in a variety of biological phenomena, including cellular metabolism, growth, and apoptosis. SIK2 is the key target for various cancers, including ovarian, breast, prostate, and lung cancers. Although potent inhibitors of SIK2 are being developed, their binding stability and functional role are not presently known. In this work, we studied the detailed interactions between SIK2 and four of its inhibitors, HG-9-91-01, KIN112, MRT67307, and MRT199665, using molecular docking, molecular dynamics simulation, binding free energy calculation, and interaction fingerprint analysis. Intermolecular interactions revealed that HG-9-91-01 and KIN112 have stronger interactions with SIK2 than those of MRT199665 and MRT67307. The key residues involved in binding with SIK2 are conserved among all four inhibitors. Our results explain the detailed interaction of SIK2 with its inhibitors at the molecular level, thus paving the way for the development of targeted efficient anti-cancer drugs.

A Blueprint for High Affinity SARS-CoV-2 Mpro Inhibitors from Activity-Based Compound Library Screening Guided by Analysis of Protein Dynamics

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein's druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.

Computational Modeling as a Tool to Investigate PPI: From Drug Design to Tissue Engineering

Protein-protein interactions (PPIs) mediate a large number of important regulatory pathways. Their modulation represents an important strategy for discovering novel therapeutic agents. However, the features of PPI binding surfaces make the use of structure-based drug discovery methods very challenging. Among the diverse approaches used in the literature to tackle the problem, linear peptides have demonstrated to be a suitable methodology to discover PPI disruptors. Unfortunately, the poor pharmacokinetic properties of linear peptides prevent their direct use as drugs. However, they can be used as models to design enzyme resistant analogs including, cyclic peptides, peptide surrogates or peptidomimetics. Small molecules have a narrower set of targets they can bind to, but the screening technology based on virtual docking is robust and well tested, adding to the computational tools used to disrupt PPI. We review computational approaches used to understand and modulate PPI and highlight applications in a few case studies involved in physiological processes such as cell growth, apoptosis and intercellular communication.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Microsecond MD Simulation and Multiple-Conformation Virtual Screening to Identify Potential Anti-COVID-19 Inhibitors Against SARS-CoV-2 Main Protease

The recent pandemic outbreak of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), raised global health and economic concerns. Phylogenetically, SARS-CoV-2 is closely related to SARS-CoV, and both encode the enzyme main protease (M^pro/3CL^pro), which can be a potential target inhibiting viral replication. Through this work, we have compiled the structural aspects of M^pro conformational changes, with molecular modeling and 1-μs MD simulations. Long-scale MD simulation resolves the mechanism role of crucial amino acids involved in protein stability, followed by ensemble docking which provides potential compounds from the Traditional Chinese Medicine (TCM) database. These lead compounds directly interact with active site residues (His41, Gly143, and Cys145) of M^pro, which plays a crucial role in the enzymatic activity. Through the binding mode analysis in the S1, S1′, S2, and S4 binding subsites, screened compounds may be functional for the distortion of the oxyanion hole in the reaction mechanism, and it may lead to the inhibition of M^pro in SARS-CoV-2. The hit compounds are naturally occurring compounds; they provide a sustainable and readily available option for medical treatment in humans infected by SARS-CoV-2. Henceforth, extensive analysis through molecular modeling approaches explained that the proposed molecules might be promising SARS-CoV-2 inhibitors for the inhibition of COVID-19, subjected to experimental validation.

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

Computational identification of potential chemoprophylactic agents according to dynamic behavior of peroxisome proliferator-activated receptor gamma

Peroxisome proliferator-activated receptor gamma (PPARγ) is an attractive target for chemoprevention of lung carcinoma, however its highly dynamic nature has plagued drug development for decades, with difficulties in receptor modeling for structure-based design. In this work, an integrated receptor-based virtual screening (VS) strategy was applied to identify PPARγ agonists as chemoprophylactic agents by using extensive docking and conformational sampling methods. Our results showed that the conformational plasticity of PPARγ, especially the H2 & S245 loop, H2' & Ω loop and AF-2 surface, is markedly affected by binding of full/partial agonists. To fully take the dynamic behavior of PPARγ into account, the VS approach effectively sorts out five commercial agents with reported antineoplastic properties. Among them, ZINC03775146 (gusperimus) and ZINC14087743 (miltefosine) might be novel PPARγ agonists with the potential for chemoprophylaxis, that simultaneously take part in a flexible switch of the AF-2 surface and state change of the Ω loop. Furthermore, the dynamic structural coupling between the H2 & S245 and H2' & Ω loops offers enticing hope for PPARγ-targeted therapeutics, by blocking kinase accessibility to PPARγ. These results might aid the development of chemopreventive drugs, and the integrated VS strategy could be conducive to drug design for highly flexible biomacromolecules.

CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations

Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.