Marine Drug Discovery through Computer-Aided Approaches.

Besides the importance of our oceans as oxygen factories, food providers, shipping pathways, and tourism enablers, oceans hide an unprecedented wealth of opportunities [...].

Advanced Methods for Natural Products Discovery: Bioactivity Screening, Dereplication, Metabolomics Profiling, Genomic Sequencing, Databases and Informatic Tools, and Structure Elucidation.

Natural Products (NP) are essential for the discovery of novel drugs and products for numerous biotechnological applications. The NP discovery process is expensive and time-consuming, having as major hurdles dereplication (early identification of known compounds) and structure elucidation, particularly the determination of the absolute configuration of metabolites with stereogenic centers. This review comprehensively focuses on recent technological and instrumental advances, highlighting the development of methods that alleviate these obstacles, paving the way for accelerating NP discovery towards biotechnological applications. Herein, we emphasize the most innovative high-throughput tools and methods for advancing bioactivity screening, NP chemical analysis, dereplication, metabolite profiling, metabolomics, genome sequencing and/or genomics approaches, databases, bioinformatics, chemoinformatics, and three-dimensional NP structure elucidation.

Computational Approaches Drive Developments in Immune-Oncology Therapies for PD-1/PD-L1 Immune Checkpoint Inhibitors.

Computational approaches in immune-oncology therapies focus on using data-driven methods to identify potential immune targets and develop novel drug candidates. In particular, the search for PD-1/PD-L1 immune checkpoint inhibitors (ICIs) has enlivened the field, leveraging the use of cheminformatics and bioinformatics tools to analyze large datasets of molecules, gene expression and protein-protein interactions. Up to now, there is still an unmet clinical need for improved ICIs and reliable predictive biomarkers. In this review, we highlight the computational methodologies applied to discovering and developing PD-1/PD-L1 ICIs for improved cancer immunotherapies with a greater focus in the last five years. The use of computer-aided drug design structure- and ligand-based virtual screening processes, molecular docking, homology modeling and molecular dynamics simulations methodologies essential for successful drug discovery campaigns focusing on antibodies, peptides or small-molecule ICIs are addressed. A list of recent databases and web tools used in the context of cancer and immunotherapy has been compilated and made available, namely regarding a general scope, cancer and immunology. In summary, computational approaches have become valuable tools for discovering and developing ICIs. Despite significant progress, there is still a need for improved ICIs and biomarkers, and recent databases and web tools have been compiled to aid in this pursuit.

A computer-aided drug design approach to discover tumour suppressor p53 protein activators for colorectal cancer therapy.

Colorectal cancer (CRC) is the third most detected cancer and the second foremost cause of cancer deaths in the world. Intervention targeting p53 provides potential therapeutic strategies, but thus far no p53-based therapy has been successfully translated into clinical cancer treatment. Here we developed a Quantitative Structure-Activity Relationships (QSAR) classification models using empirical molecular descriptors and fingerprints to predict the activity against the p53 protein, using the potency value with the active or inactive label, were developed. These models were built using in total 10,505 molecules that were extracted from the ChEMBL, ZINC and Reaxys® databases, and recent literature. Three machine learning (ML) techniques e.g., Random Forest, Support Vector Machine, Convolutional Neural Network were explored to build models for p53 inhibitor prediction. The performances of the models were successfully evaluated by internal and external validation. Moreover, based on the best in silico p53 model, a virtual screening campaign was carried out using 1443 FDA-approved drugs that were extracted from the ZINC database. A list of virtual screening hits was assented on base of some limits established in this approach, such as: (1) probability of being active against p53; (2) applicability domain; (3) prediction of the affinity between the p53, and ligands, through molecular docking. The most promising according to the limits established above was dihydroergocristine. This compound revealed cytotoxic activity against a p53-expressing CRC cell line with an IC50 of 56.8 µM. This study demonstrated that the computer-aided drug design approach can be used to identify previously unknown molecules for targeting p53 protein with anti-cancer activity and thus pave the way for the study of a therapeutic solution for CRC.

