Accurate proton affinity and gas-phase basicity values for molecules important in biocatalysis

Software Infrastructure for Next-Generation QM/MM-ΔMLP Force Fields

We present software infrastructure for the design and testing of new quantum mechanical/molecular mechanical and machine-learning potential (QM/MM-ΔMLP) force fields for a wide range of applications. The software integrates Amber's molecular dynamics simulation capabilities with fast, approximate quantum models in the xtb package and machine-learning potential corrections in DeePMD-kit. The xtb package implements the recently developed density-functional tight-binding QM models with multipolar electrostatics and density-dependent dispersion (GFN2-xTB), and the interface with Amber enables their use in periodic boundary QM/MM simulations with linear-scaling QM/MM particle-mesh Ewald electrostatics. The accuracy of the semiempirical models is enhanced by including machine-learning correction potentials (ΔMLPs) enabled through an interface with the DeePMD-kit software. The goal of this paper is to present and validate the implementation of this software infrastructure in molecular dynamics and free energy simulations. The utility of the new infrastructure is demonstrated in proof-of-concept example applications. The software elements presented here are open source and freely available. Their interface provides a powerful enabling technology for the design of new QM/MM-ΔMLP models for studying a wide range of problems, including biomolecular reactivity and protein-ligand binding.

Unveiling the Unusual Mn(CO)3 Migration in a Manganese Cyclohexenyl Complex by DFT Computations

Homogeneous catalysis involving a transition metal agostic interaction (TM…H…C) is an attractive strategy for C-H bond activation, in which the transition metal agostic intermediates serve as the critical component. To investigate the roles of manganese agostic intermediates in the unusual migration of the Mn(CO)3 fragment in the (exo-phenyl)(η3-cyclohexenyl)manganese tricarbonyl [(Ph)(η3-C6H8)Mn(CO)3] (complex 1) under the protonation of tetrafluoroboric acid-diethyl ether (HBF4.Et2O), a comprehensive density functional theory (DFT) theoretical study was performed. The computational results showed that formation of the [(cyclohex-3-enyl)-η6-benzene]manganese tricarbonyl complex [(C6H9)(η6-Ph)Mn(CO)3+][BF4] (complex 2) was achieved via a series of mono-agostic and di-agostic intermediates. The overall rate-limiting step for this unusual migration of the Mn(CO)3 fragment is the formation of the di-agostic (η2-phenyl)manganese complex 8 (4 → 5 → 8) with a Gibbs barrier of 15.4 kcal mol-1. The agostic intermediates with TM…H…C agostic interactions were well-characterized by geometry parameters, Atoms-In-Molecules (AIM) analyses, and the Natural Adaptive Orbitals (NAdOs). The located pathways in the current study successfully explained the experimental observations, and the findings on the TM…H…C agostic interaction provided a new aspect of the catalytic reaction with the manganese complex.

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

Unraveling redox pathways of the disulfide bond in dimethyl disulfide: Ab initio modeling

CONTEXT

In cellular environments, the reduction of disulfide bonds is pivotal for protein folding and synthesis. However, the intricate enzymatic mechanisms governing this process remain poorly understood. This study addresses this gap by investigating a disulfide bridge reduction reaction, serving as a model for comprehending electron and proton transfer in biological systems. Six potential mechanisms for reducing the dimethyl disulfide (DMDS) bridge through electron and proton capture were explored. Thermodynamic and kinetic analyses elucidated the sequence of proton and electron addition. MD-PMM, a method that combines molecular dynamics simulations and quantum-chemical calculations, was employed to compute the redox potential of the mechanism. This research provides valuable insights into the mechanisms and redox potentials involved in disulfide bridge reduction within proteins, offering an understanding of phenomena that are challenging to explore experimentally.

METHODS

All calculations used the Gaussian 09 software package at the MP2/6-311 + g(d,p) theory level. Visualization of the molecular orbitals and electron densities was conducted using Gaussview6. Molecular dynamics simulations were performed using GROMACS with the CHARMM36 force field. The PyMM program (Python Program for QM/MM Simulations Based on the Perturbed Matrix Method) is used to apply the Perturbed Matrix Method to MD simulations.

Tautomerization of H+KPGG: Entropic Consequences of Strong Hydrogen-Bond Networks in Peptides

Ion mobility spectrometry-mass spectrometry and quantum chemical calculations are used to determine the structures and stabilities of the singly protonated peptide H+KPGG. The two peaks making up the IMS distribution are shown to be tautomers differing by the location of the extra proton on either the lysine side chain or the N-terminus. The lysine-protonated tautomer is strongly preferred entropically while being disfavored in terms of the electronic energy and enthalpy. This relationship is shown, through comparison of all low-lying conformers of both tautomers, to be related to the strong hydrogen-bond network of the N-terminally protonated tautomer. A general relationship is demonstrated wherein stronger cross-peptide hydrogen-bond networks result in entropically disfavored conformers. Further effects of the H+KPGG hydrogen-bond network are probed by computationally examining singly and doubly methylated analogues. These results demonstrate the importance of the entropic consequences of hydrogen bonds to peptide stability as well as techniques for perturbing the hydrogen-bond network and folding preferences of peptides via minimal chemical modification.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4  kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins

