Strain in nonclassical silicon hydrides: An insight into the “ultrastability” of sila‐bi[6]prismane (Si18H12) cluster with the endohedrally trapped silicon atom,Si19H12

The recently postulated concept of “ultrastability” and “electron‐deficient aromaticity” (Vach, Nano Lett 2011, 11, 5477; Vach, J Chem Theory Comput 2012, 8, 2088) in a sila‐bi[6]prismane having an additional entrapped silicon atom, Si₁₉H₁₂, has been disproved on the basis of a careful analysis of the energetic characteristics related to the formation of this and other silicon hydrides. The central silicon atom in Si₁₉H₁₂ is weaker bound to other silicon atoms than in conventional tetrahedral silanes; moreover, Si₁₉H₁₂ possesses a significant amount of strain. The role of strain in the formation of the title compounds has been further rationalized by calculating the relative energies for the transformation to a half‐planar conformation in methane and in silane and by calculating the respective strain energies. The strain energy value in Si₁₈H₁₂ is equal to 9.93 eV whereas the same property for Si₁₉H₁₂ lies in range of 6.42–8.85 eV. Two low‐energy isomers of Si₁₉H₁₂ which lie by 2.77 and 3.42 eV (!) lower in energy than the originally considered sila‐bi[6]prismane‐based structure have been proposed. © 2015 Wiley Periodicals, Inc.     

B(13) (+) : a photodriven molecular wankel engine

Revved-up rotary: A molecular Wankel motor, the dual-ring structure B(13) (+) , is driven by circularly-polarized infrared electromagnetic radiation. Calculations show that this illumination leads to a guided unidirectional rotation of the outer ring, which is achieved with rotational frequency of the order of 300?GHz.     

Directional hydrogen bonding in the MM3 force field: II

Extensive calculations on hydrogen bonded systems were carried out using the improved MM3 directional hydrogen bond potential. The resulting total function was reoptimized. Comparisons of the hydrogen bonding potential function from ab initio calculations (MP2/6-31G**); the original MM3(89); and the reoptimized MM3 force field MM3(96), for a variety of C, N, O, and Cl systems including the formamide dimer and formamide–water complex, are described herein. Hydrogen bonding is shown to be a far more complicated and ubiquitous phenomenon than is generally recognized. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1001–1016, 1998     

Towards the design of neutral molecular tweezers for anion recognition

Molecular tweezers are simple molecular receptors that can be characterized by the presence of two flat pincers separated by a more or less rigid tether. They have the ability to form complexes with a substrate molecule by gripping the substrate between the tips of the tweezers in a similar manner to that of mechanical tweezers. Kl?rner et al. synthesized one of the structurally simplest molecular tweezers, which is reported to bind electrodeficient aromatic and aliphatic substrates as well as organic cations. Complexes between these molecular tweezers and electron-rich aromatic, aliphatic, or anionic substrates have not been observed. Inspired by several recent reports that describe the interaction of hexafluorobenzene with electron-rich sites of molecules, we conducted a theoretical study to show the possibility of building molecular tweezers, based on those synthesized by Kl?rner, which were able to bind to anions and thus increase their potential as molecular receptors. We characterized complexes formed between several fluorinated derivatives of simple tweezers and an iodine anion, and analyzed the nature of the intermolecular interactions as well as the energetics for the process of complexation. The stabilization trend reflected by the energetic results when fluorine substituents were added to benzene rings confirms our hypothesis about the possibility of obtaining neutral tweezers composed of aromatic rings that can bind anions.     

Quantum molecular mechanics—a noniterative procedure for the fast Ab Initio calculation of closed shell systems

