Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

Biomolecular dynamics with machine-learned quantum-mechanical force fields trained on diverse chemical fragments

The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.

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.

Impact of anion polarizability on ion pairing in microhydrated salt clusters

Despite longstanding interest in the mechanism of salt dissolution in aqueous media, a molecular level understanding remains incomplete. Here, cryogenic ion trap vibrational action spectroscopy is combined with electronic structure calculations to track salt hydration in a gas phase model system one water molecule at a time. The infrared photodissociation spectra of microhydrated lithium dihalide anions [LiXX'(H2O) n ]- (XX' = I2, ClI and Cl2; n = 1-3) in the OH stretching region (3800-2800 cm-1) provide a detailed picture of how anion polarizability influences the competition among ion-ion, ion-water and water-water interactions. While exclusively contact ion pairs are observed for n = 1, the formation of solvent-shared ion pairs, identified by markedly red-shifted OH stretching bands (<3200 cm-1), originating from the bridging water molecules, is favored already for n = 2. For n = 3, Li+ reaches its maximum coordination number of four only in [LiI2(H2O)3]-, in accordance with the hard and soft Lewis acid and base principle. Water-water hydrogen bond formation leads to a different solvent-shared ion pair motif in [LiI2(H2O)3]- and network formation even restabilizes the contact ion pair motif in [LiCl2(H2O)3]-. Structural assignments are exclusively possible after the consideration of anharmonic effects. Molecular dynamics simulations confirm that the significance of large amplitude motion (of the water molecules) increases with increasing anion polarizability and that needs to be considered already at cryogenic temperatures.

Mechanistic Investigations on N-Cα Bond Cleavages in Dibasic Peptides Containing Internal Lys and Arg Residues

The model nonapeptide AAARAAKAG* (* indicates amide) is used to explore N-C^α bond fragmentation under CID-MS conditions. Neighboring group participation and the effect of positioning of Lys and Arg residues on N-C^α bond cleavage is established using a library of synthetic peptide analogues. The importance of the Arg residue at position 4 and the i to i+3 spacing between Arg and Lys residues in determining the formation of the N-C^α bond cleaved product ions (c_n) is demonstrated by a comparative MS study of positional variants in analogue peptides. The effect of shortening of the Lys side chain has been established using ornithine (Orn) and diaminobutyric acid (Dab) analogues. The involvement of the Lys residue in mediating the N-C^α bond cleavage is further probed using N^ε-dimethyl and isotopically labeled ¹⁵N^α, ε lysine residues. MSⁿ experiments reveal that the c₆ ion originates from a doubly charged dehydrated b₈ ion [b₈-18]²⁺. The mechanism of this unusual fragmentation process has been probed by using position 8 analogues (Gly, Ala, and Aib). A plausible mechanism is proposed for the origin of the c₆ ion, which involves C-terminus lactam formation followed by transannular cyclization and dehydration. The results presented in this study highlight the role of reactive side chain functionalities in promoting noncanonical fragmentation pathways.

Conformational analysis of tannic acid: Environment effects in electronic and reactivity properties

Polyphenols are natural molecules of crucial importance in many applications, of which tannic acid (TA) is one of the most abundant and established. Most high-value applications require precise control of TA interactions with the system of interest. However, the molecular structure of TA is still not comprehended at the atomic level, of which all electronic and reactivity properties depend. Here, we combine an enhanced sampling global optimization method with density functional theory (DFT)-based calculations to explore the conformational space of TA assisted by unsupervised machine learning visualization and then investigate its lowest energy conformers. We study the external environment's effect on the TA structure and properties. We find that vacuum favors compact structures by stabilizing peripheral atoms' weak interactions, while in water, the molecule adopts more open conformations. The frontier molecular orbitals of the conformers with the lowest harmonic vibrational free energy have a HOMO-LUMO energy gap of 2.21 (3.27) eV, increasing to 2.82 (3.88) eV in water, at the DFT generalized gradient approximation (and hybrid) level of theory. Structural differences also change the distribution of potential reactive sites. We establish the fundamental importance of accurate structural consideration in determining TA and related polyphenol interactions in relevant technological applications.

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

