Browse
Submit Data
Databases
API
Help

Dataset Information

31 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Quantum-Chemical Quasi-Docking for Molecular Dynamics Calculations.

ABSTRACT: 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.

SUBMITTER: Sulimov A

PROVIDER: S-EPMC8781293 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Understanding the Phosphorylation Mechanism by Using Quantum Chemical Calculations and Molecular Dynamics Simulations.

Project description:Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (Egap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design.

| S-EPMC6138447 | biostudies-literature

Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations.

Project description:Haspin, an atypical serine/threonine protein kinase, is a potential target for cancer therapy. 5-iodotubercidin (5-iTU), an adenosine derivative, has been identified as a potent Haspin inhibitor in vitro. In this paper, quantum chemical calculations and molecular dynamics (MD) simulations were employed to identify and quantitatively confirm the presence of halogen bonding (XB), specifically halogen∙∙∙π (aromatic) interaction between halogenated tubercidin ligands with Haspin. Consistent with previous theoretical finding, the site specificity of the XB binding over the ortho-carbon is identified in all cases. A systematic increase of the interaction energy down Group 17, based on both quantum chemical and MD results, supports the important role of halogen bonding in this series of inhibitors. The observed trend is consistent with the experimental observation of the trend of activity within the halogenated tubercidin ligands (F < Cl < Br < I). Furthermore, non-covalent interaction (NCI) plots show that cooperative non-covalent interactions, namely, hydrogen and halogen bonds, contribute to the binding of tubercidin ligands toward Haspin. The understanding of the role of halogen bonding interaction in the ligand-protein complexes may shed light on rational design of potent ligands in the future.

| S-EPMC8840108 | biostudies-literature

Impact of Doxorubicin on Self-Organization of Congo Red: Quantum Chemical Calculations and Molecular Dynamics Simulations.

Project description:Quantum-chemical calculations and molecular dynamics simulation were applied to a model self-organization process of Congo red (CR) molecules in aqueous solution and the impact of doxorubicin (DOX) molecules on such a process. It was demonstrated that both pure CR/CR and mixed CR/DOX dimers were stable. Van der Waals interactions between aromatic units were responsible for a stacked dimer formation. An important source of stabilization in the CR/CR dimer was the polarization energy. In the CR/DOX mixed dimer long range, electrostatic interactions were the main driving force leading to complexation. An implicit solvent model showed that the formation of the CR/CR dimer was favored over the CR/DOX one. Molecular dynamics simulations demonstrated rapid complexation. In the pure CR system, short sequences of ribbon-like structures were formed. Such structures might be glued by hydrogen bonds to form bigger complexes. It was shown that the aromatic part of the DOX molecule enters CR ribbons with the sugar part covering the CR ribbons. These findings demonstrated that CR may find applications as a carrier in delivering DOX molecules; however, further more extensive investigations are required.

| S-EPMC7424579 | biostudies-literature

Stretched DNA investigated using molecular-dynamics and quantum-mechanical calculations.

Project description:We combined atomistic molecular-dynamics simulations with quantum-mechanical calculations to investigate the sequence dependence of the stretching behavior of duplex DNA. Our combined quantum-mechanical/molecular-mechanical approach demonstrates that molecular-mechanical force fields are able to describe both the backbone and base-base interactions within the highly distorted nucleic acid structures produced by stretching the DNA from the 5' ends, which include conformations containing disassociated basepairs, just as well as these force fields describe relaxed DNA conformations. The molecular-dynamics simulations indicate that the force-induced melting pathway is sequence-dependent and is influenced by the availability of noncanonical hydrogen-bond interactions that can assist the disassociation of the DNA basepairs. The biological implications of these results are discussed.

| S-EPMC2800961 | biostudies-literature

Using molecular dynamics and quantum mechanics calculations to model fluorescence observables.

Project description:We provide a critical examination of two different methods for generating a donor-acceptor electronic coupling trajectory from a molecular dynamics (MD) trajectory and three methods for sampling that coupling trajectory, allowing the modeling of experimental observables directly from the MD simulation. In the first coupling method we perform a single quantum-mechanical (QM) calculation to characterize the excited state behavior, specifically the transition dipole moment, of the fluorescent probe, which is then mapped onto the configuration space sampled by MD. We then utilize these transition dipoles within the ideal dipole approximation (IDA) to determine the electronic coupling between the probes that mediates the transfer of energy. In the second method we perform a QM calculation on each snapshot and use the complete transition densities to calculate the electronic coupling without need for the IDA. The resulting coupling trajectories are then sampled using three methods ranging from an independent sampling of each trajectory point (the independent snapshot method) to a Markov chain treatment that accounts for the dynamics of the coupling in determining effective rates. The results show that the IDA significantly overestimates the energy transfer rate (by a factor of 2.6) during the portions of the trajectory in which the probes are close to each other. Comparison of the sampling methods shows that the Markov chain approach yields more realistic observables at both high and low FRET efficiencies. Differences between the three sampling methods are discussed in terms of the different mechanisms for averaging over structural dynamics in the system. Convergence of the Markov chain method is carefully examined. Together, the methods for estimating coupling and for sampling the coupling provide a mechanism for directly connecting the structural dynamics modeled by MD with fluorescence observables determined through FRET experiments.

| S-EPMC3087818 | biostudies-literature

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations.

Project description:Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

| S-EPMC8431471 | biostudies-literature

Mechanism of titin unfolding by force: insight from quasi-equilibrium molecular dynamics calculations.

Project description:We have studied the unfolding by force of one of the immunoglobulin domains of the muscle protein titin using molecular dynamics simulations at 300 K. Previous studies, done at constant pulling rates, showed that under the effect of the force two strands connected to each other by six backbone H-bonds are pulled apart. No details about the mechanism of H-bond breaking were provided. Our simulation protocol "pull and wait" was designed to correspond to very slow pulling, more similar to the rates used in experiments than are the protocols used in previous computational studies. Under these conditions interstrand backbone H-bonds are not "ripped apart" by the application of the force. Instead, small elongations produced by the force weaken specific backbone H-bonds with respect to water-backbone H-bonds. These weakened bonds allow a single water molecule to make H-bonds to the CO and the NH of the same backbone H-bond while they are still bound to each other. The backbone H-bond then breaks (distance > 3.6 A), but its donor and acceptor atoms remain bound to the same water molecule. Further separation of the chains takes place when a second water molecule makes an H-bond with either the protein backbone donor or acceptor atom. Thus, the force does not directly break the main chain H-bonds: it destabilizes them in such a way that they are replaced by H-bonds to water. With this mechanism, the force necessary to break all the H-bonds required to separate the two strands will be strongly dependent on the pulling speed. Further simulations carried out at low forces but long waiting times (> or = 500 ps, < or = 10 ns) show that, given enough time, even a very small pulling force (< 400 pN) is sufficient to destabilize the interstrand H-bonds and allow them to be replaced by H-bonds to two water molecules. As expected, increasing the temperature to 350 K allows the interstrand H-bonds to break at lower forces than those required at 300 K.

| S-EPMC1483073 | biostudies-literature

Quantum chemical studies on molecular structure, AIM, ELF, RDG and antiviral activities of hybrid hydroxychloroquine in the treatment of COVID-19: Molecular docking and DFT calculations.

Project description:Structure-activity relationships for hydroxychloroquine compound and its derivatives resulted in a potent antiviral activity. Where hydroxychloroquine derivatives showed an apparent efficacy against coronavirus related pneumonia. For this reason, the current study is focused on the structural properties of hydroxychloroquine and hydroxychloroquine sulfate. Optimized structures of these molecules have been reported by using DFT method at B3LYP/6-31G* level of theory. The geometric were determined and compared with the experimental crystal structure. The intra and intermolecular interactions which exist within these compounds are analyzed by different methods namely the topological analysis AIM, ELF and the reduced gradient of the density. These approaches make it possible in particular to study the properties of hydrogen bonds. The highest occupied molecular orbital and the lowest unoccupied molecular orbital energy levels are constructed and the corresponding frontier energy gaps are determined to realize the charge transfer within the molecule. The densities of state diagrams were determined to calculate contributions to the molecular orbitals. The molecular electrostatic potential surfaces are determined to give a visual representation of charge distribution of these ligands and to provide information linked to electrophilic and nucleophilic sites localization. Finally, these derivatives were evaluated for the inhibition of COVID-19 activity by using the molecular docking method.

| S-EPMC7787522 | biostudies-literature

A quantum chemical molecular dynamics repository of solvated ions.

Project description:The importance of ion-solvent interactions in predicting specific ion effects in contexts ranging from viral activity through to electrolyte viscosity cannot be underestimated. Moreover, investigations of specific ion effects in nonaqueous systems, highly relevant to battery technologies, biochemical systems and colloid science, are severely limited by data deficiency. Here, we report IonSolvR - a collection of more than 3,000 distinct nanosecond-scale ab initio molecular dynamics simulations of ions in aqueous and non-aqueous solvent environments at varying effective concentrations. Density functional tight binding (DFTB) is used to detail the solvation structure of up to 55 solutes in 28 different protic and aprotic solvents. DFTB is a fast quantum chemical method, and as such enables us to bridge the gap between efficient computational scaling and maintaining accuracy, while using an internally-consistent simulation technique. We validate the database against experimental data and provide guidance for accessing individual IonSolvR records.

| S-EPMC9304403 | biostudies-literature

Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations.

Project description:Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

| S-EPMC9554917 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data