ConfID: an analytical method for conformational characterization of small molecules using molecular dynamics trajectories.

MOTIVATION: The conformational space of small molecules can be vast and difficult to assess. Molecular dynamics (MD) simulations of free ligands in solution have been applied to predict conformational populations, but their characterization is often based on clustering algorithms or manual efforts. RESULTS: Here, we introduce ConfID, an analytical tool for conformational characterization of small molecules using MD trajectories. The evolution of conformational sampling and population frequencies throughout trajectories is calculated to check for sampling convergence while allowing to map relevant conformational transitions. The tool is designed to track conformational transition events and calculate time-dependent properties for each conformational population detected. AVAILABILITY AND IMPLEMENTATION: Toolkit and documentation are freely available at http://sbcb.inf.ufrgs.br/confid. CONTACT: marcelo.poleto@ufv.br or bigrisci@inf.ufrgs.br. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Development of GROMOS-Compatible Parameter Set for Simulations of Chalcones and Flavonoids.

Chalcones and flavonoids constitute a large family of plant secondary metabolites that have been explored as a potential source of novel pharmaceutical products. While the simulation of these compounds by molecular dynamics (MD) can be a valuable strategy to assess their conformational properties and so further develop their role in drug discovery, there are no set of force field parameters specifically designed and experimentally validated for their conformational description in condensed phase. So the current work developed a new parameter set for MD simulations of these compounds' main scaffolds under GROMOS force field. We employed a protocol adjusting the atomic charges and torsional parameters to the respective quantum mechanical derived dipole moments and dihedrals rotational profiles, respectively. Experimental properties of organic liquids were used as references to the calculated values to validate the parameters. Additionally, metadynamics simulations were performed to evaluate the conformational space of complex chalcones and flavonoids, while NOE contacts during simulations were measured and compared to experimental data. Accordingly, the employed protocol allowed us to obtain force field parameters that reproduce well the target data and may be expected to contribute in more accurate computational studies on the biological/therapeutical role of such molecules.

Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities.

The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.