Conformations of the thiathiophthenes and related molecules

A conformational ab initio MO study has been carried out for the thiathiophthene molecule (TTP) and two related model compounds, thiomalonaldehyde (TMA) and its conjugate base (TMA(-)). The conformational energy surfaces for TMA, TMA(-) and TTP were generated using a least squares fit to the calculated data and plotted on a CALCOMP plotter. The results of the calculations showed that the cis-cis planar conformation of TTP is the most stable in agreement with experimental findings. For TMA and TMA(-) the cis-cis planar conformation is not the most stable. Contour plots of the five occupied -MO's of TTP show great similarity to those of naphthalene.Less detailed calculations were carried out for 3-hydroxy-prop-2-en-1-thione (HPT) and 3-mercapto-prop-2-en-1-thione (MPT). HPT was shown to be most stable in the cis planar hydrogen bonded conformation in agreement with the experimentally obtained results. For MPT the non-hydrogen bonded planar structure was found to be the most stable.     

Molecular orbital studies on the conformation of GABA (γ-aminobutyric acid)

In distinction to Extended Hückel Theory which predicts as the most stable conformation of free zwitterionic GABA a totally extended form, PCILO and SCF ab initio studies show that the intrinsically preferred conformation of the isolated molecule is a highly folded one, resulting from strong interactions between the two charged ends. Computations are also carried out for hydrated GABA in the supermolecule approach allowing moreover for the flexibility of binding of some of the water molecules of the first hydration shell. They predict the coexistence in solution of a large number of conformations showing different degrees of folding (or extension), a result confirmed by recent NMR studies. This and a number of similar results show that we have to adapt our thinking on the role of conformations in pharmacological activity to this situation, which was frequently obscured by the more abundant results of X-ray crystallography yielding a single conformation.     

Excited states of phenyl carbonyl compounds

In an attempt to clarify the origin of the dual phosphorescence in phenyl alkyl ketones, we have made some calculation (within the C.I.P.S.I. method in an excitonic scheme) to elucidate the conformation of both ground states and excited states of propiophenone. Our calculations have shown the presence of two stable isomers in the ground state, first n ^* state, and first ^* singlet and triplet states. So our work suggests that the origin of the dual phosphorescence of propiophenone could be related to the conformational change of the molecule in the n ^* state, because the most stable conformations in the n ^* state and in the ground state are different.     

Band structure of even and ions of odd polyenes

The analysis of experimental data for singlet transitions (E _n) of even polyenes (I), cations (II) and anions (III) of odd polyenes show that for infinite chains E (I)/E (II)=E (I)/E (III) = 2:1. It is shown that the energy gap is equal for the three systems. In cases (II) and (III) there is a level (NBMO) in the gap which is vacant in (II) and occupied in (III). That is why the first optical transition in (II) and (III) depends on the semiwidth of the gap.     

Bilirubin conformational enantiomer selection in sodium deoxycholate chiral micelles

Intramolecularly hydrogen-bonded, bichromophoric tetrapyrrole pigments, bilirubin-IX and mesobilirubin-XIII, adopt either of two enantiomeric conformations which are in dynamic equilibrium in solution. InpH 8 aqueous sodium deoxycholate solutions, chiral micelles preferentially select one conformational enantiomer, and the solutions exhibit a bisignate circular dichroism Cotton effect in the vicinity of the bilirubin long wavelength electronic transition. Exciton coupling theory indicates a predominance of the left-handed (or negative) chirality bilirubin conformational enantiomer.     

Ab initio study of hydrazoic acid

SCF and CI calculations were carried out on the ground^1A state of HN₃. The equilibrium geometry and vibration frequencies were computed. The results point to a planar structure (groupC _s) but to a non-linear (170 °) N-N-N conformation. The calculated vibration frequencies are in fair agreement with experimental assignments.The dissociation path of the molecule to NH and N₂ products was investigated and compared to the isoelectronic reaction of diazomethane. The dissociation energy of hydrazoic acid is estimated to be about –8 kcal/mole, with a potential barrier to dissociation of about 30 kcal/mole.Boursier IRSIA     

Cation binding to biomolecules

SCF ab initio computations are carried out on the binding of alkali and alkaline-earth cations to the phosphate monoanion. The effect of the binding on the conformational properties of the phosphodiester linkage and of the polar head of phospholipids is investigated. The results indicate that following the nature of the cation and the site of its binding, the interaction may have a profound influence on the conformation of the ligand. The consequences of this situation on the use of lanthanide probes in NMR studies are considered.     

Studies of solvent effects

