Transferable scoring function based on semiempirical quantum mechanical PM6-DH2 method: CDK2 with 15 structurally diverse inhibitors

A semiempirical quantum mechanical PM6-DH2 method accurately covering the dispersion interaction and H-bonding was used to score fifteen structurally diverse CDK2 inhibitors. The geometries of all the complexes were taken from the X-ray structures and were reoptimised by the PM6-DH2 method in continuum water. The total scoring function was constructed as an estimate of the binding free energy, i.e., as a sum of the interaction enthalpy, interaction entropy and the corrections for the inhibitor desolvation and deformation energies. The applied scoring function contains a clear thermodynamical terms and does not involve any adjustable empirical parameter. The best correlations with the experimental inhibition constants (ln K _i) were found for bare interaction enthalpy (r ² = 0.87) and interaction enthalpy corrected for ligand desolvation and deformation energies (r ² = 0.77); when the entropic term was considered, however, the correlation becomes worse but still acceptable (r ² = 0.52). The resulting correlation based on the PM6-DH2 scoring function is better than previously published function based on various docking/scoring, SAR studies or advanced QM/MM approach, however, the robustness is limited by number of available experimental data used in the correlation. Since a very similar correlation between the experimental and theoretical results was found also for a different system of the HIV-1 protease, the suggested scoring function based on the PM6-DH2 method seems to be applicable in drug design, even if diverse protein–ligand complexes have to be ranked.     

Interaction and protection mechanism between Li@C(60) and nucleic acid bases (NABs): performance of PM6-DH2 on noncovalent interaction of NABs-Li@C60

To discuss the protection mechanism of DNA from radiation as well as assess the performance of PM6-DH2 on noncovalent interactions, the interaction of four nucleic acid bases (NABs) such as adenine (A), cytosine (C), guanine (G), and thymine (T), with Li@C(60) was extensively investigated with the-state-of-art theoretical methods describing noncovalent systems, like M06-2x, PBE-D, and PM6-DH2 methods. In the gas phase, the binding strength of NABs to Li@C(60) from M06-2x decreases in the sequence, G>C>A>T. As dispersion was explicitly included, PBE-D relatively enhances the binding of A and T and corrects the sequence to, G>A>C～T. PM6-DH2 predicted similar binding energies to those from PBE-D within 0.5 kcal/mol and the same binding sequence, suggesting that the PM6-DH2 method is promising for nano-scale systems. In the aqueous solution, binding of NABs-Li@C(60) is considerably decreased, and the M06-2X and PM6-D methods yield a different sequence from the gas phase, G>A>T>C. The encapsulation of Li atom results in a lower IP for Li@C(60) than those of NABs, and the dominant localization of single-occupied molecular orbital on Li@C(60) moiety of the complexes NABs-Li@C(60) further indicates that an electron would be ejected from Li@C(60) upon radiation and Li@C(60) is therefore able to protect DNA bases from radiation. In addition, it was revealed that Li prefers coordination with the hexagonal ring at Li@C(60) , which clarifies the existing controversy in this respect. Finally, Yang's reduced density gradient approach clearly shows that the weak and strong noncovalent interaction regions in the complexes, NABs-Li@C(60) and (NABs-Li@C(60) )(+).     

Multiple peptide conformations give rise to similar binding affinities: molecular simulations of p53-MDM2

Molecular dynamics simulations, guided by experimental information (Zondlo et al. Biochemistry 2006, 45, 11945-11957) have been used successfully to reproduce experimental trends in binding affinities of variant p53 peptides with MDM2. Simulations reveal how the conformations of the peptides and the receptor modulate each other to optimize interactions. The conformations of the uncomplexed peptides are governed by a combination of helix and intrinsic disorder (in agreement with experiments), while in the complexed state two very different conformations can coexist. This yields very similar binding affinities, driven by either enthalpy or entropy.     

Solution-state conformations of 2-fluoro-, 2-chloro- and 2-bromo-acetophenone: a dipole moment and kerr effect study

Experimental dipole moments and infinite-dilution Kerr constants for 2-fluoro-, 2-chloro and 2-bromo-acetophenone (CH₃COC₆H₄X; X = F, Cl, Br) as solutes in CCl₄ at 25°C are analysed, yielding the following effective dihedral angles and percentage abundances of the less stable XO-cis conformers: X = F, 10 ± 10°, 5 ± 5%; X = Cl, 40 ± 5°, 10 ± 5%; and X = Br, 65 ± 10°, 25 ± 10%.     

