### Journal: Journal of computational chemistry

#### 28

##### Conventional strain energies of azetidine and phosphetane: Can density functional theory yield reliable results?

- Journal of computational chemistry
- Published over 7 years ago
- Discuss

The conventional strain energies for azetidine and phosphetane are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero-point vibrational energies are computed for all pertinent molecular systems using self-consistent field theory, second-order perturbation theory, and density functional theory and using the correlation consistent basis sets cc-pVDZ, cc-pVTZ, and cc-pVQZ. Single point fourth-order perturbation theory, CCSD, and CCSD(T) calculations using the cc-pVTZ and the cc-pVQZ basis sets are computed using the MP2/cc-pVTZ and MP2/cc-pVQZ optimized geometries, respectively, to ascertain the contribution of higher order correlation effects and to determine if the quadruple-zeta valence basis set is needed when higher order correlation is included. In the density functional theory study, eight different functionals are used including B3LYP, wB97XD, and M06-2X to determine if any functional can yield results similar to those obtained at the CCSD(T) level. © 2012 Wiley Periodicals, Inc.

#### 28

##### Binding structures of tri-N-acetyl-β-glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

- Journal of computational chemistry
- Published over 7 years ago
- Discuss

Lysozyme is a well-studied enzyme that hydrolyzes the β-(1,4)-glycosidic linkage of N-acetyl-β-glucosamine (NAG)(n) oligomers. The active site of hen egg-white lysozyme (HEWL) is believed to consist of six subsites, A-F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all-atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)(3) . MD trajectories are applied to analyze structures and conformation of the complex as well as protein-ligand interactions, including the hydrogen-bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)(3) are investigated independently based on a fixed-charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM-GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root-mean-squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed-charge models. © 2012 Wiley Periodicals, Inc.

#### 28

##### An analytical method for computing atomic contact areas in biomolecules

- Journal of computational chemistry
- Published almost 8 years ago
- Discuss

We propose a new analytical method for detecting and computing contacts between atoms in biomolecules. It is based on the alpha shape theory and proceeds in three steps. First, we compute the weighted Delaunay triangulation of the union of spheres representing the molecule. In the second step, the Delaunay complex is filtered to derive the dual complex. Finally, contacts between spheres are collected. In this approach, two atoms i and j are defined to be in contact if their centers are connected by an edge in the dual complex. The contact areas between atom i and its neighbors are computed based on the caps formed by these neighbors on the surface of i; the total area of all these caps is partitioned according to their spherical Laguerre Voronoi diagram on the surface of i. This method is analytical and its implementation in a new program BallContact is fast and robust. We have used BallContact to study contacts in a database of 1551 high resolution protein structures. We show that with this new definition of atomic contacts, we generate realistic representations of the environments of atoms and residues within a protein. In particular, we establish the importance of nonpolar contact areas that complement the information represented by the accessible surface areas. This new method bears similarity to the tessellation methods used to quantify atomic volumes and contacts, with the advantage that it does not require the presence of explicit solvent molecules if the surface of the protein is to be considered. © 2012 Wiley Periodicals, Inc.

#### 28

##### Accelerating VASP electronic structure calculations using graphic processing units

- Journal of computational chemistry
- Published almost 8 years ago
- Discuss

We present a way to improve the performance of the electronic structure Vienna Ab initio Simulation Package (VASP) program. We show that high-performance computers equipped with graphics processing units (GPUs) as accelerators may reduce drastically the computation time when offloading these sections to the graphic chips. The procedure consists of (i) profiling the performance of the code to isolate the time-consuming parts, (ii) rewriting these so that the algorithms become better-suited for the chosen graphic accelerator, and (iii) optimizing memory traffic between the host computer and the GPU accelerator. We chose to accelerate VASP with NVIDIA GPU using CUDA. We compare the GPU and original versions of VASP by evaluating the Davidson and RMM-DIIS algorithms on chemical systems of up to 1100 atoms. In these tests, the total time is reduced by a factor between 3 and 8 when running on n (CPU core + GPU) compared to n CPU cores only, without any accuracy loss. © 2012 Wiley Periodicals, Inc.

#### 28

##### Fast force field-based optimization of protein-ligand complexes with graphics processor

- Journal of computational chemistry
- Published almost 8 years ago
- Discuss

Usually based on molecular mechanics force fields, the post-optimization of ligand poses is typically the most time-consuming step in protein-ligand docking procedures. In return, it bears the potential to overcome the limitations of discretized conformation models. Because of the parallel nature of the problem, recent graphics processing units (GPUs) can be applied to address this dilemma. We present a novel algorithmic approach for parallelizing and thus massively speeding up protein-ligand complex optimizations with GPUs. The method, customized to pose-optimization, performs at least 100 times faster than widely used CPU-based optimization tools. An improvement in Root-Mean-Square Distance (RMSD) compared to the original docking pose of up to 42% can be achieved. © 2012 Wiley Periodicals, Inc.

#### 28

##### C library for topological study of the electronic charge density.

- Journal of computational chemistry
- Published almost 8 years ago
- Discuss

The topological study of the electronic charge density is useful to obtain information about the kinds of bonds (ionic or covalent) and the atom charges on a molecule or crystal. For this study, it is necessary to calculate, at every space point, the electronic density and its electronic density derivatives values up to second order. In this work, a grid-based method for these calculations is described. The library, implemented for three dimensions, is based on a multidimensional Lagrange interpolation in a regular grid; by differentiating the resulting polynomial, the gradient vector, the Hessian matrix and the Laplacian formulas were obtained for every space point. More complex functions such as the Newton-Raphson method (to find the critical points, where the gradient is null) and the Cash-Karp Runge-Kutta method (used to make the gradient paths) were programmed. As in some crystals, the unit cell has angles different from 90°, the described library includes linear transformations to correct the gradient and Hessian when the grid is distorted (inclined). Functions were also developed to handle grid containing files (grd from DMol® program, CUBE from Gaussian® program and CHGCAR from VASP® program). Each one of these files contains the data for a molecular or crystal electronic property (such as charge density, spin density, electrostatic potential, and others) in a three-dimensional (3D) grid. The library can be adapted to make the topological study in any regular 3D grid by modifying the code of these functions. © 2012 Wiley Periodicals, Inc.