Predicting Antifouling Activity and Acetylcholinesterase Inhibition of Marine-Derived Compounds Using a Computer-Aided Drug Design Approach.

Biofouling is the undesirable growth of micro- and macro-organisms on artificial water-immersed surfaces, which results in high costs for the prevention and maintenance of this process (billion €/year) for aquaculture, shipping and other industries that rely on coastal and off-shore infrastructure. To date, there are still no sustainable, economical and environmentally safe solutions to overcome this challenging phenomenon. A computer-aided drug design (CADD) approach comprising ligand- and structure-based methods was explored for predicting the antifouling activities of marine natural products (MNPs). In the CADD ligand-based method, 141 organic molecules extracted from the ChEMBL database and literature with antifouling screening data were used to build the quantitative structure-activity relationship (QSAR) classification model. An overall predictive accuracy score of up to 71% was achieved with the best QSAR model for external and internal validation using test and training sets. A virtual screening campaign of 14,492 MNPs from Encinar's website and 14 MNPs that are currently in the clinical pipeline was also carried out using the best QSAR model developed. In the CADD structure-based approach, the 125 MNPs that were selected by the QSAR approach were used in molecular docking experiments against the acetylcholinesterase enzyme. Overall, 16 MNPs were proposed as the most promising marine drug-like leads as antifouling agents, e.g., macrocyclic lactam, macrocyclic alkaloids, indole and pyridine derivatives.

Systematic Review: Drug Repositioning for Congenital Disorders of Glycosylation (CDG).

Advances in research have boosted therapy development for congenital disorders of glycosylation (CDG), a group of rare genetic disorders affecting protein and lipid glycosylation and glycosylphosphatidylinositol anchor biosynthesis. The (re)use of known drugs for novel medical purposes, known as drug repositioning, is growing for both common and rare disorders. The latest innovation concerns the rational search for repositioned molecules which also benefits from artificial intelligence (AI). Compared to traditional methods, drug repositioning accelerates the overall drug discovery process while saving costs. This is particularly valuable for rare diseases. AI tools have proven their worth in diagnosis, in disease classification and characterization, and ultimately in therapy discovery in rare diseases. The availability of biomarkers and reliable disease models is critical for research and development of new drugs, especially for rare and heterogeneous diseases such as CDG. This work reviews the literature related to repositioned drugs for CDG, discovered by serendipity or through a systemic approach. Recent advances in biomarkers and disease models are also outlined as well as stakeholders' views on AI for therapy discovery in CDG.

A Computer-Aided Drug Design Approach to Predict Marine Drug-Like Leads for SARS-CoV-2 Main Protease Inhibition.

The investigation of marine natural products (MNPs) as key resources for the discovery of drugs to mitigate the COVID-19 pandemic is a developing field. In this work, computer-aided drug design (CADD) approaches comprising ligand- and structure-based methods were explored for predicting SARS-CoV-2 main protease (Mpro) inhibitors. The CADD ligand-based method used a quantitative structure-activity relationship (QSAR) classification model that was built using 5276 organic molecules extracted from the ChEMBL database with SARS-CoV-2 screening data. The best model achieved an overall predictive accuracy of up to 67% for an external and internal validation using test and training sets. Moreover, based on the best QSAR model, a virtual screening campaign was carried out using 11,162 MNPs retrieved from the Reaxys® database, 7 in-house MNPs obtained from marine-derived actinomycetes by the team, and 14 MNPs that are currently in the clinical pipeline. All the MNPs from the virtual screening libraries that were predicted as belonging to class A were selected for the CADD structure-based method. In the CADD structure-based approach, the 494 MNPs selected by the QSAR approach were screened by molecular docking against Mpro enzyme. A list of virtual screening hits comprising fifteen MNPs was assented by establishing several limits in this CADD approach, and five MNPs were proposed as the most promising marine drug-like leads as SARS-CoV-2 Mpro inhibitors, a benzo[f]pyrano[4,3-b]chromene, notoamide I, emindole SB beta-mannoside, and two bromoindole derivatives.