The effect of water on the conformational preferences of acetylcholine has been studied within the discrete, the continuum and the combined discrete-continuum models described in parts I and II of this series. All the models lead to the conclusion that the trans-gauche form which is, following refined quantum-mechanical computations, the intrinsically preferred one and the one observed in the crystal of acetylcholine and of a number of analogues should remain also the preferred conformation in water. This result agrees with NMR studies. The results of the empirical discrete model used here compare favorably to those obtained by an ab initio super-molecule treatment. The continuum model utilized here represents a net improvement above such models utilized in other works.     

On the mathematical formulation and interpretation of LACO MO theory

The mathematical basis of LCAO MO theory is studied, both within the Hartree-Fock approximation and in more exact formulations. The basic LCAO expansion for molecular orbitals ¦> in terms of atomic orbitals ¦x> is conveniently written ¦> = ¦x> S ^–1 B where S is the overlap matrix for atomic orbitals and B is the matrix of atomic orbital-molecular orbital overlaps. It is suggested that matrices P and Q, defined by P=B B and Q=BnB where n is the matrix of molecular orbital occupation numbers, are appropriate to the interpretation of molecular calculations in terms of atomic orbital components, electronic populations and the degree of bonding. Implications for Hartree-Fock calculations are investigated.     

Theory of isomer shift in hemin

Using the self-consistent charge extended Hückel procedure, the charge density difference at Fe⁵⁷ nucleus, between hemin and Fe⁺³ ion is calculated. This is combined with the recent value of the calibration constant, –0.23±0.02a ₀ ³ mm/sec to obtain an isomer shift of –0.374 mm/sec between hemin and Fe⁺³ in good agreement with the value –0.392 mm/sec derived from experimental data and the calculated value of the isomer shift of Fe⁺³ with respect to K₃FeF₆ from first principle covalency investigations in the latter compound. is composed of contributions from core and valence electrons of the same order of magnitude, with the latter being more than one-half of the former. The core contribution is composed of a number of terms of comparable magnitude and differing signs, whose significance is discussed.Supported by grant HL 15196-02 from the Heart and Lung Institute of National Institute of Health.     

Bonding criteria for diatomic molecular orbitals and inter-relations among them

Bonding criteria for molecular orbitals in diatomic molecules are discussed. An orbital force criterion is shown to have several conceptual and practical advantages, providing a basis for the investigation of inter-relations among many of the commonly employed criteria.It is found that interconsistency among those criteria is guaranteed, within the framework of Koopmans' Theorem, if the orbital energies are monotonic in the range (R _e, ).The application of the orbital force criterion to the second row homonuclear diatomics exhibits reasonable chemical trends concerning the valence-shell orbitals, as well as indications of a slightly antibonding nature of the inner orbitals.Based on a section of a thesis to be submitted by Y.T. to the Senate of the Technion-Israel Institute of Technology, in partial fulfilment of the requirements for the D.Sc. degree     

Theoretical studies on the conformation of saccharides

Conformations of 2-methoxytetrahydropyran as a model for the six-membered ring in aldopyranosides have been calculated by the PCILO method using the algorithm of the conjugated gradient to optimize the geometry. The calculated geometry of the fourteen basic forms of 2-methoxytetrahydropyran was found to be in agreement with the available data obtained by X-ray diffraction of pyranosides. The results indicate differences in the geometry of 2-methoxytetrahydropyran resulting from the change of the axial vs. equatorial position of the methoxyl group. These changes are particularly meaningful in the values of bond angles and they are in agreement with the anomeric and exoanomeric effects. The experimentally found differences in the energies of an axial (⁴ C ₁) and equatorial (¹ C ₄) conformer, G = 2.9–3.7 kJ/mol, and the dipole moment, = 1.20 ± 0.05 D (1D = 3.33 10^–30mAs) agree well with the calculated values E = 3.18 kJ/mol and <> = 1.18 D which, in turn, suggest that the axial conformer is preferred over the equatorial one by a ratio a:e = 78:22.     

An ab initio study of formamide

Calculated energy and molecular properties of the ground and low-energy excited states of formamide are presented at the ground state geometry. Satisfactory results are obtained except for the ¹* energy which remains too high by 1 eV (which is nevertheless a large improvement over previous calculations). The predicted triplet energies lie at 5.4 eV (³ n*) and 5.8 eV (³*).     

Charge-transfer complexes betweenp-toluidine and iodine in solution: A kinetic study

The kinetics of charge-transfer interaction betweenp-toluidine and iodine in methylene chloride was investigated in depth. Thethermal process of formation of theinner complex was found to proceed to an equilibrium. Thephotochemical process follows a different reaction coordinate, going through the formation of an exciplex between the excitedouter complex and the amine ground state. In both cases the same ionic complex (Am ₂I⁺I ₃ ^– , whereAm stands forp-toluidine) was detected as the final product.

