Wittig reactions on photoprotoporphyrin IX: new synthetic models for the special pair of the photosynthetic reaction center

A first example of spirochlorin-chlorin dimer with fixed distances and orientations as potential model for the "special pair" of the photosynthetic reaction center is discussed. For the preparation of such a novel structure, the Wittig reagent of the desired "spacer" 5 was reacted with photoprotoporphyrin IX dimethyl ester 3 to produce the intermediate dimer 6, which on intramolecular [4 + 2] Diels-Alder cycloaddition gave an unexpected spirochlorin-chlorin dimer 9. Dehydration of dimer 6 under acid-catalyzed conditions generated the corresponding spirochlorin-porphyrin dimer 16 in quantitative yield. The asymmetry in dimer 6 caused by the biphenyl-type anisotropic effect was confirmed by NMR and model studies. The formation of dihydrobenzoporphyrin 14 by reacting chlorin 3 with the phosphonium salt of p-methylbenzylbromide 10 and isolation of 8-phenanthrenevinylporphyrin 19 from chlorin 7 further confirmed our proposed mechanism for the formation of a spirochlorin-chlorin dimer 9. Following a similar approach, chlorin 3 on reacting with bis-phosphonium salt of 4, 4'-bischloromethylbiphenyl produced conjugated chlorin dimer 25. The spectroscopic data obtained from these dimers suggest that, in these compounds, the individual chromophores are not behaving as an individual molecule, but as a single macrocycle. To examine whether the pi-pi interaction exhibited by dimer 9 resembles the structural arrangement of bacteriochlorophylls in reaction center (RC), we investigated the geometrical parameters used to characterize the pi-pi interactions in tetrapyrrolic macrocycles. Starting from the crystallographic coordinates of 9, the molecular mechanics energy minimization was performed to obtain the model dimer structure. The geometrical parameters that measure the single pyrrole ring overlap were used to compare the model structure with the crystallographic coordinates of the special pair in photosynthetic reaction center. The results indicated that the ring A of spirochlorin and the ring C of chlorin in our model dimer 9 mimic the ring A-ring A interaction found in the crystallographic special pairs, which are strategically placed by the surrounding photosynthetic reaction center protein matrix.     

Kinetics and mechanism of electron-induced splitting of a cyclobutane pyrimidine dimer with or without an electron acceptor

DNA repair has received heightened attention in recent years as ozone depletion threatens to significantly increase DNA damage by UVB radiation[1—6]. The major lesions formed in DNA by this radiation are cis-syn cyclobutane pyrimidine dimers, which are created by the linkage of two neighboring pyrimidine bases in DNA via C5-C5 and C6-C6 atoms by [2+2] cycloaddition[2,5—8]. This potentially lethal or mutagenic damage can be repaired either by the removal of the damaged bases by excisio…     

Parallel iterative reaction path optimization in ab initio quantum mechanical/molecular mechanical modeling of enzyme reactions

The determination of reaction paths for enzyme systems remains a great challenge for current computational methods. In this paper we present an efficient method for the determination of minimum energy reaction paths with the ab initio quantum mechanical/molecular mechanical approach. Our method is based on an adaptation of the path optimization procedure by Ayala and Schlegel for small molecules in gas phase, the iterative quantum mechanical/molecular mechanical (QM/MM) optimization method developed earlier in our laboratory and the introduction of a new metric defining the distance between different structures in the configuration space. In this method we represent the reaction path by a discrete set of structures. For each structure we partition the atoms into a core set that usually includes the QM subsystem and an environment set that usually includes the MM subsystem. These two sets are optimized iteratively: the core set is optimized to approximate the reaction path while the environment set is optimized to the corresponding energy minimum. In the optimization of the core set of atoms for the reaction path, we introduce a new metric to define the distances between the points on the reaction path, which excludes the soft degrees of freedom from the environment set and includes extra weights on coordinates describing chemical changes. Because the reaction path is represented by discrete structures and the optimization for each can be performed individually with very limited coupling, our method can be executed in a natural and efficient parallelization, with each processor handling one of the structures. We demonstrate the applicability and efficiency of our method by testing it on two systems previously studied by our group, triosephosphate isomerase and 4-oxalocrotonate tautomerase. In both cases the minimum energy paths for both enzymes agree with the previously reported paths.     