PSIQUE: Protein Secondary Structure Identification on the Basis of Quaternions and Electronic Structure Calculations

The secondary structure is important in protein structure analysis, classification, and modeling. We have developed a novel method for secondary structure assignment, termed PSIQUE, based on the potential energy surface (PES) of polyalanine obtained using an infinitely long chain model and density functional theory calculations. First, uniform protein segments are determined in terms of a difference of quaternions between neighboring amino acids along the protein backbone. Then, the identification of the secondary structure motifs is carried out based on the minima found in the PES. PSIQUE shows good agreement with other secondary structure assignment methods. However, it provides better discrimination of subtle secondary structures (e.g., helix types) and termini and produces more uniform segments while also accounting for local distortions. Overall, PSIQUE provides a precise and reliable assignment of secondary structures, so it should be helpful for the detailed characterization of the protein structure.

Frenkel-Holstein Hamiltonian applied to absorption spectra of quaterthiophene-based 2D hybrid organic-inorganic perovskites

For the prototypical two-dimensional hybrid organic-inorganic perovskites (2D HOIPs) (AE4T)PbX₄ (X = Cl, Br, and I), we demonstrate that the Frenkel-Holstein Hamiltonian (FHH) can be applied to describe the absorption spectrum arising from the organic component. We first model the spectra using only the four nearest neighbor couplings between translationally inequivalent molecules in the organic herringbone lattice as fitting parameters in the FHH. We next use linear-response time-dependent density functional theory (LR-TDDFT) to calculate molecular transition densities, from which extended excitonic couplings are evaluated based on the atomic positions within the 2D HOIPs. We find that both approaches reproduce the experimentally observed spectra, including changes in their shape and peak positions. The spectral changes are correlated with a decrease in excitonic coupling from X = Cl to X = I. Importantly, the LR-TDDFT-based approach with extended excitonic couplings not only gives better agreement with the experimental absorption line shape than the approach using a restricted set of fitted parameters but also allows us to relate the changes in excitonic coupling to the underlying geometry. We accordingly find that the decrease in excitonic coupling from X = Cl to Br to I is due to an increase in molecular separation, which in turn can be related to the increasing Pb-X bond length from Cl to I. Our research opens up a potential pathway to predicting optoelectronic properties of new 2D HOIPs from ab initio calculations and to gain insight into structural relations from 2D HOIP absorption spectra.

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.

Quantum mechanics of proteins in explicit water: The role of plasmon-like solute-solvent interactions

Quantum-mechanical van der Waals dispersion interactions play an essential role in intraprotein and protein-water interactions-the two main factors affecting the structure and dynamics of proteins in water. Typically, these interactions are only treated phenomenologically, via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach of density-functional tight-binding combined with the many-body dispersion formalism and demonstrate the relevance of many-body van der Waals forces both to protein energetics and to protein-water interactions. In contrast to commonly used pairwise approaches, many-body effects substantially decrease the relative stability of native states in the absence of water. Upon solvation, the protein-water dispersion interaction counteracts this effect and stabilizes native conformations and transition states. These observations arise from the highly delocalized and collective character of the interactions, suggesting a remarkable persistence of electron correlation through aqueous environments and providing the basis for long-range interaction mechanisms in biomolecular systems.

Tailored photocleavable peptides: fragmentation and neutralization pathways in high vacuum

