Clustering molecular dynamics trajectories for optimizing docking experiments.
Ontology highlight
ABSTRACT:
SUBMITTER: De Paris R
PROVIDER: S-EPMC4385651 | biostudies-literature | 2015
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Determining the optimal number and identity of structural clusters from an ensemble of molecular configurations continues to be a challenge. Recent structural clustering methods have focused on the use of internal coordinates due to the innate rotational and translational invariance of these features. The vast number of possible internal coordinates necessitates a feature space supervision step to make clustering tractable but yields a protocol that can be system type-specific. Particle positions offer an appealing alternative to internal coordinates but suffer from a lack of rotational and translational invariance, as well as a perceived insensitivity to regions of structural dissimilarity. Here, we present a method, denoted shape-GMM, that overcomes the shortcomings of particle positions using a weighted maximum likelihood alignment procedure. This alignment strategy is then built into an expectation maximization Gaussian mixture model (GMM) procedure to capture metastable states in the free-energy landscape. The resulting algorithm distinguishes between a variety of different structures, including those indistinguishable by root-mean-square displacement and pairwise distances, as demonstrated on several model systems. Shape-GMM results on an extensive simulation of the fast-folding HP35 Nle/Nle mutant protein support a four-state folding/unfolding mechanism, which is consistent with previous experimental results and provides kinetic details comparable to previous state-of-the art clustering approaches, as measured by the VAMP-2 score. Currently, training of shape-GMMs is recommended for systems (or subsystems) that can be represented by ≲200 particles and ≲100k configurations to estimate high-dimensional covariance matrices and balance computational expense. Once a shape-GMM is trained, it can be used to predict the cluster identities of millions of configurations.
| S-EPMC9228201 | biostudies-literature
Project description:We introduce automatic clustering as a computationally efficient tool for classifying and interpreting trajectories from simulations of photo-excited dynamics. Trajectories are treated as time-series data, with the features for clustering selected by variance mapping of normalized data. The L2-norm and dynamic time warping are proposed as suitable similarity measures for calculating the distance matrices, and these are clustered using the unsupervised density-based DBSCAN algorithm. The silhouette coefficient and the number of trajectories classified as noise are used as quality measures for the clustering. The ability of clustering to provide rapid overview of large and complex trajectory data sets, and its utility for extracting chemical and physical insight, is demonstrated on trajectories corresponding to the photochemical ring-opening reaction of 1,3-cyclohexadiene, noting that the clustering can be used to generate reduced dimensionality representations in an unbiased manner.
| S-EPMC10536988 | biostudies-literature
Project description:BACKGROUND:In the rational drug design process, an ensemble of conformations obtained from a molecular dynamics simulation plays a crucial role in docking experiments. Some studies have found that Fully-Flexible Receptor (FFR) models predict realistic binding energy accurately and improve scoring to enhance selectiveness. At the same time, methods have been proposed to reduce the high computational costs involved in considering the explicit flexibility of proteins in receptor-ligand docking. This study introduces a novel method to optimize ensemble docking-based experiments by reducing the size of an InhA FFR model at docking runtime and scaling docking workflow invocations on cloud virtual machines. RESULTS:First, in order to find the most affordable cost-benefit pool of virtual machines, we evaluated the performance of the docking workflow invocations in different configurations of Azure instances. Second, we validated the gains obtained by the proposed method based on the quality of the Reduced Fully-Flexible Receptor (RFFR) models produced using AutoDock4.2. The analyses show that the proposed method reduced the model size by approximately 50% while covering at least 86% of the best docking results from the 74 ligands tested. Third, we tested our novel method using AutoDock Vina, a different docking software, and showed the positive accuracy achieved in the resulting RFFR models. Finally, our results demonstrated that the method proposed optimized ensemble docking experiments and is applicable to different docking software. In addition, it detected new binding modes, which would be unreachable if employing only the rigid structure used to generate the InhA FFR model. CONCLUSIONS:Our results showed that the selective method is a valuable strategy for optimizing ensemble docking-based experiments using different docking software. The RFFR models produced by discarding non-promising snapshots from the original model are accurately shaped for a larger number of ligands, and the elapsed time spent in the ensemble docking experiments are considerably reduced.
| S-EPMC6013854 | biostudies-literature
Project description:Recent experiments with undersaturated aqueous glycine solutions have repeatedly exhibited the presence of giant liquid-like clusters or nanodroplets around 100 nm in diameter. These nanodroplets re-appear even after careful efforts for their removal and purification of the glycine solution. The composition of these clusters is not clear, although it has been suggested that they are mainly composed of glycine, a small and very soluble amino acid. To gain insights into this phenomenon, we study the aggregation of glycine in aqueous solutions at concentrations below the experimental solubility limit using large-scale molecular dynamics simulations under ambient conditions. Three protonation states of glycine (zwitterion = GLZ, anion = GLA, and cation = GLC) are simulated using molecular force fields based on the 1.14*CM1A partial charge scheme, which incorporates the OPLS all-atom force field and TIP3P water. When initiated from dispersed states, we find that giant clusters do not form in our simulations unless salt impurities are present. Moreover, if simulations are initiated from giant cluster states, we find that they tend to dissolve in the absence of salt impurities. Therefore, the simulation results provide little support for the possibility that the giant clusters seen in experiments are composed purely of glycine (and water). Considering that strenuous efforts are made in experiments to remove impurities such as salt, we propose that the giant clusters observed might instead result from the aggregation of reaction products of aqueous glycine, such as diketopiperazine or other oligoglycines which may be difficult to separate from glycine using conventional methods, or their co-aggregation with glycine.
| S-EPMC9251761 | biostudies-literature
Project description:Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.
| S-EPMC8041884 | biostudies-literature
Project description:The accurate prediction of a ligand-protein complex structure is important for computer-assisted drug development. Although many docking methods have been developed over the last three decades, the success of binding structure prediction remains greatly limited. The purpose of this study was to demonstrate the usefulness of molecular dynamics (MD) simulation in assessing a docking pose predicted using a docking program. If the predicted pose is not unstable in an aqueous environment, MD simulation equilibrates the system and removes the ligand from the predicted position. Here we investigated two proteins that are important potential therapeutic targets: ?2 adrenergic receptor (?2AR) and PR-Set7. While ?2AR is rigid and its ligands are very similar to the template ligand (carazolol), PR-Set7 is very flexible and its ligands vary greatly from the template ligand (histone H4 tail peptide). On an empirical basis, we usually expect that the docking prediction is accurate when the protein is rigid and its ligands are similar to the template ligand. The MD analyses in this study clearly suggested such a tendency. Furthermore, we discuss the possibility that the MD simulation can predict the binding pose of a ligand.
| S-EPMC5042163 | biostudies-literature
Project description:The quantum quasi-docking procedure is used to compare the docking accuracies of two quantum-chemical semiempirical methods, namely, PM6-D3H4X and PM7. Quantum quasi-docking is an approximation to quantum docking. In quantum docking, it is necessary to search directly for the global minimum of the energy of the protein-ligand complex calculated by the quantum-chemical method. In quantum quasi-docking, firstly, we look for a wide spectrum of low-energy minima, calculated using the MMFF94 force field, and secondly, we recalculate the energies of all these minima using the quantum-chemical method, and among these recalculated energies we determine the lowest energy and the corresponding ligand position. Both PM6-D3H4X and PM7 are novel methods that describe well-dispersion interactions, hydrogen and halogen bonds. The PM6-D3H4X and PM7 methods are used with the COSMO implicit solvent model as it is implemented in the MOPAC program. The comparison is made for 25 high quality protein-ligand complexes. Firstly, the docking positioning accuracies have been compared, and we demonstrated that PM7+COSMO provides better positioning accuracy than PM6-D3H4X. Secondly, we found that PM7+COSMO demonstrates a much higher correlation between the calculated and measured protein-ligand binding enthalpies than PM6-D3H4X. For future quantum docking PM7+COSMO is preferable, but the COSMO model must be improved.
| S-EPMC8781293 | biostudies-literature
Project description:Incorporating receptor flexibility into molecular docking should improve results for flexible proteins. However, the incorporation of explicit all-atom flexibility with molecular dynamics for the entire protein chain may also introduce significant error and "noise" that could decrease docking accuracy and deteriorate the ability of a scoring function to rank native-like poses. We address this apparent paradox by comparing the success of several flexible receptor models in cross-docking and multiple receptor ensemble docking for p38? mitogen-activated protein (MAP) kinase. Explicit all-atom receptor flexibility has been incorporated into a CHARMM-based molecular docking method (CDOCKER) using both molecular dynamics (MD) and torsion angle molecular dynamics (TAMD) for the refinement of predicted protein-ligand binding geometries. These flexible receptor models have been evaluated, and the accuracy and efficiency of TAMD sampling is directly compared to MD sampling. Several flexible receptor models are compared, encompassing flexible side chains, flexible loops, multiple flexible backbone segments, and treatment of the entire chain as flexible. We find that although including side chain and some backbone flexibility is required for improved docking accuracy as expected, docking accuracy also diminishes as additional and unnecessary receptor flexibility is included into the conformational search space. Ensemble docking results demonstrate that including protein flexibility leads to to improved agreement with binding data for 227 active compounds. This comparison also demonstrates that a flexible receptor model enriches high affinity compound identification without significantly increasing the number of false positives from low affinity compounds.
| S-EPMC2772076 | biostudies-literature
Project description:Amidase and laccase play a key role in the degradation process of anionic polyacrylamide (HPAM). However, the largest challenge of HPAM enzymatic degradation is whether the enzyme can bind with a substrate for a period of time. Here, the most suitable complexes, namely, Rh Amidase-HPAM-2 and Bacillus subtilis (B. subtilis) laccase-HPAM-3, were obtained by docking, and they were carried out for molecular dynamics simulation (MDS) under 298, 303, and 308 K. MDS result analysis showed that Rh Amidase-HPAM-2 was the most stable at 298 K mainly due to a salt bridge and a hydrogen bond, and B. subtilis laccase-HPAM-3 was the most stable at 298 K mainly due to two electrostatic and hydrogen bonds. The LYS96 in Rh Amidase-HPAM-2 and LYS135 in B. subtilis laccase-HPAM-3 had been the most important in their binding process. The binding of Rh Amidase-HPAM-2 and B. subtilis laccase-HPAM-3 was optimal at 303 and 298 K, respectively. HPAM was degraded by mixed bacteria, and the optimal conditions were determined to be 308 K, initial pH = 7, and an inoculated dosage of 2 mL. Under these conditions, the degradation ratio reached 39.24%. The effect of parameters on the HPAM degradation ratio followed a decreasing order of temperature > initial pH > inoculated dosage. The HPAM codegradation mechanism was supposed by mixed bacteria according to test data. The mixed bacteria secreted both amidase and laccase, and they interacted jointly with HPAM. These results lay a theoretical foundation to design and modify the enzyme through mutation experiments in the future.
| S-EPMC10034835 | biostudies-literature
Project description:The phytoconstituents of the aqueous extract from Syzygium jambos L. (Alston) leaves were defined using HPLC-PDA-MS/MS and the antioxidant, anti-aging, antibacterial, and anti-biofilm activities of the extract were in silico and in vitro investigated. The antioxidant activities were performed using in vitro DPPH and FRAP assays as well as H2-DCFDA assay in HaCaT cells in which oxidative stress was induced by UVA radiation. Anti-aging activity was tested in vitro, using aging-related enzymes. The antibacterial, anti-biofilm and inhibitory effects on bacterial mobilities (swarming and swimming) were assessed against Pseudomonas aeruginosa. Results showed that S. jambos aqueous extract contained 28 phytochemicals belonging to different metabolite classes, mainly phenolic acids, gallic acid derivatives, flavonoids, and ellagitannins. Mineral content analysis showed that S. jambos leaves contained moderate amounts of nitrogen, potassium, manganese, magnesium, and zinc, relatively low amounts of phosphorus and copper, and high concentration of calcium and iron. The extract displayed strong antioxidant activities in vitro and inhibited UVA-induced oxidative stress in HaCaT cells. Docking the major compounds identified in the extract into the four main protein targets involved in skin aging revealed an appreciable inhibitory potential of these compounds against tyrosinase, elastase, hyaluronidase, and collagenase enzymes. Moreover, molecular dynamic simulations were adopted to confirm the binding affinity of some selected compounds towards the target enzymes. The extract exhibited pronounced in vitro anti-aging effects, compared to kojic acid and quercetin (the reference compounds). It also inhibited the growth of P. aeruginosa, counteracted its ability to form biofilm, and impeded its swarming and swimming mobilities. Altogether, these findings strongly propose S. jambos leaves as a promising source of bioactive metabolites for the development of natural cosmeceutical and dermatological agents.
| S-EPMC10797028 | biostudies-literature