Molecular docking is the most frequently used computational method for studying the interactions between organic molecules and biological macromolecules. In this context, docking allows predicting the preferred pose of a ligand inside a receptor binding site. However, the selection of the “best” solution is not a trivial task, despite the widely accepted selection criterion that the best pose corresponds to the best energy score. Here, several rigid-target docking methods were evaluated on the same dataset with respect to their ability to reproduce crystallographic binding orientations, to test if the best energy score is a reliable criterion for selecting the best solution. For this, two experiments were performed: (A) to reconstruct the ligand-receptor complex by performing docking of the ligand in its own crystal structure receptor (defined as self-docking), and (B) to reconstruct the ligand-receptor complex by performing docking of the ligand in a crystal structure receptor that contains other ligand (defined as cross-docking). Root-mean square deviation (RMSD) was used to evaluate how different the obtained docking orientation is from the corresponding co-crystallized pose of the same ligand molecule. We found that docking score function is capable of predicting crystallographic binding orientations, but the best ranked solution according to the docking energy is not always the pose that reproduces the experimental binding orientation. This happened when self-docking was achieved, but it was critical in cross-docking. Taking into account that docking is typically used with predictive purposes, during cross-docking experiments, our results indicate that the best energy score is not a reliable criterion to select the best solution in common docking applications. It is strongly recommended to choose the best docking solution according to the scoring function along with additional structural criteria described for analogue ligands to assure the selection of a correct docking solution.

Although aromaticity is a concept in chemistry, in the last years, special efforts have been carried out in order to propose theoretical strategies to quantify it as a property of molecular rings. Among them, perhaps the computation of nucleus independent chemical shifts (NICSs) is the most commonly used, since it is possible to calculate it in an easy and fast way with most used quantum chemistry software. However, contradicting assignments of aromaticity by NICS and other methods have been reported in the literature, especially in studies concerning inorganic chemistry. In this Article is proposed a new and simple strategy to use the NICS information to assess aromaticity, identifying the point along the axis perpendicular to the molecular plane where the in-plane component of NICS becomes zero; it is called free of in-plane component NICS (FiPC-NICS). This spatial point is proposed as secure to consider NICS as an aromaticity descriptor; this simple proposal is evaluated in borazine and cyclotriphosphazenes. The results are compared with carefully examined aromatic stabilization energies and magnetically induced current-density analysis.

Molecular docking is a computational chemistry method which has become essential for the rational drug design process. In this context, it has had great impact as a successful tool for the study of ligand–receptor interaction modes, and for the exploration of large chemical datasets through virtual screening experiments. Despite their unquestionable merits, docking methods are not reliable for predicting binding energies due to the simple scoring functions they use. However, comparisons between two or three complexes using the predicted binding energies as a criterion are commonly found in the literature. In the present work we tested how wise is it to trust the docking energies when two complexes between a target protein and enantiomer pairs are compared. For this purpose, a ligand library composed by 141 enantiomeric pairs was used, including compounds with biological activities reported against seven protein targets. Docking results using the software Glide (considering extra precision (XP), standard precision (SP), and high-throughput virtual screening (HTVS) modes) and AutoDock Vina were compared with the reported biological activities using a classification scheme. Our test failed for all modes and targets, demonstrating that an accurate prediction when binding energies of enantiomers are compared using docking may be due to chance. We also compared pairs of compounds with different molecular weights and found the same results.

Motivation Root mean square deviation (RMSD) is one of the most useful and straightforward features for structural comparison between different conformations of the same molecule. Commonly, protein-ligand docking programs have included some utilities that allow the calculation of this value; however, they only work efficiently when exists a complete atom label equivalence between the evaluated conformations. Results We present LigRMSD, a free web-server for the automatic matching and RMSD calculations among identical or similar chemical compounds. This server allows the user to submit only a pair of identical or similar molecules or dataset of similar compounds to compare their three-dimensional conformations. Availability and implementation LigRMSD can be freely accessed at https://ligrmsd.appsbio.utalca.cl. Supplementary information Supplementary data are available at Bioinformatics online.

Liquid chromatography coupled with mass spectrometry is an outstanding methodology for fast analysis of phenolic compounds in biological samples. Twenty two compounds were quickly and accurately identified in the methanolic extract of the Antarctic lichen Ramalina terebrata for the first time using ultra high pressure liquid chromatography coupled with photodiode array detector and high resolution mass spectrometry (UHPLC-PDA-Q/Orbitrap/MS/MS). In addition, the extract and the four compounds isolated from this species were tested for the inhibitory activity of tau protein aggregation, which is a protein involved in Alzheimer’s disease (AD). All compounds showed null activity with the exception of parietin, which it was able to inhibit aggregation process of tau in a concentration range between 3 µg/mL (10 µM) to 28 µg/mL (100 µM). In addition, we show how parietin interact with tau 306VQIVYK311 hexapeptide inside of the microtubule binding domain (4R) with the help of molecular docking experiments. Finally, the constituents present in the methanolic extract could possibly contribute to the established anti-aggregation activity for this extract and this in-depth analysis of the chemical composition of R. terebrata could guide further research into its medicinal properties and potential uses.