Modeling the Effects of O-sulfonation on the CID of Serine

Author(s):  
Kenneth Lucas ◽  
George Barnes

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of \textit{O}-sulfonation on the collision induced dissociation for serine. Towards this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m\z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m\z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The \mz 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also seen that the RM1 semi empirical method was not able to obtain all of the major peaks seen in experiment for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.

2020 ◽  
Author(s):  
Kenneth Lucas ◽  
George Barnes

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of \textit{O}-sulfonation on the collision induced dissociation for serine. Towards this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m\z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m\z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The \mz 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also seen that the RM1 semi empirical method was not able to obtain all of the major peaks seen in experiment for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.


2020 ◽  
Author(s):  
Kenneth Lucas ◽  
George Barnes

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of \textit{O}-sulfonation on the collision induced dissociation for serine. Towards this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m\z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m\z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The \mz 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also seen that the RM1 semi empirical method was not able to obtain all of the major peaks seen in experiment for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.


2020 ◽  
Vol 633 ◽  
pp. A103
Author(s):  
Tao Chen ◽  
Yi Luo ◽  
Aigen Li

Aims. Due to the limitations of current computational technology, the fragmentation and isomerization products of vibrationally-excited polycyclic aromatic hydrocarbon (PAH) molecules and their derivatives have been poorly studied. In this work, we investigate the intermediate products of PAHs and their derivatives as well as the gas-phase reactions relevant to the interstellar medium, with coronene as a case study. Methods. Based on the semi-empirical method of PM3 as implemented in the CP2K program, molecular dynamics simulations were performed to model the major processes (e.g., vibrations, fragmentations, and isomerizations) of coronene and its derivatives (e.g., methylated coronene, hydrogenated coronene, dehydrogenated coronene, nitrogen-substituted coronene, and oxygen-substituted coronene) at temperatures of 3000 K and 4000 K. Results. We find that the anharmonic effects are crucial for the simulation of vibrational excitation. For the molecules studied here, H2, CO, HCN, and CH2 are the major fragments. Following the dissociation of these small units, most of the molecules could maintain their ring structures, but a few molecules would completely break into carbon chains. The transformation from a hexagon to a pentagon or a heptagon may occur and the heteroatomic substitutions (e.g., N- or O-substitutions) would facilitate the transformation.


2020 ◽  
Author(s):  
Julian Keupp ◽  
Johannes P. Dürholt ◽  
Rochus Schmid

The prototypical pillared layer MOFs, formed by a square lattice of paddle-<br>wheel units and connected by dinitrogen pillars, can undergo a breathing phase<br>transition by a “wine-rack” type motion of the square lattice. We studied this not<br>yet fully understood behavior using an accurate first principles parameterized force<br>field (MOF-FF) for larger nanocrystallites on the example of Zn 2 (bdc) 2 (dabco) [bdc:<br>benzenedicarboxylate, dabco: (1,4-diazabicyclo[2.2.2]octane)] and found clear indi-<br>cations for an interface between a closed and an open pore phase traveling through<br>the system during the phase transformation [Adv. Theory Simul. 2019, 2, 11]. In<br>conventional simulations in small supercells this mechanism is prevented by periodic<br>boundary conditions (PBC), enforcing a synchronous transformation of the entire<br>crystal. Here, we extend this investigation to pillared layer MOFs with flexible<br>side-chains, attached to the linker. Such functionalized (fu-)MOFs are experimen-<br>tally known to have different properties with the side-chains acting as fixed guest<br>molecules. First, in order to extend the parameterization for such flexible groups,<br>1a new parametrization strategy for MOF-FF had to be developed, using a multi-<br>structure force based fit method. The resulting parametrization for a library of<br>fu-MOFs is then validated with respect to a set of reference systems and shows very<br>good accuracy. In the second step, a series of fu-MOFs with increasing side-chain<br>length is studied with respect to the influence of the side-chains on the breathing<br>behavior. For small supercells in PBC a systematic trend of the closed pore volume<br>with the chain length is observed. However, for a nanocrystallite model a distinct<br>interface between a closed and an open pore phase is visible only for the short chain<br>length, whereas for longer chains the interface broadens and a nearly concerted trans-<br>formation is observed. Only by molecular dynamics simulations using accurate force<br>fields such complex phenomena can be studied on a molecular level.


2019 ◽  
Author(s):  
Chem Int

The full conformational space of N-formyl-L-alanine-amide was explored by the semi-empirical method AM1 coupled to the Multi Niche Crowding (MNC) genetic algorithm implemented in a package of programs developed in our laboratory. The structural and energy analysis of the resulting conformational space E(,ψ) exhibits 5 regions or minima ɣL, ɣD, ɛL, D and αD. The technique provides better detection of local and global minima within a reasonable time.


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


Sign in / Sign up

Export Citation Format

Share Document