#### 28

##### Complexes of 4-substituted phenolates with HF and HCN: Energy decomposition and electronic structure analyses of hydrogen bonding.

- Journal of computational chemistry
- Published over 7 years ago
- Discuss

We have computationally studied para-X-substituted phenols and phenolates (X = NO, NO(2) , CHO, COMe, COOH, CONH(2) , Cl, F, H, Me, OMe, and OH) and their hydrogen-bonded complexes with B(-) and HB (B = F and CN), respectively, at B3LYP/6-311++G** and BLYP-D/QZ4P levels of theory. Our purpose is to explore the structures and stabilities of these complexes. Moreover, to understand the emerging trends, we have analyzed the bonding mechanisms using the natural bond orbital scheme as well as Kohn-Sham molecular orbital (MO) theory in combination with quantitative energy decomposition analyses [energy decomposition analysis (EDA), extended transition state-natural orbitals for chemical valence (ETS-NOCV)]. These quantitative analyses allow for the construction of a simple physical model that explains all computational observations. © 2012 Wiley Periodicals, Inc.

#### 28

##### Electronic fluxes during diels-alder reactions involving 1,2-benzoquinones: mechanistic insights from the analysis of electron localization function and catastrophe theory.

- Journal of computational chemistry
- Published almost 8 years ago
- Discuss

By means of the joint use of electron localization function (ELF) and Thom’s catastrophe theory, a theoretical analysis of the energy profile for the hetero-Diels-Alder reaction of 4-methoxy-1,2-benzoquinone 1 and methoxyethylene 2 has been carried out. The 12 different structural stability domains obtained by the bonding evolution theory have been identified as well as the bifurcation catastrophes (fold and cusp) responsible for the changes in the topology of the system. This analysis permits finding a relationship between the ELF topology and the evolution of the bond breaking/forming processes and electron pair rearrangements through the reaction progress in terms of the different ways of pairing up the electrons. The reaction mechanism corresponds to an asynchronous electronic flux; first, the O1C5 bond is formed by the nucleophilic attack of the C5 carbon of the electron rich ethylene 2 on the most electrophilically activated carbonyl O1 oxygen of 1, and once the σ bond has been completed, the formation process of the second O4C6 bond takes place. In addition, the values of the local electrophilicity and local nucleophilcity indices in the framework of conceptual density functional theory accounts for the asychronicity of the process as well as for the observed regioselectivity. © 2012 Wiley Periodicals, Inc.

#### 27

##### VinaMPI: Facilitating multiple receptor high-throughput virtual docking on high-performance computers

- Journal of computational chemistry
- Published about 7 years ago
- Discuss

The program VinaMPI has been developed to enable massively large virtual drug screens on leadership-class computing resources, using a large number of cores to decrease the time-to-completion of the screen. VinaMPI is a massively parallel Message Passing Interface (MPI) program based on the multithreaded virtual docking program AutodockVina, and is used to distribute tasks while multithreading is used to speed-up individual docking tasks. VinaMPI uses a distribution scheme in which tasks are evenly distributed to the workers based on the complexity of each task, as defined by the number of rotatable bonds in each chemical compound investigated. VinaMPI efficiently handles multiple proteins in a ligand screen, allowing for high-throughput inverse docking that presents new opportunities for improving the efficiency of the drug discovery pipeline. VinaMPI successfully ran on 84,672 cores with a continual decrease in job completion time with increasing core count. The ratio of the number of tasks in a screening to the number of workers should be at least around 100 in order to have a good load balance and an optimal job completion time. The code is freely available and downloadable. Instructions for downloading and using the code are provided in the Supporting Information. © 2013 Wiley Periodicals, Inc.

#### 27

##### Message passing interface and multithreading hybrid for parallel molecular docking of large databases on petascale high performance computing machines

- Journal of computational chemistry
- Published over 7 years ago
- Discuss

A mixed parallel scheme that combines message passing interface (MPI) and multithreading was implemented in the AutoDock Vina molecular docking program. The resulting program, named VinaLC, was tested on the petascale high performance computing (HPC) machines at Lawrence Livermore National Laboratory. To exploit the typical cluster-type supercomputers, thousands of docking calculations were dispatched by the master process to run simultaneously on thousands of slave processes, where each docking calculation takes one slave process on one node, and within the node each docking calculation runs via multithreading on multiple CPU cores and shared memory. Input and output of the program and the data handling within the program were carefully designed to deal with large databases and ultimately achieve HPC on a large number of CPU cores. Parallel performance analysis of the VinaLC program shows that the code scales up to more than 15K CPUs with a very low overhead cost of 3.94%. One million flexible compound docking calculations took only 1.4 h to finish on about 15K CPUs. The docking accuracy of VinaLC has been validated against the DUD data set by the re-docking of X-ray ligands and an enrichment study, 64.4% of the top scoring poses have RMSD values under 2.0 Å. The program has been demonstrated to have good enrichment performance on 70% of the targets in the DUD data set. An analysis of the enrichment factors calculated at various percentages of the screening database indicates VinaLC has very good early recovery of actives. © 2013 Wiley Periodicals, Inc.