Multiwfn is a multifunctional program for wavefunction analysis. Its main functions are: (1) Calculating and visualizing real space function, such as electrostatic potential and electron localization function at point, in a line, in a plane or in a spatial scope. (2) Population analysis. (3) Bond order analysis. (4) Orbital composition analysis. (5) Plot density-of-states and spectrum. (6) Topology analysis for electron density. Some other useful utilities involved in quantum chemistry studies are also provided. The built-in graph module enables the results of wavefunction analysis to be plotted directly or exported to high-quality graphic file. The program interface is very user-friendly and suitable for both research and teaching purpose. The code of Multiwfn is substantially optimized and parallelized. Its efficiency is demonstrated to be significantly higher than related programs with the same functions. Five practical examples involving a wide variety of systems and analysis methods are given to illustrate the usefulness of Multiwfn. The program is free of charge and open-source. Its precompiled file and source codes are available from http://multiwfn.codeplex.com.

Multiwfn: a multifunctional wavefunction analyzer.

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.

Noncovalent interactions hold the key to understanding many chemical, biological, and technological problems. Describing these noncovalent interactions accurately, including their positions in real space, constitutes a first step in the process of decoupling the complex balance of forces that define noncovalent interactions. Because of the size of macromolecules, the most common approach has been to assign van der Waals interactions (vdW), steric clashes (SC), and hydrogen bonds (HBs) based on pairwise distances between atoms according to their vdW radii. We recently developed an alternative perspective, derived from the electronic density: the non-covalent interactions (NCI) index [J. Am. Chem. Soc. 2010, 132, 6498]. This index has the dual advantages of being generally transferable to diverse chemical applications and being very fast to compute, since it can be calculated from promolecular densities. Thus, NCI analysis is applicable to large systems, including proteins and DNA, where analysis of noncova...

/pdf/nciplot-a-program-for-plotting-noncovalent-interaction-5bhkffn4zn.pdf

NCIPLOT: A Program for Plotting Noncovalent Interaction Regions

Four Cu(I) complexes were synthesized with a family of pyridylmethylamide ligands, HLR [HLR = N-(2-pyridylmethyl)acetamide, R = null; 2,2-dimethyl-N-(2-pyridylmethyl)propionamide, R = Me3; 2,2,2-triphenyl-N-(2-pyridylmethyl)acetamide, R = Ph3)]. Complexes 1–3 were synthesized from the respective ligand and [Cu(CH3CN)4]PF6 in a 2 : 1 molar ratio: [Cu(HL)2]PF6 (1), [Cu2(HLMe3)4](PF6)2 (2), [Cu(HLPh3)2]PF6 (3). Complex 4, [Cu(HL)(CH3CN)(PPh3)]PF6, was synthesized from the reaction of HL with [Cu(CH3CN)4]PF6 and PPh3 in a 1 : 1 : 1 molar ratio. X-Ray crystal structures reveal that complexes 1, 3 and 4 are mononuclear Cu(I) species, while complex 2 is a Cu(I) dimer. The copper ions are four-coordinate with geometries ranging from distorted tetrahedral to seesaw in 1, 2, and 4. Complexes 1 and 2 are very air sensitive and they display similar electrochemical properties. The coordination geometry of complex 3 is nearly linear, two-coordinate. Complex 3 is exceptionally stable with respect to oxidation in the air, and its cyclic voltammetry shows no oxidation wave in the range of 0–1.5 V. The unusual inertness of complex 3 towards oxidation is attributed to the protection from bulky triphenyl substituent of the HLPh3 ligand. A new geometric parameter for four-coordinate compounds, τ4, is proposed as an improved, simple metric for quantitatively evaluating the geometry of four-coordinate complexes and compounds.

Structural variation in copper(I) complexes with pyridylmethylamide ligands: structural analysis with a new four-coordinate geometry index, tau4.

