Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

Fantasy versus reality in fragment-based quantum chemistry

Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.

Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective

Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the φ and ψ torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the φ/ψ scans in Ramachandran plots. At the molecular level, it is surprisingly electrostatic destabilization that causes the high-energy regions in the Ramachandran plot, not molecular steric hindrance (related to the intra-atomic energy). At the functional group and atomic levels, the importance of key peptide atoms (Oi-1 , Ci , Ni , Ni+1 ) and some sidechain hydrogen atoms (Hγ ) are identified as responsible for the destabilization seen in the energetically disfavored Ramachandran regions. Consistently, the Oi-1 atoms are particularly important for the explanation of dipeptide intrinsic behavior, where electrostatic and steric destabilization unusually complement one another. The findings suggest that, at least for these dipeptides, it is the peptide group atoms that dominate the intrinsic behavior, more so than the sidechain atoms. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Comparative Assessment of Different RNA Tetranucleotides from the DFT-D3 and Force Field Perspective

Classical force field (FF) molecular dynamics (MD) simulations of RNA tetranucleotides have substantial problems in reproducing conformer populations indicated by NMR experiments. To provide more information about the possible sources of errors, we performed quantum mechanical (QM, TPSS-D3/def2-TZVP) and molecular mechanics (MM, AMBER parm99bsc0+χ_OL3) calculations of different r(CCCC), r(GACC), and r(UUUU) conformers obtained from explicit solvent MD simulations. Solvent effects in the static QM and MM calculations were mimicked using implicit solvent models (COSMO and Poisson-Boltzmann, respectively). The comparison of QM and MM geometries and energies revealed that the two methodologies provide qualitatively consistent results in most of the cases. Even though we found some differences, these were insufficient to indicate any systematic corrections of the RNA FF terms that could improve the performance of classical MD in simulating tetranucleotides. On the basis of these findings, we inferred that the overpopulation of intercalated conformers in the MD simulations of RNA tetramers, which were not observed experimentally, might be predominantly caused by imbalanced water-solvent and water-water interactions. Apart from the large-scale QM calculations performed to assess the performance of the AMBER FF, a representative spectrum of faster QM methods was tested.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

Gas-phase basicities of polyfunctional molecules. Part 4: Carbonyl groups as basic sites

This article constitutes the fourth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435). This fourth part is devoted to carbonyl containing polyfunctional molecules. After a short reminder of the methods of determination of gas-phase basicity and the underlying physicochemical concepts, specific examples are examined under two major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) ketones, diketones, ketoalcohols, and ketoethers are considered. A second chapter describes the protonation energetic of gaseous acids and derivatives including diacids, diesters, diamides, anhydrides, imides, ureas, carbamates, amino acid derivatives, and peptides. Experimental data were re-evaluated according to the presently adopted basicity scale. Structural and energetic information given by G3 and G4 quantum chemistry computations on typical systems are presented.

Towards an understanding of the nucleation of alpha-para amino benzoic acid from ethanolic solutions: a multi-scale approach

