Direct Dynamics Simulations of Fragmentation of a Zn(II)-2Cys-2His Oligopeptide. Comparison with Mass Spectrometry Collision-Induced Dissociation

CH vs. HC-Promiscuous Metal Sponges in Antimicrobial Peptides and Metallophores

Histidine and cysteine residues, with their imidazole and thiol moieties that deprotonate at approximately physiological pH values, are primary binding sites for Zn(II), Ni(II) and Fe(II) ions and are thus ubiquitous both in peptidic metallophores and in antimicrobial peptides that may use nutritional immunity as a way to limit pathogenicity during infection. We focus on metal complex solution equilibria of model sequences encompassing Cys-His and His-Cys motifs, showing that the position of histidine and cysteine residues in the sequence has a crucial impact on its coordination properties. CH and HC motifs occur as many as 411 times in the antimicrobial peptide database, while similar CC and HH regions are found 348 and 94 times, respectively. Complex stabilities increase in the series Fe(II) < Ni(II) < Zn(II), with Zn(II) complexes dominating at physiological pH, and Ni(II) ones-above pH 9. The stabilities of Zn(II) complexes with Ac-ACHA-NH₂ and Ac-AHCA-NH₂ are comparable, and a similar tendency is observed for Fe(II), while in the case of Ni(II), the order of Cys and His does matter-complexes in which the metal is anchored on the third Cys (Ac-AHCA-NH₂) are thermodynamically stronger than those where Cys is in position two (Ac-ACHA-NH₂) at basic pH, at which point amides start to take part in the binding. Cysteine residues are much better Zn(II)-anchoring sites than histidines; Zn(II) clearly prefers the Cys-Cys type of ligands to Cys-His and His-Cys ones. In the case of His- and Cys-containing peptides, non-binding residues may have an impact on the stability of Ni(II) complexes, most likely protecting the central Ni(II) atom from interacting with solvent molecules.

Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations

In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.

Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories

In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.

Unimolecular Fragmentation Properties of Thermometer Ions from Chemical Dynamics Simulations

Thermometer ions are widely used to calibrate the internal energy of the ions produced by electrospray ionization in mass spectrometry. Typically, benzylpyridinium ions with different substituents are used. More recently, benzhydrylpyridinium ions were proposed for their lower bond dissociation energies. Direct dynamics simulations using M06-2X/6-31G(d), DFTB, and PM6-D3 are performed to characterize the activation energies of two representative systems: para-methylbenzylpyridinium ion (p-Me-BnPy⁺) and methyl,methylbenzhydrylpyridinium ion (Me,Me-BhPy⁺). Simulation results are used to calculate rate constants for the two systems. These rate constants and their uncertainties are used to find the Arrhenius activation energies and RRK fitted threshold energies which give reasonable agreement with calculated bond dissociation energies at the same level of theory. There is only one fragmentation mechanism observed for both systems, which involves C-N bond dissociation via a loose transition state, to generate either benzylium or benzhydrylium ion and a neutral pyridine molecule. For p-Me-BnPy⁺ using DFTB and PM6-D3 the formation of tropylium ion, from rearrangement of benzylium ion, was observed but only at higher excitation energies and for longer simulation times. These observations suggest that there is no competition between reaction pathways that could affect the reliability of internal energy calibrations. Finally, we suggest using DFTB with a modified-Arrhenius model in future studies.

Collisional dynamics simulations revealing fragmentation properties of Zn(ii)-bound poly-peptide

Collisional simulations show how peptide fragmentation is modified by the presence of Zn(ii).