New insights into the ultrafast photophysics of oxidized and reduced FAD in solution

The ultrafast photophysics of oxidized and reduced flavin adenine dinucleotide (FAD) in aqueous solution was studied by broadband UV-vis femtosecond transient absorption spectroscopy. We observed that oxidized FAD (FAD(ox)) in solution readily aggregates at submillimolar concentration. Upon excitation of FAD(ox), three excited-state lifetimes were found and assigned to three different species: the closed (stacked) conformation of the monomer (～5.4 ps), the open (extended) conformation of the monomer (～2.8 ns), and the dimer (～27 ps). In the case of the stacked conformation of the monomer, we show that intramolecular electron transfer from the adenine to the isoalloxazine ring occurs with a time constant of 5.4 ps and is followed by charge recombination on a faster time scale, namely, 390 fs. We additionally demonstrate that deprotonated reduced flavin (FADH(-)) undergoes biphotonic ionization under high excitation fluence and dissociates into a hydrated electron and the neutral semiquinone radical FADH(?).     

Coordinate reduction for exploring chemical reaction paths

The potential energy surface for the reaction of a typical molecular system composed of N atoms is defined uniquely by 3N-6 coordinates. These coordinates can be defined by the Cartesian coordinates of the atomic centers (minus overall translation and rotation), or a set of internally defined coordinates such as bond stretches, angle bends, and torsions. By applying principal component analysis to the geometries along a reaction path, a reduced set of coordinates, d ≪ 3N-6, can be obtained. This reduced set of coordinates can reproduce the changes in geometry along the reaction path with chemical accuracy and may help improve the efficiency of reaction path optimization algorithms.     

A new module for constrained multi‐fragment geometry optimization in internal coordinates implemented in the MOLCAS package

A parallel procedure for an effective optimization of relative position and orientation between two or more fragments has been implemented in the MOLCAS program package. By design, the procedure does not perturb the electronic structure of a system under the study. The original composite system is divided into frozen fragments and internal coordinates linking those fragments are the only optimized parameters. The procedure is capable to handle fully independent (no border atoms) fragments as well as fragments connected by covalent bonds. In the framework of the procedure, the optimization of relative position and orientation of the fragments are carried out in the internal “Z‐matrix” coordinates using numerical derivatives. The total number of required single points energy evaluations scales with the number of fragments rather than with the total number of atoms in the system. The accuracy and the performance of the procedure have been studied by test calculations for a representative set of two‐ and three‐fragment molecules with artificially distorted structures. The developed approach exhibits robust and smooth convergence to the reference optimal structures. As only a few internal coordinates are varied during the procedure, the proposed constrained fragment geometry optimization can be afforded even for high level ab initio methods like CCSD(T) and CASPT2. This capability has been demonstrated by applying the method to two larger cases, CCSD(T) and CASPT2 calculations on a positively charged benzene lithium complex and on the oxygen molecule interacting to iron porphyrin molecule, respectively. © 2013 Wiley Periodicals, Inc.     

Quantitative fluorescence kinetic analysis of NADH and FAD in human plasma using three- and four-way calibration methods capable of providing the second-order advantage

The metabolic coenzymes reduced nicotinamide adenine dinucleotide (NADH) and flavin adenine dinucleotide (FAD) are the primary electron donor and acceptor respectively, participate in almost all biological metabolic pathways. This study develops a novel method for the quantitative kinetic analysis of the degradation reaction of NADH and the formation reaction of FAD in human plasma containing an uncalibrated interferent, by using three-way calibration based on multi-way fluorescence technique. In the three-way analysis, by using the calibration set in a static manner, we directly predicted the concentrations of both analytes in the mixture at any time after the start of their reactions, even in the presence of an uncalibrated spectral interferent and a varying background interferent. The satisfactory quantitative results indicate that the proposed method allows one to directly monitor the concentration of each analyte in the mixture as the function of time in real-time and nondestructively, instead of determining the concentration after the analytical separation. Thereafter, we fitted the first-order rate law to their concentration data throughout their reactions. Additionally, a four-way calibration procedure is developed as an alternative for highly collinear systems. The results of the four-way analysis confirmed the results of the three-way analysis and revealed that both the degradation reaction of NADH and the formation reaction of FAD in human plasma fit the first-order rate law. The proposed methods could be expected to provide promising tools for simultaneous kinetic analysis of multiple reactions in complex systems in real-time and nondestructively.     