---

Kinetische Untersuchung der Charge-Transfer-Komplexe zwischenp-Toluidin und Jod

Zusammenfassung Die Kinetik der Charge-Transfer-Wechselwirkung zwischenp-Toluidin und Jod in Methylenchlorid wurde ausführlich untersucht. Derthermische Prozeß, der zur Bildung desinner-Komplexes führt, geht bis zu einem Gleichgewicht. Derphotochemische Prozeß folgt einer unterschiedlichen Reaktionskoordinate und verläuft über die Bildung eines Exziplexes zwischen dem angeregtenouter-Komplex und dem Amin im Grundzustand. In beiden Fällen wurde derselbe ionische Komplex (Am ₂I⁺I ₃ ^– , wobeiAm fürp-Toluidin steht) als Endprodukt festgestellt.

     

Linear combination of hybrid orbitals: Acyclic alkanes

A LCHO SCF-CI semiempirical scheme is formulated for compounds of carbon and hydrogen, and far-UV transition energies computed for a number of acyclic alkanes. The type of spectral bands (CH CH*, etc.) and symmetries of subbands are assigned, and their dependence on conformation and branching discussed. It is argued that, for pentanes and higher alkanes, spectral quantities can be estimated from data on the lower alkanes.Multi-Conformational Compounds, part X. For part IX, see Ref. [1].     

SCF perturbation theory for unimolecular reactions

A CNDO/2 SCF perturbation theory is presented for interpreting the form of CNDO/2 potential energy surfaces of unimolecular reactions. The analysis is performed by calculating the energy change E arising from a distortion of the molecular geometry along the reaction coordinate. E is decomposed into different perturbational contributions which are appropriate for an interpretation of the perturbation energy E. Moreover, E is resolved into energy parts arising from a single occupied orbital and contributions due to pairwise orbital interactions. In this way one evaluates numerically how the form of the occupied and unoccupied orbitals determines the magnitude of E. If the distortion occurs along a definite symmetry coordinate, group-theoretical arguments can be applied to discuss the magnitude of characteristic components of the perturbation energy. The SCF perturbation theory is used to analyze the isomerization of ethylene, cis-2-butene and cis-2-butenenitrile.This work was partially supported by Nato-Grant No. 1072     

Ab initio SCF and CI study of the electronic spectrum of pyridine N-oxide

Configuration interaction (CI) studies of ground, n *, * * electronically excited states are reported for pyridine N-oxide. The transition energy to the lowest * excited ¹ B ₂ state is calculated at 4.35 eV, compared to the experimental spectrum range of 3.67–4.0 eV. This state lies below the lowest n * excited ¹ A ₂ state calculated at 4.81 eV above the ground state. The only experimentally reported triplet state at 2.92 eV above the ground state is predicted to be the ³ A ₁ (*) state. The calculated energy lies at 3.27 eV. Numerous other high-lying singlet states as well as the triplet states have also been calculated. The intramolecular charge transfer character of the ground and the excited states have been studied in terms of the calculated dipole moment and other physical properties.     

Conformational analysis of dimethylphosphate with quantum-mechanical and classical methods

In order to shed light on the conformational behavior of polynucleotide chains, and in particular to clarify the origins of the barriers to internal rotation in the phosphodiester linkage, we computed, with a quantum-mechanical ab initio procedure, the energies associated to 86 combinations of the two torsion angles in the dimethylphosphate anion (CH₃O)₂PO₂ ^–, and then we sought for an analytical expression apt to reproduce these energies with the highest possible accuracy. An excellent agreement (standard deviation of the fitted energies from the ab initio energies 0.28 kcal/mole) with the quantum-mechanical calculations was reached with a potential consisting of four terms: 1) a 6–12 Lennard-Jones contribution, in which different parameters are used to describe the interactions of methyls with the ester oxygens and with the anionic oxygens; 2) a contribution with twofold periodicity, accounting for the anomeric effects connected to the interactions between the lone pair electrons and the polar bonds of phosphorus with the anionic oxygens; 3) a contribution with threefold periodicity, representing the usual bond-staggering term; and 4) a Coulombic contribution, arising from electrostatic interactions between partially charged atoms.     

Acidity scale for strong Bronsted acids referred to their complexes with water and in the Gaseous phase

CNDO/2 calculations on AH...B complexes in the gaseous phase (AH = strong acids, B = water) gives the following order H₂O < HF < HCl < HSO₃CF₃ < H₂SO₄ < HSO₃F < HClO₄ < H₃O⁺ if we take as reference the Wiberg residual indexp(H... O). Other properties and in particular the charge transfer behave in the same fashion. Charge variations for each atom give an explanation for the Pauling rule for the oxygenated acids.