Current experiments only provide low-resolution information on gaseous protein ions generated by electrospray ionization (ESI). Molecular dynamics (MD) simulations can yield complementary insights. Unfortunately, conventional MD does not capture the mobile nature of protons in gaseous proteins. Mobile-proton MD (MPMD) overcomes this limitation. Earlier MPMD data at 300 K indicated that protein ions generated by "native" ESI retain solution-like structures with a hydrophobic core and zwitterionic exterior [Bakhtiari, M.; Konermann, L. J. Phys. Chem. B 2019, 123, 1784-1796]. MPMD redistributes protons using electrostatic and proton affinity calculations. The robustness of this approach has never been scrutinized. Here, we close this gap by benchmarking MPMD against density functional theory (DFT) at the B3LYP/6-31G* level, which is well suited for predicting proton affinities. The computational cost of DFT necessitated the use of small peptides. The MPMD energetic ranking of proton configurations was found to be consistent with DFT single-point energies, implying that MPMD can reliably identify favorable protonation sites. Peptide MPMD runs converged to DFT-optimized structures only when applying 300-500 K temperature cycling, which was necessary to prevent trapping in local minima. Temperature cycling MPMD was then applied to gaseous protein ions. Native ubiquitin converted to slightly expanded structures with a zwitterionic core and a nonpolar exterior. Our data suggest that such inside-out protein structures are intrinsically preferred in the gas phase, and that they form in ESI experiments after moderate collisional excitation. This is in contrast to native ESI (with minimal collisional excitation, simulated by MPMD at 300 K), where kinetic trapping promotes the survival of solution-like structures. In summary, this work validates the MPMD approach for simulations on gaseous peptides and proteins.

Oxidation Products from the Neolignan Licarin A by Biomimetic Reactions and Assessment of in vivo Acute Toxicity

Licarin A, a dihydrobenzofuranic neolignan presents in several medicinal plants and seeds of nutmeg, exhibits strong activity against protozoans responsible for Chagas disease and leishmaniasis. From biomimetic reactions by metalloporphyrin and Jacobsen catalysts, seven products were determined: four isomeric products yielded by epoxidation from licarin A, besides a new product yielded by a vicinal diol, a benzylic aldehyde, and an unsaturated aldehyde in the structure of the licarin A. The incubation with rat and human liver microsomes partially reproduced the biomimetic reactions by the production of the same epoxidized product of m/z 343 [M + H]⁺. In vivo acute toxicity assays of licarin A suggested liver toxicity based on biomarker enzymatic changes. However, microscopic analysis of tissues sections did not show any tissue damage as indicative of toxicity after 14 days of exposure. New metabolic pathways of the licarin A were identified after in vitro biomimetic oxidation reaction and in vitro metabolism by rat or human liver microsomes.

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

We report QDπ-v1.0 for modeling the internal energy of drug molecules containing H, C, N, and O atoms. The QDπ model is in the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP) that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model that is corrected to a quantitatively high-level of accuracy through a deep-learning potential (DeepPot-SE). The model has the advantage that it is able to properly treat electrostatic interactions and handle changes in charge/protonation states. The model is trained against reference data computed at the ωB97X/6-31G* level (as in the ANI-1x data set) and compared to several other approximate semiempirical and machine learning potentials (ANI-1x, ANI-2x, DFTB3, MNDO/d, AM1, PM6, GFN1-xTB, and GFN2-xTB). The QDπ model is demonstrated to be accurate for a wide range of intra- and intermolecular interactions (despite its intended use as an internal energy model) and has shown to perform exceptionally well for relative protonation/deprotonation energies and tautomers. An example application to model reactions involved in RNA strand cleavage catalyzed by protein and nucleic acid enzymes illustrates QDπ has average errors less than 0.5 kcal/mol, whereas the other models compared have errors over an order of magnitude greater. Taken together, this makes QDπ highly attractive as a potential force field model for drug discovery.

Benchmark Ab Initio Determination of the Conformers, Proton Affinities, and Gas-Phase Basicities of Cysteine

A systematic conformational mapping combined with literature data leads to 85 stable neutral cysteine conformers. The implementation of the same mapping process for the protonated counterparts reveals 21 N-(amino-), 64 O-(carbonyl-), and 37 S-(thiol-)protonated cysteine conformers. Their relative energies and harmonic vibrational frequencies are given at the MP2/aug-cc-pVDZ level of theory. Further benchmark ab initio computations are performed for the 10 lowest-lying neutral and protonated amino acid conformers (for each type) such as CCSD(T)-F12a/cc-pVDZ-F12 geometry optimizations (and frequency computations for cysteine) as well as auxiliary correction computations of the basis set effects up to CCSD(T)-F12b/cc-pVQZ-F12, electron correlation effects up to CCSDT(Q), core correlation effects, second-order Douglass-Kroll relativistic effects, and zero-point energy contributions. Boltzmann-averaged 0 (298.15) K proton affinity and [298.15 K gas-phase basicity] values of cysteine are predicted to be 214.96 (216.39) [208.21], 201.83 (203.55) [194.16], and 193.31 (194.74) [186.40] kcal/mol for N-, O-, and S-protonation, respectively, also considering the previously described auxiliary corrections.

