Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Spyrmsd: symmetry-corrected RMSD calculations in Python.

ABSTRACT: Root mean square displacement (RMSD) calculations play a fundamental role in the comparison of different conformers of the same ligand. This is particularly important in the evaluation of protein-ligand docking, where different ligand poses are generated by docking software and their quality is usually assessed by RMSD calculations. Unfortunately, many RMSD calculation tools do not take into account the symmetry of the molecule, remain difficult to integrate flawlessly in cheminformatics and machine learning pipelines-which are often written in Python-or are shipped within large code bases. Here we present a new open-source RMSD calculation tool written in Python, designed to be extremely lightweight and easy to integrate into existing software.

SUBMITTER: Meli R

PROVIDER: S-EPMC7457508 | biostudies-literature | 2020 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

spyrmsd: symmetry-corrected RMSD calculations in Python.

Meli Rocco R Biggin Philip C PC

Journal of cheminformatics 20200831 1

Root mean square displacement (RMSD) calculations play a fundamental role in the comparison of different conformers of the same ligand. This is particularly important in the evaluation of protein-ligand docking, where different ligand poses are generated by docking software and their quality is usually assessed by RMSD calculations. Unfortunately, many RMSD calculation tools do not take into account the symmetry of the molecule, remain difficult to integrate flawlessly in cheminformatics and mac ...[more]

PMID: 33431033

Similar Datasets

A Python tool to set up relative free energy calculations in GROMACS.

Project description:Free energy calculations based on molecular dynamics (MD) simulations have seen a tremendous growth in the last decade. However, it is still difficult and tedious to set them up in an automated manner, as the majority of the present-day MD simulation packages lack that functionality. Relative free energy calculations are a particular challenge for several reasons, including the problem of finding a common substructure and mapping the transformation to be applied. Here we present a tool, alchemical-setup.py, that automatically generates all the input files needed to perform relative solvation and binding free energy calculations with the MD package GROMACS. When combined with Lead Optimization Mapper (LOMAP; Liu et al. in J Comput Aided Mol Des 27(9):755-770, 2013), recently developed in our group, alchemical-setup.py allows fully automated setup of relative free energy calculations in GROMACS. Taking a graph of the planned calculations and a mapping, both computed by LOMAP, our tool generates the topology and coordinate files needed to perform relative free energy calculations for a given set of molecules, and provides a set of simulation input parameters. The tool was validated by performing relative hydration free energy calculations for a handful of molecules from the SAMPL4 challenge (Mobley et al. in J Comput Aided Mol Des 28(4):135-150, 2014). Good agreement with previously published results and the straightforward way in which free energy calculations can be conducted make alchemical-setup.py a promising tool for automated setup of relative solvation and binding free energy calculations.

| S-EPMC4749424 | biostudies-literature

FragBuilder: an efficient Python library to setup quantum chemistry calculations on peptides models.

Project description:We present a powerful Python library to quickly and efficiently generate realistic peptide model structures. The library makes it possible to quickly set up quantum mechanical calculations on model peptide structures. It is possible to manually specify a specific conformation of the peptide. Additionally the library also offers sampling of backbone conformations and side chain rotamer conformations from continuous distributions. The generated peptides can then be geometry optimized by the MMFF94 molecular mechanics force field via convenient functions inside the library. Finally, it is possible to output the resulting structures directly to files in a variety of useful formats, such as XYZ or PDB formats, or directly as input files for a quantum chemistry program. FragBuilder is freely available at https://github.com/jensengroup/fragbuilder/ under the terms of the BSD open source license.

| S-EPMC3961104 | biostudies-literature

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers.

Project description:Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

| S-EPMC4612423 | biostudies-literature

Validation of experimental molecular crystal structures with dispersion-corrected density functional theory calculations.

Project description:This paper describes the validation of a dispersion-corrected density functional theory (d-DFT) method for the purpose of assessing the correctness of experimental organic crystal structures and enhancing the information content of purely experimental data. 241 experimental organic crystal structures from the August 2008 issue of Acta Cryst. Section E were energy-minimized in full, including unit-cell parameters. The differences between the experimental and the minimized crystal structures were subjected to statistical analysis. The r.m.s. Cartesian displacement excluding H atoms upon energy minimization with flexible unit-cell parameters is selected as a pertinent indicator of the correctness of a crystal structure. All 241 experimental crystal structures are reproduced very well: the average r.m.s. Cartesian displacement for the 241 crystal structures, including 16 disordered structures, is only 0.095 Å (0.084 Å for the 225 ordered structures). R.m.s. Cartesian displacements above 0.25?A either indicate incorrect experimental crystal structures or reveal interesting structural features such as exceptionally large temperature effects, incorrectly modelled disorder or symmetry breaking H atoms. After validation, the method is applied to nine examples that are known to be ambiguous or subtly incorrect.

| S-EPMC2940256 | biostudies-literature

On the use of orientational restraints and symmetry corrections in alchemical free energy calculations.

Project description:Alchemical free energy calculations are becoming a useful tool for calculating absolute binding free energies of small molecule ligands to proteins. Here, we find that the presence of multiple metastable ligand orientations can cause convergence problems when distance restraints alone are used. We demonstrate that the use of orientational restraints can greatly accelerate the convergence of these calculations. However, even with this acceleration, we find that sufficient sampling requires substantially longer simulations than are used in many published protocols. To further accelerate convergence, we introduce a new method of configuration space decomposition by orientation which reduces required simulation lengths by at least a factor of 5 in the cases examined. Our method is easily parallelizable, well suited for cases where a ligand cocrystal structure is not available, and can utilize initial orientations generated by docking packages.

| S-EPMC3583553 | biostudies-literature

CARON--average RMSD of NMR structure ensembles.

Project description:The NMR protein structures are often deposited in the Protein Data Bank as ensembles of models that agree with the experimental restraints. Information about stereochemical variability and the molecular flexibility can be obtained by systematic comparison of all models. Here we describe CARON, a software that allows the computation of the root-mean-square-distances between equivalent atoms and residues in all models and introduces these values into the occupancy and the B-factor fields of PDB-formatted files. This tool allows the user to both get a quantitative estimation of the conformational homogeneity of the models and to exploit this information in common computer graphics programs.

| S-EPMC2828896 | biostudies-literature

Relaxation estimation of RMSD in molecular dynamics immunosimulations.

Project description:Molecular dynamics simulations have to be sufficiently long to draw reliable conclusions. However, no method exists to prove that a simulation has converged. We suggest the method of "lagged RMSD-analysis" as a tool to judge if an MD simulation has not yet run long enough. The analysis is based on RMSD values between pairs of configurations separated by variable time intervals Δt. Unless RMSD(Δt) has reached a stationary shape, the simulation has not yet converged.

| S-EPMC3457668 | biostudies-literature