Design,synthesis and computational approach to study novel pyrrole scaffolds as active inhibitors of enoyl ACP reductase (InhA) and Mycobacterium tuberculosis antagonists

Novel 54 pyrrolyl acetohydrazide analogues were designed, synthesized and screened for antitubercular activity against InhA. Enoyl-ACP reductase/InhA one of the significant enzymes implicated in type II FAS (fatty acid synthase) biosynthetic pathway of bacterial outer cell membrane, in addition Mycobacterium tuberculosis H37Rv inhibition potency has proven to be one of the most promising drug target used for designing and testing against TB. In silico molecular modeling was achieved using Surflex-docking method to recognize important binding sites of the enoyl-ACP reductase. 3D-QSAR studies like CoMFA and CoMSIA approaches were studied to create 3D-QSAR depictions for InhA inhibitors. Based on docking results the synthesized molecules are oriented towards the core of the active site. The pyrolyl acetohydrazides exhibited one or two H-bonding connections with InhA enzyme. Molecule 8l (MIC 0.4 μg/mL, InhA- 70% at 50 μM) showed H-bonding connections with Tyr158 and NAD⁺ in a similar mode to that of ligand pyrrolidine carboxamide. The tested molecules also showed good antibacterial (Escherichia coli, Staphylococcus aureus) activity (MIC 0.4–25 μg/mL), while the study against A549 lung adenocarcinoma cell line confirms nontoxic nature of the reported molecules. The QSAR model from CoMFA and CoMSIA through the database configuration showed the most excellent data. The prognostic ability of CoMFA and CoMSIA representations was developed by a test set of 15 molecules that produced cross validated correlation coefficients (q²) of 0.642 and 0.701, respectively. This study gives insight into structural requirement needed for the development of more active InhA inhibitors through in silico approach.     

Opioid conformations: The remarkable effect of C-14 substituents on ring c conformations of 6α and 6β epimers of 6-substituted epoxymorphinan opioid ligands

High field ¹H nmr studies have shown that the nature of the C -14 substituents has a remarkable influence on ring C conformation of epoxymorphinan opioids which have a 6α-hydroxyl group. In sharp contrast to the ring C twist-boat conformation observed in those 6α-hydroxy compounds which also have a 14-hydroxyl group (α-naltrexol; α-oxymorphol), ring C exists predominantly as a chair conformer in 6α-hydroxy compounds which have a proton bound to C -14 (dihydromorphine; dihydrocodeine). The 6β-hydroxy compounds (β-naltrexol; β-oxymorphol) have a ring C chair conformation in agreement with earlier studies.     

The use of multidimensional Franck-Condon simulations to assess model chemistries: a case study on phenol

Multidimensional Franck-Condon simulations of the dispersed fluorescence spectra of phenol generated with geometries obtained from the highly correlated post-Hartree-Fock methods CASSCF, MRCI, and SACCI are presented. While the simulations based on CASSCF and MRCI optimized geometries are very similar to each other and fail to reproduce the experimentally measured intensities faithfully, the simulations obtained from SACCI optimized geometries are very close to the experimental spectra. The code developed for the multidimensional Franck-Condon simulations is described. It is shown that the integral storage problem common to the evaluation of multidimensional Franck-Condon integrals can be overcome by saving all quantities needed to disk. This strategy allows the code to run on computers with limited resources and is very well suited for application to molecules with a very large number of vibrational modes.     

Binding of phosphinate and phosphonate inhibitors to aspartic proteases: a first-principles study

Phosphinate and phosphonate derivatives are potent inhibitors of aspartic proteases (APs). The affinity for the enzyme might be caused by the presence of low barrier hydrogen bonds between the ligand and the catalytic Asp dyad in the cleavage site. We have used density functional theory calculations along with hybrid quantum mechanics/molecular mechanics Car-Parrinello molecular dynamics simulations to investigate the hydrogen-bonding pattern at the binding site of the complexes of human immunodeficiency virus type-1 AP and the eukaryotic endothiapepsin and penicillopepsin. Our calculations are in fair agreement with the NMR data available for endothiapepsin (Coates et al. J. Mol. Biol. 2002, 318, 1405-1415) and show that the most stable active site configuration is the diprotonated, negatively charged form. In the viral complex both protons are located at the catalytic Asp dyad, while in the eukaryotic complexes the proton shared by the closest oxygen atoms is located at the phosphinic/phosphonic group.     

