M. Waroquier

Normal modes for large molecules with arbitrary link constraints in the mobile block Hessian approach

A. Ghysels, D. Van Neck, B.R. Brooks, V. Van Speybroeck, M. Waroquier
Journal of Chemical Physics
130 (8), 084107
2009
A1

Abstract 

In a previous paper [ Ghysels et al., J. Chem. Phys. 126, 224102 (2007) ] the mobile block Hessian (MBH) approach was presented. The method was designed to accurately compute vibrational modes of partially optimized molecular structures. The key concept was the introduction of several blocks of atoms, which can move as rigid bodies with respect to a local, fully optimized subsystem. The choice of the blocks was restricted in the sense that none of them could be connected, and also linear blocks were not taken into consideration. In this paper an extended version of the MBH method is presented that is generally applicable and allows blocks to be adjoined by one or two common atoms. This extension to all possible block partitions of the molecule provides a structural flexibility varying from very rigid to extremely relaxed. The general MBH method is very well suited to study selected normal modes of large macromolecules (such as proteins and polymers) because the number of degrees of freedom can be greatly reduced while still keeping the essential motions of the molecular system. The reduction in the number of degrees of freedom due to the block linkages is imposed here directly using a constraint method, in contrast to restraint methods where stiff harmonic couplings are introduced to restrain the relative motion of the blocks. The computational cost of this constraint method is less than that of an implementation using a restraint method. This is illustrated for the α-helix conformation of an alanine-20-polypeptide. © 2009 American Institute of Physics

Mobile Block Hessian Approach with Adjoined Blocks: An Efficient Approach for the Calculation of Frequencies in Macromolecules

A. Ghysels, V. Van Speybroeck, E. Pauwels, D. Van Neck, B.R. Brooks, M. Waroquier
Journal of Chemical Theory and Computation (JCTC)
5 (5), 1203-1215
2009
A1

Abstract 

In an earlier work, the authors developed a new method, the mobile block Hessian (MBH) approach, to accurately calculate vibrational modes for partially optimized molecular structures [ J. Chem. Phys. 2007, 126 (22), 224102.]. It is based on the introduction of blocks, consisting of groups of atoms, that can move as rigid bodies. The internal geometry of the blocks need not correspond to an overall optimization state of the total molecular structure. The standard MBH approach considers free blocks with six degrees of freedom. In the extended MBH approach introduced herein, the blocks can be connected by one or two adjoining atoms, which further reduces the number of degrees of freedom. The new approach paves the way for the normal-mode analysis of biomolecules such as proteins. It rests on the hypothesis that low-frequency modes of proteins can be described as pure rigid-body motions of blocks of consecutive amino acid residues. The method is validated for a series of small molecules and further applied to alanine dipeptide as a prototype to describe vibrational interactions between two peptide units; to crambin, a small protein with 46 amino acid residues; and to ICE/caspase-1, which contains 518 amino acid residues.

Nucleophile-dependent regioselective ring opening of 2-substituted N,N-dibenzylaziridinium ions: bromideversushydride

S. Young Yun, S. Catak, W. Koo Lee, M. D'Hooghe, N. De Kimpe, V. Van Speybroeck, M. Waroquier, Y. Kim, H-J. Ha
Chemical Communications
(18), 2508-2510
2009
A1

Abstract 

The ring opening of 2-substituted N,N-dibenzylaziridinium ions by bromide exclusively occurs at the substituted aziridine carbon atom in a stereospecific way, whereas the opposite regioselectivity was observed for hydride-induced ring opening at the unsubstituted position; furthermore, this unprecedented hydride-promoted reactivity was validated by means of Density Functional Theory (DFT) calculations.

ENDOR and HYSCORE analysis and DFT-assisted identification of the third major stable radical in sucrose single crystals X-irradiated at room temperature

H. De Cooman, E. Pauwels, H. Vrielinck, E. Sagstuen, S. Van Doorslaer, F. Callens, M. Waroquier
Physical Chemistry Chemical Physics (PCCP)
11 (7), 1105-1114
2009
A1

Abstract 