The electronic structure of the flavin cofactor in DNA photolyase.

Density functional theory is used to calculate the electronic structure of the neutral flavin radical, FADH(*), formed in the light-induced electron-transfer reaction of DNA repair in cis,syn-cyclobutane pyrimidine dimer photolyases. Using the hybrid B3LYP functional together with the double-zeta basis set EPR-II, (1)H, (13)C, (15)N, and (17)O isotropic and anisotropic hyperfine couplings are calculated and explained by reference to the electron densities of the highest occupied molecular orbital and of the unpaired spin distribution on the radical. Comparison of calculated and experimental hyperfine couplings obtained from EPR and ENDOR/TRIPLE resonance leads to a refined structure for the FAD cofactor in Escherichia coli DNA photolyase. Hydrogen bonding at N3H, O4, and N5H results in significant changes in the unpaired spin density distribution and hyperfine coupling constants. The calculated electronic structure of FADH(*) provides evidence for a superexchange-mediated electron transfer between the cyclobutane pyrimidine dimer lesion and the 7,8-dimethyl isoalloxazine moiety of the flavin cofactor via the adenine moiety.     

Calculation of the O-H stretching vibrational overtone spectrum of the water dimer

The O-H stretching vibrational overtone spectrum of the water dimer has been calculated with the dimer modeled as two individually vibrating monomer units. Vibrational term values and absorption intensities have been obtained variationally with a computed dipole moment surface and an internal coordinate Hamiltonian, which consists of exact kinetic energy operators within the Born-Oppenheimer approximation of the monomer units. Three-dimensional ab initio potential energy and dipole moment surfaces have been calculated using the internal coordinates of the monomer units using the coupled cluster method including single, double, and perturbative triple excitations [CCSD(T)] with the augmented correlation consistent valence triple zeta basis set (aug-cc-pVTZ). The augmented correlation consistent valence quadruple zeta basis set (aug-cc-pVQZ), counterpoise correction, basis set extrapolation to the complete basis set limit, relativistic corrections, and core and valence electron correlations effects have been included in one-dimensional potential energy surface cuts. The aim is both to investigate the level of ab initio and vibrational calculations necessary to produce accurate results when compared with experiment and to aid the detection of the water dimer under atmospheric conditions.     

Counterpoise-corrected potential energy surfaces of simple H-bonded systems

Geometries and stabilization energies of various simple H-bonded complexes (water dimer, hydrogen fluoride dimer, formamide dimer, formic acid dimer) have been determined by a gradient optimization that eliminates the basis set superposition error (BSSE) by the counterpoise (CP) method in each gradient cycle as well as by the standard gradient optimization. Both optimization methods lead to different potential energy surfaces (PES). The difference depends on the theoretical level used and is larger if correlation energy is considered. Intermolecular distances from the CP-corrected PES are consistently longer, and this difference might be significant (∼0.1 ?); also angular characteristics determined from both surfaces differ significantly. Different geometries were obtained even when passing to larger basis sets (aug-cc-pVDZ). The standard optimization procedure can result in a completely wrong structure. For example, the “quasi-linear” structure of the (HF)₂ (global minimum) does not exist at the standard MP2/ 6-31G** PES (where only cyclic structure was detected) and is found only at the CP-corrected PES. Stabilization energies obtained from the CP-corrected PES are always larger than these from the standard PES where the BSSE is added only a posteriori for the final optimized structure; both energies converge only when passing to a larger basis set (aug-cc-pVDZ). Received: 11 March 1998 / Accepted: 19 June 1998 / Published online: 4 September 1998 RID=" ID=" <E6>Acknowledgements.</E6> The project was supported by the Grant Agency of the Czech Republic (Grant No. 203/98/1166). RID=" ID=" <E5>Correspondence to</E5>: P. Hobza     

