Copper(II) diamino acid complexes: quantum chemical computations regarding diastereomeric effects on the energy of complexation

[reaction: see text] Quantum chemical calculations were used to rationalize the observed enantiodifferentiation in the complexation of alpha-amino acids to chiral Cu(II) complexes. Apart from Cu(II)[bond]pi interactions and steric repulsions between the anchoring cholesteryl-Glu moiety and an aromatic amino acid R group, hydrogen bonding also plays a role. In fact, in the case of tryptophan, C[double bond]O...H[bond]N hydrogen bonding between the glutamate moiety and the tryptophan N[bond]H group compensates for the loss of intramolecular hydrogen-bonding and diminished Cu(II)[bond]pi interactions.     

Structure,stability, and electronic and NMR properties of various oxo- and nitrido-derivatives of [L(Salen)Mn(III)]+, where L = none and imidazole. A density functional study

Structure, stability, and electronic and NMR properties of [(Salen)Mn(III)](+)-derived intermediates/reactants in the epoxidation/amination of unfunctionalized olefins, namely [(Salen)Mn(V)O](+) (1-oxo), [(Salen)Mn(IV)O] (2-oxo), and [(Salen)Mn(V)N] (3), have been studied with the B3LYP density functional method. It has been shown that the (1)A, (3)A, and (5)A states of cationic 1-oxo species are virtually degenerate, while for the neutral 2-oxo species the ground (4)A state lies 6.4 kcal/mol lower than (2)A. In the nitrido species 3, the (1)A state has been shown to be the ground state in agreement with experiment. We have investigated isomerization of 1-oxo and 2-oxo species into unusual [(OSalen)Mn(III)](+) (1-N-oxo and 1-peroxo) and [(OSalen)Mn(II)] (2-N-oxo and 2-peroxo) species, respectively. For cationic species 1, the 1-N-oxo isomers are more stable (by 10-12 kcal/mol) than the 1-oxo isomer and are separated from the latter by 21-22 kcal/mol barriers. On the other hand, 1-peroxo isomers are calculated to be 14-16 kcal/mol higher than the 1-oxo isomer. For neutral species 2, however, both 2-N-oxo and 2-peroxo isomers lie significantly higher in energy than the 2-oxo isomer. It has been shown that coordination of axial imidazole ligand alters relative energies of spin states for 1- and 2-oxo species, destabilizing low-spin states. For singlet states of H(2)Salen, 1-oxo, and 3, we have calculated (1)H, (13)C, (15)N, and (17)O NMR chemical shifts using the gauge-independent-atomic orbital (GIAO) approach.     

Semiclassical study of collision-induced predissociation: Comparison of the Landau-Zener model with the method of analytic continuation

Semiclassical calculations are carried out by two methods for the problem of collision-induced predissociation of electronically excited I₂. The first method is that of surface-hopping with the Landau-Zener model. The second method is similar to surface-hopping, except that analytic continuation of the adiabatic potential energy surfaces replaces the Landau-Zener model. Results of the calculations by the two methods compare favorably with each other and with experiment. The possible advantages of the second method are discussed.Camille and Henry Dreyfus Teacher-Scholar, Alfred P. Sloan Research Fellow.     

A comparative study of various computational approaches in calculating the structure of pyridoxal 5'-phosphate (PLP)-dependent beta-lyase protein. The importance of protein environment

Various computational approaches, using molecular mechanics (Amber), semiempirical (AM1), density functional (B3LYP), and various ONIOM methods, have been comparatively investigated for the structure of Escherichia coli NifS CsdB protein. The structure of the entire monomer containing 407 amino acid residues and 579 surrounding water molecules has been optimized. The full geometry optimization in the "active site-only" approach (including only active site atoms) has been found to give the largest root-mean-square (RMS) deviation from the X-ray structure; a much better agreement has been achieved by keeping the atoms leading to the backbones of some amino acids frozen in their positions in the X-ray structure. The best agreement has been attained by including the surrounding protein in the calculations using the two-layer ONIOM (B3LYP:Amber) approach. The results presented in this study conclusively demonstrate the importance of the protein/active-site interaction on the active-site structure of the enzyme. The present theoretical study represents the largest system studied at the ONIOM level to date, containing 7992 atoms, including 84 atoms in the QM region and rest in the MM region.     

Ab initio potential energy surfaces of the ion-molecule reaction: C2H2 + O+

High level ab initio calculations using complete active space self-consistent field and multi reference single and double excitation configuration interaction methods with cc-pVDZ (correlation consistent polarized valence double zeta) and cc-pVTZ (triple zeta) basis sets have been performed to elucidate the reaction mechanism of the ion-molecule reaction, C2H2(1Sigmag+) + O+(4S), for which collision experiment has been performed by Chiu et al. [J. Chem. Phys. 109, 5300 (1998)]. The minor low-energy process leading to the weak spin-forbidden product C2H2+ (2Piu) + O(1D) has been studied previously and will not be discussed here. The major pathways to form charge-transfer (CT) products, C2H2+ (2Piu) + O(3P) (CT1) and C2H2+ (4A2) + O(3P) (CT2), and the covalently bound intermediates are investigated. The approach of the oxygen atom cation to acetylene goes over an energy barrier TS1 of 29 kcal/mol (relative to the reactant) and adiabatically leads the CT2 product or a weakly bound intermediate Int1 between CT2 products. This transition state TS1 is caused by the avoided crossing between the reactant and CT2 electronic states. As the C-O distance becomes shorter beyond the above intermediate, the C1 reaction pathway is energetically more favorable than the Cs pathway and goes over the second transition state TS2 of a relative energy of 39 kcal/mol. Although this TS connects diabatically to the covalent intermediate Int2, there are many states that interact adiabatically with this diabatic state and these lead to the other charge-transfer product CT1 via either of several nonadiabatic transitions. These findings are consistent with the experiment, in which charge transfer and chemical reaction products are detected above 35 and 39 kcal/mol collision energies, respectively.     

