Search

de novo Drug Design (dnDD) aims to create new molecules that satisfy multiple conflicting objectives. Since several desired properties can be considered in the optimization process, dnDD is naturally categorized as a many-objective optimization problem (ManyOOP), where more than three objectives must be simultaneously optimized. However, a large number of objectives typically pose several challenges that affect the choice and the design of optimization methodologies. Herein, we cover the application of multi- and many-objective optimization methods, particularly those based on Evolutionary Computation and Machine Learning techniques, to enlighten their potential application in dnDD. Additionally, we comprehensively analyze how molecular properties used in the optimization process are applied as either objectives or constraints to the problem. Finally, we discuss future research in many-objective optimization for dnDD, highlighting two important possible impacts: i) its integration with the development of multi-target approaches to accelerate the discovery of innovative and more efficacious drug therapies and ii) its role as a catalyst for new developments in more fundamental and general methodological frameworks in the field.

Abstract The COVID-19 caused by the SARS-CoV-2 virus was declared a pandemic disease in March 2020 by the World Health Organization (WHO). Structure-Based Drug Design strategies based on docking methodologies have been widely used for both new drug development and drug repurposing to find effective treatments against this disease. In this work, we present the developments implemented in the DockThor-VS web server to provide a virtual screening (VS) platform with curated structures of potential therapeutic targets from SARS-CoV-2 incorporating genetic information regarding relevant non-synonymous variations. The web server facilitates repurposing VS experiments providing curated libraries of currently available drugs on the market. At present, DockThor-VS provides ready-for-docking 3D structures for wild type and selected mutations for Nsp3 (papain-like, PLpro domain), Nsp5 (Mpro, 3CLpro), Nsp12 (RdRp), Nsp15 (NendoU), N protein, and Spike. We performed VS experiments of FDA-approved drugs considering the therapeutic targets available at the web server to assess the impact of considering different structures and mutations to identify possible new treatments of SARS-CoV-2 infections. The DockThor-VS is freely available at www.dockthor.lncc.br .

Abstract Scoring functions are essential for modern in silico drug discovery. However, the accurate prediction of binding affinity by scoring functions remains a challenging task. The performance of scoring functions is very heterogeneous across different target classes. Scoring functions based on precise physics-based descriptors better representing protein–ligand recognition process are strongly needed. We developed a set of new empirical scoring functions, named DockTScore, by explicitly accounting for physics-based terms combined with machine learning. Target-specific scoring functions were developed for two important drug targets, proteases and protein–protein interactions, representing an original class of molecules for drug discovery. Multiple linear regression (MLR), support vector machine and random forest algorithms were employed to derive general and target-specific scoring functions involving optimized MMFF94S force-field terms, solvation and lipophilic interactions terms, and an improved term accounting for ligand torsional entropy contribution to ligand binding. DockTScore scoring functions demonstrated to be competitive with the current best-evaluated scoring functions in terms of binding energy prediction and ranking on four DUD-E datasets and will be useful for in silico drug design for diverse proteins as well as for specific targets such as proteases and protein–protein interactions. Currently, the MLR DockTScore is available at www.dockthor.lncc.br .

Abstract Background There are limited long-term outcome data in eculizumab-treated patients with atypical hemolytic uremic syndrome (aHUS). We report final results from the largest prospective, observational, multicenter study of patients with aHUS treated with eculizumab. Methods Patients with aHUS who participated in any of five parent eculizumab trials and received at least one eculizumab infusion were eligible for enrollment in a long-term follow-up study. Rates of thrombotic microangiopathy (TMA) manifestations off versus on eculizumab were evaluated. Additional endpoints included change from baseline estimated glomerular filtration rate (eGFR), long-term renal outcomes, and serious targeted treatment-emergent adverse events. Results Among 93 patients (0–80 years of age), 51 (55%) remained on eculizumab and 42 (45%) discontinued; for those who discontinued, 21 (50%) reinitiated therapy. Patients who reinitiated eculizumab had similar baseline clinical characteristics to patients who remained on eculizumab, with higher likelihood of genetic/autoimmune complement abnormalities, more prior TMAs, and longer disease course versus those who did not reinitiate. Mean eGFR improved rapidly and remained stable for up to 6 years on eculizumab. In patients who discontinued, there was a trend toward decreasing renal function over time from discontinuation. Additionally, off-treatment TMA manifestation rates were higher in those aged

Structure-based virtual screening (VS) is a widely used approach that employs the knowledge of the three-dimensional structure of the target of interest in the design of new lead compounds from large-scale molecular docking experiments. Through the prediction of the binding mode and affinity of a small molecule within the binding site of the target of interest, it is possible to understand important properties related to the binding process. Empirical scoring functions are widely used for pose and affinity prediction. Although pose prediction is performed with satisfactory accuracy, the correct prediction of binding affinity is still a challenging task and crucial for the success of structure-based VS experiments. There are several efforts in distinct fronts to develop even more sophisticated and accurate models for filtering and ranking large libraries of compounds. This paper will cover some recent successful applications and methodological advances, including strategies to explore the ligand entropy and solvent effects, training with sophisticated machine-learning techniques, and the use of quantum mechanics. Particular emphasis will be given to the discussion of critical aspects and further directions for the development of more accurate empirical scoring functions.