Jaguar is an ab initio quantum chemical program that specializes in fast electronic structure predictions for molecular systems of medium and large size. Jaguar focuses on computational methods with reasonable computational scaling with the size of the system, such as density functional theory (DFT) and local second-order Moller–Plesset perturbation theory. The favorable scaling of the methods and the high efficiency of the program make it possible to conduct routine computations involving several thousand molecular orbitals. This performance is achieved through a utilization of the pseudospectral approximation and several levels of parallelization. The speed advantages are beneficial for applying Jaguar in biomolecular computational modeling. Additionally, owing to its superior wave function guess for transition-metal-containing systems, Jaguar finds applications in inorganic and bioinorganic chemistry. The emphasis on larger systems and transition metal elements paves the way toward developing Jaguar for its use in materials science modeling. The article describes the historical and new features of Jaguar, such as improved parallelization of many modules, innovations in ab initio pKa prediction, and new semiempirical corrections for nondynamic correlation errors in DFT. Jaguar applications in drug discovery, materials science, force field parameterization, and other areas of computational research are reviewed. Timing benchmarks and other results obtained from the most recent Jaguar code are provided. The article concludes with a discussion of challenges and directions for future development of the program. © 2013 Wiley Periodicals, Inc.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/qua.24481

Jaguar: A high-performance quantum chemistry software program with strengths in life and materials sciences

Electron transfer (ET) reactions provide a nexus among chemistry, biochemistry, and physics. These reactions underpin the “power plants” and “power grids” of bioenergetics, and they challenge us to understand how evolution manipulates structure to control ET kinetics. Ball-and-stick models for the machinery of electron transfer, however, fail to capture the rich electronic and nuclear dynamics of ET molecules: these static representations disguise, for example, the range of thermally accessible molecular conformations. The influence of structural fluctuations on electron-transfer kinetics is amplified by the exponential decay of electron tunneling probabilities with distance, as well as the delicate interference among coupling pathways. Fluctuations in the surrounding medium can also switch transport between coherent and incoherent ET mechanisms − and may gate ET so that its kinetics is limited by conformational interconversion times, rather than by the intrinsic ET time scale. Moreover, preparation of a ...

Steering Electrons on Moving Pathways

Formation of a transient [myoglobin (Mb), cytochrome b5 (cyt b5)] complex is required for the reductive repair of inactive ferri-Mb to its functional ferro-Mb state. The [Mb, cyt b5] complex exhibits dynamic docking (DD), with its cyt b5 partner in rapid exchange at multiple sites on the Mb surface. A triple mutant (Mb(3M)) was designed as part of efforts to shift the electron-transfer process to the simple docking (SD) regime, in which reactive binding occurs at a restricted, reactive region on the Mb surface that dominates the docked ensemble. An electrostatically-guided Brownian dynamics (BD) docking protocol was used to generate an initial ensemble of reactive configurations of the complex between unrelaxed partners. This ensemble samples a broad and diverse array of heme–heme distances and orientations. These configurations seeded all-atom constrained molecular dynamics simulations (MD) to generate relaxed complexes for the calculation of electron tunneling matrix elements (TDA) through tunneling-pathway analysis. This procedure for generating an ensemble of relaxed complexes combines the ability of BD calculations to sample the large variety of available conformations and interprotein distances, with the ability of MD to generate the atomic level information, especially regarding the structure of water molecules at the protein–protein interface, that defines electron-tunneling pathways. We used the calculated TDA values to compute ET rates for the [Mb(wt), cyt b5] complex and for the complex with a mutant that has a binding free energy strengthened by three D/E → K charge-reversal mutations, [Mb(3M), cyt b5]. The calculated rate constants are in agreement with the measured values, and the mutant complex ensemble has many more geometries with higher TDA values than does the wild-type Mb complex. Interestingly, water plays a double role in this electron-transfer system, lowering the tunneling barrier as well as inducing protein interface remodeling that screens the repulsion between the negatively-charged propionates of the two hemes.