Multi-configuration electron-hole potential method for excited states

The recently proposed electron-hole potential (EHP) method for excited states is extended to the multi-configurational case. The variation equation is solved using the quadratic convergence method. The EHP methods are shown to be approximations to the complete singly excited configuration interaction (CSECI) in the variational sense. Extended Brillouin theorems are proved for the EHP methods. The excitation energies and wave functions obtained by one and two configurational EHP methods agree well with those of the CSECI method. The EHP methods have clear advantage in the computer time requirement over the CI method and are especially suited for a calculation of approximate excited states of large molecules. The EHP methods are applicable to excited states which belong to the same irreducible representation as the ground state.     

Computational study of the aminolysis of esters. The reaction of methylformate with ammonia

The aminolysis of esters is a basic organic reaction considered as a model for the interaction of carbonyl group with nucleophiles. In the present computational study the different possible mechanistic pathways of the reaction are reinvestigated by applying higher level electronic structure theory, examining the general base catalysis by the nucleophile, and a more comprehensive study the solvent effect. Both the ab initio QCISD/6-31(d,p) method and density functional theory at the B3LYP/6-31G(d) level were employed to calculate the reaction pathways for the simplest model aminolysis reaction between methylformate and ammonia. Solvent effects were assessed by the PCM method. The results show that in the case of noncatalyzed aminolysis the addition/elimination stepwise mechanism involving two transition states and the concerted mechanism have very similar activation energies. However, in the case of catalyzed aminolysis by a second ammonia molecule the stepwise mechanism has a distinctly lower activation energy. All transition states in the catalyzed aminolysis are 10-17 kcal/mol lower than those for the uncatalyzed process.     

Can chlorine anion catalyze the reaction of HOCl with HCl?

The reaction of HOCl + HCl → Cl₂ + H₂O in the presence of chlorine anion Cl⁻ has been studied using ab initio methods. The overall exothermicity is 15.5 kcal mol⁻¹ and this reaction has been shown to have a high activation barrier of 46.5 kcal mol⁻¹. Cl⁻ is found to catalyze the reaction via the formation of HOCl·Cl⁻, ClH·HOCl·Cl⁻ and Cl⁻·H₂) intermediate ion-molecule complexes or by interacting with a concerted four-center transition state of the reaction of HOCl + HCl.     

Regulation of geometry around the ruthenium center of bis(2-pyridinecarboxylato) complexes by the nitrosyl moiety: syntheses, structures, and theoretical studies

cis-[Ru(NO)Cl(pyca)(2)] (pyca = 2-pyridinecarboxylato), in which the two pyridyl nitrogen atoms of the two pyca ligands coordinate at the trans position to each other and the two carboxylic oxygen atoms at the trans position to the nitrosyl ligand and the chloro ligand, respectively (type I shown as in Chart 1), reacted with NaOCH(3) to generate cis-[Ru(NO)(OCH(3))(pyca)(2)] (type I). The geometry of this complex was confirmed to be the same as the starting complex by X-ray crystallography: C(13.5)H(13)N(3)O(6.5)Ru; monoclinic, P2(1)/n; a = 8.120(1), b = 16.650(1), c = 11.510(1) A; beta = 99.07(1) degrees; V = 1536.7(2) A(3); Z = 4. The cis-trans geometrical change reaction occurred in the reactions of cis-[Ru(NO)(OCH(3))(pyca)(2)] (type I) in water and alcohol (ROH, R = CH(3), C(2)H(5)) to form [[trans-Ru(NO)(pyca)(2)](2)(H(3)O(2))](+) (type V) and trans-[Ru(NO)(OR)(pyca)(2)] (type V). The reactions of the trans-form complexes, trans-[Ru(NO)(H(2)O)(pyca)(2)](+) (type V) and trans-[Ru(NO)(OCH(3))(pyca)(2)] (type V), with Cl(-) in hydrochloric acid solution afforded the cis-form complex, cis-[Ru(NO)Cl(pyca)(2)] (type I). The favorable geometry of [Ru(NO)X(pyca)(2)](n)(+) depended on the nature of the coexisting ligand X. This conclusion was confirmed by theoretical, synthetic, and structural studies. The mono-pyca-containing nitrosylruthenium complex (C(2)H(5))(4)N[Ru(NO)Cl(3)(pyca)] was synthesized by the reaction of [Ru(NO)Cl(5)](2)(-) with Hpyca and characterized by X-ray structural analysis: C(14)H(24)N(3)O(3)Cl(3)Ru; triclinic, Ponemacr;, a = 7.631(1), b = 9.669(1), c = 13.627(1) A; alpha = 83.05(2), beta = 82.23(1), gamma = 81.94(1) degrees; V = 981.1(1) A(3); Z = 2. The type II complex of cis-[Ru(NO)Cl(pyca)(2)] was synthesized by the reaction of [Ru(NO)Cl(3)(pyca)](-) or [Ru(NO)Cl(5)](2)(-) with Hpyca and isolated by column chromatography. The structure was determined by X-ray structural analysis: C(12)H(8)N(3)O(5)ClRu; monoclinic, P2(1)/n; a = 10.010(1), b = 13.280(1), c = 11.335(1) A; beta = 113.45(1) degrees; V = 1382.4(2) A(3); Z = 4.