skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Collisional dynamics simulations revealing fragmentation properties of Zn( ii )-bound poly-peptide
Chemical dynamics simulations are performed to study the collision induced gas phase unimolecular fragmentation of a model peptide with the sequence acetyl-His 1 -Cys 2 -Gly 3 -Pro 4 -Tyr 5 -His 6 -Cys 7 (analogue methanobactin peptide-5, amb 5 ) and in particular to explore the role of zinc binding in reactivity. Fragmentation pathways, their mechanisms, and collision energy transfer are discussed. The probability distributions of the pathways are compared with the results of the experimental IM-MS, MS/MS spectrum and previous thermal simulations. Collisional activation gives both statistical and non-statistical fragmentation pathways with non-statistical shattering mechanisms accounting for a relevant percentage of reactive trajectories, becoming dominant at higher energies. The tetra-coordination of zinc changes qualitative and quantitative fragmentation, in particular the shattering. The collision energy threshold for the shattering mechanism was found to be 118.9 kcal mol −1 which is substantially higher than the statistical Arrhenius activation barrier of 35.8 kcal mol −1 identified previously during thermal simulations. This difference can be attributed to the tetra-coordinated zinc complex that hinders the availability of the sidechains to undergo direct collision with the Ar projectile.  more » « less
Award ID(s):
1764436
PAR ID:
10248603
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
26
ISSN:
1463-9076
Page Range / eLocation ID:
14551 to 14559
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT This study focuses on investigating the conformational structure and zinc(II) affinity of a zinc finger‐like motif (ZFM) peptide with the sequence acetyl‐His1‐Cys2‐Gly3‐Pro4‐Gly5‐His6‐Cys7, where bold highlights the potential zinc(II) binding sites. Zinc fingers are crucial protein motifs known for their high specificity and affinity for zinc ions. The ZFM peptide's sequence contains the 2His‐2Cys zinc‐binding sites similar to those in natural zinc finger proteins but without the hydrophobic core, making it a valuable model for studying zinc(II)–peptide interactions. Previous research on related peptides showed that collision cross sections and B3LYP modeling predicted that the His‐2Cys‐carboxyl terminus coordination of zinc(II) was more stable than the 2His‐2Cys. Employing a comprehensive approach integrating ion mobility–mass spectrometry and theoretical modeling techniques, various zinc(II) binding modes of the ZFM have been thoroughly compared to ascertain their influence on the competitive threshold collision‐induced dissociation method for measuring the relative gas‐phase Zn(II) affinity of the ZFM peptide. The measured Zn(II) affinity of ZFM is greater than those measured recently for two peptides with similar primary structures, acetyl‐His1‐Cys2‐Gly3‐Pro4‐Gly5‐Gly6‐Cys7and acetyl‐Asp1‐His2‐Gly3‐Pro4‐Gly5‐Gly6‐Cys7, indicating the preference for the His1‐Cys2‐His6‐Cys7side groups for coordinating zinc(II) over the His‐2Cys‐carboxyl terminus or Asp‐His‐Cys‐carboxyl terminus in these related heptapeptides. 
    more » « less
  2. Abstract The analog methanobactin (amb) peptide with the sequence ac‐His1‐Cys2‐Gly3‐Pro4‐Tyr5‐His6‐Cys7(amb5A) will bind the metal ions of zinc, nickel, and copper. To further understand how amb5Abinds these metals, we have undertaken a series of studies of structurally related heptapeptides where one or two of the potential His or Cys binding sites have been replaced by Gly, or the C‐terminus has been blocked by amidation. The studies were designed to compare how these metals bind to these sequences in different pH solutions of pH 4.2 to 10 and utilized native electrospray ionization (ESI) with ion mobility‐mass spectrometry (IM‐MS) which allows for the quantitative analysis of the charged species produced during the reactions. The native ESI conditions were chosen to conserve as much of the solution‐phase behavior of the amb peptides as possible and an analysis of how the IM‐MS results compare with the expected solution‐phase behavior is discussed. The oligopeptides studied here have applications for tag‐based protein purification methods, as therapeutics for diseases caused by elevated metal ion levels or as inhibitors for metal‐protein enzymes such as matrix metalloproteinases. 
    more » « less
  3. With dual goals of efficient and accurate modeling of solvation thermodynamics in molten salt liquids, we employ ab initio molecular dynamics (AIMD) simulations, deep neural network interatomic potentials (NNIP), and quasichemical theory (QCT) to calculate the excess chemical potentials for the solute ions Na + and Cl − in the molten NaCl liquid. NNIP-based molecular dynamics simulations accelerate the calculations by 3 orders of magnitude and reduce the uncertainty to 1 kcal mol −1 . Using the Density Functional Theory (DFT) level of theory, the predicted excess chemical potential for the solute ion pair is −178.5 ± 1.1 kcal mol −1 . A quantum correction of 13.7 ± 1.9 kcal mol −1 is estimated via higher-level quantum chemistry calculations, leading to a final predicted ion pair excess chemical potential of −164.8 ± 2.2 kcal mol −1 . The result is in good agreement with a value of −163.5 kcal mol −1 obtained from thermo-chemical tables. This study validates the application of QCT and NNIP simulations to the molten salt liquids, allowing for significant insights into the solvation thermodynamics crucial for numerous molten salt applications. 
    more » « less
  4. Laskin, J; Ouyang, Z (Ed.)
    Chirality effects on the intrinsic gas-phase acidity of oligopeptides have been studied using a pair of stereoisomeric tripeptides consisting of a D/L-cysteine (C) and two residues of alanine (A): CAA and dCAA, where the C-terminus is amidated. Mass spectrometry measurements through bracketing via collision-induced dissociation clearly show that CAA is a stronger gas-phase acid than dCAA. Quantitative values of the acidity were determined using the extended Cooks kinetic method. The resulting deprotonation enthalpy (∆acidH) for CAA is 326.2 kcal/mol (1364.7 kJ/mol) and for dCAA it is 326.8 kcal/mol (1367.6 kJ/mol). The corresponding gas-phase acidity (∆acidG) for CAA is 321.3 kcal/mol (1344.2 kJ/mol) and for dCAA it is 322.0 kcal/mol (1347.3 kJ/mol). Changing the N-terminal cysteine from the L-form to the D-form reduces the gas-phase acidity by about 0.6 kcal/mol (2.5 kJ/mol). Extensive conformational searches followed by quantum chemical calculations at the ωB97X-D/6-311+G(d,p) level of theory yielded a set of lowest energy conformations for each peptide species. Theoretical gas-phase acidities calculated using the Boltzmann averaged conformational contributions are in good agreement with the experimental data. The shift in the acidity is likely due to the conformational effect induced by D-cysteine, which increases the stability of the neutral dCAA, and hence reduces its acidity. A chirality change on a single amino acid can have a noticeable effect on the biochemical properties of peptides and proteins. 
    more » « less
  5. ABSTRACT Access to benzofuran‐2(3H)‐one derivatives from readily available substrates under mild conditions is crucial in the pharmaceutical and plastics industries. We identified (Z)‐3‐(2‐phenylhydrazineylidene)benzofuran‐2(3H)‐one (P) during the recrystallization of (E)‐2‐(2,2‐dichloro‐1‐(phenyldiazenyl)vinyl)phenol using a 96% ethanol solution. The mechanism of the unexpected substrate conversion leading toPis investigated using density functional calculations. The computations revealed that ethanol is required to initiate the reaction viaTS1E, which involves a concerted deprotonation of ethanol by the basic diaza group of the substrate and an ethoxy group attacking the electrophilic center (Cl2C), with an energy barrier of 28.3 kcal/mol. The resulting intermediate (I1E) is calculated to be unstable and can yield a cyclic chloroacetal adduct with a lower energy barrier of 2.2 kcal/mol via the ring‐closure transition state (TS2E). In the absence of water, the next steps are impossible because water is required to cleave the ether bond, yieldingP. A small amount of water (4% of the recrystallization solvent) can promote further transformation ofI2Evia the transition statesTS3E(∆G = 11.1 kcal/mol) andTS4E(∆G = 10.5 kcal/mol). A comparison of the ethanol/water‐ and only water‐promoted free energy profiles shows that the presence of ethanol is crucial for lowering the energy barriers (by about 5 kcal/mol) for the initial two steps leading to the cyclic chloroacetal (I2E), whereas water is then required to initiate product formation. 
    more » « less