The molecular assembly and subsequent nucleation of para-amino benzoic acid (PABA) from ethanolic solutions is probed using a multi-scale and multi-technique approach. This is applied by examining and interrelating information regarding the molecular, solution-state, cluster, solid-state and surface structures to understand why the alpha form of PABA is crystallised in preference to its low temperature beta form. Calculations suggest that conformational changes within the solute molecule play little or no role in directing the nucleation of either the alpha or beta crystal forms. Combined ab initio and molecular dynamics calculations of the stability of small clusters in solution suggests that the hydrogen-bonded carboxylic acid dimers, present in the alpha structure, are the most stable in solution and play a major role in the self-assembly and polymorphic expression of the alpha form in ethanol in preference to the beta form. These calculations are in good agreement with X-ray small-angle scattering analysis which reveals the presence of PABA clusters in ethanol which are consistent with the size and shape of a carboxylic acid dimer. SAXS studies also reveal the presence of larger cluster structures in a size range 10-40 nm which appear to grow, perhaps reflecting a change in the balance between monomers and dimers within the solution during the nucleation process. The results of crystallisation-kinetics experiments indicate an instantaneous nucleation mechanism where the number of instantaneously nucleated crystallites is calculated to be 1360-660 nuclei per ml and the subsequent growth is found to be only rate limited by diffusion of the growth unit to the crystallite surface. A linear dependence of growth rate with respect to supersaturation is observed for the (0 1 -1) capping face, which is associated with strong π-π stacking interactions. This is consistent with a solid-on-solid mechanism associated with surface roughened growth and concomitant poor lattice-perfection. Conversely, the side (1 0 -1) surface has a growth mechanism consistent with a 2D nucleation birth and spread mechanism. Hence, these mechanisms result in very fast growth along the b-axis and the needle-like morphology that is observed for alpha-PABA.

Comprehensive analysis of energy minima of the 20 natural amino acids