In this article, we advance the foundations of a strategy to develop a molecular mechanics method based not on classical mechanics and force fields but entirely on quantum mechanics and localized electron‐pair orbitals, which we call quantum molecular mechanics (QMM). Accordingly, we introduce a new manner of calculating Hartree–Fock ab initio wavefunctions of closed shell systems based on variationally preoptimized nonorthogonal electron pair orbitals constructed by linear combinations of basis functions centered on the atoms. QMM is noniterative and requires only one extremely fast inversion of a single sparse matrix to arrive to the one‐particle density matrix, to the electron density, and consequently, to the ab initio electrostatic potential around the molecular system, or cluster of molecules. Although QMM neglects the smaller polarization effects due to intermolecular interactions, it fully takes into consideration polarization effects due to the much stronger intramolecular geometry distortions. For the case of methane, we show that QMM was able to reproduce satisfactorily the energetics and polarization effects of all distortions of the molecule along the nine normal modes of vibration, well beyond the harmonic region. We present the first practical applications of the QMM method by examining, in detail, the cases of clusters of helium atoms, hydrogen molecules, methane molecules, as well as one molecule of HeH⁺ surrounded by several methane molecules. We finally advance and discuss the potentialities of an exact formula to compute the QMM total energy, in which only two center integrals are involved, provided that the fully optimized electron‐pair orbitals are known. © 2012 Wiley Periodicals, Inc.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution.

We have studied the keto-enol interconversion of acetone to understand the mechanism of tautomerism relevant to numerous organic and biochemical processes. Applying the ab initio metadynamics method, we simulated the keto-enol isomerism both in the gas phase and in the presence of water. For the gas-phase intramolecular mechanism we show that no other hydrogen-transfer reactions can compete with the simple keto-enol tautomerism. We obtain an intermolecular mechanism and remarkable participation of water when acetone is solvated by neutral water. The simulations reveal that C deprotonation is the kinetic bottleneck of the keto-enol transformation, in agreement with experimental observations. The most interesting finding is the formation of short H-bonded chains of water molecules that provide the route for proton transfer from the carbon to the oxygen atom of acetone. The mechanistic picture that emerged from the present study involves proton migration and emphasizes the importance of active solvent participation in tautomeric interconversion.     

Use of symmetric rank-one Hessian update in molecular geometry optimization

The use of the symmetric rank-one Hessian update and the Broyden–Fletcher–Goldfarb–Shano (BFGS) update formula are considered in an ab initio molecular geometry optimization algorithm. It is noted that the symmetric rank-one Hessian update has an advantage when compared with the BFGS update formula and this advantage must be more evident in the optimization of molecular geometry, because the total energy surface is a near-quadratic function with a small nonlinearity close to a minimum point. The results obtained in geometry optimization of a test group of molecules support this proposal and show that the use of the symmetric rank-one Hessian update formula permits reduction of the number of energy and gradient evaluations needed to locate a minimum on the energy surface. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1877–1886, 1998     

Ab initio molecular dynamics simulation of a water-hydrogen fluoride equimolar mixture.

Hydrogen fluoride and water can be mixed in any proportion. The resulting solutions have unique acidic properties. In particular, hydrogen fluoride undergoes a weak-to-strong acidity transition with increasing concentration of HF To supplement the knowledge already obtained on dilute or moderately concentrated solutions and gas-phase aggregates, an equimolar mixture is studied here by Car-Parrinello molecular dynamics. The natures of the ions and of the complexes formed in the equimolar liquid were determined. Specifically, H3O+, H5O2+, FH-OH2, and HF2- were spontaneously obtained while only hydronium and fluoride ions pre-exist in the equimolar crystal. The behaviour of the proton in the equimolar liquid was compared with mixtures of other proportions simulated previously in an attempt to relate proton dynamics to acidity. In the same way, the behaviour of HF2- was also examined. In this case, proton localization and transfer appeared to be driven by the fluctuating environment of the solvated ion.     

Berry pseudorotation mechanism for the interpretation of the 19F NMR spectrum in PF5 by ab initio molecular dynamics simulations.