Photocleavable tags (PCTs) have the potential for excellent spatio-temporal control over the release of subunits of complex molecules. Here, we show that electrosprayed oligopeptides, functionalized by a tailored ortho-nitroarylether can undergo site-specific photo-activated cleavage under UV irradiation (266 nm) in high vacuum. The comparison of UV photodissociation (UVPD) and collision-induced dissociation (CID) points to the thermal nature of the cleavage mechanism, a picture corroborated by the temperature dependence of the process. Two competing photodissociation pathways can be identified. In one case a phenolate anion is separated from a neutral zwitterion. In the other case a neutral phenol derivative leaves a negatively charged peptide behind. To understand the factors favoring one channel over the other, we investigate the influence of the peptide length, the nature of the phenolic group and the position of the nitro-group (ortho vs. para). The observed gas phase cleavage of a para-nitro benzylic ether markedly differs from the established behavior in solution.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

Mapping and classifying molecules from a high-throughput structural database

High-throughput computational materials design promises to greatly accelerate the process of discovering new materials and compounds, and of optimizing their properties. The large databases of structures and properties that result from computational searches, as well as the agglomeration of data of heterogeneous provenance leads to considerable challenges when it comes to navigating the database, representing its structure at a glance, understanding structure-property relations, eliminating duplicates and identifying inconsistencies. Here we present a case study, based on a data set of conformers of amino acids and dipeptides, of how machine-learning techniques can help addressing these issues. We will exploit a recently-developed strategy to define a metric between structures, and use it as the basis of both clustering and dimensionality reduction techniques-showing how these can help reveal structure-property relations, identify outliers and inconsistent structures, and rationalise how perturbations (e.g. binding of ions to the molecule) affect the stability of different conformers.

Trends for isolated amino acids and dipeptides: Conformation, divalent ion binding, and remarkable similarity of binding to calcium and lead

We derive structural and binding energy trends for twenty amino acids, their dipeptides, and their interactions with the divalent cations Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+. The underlying data set consists of more than 45,000 first-principles predicted conformers with relative energies up to ~4 eV (~400 kJ/mol). We show that only very few distinct backbone structures of isolated amino acids and their dipeptides emerge as lowest-energy conformers. The isolated amino acids predominantly adopt structures that involve an acidic proton shared between the carboxy and amino function. Dipeptides adopt one of two intramolecular-hydrogen bonded conformations C5 or . Upon complexation with a divalent cation, the accessible conformational space shrinks and intramolecular hydrogen bonding is prevented due to strong electrostatic interaction of backbone and side chain functional groups with cations. Clear correlations emerge from the binding energies of the six divalent ions with amino acids and dipeptides. Cd2+ and Hg2+ show the largest binding energies-a potential correlation with their known high acute toxicities. Ca2+ and Pb2+ reveal almost identical binding energies across the entire series of amino acids and dipeptides. This observation validates past indications that ion-mimicry of calcium and lead should play an important role in a toxicological context.

Nuclear Quantum Effects in H(+) and OH(-) Diffusion along Confined Water Wires

The diffusion of protons and hydroxide ions along water wires provides an efficient mechanism for charge transport that is exploited by biological membrane channels and shows promise for technological applications such as fuel cells. However, what is lacking for a better control and design of these systems is a thorough theoretical understanding of the diffusion process at the atomic scale. Here we focus on two aspects of this process that are often disregarded because of their high computational cost: the use of first-principles potential energy surfaces and the treatment of the nuclei as quantum particles. We consider proton and hydroxide ions in finite water wires using density functional theory augmented with an apolar cylindrical confining potential. We employ machine learning techniques to identify the charged species, thus obtaining an agnostic definition that takes explicitly into account the delocalization of the charge in the Grotthus-like mechanism. We include nuclear quantum effects (NQEs) through the thermostated ring polymer molecular dynamics method and model finite system size effects by considering Langevin dynamics on the potential of mean force of the charged species, allowing us to extract the same "universal" diffusion coefficient from simulations with different wire sizes. In the classical case, diffusion coefficients depend significantly on the potential energy surface, in particular on how dispersion forces modulate water-water distances. NQEs, however, make the diffusion less sensitive to the underlying potential and geometry of the wire.