Structures and IR Spectra of the Gramicidin S Peptide: Pushing the Quest for Low-Energy Conformations

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Extended Lagrangian Born-Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

Diserinol Isophthalamide: A Novel Reagent for Complexation with Biomolecular Anions in Electrospray Ionization Mass Spectrometry

Transferring biomolecules from solution to vacuum facilitates a detailed analysis of molecular structure and dynamics by isolating molecules of interest from a complex environment. However, inherent in the ion desolvation process is the loss of solvent hydrogen bonding partners, which are critical for the stability of a condensed-phase structure. Thus, transfer of ions to vacuum can favor structural rearrangement, especially near solvent-accessible charge sites, which tend to adopt intramolecular hydrogen bonding motifs in the absence of solvent. Complexation of monoalkylammonium moieties (e.g., lysine side chains) with crown ethers such as 18-crown-6 can disfavor structural rearrangement of protonated sites, but no equivalent ligand has been investigated for deprotonated groups. Herein we describe diserinol isophthalamide (DIP), a novel reagent for the gas-phase complexation of anionic moieties within biomolecules. Complexation is observed to the C-terminus or side chains of the small model peptides GD, GE, GG, DF-OMe, VYV, YGGFL, and EYMPME in electrospray ionization mass spectrometry (ESI-MS) studies. In addition, complexation is observed with the phosphate and carboxylate moieities of phosphoserine and phosphotyrosine. DIP performs favorably in comparison to an existing anion recognition reagent, 1,1'-(1,2-phenylene)bis(3-phenylurea), that exhibits moderate carboxylate binding in organic solvent. This improved performance in ESI-MS experiments is attributed to reduced steric constraints to complexation with carboxylate groups of larger molecules. Overall, diserinol isophthalamide is an effective complexation reagent that can be applied in future work to study retention of solution-phase structure, investigate intrinsic molecular properties, and examine solvation effects.

Reconstructing the infrared spectrum of a peptide from representative conformers of the full canonical ensemble

Leucine enkephalin (LeuEnk), a biologically active endogenous opioid pentapeptide, has been under intense investigation because it is small enough to allow efficient use of sophisticated computational methods and large enough to provide insights into low-lying minima of its conformational space. Here, we reproduce and interpret experimental infrared (IR) spectra of this model peptide in gas phase using a combination of replica-exchange molecular dynamics simulations, machine learning, and ab initio calculations. In particular, we evaluate the possibility of averaging representative structural contributions to obtain an accurate computed spectrum that accounts for the corresponding canonical ensemble of the real experimental situation. Representative conformers are identified by partitioning the conformational phase space into subensembles of similar conformers. The IR contribution of each representative conformer is calculated from ab initio and weighted according to the population of each cluster. Convergence of the averaged IR signal is rationalized by merging contributions in a hierarchical clustering and the comparison to IR multiple photon dissociation experiments. The improvements achieved by decomposing clusters containing similar conformations into even smaller subensembles is strong evidence that a thorough assessment of the conformational landscape and the associated hydrogen bonding is a prerequisite for deciphering important fingerprints in experimental spectroscopic data.

Manifolds of low energy structures for a magic number of hydrated sulfate: SO42-(H2O)24

Low energy structures of SO₄^2-(H₂O)₂₄ have been obtained using a combination of classical molecular dynamics simulations and refinement of structures and energies by quantum chemical calculations. Extensive exploration of the potential energy surface led to a number of low-energy structures, confirmed by accurate calibration calculations. An overall analysis of this large set was made after devising appropriate structural descriptors such as the numbers of cycles and their combinations. Low energy structures bear common motifs, the most prominent being fused cycles involving alternatively four and six water molecules. The latter adopt specific conformations which ensure the appropriate surface curvature to form a closed cage without dangling O-H bonds and at the same time provide 12-coordination of the sulfate ion. A prominent feature to take into account is isomerism via inversion of hydrogen bond orientations along cycles. This generates large families of ca. 100 isomers for this cluster size, spanning energy windows of 10-30 kJ mol^-1. This relatively ignored isomerism must be taken into account to identify reliably the lowest energy minima. The overall picture is that the magic number cluster SO₄^2-(H₂O)₂₄ does not correspond to formation of a single, remarkable structure, but rather to a manifold of structural families with similar stabilities. Extensive calculations on isomerization mechanisms within a family indicate that large barriers are associated to direct inversion of hydrogen bond networks. Possible implications of these results for magic number clusters of other anions are discussed.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

A Decapeptide Hydrated by Two Waters: Conformers Determined by Theory and Validated by Cold Ion Spectroscopy

The intrinsic structures of biomolecules in the gas phase may not reflect their native solution geometries. Microsolvation of the molecules bridges the two environments, enabling a tracking of molecular structural changes upon hydration at the atomistic level. We employ density functional calculations to compute a large pool of structures and vibrational spectra for a gas-phase complex, in which a doubly protonated decapeptide, gramicidin S, is solvated by two water molecules. Though most vibrations of this large complex are treated in a harmonic approximation, the water molecules and the vibrations of the host ion coupled to them are locally described by a quantum mechanical vibrational self-consistent field theory with second-order perturbation correction (VSCF-PT2). Guided and validated by the available cold ion spectroscopy data, the computational analysis identifies structures of the three experimentally observed conformers of the complex. They, mainly, differ by the hydration sites, of which the one at the Orn side chain is the most important for reshaping the peptide toward its native structure. The study demonstrates the ability of a quantum chemistry approach that intelligently combines the semiempirical and ab initio computations to disentangle a complex interplay of intra- and intermolecular hydrogen bonds in large molecular systems.

Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach

