On modeling the pressure-dependent photoisomerization of trans-stilbene by including slow intramolecular vibrational energy redistribution

Experimental data for the photoisomerization of trans-stilbene (S(1)) in thermal bath gases at pressures up to 20 bar obtained previously by Meyer, Schroeder, and Troe (J. Phys. Chem. A 1999, 103, 10528-10539) are modeled by using a full collisional-reaction master equation that includes non-RRKM (Rice-Ramsperger-Kassel-Marcus) effects due to slow intramolecular vibrational energy redistribution (IVR). The slow IVR effects are modeled by incorporating the theoretical results obtained recently by Leitner et al. (J. Phys. Chem. A 2003, 107, 10706-10716), who used the local random matrix theory. The present results show that the experimental rate constants of Meyer et al. are described to within about a factor of 2 over much of the experimental pressure range. However, a number of assumptions and areas of disagreement will require further investigation. These include a discrepancy between the calculated and experimental thermal rate constants near zero pressure, a leveling off of the experimental rate constants that is not predicted by theory and which depends on the identity of the collider gas, the need to use rate constants for collision-induced IVR that are larger than the estimated total collision rate constants, and the choice of barrier-crossing frequency. Despite these unsettled issues, the theory of Leitner et al. shows great promise for accounting for possible non-RRKM effects in an important class of reactions.     

Estimates of the energy of intramolecular hydrogen bonds

A method for the estimation of the energy of intramolecular hydrogen bonds in conjugated systems existing in a variety of conformations is presented. The method is applied to determine the intramolecular hydrogen bond energy in 3-aminopropenal and 3-aminopropenthial. According to the proposed estimation scheme, the intramolecular H-bond energies are found to be of the order of 5-7 kcal/mol. These results are compared with those obtained by using other estimation schemes as well as with the recent results by other authors. Also, the H-bond energies in dimers and trimers of the two molecules are calculated and compared with the corresponding data for internally hydrogen-bonded monomers. This comparison shows that the bond equalization effect is primarily due to proton donor-proton acceptor proximity. In comparison with intermolecular hydrogen bonds, the rigidity of the chelate skeleton enhances this proximity effect. The same effect can be seen in systems with intermolecular hydrogen bonds, although its magnitude is diminished because of the absence of additional forces which pull the proton donor and proton acceptor groups toward each other. No specific resonance-assisted origin of the intramolecular hydrogen bond energy seems to be needed to elucidate the energetics of these bonds.     

On the electrostatic energy of charged alkali clusters

In the framework of the jellium model an algorithm is presented to calculate the electrostatic energy of charged metal clusters of arbitrary shape. Explicit expressions are given for some special shapes.     

Quantum effects in intramolecular energy transfer

Anharmonic potentials which are presumed to be classically quasi-ergodic, but which have symmetries leading to degenerate quantum states, fail in the quantum case to transfer energy equivalently among rigorously equivalent phase space locations. This is shown using simple group theory and is illustrated for the case of the Henon-Heiles potential     

Non-covalent switch for intramolecular energy transfer

The axially coordinated complex of phenylazopyridine and a Zn-porphyrin/free-base porphyrin conjugate provides a switch for intramolecular energy transfer, with reversible complexation/decomplexation as a switching protocol.     

A repertoire of pyridinium-phenyl-methyl cross-talk through a cascade of intramolecular electrostatic interactions

Direct intramolecular cation-pi interaction between phenyl and pyridinium moieties in 1a(+) has been experimentally evidenced through pH-dependent (1)H NMR titration. The basicity of the pyridinyl group (pK(a) 2.9) in 1a can be measured both from the pH-dependent chemical shifts of the pyridinyl protons as well as from the protons of the neighboring phenyl and methyl groups as a result of electrostatic interaction between the phenyl and the pyridinium ion in 1a(+) at the ground state. The net result of this nearest neighbor electrostatic interaction is that the pyridinium moiety in 1a becomes more basic (pK(a) 2.92) compared to that in the standard 2a (pK(a) 2.56) as a consequence of edge-to-face cation (pyridinium)-pi (phenyl) interaction, giving a free energy of stabilization (DeltaDeltaG(o)pKa) of -2.1 kJ mol(-1). The fact that the pH-dependent downfield shifts of the phenyl and methyl protons give the pK(a) of the pyridine moiety of 1a also suggests that the nearest neighbor cation (pyridinium)-pi (phenyl) interaction also steers the CH (methyl)-pi (phenyl) interaction in tandem. This means that the whole pyridine-phenyl-methyl system in 1a(+) is electronically coupled at the ground state, cross-modulating the physicochemical property of the next neighbor by using the electrostatics as the engine, and the origin of this electrostatics is a far away point in the molecule-the pyridinyl-nitrogen. The relative chemical shift changes and the pK(a) differences show that the cation (pyridinium)-pi (phenyl) interaction is indeed more stable (DeltaDeltaG(o)pKa = -2.1 kJ mol(-1)) than that of the CH (methyl)-pi (phenyl) interaction (DeltaDeltaG(o)pKa = -0.8 kJ mol(-1)). Since the pK(a) of the pyridine moiety in 1a is also obtained through the pH-dependent shifts of both phenyl and methyl protons, it suggests that the net electrostatic mediated charge transfer from the phenyl to the pyridinium and its effect on the CH (methyl)-pi (phenyl) interaction corresponds to DeltaG(o)pKa of the pyridinium ion (approximately 17.5 kJ mol(-1)), which means that the aromatic characters of the phenyl and the pyridinium rings in 1a(+) have been cross-modulated owing to the edge-to-face interaction proportional to this DeltaG(o)pKa change.     