CO2 diffusivity in LiY and NaY faujasite systems: a combination of molecular dynamics simulations and quasi-elastic neutron scattering experiments

Quasi-elastic Neutron Scattering combined with Molecular Dynamics simulations have been carried out to gain further insight into the CO₂ dynamics in LiY and NaY Faujasites. In both materials, it was pointed out that the transport diffusivity (D_T) increases with the loading whereas the self diffusivity (D_S) decreases. In addition, it was shown that LiY exhibits a significant slower CO₂ self diffusivity process due to a strong interaction between the Li⁺ cation and the adsorbate molecules at the initial stage of diffusion. This result is consistent with higher simulated activation energy in this cation exchanged faujasite form. By contrast, the transport diffusivity is revealed to be slightly faster in LiY than in NaY.     

The NMR spectra and conformations of dihydropyran derivatives—I: The conformations of 2-alkoxy-5,6-dihydro-α-pyran-6-carboxylic esters

On the basis of NMR spectra of trans- and cis-2-alkoxy-5,6-α-pyran-6-carboxylic esters it was found that, at room temperature, the trans-compounds exist exclusively in the conformation with equatorial carbalkoxy and pseudoaxial alkoxy groups. The cis-isomers appear to be in conformational equilibrium between a form with equatorial carbalkoxy and pseudoequatorial alkoxy groups and that with axial and pseudoaxial substituents. In the latter case, the axial carbalkoxy group is bent out off its normal position by about 15°.     

Investigation of binding features: Effects on the interaction between CYP2A6 and inhibitors

A computational investigation has been carried out on CYP2A6 and its naphthalene inhibitors to explore the crucial molecular features contributing to binding specificity. The molecular bioactive orientations were obtained by docking (FlexX) these compounds into the active site of the enzyme. And the density functional theory method was further used to optimize the molecular structures with the subsequent analysis of molecular lipophilic potential (MLP) and molecular electrostatic potential (MEP). The minimal MLPs, minimal MEPs, and the band gap energies (the energy difference between the highest occupied molecular orbital and lowest unoccupied molecular orbital) showed high correlations with the inhibition activities (pIC₅₀s), illustrating their significant roles in driving the inhibitor to adopt an appropriate bioactive conformation oriented in the active site of CYP2A6 enzyme. The differences in MLPs, MEPs, and the orbital energies have been identified as key features in determining the binding specificity of this series of compounds to CYP2A6 and the consequent inhibitory effects. In addition, the combinational use of the docking, MLP and MEP analysis is also demonstrated as a good attempt to gain an insight into the interaction between CYP2A6 and its inhibitors. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Eclipsed and staggered conformations of (SIH3)2F+: An ab initio study

Ab initio calculations at 6–31G**, 6–31++G**, and MP2/6–31G** levels were performed on disilyl–fluoronium, (SiH₃)₂F⁺, with the SiH₃ group eclipsed or staggered. Optimized geometries, total energies, dipole moments, atomic charges, electronic density, and vibrational frequencies were computed. The results were compared with calculated structural parameters and vibrational frequencies of H₃SiF, H₂SiF⁺, H₂SiF^?, and H₄SiF⁺ ions. The basis-set effects were studied. Several thermochemistry parameters—ZPE, thermal energy, rotational constants, and entropies—were also calculated. © 1994 John Wiley & Sons, Inc.     

Green corrosion inhibitors of steel based on peptides and their constituents: a combination of experimental and computational approach

Journal of Solid State Electrochemistry - In this work, six groups of non-toxic compounds were tested as inhibitors for mild steel in 1 M HCl solution: glycine (Gly), glutamic acid (Glu),...     

Flexible matching of test ligands to a 3D pharmacophore using a molecular superposition force field: Comparison of predicted and experimental conformations of inhibitors of three enzymes

