Accuracy evaluation and addition of improved dihedral parameters for the MMFF94s

PheSA: An Open-Source Tool for Pharmacophore-Enhanced Shape Alignment

PheSA is an open-source pharmacophore- and shape-based screening and molecular alignment tool that is fully open-source as part of OpenChemLib. Supporting standard ligand-based screening, flexible refinement of alignments, and receptor-guided shape docking, PheSA is a very flexible tool and can be used for different use cases in structure-based drug design. We present the algorithm and different benchmark studies that investigate the screening performance and also the quality of the generated alignments and the pose prediction performance of the receptor-guided PheSA algorithm. An important finding is the effect of the type of similarity metric used for measuring screening enrichment (symmetric Tanimoto versus asymmetric Tversky), whereby we could observe improved enrichment rates by using Tversky. PheSA exhibits enrichments on the DUD-E that are on par with commercial methods.

Ensemble docking based virtual screening of SARS-CoV-2 main protease inhibitors

During the first years of COVID-19 pandemic, X-ray structures of the coronavirus drug targets were acquired at an unprecedented rate, giving hundreds of PDB depositions in less than a year. The main protease (Mpro) of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) is the primary validated target of direct-acting antivirals. The selection of the optimal ensemble of structures of Mpro for the docking-driven virtual screening campaign was thus non-trivial and required a systematic and automated approach. Here we report a semi-automated active site RMSD based procedure of ensemble selection from the SARS-CoV-2 Mpro crystallographic data and virtual screening of its inhibitors. The procedure was compared with other approaches to ensemble selection and validated with the help of hand-picked and peer-reviewed activity-annotated libraries. Prospective virtual screening of non-covalent Mpro inhibitors resulted in a new chemotype of thienopyrimidinone derivatives with experimentally confirmed enzyme inhibition.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Synthesis and in silico study of chenodeoxycholic acid and its analogues as an alternative inhibitor of spike glycoprotein of SARS-CoV-2

COVID-19, caused by SARS-CoV-2, is a viral infection that has generated one of the most significant health problems in the world. Spike glycoprotein is a crucial enzyme in viral replication and transcription mediation. There are reports in the literature on using bile acid in the fight against this virus through in vitro tests. This work presents the synthesis of nine chenodeoxycholic acid derivatives (1-9), which were prepared by oxidation, acetylation, formylation, and esterification reactions, and the analogs 6-9 have not yet been reported in the literature and the possibility of conducting an in silico study of bile acid derivatives as a therapeutic alternative to combat the virus using glycoprotein as a macromolecular target. As a result, five compounds (1, 6-9) possessed favorable competitive interactions with the lowest energies compared to the native ligand (BLA), and the highlighted compound 9 got the best scores. At the same time, analog 1 presented the best ADME filter result. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. Physical and spectroscopic data have fully characterized all the compounds.Communicated by Ramaswamy H. Sarma.

High-Quality Conformer Generation with CONFORGE: Algorithm and Performance Assessment

Knowledge of the putative bound-state conformation of a molecule is an essential prerequisite for the successful application of many computer-aided drug design methods that aim to assess or predict its capability to bind to a particular target receptor. An established approach to predict bioactive conformers in the absence of receptor structure information is to sample the low-energy conformational space of the investigated molecules and derive representative conformer ensembles that can be expected to comprise members closely resembling possible bound-state ligand conformations. The high relevance of such conformer generation functionality led to the development of a wide panel of dedicated commercial and open-source software tools throughout the last decades. Several published benchmarking studies have shown that open-source tools usually lag behind their commercial competitors in many key aspects. In this work, we introduce the open-source conformer ensemble generator CONFORGE, which aims at delivering state-of-the-art performance for all types of organic molecules in drug-like chemical space. The ability of CONFORGE and several well-known commercial and open-source conformer ensemble generators to reproduce experimental 3D structures as well as their computational efficiency and robustness has been assessed thoroughly for both typical drug-like molecules and macrocyclic structures. For small molecules, CONFORGE clearly outperformed all other tested open-source conformer generators and performed at least equally well as the evaluated commercial generators in terms of both processing speed and accuracy. In the case of macrocyclic structures, CONFORGE achieved the best average accuracy among all benchmarked generators, with RDKit's generator coming close in second place.

De novo design of bioactive phenol and chromone derivatives for inhibitors of Spike glycoprotein of SARS-CoV-2 in silico