Antifouling Napyradiomycins from Marine-Derived Actinomycetes Streptomyces aculeolatus.

The undesired attachment of micro and macroorganisms on water-immersed surfaces, known as marine biofouling, results in severe prevention and maintenance costs (billions €/year) for aquaculture, shipping and other industries that rely on coastal and off-shore infrastructures. To date, there are no sustainable, cost-effective and environmentally safe solutions to address this challenging phenomenon. Therefore, we investigated the antifouling activity of napyradiomycin derivatives that were isolated from actinomycetes from ocean sediments collected off the Madeira Archipelago. Our results revealed that napyradiomycins inhibited ≥80% of the marine biofilm-forming bacteria assayed, as well as the settlement of Mytilus galloprovincialis larvae (EC50 < 5 µg/ml and LC50/EC50 >15), without viability impairment. In silico prediction of toxicity end points are of the same order of magnitude of standard approved drugs and biocides. Altogether, napyradiomycins disclosed bioactivity against marine micro and macrofouling organisms, and non-toxic effects towards the studied species, displaying potential to be used in the development of antifouling products.

Intra-clade metabolomic profiling of MAR4 Streptomyces from the Macaronesia Atlantic region reveals a source of anti-biofilm metabolites.

The search for new and effective strategies to reduce bacterial biofilm formation is of utmost importance as bacterial resistance to antibiotics continues to emerge. The use of anti-biofilm agents that can disrupt recalcitrant bacterial communities can be an advantageous alternative to antimicrobials, as their use does not lead to the development of resistance mechanisms. Six MAR4 Streptomyces strains isolated from the Madeira Archipelago, at the unexplored Macaronesia Atlantic ecoregion, were used to study the chemical diversity of produced hybrid isoprenoids. These marine actinomycetes were investigated by analysing their crude extracts using LC-MS/MS and their metabolomic profiles were compared using multivariate statistical analysis (principal component analysis), showing a separation trend closely related to their phylogeny. Molecular networking unveiled the presence of a class of metabolites not previously described from MAR4 strains and new chemical derivatives belonging to the napyradiomycin and marinone classes. Furthermore, these MAR4 strains produce metabolites that inhibit biofilm formation of Staphylococcus aureus and Marinobacter hydrocarbonoclasticus. The anti-biofilm activity of napyradiomycin SF2415B3 (1) against S. aureus was confirmed.

NavMol 3.0: enabling the representation of metabolic reactions by blind users.

Summary: The representation of metabolic reactions strongly relies on visualization, which is a major barrier for blind users. The NavMol software renders the communication and interpretation of molecular structures and reactions accessible by integrating chemoinformatics and assistive technology. NavMol 3.0 provides a molecular editor for metabolic reactions. The user can start with templates of reactions and build from such cores. Atom-to-atom mapping enables changes in the reactants to be reflected in the products (and vice-versa) and the reaction centres to be automatically identified. Blind users can easily interact with the software using the keyboard and text-to-speech technology. Availability and implementation: NavMol 3.0 is free and open source under the GNU general public license (GPLv3), and can be downloaded at http://sourceforge.net/projects/navmol as a JAR file. Contact: joao@airesdesousa.com.

Computational Methodologies in the Exploration of Marine Natural Product Leads.

Computational methodologies are assisting the exploration of marine natural products (MNPs) to make the discovery of new leads more efficient, to repurpose known MNPs, to target new metabolites on the basis of genome analysis, to reveal mechanisms of action, and to optimize leads. In silico efforts in drug discovery of NPs have mainly focused on two tasks: dereplication and prediction of bioactivities. The exploration of new chemical spaces and the application of predicted spectral data must be included in new approaches to select species, extracts, and growth conditions with maximum probabilities of medicinal chemistry novelty. In this review, the most relevant current computational dereplication methodologies are highlighted. Structure-based (SB) and ligand-based (LB) chemoinformatics approaches have become essential tools for the virtual screening of NPs either in small datasets of isolated compounds or in large-scale databases. The most common LB techniques include Quantitative Structure⁻Activity Relationships (QSAR), estimation of drug likeness, prediction of adsorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, similarity searching, and pharmacophore identification. Analogously, molecular dynamics, docking and binding cavity analysis have been used in SB approaches. Their significance and achievements are the main focus of this review.