The present work reports the calculation of vibrational infrared (IR) and Raman spectra of large molecular systems employing molecular tailoring approach (MTA). Further, it extends the grafting procedure for the accurate evaluation of IR and Raman spectra of large molecular systems, employing a new methodology termed as Fragments-in-Fragments (FIF), within MTA. Unlike the previous MTA-based studies, the accurate estimation of the requisite molecular properties is achieved without performing any full calculations (FC). The basic idea of the grafting procedure is implemented by invoking the nearly basis-set-independent nature of the MTA-based error vis-à-vis the respective FCs. FIF has been tested out for the estimation of the above molecular properties for three isomers, viz., β-strand, 310- and α-helix of acetyl(alanine)nNH2 (n = 10, 15) polypeptides, three conformers of doubly protonated gramicidin S decapeptide and trpzip2 protein (PDB id: 1LE1), respectively, employing BP86/TZVP, M06/6-311G**, and M05-2X/6-31G** levels of theory. For most of the cases, a maximum difference of 3 cm(-1) is achieved between the grafted-MTA frequencies and the corresponding FC values. Further, a comparison of the BP86/TZVP level IR and Raman spectra of α-helical (alanine)20 and its N-deuterated derivative shows an excellent agreement with the existing experimental spectra. In view of the requirement of only MTA-based calculations and the ability of FIF to work at any level of theory, the current methodology provides a cost-effective solution for obtaining accurate spectra of large molecular systems.

Structural analyses of isolated cyclic tetrapeptides with varying amino acid residues

Structural analyses of isolated cyclic tetrapeptides with varying amino acid residues were performed by applying combined IR/UV spectroscopy in the molecular beam and DFT calculations. The intrinsic structural properties especially with regard to the influence of different amino acid residues are fundamental for optimizing their binding ability.

Mapping gas phase dipeptide motions in the far-infrared and terahertz domain

Vibrational signatures of Ac-Phe-AA-NH₂ dipeptides are recorded and analysed in the far IR/THz spectral domain (100–800 cm⁻¹, 3–24 THz), with the ‘AA’ amino acid chosen within the series ‘AA’ = Gly, Ala, Pro, Cys, Ser, Val. Phe stands for phenylalanine.

Complexes of Ni(ii) and Cu(ii) with small peptides: deciding whether to deprotonate

Infrared multiple photon dissociation (IRMPD) spectroscopy differentiates two binding modes (iminol versus charge solvated) for Ni(ii) bound to model peptides.

Aminophenol isomers unraveled by conformer-specific far-IR action spectroscopy

Far-infrared action spectroscopy of aminophenol in the gas-phase revealed isomer- and conformer-specific vibrational signatures and provided the heights of NH₂ inversion barrier.

Microhydration of Protonated Nα-Acetylhistidine: A Theoretical Approach

Extensive exploration of the potential energy surfaces of protonated Nα-acetylhistidine hydrated by 0-3 molecules of water was performed. The methodology combined hierarchical and genealogical (Darwin family tree) approaches using polarizable AMOEBA force field and M06 functional. It is demonstrated that this mixed approach allows recovering a larger number of conformers than the number recovered by using any one of the two methods alone. Hydration enthalpies of protonated Nα-acetylhistidine and of model compounds have been computed using higher theoretical methods, up to the G4MP2 procedure. Excellent agreement with experiment is observed for successive hydration of methylamonium and imidazolium cations using MP2/6-311++G(2d,2p)//M06/6-311++G(d,p) and G4MP2 methods, thereby validating the theory levels used for hydrated protonated Nα-acetylhistidine. It is found that the first hydration enthalpy of protonated Nα-acetylhistidine is ca. 10 kJ mol(-1) lower than that of imidazolium, a result explained by the local environment of the positively charged imidazolium moiety.

Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field

The Driven Molecular Dynamics approach has been adapted and associated with the AMOEBA polarizable force field to assign and visualize vibrational modes in infrared spectra obtained by molecular dynamics simulations.

Theoretical Sum Frequency Generation Spectroscopy of Peptides

Vibrational sum frequency generation (SFG) has become a very promising technique for the study of proteins at interfaces, and it has been applied to important systems such as anti-microbial peptides, ion channel proteins, and human islet amyloid polypeptide. Moreover, so-called "chiral" SFG techniques, which rely on polarization combinations that generate strong signals primarily for chiral molecules, have proven to be particularly discriminatory of protein secondary structure. In this work, we present a theoretical strategy for calculating protein amide I SFG spectra by combining line-shape theory with molecular dynamics simulations. We then apply this method to three model peptides, demonstrating the existence of a significant chiral SFG signal for peptides with chiral centers, and providing a framework for interpreting the results on the basis of the dependence of the SFG signal on the peptide orientation. We also examine the importance of dynamical and coupling effects. Finally, we suggest a simple method for determining a chromophore's orientation relative to the surface using ratios of experimental heterodyne-detected signals with different polarizations, and test this method using theoretical spectra.