The current drug options for the treatment of chronic Chagas disease have not been sufficient and high hopes have been placed on the use of genomic data from the human parasite Trypanosoma cruzi to identify new drug targets and develop appropriate treatments for both acute and chronic Chagas disease. However, the lack of a complete assembly of the genomic sequence and the presence of many predicted proteins with unknown or unsure functions has hampered our complete view of the parasite's metabolic pathways. Moreover, pinpointing new drug targets has proven to be more complex than anticipated and has revealed large holes in our understanding of metabolic pathways and their integrated regulation, not only for this parasite, but for many other similar pathogens. Using an in silicocomparative study on pathway annotation and searching for analogous and specific enzymes, we have been able to predict a considerable number of additional enzymatic functions in T. cruzi. Here we focus on the energetic pathways, such as glycolysis, the pentose phosphate shunt, the Krebs cycle and lipid metabolism. We point out many enzymes that are analogous to those of the human host, which could be potential new therapeutic targets.

We analyzed the performance of a real coded "steady-state" genetic algorithm (SSGA) using a grid-based methodology in docking five HIV-1 protease-ligand complexes having known three-dimensional structures. All ligands tested are highly flexible, having more than 10 conformational degrees of freedom. The SSGA was tested for the rigid and flexible ligand docking cases. The implemented genetic algorithm was able to dock successfully rigid and flexible ligand molecules, but with a decreasing performance when the number of ligand conformational degrees of freedom increased. The docked lowest-energy structures have root mean square deviation (RMSD) with respect to the corresponding experimental crystallographic structure ranging from 0.037 Å to 0.090 Å in the rigid docking, and 0.420 Å to 1.943 Å in the flexible docking. We found that not only the number of ligand conformational degrees of freedom is an important aspect to the algorithm performance, but also that the more internal dihedral angles are critical. Furthermore, our results showed that the initial population distribution can be relevant for the algorithm performance.

An approach to the hydrophobic-polar (HP) protein folding model was developed using a genetic algorithm (GA) to find the optimal structures on a 3D cubic lattice. A modification was introduced to the scoring system of the original model to improve the model's capacity to generate more natural-like structures. The modification was based on the assumption that it may be preferable for a hydrophobic monomer to have a polar neighbor than to be in direct contact with the polar solvent. The compactness and the segregation criteria were used to compare structures created by the original HP model and by the modified one. An islands' algorithm, a new selection scheme and multiple-points crossover were used to improve the performance of the algorithm. Ten sequences, seven with length 27 and three with length 64 were analyzed. Our results suggest that the modified model has a greater tendency to form globular structures. This might be preferable, since the original HP model does not take into account the positioning of long polar segments. The algorithm was implemented in the form of a program with a graphical user interface that might have a didactical potential in the study of GA and on the understanding of hydrophobic core formation.

The main goal of this study is to find the most effective set of parameters for the Simplified Generalized Simulated Annealing algorithm, SGSA, when applied to distinct cost function as well as to find a possible correlation between the values of these parameters sets and some topological characteristics of the hypersurface of the respective cost function. The SGSA algorithm is an extended and simplified derivative of the GSA algorithm, a Markovian stochastic process based on Tsallis statistics that has been used in many classes of problems, in particular, in biological molecular systems optimization. In all but one of the studied cost functions, the global minimum was found in 100% of the 50 runs. For these functions the best visiting parameter, qV, belongs to the interval [1.2, 1.7]. Also, the temperature decaying parameter, qT, should be increased when better precision is required. Moreover, the similarity in the locus of optimal parameter sets observed in some functions indicates that possibly one could extract topological information about the cost functions from these sets.

We report the first successful extraction of oyster DNA from a pearl and use it to identify the source oyster species for the three major pearl-producing oyster species Pinctada margaritifera, P. maxima and P. radiata. Both mitochondrial and nuclear gene fragments could be PCR-amplified and sequenced. A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay in the internal transcribed spacer (ITS) region was developed and used to identify 18 pearls of unknown origin. A micro-drilling technique was developed to obtain small amounts of DNA while maintaining the commercial value of the pearls. This DNA fingerprinting method could be used to document the source of historic pearls and will provide more transparency for traders and consumers within the pearl industry.

Emeralds from the Musakashi area of Zambia entered the gem trade in early 2002 in rather small quantities. Faceted stones are mostly available in smaller sizes (1-5 ct) but occasionally reach up to 20 ct or more. They can be of very fine gem quality and exhibit similarities to Colombian emeralds, not only in appearance, but also in their gemmological properties, inclusion features and chemical properties. Seventy-eight Musakashi emeralds were characterised gemmologically and chemically for this report. They can be distinguished from Colombian emeralds (as well as other similar emeralds from Afghanistan's Panjshir Valley) based on certain diagnostic internal features (e.g. 'sawtooth'-outlined multiphase fluid inclusions) and their chemical composition (e.g. plotting Fe vs V/Cr, and cluster visualisation using a t-SNE machine-learning algorithm). Separating them from the classic schist-hosted emeralds of the Kafubu area in Zambia is straightforward based on differences in their gemmological properties, inclusions, UV-Vis-NIR and Raman spectra, and chemical composition. [ABSTRACT FROM AUTHOR]