Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock

<p>Abstract</p> <p>Background</p> <p>Molecular docking methods are commonly used for predicting binding modes and energies of ligands to proteins. For accurate complex geometry and binding energy estimation, an appropriate method for calculating partial charges is essential. AutoDockTools software, the interface for preparing input files for one of the most widely used docking programs AutoDock 4, utilizes the Gasteiger partial charge calculation method for both protein and ligand charge calculation. However, it has already been shown that more accurate partial charge calculation - and as a consequence, more accurate docking- can be achieved by using quantum chemical methods. For docking calculations quantum chemical partial charge calculation as a routine was only used for ligands so far. The newly developed Mozyme function of MOPAC2009 allows fast partial charge calculation of proteins by quantum mechanical semi-empirical methods. Thus, in the current study, the effect of semi-empirical quantum-mechanical partial charge calculation on docking accuracy could be investigated.</p> <p>Results</p> <p>The docking accuracy of AutoDock 4 using the original AutoDock scoring function was investigated on a set of 53 protein ligand complexes using Gasteiger and PM6 partial charge calculation methods. This has enabled us to compare the effect of the partial charge calculation method on docking accuracy utilizing AutoDock 4 software. Our results showed that the docking accuracy in regard to complex geometry (docking result defined as accurate when the RMSD of the first rank docking result complex is within 2 Å of the experimentally determined X-ray structure) significantly increased when partial charges of the ligands and proteins were calculated with the semi-empirical PM6 method.</p> <p>Out of the 53 complexes analyzed in the course of our study, the geometry of 42 complexes were accurately calculated using PM6 partial charges, while the use of Gasteiger charges resulted in only 28 accurate geometries. The binding affinity estimation was not influenced by the partial charge calculation method - for more accurate binding affinity prediction development of a new scoring function for AutoDock is needed.</p> <p>Conclusion</p> <p>Our results demonstrate that the accuracy of determination of complex geometry using AutoDock 4 for docking calculation greatly increases with the use of quantum chemical partial charge calculation on both the ligands and proteins.</p

Electron transfer mechanism for periodic acid oxidation of aromatic substrates

Oxidn. of phenolic ethers, e.g. $I (R = Me, R1 = H)$ with 1 equiv of $H_5IO$6 gave the iodinated product $I (R = Me, R1 = iodo)$. Use of 2 equiv of oxidant gave $I (R = CHO, R1 = iodo)$in which the side chain was oxidized. Oxidn. of binaphthyl II with either $H_5IO_6$ or ceric ammonium nitrate in MeOH gave the same cyclized product III. ESR signals were found for the $H_5IO_6$ oxidn. of 11 arom. substrates. The ESR data and addnl. chem. evidence suggest a radical cation mechanism for the oxidation

Structural Crystallography and Crystal Chemistry

CaH406P-.K +, M r = 206.10, is orthorhombic, space group Pbca (from systematic absences), a = 14.538(4), b = 13.364(5), c = 6.880 (6)A, U = 1383.9 A 3, D x = 2.07 Mg m -a, Z = 8, ~.(Mo Ka) = 0.7107/~, p(MO Ka) = 1.015 mm -1. The final R value is 0.042 for a total of 1397 reflections. The high energy P-O(13) and the enolic C(1)-O(13) bonds are 1.612 and 1.374 A respectively. The enolpyruvate moiety is essentially planar. The orientation of the phosphate with respect to the pyruvate group in PEP.K is distinctly different from that in the PEP-cyclohexylammonium salt, the torsion angle C (2)-C (1)-O(13)- P being -209.1 in the former and -90 ° in the latter. The K + ion binds simultaneously to both the phosphate and carboxyl ends of the same PEP molecule. The ester O(13) is also a binding site for the cation. The K + ion is coplanar with the pyruvate moiety and binds to 0(22) and O(13) almost along their lone-pair directions. The carbonyl 0(22) prefers to bind to the K + ion rather than take part in the formation of hydrogen bonds usually observed in carboxylic acid structures

Ultrasonic evaluation of delamination in quasi-isotropic CFRP laminates subjected to low-velocity impact