For the first time, theoretical evidence that confirms the importance of the Berry pseudorotation process in the interpretation of the 19F NMR spectrum of phosphorus pentafluoride (PF5) is presented. Ab initio molecular dynamics simulations have been performed to generate a large number of configurations used for NMR parameter computations at the density functional theory level. Two different temperatures were set to highlight the effect of pseudorotation process on the NMR spectrum. Average 19F chemical shifts and spin-spin coupling constants calculated for the five fluorine atoms converge towards the NMR equivalence of the five atoms when the Berry pseudorotation mechanism is accounted for.     

Spectroscopic properties in the liquid phase: combining high-level ab initio calculations and classical molecular dynamics.

We present an integrated computational tool, rooted in density functional theory, the polarizable continuum model, and classical molecular dynamics employing spherical boundary conditions, to study the spectroscopic observables of molecules in solution. As a test case, a modified OPLS-AA force field has been developed and used to compute the UV and NMR spectra of acetone in aqueous solution. The results show that provided the classical force fields are carefully reparameterized and validated, the proposed approach is robust and effective, and can also be used by nonspecialists to provide a general and powerful complement to experimental techniques.     

Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations

Computational studies of proteins based on empirical force fields represent a powerful tool to obtain structure-function relationships at an atomic level, and are central in current efforts to solve the protein folding problem. The results from studies applying these tools are, however, dependent on the quality of the force fields used. In particular, accurate treatment of the peptide backbone is crucial to achieve representative conformational distributions in simulation studies. To improve the treatment of the peptide backbone, quantum mechanical (QM) and molecular mechanical (MM) calculations were undertaken on the alanine, glycine, and proline dipeptides, and the results from these calculations were combined with molecular dynamics (MD) simulations of proteins in crystal and aqueous environments. QM potential energy maps of the alanine and glycine dipeptides at the LMP2/cc-pVxZ//MP2/6-31G* levels, where x = D, T, and Q, were determined, and are compared to available QM studies on these molecules. The LMP2/cc-pVQZ//MP2/6-31G* energy surfaces for all three dipeptides were then used to improve the MM treatment of the dipeptides. These improvements included additional parameter optimization via Monte Carlo simulated annealing and extension of the potential energy function to contain peptide backbone phi, psi dihedral crossterms or a phi, psi grid-based energy correction term. Simultaneously, MD simulations of up to seven proteins in their crystalline environments were used to validate the force field enhancements. Comparison with QM and crystallographic data showed that an additional optimization of the phi, psi dihedral parameters along with the grid-based energy correction were required to yield significant improvements over the CHARMM22 force field. However, systematic deviations in the treatment of phi and psi in the helical and sheet regions were evident. Accordingly, empirical adjustments were made to the grid-based energy correction for alanine and glycine to account for these systematic differences. These adjustments lead to greater deviations from QM data for the two dipeptides but also yielded improved agreement with experimental crystallographic data. These improvements enhance the quality of the CHARMM force field in treating proteins. This extension of the potential energy function is anticipated to facilitate improved treatment of biological macromolecules via MM approaches in general.     

Exploring the hydration of Pb2+: ab initio studies and first-principles molecular dynamics

Even though lead is a well-known toxicant widely scattered throughout the world since antiquity, its chemistry is poorly documented at the molecular level. Here we investigate the hydration of the Pb(2+) ion by means of first-principles molecular dynamics (Car-Parrinello molecular dynamics, CPMD). We found that the hydrated cation is heptacoordinated in a dynamically holodirected arrangement roughly corresponding to a fluxional distorted pentagonal bipyramid. The time-averaged Pb-O bond length is especially large and amounts to 2.70 A with an associated root-mean-square deviation of 0.26 A. This results from a dynamic exchange between short (<2.6 A), intermediate (2.6-3.0 A) and long (>3.0 A) Pb-O bonds. The latter very long Pb-O distance implies that the determination of the coordination number n(c) from experimental work may not necessarily yield values directly comparable to the theoretical value of n(c)=7, since not all experimental techniques would recognize such a long distance as a bond to the metal cation. Pronounced disorders are evidenced in the second shell, characteristic of a chaotropic cation, and exchanges between the first and second shells cannot be excluded on a timescale of a few tens of picoseconds.