This work presents the synthesis of 12 phenol and chromone derivatives, prepared by the analogs, and the possibility of conducting an in silico study of its derivatives as a therapeutic alternative to combat the SARS-CoV-2, pathogen responsible for COVID-19 pandemic, using its S-glycoprotein as a macromolecular target. After the initial screening for the ranking of the products, it was chosen which structure presented the best energy bond with the target. As a result, derivative 4 was submitted to a molecular growth study using artificial intelligence, where 8436 initial structures were obtained that passed through the interaction filters and similarity to the active glycoprotein pocket through the MolAICal computational package. Thus, 557 Hits with active configuration were generated, which is very promising compared to the BLA reference link for inhibiting the biological target. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. The Hit 48 and 250 are the most active compounds against SARS-CoV-2. In summary, the results show that the Hit 250 would be more active than the natural compound, which could be further developed for further testing against SARS-CoV-2. The study employs the de novo approach to design new drugs, combining artificial intelligence and molecular dynamics simulations to create efficient molecular structures. This research aims to contribute to the development of effective therapeutic strategies against the pandemic.

Phomactinine, the First Nitrogen-Bearing Phomactin, Produced by Biatriospora sp. CBMAI 1333

Metabolomics analyses and improvement of growth conditions were applied toward diversification of phomactin terpenoids by the fungus Biatriospora sp. CBMAI 1333. Visualization of molecular networking results on Gephi assisted the observation of phomactin diversification and guided the isolation of new phomactin variants by applying a modified version of chemometrics based on a fractional factorial design. Consequentially, the first nitrogen-bearing phomactin, phomactinine (1), with a new rearranged carbon skeleton, was isolated and identified. The strategy combining metabolomics and chemometrics can be extended to include bioassay potency, structure novelty, and metabolic diversification connected or not to genomic analyses.

Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R²) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.

Semisynthesis, in silico study and in vitro antibacterial evaluation of fucosterol derivatives

Thirteen fucosterol derivatives were prepared by structural modification at the hydroxyl group in C-3 and catalytic hydrogenation at the carbon-carbon double bond in C-5(6) and C-24(28). The structures of all compounds were established based on their spectral data (IR, MS, and NMR). Fucosterol (1) and its derivatives (2-12, and a mixture of 13a and 13b) were evaluated for their in vitro antibacterial activity against Klebsiella pneumoniae (ATCC 10031), Escherichia coli (ATCC 10536), Pseudomonas aeruginosa (ATCC 15442), Streptococcus mutans (ATCC 0046) and Staphylococcus aureus using the microdilution method. Among them, 1, 8, 9, 10, and a mixture of 13a and 13b exhibited the best antibacterial activity. The derivative 7 was inactive against all bacterial strains evaluated (MIC ≥ 2.327 mM). In addition, the investigation of binding interactions of more active compounds (1, 8, 9, 10, and mixture of 13a and 13b) to appropriate proteins was performed using molecular docking. This paper registers for the first time the in silico studies on the antibacterial activity of compounds 1, 8, 9, 10, and mixture of 13a/13b, and the spectral data of compounds 4, 6, and 7.

Discovery and biosynthesis of macrophasetins from the plant pathogen fungus Macrophomina phaseolina

3-Decalinoyltetramic acids (DTAs) are a class of natural products with chemical diversity and potent bioactivities. In fungal species there is a general biosynthetic route to synthesize this type of compounds, which usually features a polyketide synthase-nonribosomal peptide synthetase (PKS-NRPS) and a lipocalin-like Diels-Alderase (LLDAse). Using a synthetic biology approach, combining the bioinformatics analysis prediction and heterologous expression, we mined a PKS-NRPS and LLDAse encoding gene cluster from the plant pathogenic fungus Macrophomina phaseolina and characterized the cluster to be responsible for the biosynthesis of novel DTAs, macrophasetins. In addition, we investigated the biosynthesis of these compounds and validated the accuracy of the phylogeny-guided bioinformatics analysis prediction. Our results provided a proof of concept example to this approach, which may facilitate the discovery of novel DTAs from the fungal kingdom.

Comparing the Performances of Force Fields in Conformational Searching of Hydrogen-Bond-Donating Catalysts

Here, we compare the relative performances of different force fields for conformational searching of hydrogen-bond-donating catalyst-like molecules. We assess the force fields by their predictions of conformer energies, geometries, low-energy, nonredundant conformers, and the maximum numbers of possible conformers. Overall, MM3, MMFFs, and OPLS3e had consistently strong performances and are recommended for conformationally searching molecules structurally similar to those in this study.

Targeting Epigenetic Regulators Using Machine Learning: Potential Sirtuin 2 Inhibitors

