Combining classical molecular docking with self-consistent charge density-functional tight-binding computations for the efficient and quality prediction of ligand binding structure

To improve the successful prediction rate of the existing molecular docking methods, a new docking approach is proposed that consists of three steps: generating an ensemble of docked poses with a conventional docking method, performing clustering analysis of the ensemble to select the representative poses, and optimizing the representative structures with a low-cost quantum mechanics method. Three quantum mechanics methods, self-consistent charge density-functional tight-binding, ONIOM(DFT:PM6), and ONIOM(SCC-DFTB:PM6), are tested on 18 ligand-receptor bio-complexes. The rate of successful binding pose predictions by the proposed self-consistent charge density-functional tight-binding docking method is the highest, at 67%. The self-consistent charge density-functional tight-binding docking method should be useful for the structure-based drug design.

Intermolecular π/π and H/π interactions in dimers researched by different computational methods

The analysis of π/π and H /π interactions in complexes are a challenging aspect of theoretical research. Due to the different approximations of different levels of theory, results tend to be inconsistent. We compared the reliabilities of HF, SVWN, M06L, PW91, BLYP, B3LYP, BHandHLYP, B97D, MP2, and DFTB-D approaches in researching π/π and H /π interactions by calculating the binding energies of five benzene-containing dimers. The effects of 6-31+G**, 6-311++G** and 6-311++G(2df,2p) basis sets on the results were analyzed too. We found that the DFTB-D and B97D methods combined with the 6-311++G** basis set perform well for dimers that contain π/π and H /π interactions. With high efficiency and satisfactory precision, DFTB-D is helpful for the calculation of complexes containing π/π and H /π stacking. We further calculated the structures and properties of phenylalanine-containing dimers using the DFTB-D and B97D methods. The properties of low energy conformers such as rotational constants, dipole moments and molecular orbitals were also analyzed. These data should be helpful for research into systems that contain π/π and H /π stacking.

Synthesis, antibacterial activities, and theoretical studies of dicoumarols

In this work, we report the synthesis, evaluation of antibacterial activity and theoretical study of dicoumarols including DC, 2-PyDC, 3-PyDC and 4-PyDC.

UNDERSTANDING THE MOLECULAR MECHANISM OF BINDING MODES OF AURORA A INHIBITORS BY LONG TIME SCALE GPU DYNAMICS

Inhibition of Aurora A kinase interaction is considered to be a promising approach for the discovery of new molecularly targeted cancer therapeutics. In this study, the binding mechanisms of two different inhibitors with a contrasting binding affinity to Aurora A were investigated by long time scale GPU molecular dynamics (MD) simulations coupled with molecular mechanics-Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA) method. The results showed that the predicted binding free energies of these two complexes were consistent with the experimental data. Through analyzing the individual energy components of binding free energy, we found that the van der Waals contribution was the main force to drive the inhibitor–protein binding and the electrostatic contribution was also a crucial factor for the inhibitor–Aurora A binding. The structural analysis demonstrated that the inhibitor HPM could produce more hydrophobic interaction contacts with Aurora A than that of 2JZ, and the loss of key hydrogen bonds between the inhibitor and residue Arg137 in the hinge region of Aurora A was another important reason for the weaker binding affinity of 2JZ to Aurora A. This study sheds more light on the development of the efficient inhibitors targeting the Aurora A.

MOLECULAR SIMULATION OF α-TOCOPHEROL PASSING ACROSS DPPC LIPID USING POTENTIAL OF MEAN FORCE AND ACCELERATED MOLECULAR DYNAMICS METHOD

In this paper the process of α-tocopherol (TCP) passing across DPPC membrane was simulated using both the potential of mean force (PMF) and the accelerated molecular dynamics (aMD) methods, respectively. Energy properties, hydrogen bonds and orientation have been compared between these two methods and several conclusions have been obtained. The results indicate that TCP tends to stay at z = 1.2 nm of lipid bilayer. The binding free energy profiles of these two methods are alike. All these show that aMD could obtain comparable results as PMF method, while needs less computation time and resources. Therefore, aMD method could be used as an alternative method for prediction of transport properties of drug-lipid system.

DFT STUDY OF THE PROTONATION AND DEPROTONATION ENTHALPIES OF BENZOXAZOLE, 1,2-BENZISOXAZOLE AND 2,1-BENZISOXAZOLE AND IMPLICATIONS FOR THE STRUCTURES AND ENERGIES OF THEIR ADDUCTS WITH EXPLICIT WATER MOLECULES

Benzoxazole, 1,2-benzisoxazole and 2,1-benzisoxazole are biologically active molecules with potential applications in drug design. Their interaction with aqueous medium in biological systems may be simulated by considering their interaction with explicit water molecules. Such studies provide information on the structures, energies and type of interactions stabilizing the resulting geometric systems. The objective of the current study was to utilize theoretical approaches to investigate the structures, stabilization energy and binding energy of benzoxazole–water, 1,2-benzisoxazole–water and 2,1-benzisoxazole–water complexes. The calculations were performed utilizing the density functional theory (DFT)/M06-2X/6-311 ++ G(d,p) method and the DFT/ωB97XD method with both the 6-311 ++ G(d,p) and the aug-cc-pVDZ basis sets. The results suggest that the stability of the different clusters depends on interrelated factors including the rings formed by intermolecular hydrogen bonds and the proton affinity (PA) or acidity of the atoms forming the intermolecular hydrogen bonds with the water molecules. A comparison across methods indicates that the results follow similar trends with different methods.

FORMIC ACID DEHYDROGENATION ON SURFACES — A REVIEW OF COMPUTATIONAL ASPECT

In this review, we have mainly shown the recent computational studies on formic acid adsorption and selective dissociation to produce hydrogen ( HCOOH → CO 2 + H 2) on several metal ( Pt , Pd , Ni , Cu , Rh and Au ) and metal oxide ( TiO 2, MgO , ZnO and NiO ) surfaces, and both thermal decomposition and electro-catalytic oxidation have been discussed. The decomposition mechanisms of formic acid have been studied by using different computational models and methods, not only interesting and exciting but also different and controversial results have been reported. It is noted that the model systems used in these studies are too simple and idealized, and they cannot represent the real catalysts or the catalytic systems, and more sophisticated computational methodologies and real model systems under the consideration of the working conditions are therefore needed.