A Computer-Driven Approach to Discover Natural Product Leads for Methicillin-Resistant Staphylococcus aureus Infection Therapy.

The risk of methicillin-resistant Staphylococcus aureus (MRSA) infection is increasing in both the developed and developing countries. New approaches to overcome this problem are in need. A ligand-based strategy to discover new inhibiting agents against MRSA infection was built through exploration of machine learning techniques. This strategy is based in two quantitative structure⁻activity relationship (QSAR) studies, one using molecular descriptors (approach A) and the other using descriptors (approach B). In the approach A, regression models were developed using a total of 6645 molecules that were extracted from the ChEMBL, PubChem and ZINC databases, and recent literature. The performance of the regression models was successfully evaluated by internal and external validation, the best model achieved R² of 0.68 and RMSE of 0.59 for the test set. In general natural product (NP) drug discovery is a time-consuming process and several strategies for dereplication have been developed to overcome this inherent limitation. In the approach B, we developed a new NP drug discovery methodology that consists in frontloading samples with 1D NMR descriptors to predict compounds with antibacterial activity prior to bioactivity screening for NPs discovery. The NMR QSAR classification models were built using 1D NMR data (¹H and 13C) as descriptors, from crude extracts, fractions and pure compounds obtained from actinobacteria isolated from marine sediments collected off the Madeira Archipelago. The overall predictability accuracies of the best model exceeded 77% for both training and test sets.

Machine Learning Methods to Predict Density Functional Theory B3LYP Energies of HOMO and LUMO Orbitals.

Machine learning algorithms were explored for the fast estimation of HOMO and LUMO orbital energies calculated by DFT B3LYP, on the basis of molecular descriptors exclusively based on connectivity. The whole project involved the retrieval and generation of molecular structures, quantum chemical calculations for a database with >111 000 structures, development of new molecular descriptors, and training/validation of machine learning models. Several machine learning algorithms were screened, and an applicability domain was defined based on Euclidean distances to the training set. Random forest models predicted an external test set of 9989 compounds achieving mean absolute error (MAE) up to 0.15 and 0.16 eV for the HOMO and LUMO orbitals, respectively. The impact of the quantum chemical calculation protocol was assessed with a subset of compounds. Inclusion of the orbital energy calculated by PM7 as an additional descriptor significantly improved the quality of estimations (reducing the MAE in >30%).

Dereplication: racing to speed up the natural products discovery process.

Covering: 1993-2014 (July)To alleviate the dereplication holdup, which is a major bottleneck in natural products discovery, scientists have been conducting their research efforts to add tools to their "bag of tricks" aiming to achieve faster, more accurate and efficient ways to accelerate the pace of the drug discovery process. Consequently dereplication has become a hot topic presenting a huge publication boom since 2012, blending multidisciplinary fields in new ways that provide important conceptual and/or methodological advances, opening up pioneering research prospects in this field.

QSAR-assisted virtual screening of lead-like molecules from marine and microbial natural sources for antitumor and antibiotic drug discovery.

A Quantitative Structure-Activity Relationship (QSAR) approach for classification was used for the prediction of compounds as active/inactive relatively to overall biological activity, antitumor and antibiotic activities using a data set of 1746 compounds from PubChem with empirical CDK descriptors and semi-empirical quantum-chemical descriptors. A data set of 183 active pharmaceutical ingredients was additionally used for the external validation of the best models. The best classification models for antibiotic and antitumor activities were used to screen a data set of marine and microbial natural products from the AntiMarin database-25 and four lead compounds for antibiotic and antitumor drug design were proposed, respectively. The present work enables the presentation of a new set of possible lead like bioactive compounds and corroborates the results of our previous investigations. By other side it is shown the usefulness of quantum-chemical descriptors in the discrimination of biologically active and inactive compounds. None of the compounds suggested by our approach have assigned non-antibiotic and non-antitumor activities in the AntiMarin database and almost all were lately reported as being active in the literature.