Proton Affinity Values of Fentanyl and Fentanyl Analogues Pertinent to Ambient Ionization and Detection

Proton affinity is a major factor in the atmospheric pressure chemical ionization of illicit drugs. The detection of illicit drugs by mass spectrometry and ion mobility spectrometry relies on the analytes having greater proton affinities than background species. Evaluating proton affinities for fentanyl and its analogues is informative for predicting the likelihood of ionization in different environments and for optimizing the compounds' ionization and detection, such as through the addition of dopant chemicals. Herein, density functional theory was used to computationally determine the proton affinity and gas-phase basicity of 15 fentanyl compounds and several relevant molecules as a reference point. The range of proton affinities for the fentanyl compounds was from 1018 to 1078 kJ/mol. Fentanyl compounds with the higher proton affinity values appeared to form a bridge between the oxygen on the amide and the protonated nitrogen on the piperidine ring based on models and calculated bond distances. Experiments with fragmentation of proton-bound clusters using atmospheric flow tube-mass spectrometry (AFT-MS) provided estimates of relative proton affinities and showed proton affinity values of fentanyl compounds >1000 kJ/mol, which were consistent with the computational results. The high proton affinities of fentanyl compounds facilitate their detection by ambient ionization techniques in complex environments. The detection limits of the fentanyl compounds with AFT-MS are in the low femtogram range, which demonstrates the feasibility of trace vapor drug detection.

A New Pt(II) Complex with Anionic s-Triazine Based NNO-Donor Ligand: Synthesis, X-ray Structure, Hirshfeld Analysis and DFT Studies

The reaction of PtCl₂ with s-triazine-type ligand (HTriaz) (1:1) in acetone under heating afforded a new [Pt(Triaz)Cl] complex. Single-crystal X-ray diffraction analysis showed that the ligand (HTriaz) is an NNO tridentate chelate via two N-atoms from the s-triazine and hydrazone moieties and one oxygen from the deprotonated phenolic OH. The coordination environment of the Pt(II) is completed by one Cl⁻¹ ion trans to the Pt-N_(hydrazone). Hirshfeld surface analysis showed that the most dominant interactions are the H···H, H···C and O···H intermolecular contacts. These interactions contributed by 60.9, 11.2 and 8.3% from the whole fingerprint area, respectively. Other minor contributions from the Cl···H, C···N, N···H and C···C contacts were also detected. Among these interactions, the most significant contacts are the O···H, H···C and H···H interactions. The amounts of the electron transfer from the ligand groups to Pt(II) metal center were predicted using NBO calculations. Additionally, the electronic spectra were assigned based on the TD-DFT calculations.

Proton affinity and gas phase basicity of diamandoid molecules: diamantane to C131H116

Calculated proton affinities (PAs) and gas phase basicities (GPBs) are reported for diamantane (C₁₄H₂₀), triamantane (C₁₈H₂₄), 'globular and planar' isomers of tetramantane (C₂₂H₂₈) and pentamantane (C₂₆H₃₂), and for one 'globular' isomer of each of the larger diamondoid molecules: C₅₁H₅₈, C₇₈H₇₂, C₁₀₂H₉₀, and C₁₃₁H₁₁₆. Assuming C_xH_y as the parent diamondoid molecule, we calculated PA and GPB values for a variety of C_xH_y+1⁺ isomers, as well as for the reaction C_xH_y + H⁺ yielding C_xH_y-1⁺ + H₂(g); the latter is slightly favored based on GPB values for diamantane through pentamantane, but less favored compared to certain C_xH_y+1⁺ isomers of C₅₁H₅₈, C₁₀₂H₉₀, and C₁₃₁H₁₁₆. Indeed, the GPB values of C₅₁H₅₈, C₁₀₂H₉₀, and C₁₃₁H₁₁₆ classifiy them as 'superbases'. Calculations that had the initial location of the proton in an interstitial site inside the diamondoid molecule always showed the H having moved to the outside of the diamondoid molecule; for this reason, we focused on testing a variety of initial configurations with the proton placed in an initial position on the surface. Additional protons were added to determine the limiting number that could be, per these calculations, taken up by the diamondoid molecules and the maximum number of protons are shown in parentheses: C₁₄H₂₀(2), C₁₈H₂₄(3), C₂₂H₂₈(3), C₂₆H₃₂(3), C₅₁H₅₈(4). Bader charge distributions obtained for C_xH_y+1⁺ isomers (for diamantane through pentamantane) suggest that the positive charge is essentially completely delocalized over all the H atoms. NMR spectra were calculated for different isomers of C₁₄H₁₉⁺, and compared to the published NMR spectrum for when diamantane was mixed with magic acid and H₂(g) was produced.