Characterization and modeling of multi-element electrostatic lens systems

Electrostatic cylinder lenses are widely used to control beams of charged particles with different energies and directions in several fields, especially in electron spectroscopy. In this study, we investigated the zoom-lens properties of the three-, four- and five-element electrostatic lens systems as a function of the lens voltages and their dimensions. Lens systems, which consist of cylindrical electrodes, each spaced 0.1 diameter apart, were designed by using the SIMION and LENSYS programs. Three different lens configurations were studied to form an image at a specific position for use in experimental studies. We also discussed the line-shape profile of a three-element lens system as a function of the focusing voltage for einzel lens configuration.     

Computational experiments on the intramolecular energy flow in macromolecules

The microscopic details of the flow of energy in a single chain of polyethylene containing 300 atoms is discussed. The intramolecular dynamics of the polyethylene molecule is studied as a function of CH stretch excitation, temperature, and pressure. The rate of energy flow from CH stretching modes is found to be very rapid and irreversible, occurring on a timescale of less than 0.5 ps at low temperatures, and increases with temperature. A general characteristic two-phase energy flow behavior is observed, where there is initially a very rapid flow (due to the decay of the initial excitation) followed by a slower flow (due to energy redistribution throughout the system). The mechanism for the initial facile energy flow is shown to involve strong resonant pathways. In particular, a CH stretch/HCH bend Fermi (1:2) resonance is shown to dominate the short-time dynamics and facilitates the overall process of energy redistribution. The increase in the rate of energy flow as a function of the backbone temperature is found to be due to the increase in the density of the bath states for energy redistribution which subsequently results in the formation of new low-order resonant interactions (1:1, etc). The long-time dynamics, associated to complete redistribution of the initial CH stretch energy with all of the 894 available vibrational modes, occurs within a time of 2 ps. This timescale corresponds to the time for intramolecular redistribution. A comparison of the intramolecular redistribution time to that of intermolecular redistribution (redistribution in the condensed or solid phase as opposed to a single chain) is also made. A preliminary study of energy flow in a crystal of polyethylene (system containing 19 polyethylene chains) shows that the energy flow exhibits two very different time behaviors. The first is for the intramolecular redistribution as in the single chain study and the second is for intermolecular (chain-to-chain) redistribution. The timescale for intermolecular redistribution is found to be on the order of 0.2 ns at room temperature and pressure, about two orders of magitude larger than the intramolecular timescale.     

Accurate calculation of the binding energy of the water dimer

The binding energy of the water dimer is calculated at the MP2 level using efficient basis sets augmented with bond functions. The intermolecular energy is determined by the supermolecular approach and the basis set superposition error is corrected by the counterpoise method. Bond functions are found useful and very effective in recovering the dispersion energy, which is traditionally achieved by polarization functions. The calculated binding energy of the water dimer is systematically converged to a value of 4.75 kcal mol⁻ as bond functions are gradually added to nucleus-centered basis sets.     

Heavy mass barrier to intramolecular energy transfer

The dynamics of a collection of seven Morse or harmonic oscillators are investigated to model a molecule in which two halves are separated by a heavy atom. The results are related to a recent experiment on intramolecular dynamics and suggest an extension of the anharmonic local mode concept to groups.     

Estimating the energy of intramolecular hydrogen bonds in chitosan oligomers

The effect the number of chitosan monomer units CTS_n (n = 1–5), the protonation of chitosan dimers, and the interaction between CTS_n (n = 1–3) and acetate ions have on the energy of intramolecular hydrogen bonds is investigated by means of QTAIM analysis and solving the vibrational problem within the cluster-continuum model. It is established that the number of H-bonds in CTS_n is 2n ? 1 and the total energy of H-bonds grows by ~20 kJ/mol. It is concluded that the hydrogen bonds between CTS and acetate ions play a major role in the stabilization of polyelectrolyte complexes in dilute acetic acid solutions of CTS.     

Accurate and interpretable computational modeling of chemical mutagenicity

We describe a method for modeling chemical mutagenicity in terms of simple rules based on molecular features. A classification model was built using a rule-based ensemble method called RuleFit, developed by Friedman and Popescu. We show how performance compares favorably against literature methods. Performance was measured through the use of cross-validation and testing on external test sets. All data sets used are publicly available. The method automatically generated transparent rules in terms of molecular structure that agree well with known toxicology. While we have focused on chemical mutagenicity in demonstrating this method, we anticipate that it may be more generally useful in modeling other molecular properties such as other types of chemical toxicity.     

