A second look at the role of the fast Fourier transform as an elliptic solver

A fast cosine transform (FCT) is coupled with a tridiagonal solver for the purpose of solving the Poisson equation on irregular and non‐uniform rectangular staggered grids. This kind of solution is required for the pressure field during the simulation of the incompressible Navier–Stokes equations when using the projection method. A new technique using the FCT–tridiagonal solver is derived for the cases where the boundaries of the flow regime do not coincide with the boundaries of the computational domain and for non‐uniform grids. The technique is based on an iterative procedure where a defect equation is solved in every iteration, followed by a relaxation procedure. The method is investigated analytically and numerically to show that the solution converges as a geometric series. The method is further investigated for the effects of the relative size of the rigid body, the grid stretching, size and aspect ratio. The new solver is incorporated with the direct numerical simulation (DNS) and large eddy simulation (LES) techniques to simulate the flows around a backward‐facing step and a 3D rectangular obstacle, yielding results that qualitatively compare well with known results. Copyright © 2005 John Wiley & Sons, Ltd.     

On Trees and Characters

A new family of trees, defined in term of Young diagrams, is introduced. Values of central characters of the symmetric group are represented as a weighted enumeration of such trees. The proof involves a new decomposition theorem for representations corresponding to general shapes.     

Can DFT methods correctly and efficiently predict the coordination number of copper(I) complexes? A case study

The coordination number of various experimentally known Cu(I) compounds is studied using density functional theory. Various basis sets are tested, aiming to establish a reliable level for prediction of the coordination number of these and other Cu(I) complexes. It is found that most levels exhibit correct trends, namely, the bulkier ligands demonstrate larger preference for coordination of two ligands. Proper absolute values are obtained when dispersion corrections are also included in the calculations. It is concluded that the fairly small modified 6‐31+G* basis set due to Pulay represents a good compromise between accuracy and efficiency, followed by Balabanov and Peterson's all‐electron aug‐cc‐pVDZ basis set. The overall energy is decomposed into various components whose relative contribution to the overall tendency of forming a complex with a particular coordination is examined. It is shown that two opposing contributions play a major role: the interaction energy of the ligand being added and the deformation energy of the copper's coordination sphere prior to the ligand addition. The former being a stabilizing contribution, leads to higher coordination numbers while the later, a destabilizing contribution, is shown to favor lower coordination numbers. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Apparent NAC effect in chorismate mutase reflects electrostatic transition state stabilization

The catalytic reaction of chorismate mutase (CM) has been the subject of major current attention. Nevertheless, the origin of the catalytic power of CM remains an open question. In particular, it has not been clear whether the enzyme works by providing electrostatic transition state stabilization (TSS), by applying steric strain, or by populating near attack conformation (NAC). The present work explores this issue by a systematic quantitative analysis. The overall catalytic effect is reproduced by the empirical valence bond (EVB) method. In addition, the binding free energy of the ground state and the transition state is evaluated, demonstrating that the enzyme works by TSS. Furthermore, the evaluation of the electrostatic contribution to the reduction of the activation energy establishes that the TSS results from electrostatic effects. It is also found that the apparent NAC effect is not the reason for the catalytic effect but the result of the TSS. It is concluded that in CM as in other enzymes the key catalytic effect is electrostatic TSS. However, since the charge distribution of the transition state and the reactant state is similar, the stabilization of the transition state leads to reduction in the distance between the reacting atoms in the reactant state.     

Density embedded VB/MM: a hybrid ab initio VB/MM with electrostatic embedding

A hybrid QM/MM method that combines ab initio valence-bond (VB) with molecular mechanics (MM) is presented. The method utilizes the ab initio VB approach to describe the reactive fragments and MM to describe the environment thus allows VB calculations of reactions in large biological systems. The method, termed density embedded VB/MM (DE-VB/MM), is an extension of the recently developed VB/MM method. It involves calculation of the electrostatic interaction between the reactive fragments and their environment using the electrostatic embedding scheme. Namely, the electrostatic interactions are represented as one-electron integrals in the ab initio VB Hamiltonian, hence taking into account the wave function polarization of the reactive fragments due to the environment. Moreover, the assumptions that were utilized in an earlier version of the method, VB/MM, to formulate the electrostatic interactions effect on the off-diagonal matrix elements are no longer required in the DE-VB/MM methodology. Using DE-VB/MM, one can calculate, in addition to the adiabatic ground state reaction profile, the energy of the diabatic VB configurations as well as the VB state correlation diagram for the reaction. The abilities of the method are exemplified on the identity SN2 reaction of a chloride anion with methyl chloride in aqueous solution. Both the VB configurations diagram and the state correlation diagram are presented. The results are shown to be in very good agreement with both experimental and other computational data, suggesting that DE-VB/MM is a proper method for application to different reactivity problems in biological systems.     

Valence Bond and Enzyme Catalysis: A Time To Break Down and a Time To Build Up

Understanding enzyme catalysis and developing ability to control of it are two great challenges in biochemistry. A few successful examples of computational‐based enzyme design have proved the fantastic potential of computational approaches in this field, however, relatively modest rate enhancements have been reported and the further development of complementary methods is still required. Herein we propose a conceptually simple scheme to identify the specific role that each residue plays in catalysis. The scheme is based on a breakdown of the total catalytic effect into contributions of individual protein residues, which are further decomposed into chemically interpretable components by using valence bond theory. The scheme is shown to shed light on the origin of catalysis in wild‐type haloalkane dehalogenase (wt‐DhlA) and its mutants. Furthermore, the understanding gained through our scheme is shown to have great potential in facilitating the selection of non‐optimal sites for catalysis and suggesting effective mutations to enhance the enzymatic rate.