Synthesis of "sugar-rods" with phytohemagglutinin cross-linking properties by using the palladium-catalyzed Sonogashira reaction

A palladium-catalyzed Sonogashira reaction has been applied for the syntheses of divalent "sugar-rods" which exhibited excellent lectin cross-linking properties. The procedure, which involves a tetrakis(triphenylphosphine)palladium-catalyzed cross-coupling reaction between an alkyne and a halogen-bearing sp2-carbon in DMF at 60 degrees C, is very efficient and the dimeric or heterobifunctional "sugar-rods" (8-13, 15-17) were isolated in 65-100% yields. Dimers 8a and 15a were both shown to form insoluble cross-linked lattices when mixed with the tetrameric plant lectin from Canavalia ensiformis (Concavalin A, Con A). Moreover, the relative inhibitory properties of the synthetic dimannosides were determined by means of the hemagglutination of rabbit erythrocytes, whereby dimer 15a was shown to be 20-fold more potent than monomeric methyl alpha-D-mannopyranoside.     

Partial optimization of large molecules and clusters

By examining the displacement coordinate metric three modes of constrained optimization for large molecules and clusters are suggested. The first method corresponds to a conventional optimization using internal coordinates. The second mode has applications with respect to both internal and cartesian coordinates. The final mode is particularly interesting because it can result in computational savings. A mixture of both internal and cartesian coordinates is specified where these coordinates are usually a subset of the molecules or clusters total coordinate set. In the optimization only a subset of the energy derivatives need be evaluated reducing the computational effort associated with the gradient calculation.     

Density-functional geometry optimization of the 150,000-atom photosystem-I trimer

We present a linear-scaling method based on the use of density-functional theory (DFT) for the system-wide optimization of x-ray structural coordinates and apply it to optimize the 150,000 atoms of the photosystem-I (PS-I) trimer. The method is based on repetitive applications of a multilevel ONIOM procedure using the PW916-31G(d) DFT calculations for the high level and PM3 for the lower level; this method treats all atoms in the structure equivalently, a structure in which the majority of the atoms can be considered as part of some internal "active site." To obtain a realistic single structure, some changes to the original protein model were necessary but these are kept to a minimum in order that the optimized structure most closely resembles the original x-ray one. Optimization has profound effects on the perceived electronic properties of the cofactors, with, e.g., optimization lowering the internal energy of the chlorophylls by on average 53 kcal mol(-1) and eliminates the enormous 115 kcal mol(-1) energy spread depicted by the original x-ray heavy-atom coordinates. A highly precise structure for PS-I results that is suitable for analysis of device function. Significant qualitative features of the structure are also improved such as correction of an error in the stereochemistry of one of the chlorophylls in the "special pair" of the reaction center, as well as the replacement of a water molecule with a metal cation in a critical region on the C3 axis. The method also reveals other unusual features of the structure, leading both to suggestions concerning device functionality and possible mutations between gene sequencing and x-ray structure determination. The optimization scheme is thus shown to augment the molecular modeling schemes that are currently used to add medium-resolution structural information to the raw scattering data in order to obtain atomically resolved structures. System-wide optimization is now a feasible process and its use within protein x-ray data refinement should be considered.     

Vibrational calculations in formaldehyde: the CH stretch system

An alternative procedure for the calculation of highly excited vibrational levels in S0 formaldehyde was developed to apply to larger molecules. It is based on a new set of symmetrized vibrational valence coordinates. The fully symmetrized vibrational kinetic energy operator is derived in these coordinates using the Handy expression [Molec. Phys. 61, 207 (1987)]. The potential energy surface is expressed as a fully symmetrized quartic expansion in the coordinates. We have performed ab initio electronic computations using GAMESS to obtain all force constants of the S0 formaldehyde quartic force field. Our large scale vibrational calculations are based on a fully symmetrized vibrational basis set, in product form. The vibrational levels are calculated one by one using an artificial intelligence search/selection procedure and subsequent Lanczos iteration, providing access to extremely high vibrational energies. In this work special attention has been given to the CH stretch system by calculating the energies up to the fifth CH stretch overtone at ∼16000 cm⁻¹, but the method has also been tested on two highly excited combination levels including other lower frequency modes.      