Analytical calculation of the electrostatic interaction energy within the CNDO framework

This work analyses the adequacy of an analytical electrostatic energy formulation within the CNDO framework to predict the stable conformations of large molecular complexes. Comparisons are made with abinitio results for small systems such as water-formamide, methanol-water-imidazole, or guanine-cytosine and with AM 1 results for two large systems: a molecular tweezer + the 9-methyladenine complex and a model active site of the α-chymotrypsin and its ligand complex. This approach is efficient in providing reliable confromers for large molecular systems in a very fast way. © 1993 John Wiley & Sons, Inc.     

On the question of heavy-atom blocking of intramolecular vibrational energy transfer

The possibility that a heavy atom can block internal energy transfer in excited species has been further investigated by chemical activation study of 4-(trimethyl lead)-2-butyl and 5-(trimethyl tin)-2-pentyl radicals. These radicals were prepared and excited by the reaction of H atoms with the corresponding 1-olefins. Energy randomization was found to occur on a subpicosecond time scale in all the radicals studied.     

Accuracy assessment of semiempirical molecular electrostatic potential of proteins

 Accurate electrostatic maps of proteins are of great importance in research of protein interaction with ligands, solvent media, drugs, and other biomolecules. The large size of real-life proteins imposes severe limitations on computational methods one can use for obtaining the electrostatic map. Well-known accurate second-order M?ller–Plesset and density functional theory methods are not routinely applicable to systems larger than several hundred atoms. Conventional semiempirical tools, as less resource demanding ones, could be an attractive solution but they do not yield sufficiently accurate calculation results with reference to protein systems, as our analysis demonstrates. The present work performs a thorough analysis of the accuracy issues of the modified neglect of differential overlap type semiempirical Hamiltonians AM1 and PM3 on example of the calculation of the molecular electrostatic potential and the dipole moment of natural amino acids. Real capabilities and limitations of these methods with application to protein modeling are discussed. Received: 26 April 2002 / Accepted: 19 September 2002 / Published online: 14 February 2003     

Controlling the DNA binding specificity of bHLH proteins through intramolecular interactions

Reversible control of the conformation of proteins was employed to probe the relationship between flexibility and specificity of the basic helix-loop-helix protein MyoD. A fusion protein (apaMyoD) was designed where the basic DNA binding helix of MyoD was stablized by an amino-terminal extension with a sequence derived from the bee venom peptide apamin. The disulfide-stabilized helix from apamin served as a nucleus for a helix that extended for a further ten residues, thereby holding apaMyoD's DNA recognition helix in a predominantly alpha-helical conformation. The thermal stability of the DNA complexes of apaMyoD was increased by 13 degrees C relative to MyoD-bHLH. Measurements of the fluorescence anisotropy change on DNA binding indicated that apaMyoD bound to E-box-containing DNA sequences with enhanced affinity relative to MyoD-bHLH. Consequently, the DNA binding specificity of apaMyoD was increased 10-fold.     

Accurate free energy calculation along optimized paths

The path‐based methods of free energy calculation, such as thermodynamic integration and free energy perturbation, are simple in theory, but difficult in practice because in most cases smooth paths do not exist, especially for large molecules. In this article, we present a novel method to build the transition path of a peptide. We use harmonic potentials to restrain its nonhydrogen atom dihedrals in the initial state and set the equilibrium angles of the potentials as those in the final state. Through a series of steps of geometrical optimization, we can construct a smooth and short path from the initial state to the final state. This path can be used to calculate free energy difference. To validate this method, we apply it to a small 10‐ALA peptide and find that the calculated free energy changes in helix‐helix and helix‐hairpin transitions are both self‐convergent and cross‐convergent. We also calculate the free energy differences between different stable states of β‐hairpin trpzip2, and the results show that this method is more efficient than the conventional molecular dynamics method in accurate free energy calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A scheme estimating the energy of intramolecular hydrogen bonds in diols

The relative energies of conformers of 1,2-ethanediol, 1,3-propanediol, and 1,4-butanediol are split into a sum of five different terms including the intramolecular OH?O interaction. This scheme allows to estimate the energy of the O-H?O intramolecular hydrogen bond of the tGG′g and gGG′g conformers of 1,3-propanediol, the g′GG′Gt and g′GG′Gg conformers of 1,4-butanediol, and the energy of the non-bonded O-H?O interaction in the g′Gt, g′Gg and g′Gg′ conformers of 1,2-ethanediol. This scheme provides pure hydrogen bond energies without assuming the geometry and/or electronic features to be constant between the conformation having a IHB and a reference conformation. The fitted energies show a perfect linear correlation with the corresponding r(H?O)⁻¹ values. QTAIM atomic electron population and energies of the donor hydrogen calculated along the H-O-C-C internal rotation are found to be linearly correlated. These linear correlations display small changes at the BCP formation in 1,3-propanediol.