A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

Molecular structure does not easily identify the intricate noncovalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect noncovalent interactions in real space, based on the electron density and its derivatives. Our approach reveals the underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of nonbonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.

/pdf/revealing-noncovalent-interactions-1oqknfvmfi.pdf

Revealing noncovalent interactions.

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

QM/MM Methods for Biomolecular Systems

Dispersion corrections to standard Kohn–Sham density functional theory (DFT) are reviewed. The focus is on computationally efficient methods for large systems that do not depend on virtual orbitals or rely on separated fragments. The recommended approaches (van der Waals density functional and DFT-D) are asymptotically correct and can be used in combination with standard or slightly modified (short-range) exchange–correlation functionals. The importance of the dispersion energy in intramolecular cases (conformational problems and thermochemistry) is highlighted. © 2011 John Wiley & Sons, Ltd. WIREs Comput Mol Sci 2011 1 211-228 DOI: 10.1002/wcms.30

Density functional theory with London dispersion corrections

Introduction
PaDEL-Descriptor is a software for calculating molecular descriptors and fingerprints. The software currently calculates 797 descriptors (663 1D, 2D descriptors, and 134 3D descriptors) and 10 types of fingerprints. These descriptors and fingerprints are calculated mainly using The Chemistry Development Kit. Some additional descriptors and fingerprints were added, which include atom type electrotopological state descriptors, McGowan volume, molecular linear free energy relation descriptors, ring counts, count of chemical substructures identified by Laggner, and binary fingerprints and count of chemical substructures identified by Klekota and Roth.

Methods
PaDEL-Descriptor was developed using the Java language and consists of a library component and an interface component. The library component allows it to be easily integrated into quantitative structure activity relationship software to provide the descriptor calculation feature while the interface component allows it to be used as a standalone software. The software uses a Master/Worker pattern to take advantage of the multiple CPU cores that are present in most modern computers to speed up calculations of molecular descriptors.

Results
The software has several advantages over existing standalone molecular descriptor calculation software. It is free and open source, has both graphical user interface and command line interfaces, can work on all major platforms (Windows, Linux, MacOS), supports more than 90 different molecular file formats, and is multithreaded.

Conclusion
PaDEL-Descriptor is a useful addition to the currently available molecular descriptor calculation software. The software can be downloaded at http://padel.nus.edu.sg/software/padeldescriptor. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011

https://onlinelibrary.wiley.com/doi/pdf/10.1002/jcc.21707

PaDEL‐descriptor: An open source software to calculate molecular descriptors and fingerprints

The pharmaceutical agent theophylline (THEO) is primarily used as a bronchodilator and is commercially available in both tablet and liquid dosage forms. THEO is highly hygroscopic, reducing its stability, overall shelf-life, and therefore usage as a drug. THEO and dicarboxylic acid cocrystals were designed by Trask et al. in an attempt to decrease the hygroscopic behaviour of THEO; cocrystallisation of THEO with malonic acid (MA) did not improve the hygroscopic stability of THEO in simulated atmospheric humidity testing. The current study employed high-resolution X-ray crystallography, and Density Functional Theory (DFT) calculations to examine the electron density distribution (EDD) changes between the cocrystal and its individual components. The EED changes identified the reasons why the THEO:MA cocrystal did not alter the hygroscopic profile of THEO. The cocrystal was equally porous, with atomic packing factors (APF) similar to those of THEO 0.73 vs. 0.71, respectively. The THEO:MA (1) cocrystal structure is held together by an array of interactions; a heterogeneous synthon between the imidazole and a carboxylic fragment stabilising the asymmetric unit, a pyrimidine-imidazole homosynthon, and an aromatic cycle stack between two THEO moieties have been identified, providing 9.7–12.9 kJ mol−1 of stability. These factors did not change the overall relative stability of the cocrystal relative to its individual THEO and MA components, as shown by cocrystal (1) and THEO being equally stable, with calculated lattice energies within 2.5 kJ mol−1 of one other. The hydrogen bond analysis and fragmented atomic charge analysis highlighted that the formation of (1) combined both the EDD of THEO and MA with no net chemical change, suggesting that the reverse reaction — (1) back to THEO and MA — is of equal potential, ultimately producing THEO hydrate formation, in agreement with the work of Trask et al. These results highlight that a review of the EDD change associated with a chemical reaction can aid in understanding cocrystal design. In addition, they indicate that cocrystal design requires further investigation before becoming a reliable process, with particular emphasis on identifying the appropriate balance of synthon engineering, weak interactions, and packing dynamics.

/pdf/an-experimental-and-theoretical-charge-density-study-of-1qf1f93t.pdf

An experimental and theoretical charge density study of theophylline and malonic acid cocrystallization

A novel quinoline-substituted pyrimidine ligand forms two different coloured complexes upon reaction with Re(CO)5Br. These compounds display distinct photophysical properties that are dictated by their stereochemistry.

/pdf/stereoelectronic-control-of-photophysics-red-and-yellow-1cmiznnbtw.pdf

Stereoelectronic control of photophysics: red and yellow axial and equatorial anomers of a rhenium-quinoline complex

Thiation of N , N -biphthalimide in a search for new axially chiral ligands gave a mixture of oligothiobiphthalimides; X-ray crystallography shows that the achiral isomer of dithiodioxobiphthalimide, 1,3-dithioxo-1H-3H-[2,2 ' ]biisoindolyl-1 ' ,3 ' -dione is formed, which may imply a change in the mechanism of thiation between mono- and di-phthalimides.

Synthesis and X-ray structure determination of thionated biphthalimides

Hartree-Fock, density functional, and MP2 methods are applied to the problem of designing neutral, bicyclic C3-symmetric cages incorporating interacting anion- and cation-binding sites which strongly bind NaCl as an ion contact pair. A large number of trial ligands L and their complexes L:NaCl are tested, with the focus on maximizing binding by (i) optimizing the cavity size and shape and (ii) varying the nature of the anion- and cation-binding functionalities. The corresponding complexes L:Cl(-) and L:Na(+) are also studied in some detail. An analysis of their structures and charge distributions helps to build a consistent picture of the requirements for a successful NaCl binding. The 'best' candidate ligand utilizes a tripodal triether-substituted amine N(CH2CH2OR-)3 to bind the sodium cation; three thiourea groups in a tripodal arrangement with a 1,3,5-trisubstituted benzyl spacer group {C6H3(CH2NHC [Formula: see text] XNH-)3 X=O,S} to bind chloride; and a -CH2CH2- spacer linking the two binding sites. A simple Quantitative Structure-Property analysis suggests that the binding cavity shape and size is near to the optimal one for this system.

James Alexis Platts

Papers

An experimental and theoretical charge density study of theophylline and malonic acid cocrystallization

Stereoelectronic control of photophysics: red and yellow axial and equatorial anomers of a rhenium-quinoline complex

Synthesis and X-ray structure determination of thionated biphthalimides

Quantum-Chemical Design of Cryptand-like Ditopic Salt Binders.

Comparison of proposed putative bioactive conformations of altered peptide ligands of myelin basic protein by spectroscopic and modelling studies