Molecular Basis for the Remarkably Different Gas-Phase Behavior of Deprotonated Thyroid Hormones Triiodothyronine (T3) and Reverse Triiodothyronine (rT3): A Clue for Their Discrimination?

Thyroid hormones are biologically active small molecules responsible for growth and development regulation, basal metabolic rate, and lipid and carbohydrate metabolism. Liquid chromatography mass spectrometry (LC–MS) can be used to quantify thyroid hormones blood level with high speed and selectivity, aiming to improve the diagnosis and treatment of the severe pathological conditions in which they are implicated, i.e., hypo- and hyperthyroidism. In this work, the gas-phase behavior of the isomeric thyroid hormones triiodothyronine (T3) and reverse triiodothyronine (rT3) in their deprotonated form was studied at a molecular level using MS-based techniques. Previously reported collision-induced dissociation experiments yielded distinct spectra despite the high structural similarity of the two compounds, suggesting different charge sites to be responsible. Infrared multiple photon dissociation spectroscopy on [T3-H]⁻ and [rT3-H]⁻ was performed, and the results were interpreted using DFT and MP2 calculations, assessing the prevalence of T3 in the carboxylate form and rT3 as a phenolate isomer. The different deprotonation sites of the two isomers were also found to drive their ion-mobility behavior. In fact, [T3-H]⁻ and [rT3-H]⁻ were successfully separated. Drift times were correlated with collisional cross section values of 209 and 215 Ų for [T3-H]⁻ and [rT3-H]⁻, respectively. Calculations suggested the charge site to be the main parameter involved in the different mobilities of the two anions. Finally, bare [T3-H]⁻ and [rT3-H]⁻ were made to react with neutral acetylacetone and trifluoroacetic acid, confirming rT3 to be more acidic than T3 in agreement with the calculated gas-phase acidities of T3 and rT3 equal to 1345 and 1326 kJ mol^–1, respectively.

Substituent Effects on the Basicity of Patriscabrin A and Lettucenin A: Evolution Favors the Aromatic?

Basicities for derivatives of patriscabrin A and lettucenin A were calculated with density functional theory. A significant correlation is observed between the basicity and Hammett σ parameters. Protonation increases the aromatic character of the cyclic moieties of each natural product. The naturally occurring structures are predicted to be the most aromatic.

Synthesis and Halochromic Properties of 1,2,6-Tri- and 1,2,3,6-Tetra-aryl Azulenes

A series of novel 2,6-functionalized azulene molecules Azu1-3 with varied fluorene substituents at the 1- and 3-positions of azulene as well as at the 5'-position of 2-thiophene group were synthesized. Their electronic absorption and emission spectra at neutral and protonated states were examined. It was found that after functionalization with fluorenyl groups, Azu1-3 exhibited absorption maxima at 445, 451 to 468 nm, respectively. In contrast, their corresponding protonated species showed much redshifted absorption maxima at 560, 582 to 643 nm, respectively, mainly due to the extension of conjugation length and the large dipole moment along the C_2v axis of 2,6-substituted azulene molecules. Azu1-3 are non-fluorescent in their neutral forms, but became emissive in their protonated states. Analysis of absorption and emission spectra shows that substitution of the 1- or 3-position of azulene led to decrease in response to trifluoroacetic acid.

Benchmark ab initio proton affinity of glycine

A systematic conformational search reveals three N- (amino) and eight O- (carbonyl) protonated glycine conformers with benchmark equilibrium(adiabatic) relative energies in the 0.00-7.51(0.00-7.37) and 25.91-31.61(24.45-30.28) kcal mol-1 ranges, respectively. Benchmark ab initio structures of the glycine conformers and its protonated species are obtained at the CCSD(T)-F12b/aug-cc-pVTZ level of theory and the relative energy computations consider basis-set effects up to aug-cc-pVQZ with CCSD(T)-F12b, electron correlation up to CCSDT(Q), core correlation corrections, scalar relativistic effects, and zero-point energy contributions. The best predictions for Boltzmann-averaged 0(298.15) K proton affinities and [298.15 K gas-phase basicities] of glycine are 211.00(212.43)[204.75] and 186.38(187.64)[180.21] kcal mol-1 for N- and O-protonation, respectively, in excellent agreement with experiments.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Detection of Pseudomonas aeruginosa quorum sensing molecules at an electrified liquid|liquid micro-interface through facilitated proton transfer

