Conformation dependence of pKa’s of the chromophores from the purple asFP595 and yellow zFP538 fluorescent proteins

Two members of the green fluorescent protein family, the purple asFP595 and yellow zFP538 proteins, are perspective fluorescent markers for use in multicolor imaging and resonance energy-transfer applications. We report the results of quantum based calculations of the solution pK_a values for selected protonation sites of the denatured asFP595 and zFP538 chromophores in the trans- and cis-conformations in order to add in the interpretation of photo-physical properties of these proteins. The pK_a values were determined from the theromodynamic cycle based on B3LYP/6-311++G(2df, 2p) calculations of the gas phase free energies of the molecules and the B3LYP/6-311++G(d, p) calculations of solvation energies. The results show that the pK_a’s of the protonation sites of the chromophore from asFP595 noticeably depend on the isomer conformation (cis- or trans-), while those of zFP538 are much less sensitive to isomerization.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

Gas electron diffraction study of the geometric structure of triallylborane molecule

The main structural parameters of the triallylborane molecule having the C ₃ symmetry were determined by gas electron diffraction and quantum-chemical calculations at the MP2/6-31G(d,p) and B3LYP/6-31G(d,p) levels. The parameters calculated by the MP2/6-31G(d,p) method are in better agreement with the experimental data than those calculated by the B3LYP/6-31G(d,p) method.__________Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 98–101, January, 2005.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

FT-IR and FT-Raman spectra and vibrational investigation of 4-chloro-2-fluoro toluene using ab initio HF and DFT (B3LYP/B3PW91) calculations

FT-IR (4000-100 cm(-1)) and FT-Raman (4000-100 cm(-1)) spectra of solid sample of 4-chloro-2-fluoro toluene (4Cl2FT) have been recorded using Bruker IFS 66 V spectrometer. Ab initio-HF (HF/6-311++G (d, p)) and DFT (B3LYP/6-311++G and B3PW91/6-311++G (d, p)) calculations have been performed giving energies, optimized structures, harmonic vibrational frequencies, depolarization ratios, IR intensities, Raman activities. The vibrational frequencies are calculated and scaled values are compared with FT-IR and FT-Raman experimental values. The isotropic HF and DFT analyses showed good agreement with experimental observations. The differences between the observed and scaled wave number values of most of the fundamentals are very small in B3LYP than HF. Comparison of the simulated spectra provides important information about the ability of the computational method (B3LYP) to describe the vibrational modes. The influences of substitutions on the geometry of molecule and its normal modes of vibrations have also been discussed. The changes made by substitutions on the benzene are much responsible for the non-linearity of the molecule. This is an attractive entity for the future studies of non-linear optics.     

Conformational and electronic properties of N-methacryloyl-(L)-glutamic acid

Conformers of N-methacryloyl-(L)-glutamic acid (MAGA, C₉H₁₃NO₅) in the gas phase, benzene or methanol environment were searched by the density functional theory (DFT) and time-dependent DFT. Selected conformational isomers based on potential energy surface analysis were studied by considering the structural and electronic properties. B3LYP functional, and HF and MP2 levels were used with the 6-31+G(d,p), 6-311+G(3df,p) and aug-cc-pVDZ basis sets. We explored the effect of medium on the conformational preferences, UV spectra and frontier molecular orbitals of the conformers. The outcomes of this research will be useful for future studies including moieties analogous to MAGA.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

Electronegativity and resonance parameters from the geometry of monosubstituted benzene rings

Reference values of the structural substituent parameters, S _E and S _R, measuring the electronegativity and resonance effects, respectively, of functional groups (Campanelli et al. J Phys Chem A 107:6429–6440, 2003) have been determined from the benzene ring geometries of 100 Ph–X species, including different conformations of the same molecule. Geometries have been obtained by quantum chemical calculations at the HF/6-31G*, HF/6-311++G**, and B3LYP/6-311++G** levels of theory. The substituent parameters from HF/6-311++G** calculations are in close agreement with those determined at the HF/6-31G* level. Using the B3LYP density functional yields S _E and S _R values which—in general—correlate well with the corresponding HF values. Exceptions occur with some charged groups, and, in the case of S _E, with a few dipolar groups having very high or low electronegativities. S _R values from B3LYP calculations are about 22% smaller than the corresponding HF values. The variations of the benzene ring geometry caused by electronegativity, resonance, and steric effects are illustrated in some detail.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Evaluation of <Superscript>13</Superscript>C NMR spectra of cyclopropenyl and cyclopropyl acetylenes by theoretical calculations

A convenient methodology was developed for a very accurate calculation of ¹³C NMR chemical shifts of the title compounds. GIAO calculations with density functional methods (B3LYP, B3PW91, PBE1PBE) and 6-311+G(2d,p) basis set predict experimental chemical shifts of 3-ethynylcyclopropene (1), 1-ethynylcyclopropane (2) and 1,1-diethynylcyclopropane (3) with high accuracy of 1–2 ppm. The present article describes in detail the effect of geometry choice, density functional method, basis set and effect of solvent on the accuracy of GIAO calculations of ¹³C NMR chemical shifts. In addition, the particular dependencies of ¹³C chemical shifts on the geometry of cyclopropane ring were investigated.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

Theoretical study and atoms in molecule analysis of hydrogen bonded clusters of ammonia and isocyanic acid

Abstract  