The optimization of MCSCF functions by application of the generalized brillouin theorem: The LiH2 potential energy surface

The generalized Brillouin theorem is used to construct an optimization procedure for MCSCF functions by iterative contracted CI calculations. Special attention is paid to the MO transformation step in each iteration. In this method the MCSCF calculation may easily be augmented by a restricted CI calculation involving a configuration set which is uniquely determined by the trial function. An application to the calculation of the potential energy surface for linear LiH₂ in the reaction LiH + HLi + H₂ leads to the conclusion that this restricted CI is necessary, in order to obtain satisfactory results for the potential energy barrier in this reaction.     

化学反应的内禀反应坐标法(Ⅲ)——甲硫醛反式异构化反应的IRC解析

本文用内禀反应坐标(IRC)方法,采用Gaussian程序,4-31G基函数,研究了甲硫醛反式异构化反应机理。     

Accurate theoretical determination of the structure of aromatic complexes is complicated: the phenol dimer and phenol...methanol cases

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.     

Estimating the hessian for gradient-type geometry optimizations

Optimization methods that use gradients require initial estimates of the Hessian or second derivative matrix; the more accurate the estimate, the more rapid the convergence. For geometry optimization, an approximate Hessian or force constant matrix is constructed from a simple valence force field that takes into account the inherent connectivity and flexibility of the molecule. Empirical rules are used to estimate the diagonal force constants for a set of redundant internal coordinates consisting of all stretches, bends, torsions and out-of-plane deformations involving bonded atoms. The force constants are transformed from the redundant internal coordinates to Cartesian coordinates, and then from Cartesian coordinates to the non-redundant internal coordinates used in the specification of the geometry and optimization. This method is especially suitable for cyclic molecules. Problems associated with the choice of internal coordinates for geometry optimization are also discussed.Fellow of the Alfred P. Sloan Foundation, 1981–83     

On the Importance of CP-corrected Gradient Optimization in the Study of Hydrogen Bonded Systems

Geometries, harmonic vibrational frequencies and interaction energies of the water-hydrogen sulfide dimer, hydrogen fluoride dimer and glycine zwitterion-water dimer were determined by the counterpoise-corrected (CP-corrected) gradient optimization that explicitly corrects for the basis set superpusition error (BSSE) and CP-uncorrected (normal) gradient opfimization respectively at the B3LYP and MP2 levels of theory, employing the popular Pople‘s standard 6-31G(d), 6-31G(d,p) and 6-311 G(d,p) basis sets in order to assess the importance of CP-corrected gradient optimiTation in the study of hydrogen bonded systems. The normal optimization of these three H-bonded systems obtained using these popular basis sets all yielded erratic results, whereas use of CP-corrected gradient optimization led to consistent results with those from larger basis sets. So this CP receipt becomes useful and necessary to correctly describe large systems, where the use of small basis sets may be necessary.     

Single‐ended transition state finding with the growing string method

Reaction path finding and transition state (TS) searching are important tasks in computational chemistry. Methods that seek to optimize an evenly distributed set of structures to represent a chemical reaction path are known as double‐ended string methods. Such methods can be highly reliable because the endpoints of the string are fixed, which effectively lowers the dimensionality of the reaction path search. String methods, however, require that the reactant and product structures are known beforehand, which limits their ability for systematic exploration of reactive steps. In this article, a single‐ended growing string method (GSM) is introduced which allows for reaction path searches starting from a single structure. The method works by sequentially adding nodes along coordinates that drive bonds, angles, and/or torsions to a desired reactive outcome. After the string is grown and an approximate reaction path through the TS is found, string optimization commences and the exact TS is located along with the reaction path. Fast convergence of the string is achieved through use of internal coordinates and eigenvector optimization schemes combined with Hessian estimates. Comparison to the double‐ended GSM shows that single‐ended method can be even more computationally efficient than the already rapid double‐ended method. Examples, including transition metal reactivity and a systematic, automated search for unknown reactivity, demonstrate the efficacy of the new method. This automated reaction search is able to find 165 reaction paths from 333 searches for the reaction of NH₃BH₃ and (LiH)₄, all without guidance from user intuition. © 2015 Wiley Periodicals, Inc.