Energy minima of the 20 natural amino acids (capped by a peptide bond at both the N and the C termini, CH3-C(═O)-N(H)-(H)Cα(R)-C(═O)-N(H)-CH3), were obtained by ab initio geometry optimization. Starting with a large number of minima, quickly generated by MarvinView, geometry optimization at the HF/6-31G(d,p) level of theory reduced the number of minima, followed by further optimization at the B3LYP/apc-1 and MP2/cc-pVDZ levels, which caused some minima to disappear and some stable minima to migrate on the Ramachandran map. There is a relation between the number of minima and the size and the flexibility of the side chain. The energy minima of the 20 amino acids are mainly located in the regions of βL, γL, δL, and αL of the Ramachandran map. Multipole moments of atoms occurring in the fragment [-NH-Cα-C(═O)-] common to all 20 amino acids were calculated at the three levels of theory mentioned above. The near parallelism in behavior of these moments between levels of theory is beneficial toward estimating moments with the more expensive B3LYP and MP2 methods from data calculated with the cheaper HF method. Finally, we explored the transferability of properties between different amino acids: the bond length and angles of the common fragment [-NH-Cα(HαCβ)-C'(═O)-] in all amino acids except Gly and Pro. All bond lengths are highly transferable between different amino acids, and the standard deviations are small.

Gas-phase enthalpies of formation and enthalpies of sublimation of amino acids based on isodesmic reaction calculations

Accurate gas-phase enthalpies of formation (ΔfH298°) of 20 common α-amino acids, seven uncommon amino acids, and three small peptides were calculated by combining G4 theory calculations with an isodesmic reaction approach. The internal consistency over a set of ΔfH298°(g) values was achieved by sequential adjustment of their values through the isodesmic reactions. Four amino acids, alanine, β-alanine, sarcosine, and glycine, with reliable internally self-consistent experimental data, were chosen as the key reference compounds. These amino acids together with about 100 compounds with reliable experimental data (their accuracy was supported by G4 calculations) were used to estimate the enthalpies of formation of remaining amino acids. All of the amino acids with the previously established enthalpies of formation were later used as the reference species in the isodesmic reactions for the other amino acids. A systematic comparison was made of 14 experimentally determined enthalpies of formation with the results of calculations. The experimental enthalpies of formation for 10 amino acids were reproduced with good accuracy, but the experimental and calculated values for 4 compounds differed by 11–21 kJ/mol. For these species, the theoretical ΔfH298°(g) values were suggested as more reliable than the experimental values. On the basis of theoretical results, the recommended values for the gas-phase enthalpies of formation were also provided for amino acids for which the experimental ΔfH298°(g) were not available. The enthalpies of sublimation were evaluated for all compounds by taking into account the literature data on the solid-phase enthalpies of formation and the ΔfH298°(g) values recommended in our work. A special attention was paid to the accurate prediction of enthalpies of formation of amino acids from the atomization reactions. The problems associated with conformational flexibility of these compounds and harmonic treatment of low frequency torsional modes were discussed. The surprisingly good agreement between the ΔfH298°(g) values calculated from the atomization and isodesmic reactions is largely the result of a fortuitous mutual compensation of various corrections used in the atomization reaction procedure.

Achieving the CCSD(T) Basis-Set Limit in Sizable Molecular Clusters: Counterpoise Corrections for the Many-Body Expansion

An efficient procedure is introduced to obtain the basis-set limit in electronic structure calculations of large molecular and ionic clusters. This approach is based on a Boys-Bernardi-style counterpoise correction for clusters containing arbitrarily many monomer units, which is rendered computationally feasible by means of a truncated many-body expansion. This affords a tractable way to apply the sequence of correlation-consistent basis sets (aug-cc-pVXZ) to large systems and thereby obtain energies extrapolated to the complete basis set (CBS) limit. A three-body expansion with three-body counterpoise corrections is shown to afford errors of ≲0.1-0.2 kcal/mol with respect to traditional MP2/CBS results, even for challenging systems such as fluoride-water clusters. A triples correction, δCCSD(T) = ECCSD(T) - EMP2, can be estimated accurately and efficiently as well. Because the procedure is embarrassingly parallelizable and requires no electronic structure calculations in systems larger than trimers, it is extendible to very large clusters. As compared to traditional CBS extrapolations, computational time is dramatically reduced even without parallelization.

Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: the challenge of the conformational equilibrium in glycine

The structures, relative stabilities, and infrared spectra of the six low-energy conformers of glycine have been characterized using a state-of-the-art quantum-mechanical approach allowing the bond distances, conformational enthalpies and vibrational frequencies to be determined well within the chemical accuracy. Transition state structures governing interconversion among the different energy minima have also been characterized. In detail, the gas-phase thermodynamic properties (at 15 K and 410 K) of the glycine conformers considered have been obtained with a 1 kJ mol(-1) accuracy, and it has been shown that the employment of DFT geometries usually reduces such accuracy by at most 0.1 kJ mol(-1). Regarding molecular structures, the use of two different composite schemes allowed us to further confirm the suitability of a rather cost-effective approach and provide geometrical parameters with an overall accuracy better than 0.002 Å for distances and 1 degree for angles. Thanks to a hybrid CC/DFT approach, the infrared spectra of all conformers considered and of several deuterated isotopologues have been reproduced (when experimental data were available) or predicted with an accuracy of 10 cm(-1). Finally, the joint thermodynamic and spectroscopic investigation allowed us to shed some light on the possible observation of elusive conformers. On the whole, the high accuracy of the computational results allows us to draw a fully consistent interpretation of the available experimental data and to obtain a more complete characterization of the potential energy surface of glycine.

Peptide bond distortions from planarity: new insights from quantum mechanical calculations and peptide/protein crystal structures

By combining quantum-mechanical analysis and statistical survey of peptide/protein structure databases we here report a thorough investigation of the conformational dependence of the geometry of peptide bond, the basic element of protein structures. Different peptide model systems have been studied by an integrated quantum mechanical approach, employing DFT, MP2 and CCSD(T) calculations, both in aqueous solution and in the gas phase. Also in absence of inter-residue interactions, small distortions from the planarity are more a rule than an exception, and they are mainly determined by the backbone ψ dihedral angle. These indications are fully corroborated by a statistical survey of accurate protein/peptide structures. Orbital analysis shows that orbital interactions between the σ system of C(α) substituents and the π system of the amide bond are crucial for the modulation of peptide bond distortions. Our study thus indicates that, although long-range inter-molecular interactions can obviously affect the peptide planarity, their influence is statistically averaged. Therefore, the variability of peptide bond geometry in proteins is remarkably reproduced by extremely simplified systems since local factors are the main driving force of these observed trends. The implications of the present findings for protein structure determination, validation and prediction are also discussed.

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.