Miniaturization of electrochemical detection methods for point-of-care-devices is ideal for their integration and use within healthcare environments. Simultaneously, the prolific pathogenic bacteria Pseudomonas aeruginosa poses a serious health risk to patients with compromised immune systems. Recognizing these two factors, a proof-of-concept electrochemical method employing a micro-interface between water and oil (w/o) held at the tip of a pulled borosilicate glass capillary is presented. This method targets small molecules produced by P. aeruginosa colonies as signalling factors that control colony growth in a pseudo-multicellular process known as quorum sensing (QS). The QS molecules of interest are 4-hydroxy-2-heptylquinoline (HHQ) and 2-heptyl-3,4-dihydroxyquinoline (PQS, Pseudomonas quinolone signal). Hydrophobic HHQ and PQS molecules, dissolved in the oil phase, were observed electrochemically to facilitate proton transfer across the w/o interface. This interfacial complexation can be exploited as a facile electrochemical detection method for P. aeruginosa and is advantageous as it does not depend on the redox activity of HHQ/PQS. Interestingly, the limit-of-linearity is reached as [H+] ≈ [ligand]. Density functional theory calculations were performed to determine the proton affinities and gas-phase basicities of HHQ/PQS, as well as elucidate the likely site of stepwise protonation within each molecule.

Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids

Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.

Mobile Protons Limit the Stability of Salt Bridges in the Gas Phase: Implications for the Structures of Electrosprayed Protein Ions

Electrosprayed protein ions can retain native-like conformations. The intramolecular contacts that stabilize these compact gas-phase structures remain poorly understood. Recent work has uncovered abundant salt bridges in electrosprayed proteins. Salt bridges are zwitterionic BH⁺/A^- contacts. The low dielectric constant in the vacuum strengthens electrostatic interactions, suggesting that salt bridges could be a key contributor to the retention of compact protein structures. A problem with this assertion is that H⁺ are mobile, such that H⁺ transfer can convert salt bridges into neutral B⁰/HA⁰ contacts. This possible salt bridge annihilation puts into question the role of zwitterionic motifs in the gas phase, and it calls for a detailed analysis of BH⁺/A^- versus B⁰/HA⁰ interactions. Here, we investigate this issue using molecular dynamics (MD) simulations and electrospray experiments. MD data for short model peptides revealed that salt bridges with static H⁺ have dissociation energies around 700 kJ mol^-1. The corresponding B⁰/HA⁰ contacts are 1 order of magnitude weaker. When considering the effects of mobile H⁺, BH⁺/A^- bond energies were found to be between these two extremes, confirming that H⁺ migration can significantly weaken salt bridges. Next, we examined the protein ubiquitin under collision-induced unfolding (CIU) conditions. CIU simulations were conducted using three different MD models: (i) Positive-only runs with static H⁺ did not allow for salt bridge formation and produced highly expanded CIU structures. (ii) Zwitterionic runs with static H⁺ resulted in abundant salt bridges, culminating in much more compact CIU structures. (iii) Mobile H⁺ simulations allowed for the dynamic formation/annihilation of salt bridges, generating CIU structures intermediate between scenarios (i) and (ii). Our results uncover that mobile H⁺ limit the stabilizing effects of salt bridges in the gas phase. Failure to consider the effects of mobile H⁺ in MD simulations will result in unrealistic outcomes under CIU conditions.

For Whom a Puddle Is the Sea? Adsorption of Organic Guests on Hydrated MCM-41 Silica

Thermal and hydration effects on the mobility of compact and branched organic molecules and a bulky pharmaceutical substance loaded in submonolayer amounts onto mesoporous silica have been elucidated using ¹H and ³¹P solid-state NMR. In all cases, the ambient hydration has a stronger effect than an increase in temperature to 370 K for water-free silica. The effect of hydration depends on the guest and ranges from complete solvation to a silica-water-guest sandwich structure to a silica-guest/silica-water pattern. The mobility of the guests under different conditions has been described. The specific structure of the MCM-41 surface allows one to study very slow surface diffusion, a diffusivity of about 10^-15-10^-16 m²/s. The data reported are relevant to any nonfunctionalized silica, while the method used is applicable to any phosphor-containing guest on any host.

Gas-Phase Protonation Thermodynamics of Biological Lipids: Experiment, Theory, and Implications

Phospholipids are important to cellular function and are a vital structural component of plasma and organelle membranes. These membranes isolate the cell from its environment, allow regulation of the internal concentrations of ions and small molecules, and host diverse types of membrane proteins. It remains extremely challenging to identify specific membrane protein-lipid interactions and their relative strengths. Native mass spectrometry, an intrinsically gas-phase method, has recently been demonstrated as a promising tool for identifying endogenous protein-lipid interactions. However, to what extent the identified interactions reflect solution- versus gas-phase binding strengths is not known. Here, the "Extended" Kinetic Method and ab initio computations at three different levels of theory are used to experimentally and theoretically determine intrinsic gas-phase basicities (GB, ΔG for deprotonation of the protonated base) and proton affinities (PA, ΔH for deprotonation of the protonated base) of six lipids representing common phospholipid types. Gas-phase acidities (ΔG and ΔH for deprotonation) of neutral phospholipids are also evaluated computationally and ranked experimentally. Intriguingly, it is found that two of these phospholipids, sphingomyelin and phosphatidylcholine, have the highest GB of any small, monomeric biomolecules measured to date and are more basic than arginine. Phosphatidylethanolamine and phosphatidylserine are found to be similar in GB to basic amino acids lysine and histidine, and phosphatidic acid and phosphatidylglycerol are the least basic of the six lipid types studied, though still more basic than alanine. Kinetic Method experiments and theory show that the gas-phase acidities of these phospholipids are high but less extreme than their GB values, with phosphatidylserine and phosphatidylglycerol being the most acidic. These results indicate that sphingomyelin and phosphatidylcholine lipids can act as charge-reducing agents when dissociated from native membrane protein-lipid complexes in the gas phase and provide a straightforward model to explain the results of several recent native mass spectrometry studies of protein-lipid complexes.

