Building nonparametric nn-body force fields using Gaussian process regression

Constructing a classical potential suited to simulate a given atomic system is a remarkably difficult task. This chapter presents a framework under which this problem can be tackled, based on the Bayesian construction of nonparametric force fields of a given order using Gaussian process (GP) priors. The formalism of GP regression is first reviewed, particularly in relation to its application in learning local atomic energies and forces. For accurate regression it is fundamental to incorporate prior knowledge into the GP kernel function. To this end, this chapter details how properties of smoothness, invariance and interaction order of a force field can be encoded into corresponding kernel properties. A range of kernels is then proposed, possessing all the required properties and an adjustable parameter $n$ governing the interaction order modelled. The order $n$ best suited to describe a given system can be found automatically within the Bayesian framework by maximisation of the marginal likelihood. The procedure is first tested on a toy model of known interaction and later applied to two real materials described at the DFT level of accuracy. The models automatically selected for the two materials were found to be in agreement with physical intuition. More in general, it was found that lower order (simpler) models should be chosen when the data are not sufficient to resolve more complex interactions. Low $n$ GPs can be further sped up by orders of magnitude by constructing the corresponding tabulated force field, here named "MFF".Comment: 31 pages, 11 figures, book chapte

Identification of New Inhibitors with Potential Antitumor Activity from Polypeptide Structures via Hierarchical Virtual Screening.

Leukemias are neoplasms that affect hematopoietic cells, which are developed by genetic alterations (mutations) that lead to the loss of proliferation control mechanisms (maturation and/or cell death). The α4β1 integrin receptor is a therapeutic target for inflammation, autoimmune diseases and lymphoid tumors. This study was carried out to search through the antagonists-based virtual screening for α4β1 receptor. Initially, seventeen (17) structures were selected (based on the inhibitory activity values, IC50) and the structure with the best value was chosen as the pivot. The pharmacophoric pattern was determined from the online PharmaGist server and resulted in a model of score value equal to 97.940 with 15 pharmacophoric characteristics that were statistically evaluated via Pearson correlations, principal component analysis (PCA) and hierarchical clustering analysis (HCA). A refined model generated four pharmacophoric hypotheses totaling 1.478 structures set of Zinc_database. After, the pharmacokinetic, toxicological and biological activity predictions were realized comparing with pivot structure that resulted in five (ZINC72088291, ZINC68842860, ZINC14365931, ZINC09588345 and ZINC91247798) structures with optimal in silico predictions. Therefore, future studies are needed to confirm antitumor potential activity of molecules selected this work with in vitro and in vivo assays

Identification of Potential New Aedes aegypti Juvenile Hormone Inhibitors from N-Acyl Piperidine Derivatives: A Bioinformatics Approach

Aedes aegypti mosquitoes transmit several human pathogens that cause millions of deaths worldwide, mainly in Latin America. The indiscriminate use of insecticides has resulted in the development of species resistance to some such compounds. Piperidine, a natural alkaloid isolated from Piper nigrum, has been used as a hit compound due to its larvicidal activity against Aedes aegypti. In the present study, piperidine derivatives were studied through in silico methods: pharmacophoric evaluation (PharmaGist), pharmacophoric virtual screening (Pharmit), ADME/Tox prediction (Preadmet/Derek 10.0®), docking calculations (AutoDock 4.2) and molecular dynamics (MD) simulation on GROMACS-5.1.4. MP-416 and MP-073 molecules exhibiting ΔG binding (MMPBSA −265.95 ± 1.32 kJ/mol and −124.412 ± 1.08 kJ/mol, respectively) and comparable to holo (ΔG binding = −216.21 ± 0.97) and pyriproxyfen (a well-known larvicidal, ΔG binding= −435.95 ± 2.06 kJ/mol). Considering future in vivo assays, we elaborated the theoretical synthetic route and made predictions of the synthetic accessibility (SA) (SwissADME), lipophilicity and water solubility (SwissADME) of the promising compounds identified in the present study. Our in silico results show that MP-416 and MP-073 molecules could be potent insecticides against the Aedes aegypti mosquitoes