Sirtuin 2 (SIRT2) is a nicotinamide adenine dinucleotide (NAD+)-dependent deacetylase that has been identified as a target for many diseases, including Parkinson’s disease (PD) and leukemia. Using 234 SIRT2 inhibitors from the ZINC15 database, we generated molecular descriptors with PaDEL and constructed a machine-learning (ML) model for the binary classification of SIRT2 inhibitors. To predict compounds with novel inhibitory mechanisms, we then applied the model on the ZINC15/FDA subset, yielding 107 potential SIRT2 inhibitors. For validation of these substances, we employed the binding analysis software AutoDock Vina to perform virtual screening, with which 43 compounds were considered best inhibitors at the [Formula: see text][Formula: see text]kcal/mol binding affinity threshold. Our results demonstrate the potential of ligand-based (LB) ML techniques in conjunction with receptor-based virtual screening (RBVS) to facilitate the drug discovery or repurposing.

Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective

Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.

Structural bioinformatics used to predict the protein targets of remdesivir and flavones in SARS-CoV-2 infection

BACKGROUND

During the current SARS-CoV-2 pandemic, the identification of effective antiviral drugs is crucial. Unfortunately, no specific treatment or vaccine is available to date.

OBJECTIVE

Here, we aimed to predict the interactions between SARS-CoV-2 proteins and protein targets from the human body for some flavone molecules (kaempferol, morin, pectolinarin, myricitrin, and herbacetin) in comparison to synthetic compounds (hydroxychloroquine, remdesivir, ribavirin, ritonavir, AMD-070, favipiravir).

METHODS

Using MOE software and advanced bioinformatics and cheminformatics portals, we conducted an extensive analysis based on various structural and functional features of compounds, such as their amphiphilic field, flexibility, and steric features. The structural similarity analysis of natural and synthetic compounds was performed using Tanimoto coefficients. The interactions of some compounds with SARS-CoV-2 3CLprotease or RNA-dependent RNA polymerase were described using 2D protein-ligand interaction diagrams based on known crystal structures. The potential targets of considered compounds were identified using the SwissTargetPrediction web tool.

RESULTS

Our results showed that remdesivir, pectolinarin, and ritonavir present a strong structural similarity which may be correlated to their similar biological activity. As common molecular targets of compounds in the human body, ritonavir, kaempferol, morin, and herbacetin can activate multidrug resistance-associated proteins, while remdesivir, ribavirin, and pectolinarin appear as ligands for adenosine receptors.

CONCLUSION

Our evaluation recommends remdesivir, pectolinarin, and ritonavir as promising anti-SARS-CoV-2 agents.

Advanced Methods for Studying Structure and Interactions of Macrolide Antibiotics

Macrolide antibiotics are macrocyclic compounds that are clinically used and prescribed for the treatment of upper and lower respiratory tract infections. They inhibit the synthesis of bacterial proteins by reversible binding to the 23S rRNA at or near the peptidyl transferase center. However, their excellent antibacterial profile was largely compromised by the emergence of bacterial resistance. Today, fighting resistance to antibiotics is one of the greatest challenges in medicinal chemistry. Considering various physicochemical properties of macrolides, understanding their structure and interactions with macromolecular targets is crucial for the design of new antibiotics efficient against resistant pathogens. The solid-state structures of some macrolide-ribosome complexes have recently been solved, throwing new light on the macrolide binding mechanisms. On the other hand, a combination of NMR spectroscopy and molecular modeling calculations can be applied to study free and bound conformations in solution. In this article, a description of advanced physicochemical methods for elucidating the structure and interactions of macrolide antibiotics in solid state and solution will be provided, and their principal advantages and drawbacks will be discussed.

Quantitative structure-property relationship of the photoelectrochemical oxidation of phenolic pollutants at modified nanoporous titanium oxide using supervised machine learning

Here we report on an advanced photoelectrochemical (PEC) oxidation of 22 phenolic pollutants based on modified nanoporous TiO2, which was directly grown on a titanium substrate electrochemically. Their degradation rate constants were experimentally determined and their physicochemical properties were computaionally calculated. The quantitative structure-property relationship (QSPR) was elucidated by employing multiple linear regression (MLR) method. A supervised machine learning approach was employed to build QSPR models. The high predictive abilities of the QSPR model were validated via leave-one-out (LOO) method and a strict regimen of statistical validation tests. The significant descriptors identified in the QSPR Model for the phenolic compounds were also assessed using a typical dye pollutant Rhodamine B, further confirming the high effectiveness and predictability of the optimized model. Our study has shown that the integrated effect of the structural, hydrophobic and topological properties along with electronic property should be considered in order to design an efficient PEC catalytic approach for environmental applications.