Charge Retention/Charge Depletion in ESI-MS: Theoretical Rationale

Gas phase modification in ESI-MS can significantly alter the charge state distribution of small peptides and proteins. The preceding paper presented a systematic experimental study on this topic using Substance P and proposed a charge retention/charge depletion mechanism, explaining different gas- and liquid-phase modifications [Thinius et al. J. Am. Soc. Mass Spec. 2020, 10.1021/jasms.9b00044]. In this work, we aim to support this rational by theoretical investigations on the proton transfer processes from (multiply) charged analytes toward solvent clusters. As model systems we use small (di)amines as analytes and methanol (MeOH) and acetonitrile (ACN) as gas phase modifiers. The calculations are supported by a set of experiments using (di)amines, to bridge the gap between the present model system and Substance P used in the preceding study. Upon calculation of the thermochemical stability as well as the proton transfer pathways, we find that both ACN and MeOH form stable adduct clusters at the protonation site. MeOH can form large clusters through a chain of H-bridges, eventually lowering the barriers for proton transfer to an extent that charge transfer from the analyte to the MeOH cluster becomes feasible. ACN, however, cannot form H-bridged structures due to its aprotic nature. Hence, the charge is retained at the original protonation site, i.e., the analyte. The investigation confirms the proposed charge retention/charge depletion model. Thus, adding aprotic solvent vapors to the gas phase of an ESI source more likely yields higher charge states than using protic compounds.

Charge Retention/Charge Depletion in ESI-MS: Experimental Evidence

The effects of liquid and gas phase additives (chemical modifiers) on the ion signal distribution for Substance P (SP), recorded with a nanoelectrospray setup, are evaluated. Depletion of the higher charge state of Substance P ([SP+3H]³⁺) is observed with polar protic gas phase modifiers. This is attributed to their ability to form larger hydrogen-bonded clusters, whose proton affinity increases with cluster size. These clusters are able to deprotonate the higher charge state. "Supercharging agents" (SCAs) as well as aprotic polar gas phase modifiers, which promote the retention of the higher charge state of Substance P, do not form such large clusters under the given conditions and are therefore not able to deprotonate Substance P. Both SCAs and aprotic modifiers form clusters with the higher charge state, leading to stabilization of the charge. Whereas supercharging agents have low vapor pressures and are therefore enriched in late-stage electrospray droplets, the gas phase modifiers are volatile organic solvents. Collision induced dissociation experiments revealed that the addition of a modifier significantly delays the droplet evaporation and ion release process. This indicates that the droplet takes up the gas phase modifier to a certain extent (accommodation). Depending on the modifier's properties either charge depletion or retention may eventually be promoted.

Rapid and accurate molecular deprotonation energies from quantum alchemy

We assess the applicability of alchemical perturbation density functional theory (APDFT) for quickly and accurately estimating deprotonation energies. We have considered all possible single and double deprotonations in one hundred small organic molecules drawn at random from QM9 [Ramakrishnan et al., JCTC, 2015]. Numerical evidence is presented for 5160 deprotonated species at both HF/def2-TZVP and CCSD/6-31G* levels of theory. We show that the perturbation expansion formalism of APDFT quickly converges to reliable results: using CCSD electron densities and derivatives, regular Hartree-Fock calculations are outperformed at the second or third order for ranking all possible doubly or singly deprotonated molecules, respectively. CCSD single deprotonation energies are reproduced within 1.4 kcal mol^-1 on average within third order APDFT. We introduce a hybrid approach where the computational cost of APDFT is reduced even further by mixing first order terms at a higher level of theory (CCSD) with higher order terms at a lower level of theory only (HF). We find that this approach reaches 2 kcal mol^-1 accuracy in absolute deprotonation energies compared to CCSD at 2% of the computational cost of third order APDFT.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

Density Functional Studies of Coenzyme NADPH and Its Oxidized Form NADP+ : Structures, UV-Vis Spectra, and the Oxidation Mechanism of NADPH

