226
|
Nachtigallová D, Hobza P, Ritze HH. Electronic splitting in the excited states of DNA base homodimers and -trimers: an evaluation of short-range and Coulombic interactions. Phys Chem Chem Phys 2008; 10:5689-97. [DOI: 10.1039/b806323k] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
|
227
|
Valdes H, Pluháčková K, Pitonák M, Řezáč J, Hobza P. Benchmark database on isolated small peptides containing an aromatic side chain: comparison between wave function and density functional theory methods and empirical force field. Phys Chem Chem Phys 2008; 10:2747-57. [DOI: 10.1039/b719294k] [Citation(s) in RCA: 135] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
228
|
Valdés H, Klusák V, Pitoňák M, Exner O, Starý I, Hobza P, Rulíšek L. Evaluation of the intramolecular basis set superposition error in the calculations of larger molecules: [n]helicenes and Phe-Gly-Phe tripeptide. J Comput Chem 2008; 29:861-70. [DOI: 10.1002/jcc.20841] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
229
|
Černý J, Tong X, Hobza P, Müller-Dethlefs K. Competition between stacking and hydrogen bonding: theoretical study of the phenol⋯Ar cation and neutral complex and comparison to experiment. Phys Chem Chem Phys 2008; 10:2780-4. [DOI: 10.1039/b801460b] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
230
|
Callahan MP, Gengeliczki Z, Svadlenak N, Valdes H, Hobza P, de Vries MS. Non-standard base pairing and stacked structures in methyl xanthine clusters. Phys Chem Chem Phys 2008; 10:2819-26. [DOI: 10.1039/b719874d] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
231
|
Řezáč J, Jurečka P, Riley KE, Černý J, Valdes H, Pluháčková K, Berka K, Řezáč T, Pitoňák M, Vondrášek J, Hobza P. Quantum Chemical Benchmark Energy and Geometry Database for Molecular Clusters and Complex Molecular Systems (www.begdb.com): A Users Manual and Examples. ACTA ACUST UNITED AC 2008. [DOI: 10.1135/cccc20081261] [Citation(s) in RCA: 129] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Our previous benchmark CCSD(T)/ complete basis set limit calculations were collected into a database named begdb - Benchmark Energy and Geometry DataBase. Web-based interface to this database was prepared and is available at www.begdb.com. Users can browse, search and plot the data online or download structures and energy tables.
Collapse
|
232
|
Wang W, Hobza P. Application of Berlin's Theorem to Bond-Length Changes in Isolated Molecules and Red- and Blue-Shifting H-Bonded Clusters. ACTA ACUST UNITED AC 2008. [DOI: 10.1135/cccc20080862] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The origin of the bond-length change in molecule or molecular cluster has been investigated at the MP2/aug-cc-pVDZ level of theory using the electrostatic potential or the electron density difference analysis method. Our results have clearly shown that the bond-length change of a chemical bond is determined mainly by the balance of the electrostatic forces exerted by electrons on the two nuclei. The factors that affect the balance of the electrostatic forces include four parts: (i) The abstraction of the electron density from Berlin's binding region between the two nuclei. (ii) The accumulation of the electron density in Berlin's antibinding regions. (iii) The accumulation of the electron density in Berlin's binding region between the two nuclei. (iv) The abstraction of the electron density from Berlin's antibinding regions. Using the change of the electron density around the two nuclei of a chemical bond, we have succeeded in explaining two important chemical phenomena: (i) breakdown of bond length-bond strength correlation; (ii) the bond-length change in the hydrogen bond.
Collapse
|
233
|
Sedlák R, Jurecka P, Hobza P. Density functional theory-symmetry adapted perturbation treatment energy decomposition of nucleic acid base pairs taken from DNA crystal geometry. J Chem Phys 2007; 127:075104. [PMID: 17718635 DOI: 10.1063/1.2759207] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
234
|
Beran S, Carsky P, Hobza P, Panciř J, Polak R, Slanina Z, Zahradnik R. Theory of Chemical Reactivity. RUSSIAN CHEMICAL REVIEWS 2007. [DOI: 10.1070/rc1978v047n11abeh002291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
235
|
Wang W, Pitonák M, Hobza P. CH Stretching Vibrational Shift of Benzene Dimer: Consistency of Experiment and Calculation. Chemphyschem 2007; 8:2107-11. [PMID: 17763495 DOI: 10.1002/cphc.200700415] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Three low-energy structures of the benzene dimer are investigated by several theoretical procedures (RI-MP2, CCSD(T), RI-DFT-D, DFT/BH&H) covering London dispersion energy. The RI-DFT-D and CCSD(T) calculations are used to verify the DFT/BH&H dimer characteristics, as only at this level can anharmonic calculations be performed. It is ascertained that the T-shaped (C(2v)) structure, in which the C-H stretching frequency of the proton donor shows a significant blue shift, is not stable at any level of theory. It is either a transition structure or a minimum which is easily transformed into a parallel-displaced structure or a T-shaped (C(s)) structure, even at low temperature. Consequently, no blue shift can be detected. On the other hand, the calculated anharmonic IR spectra of the two most stable structures of benzene dimer, namely, the T-shaped (C(s)) and the parallel-displaced ones, give rise to a small red (and no blue) shift of the C-H stretching vibration. This finding is fully consistent with the experimental results.
Collapse
|
236
|
Riley KE, Vondrásek J, Hobza P. Performance of the DFT-D method, paired with the PCM implicit solvation model, for the computation of interaction energies of solvated complexes of biological interest. Phys Chem Chem Phys 2007; 9:5555-60. [PMID: 17957311 DOI: 10.1039/b708089a] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work we investigate the performance of the DFT method, augmented with an empirical dispersion function (DFT-D), paired with the PCM implicit solvation model, for the computation of noncovalent interaction energies of biologically-relevant, solvated model complexes. It is found that this method describes intermolecular interactions within water and ether (protein-like) environments with roughly the same accuracy as in the gas phase. Another important finding is that, when environmental effects are taken into account, the empirical dispersion term associated with the DFT-D method need be modified very little (or not at all), in order to obtain the optimum, most well balanced, performance.
Collapse
|
237
|
Bendová L, Jurecka P, Hobza P, Vondrásek J. Model of Peptide Bond−Aromatic Ring Interaction: Correlated Ab Initio Quantum Chemical Study. J Phys Chem B 2007; 111:9975-9. [PMID: 17672495 DOI: 10.1021/jp072859+] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Aromatic ring-peptide bond interactions (modeled as benzene and formamide, N-methylformamide and N-methylacetamide) are studied by means of advanced computational chemistry methods: second-order Möller-Plesset (MP2), coupled-cluster single and double excitation model [CCSD(T)], and density functional theory with dispersion (DFT-D). The geometrical preferences of these interactions as well as their interaction energy content, in both parallel and T-shaped arrangements, are investigated. The stabilization energy reaches a value of over 5 kcal/mol for the N-methylformamide-benzene complex at the CCSD(T)/complete basis set (CBS) level. Decomposition of interaction energy by the DFT-symmetry-adapted perturbation treatment (SAPT) technique shows that the parallel and T-shaped arrangements, although similar in their total interaction energies, differ significantly in the proportion of electrostatic and dispersion terms.
Collapse
|
238
|
Riley KE, Hobza P. Assessment of the MP2 method, along with several basis sets, for the computation of interaction energies of biologically relevant hydrogen bonded and dispersion bound complexes. J Phys Chem A 2007; 111:8257-63. [PMID: 17649987 DOI: 10.1021/jp073358r] [Citation(s) in RCA: 146] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In the past several years the MP2 method has been used extensively in studies of noncovalent interactions within biological systems such as proteins, DNA/RNA, and protein-ligand complexes. In this work we assess the performance that can be expected of this method, when paired with several different medium and extended basis sets, for the accurate computation of binding energies of hydrogen bonded and dispersion bound biologically derived complexes. It is found that, overall, the MP2/cc-pVTZ method produces the best, most well balanced, description of noncovalent interactions. Another interesting observation made in this study is that generally the MP2 technique, when paired with any basis set, does not yield reliable results for cyclic hydrogen bonds such as those found in nucleic acid base pairs.
Collapse
|
239
|
Reha D, Hocek M, Hobza P. Exceptional thermodynamic stability of DNA duplexes modified by nonpolar base analogues is due to increased stacking interactions and favorable solvation: Correlated ab initio calculations and molecular dynamics simulations. Chemistry 2007; 12:3587-95. [PMID: 16502452 DOI: 10.1002/chem.200501126] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The geometries of DNA hexamer (5'-GGAACC-3') and DNA 13-mer (5'-GCGTACACATGCG-3') have been determined by molecular dynamics (MD) simulations using an empirical force field. The central canonical base pair was replaced by a pair of nonpolar base analogues, 2,2'-bipyridyl and 3-methylisocarbostyril. The stabilization energy of the model system (model A) consisting of a central base pair (base-analogue pair) and two neighboring base pairs was determined by the RI-MP2 method using an extended aug-cc-pVDZ basis set. The geometry of the model was averaged from structures determined by MD simulations. The role of the solvent was covered by the COSMO continuum solvent model and calculations were performed for a larger model system (model B) which also contained a sugar-phosphate backbone. The total stabilization energies of the unperturbed system and the system perturbed by a base-analogue pair (model A) were comparable to the stability of both duplexes experimentally determined. This is due to large stacking interaction energy of the base-analogue self-pair which compensates for the missing hydrogen-bonding energy of the replaced adenine...thymine base pair. The selectivity of the base-analogue pair was reproduced (model B) when their desolvation energy was included with the interaction energy of both strands determined by the approximate SCC-DFTB-D method.
Collapse
|
240
|
Dabkowska I, Gonzalez HV, Jurecka P, Hobza P. Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels. J Phys Chem A 2007; 109:1131-6. [PMID: 16833422 DOI: 10.1021/jp046738a] [Citation(s) in RCA: 114] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.
Collapse
|
241
|
Abo-riziq A, Crews BO, Compagnon I, Oomens J, Meijer G, Von Helden G, Kabelac M, Hobza P, de Vries MS. The mid-IR spectra of 9-ethyl guanine, guanosine, and 2-deoxyguanosine. J Phys Chem A 2007; 111:7529-36. [PMID: 17625809 DOI: 10.1021/jp072183i] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present the mid-IR (400-1800 cm(-1)) spectra of 9-ethyl guanine, guanosine, and 2-deoxyguanosine measured by IR-UV double-resonance spectroscopy. We compare the recorded mid-IR spectra with the spectra of the most stable structures obtained from RI-MP2 and RI-DFT-D calculations. The results confirm the enol form for all structures and demonstrate the efficacy of a new approach to DFT calculations that includes dispersion interactions.
Collapse
|
242
|
Kabelác M, Valdes H, Sherer EC, Cramer CJ, Hobza P. Benchmark RI-MP2 database of nucleic acid base trimers: performance of different density functional models for prediction of structures and binding energies. Phys Chem Chem Phys 2007; 9:5000-8. [PMID: 17851596 DOI: 10.1039/b707182e] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new database of nucleic acid base trimers has been developed that includes 141 geometries and stabilization energies obtained at the RI-MP2 level of theory with the TZVPP basis set. Compared to previously compiled biologically oriented databases, this new construct includes considerably more complicated structures; the various intermolecular interactions in the trimers are quite heterogeneous and in particular include simultaneous hydrogen bonding and stacking interactions, which is similar to the situation in actual biopolymers. Validation against these benchmark data is therefore a more demanding task for approximate models, since correct descriptions of all energy terms are unlikely to be accomplished by fortuitous cancellations of systematic errors. The density functionals TPSS (both with and without an empirical dispersion term), PWB6K, M05-2X, and BH&H, and the self-consistent charge density functional tight binding method augmented with an empirical dispersion term (SCC-DFTB-D) were assessed for their abilities accurately to compute structures and energies. The best reproduction of the BSSE corrected RI-MP2 stabilization energies was achieved by the TPSS functional (TZVPP basis set) combined with empirical dispersion; removal of the dispersion correction leads to significantly degraded performance. The M05-2X and PWB6K functionals performed very well in reproducing the RI-MP2 geometries, but showed a systematic moderate underestimation of the magnitude of base stacking interactions. The SCC-DFTB-D method predicts geometries in fair agreement with RI-MP2; given its computational efficiency it represents a good option for initial scanning of analogous biopolymeric potential energy surfaces. BH&H gives geometries of comparable quality to the other functionals but significantly overestimates interaction energies other than stacking.
Collapse
|
243
|
Kolár M, Hobza P. Accurate Theoretical Determination of the Structure of Aromatic Complexes Is Complicated: The Phenol Dimer and Phenol···Methanol Cases. J Phys Chem A 2007; 111:5851-4. [PMID: 17566991 DOI: 10.1021/jp071486+] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The structure of the phenol dimer and phenol...methanol complexes was determined by gradient optimization using the Hartree-Fock (HF), MP2, DFT, and RI-DFT-D methods with various basis sets. Theoretical rotational constants were compared with experimental values and the following conclusions were made: (1) HF and DFT methods fail to predict cluster geometries; (2) MP2 with a medium basis set yields reliable cluster geometries but only because of a compensation for errors; (3) when the AO basis set is enlarged, the geometry becomes incorrect, and the theoretical geometry becomes reliable only when the higher correlation energy contributions (CCSD(T)) are included; and (4) the RI-DFT-D procedure covering the dispersion energy provides excellent geometries.
Collapse
|
244
|
Zendlová L, Hobza P, Kabelác M. Potential energy surfaces of the microhydrated guanine...cytosine base pair and its methylated analogue. Chemphyschem 2007; 7:439-47. [PMID: 16463334 DOI: 10.1002/cphc.200500311] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
A complete scan of the potential and free-energy surfaces of monohydrated and dihydrated guanine...cytosine and 9-methylguanine...1-methylcytosine base pairs was realized by the molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The most stable and populated structures localized were further fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method with a large basis set. A systematic study of microhydration of these systems using a high-level correlated ab initio approach is presented for the first time. The different behavior of guanine...cytosine and adenine...thymine complexes is also discussed. These studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for better understanding of the influence of the solvent on the stability of the structure of DNA.
Collapse
|
245
|
Pluhácková K, Hobza P. On the Nature of the Surprisingly Small (Red) Shift in the Halothane⋅⋅⋅Acetone Complex. Chemphyschem 2007; 8:1352-6. [PMID: 17487895 DOI: 10.1002/cphc.200700153] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The halothaneacetone and fluoroformacetone complexes are studied using the second-order Møller-Plesset (MP2) method with a cc-pVTZ basis set and the density functional theory (DFT) method with a TZVP basis set. Whereas halothane exhibits a small red shift upon complexation, fluoroform shows a pronounced blue shift. To explain this difference in behavior, we perform symmetry-adapted perturbation theory (SAPT) and natural bond orbital (NBO) analyses. Although the composition of the total stabilization energy of each complex is different, that alone does not provide a satisfactory explanation for the difference in the spectral shifts. This difference is interpreted as a result of the interplay of the hyperconjugation and rehybridization mechanisms. The small and surprising red shift of the C-H stretching frequency of halothane, which resulted from the complexation of this species with acetone,is explained by the compensation of the two above-mentioned mechanisms. On the other hand, the fluoroformacetone complex exhibits a blue shift of the C--H stretching frequency upon complexation, the most likely reason for this shift being a concerted occurrence of the hyperconjugation and rehybridization mechanisms. The calculated shift of the C-H stretching vibration frequencies of halothane (+27 cm(-1)) agree with the experimental value of +5 cm(-1).
Collapse
|
246
|
Abstract
Non-covalent interactions play an important role in chemistry, physics and especially in biodisciplines. They determine the structure of biomacromolecules such as DNA and proteins and are responsible for the molecular recognition process. Theoretical evaluation of interaction energies is difficult; however, perturbation as well as variation (supermolecular) methods are briefly described. Accurate interaction energies can be obtained by complete basis set limit calculations providing a large portion of correlation energy is covered (e.g. by performing CCSD(T) calculations). The role of H-bonding and stacking interactions in the stabilisation of DNA, oligopeptides and proteins is described, and the importance of London dispersion energy is shown.
Collapse
|
247
|
Kubar T, Hanus M, Ryjácek F, Hobza P. Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations. Chemistry 2007; 12:280-90. [PMID: 16294358 DOI: 10.1002/chem.200500725] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5-ethyl-6-phenylphenanthridinium and uncharged 3,8-diamino-6-phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy. In the case of the cationic intercalators, the electrostatic and charge-transfer terms are also important. The DeltaG of ethidium intercalation to DNA was estimated at -4.5 kcal mol(-1) and this value agrees well with the experimental result. Of six contributions to the final free energy, the interaction energy value is crucial. The intercalation process is governed by the non-covalent stacking (including charge-transfer) interaction while the hydrogen bonding between the ethidium amino groups and the DNA backbone is less important. This is confirmed by the evaluation of the interaction energy as well as by the calculation of the free energy change. The intercalation affects the macroscopic properties of DNA in terms of its flexibility. This explains the easier entry of another intercalator molecule in the vicinity of an existing intercalation site.
Collapse
|
248
|
Muchová E, Slavícek P, Sobolewski AL, Hobza P. Glycine in an Electronically Excited State: Ab Initio Electronic Structure and Dynamical Calculations. J Phys Chem A 2007; 111:5259-69. [PMID: 17530834 DOI: 10.1021/jp071483x] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The goal of this study is to explore the photochemical processes following optical excitation of the glycine molecule into its two low-lying excited states. We employed electronic structure methods at various levels to map the PES of the ground state and the two low-lying excited states of glycine. It follows from our calculations that the photochemistry of glycine can be regarded as a combination of photochemical behavior of amines and carboxylic acid. The first channel (connected to the presence of amino group) results in ultrafast decay, while the channels characteristic for the carboxylic group occur on a longer time scale. Dynamical calculations provided the branching ratio for these channels. We also addressed the question whether conformationally dependent photochemistry can be observed for glycine. While electronic structure calculations favor this possibility, the ab initio multiple spawning (AIMS) calculations showed only minor relevance of the reaction path resulting in conformationally dependent dynamics.
Collapse
|
249
|
Abstract
Gas-phase spectroscopy lends itself ideally to the study of isolated molecules and provides important data for comparison with theory. In recent years, we have seen enormous progress in the study of biomolecular building blocks in the gas phase. The motivation for such work is threefold: (a) It is important to distinguish between intrinsic molecular properties and properties that result from the biological environment. (b) Gas-phase spectroscopy of clusters provides insights into fundamental interactions and into microsolvation. (c) Gas-phase data support quantum-chemical calculations. This review focuses on the current status of (poly)amino acids and DNA bases. Recent results help elucidate structure and hydrogen-bonded interactions, as well as showcase a successful interplay between theory and experiment.
Collapse
|
250
|
Lee EC, Kim D, Jurecka P, Tarakeshwar P, Hobza P, Kim KS. Understanding of Assembly Phenomena by Aromatic−Aromatic Interactions: Benzene Dimer and the Substituted Systems. J Phys Chem A 2007; 111:3446-57. [PMID: 17429954 DOI: 10.1021/jp068635t] [Citation(s) in RCA: 512] [Impact Index Per Article: 30.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Interactions involving aromatic rings are important in molecular/biomolecular assembly and engineering. As a consequence, there have been a number of investigations on dimers involving benzene or other substituted pi systems. In this Feature Article, we examine the relevance of the magnitudes of their attractive and repulsive interaction energy components in governing the geometries of several pi-pi systems. The geometries and the associated binding energies were evaluated at the complete basis set (CBS) limit of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using a least biased scheme for the given data set. The results for the benzene dimer indicate that the floppy T-shaped structure (center-to-center distance: 4.96 A, with an axial benzene off-centered above the facial benzene) is isoenergetic in zero-point-energy (ZPE) corrected binding energy (D0) to the displaced-stacked structure (vertical interplanar distance: 3.54 A). However, the T-shaped structure is likely to be slightly more stable (D0 approximately equal to 2.4-2.5 kcal/mol) if quadruple excitations are included in the coupled cluster calculations. The presence of substituents on the aromatic ring, irrespective of their electron withdrawing or donating nature, leads to an increase in the binding energy, and the displaced-stacked conformations are more stabilized than the T-shaped conformers. This explains the wide prevalence of displaced stacked structures in organic crystals. Despite that the dispersion energy is dominating, the substituent as well as the conformational effects are correlated to the electrostatic interaction. This electrostatic origin implies that the substituent effect would be reduced in polar solution, but important in apolar media, in particular, for assembling processes.
Collapse
|