Binding free energy calculation for duocarmycin/DNA complex based on the QPLD-derived partial charge model

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

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.

Binding affinities of the farnesoid X receptor in the D3R Grand Challenge 2 estimated by free-energy perturbation and docking

We have studied the binding of 102 ligands to the farnesoid X receptor within the D3R Grand Challenge 2016 blind-prediction competition. First, we employed docking with five different docking software and scoring functions. The selected docked poses gave an average root-mean-squared deviation of 4.2 Å. Consensus scoring gave decent results with a Kendall's τ of 0.26 ± 0.06 and a Spearman's ρ of 0.41 ± 0.08. For a subset of 33 ligands, we calculated relative binding free energies with free-energy perturbation. Five transformations between the ligands involved a change of the net charge and we implemented and benchmarked a semi-analytic correction (Rocklin et al., J Chem Phys 139:184103, 2013) for artifacts caused by the periodic boundary conditions and Ewald summation. The results gave a mean absolute deviation of 7.5 kJ/mol compared to the experimental estimates and a correlation coefficient of R 2 = 0.1. These results were among the four best in this competition out of 22 submissions. The charge corrections were significant (7-8 kJ/mol) and always improved the results. By employing 23 intermediate states in the free-energy perturbation, there was a proper overlap between all states and the precision was 0.1-0.7 kJ/mol. However, thermodynamic cycles indicate that the sampling was insufficient in some of the perturbations.

In silico and in vitro studies of cinnamaldehyde and their derivatives against LuxS in Streptococcus pyogenes: effects on biofilm and virulence genes

The LuxS-based signalling pathway has an important role in physiological and pathogenic functions that are capable of causing different infections. In the present study, cinnamaldehyde (CN) and their derivatives were evaluated for their inhibitory efficiency against LuxS by molecular modelling, docking, dynamics and free-energy calculations. Sequence and structure-similarity analysis of LuxS protein, five different amino acids were found to be highly conserved, of which GLY128 was identified as the key residue involved in the effective binding of the ligands. Quantum-polarized ligand docking protocol showed that 2nitro and 4nitro CN has a higher binding efficiency than CN, which very well corroborates with the in vitro studies. COMSTAT analysis for the microscopic images of the S. pyogenes biofilm showed that the ligands have antibiofilm potential. In addition, the results of quantitative polymerase chain reaction (qPCR) analysis revealed that the transcripts treated with the compounds showed decrease in luxS expression, which directly reflects with the reduction in expression of speB. No substantial effect was observed on the virulence regulator (srv) transcript. These results confirm that speB is controlled by the regulation of luxS. The decreased rate of S. pyogenes survival in the presence of these ligands envisaged the fact that the compounds could readily enhance opsonophagocytosis with the reduction of virulence factor secretion. Thus, the overall data supports the use of CN derivatives against quorum sensing-mediated infections caused by S. pyogenes.

Designing a ligand for pharmaceutical purposes

BACKGROUND

Drug approval applications to the FDA have shown a remarkably small increment compared with what was expected. In the last few years several efforts have been made to improve the results of rational drug design approaches and in particular to predict inhibitor-target structure and to evaluate the free energy of binding. Virtual database screening, combined with other computational methods, is one of the most promising methods to overcome this key issue.

OBJECTIVE

It is possible to understand how computational medicinal chemistry is changing, improving from its errors and moving towards becoming a more important tool for drug development.

METHODS

Some of the most recent modeling techniques have been presented and in particular the benefits of combining these techniques are highlighted.

RESULTS/CONCLUSION

At present computational chemists can understand the peculiar problems associated with the study of biological systems and on this basis they can choose the right collection of complementary in silico approaches to solve the medicinal chemistry problem in a better manner.

Biological evaluation and docking studies of recently identified inhibitors of phosphoinositide-3-kinases

The alpha isoform of the phosphatidylinositol-3-kinases (PI3Kα) is often mutated, amplified and overexpressed in human tumors. In an effort to develop new inhibitors targeting this enzyme, we carried out a pharmacophore model study based on six PI3Kα-selective compounds. The pharmacophore searching identified three structurally novel inhibitors of PI3Kα and its H1047R mutant. Our biological studies show that two of our hit molecules suppressed the formation of pAKT, a downstream effector of PI3Kα, and induced apoptosis in the HCT116 colon cancer cell line. QPLD-based docking showed that residues Asp933, Glu849, Val851, and Gln859 appeared to be key binding residues for active inhibitors.

Molecular modelling methods to quantitate drug-DNA interactions

We describe a molecular modelling method for calculating the binding affinity of ligands for DNA. Though theoretically applicable to any form of noncovalent interaction, we concentrate on the case of predicting the sequence selectivity of a minor-groove binding ligand. The method is based on performing molecular dynamics (MD) simulations on DNA sequences, with and without the ligand bound, and postprocessing the molecular dynamics trajectory data to obtain approximate free energies of binding. We discuss issues relating to the preparation of the structures for simulation, choices for the molecular dynamics simulation method itself, methods for evaluating the reliability and stability of the simulation data, and finally alternative approaches to postprocessing the data to extract approximate free energies of binding.

A theoretical study of repeating sequence in HRP II: a combination of molecular dynamics simulations and (17)O quadrupole coupling tensors

Histidine rich protein II derived peptide (HRP II 169-182) was investigated by molecular dynamics, MD, simulation and (17)O electric field gradient, EFG, tensor calculations. MD simulation was performed in water at 300 K with alpha-helix initial structure. It was found that peptide loses its initial alpha-helix structure rapidly and is converted to random coil and bent secondary structures. To understand how peptide structure affects EFG tensors, initial structure and final conformations resulting from MD simulations were used to calculate (17)O EFG tensors of backbone carbonyl oxygens. Calculations were performed using B3LYP method and 6-31+G basis set. Calculated (17)O EFG tensors were used to evaluate quadrupole coupling constants, QCC, and asymmetry parameters, eta(Q). Difference between the calculated QCC and eta(Q) values revealed how hydrogen-bonding interactions affect EFG tensors at the sites of each oxygen nucleus.