Ultrasonic C-Scan is used very often to detect flaws and defects in the composite components resulted during fabrication and damages resulting from service conditions. Evaluation and characterization of defects and damages of composites require experience and good understanding of the material as they are distinctly different in composition and behavior as compared to conventional metallic materials. The failure mechanisms in composite materials are quite complex. They involve the interaction of matrix cracking, fiber matrix interface debonding, fiber pullout, fiber fracture and delamination. Generally all of them occur making the stress and failure analysis very complex. Under low-velocity impact loading delamination is observed to be a major failure mode. In composite materials the ultrasonic waves suffer high acoustic attenuation and scattering effect, thus making data interpretation difficult. However these difficulties can be overcome to a greater extent by proper selection of probe, probe parameter settings like pulse width, pulse amplitude, pulse repetition rate, delay, blanking, gain etc., and data processing which includes image processing done on the image obtained by the C-Scan

FAILURE BEHAVIOR OF WOVEN FABRIC COMPOSITES

Woven fabrics are now considered important reinforcing materials for use as unidirectional (UD) tapes in composites technology. The understanding of the failure mechanisms in woven fabric composites, however, is not complete. The tensile behavior of on-axis and off-axis specimens of E-glass/epoxy and carbon/epoxy plain weave fabric composites with and without a hole is studied. The failure behavior of woven fabric composites differs from UD tape composites. Because of the interlaced nature of warp and weft fibers, their reorientation along the loading direction is possible in off-axis specimens. The reorientation of fibers may increase the notched strength over the unnotched strength. The failure modes in angle-ply laminates of UD tape are mostly matrix cracking, fiber pullout, delamination, and to some extent, fiber breakage, whereas in woven fabric composites, the ultimate failure is always due to fiber breakage which, in general, is a preferred failure mode

Compression after impact testing of carbon fiber reinforced plastic laminates

Carbon fiber reinforced plastic (CFRP) composite laminates are susceptible for multiple delaminations, which reduce the compressive strength significantly even when the impact load is insufficient to cause visible damage. The post-impact compression test is used widely to assess the relative performance of different composite laminates with different fiber matrix combinations. In the present work, CFRP (T300/914) laminates were subjected to low-velocity impact loading simulating tool drops and runway debris. Compression after impact (CAI) tests were performed using fixtures as per NASA specifications. Results obtained indicate definite correlation between impact energy, failure strain, and residual compressive strength

The Structure of Monosodium Phosphoenolpyruvate

CDH406P-.Na +.H20 , M r = 208.0, is monoclinic, Cc, a = 11.423 (2), b = 23.253 (5), c - 6.604 (1) A, fl = 123.63 (1) °, U = 1460.6 A 3, D x =. 1.89 Mg m -a, Z = 8, 2(Mo Ka) = 0.7107 A, p(Mo Ka) = 0.44 mm -~, F(000) = 840. Final R = 0.063 for 1697 reflections.The two crystallographically independent molecules of phosphoenolpyruvate (PEP) (A and B) are almost mirror images of each other, the mirror being the planar enolpyruvate group. The torsion angle C(3)-C(2)- O(1)-P(1) is 122.6 in A and -112.0 ° in B, in contrast to -209.1 ° in PEP.K. The enolic C(2)-O(1) has a partial double-bond character [1.401 (A), 1.386A (B)]. The high-energy P~O bond (1.595 and 1.610A) is comparable to that in PEP.K (1.612 A). Na(1) has six nearest neighbours while Na(2) has only five. The Na + ions are involved in binding only the phosphates of different molecules, in contrast to the K ÷ ion in PEP. K, which binds to both the phosphate and carboxyl ends of the same molecule. The planar carboxyl groups stack on each other at an average distance of 3.2 A instead of forming hydrogen-bonded dimers usually found in carboxylate structures

Rapid screening for HIV-1 protease inhibitor leads through X-ray diffraction

Knowledge of the three-dimensional structures of HIV-1 protease and of its complexes with various inhibitors has played a key role in development of drugs against AIDS. Hexagonal crystals of unliganded tethered HIV-1 protease in which the enzyme conformation is identical to its ligand-bound state can be used in combination with the soaking method in order to identify potential inhibitor leads via X-ray diffraction. The advantages of the soaking method are the generality of application and the rapidity of structure determination for iterative structure-based drug design. Structures of two ligand complexes with HIV-1 protease determined using this method are shown to be very similar to the structures obtained earlier via co-crystallization