Docking rigid macrocycles using Convex-PL, AutoDock Vina, and RDKit in the D3R Grand Challenge 4

SAnDReS 2.0: Development of machine-learning models to explore the scoring function space

Classical scoring functions may exhibit low accuracy in determining ligand binding affinity for proteins. The availability of both protein-ligand structures and affinity data make it possible to develop machine-learning models focused on specific protein systems with superior predictive performance. Here, we report a new methodology named SAnDReS that combines AutoDock Vina 1.2 with 54 regression methods available in Scikit-Learn to calculate binding affinity based on protein-ligand structures. This approach allows exploration of the scoring function space. SAnDReS generates machine-learning models based on crystal, docked, and AlphaFold-generated structures. As a proof of concept, we examine the performance of SAnDReS-generated models in three case studies. For all three cases, our models outperformed classical scoring functions. Also, SAnDReS-generated models showed predictive performance close to or better than other machine-learning models such as KDEEP, CSM-lig, and ΔVinaRF20. SAnDReS 2.0 is available to download at https://github.com/azevedolab/sandres.

Predicting blood-brain barrier permeability of molecules with a large language model and machine learning

Predicting the blood-brain barrier (BBB) permeability of small-molecule compounds using a novel artificial intelligence platform is necessary for drug discovery. Machine learning and a large language model on artificial intelligence (AI) tools improve the accuracy and shorten the time for new drug development. The primary goal of this research is to develop artificial intelligence (AI) computing models and novel deep learning architectures capable of predicting whether molecules can permeate the human blood-brain barrier (BBB). The in silico (computational) and in vitro (experimental) results were validated by the Natural Products Research Laboratories (NPRL) at China Medical University Hospital (CMUH). The transformer-based MegaMolBART was used as the simplified molecular input line entry system (SMILES) encoder with an XGBoost classifier as an in silico method to check if a molecule could cross through the BBB. We used Morgan or Circular fingerprints to apply the Morgan algorithm to a set of atomic invariants as a baseline encoder also with an XGBoost classifier to compare the results. BBB permeability was assessed in vitro using three-dimensional (3D) human BBB spheroids (human brain microvascular endothelial cells, brain vascular pericytes, and astrocytes). Using multiple BBB databases, the results of the final in silico transformer and XGBoost model achieved an area under the receiver operating characteristic curve of 0.88 on the held-out test dataset. Temozolomide (TMZ) and 21 randomly selected BBB permeable compounds (Pred scores = 1, indicating BBB-permeable) from the NPRL penetrated human BBB spheroid cells. No evidence suggests that ferulic acid or five BBB-impermeable compounds (Pred scores < 1.29423E-05, which designate compounds that pass through the human BBB) can pass through the spheroid cells of the BBB. Our validation of in vitro experiments indicated that the in silico prediction of small-molecule permeation in the BBB model is accurate. Transformer-based models like MegaMolBART, leveraging the SMILES representations of molecules, show great promise for applications in new drug discovery. These models have the potential to accelerate the development of novel targeted treatments for disorders of the central nervous system.

FaissMolLib: An efficient and easy deployable tool for ligand-based virtual screening

Virtual screening-based molecular similarity and fingerprint are crucial in drug design, target prediction, and ADMET prediction, aiding in identifying potential hits and optimizing lead compounds. However, challenges such as lack of comprehensive open-source molecular fingerprint databases and efficient search methods for virtual screening are prevalent. To address these issues, we introduce FaissMolLib, an open-source virtual screening tool that integrates 2.8 million compounds from ChEMBL and ZINC databases. Notably, FaissMolLib employs the highly efficient Faiss search algorithm, outperforming the Tanimoto algorithm in identifying similar molecules with its tighter clustering in scatter plots and lower mean, standard deviation, and variance in key molecular properties. This feature enables FaissMolLib to screen 2.8 million compounds in just 0.05 seconds, offering researchers an efficient, easily deployable solution for virtual screening on laptops and building unique compound databases. This significant advancement holds great potential for accelerating drug discovery efforts and enhancing chemical data analysis. FaissMolLib is freely available at http://liuhaihan.gnway.cc:80. The code and dataset of FaissMolLib are freely available at https://github.com/Superhaihan/FiassMolLib.

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.

Excavating phytochemicals from plants possessing antiviral activities for identifying SARS-CoV hemagglutinin-esterase inhibitors by diligent computational workflow

Coronaviruses (CoVs) belong to a group of RNA viruses that cause diseases in vertebrates including. Newer and deadlier than SARS CoV-2 are sought to appear in future for which the scientific community must be prepared with the strategies for their control. Spike protein (S-protein) of all the CoVs require angiotensin-converting enzyme2 (ACE2), while CoVs also require hemagglutinin-acetylesterase (HE) glycoprotein receptor to simultaneously interact with O-acetylated sialic acids on host cells, both these interactions enable viral particle to enter host cell leading to its infection. Target inhibition of viral S-protein and HE glycoprotein receptor can lead to a development of therapy against the SARS CoV-2. The proposition is to recognize molecules from the bundle of phytochemicals of medicinal plants known to possess antiviral potentials as a lead that could interact and mask the active site of, HE glycoprotein which would ideally bind to O-acetylated sialic acids on human host cells. Such molecules can be addressed as 'HE glycoprotein blockers'. A library of 110 phytochemicals from Withania somnifera, Asparagus racemosus, Zinziber officinalis, Allium sativum, Curcuma longa and Adhatoda vasica was constructed and was used under present study. In silico analysis was employed with plant-derived phytochemicals. The molecular docking, molecular dynamics simulations over the scale of 1000 ns (1 μs) and ADMET prediction revealed that the Withania somnifera (ashwagandha) and Asparagus racemosus (shatavari) plants possessed various steroidal saponins and alkaloids which could potentially inhibit the COVID-19 virus and even other CoVs targeted HE glycoprotein receptor.

Hydroxylation of Antitubercular Drug Candidate, SQ109, by Mycobacterial Cytochrome P450

Spreading of the multidrug-resistant (MDR) strains of the one of the most harmful pathogen Mycobacterium tuberculosis (Mtb) generates the need for new effective drugs. SQ109 showed activity against resistant Mtb and already advanced to Phase II/III clinical trials. Fast SQ109 degradation is attributed to the human liver Cytochrome P450s (CYPs). However, no information is available about interactions of the drug with Mtb CYPs. Here, we show that Mtb CYP124, previously assigned as a methyl-branched lipid monooxygenase, binds and hydroxylates SQ109 in vitro. A 1.25 Å-resolution crystal structure of the CYP124-SQ109 complex unambiguously shows two conformations of the drug, both positioned for hydroxylation of the ω-methyl group in the trans position. The hydroxylated SQ109 presumably forms stabilizing H-bonds with its target, Mycobacterial membrane protein Large 3 (MmpL3). We anticipate that Mtb CYPs could function as analogs of drug-metabolizing human CYPs affecting pharmacokinetics and pharmacodynamics of antitubercular (anti-TB) drugs.