Modeling the cluster formation during infrared and ultraviolet matrix-assisted laser desorption ionization of oligonucleotides in succinic acid matrix with molecular mechanics

A large succinic acid (HOOC(CH2)(2)COOH) matrix containing 7 x 7 x 7 unit cells with guest oligonucleotide AGCAGCT was modeled with molecular dynamics simulation for infrared matrix-assisted laser desorption ionization. We simulated the laser heating of the succinic acid (this data is still missing from the literature) with lambda = 2.94 mum infrared laser pulses and compared it to ultraviolet excitation. We did this in order to elucidate the cluster formation of succinic acid in the gas phase in itself and around the analyte. At this wavelength, the laser energy is coupled into the matrix through the OH vibrations of the carboxyl groups. The most pronounced difference we observed at 1,500 K simulation is that infrared heating generates about 10-15 more succinic acid molecules bound to the analyte in noncovalent complex form than the ultraviolet mode, which generates about 2 molecules thus bound. We report energy redistribution within the matrix between the host and guest species as well as other dynamical properties. The parameter and topology data for succinic acid that we used are reported and ready for use in CHARMM computer code environment for simulation.