Summary A computer procedure TFIT, which uses a molecular superposition force field to flexibly match test compounds to a 3D pharmacophore, was evaluated to find out whether it could reliably predict the bioactive conformations of flexible ligands. The program superposition force field optimizes the overlap of those atoms of the test ligand and template that are of similar chemical type, by applying an attractive force between atoms of the test ligand and template which are close together and of similar type (hydrogen bonding, charge, hydrophobicity). A procedure involving Monte Carlo torsion perturbations, followed by torsional energy minimization, is used to find conformations of the test ligand which cominimize the internal energy of the ligand and the superposition energy of ligand and template. The procedure was tested by applying it to a series of flexible ligands for which the bioactive conformation was known experimentally. The 15 molecules tested were inhibitors of thermolysin, HIV-1 protease or endothiapepsin for which X-ray structures of the bioactive conformation were available. For each enzyme, one of the molecules served as a template and the others, after being conformationally randomized, were fitted. The fitted conformation was then compared to the known binding geometry. The matching procedure was successful in predicting the bioactive conformations of many of the structures tested. Significant deviation from experimental results was found only for parts of molecules where it was readily apparent that the template did not contain sufficient information to accurately determine the bioactive conformation.     

Laser spectroscopic study on the conformations and the hydrated structures of benzo-18-crown-6-ether and dibenzo-18-crown-6-ether in supersonic jets

The laser-induced fluorescence spectra of jet-cooled benzo-18-crown-6 (B18C6) and dibenzo-18-crown-6 (DB18C6) exhibit a number of vibronic bands in the 35 000-37 000 cm(-1) region. We attribute these bands to monomers and hydrated clusters by fluorescence-detected IR-UV and UV-UV double resonance spectroscopy. We found four and two conformers for bare B18C6 and DB18C6, and the hydration of one water molecule reduces the number of isomers to three and one for B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1), respectively. The IR-UV spectra of B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1) suggest that all isomers of the monohydrated clusters have a double proton-donor type (bidentate) hydration. That is, the water molecule is bonded to B18C6 or DB18C6 via two O-H[dot dot dot]O hydrogen bonds. The blue shift of the electronic origin of the monohydrated clusters and the quantum chemical calculation suggest that the water molecule in B18C6-(H(2)O)(1) and DB18C6-(H(2)O)(1) prefers to be bonded to the ether oxygen atoms near the benzene ring.     

Molecular conformations of a partially halogenated ether: A study based on infrared spectroscopy and density functional theory calculations

A new partially halogenated ether (ClCF₂CF(CF₃)OCF₂CH₃) has been synthesized and characterized using DSC, GC, ¹H and ¹⁹F NMR, IR. The experimental infrared spectra of this “flexible” molecule have been successfully interpreted on the basis of reliable Density Functional Theory calculations. An efficient method useful for the identification of the many stable conformers has been developed and applied. Infrared spectra of the stable conformers have been simulated after full geometry optimization. The results obtained allow detection of conformation-sensitive bands, making possible the interpretation of fine details in the spectra.     

Analysis of the RGD sequence in protein structures: comparison to the conformations of the RGDW and DRGDW peptides determined by molecular dynamics simulations

Geometric properties of the RGD sequence in a data set of protein crystal and NMR structures deposited in the Protein Data Bank were examined to identify structural characteristics that are related to cell adhesion activity. Interatomic distances and dihedral angles are examined. These geometric measures are then used in an analysis of the conformations of the RGDW and _DRGDW peptides obtained from molecular dynamics simulations (Stote RH, et al. (2000) J Phys ChemB 104:1624). This analysis leads to the suggestion that differences in the accessible conformations contribute to the difference in biological activity between the RGDW and the _DRGDW peptides. Received: 15 August 2000 / Accepted: 4 October 2000 / Published online: 21 March 2001     

Boosted feature selectors: a case study on prediction P-gp inhibitors and substrates

Feature selection is commonly used as a preprocessing step to machine learning for improving learning performance, lowering computational complexity and facilitating model interpretation. This paper proposes the application of boosting feature selection to improve the classification performance of standard feature selection algorithms evaluated for the prediction of P-gp inhibitors and substrates. Two well-known classification algorithms, decision trees and support vector machines, were used to classify the chemical compounds. The experimental results showed better performance for boosting feature selection with respect to the standard feature selection algorithms while maintaining the capability for feature reduction.