Density functional theory has been used to study the biologically important coenzyme NADPH and its oxidized form NADP⁺ . It was found that free NADPH prefers a compact structure in gas phase and exists in more extended geometries in aqueous solution. Ultraviolet-visible absorption spectra in aqueous solution were calculated for NADPH with an explicit treatment of 100 surrounding water molecules in combination with the COSMO solvation model for bulk hydration effects. The obtained spectra using the B3LYP hybrid density functional agree quite well with experimental data. The changes of Gibbs free energies ΔG in reactions of NADPH with O₂ observed experimentally in cardiovascular and in chemical systems, that is, NADPH + 2 ³ O₂  → NADP⁺  + 2 O₂ ^-  + H⁺ and NADPH + ¹ O₂  + H⁺  → NADP⁺  + H₂ O₂ , respectively, were calculated. The NADPH oxidation reaction in the cardiovascular system cannot proceed without activation since the obtained ΔG is positive. The reaction of NADPH in the chemical system with singlet oxygen was found to proceed in two ways, each consisting of two steps, that is, NADPH firstly reacts with ¹ O₂ barrierlessly to form NADP⁺ and HO₂ ^- , from which H₂ O₂ is formed in a spontaneous reaction with H⁺ , or ¹ O₂ and H⁺ initially form ¹ HO₂ ⁺ , which further reacts with NADPH to yield NADP⁺ and H₂ O₂ . © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Alkali Metal Cation Affinities of Neutral Maingroup-Element Hydrides across the Periodic Table

We have carried out an extensive quantum chemical exploration of gas-phase alkali metal cation affinities (AMCAs) of archetypal neutral bases across the periodic system using relativistic density functional theory. One objective of this work is to provide an intrinsically consistent set of values of the 298 K AMCAs of all neutral maingroup-element hydrides XH_n of groups 15–18 along the periods 1–6. Our main purpose is to understand these trends in terms of the underlying bonding mechanism using Kohn–Sham molecular orbital theory together with a canonical energy decomposition analysis (EDA). We compare the trends in XH_n AMCAs with the trends in XH_n proton affinities (PAs). We also examine the differences between the trends in AMCAs of the neutral XH_n bases with those in the corresponding anionic XH_n–1^– bases. Furthermore, we analyze how the cation affinity of our neutral Lewis bases changes along the group-1 cations H⁺, Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺.

On the Use of Popular Basis Sets: Impact of the Intramolecular Basis Set Superposition Error

The magnitude of intramolecular basis set superposition error (BSSE) is revealed via computing systematic trends in molecular properties. This type of error is largely neglected in the study of the chemical properties of small molecules and it has historically been analyzed just in the study of large molecules and processes dominated by non-covalent interactions (typically dimerization or molecular complexation and recognition events). In this work we try to provide proof of the broader prevalence of this error, which permeates all types of electronic structure calculations, particularly when employing insufficiently large basis sets.

Alternated Branching Ratios by Anomaly in Collision-Induced Dissociation of Proton-Bound Hoogsteen Base Pairs of 1-Methylcytosine with 1-Methylguanine and 9-Methylguanine

A comparative study on the proton-bound complexes of 1-methylcytosine (1-mC) with 1-methylguanine (1-mG) and 9-methylguanine (9-mG), [1-mC:1-mG:H]⁺ and [1-mC:9-mG:H]⁺, respectively, was carried out using energy-resolved collision-induced dissociation (ER-CID) experiments in combination with quantum chemical calculations. In ER-CID experiments, the measured survival yields indicated an essentially identical stability for the two proton-bound complexes. In comparison with the lowest-energy structures and base-pairing energetics predicted at the B3LYP/6-311+G(2d,2p) theory level, both complexes produced in this study were suggested to be proton-bound Hoogsteen base pairs. Curiously, despite the similarity in structures, binding energetics, and potential energy surfaces predicted by the B3LYP theory, the fragment branching ratios exhibited an intriguing alternation between the two proton-bound Hoogsteen base pairs. The CID of [1-mC:1-mG:H]⁺ produced protonated cytosines, [1-mC:H]⁺, more abundantly than [1-mG:H]⁺, whereas that of [1-mC:9-mG:H]⁺ gave rise to a more pronounced production of protonated guanines, [9-mG:H]⁺. However, using the proton affinities of moieties predicted by the high-accuracy methods, including CBS-QB3 and the Guassian-4 theory, the anomaly known for [Cytosine:Guanine:H]⁺ (J. Am. Soc. Mass Spectrom. 29, 2368-2379 (2018)) successfully accounted for the alternated branching ratios. Thereby, the anomaly, more specifically, the production of proton-transferred fragments of O-protonated cytosines in the CID of proton-bound Hoogsteen base pairs, is indeed real, which is disclosed as the alternated branching ratios in the CID spectra of [1-mC:1-mG:H]⁺ and [1-mC:9-mG:H]⁺ in this study. Graphical Abstract .

Thiaborane clusters with an exoskeletal B-H group

The thiaboranes closo-1-SB11H11 (1a) and 12-I-closo-1-SB11H10 (1b) react with 4-(dimethylamino)pyridine under inert conditions upon the formation of the nido-type thiaboranes 9-B{(4-Me2N)C5NH4}2(H)-7-SB10H11 and 9-B{(4-Me2N)C5NH4}2(H)-5-I-7-SB10H10 containing an exoskeletal B-H group. The same type of B-H moiety is also stabilised by one bipyridine molecule in a chelating fashion. These complexes are unstable in solution, and in air and hydrolyse to monodeboronated ionic compounds having [nido-7-SB10H11]- or [5-I-nido-7-SB10H10]- anions which are also products of the reactions of 1a and 1b with other N-bases such as pyridine, ammonia and DABCO. The extrusion of one boron and one sulphur atom takes place when 1a reacts with 2,6-di-tert-butylpyridine to yield decaborane.