Recently, the chemical structure of two of the three major stable radicals (T2 and T3) produced in sucrose single crystals by X-irradiation at room temperature was identified by comparing Density Functional Theory (DFT) calculations of Electron Magnetic Resonance parameters with experimental results [H. De Cooman, E. Pauwels, H. Vrielinck, E. Sagstuen, F. Callens and M. Waroquier, J. Phys. Chem. B, 2008, 112, 7298–7307]. Ambiguities concerning an unusual proton hyperfine coupling (HFC) tensor prevented the identification of the third major stable radical (T1). In the present work, experimental results of continuous wave Electron Nuclear Double Resonance experiments on sucrose single crystals and Hyperfine Sublevel Correlation Spectroscopy experiments on sucrose powder are presented that lift these remaining ambiguities. Using the final set of experimental HFC tensors and employing advanced DFT calculations, the chemical structure of the T1 radical is established: an allylic-type radical with approximately half of the spin density localised on the C2′ carbon of the fructose unit, involving glycosidic bond cleavage at the fructose side and a concerted formation of a carbonyl group at the C1′ carbon. The electronic structure of the T1 radical is discussed in more detail by means of additional DFT calculations, yielding a better understanding of the peculiar properties of the unusual proton HFC tensor mentioned above.

Bond Dissociation Enthalpies of Large Aromatic Carbon-Centered Radicals

K. Hemelsoet, V. Van Speybroeck, M. Waroquier
Journal of Physical Chemistry A
112 (51), 13566-13573
2008
A1

Abstract 

Carbon−hydrogen bond dissociation enthalpy (BDE) values are computed for the class of benzylic radicals. An extended and representative set of large methylated polyaromatics has been submitted to an accurate computational study using various levels of theory. The hybrid B3P86 as well as two contemporary functionals (BMK and M05-2X) are applied. For a selection of species, the suitability of the DFT methods is validated through comparison with high-level G3(MP2)-RAD and SCS-ROMP2 results. The influence of the polyaromatic environment on the BDE results is thoroughly discussed. The results are compared with other hydrocarbon radical types in order to obtain a generalized radical stability scale. In order to complete this investigation, also carbon−carbon BDE values have been calculated, giving information about the influence of the local environment on removing the methyl group from the polyaromatic.

MD-TRACKS: A Productive Solution for the Advanced Analysis of Molecular Dynamics and Monte Carlo simulations

T. Verstraelen, M. Van Houteghem, V. Van Speybroeck, M. Waroquier
Journal of Chemical Information and Modeling (JCIM)
48 (12), 2414–2424
2008
A1

Abstract 

In this paper, we present MD-TRACKS, an advanced statistical analysis toolkit for Molecular Dynamics and Monte Carlo simulations. The program is compatible with different molecular simulation codes, and the analysis results can be loaded into spreadsheet software and plotting tools. The analysis is performed with commands that operate on a binary trajectory database. These commands process not only plain trajectory data but also the output of other MD-TRACKS commands, which enables complex analysis work flows that are easily programmed in shell scripts. The applicability, capabilities, and ease of use of MD-TRACKS are illustrated by means of examples, that is, the construction of vibrational spectra and radial distribution functions from a molecular dynamics run is discussed in the case of tetrahydrofuran. These properties are compared with the experimental data available in the literature. MD-TRACKS is open-source software distributed at http://molmod.ugent.be/code/.

Schonland ambiguity in the electron nuclear double resonance analysis of hyperfine interactions: Principles and practice

H. Vrielinck, H. De Cooman, M.A. Tarpan, E. Sagstuen, M. Waroquier, F. Callens
Journal of Magnetic Resonance
195 (2), 196-205
2008
A1

Abstract 

For the analysis of the angular dependence of electron paramagnetic resonance (EPR) spectra of low-symmetry centres with S = 1/2 in three independent planes, it is well-established—but often overlooked—that an ambiguity may arise in the best-fit tensor result. We investigate here whether a corresponding ambiguity also arises when determining the hyperfine coupling (HFC) tensor for nuclei with I = 1/2 from angular dependent electron nuclear double resonance (ENDOR) measurements. It is shown via a perturbation treatment that for each set of MS ENDOR branches two best-fit tensors can be derived, but in general only one unique solution simultaneously fits both. The ambiguity thus only arises when experimental data of only one MS multiplet are used in analysis or in certain limiting cases. It is important to realise that the ambiguity occurs in the ENDOR frequencies and therefore the other best-fit result for an ENDOR determined tensor depends on various details of the ENDOR experiment: the MS state of the fitted transitions, the microwave frequency (or static magnetic field) in the ENDOR measurements and the rotation planes in which data have been collected. The results are of particular importance in the identification of radicals based on comparison of theoretical predictions of HFCs with published literature data. A procedure for obtaining the other best-fit result for an ENDOR determined tensor is outlined.