/pdf/interfacial-hydration-dynamics-and-electron-transfer-multi-5kxge1dd1q.pdf

Interfacial hydration, dynamics and electron transfer: multi-scale ET modeling of the transient [myoglobin, cytochrome b5] complex.

In inflammatory bowel disease (IBD) patients, the enzyme phospholipase A2 (PLA2) is overexpressed in the inflamed intestinal tissue, and hence may be exploited as a prodrug-activating enzyme allowing drug targeting to the site(s) of gut inflammation. The purpose of this work was to develop powerful modern computational approaches, to allow optimized a-priori design of phospholipid (PL) based prodrugs for IBD drug targeting. We performed simulations that predict the activation of PL-drug conjugates by PLA2 with both human and bee venom PLA2. The calculated results correlated well with in-vitro experimental data. In conclusion, a-priori drug design using a computational approach complements and extends experimentally derived data, and may improve resource utilization and speed drug development.

Phospholipid-Based Prodrugs for Drug Targeting in Inflammatory Bowel Disease: Computational Optimization and In-Vitro Correlation.

The monosila analogue of the C6Me7+ cation can be accessed by hydride abstraction from (pentamethylcyclopentadienyl)dimethylsilane, Cp*SiMe2H. Treatment of this material with triphenylmethylium tetrakis(pentafluorophenyl)borate in dichloromethane at −50 °C produced a single cationic species stable for a day or longer. The presence of a single resonance for the ring methyls and a single resonance for the silyl methyls in both the 1H and the 13C spectra indicated either a static structure with fivefold symmetry or a dynamic structure in which the dimethylsilyl group equilibrates rapidly among positions. Density function theory calculations found a minimum whose structure had the unsymmetrical silabicyclo[3.1.0]hexenyl structure. Calculated 1H, 13C, and 29Si chemical shifts for this structure were in good agreement with the observed values for this structure. Another calculated minimum was the silacyclohexadienyl structure, but its calculated chemical shifts were very different from those observed. We conclude that the cation has the bicyclic structure but rapidly interconverts the ring positions to produce the very simple NMR spectra.

The C5SiMe7+ cation: pyramidal, bicyclic, or cyclohexadienyl?

Targeting drugs to the inflamed intestinal tissue(s) represents a major advancement in the treatment of inflammatory bowel disease (IBD). In this work we present a powerful in-silico modeling approach to guide the molecular design of novel prodrugs targeting the enzyme PLA2, which is overexpressed in the inflamed tissues of IBD patients. The prodrug consists of the drug moiety bound to the sn-2 position of phospholipid (PL) through a carbonic linker, aiming to allow PLA2 to release the free drug. The linker length dictates the affinity of the PL-drug conjugate to PLA2, and the optimal linker will enable maximal PLA2-mediated activation. Thermodynamic integration and Weighted Histogram Analysis Method (WHAM)/Umbrella Sampling method were used to compute the changes in PLA2 transition state binding free energy of the prodrug molecule (∆∆Gtr) associated with decreasing/increasing linker length. The simulations revealed that 6-carbons linker is the optimal one, whereas shorter or longer linkers resulted in decreased PLA2-mediated activation. These in-silico results were shown to be in excellent correlation with experimental in-vitro data. Overall, this modern computational approach enables optimization of the molecular design of novel prodrugs, which may allow targeting the free drug specifically to the diseased intestinal tissue of IBD patients.

Shahar Keinan

Papers

Steering Electrons on Moving Pathways

Interfacial hydration, dynamics and electron transfer: multi-scale ET modeling of the transient [myoglobin, cytochrome b5] complex.

Phospholipid-Based Prodrugs for Drug Targeting in Inflammatory Bowel Disease: Computational Optimization and In-Vitro Correlation.

The C5SiMe7+ cation: pyramidal, bicyclic, or cyclohexadienyl?

Computational modeling and in-vitro/in-silico correlation of phospholipid-based prodrugs for targeted drug delivery in inflammatory bowel disease