Inside-protonated 1-azaadamantane: computational studies on the structure, stability, and generation

Post Hartree–Fock and density functional theory methods have been employed to study inside-protonated 1-azaadamantane 7 and its complexes with the fluoride counterion (contact ion pairs) 10 and 11. The study also involved 1-azaadamantane 4, its outside-protonated form 8, and 1-azaadamantane radical cation 17. Inside-protonated 1-azaadamantane 7 is more than 82 kcal mol−1less stable than out-isomer 8. The repulsive interaction between the internal N+–H group and the azaadamantane cage and a substantial deformation of this cage greatly weaken the C–N and C–C bonds and, consequently, lead to a low kinetic stability of in-ion 7 in the studied unimolecular and bimolecular reactions involving the removal of the encapsulated proton from the cage. Among these reactions, a 7 → 8 rearrangement through a reversible cage opening at the C–N bond was found to be the main transformation channel ([Formula: see text] < 16 kcal mol−1) for in-ion 7. This rearrangement can be catalyzed by an external base, e.g., the fluoride anion. A 1,4-hydrogen migration in 1-azaadamantane radical cation 17 as a possible pathway to the inside-protonated 1-azaadamantane 7 was explored. It was found that this process has a prohibitively high activation barrier, [Formula: see text] > 104 kcal mol−1, and is not able to compete with the α-C–C cleavage of the azaadamantane cage ([Formula: see text] < 26 kcal mol−1).

Relative glycosidic bond stabilities of naturally occurring methylguanosines: 7-methylation is intrinsically activating

The frequency and diversity of posttranscriptional modifications add an additional layer of chemical complexity beyond canonical nucleic acid sequence. Methylations are particularly frequently occurring and often highly conserved throughout the kingdoms of life. However, the intricate functions of these modified nucleic acid constituents are often not fully understood. Systematic foundational research that reduces systems to their minimum constituents may aid in unraveling the complexities of nucleic acid biochemistry. Here, we examine the relative intrinsic N-glycosidic bond stabilities of guanosine and five naturally occurring methylguanosines (O2'-, 1-, 7-, N2,N2-di-, and N2,N2,O2'-trimethylguanosine) probed by energy-resolved collision-induced dissociation tandem mass spectrometry and complemented with quantum chemical calculations. Apparent glycosidic bond stability is generally found to increase with increasing methyl substitution (canonical < mono- < di- < trimethylated). Many biochemical transformations, including base excision repair mechanisms, involve protonation and/or noncovalent interactions to increase nucleobase leaving-group ability. The protonated gas-phase methylguanosines require less activation energy for glycosidic bond cleavage than their sodium cationized forms. However, methylation at the N7 position intrinsically weakens the glycosidic bond of 7-methylguanosine more significantly than subsequent cationization, and thus 7-methylguanosine is suggested to be under perpetually activated conditions. N7 methylation also alters the nucleoside geometric preferences relative to the other systems, including the nucleobase orientation in the neutral form, sugar puckering in the protonated form, and the preferred protonation and sodium cation binding sites. All of the methylated guanosines examined here are predicted to have proton affinities and gas-phase basicities that exceed that of canonical guanosine. Additionally, the proton affinity and gas-phase basicity trends exhibit a roughly inverse correlation with the apparent glycosidic bond stabilities.

Binding preference of nitroimidazolic radiosensitizers to nucleobases and nucleosides probed by electrospray ionization mass spectrometry and density functional theory

Nitroimidazolic radiosensitizers are used in radiation therapy to selectively sensitize cancer cells deprived of oxygen, and the actual mechanism of radiosensitization is still not understood. Selecting five radiosensitizers (1-methyl-5-nitroimidazole, ronidazole, ornidazole, metronidazole, and nimorazole) with a common 5-nitroimidazolic ring with different substitutions at N1 and C2 positions of the imidazole moiety, we investigate here their binding to nucleobases (A, T, G, and C) and nucleosides (As, Td, Gs, and Cd) via the positive electrospray ionization mass spectrometry experiments. In addition, quantum chemical calculations at the M062x/6-311+G(d,p) level of theory and basis set were used to determine binding energies of the proton bound dimers of a radiosensitizer and a nucleobase. The positive electrospray ionization leads to the formation of proton bound dimers of all radiosensitizers except 1-methyl-5-nitroimidazole in high abundance with C and smaller abundance with G. Ronidazole and metronidazole formed less abundant dimers also with A, while no dimers were observed to be formed at all with T. In contrast to the case of the nucleoside Td, the dimer intensity is as high as that with Cd, while the abundance of the dimer with Gs is smaller than that of the former. The experimental results are consistent with the calculations of binding energies suggesting proton bound dimers with C and G to be the strongest bound ones. Finally, a barrier-free proton transfer is observed when protonated G or C approaches the nitroimidazole ring.