Radiation-Induced Radicals in Glucose-1-phosphate. I. Electron Paramagnetic Resonance and Electron Nuclear Double Resonance Analysis of in situ X-Irradiated Single Crystals at 77 K

H. De Cooman, G. Vanhaelewyn, E. Pauwels, E. Sagstuen, M. Waroquier
Journal of Physical Chemistry B
112 (47), 15045-15053
2008
A1

Abstract 

Electron magnetic resonance analysis of radiation-induced defects in dipotassium glucose-1-phosphate dihydrate single crystals in situ X-irradiated and measured at 77 K shows that at least seven different carbon-centered radical species are trapped. Four of these (R1−R4) can be fully or partly characterized in terms of proton hyperfine coupling tensors. The dominant radical (R2) is identified as a C1-centered species, assumedly formed by a scission of the sugar−phosphate junction and the concerted formation of a carbonyl group at the neighboring C2 carbon. This structure is chemically identical to a radical recently identified in irradiated sucrose single crystals. Radical species R1 and R4 most likely are C3- and C6-centered species, respectively, both formed by a net hydrogen abstraction. R3 is suggested to be chemically similar to but geometrically different from R4. Knowledge of the identity of the sugar radicals present at 77 K provides a first step in elucidating the formation mechanism of the phosphoryl radicals previously detected after X-irradiation at 280 K. In paper II, the chemical identity, precise conformation, and possible formation mechanisms of these radical species are investigated by means of DFT calculations and elementary insight into the radiation chemistry of sugar and sugar derivatives is obtained.

Radiation-Induced Radicals in Glucose-1-phosphate. II. DFT Analysis of Structures and Possible Formation Mechanisms

E. Pauwels, H. De Cooman, G. Vanhaelewyn, E. Sagstuen, F. Callens, M. Waroquier
Journal of Physical Chemistry B
112 (47), 15054-15063
2008
A1

Abstract 

Four radiation-induced carbon-centered radicals in dipotassium glucose-1-phosphate dihydrate single crystals are examined with DFT methods, consistently relying on a periodic computational scheme. Starting from a set of plausible radical models, EPR hyperfine coupling tensors are calculated for optimized structures and compared with data obtained from EPR/ENDOR measurements, which are described in part I of this work. In this way, an independent structural identification is made of all the radicals that were observed in the experiments (R1−R4) and tentative reaction schemes are proposed. Also, the first strong evidence for conformational freedom in sugar radicals is established: two species are found to have the same chemical composition but different conformations and consequently different hyperfine coupling tensors. Analysis of the calculated energies for all model compounds suggests that the radiation chemistry of sugars, in general, is kinetically and not necessarily thermodynamically controlled.

An Intrinsic Radical Stability Scale from the Perspective of Bond Dissociation Enthalpies: A Companion to Radical Electrophilicities

F. De Vleeschouwer, V. Van Speybroeck, M. Waroquier, P. Geerlings, F. De Proft
Journal of Organic Chemistry
73 (22) 9109-9120
2008
A1

Abstract 

Bond dissociation enthalpies (BDEs) of a large series of molecules of the type A−B, where a series of radicals A ranging from strongly electrophilic to strongly nucleophilic are coupled with a series of 8 radicals (CH2OH, CH3, NF2, H, OCH3, OH, SH, and F) also ranging from electrophilic to nucleophilic, are computed and analyzed using chemical concepts emerging from density functional theory, more specifically the electrophilicities of the individual radical fragments A and B. It is shown that, when introducing the concept of relative radical electrophilicity, an (approximately) intrinsic radical stability scale can be developed, which is in good agreement with previously proposed stability scales. For 47 radicals, the intrinsic stability was estimated from computed BDEs of their combinations with the strongly nucleophilic hydroxymethyl radical, the neutral hydrogen atom, and the strongly electrophilic fluorine atom. Finally, the introduction of an extra term containing enhanced Pauling electronegativities in the model improves the agreement between the computed BDEs and the ones estimated from the model, resulting in a mean absolute deviation of 16.4 kJ mol−1. This final model was also tested against 82 experimental values. In this case, a mean absolute deviation of 15.3 kJ mol−1 was found. The obtained sequences for the radical stabilities are rationalized using computed spin densities for the radical systems.

Pages

Subscribe to RSS - M. Waroquier