Ab initio and density functional calculations were used to analyze the interaction between a molecule of the isocyanic acid with 1 up to 4 molecules of ammonia at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) computational levels. The cooperative effect is increased with the increasing size of studied clusters. Red shifts of the H–N stretching frequency for complexes involving the isocyanic acid as an H-donor were predicted. Atom in molecules was used to analyze cooperative effects on topological parameters.     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Density functional theory and Hartree–Fock-density functional theory calculations of O, S, and Ge quadrupole coupling constants

The performance of several density functional theory and Hartree–Fock density functional theory methods in conjunction with Pople type bases for the calculation of ¹⁷O, ³³S, and ⁷³Ge quadrupole coupling constants in gaseous state molecules was investigated.Assessment of the several models was made by linear regression analysis of the calculated gradient of the molecular electric field versus the experimental nuclear quadrupole coupling constant (NQCCs). Calculations for oxygen on six molecules with the B3LYP/6-311++G(3df,3p) model yield a residual standard deviation of 0.057 MHz (1.4%); for sulfur on 12 molecules with the B3LYP/6-311G(3df,3p) model, 0.42 MHz (1.8%); and for germanium on nine molecules with the B3P86/6-311G(2d) model, 0.83 MHz (1.0%).In the case of germyl acetylene, our calculations indicate that the experimental NQCC reported some time ago by Thomas and Laurie [J. Chem. Phys. 44 (1966) 2602] was incorrectly assigned with respect to algebraic sign.Predictions are made of the ¹⁷O and ³³S NQCCs in furan, 4H-pyran-4-one, 4H-pyran-4-thione, and 4H-thiapyran-4-thione; and of the ⁷³Ge NQCC in germyl bromide.     

Quantum chemical calculation studies on 4-phenyl-1-(1-phenylethylidene)thiosemicarbazide

Ab initio calculations of the structure, atomic charges, natural bond orbital and thermodynamic functions have been performed at HF/6-311G^** and B3LYP/6-311G^** levels of theory for the title compound of 4-phenyl-1-(1-phenylethylidene) thiosemicarbazide. The calculated results show that the sulfur atom and all of the nitrogen atoms have bigger negative charges and that they are the potential sites to react with the metallic ions, which make the title compound a multidentate ligand. The coordination ability of the sulfur atom and the nitrogen atom of C=N double bond increases with the increase of polarity of the solvent. Electronic absorption spectra have been calculated by time-dependent density functional theory (TD-DFT) method. The calculation of the second-order optical nonlinearity has also been carried out with the PM3 semi-empirical method, resulting in the molecular hyperpolarisability is 5.477×10⁻³⁰ esu.     

Synthese et Etude Conformationnelle Par RMN (1H, 13C, 31P) et DFT des Cycloalcoxyphosphinallenes et des Hydrazones β-Cycloalcoxyphosphonatees

Abstract

The addition of hydrazine and its derivatives to cycloalkoxyphosphinallenes leads to β-(5,5-dialkyl-2-oxo-1,3,2-dioxaphosphoranyl)-hydrazones in good yields. The structure of the obtained compounds were elucidated by the NMR (¹H, ¹³C, ³¹P) spectroscopy and density functional theory (DFT) calculations at B3LYP/6-311++G (2d, 2p) level of theory.     

Further evidence for the inseparability of the theoretical enthalpic terminal and penultimate unit effects in the radical copolymerization of styrene with acrylonitrile

The theoretical enthalpies of propagation reactions at 0 K without zero‐point vibrational energy corrections according to terminal and penultimate models of the radical copolymerization of styrene with acrylonitrile are reported from molecular orbital calculations at the following levels of theory and basis sets: HF/6‐31G(d); B3‐LYP/6‐31G(d); B3‐LYP/6‐311G(d,p) and B3‐LYP/6‐311+G(3df)//6‐311G(d,p). Both the enthalpic terminal and penultimate unit effects, determined according to the theoretical thermochemistry, depend on the level of theory and basis set used for the molecular orbital calculations. The best performing B3LYP/6‐311+G(3df)//B3LYP/6‐311G(d,p) procedure gives theoretical enthalpies for the addition of styrene and acrylonitrile to CH that differ from experimental values by 0.6 and 1.6 kcal mol^?1, respectively. An analysis of the results obtained here leads to the conclusion that at least for the styrene–acrylonitrile monomer system, that is, a monomer system known as one of the few systems that do not conform to terminal model composition and microstructure equations, the enthalpic terminal unit effects seem to depend on the penultimate units of the growing radical. This finding, together with the outcome from our previous work on the dependence of the penultimate effects on the terminal units in a growing macroradical, indicates the inseparability of the enthalpic terminal (implicit) and explicit penultimate unit effects on the radical copolymerization. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 1778–1787, 2003     

Mechanism of formation and monomolecular decomposition of aci-nitromethanes: a quantum-chemical study

A quantum-chemical study of the reactions of formation of aci-nitromethane (aci-NM) and aci-dinitromethane (aci-DNM) and their decomposition with elimination of water was carried out. The methods employed were the ab initio RHF method with inclusion of electron correlation at the MP2 level of theory and the Dunning—Hay double zeta basis set augmented with polarization d-functions on heavy-element atoms, the DFT approach at the B3LYP level, and the semiempirical PM3 method. The formation of aci-NM and aci-DNM was found to be the limiting stage of the mechanism under study. For DNM, the barrier to reaction is substantially lower than for NM. The estimates of the heights of the barriers to formation found from density functional calculations at the B3LYP/6-311++G(df,p) level (258 kJ mol^–1 for aci-NM and 218.5 kJ mol^–1 for aci-DNM) are thought to be the most reliable.