A chemoinformatics approach to the discovery of lead-like molecules from marine and microbial sources en route to antitumor and antibiotic drugs.

The comprehensive information of small molecules and their biological activities in the PubChem database allows chemoinformatic researchers to access and make use of large-scale biological activity data to improve the precision of drug profiling. A Quantitative Structure-Activity Relationship approach, for classification, was used for the prediction of active/inactive compounds relatively to overall biological activity, antitumor and antibiotic activities using a data set of 1804 compounds from PubChem. Using the best classification models for antibiotic and antitumor activities a data set of marine and microbial natural products from the AntiMarin database were screened-57 and 16 new lead compounds for antibiotic and antitumor drug design were proposed, respectively. All compounds proposed by our approach are classified as non-antibiotic and non-antitumor compounds in the AntiMarin database. Recently several of the lead-like compounds proposed by us were reported as being active in the literature.

Investigating the antiviral therapeutic potentialities of marine polycyclic lamellarin pyrrole alkaloids as promising inhibitors for SARS-CoV-2 and Zika main proteases (Mpro).

The new coronavirus variant (SARS-CoV-2) and Zika virus are two world-wide health pandemics. Along history, natural products-based drugs have always crucially recognized as a main source of valuable medications. Considering the SARS-CoV-2 and Zika main proteases (Mpro) as the re-production key element of the viral cycle and its main target, herein we report an intensive computer-aided virtual screening for a focused list of 39 marine lamellarins pyrrole alkaloids, against SARS-CoV-2 and Zika main proteases (Mpro) using a set of combined modern computational methodologies including molecular docking (MDock), molecule dynamic simulations (MDS) and structure-activity relationships (SARs) as well. Indeed, the molecular docking studies had revealed four promising marine alkaloids including [lamellarin H (14)/K (17)] and [lamellarin S (26)/Z (39)], according to their notable ligand-protein energy scores and relevant binding affinities with the SARS-CoV-2 and Zika (Mpro) pocket residues, respectively. Consequentially, these four chemical hits were further examined thermodynamically though investigating their MD simulations at 100 ns, where they showed prominent stability within the accommodated (Mpro) pockets. Moreover, in-deep SARs studies suggested the crucial roles of the rigid fused polycyclic ring system, particularly aromatic A- and F- rings, position of the phenolic -OH and δ-lactone functionalities as essential structural and pharmacophoric features. Finally, these four promising lamellarins alkaloids were investigated for their in-silico ADME using the SWISS ADME platform, where they displayed appropriated drug-likeness properties. Such motivating outcomes are greatly recommending further in vitro/vivo examinations regarding those lamellarins pyrrole alkaloids (LPAs).Communicated by Ramaswamy H. Sarma.

Investigating the hepatoprotective potentiality of marine-derived steroids as promising inhibitors of liver fibrosis.

It has been reported that organic extracts derived from soft corals belonging to the genus Sarcophyton have exhibited a wide range of therapeutic characteristics. Based on biochemical and histological techniques, we aimed to assess the hepatoprotective role of the organic extract and its principal steroidal contents derived from the Red Sea soft coral Sarcophyton glaucum on acetaminophen-induced liver fibrosis in rats. Serum liver function parameters (ALT, AST, ALP and total bilirubin) were quantified using a spectrophotometer, and both alpha-fetoprotein (AFP) and carcinoembryonic antigen (CEA) levels were determined by using enzyme-linked immunosorbent assay (ELISA) kits while transformed growth factor beta (TGF-ß) and tumor necrosis factor α (TNF-α) in liver tissue homogenate were determined using ELISA, and TGF-ß and TNF-α gene expression in liver tissue was determined using real-time PCR following extraction and purification. Histopathological alterations in hepatic tissue were also examined under a microscope. In order to prioritize the isolation and characterization of the most promising marine steroids from the organic extract of the Red Sea soft coral Sarcophyton glaucum as hepatoprotective agents, a computational approach was employed. This approach involved molecular docking (MDock) and analysis of the structure-activity relationship (SAR) against glutathione-S-transferase (GST) and Cu-Zn human superoxide dismutase (Cu-ZnSOD) enzymes. Although the major role in the detoxification of foreign chemicals and toxic metabolites of GST and SOD enzymes is known, there is a lack of knowledge about the mode of action of the hepatoprotective process and those of the targets involved. The present study investigated the multiple interactions of a series of marine steroids with the GST and SOD enzymes, in order to reveal insights into the process of hepatoprotection.

1D 13C-NMR data as molecular descriptors in spectra--structure relationship analysis of oligosaccharides.

Spectra-structure relationships were investigated for estimating the anomeric configuration, residues and type of linkages of linear and branched trisaccharides using 13C-NMR chemical shifts. For this study, 119 pyranosyl trisaccharides were used that are trimers of the α or ß anomers of D-glucose, D-galactose, D-mannose, L-fucose or L-rhamnose residues bonded through a or b glycosidic linkages of types 1→2, 1→3, 1→4, or 1→6, as well as methoxylated and/or N-acetylated amino trisaccharides. Machine learning experiments were performed for: (1) classification of the anomeric configuration of the first unit, second unit and reducing end; (2) classification of the type of first and second linkages; (3) classification of the three residues: reducing end, middle and first residue; and (4) classification of the chain type. Our previously model for predicting the structure of disaccharides was incorporated in this new model with an improvement of the predictive power. The best results were achieved using Random Forests with 204 di- and trisaccharides for the training set-it could correctly classify 83%, 90%, 88%, 85%, 85%, 75%, 79%, 68% and 94% of the test set (69 compounds) for the nine tasks, respectively, on the basis of unassigned chemical shifts.

Investigating the structure-activity relationship of marine polycyclic batzelladine alkaloids as promising inhibitors for SARS-CoV-2 main protease (Mpro).

Over a span of two years ago, since the emergence of the first case of the novel coronavirus (SARS-CoV-2) in China, the pandemic has crossed borders causing serious health emergencies, immense economic crisis and impacting the daily life worldwide. Despite the discovery of numerous forms of precautionary vaccines along with other recently approved orally available drugs, yet effective antiviral therapeutics are necessarily needed to hunt this virus and its variants. Historically, naturally occurring chemicals have always been considered the primary source of beneficial medications. Considering the SARS-CoV-2 main protease (Mpro) as the duplicate key element of the viral cycle and its main target, in this paper, an extensive virtual screening for a focused chemical library of 15 batzelladine marine alkaloids, was virtually examined against SARS-CoV-2 main protease (Mpro) using an integrated set of modern computational tools including molecular docking (MDock), molecule dynamic (MD) simulations and structure-activity relationships (SARs) as well. The molecular docking predictions had disclosed four promising compounds including batzelladines H-I (8-9) and batzelladines F-G (6-7), respectively according to their prominent ligand-protein energy scores and relevant binding affinities with the (Mpro) pocket residues. The best two chemical hits, batzelladines H-I (8-9) were further investigated thermodynamically though studying their MD simulations at 100 ns, where they showed excellent stability within the accommodated (Mpro) pocket. Moreover, SARs studies imply the crucial roles of the fused tricyclic guanidinic moieties, its degree of unsaturation, position of the N-OH functionality and the length of the side chain as a spacer linking between two active sites, which disclosed fundamental structural and pharmacophoric features for efficient protein-ligand interaction. Such interesting findings are greatly highlighting further in vitro/vivo examinations regarding those marine natural products (MNPs) and their synthetic equivalents as promising antivirals.