Acidity constants and protonation sites of cyclic dinucleotides determined by capillary electrophoresis, quantum chemical calculations, and NMR spectroscopy

Cyclic dinucleotides (CDNs) are important second messengers in bacteria and eukaryotes. Detailed characterization of their physicochemical properties is a prerequisite for understanding their biological functions. Herein, we examine acid-base and electromigration properties of selected CDNs employing capillary electrophoresis (CE), density functional theory (DFT), and nuclear magnetic resonance (NMR) spectroscopy to provide benchmark pKa values, as well as to unambiguously determine the protonation sites. Acidity constants (pKa) of the NH+ moieties of adenine and guanine bases and actual and limiting ionic mobilities of CDNs were determined by nonlinear regression analysis of the pH dependence of their effective electrophoretic mobilities measured by CE in aqueous background electrolytes in a wide pH range (0.98-11.48), at constant temperature (25°C), and constant ionic strength (25 mM). The thermodynamic pKa values were found to be in the range 3.31-4.56 for adenine and 2.28-3.61 for guanine bases, whereas the pKa of enol group of guanine base was in the range 10.21-10.40. Except for systematic shifts of ∼2 pKa, the pKa values calculated by the DFT-D3//COSMO-RS composite protocol that included large-scale conformational sampling and "cross-morphing" were in a relatively good agreement with the pKas determined by CE and predict N1 atom of adenine and N7 atom of guanine as the protonation sites. The protonation of the N1 atom of adenine and N7 atom of guanine in acidic background electrolytes (BGEs) and the dissociation of the enol group of guanine in alkaline BGEs was confirmed also by NMR spectroscopy.

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

What are the minimal folding seeds in proteins? Experimental and theoretical assessment of secondary structure propensities of small peptide fragments

Certain peptide sequences, some of them as short as amino acid triplets, are significantly overpopulated in specific secondary structure motifs in folded protein structures. For example, 74% of the EAM triplet is found in α-helices, and only 3% occurs in the extended parts of proteins (typically β-sheets). In contrast, other triplets (such as VIV and IYI) appear almost exclusively in extended parts (79% and 69%, respectively). In order to determine whether such preferences are structurally encoded in a particular peptide fragment or appear only at the level of a complex protein structure, NMR, VCD, and ECD experiments were carried out on selected tripeptides: EAM (denoted as pro-'α-helical' in proteins), KAM(α), ALA(α), DIC(α), EKF(α), IYI(pro-β-sheet or more generally, pro-extended), and VIV(β), and the reference α-helical CATWEAMEKCK undecapeptide. The experimental data were in very good agreement with extensive quantum mechanical conformational sampling. Altogether, we clearly showed that the pro-helical vs. pro-extended propensities start to emerge already at the level of tripeptides and can be fully developed at longer sequences. We postulate that certain short peptide sequences can be considered minimal "folding seeds". Admittedly, the inherent secondary structure propensity can be overruled by the large intramolecular interaction energies within the folded and compact protein structures. Still, the correlation of experimental and computational data presented herein suggests that the secondary structure propensity should be considered as one of the key factors that may lead to understanding the underlying physico-chemical principles of protein structure and folding from the first principles.

Experimental Evidence and Mechanistic Description of the Phenolic H-Transfer to the Cu2O2 Active Site of oxy-Tyrosinase

Tyrosinase is a ubiquitous coupled binuclear copper enzyme that activates O2 toward the regioselective monooxygenation of monophenols to catechols via a mechanism that remains only partially defined. Here, we present new mechanistic insights into the initial steps of this monooxygenation reaction by employing a pre-steady-state, stopped-flow kinetics approach that allows for the direct measurement of the monooxygenation rates for a series of para-substituted monophenols by oxy-tyrosinase. The obtained biphasic Hammett plot and the associated solvent kinetic isotope effect values provide direct evidence for an initial H-transfer from the protonated phenolic substrate to the Cu2O2 core of oxy-tyrosinase. The correlation of these experimental results to quantum mechanics/molecular mechanics calculations provides a detailed mechanistic description of this H-transfer step. These new mechanistic insights revise and expand our fundamental understanding of Cu2O2 active sites in biology.

Discovery of a Druggable, Cryptic Pocket in SARS-CoV-2 nsp16 Using Allosteric Inhibitors

A collaborative, open-science team undertook discovery of novel small molecule inhibitors of the SARS-CoV-2 nsp16-nsp10 2'-O-methyltransferase using a high throughput screening approach with the potential to reveal new inhibition strategies. This screen yielded compound 5a, a ligand possessing an electron-deficient double bond, as an inhibitor of SARS-CoV-2 nsp16 activity. Surprisingly, X-ray crystal structures revealed that 5a covalently binds within a previously unrecognized cryptic pocket near the S-adenosylmethionine binding cleft in a manner that prevents occupation by S-adenosylmethionine. Using a multidisciplinary approach, we examined the mechanism of binding of compound 5a to the nsp16 cryptic pocket and developed 5a derivatives that inhibited nsp16 activity and murine hepatitis virus replication in rat lung epithelial cells but proved cytotoxic to cell lines canonically used to examine SARS-CoV-2 infection. Our study reveals the druggability of this newly discovered SARS-CoV-2 nsp16 cryptic pocket, provides novel tool compounds to explore the site, and suggests a new approach for discovery of nsp16 inhibition-based pan-coronavirus therapeutics through structure-guided drug design.

Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6)

Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG^‡ ≈ 21 kcal·mol^-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF₃ or -CHF₂ substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.

Can Copper(I) and Silver(I) be Hydrogen Bond Acceptors?

Gold(I) centers can form moderately strong (Au⋅⋅⋅H) hydrogen bonds with tertiary ammonium groups, as has been demonstrated in the 3AuCl⁺ (3⁺ =1-(tert-butyl)-3-phenyl-4-(2-((dimethylammonio)methyl)phenyl)-1,2,4-triazol-5-ylidene) complex. However, similar hydrogen bonding interactions with isoelectronic silver(I) or copper(I) centers are unknown. Herein, we first explored whether the Au⋅⋅⋅H bond originally observed in 3AuCl⁺ can be strengthened by replacing Cl with Br or I. Experimental gas-phase IR spectra in the ∼3000 cm^-1 region showed only a small effect of the halogen on the Au⋅⋅⋅H bond. Next, we measured the spectra of 3AgCl⁺ , which exhibited significant differences compared to its 3AuX⁺ counterparts. The difference has been explained by DFT calculations which indicated that the Ag⋅⋅⋅H interaction is only marginal in this complex, and a Cl⋅⋅⋅H hydrogen bond is formed instead. Calculations predicted the same for the 3CuCl⁺ complex. However, we noticed that for Ag and Cu complexes with less flexible ligands, such as dimethyl(2-(dimethylammonio)phenyl)phosphine (L⁷ H⁺ ), the computations predict the presence of the respective Ag⋅⋅⋅H and Cu⋅⋅⋅H hydrogen bonds, with a strength similar to the Au⋅⋅⋅H bond in 3AuCl⁺ . We, therefore, propose possible complexes where the presence of (Ag/Cu)⋅⋅⋅H bonds could be experimentally verified to broaden our understanding of these unusual interactions.

Redox Potentials with COSMO-RS: Systematic Benchmarking with Different Databases

Recent techniques of computational electrochemistry can yield redox potentials with accuracy as good as 0.1 V. Yet, many such methods are not universal, easy to use, or computationally efficient. Herein, we provide a systematic benchmarking of a relatively cheap and straightforward computational approach for fairly accurate computations of redox potentials. It is based on a combination of the conductor-like screening model for real solvents (COSMO-RS) and the density functional theory (DFT). The benchmarking is done with databases covering diverse redox systems, including transition-metal aquacomplexes and various organic and inorganic compounds. In addition, we also present our own test set aiming at maximum chemical diversity and maximum range of redox potential values. We assess the performance of the fairly efficient computational protocol combining the COSMO-RS with the BP86 DFT functional. This is done by calibrating it against different high-level state-of-the-art techniques, in particular, polarizable continuum model (PCM) connected to composite G3(MP2,CC)(+) method, domain-based pair natural orbital implementation of coupled cluster theory, or complete basis set CBS-QB3 method. We also elaborate on the absolute reduction potential value of standard hydrogen electrode to be used with COSMO-RS, and we propose the value of approx. 4.4 V. The COSMO-RS/BP86-D3(BJ) combination outperforms other methods on a wide range of redox systems. However, we show that its accuracy is based on a balanced error cancelation and, therefore, it cannot be further systematically improved. As a result, the proposed procedure represents a pragmatic choice for large-scale screening, yet it could be combined with more advanced methods for detailed studies.

Molecular dynamics simulations provide structural insight into binding of cyclic dinucleotides to human STING protein

Human stimulator of interferon genes (hSTING) is a signaling adaptor protein that triggers innate immune system by response to cytosolic DNA and second messenger cyclic dinucleotides (CDNs). Natural CDNs contain purine nucleobase with different phosphodiester linkage types (3'-3', 2'-2' or mixed 2'-3'-linkages) and exhibit different binding affinity towards hSTING, ranging from micromolar to nanomolar. High-affinity CDNs are considered as suitable candidates for treatment of chronic hepatitis B and cancer. We have used molecular dynamics simulations to investigate dynamical aspects of binding of natural CDNs (specifically, 2'-2'-cGAMP, 2'-3'-cGAMP, 3'-3'-cGAMP, 3'-3'-c-di-AMP, and 3'-3'-c-di-GMP) with hSTING^wt protein. Our results revealed that CDN/hSTING^wt interactions are controlled by the balance between fluctuations (conformational changes) in the CDN ligand and the protein dynamics. Binding of different CDNs induces different degrees of conformational/dynamics changes in hSTING^wt ligand binding cavity, especially in α₁-helices, the so-called lid region and α₂-tails. The ligand residence time in hSTING^wt protein pocket depends on different contribution of R232 and R238 residues interacting with oxygen atoms of phosphodiester groups in ligand, water distribution around interacting charged centers (in protein residues and ligand) and structural stability of closed conformation state of hSTING^wt protein. These findings may perhaps guide design of new compounds modulating hSTING activity.Communicated by Ramaswamy H. Sarma.

Cationic Gold(II) Complexes: Experimental and Theoretical Study

Gold(II) complexes are rare, and their application to the catalysis of chemical transformations is underexplored. The reason is their easy oxidation or reduction to more stable gold(III) or gold(I) complexes, respectively. We explored the thermodynamics of the formation of [Au^II (L)(X)]⁺ complexes (L=ligand, X=halogen) from the corresponding gold(III) precursors and investigated their stability and spectral properties in the IR and visible range in the gas phase. The results show that the best ancillary ligands L for stabilizing gaseous [Au^II (L)(X)]⁺ complexes are bidentate and tridentate ligands with nitrogen donor atoms. The electronic structure and spectral properties of the investigated gold(II) complexes were correlated with quantum chemical calculations. The results show that the molecular and electronic structure of the gold(II) complexes as well as their spectroscopic properties are very similar to those of analogous stable copper(II) complexes.

Elucidation of the tyrosinase/O2/monophenol ternary intermediate that dictates the monooxygenation mechanism in melanin biosynthesis

Melanins are highly conjugated biopolymer pigments that provide photoprotection in a wide array of organisms, from bacteria to humans. The rate-limiting step in melanin biosynthesis, which is the ortho-hydroxylation of the amino acid L-tyrosine to L-DOPA, is catalyzed by the ubiquitous enzyme tyrosinase (Ty). Ty contains a coupled binuclear copper active site that binds O2 to form a μ:η2:η2-peroxide dicopper(II) intermediate (oxy-Ty), capable of performing the regioselective monooxygenation of para-substituted monophenols to catechols. The mechanism of this critical monooxygenation reaction remains poorly understood despite extensive efforts. In this study, we have employed a combination of spectroscopic, kinetic, and computational methods to trap and characterize the elusive catalytic ternary intermediate (Ty/O2/monophenol) under single-turnover conditions and obtain molecular-level mechanistic insights into its monooxygenation reactivity. Our experimental results, coupled with quantum-mechanics/molecular-mechanics calculations, reveal that the monophenol substrate docks in the active-site pocket of oxy-Ty fully protonated, without coordination to a copper or cleavage of the μ:η2:η2-peroxide O-O bond. Formation of this ternary intermediate involves the displacement of active-site water molecules by the substrate and replacement of their H bonds to the μ:η2:η2-peroxide by a single H bond from the substrate hydroxyl group. This H-bonding interaction in the ternary intermediate enables the unprecedented monooxygenation mechanism, where the μ-η2:η2-peroxide O-O bond is cleaved to accept the phenolic proton, followed by substrate phenolate coordination to a copper site concomitant with its aromatic ortho-hydroxylation by the nonprotonated μ-oxo. This study provides insights into O2 activation and reactivity by coupled binuclear copper active sites with fundamental implications in biocatalysis.

Exhaustive Mapping of the Conformational Space of Natural Dipeptides by the DFT-D3//COSMO-RS Method

We extensively mapped energy landscapes and conformations of 22 (including three His protonation states) proteinogenic α-amino acids in trans configuration and the corresponding 484 (22²) dipeptides. To mimic the environment in a protein chain, the N- and C-termini of the studied systems were capped with acetyl and N-methylamide groups, respectively. We systematically varied the main chain dihedral angles (ϕ, ψ) by 40° steps and all side chain angles by 90° or 120° steps. We optimized the molecular geometries with the GFN2-xTB semiempirical (SQM) method and performed single point density functional theory calculations at the BP86-D3/DGauss-DZVP//COSMO-RS level in water, 1-octanol, N,N-dimethylformamide, and n-hexane. For each restrained (nonequilibrium) structure, we also calculated energy gradients (in water) and natural atomic charges. The exhaustive and unprecedented QM-based sampling enabled us to construct Ramachandran plots of quantum mechanical (QM(BP86-D3)//COSMO-RS) energies calculated on SQM structures, for all 506 (484 dipeptides and 22 amino acids) studied systems. We showed how the character of an amino acid side chain influences the conformational space of single amino acids and dipeptides. With clustering techniques, we were able to identify unique minima of amino acids and dipeptides (i.e., minima on the GFN2-xTB potential energy surfaces) and analyze the distribution of their BP86-D3//COSMO-RS conformational energies in all four solvents. We also derived an empirical formula for the number of unique minima based on the overall number of rotatable bonds within each peptide. The final peptide conformer data set (PeptideCs) comprises over 400 million structures, all of them annotated with QM(BP86-D3)//COSMO-RS energies. Thanks to its completeness and unbiased nature, the PeptideCs can serve, inter alia, as a data set for the validation of new methods for predicting the energy landscapes of protein structures. This data set may also prove to be useful in the development and reparameterization of biomolecular force fields. The data set is deposited at Figshare (10.25452/figshare.plus.19607172) and can be accessed using a simple web interface at http://peptidecs.uochb.cas.cz.

Cyclic Octapeptides Composed of Two Glutathione Units Outperform the Monomer in Lead Detoxification

A rationally-designed scaffold of cyclic octapeptides composed of two units of the natural tripeptide glutathione (GSH) was optimized to strongly and selectively capture toxic lead ions (Pb(II)). Using state-of-the-art computational tools, a list of eleven plausible peptides was shortened to five analogs based on their calculated affinity to Pb(II) ions. We then synthesized and investigated them for their abilities to recover Pb-poisoned human cells. A clear pattern was observed from the in vitro detoxification results, indicating the importance of cavity size and polar moieties to enhance metal capturing. These, together with the apparent benefit in cyclizing the peptides, improved the detoxification of the two lead peptides by approximately two folds compared to GSH and the benchmark chelating agents against Pb poisoning. Moreover, the two peptides did not show any toxicity and, therefore, were thoroughly investigated to determine their potential as next-generation remedies for Pb poisoning.

Comprehensive Theoretical View of the [Cu2O2] Side-on-Peroxo-/Bis-μ-Oxo Equilibria

Coupled binuclear copper (CBC) sites are employed by many metalloenzymes to catalyze a broad set of biochemical transformations. Typically, the CBC catalytic sites are activated by the O2 molecule to form various [Cu2O2] reactive species. This has also inspired synthesis and development of various biomimetic inorganic complexes featuring the CBC core. From theoretical perspective, the [Cu2O2] reactivity often hinges on the side-on-peroxo-dicopper(II) (P) vs. bis-μ-oxo-dicopper(III) (O) isomerism - an equilibrium that has become almost iconic in theoretical bioinorganic chemistry. Herein, we present a comprehensive calibration and evaluation of the performance of various composite computational protocols available in contemporary computational chemistry, involving coupled-cluster and multi-reference (relativistic) wave function methods, popular density functionals and solvation models. Starting with the well-studied reference [Cu2O2(NH3)6]2+ system, we compared the performance of electronic structure methods and discussed the relativistic effects. This allowed us to select several 'calibrated' DFT functionals that can be conveniently employed to study ten experimentally well-characterized [Cu2O2] inorganic systems. We mostly correctly predicted the lowest-energy structures (P vs. O) of the studied systems. In addition, we present calibration of the used electronic structure methods for prediction of the spectroscopic features of the [Cu2O2] core.

Evidence for H-bonding interactions to the μ-η2:η2-peroxide of oxy-tyrosinase that activate its coupled binuclear copper site

The factors that control the diverse reactivity of the μ-η²:η²-peroxo dicopper(II) oxy-intermediates in the coupled binuclear copper proteins remain elusive. Here, spectroscopic and computational methods reveal H-bonding interactions between active-site waters and the μ-η²:η²-peroxide of oxy-tyrosinase, and define their effects on the Cu(II)₂O₂ electronic structure and O₂ activation.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Understanding desaturation/hydroxylation activity of castor stearoyl Δ9-Desaturase through rational mutagenesis

•

Rationally designed mutations in the Δ⁹ desaturase promoted hydroxylation activity.

•

Proton and electron transfer to the active site is crucial for the Δ⁹D to desaturate

•

Detailed analysis of all enzymatic products of the Δ⁹D was carried out

•

Insight into the chemo-, and stereoselectivity of non-heme diiron enzymes was obtained

A recently proposed reaction mechanism of soluble Δ⁹ desaturase (Δ⁹D) allowed us to identify auxiliary residues His203, Asp101, Thr206 and Cys222 localized near the di-iron active site that are supposedly involved in the proton transfer (PT) to and from the active site. The PT, along with the electron transfer (ET), seems to be crucial for efficient desaturation. Thus, perturbing the major PT chains is expected to impair the native reaction and (potentially) amplify minor reaction channels, such as the substrate hydroxylation. To verify this hypothesis, we mutated the four residues mentioned above into their counterparts present in a soluble methane monooxygenase (sMMO), and determined the reaction products of mutants. We found that the mutations significantly promote residual monohydroxylation activities on stearoyl-CoA, often at the expense of native desaturation activity. The favored hydroxylation positions are C₉, followed by C₁₀ and C₁₁. Reactions with unsaturated substrate, oleoyl-CoA, yield erythro-9,10-diol, cis-9,10-epoxide and a mixture of allylic alcohols. Additionally, using 9- and 11-hydroxystearoyl-CoA, we showed that the desaturation reaction can proceed only with the hydroxyl group at position C₁₁, whereas the hydroxylation reaction is possible in both cases, i.e. with hydroxyl at position C₉ or C₁₁. Despite the fact that the overall outcome of hydroxylation is rather modest and that it is mostly the desaturation/hydroxylation ratio that is affected, our results broaden understanding of the origin of chemo- and stereoselectivity of the Δ⁹D and provide further insight into the catalytic action of the NHFe₂ enzymes.

Thiolation and Carboxylation of Glutathione Synergistically Enhance Its Lead-Detoxification Capabilities

The natural tripeptide glutathione (GSH) is a ubiquitous compound harboring various biological tasks, among them interacting with essential and toxic metal ions. Yet, although weakly binding the poisonous metal lead (Pb), GSH poorly detoxifies it. β-Mercaptoaspartic acid is a new-to-nature novel amino acid that was found to enhance the Pb-detoxification capability of a synthetic cyclic tetrapeptide. Aiming to explore the advantages of noncanonical amino acids (ncAAs) of this nature, we studied the detoxification capabilities of GSH and three analogue peptides, each of which contains at least one ncAA that harbors both free carboxylate and thiolate groups. A thorough investigation that includes in vitro detoxification and mechanistic evaluations, metal-binding affinity, metal selectivity, and computational studies shows that these ncAAs are highly beneficial in additively enhancing Pb binding and reveals the importance of both high affinity and metal selectivity in synergistically reducing Pb toxicity in cells. Hence, such ncAAs join the chemical toolbox against Pb poisoning and pollution, enabling peptides to strongly and selectively bind the toxic metal ion.

Enzymatic Synthesis of 3'-5', 3'-5' Cyclic Dinucleotides, Their Binding Properties to the Stimulator of Interferon Genes Adaptor Protein, and Structure/Activity Correlations

The 3'-5', 3'-5' cyclic dinucleotides (3'3'CDNs) are bacterial second messengers that can also bind to the stimulator of interferon genes (STING) adaptor protein in vertebrates and activate the host innate immunity. Here, we profiled the substrate specificity of four bacterial dinucleotide synthases from Vibrio cholerae (DncV), Bacillus thuringiensis (btDisA), Escherichia coli (dgcZ), and Thermotoga maritima (tDGC) using a library of 33 nucleoside-5'-triphosphate analogues and then employed these enzymes to synthesize 24 3'3'CDNs. The STING affinity of CDNs was evaluated in cell-based and biochemical assays, and their ability to induce cytokines was determined by employing human peripheral blood mononuclear cells. Interestingly, the prepared heterodimeric 3'3'CDNs bound to the STING much better than their homodimeric counterparts and showed similar or better potency than bacterial 3'3'CDNs. We also rationalized the experimental findings by in-depth STING-CDN structure-activity correlations by dissecting computed interaction free energies into a set of well-defined and intuitive terms. To this aim, we employed state-of-the-art methods of computational chemistry, such as quantum mechanics/molecular mechanics (QM/MM) calculations, and complemented the computed results with the {STING:3'3'c-di-ara-AMP} X-ray crystallographic structure. QM/MM identified three outliers (mostly homodimers) for which we have no clear explanation of their impaired binding with respect to their heterodimeric counterparts, whereas the R² = 0.7 correlation between the computed ΔG'_{int_rel} and experimental ΔT_m's for the remaining ligands has been very encouraging.

Ir-Catalyzed Cycloaddition of Tribenzocyclyne with Biphenylenes

We demonstrate that Ir-catalyzed C-C bond activation in biphenylenes followed by a reaction with tribenzocyclyne is a suitable method for synthesizing strained and unknown monoadducts with the tetradehydrotetrabenzo[a,c,e,i]cyclododecene scaffold ([12]annulenes). Modification of reaction conditions also furnished [12]annulene products with cis and/or trans double bonds formed by hydrogen transfer. The [9]annulene side product was formed upon the reaction of the benzyl radical with tribenzocyclyne during the Bergman cyclization. All isolated compounds were fully characterized by HRMS, NMR, and X-ray diffraction analysis.

Potent Cyclic Tetrapeptide for Lead Detoxification

Lead (Pb) is a ubiquitous poisonous metal, affecting the health of vast populations worldwide. Medications to treat Pb poisoning suffer from various limitations and are often toxic owing to insufficient metal selectivity. Here, we report a cyclic tetrapeptide that selectively binds Pb and eradicates its toxic effect on the cellular level, with superior potency than state-of-the-art drugs. The Pb-peptide complex is remarkably strong and was characterized experimentally and computationally. Accompanied by the lack of toxicity and enhanced stability of this peptide, these qualities indicate its merit as a potential remedy for Pb poisoning.

Synthesis and Biological Evaluation of Phosphoester and Phosphorothioate Prodrugs of STING Agonist 3',3'-c-Di(2'F,2'dAMP)

Cyclic dinucleotides (CDNs) are second messengers that bind to the stimulator of interferon genes (STING) and trigger the expression of type I interferons and proinflammatory cytokines. Here we evaluate the activity of 3',3'-c-di(2'F,2'dAMP) and its phosphorothioate analogues against five STING allelic forms in reporter-cell-based assays and rationalize our findings with X-ray crystallography and quantum mechanics/molecular mechanics calculations. We show that the presence of fluorine in the 2' position of 3',3'-c-di(2'F,2'dAMP) improves its activity not only against the wild type (WT) but also against REF and Q STING. Additionally, we describe the synthesis of the acyloxymethyl and isopropyloxycarbonyl phosphoester prodrugs of CDNs. Masking the negative charges of the CDNs results in an up to a 1000-fold improvement of the activities of the prodrugs relative to those of their parent CDNs. Finally, the uptake and intracellular cleavage of pivaloyloxymethyl prodrugs to the parent CDN is rapid, reaching a peak intracellular concentration within 2 h.

Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein

STING (stimulator of interferon genes) is a key regulator of innate immunity that has recently been recognized as a promising drug target. STING is activated by cyclic dinucleotides (CDNs) which eventually leads to expression of type I interferons and other cytokines. Factors underlying the affinity of various CDN analogues are poorly understood. Herein, we correlate structural biology, isothermal calorimetry (ITC) and computational modeling to elucidate factors contributing to binding of six CDNs-three pairs of natural (ribo) and fluorinated (2'-fluororibo) 3',3'-CDNs. X-ray structural analyses of six {STING:CDN} complexes did not offer any explanation for the different affinities of the studied ligands. ITC showed entropy/enthalpy compensation up to 25 kcal mol-1 for this set of similar ligands. The higher affinities of fluorinated analogues are explained with help of computational methods by smaller loss of entropy upon binding and by smaller strain (free) energy.

Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments

Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.

Plasmon-assisted click chemistry at low temperature: an inverse temperature effect on the reaction rate

Plasmon assistance promotes a range of chemical transformations by decreasing their activation energies. In a common case, thermal and plasmon assistance work synergistically: higher temperature results in higher plasmon-enhanced catalysis efficiency. Herein, we report an unexpected tenfold increase in the reaction efficiency of surface plasmon-assisted Huisgen dipolar azide-alkyne cycloaddition (AAC) when the reaction mixture is cooled from room temperature to -35 °C. We attribute the observed increase in the reaction efficiency to complete plasmon-induced annihilation of the reaction barrier, prolongation of plasmon lifetime, and decreased relaxation of plasmon-excited-states under cooling. Furthermore, control quenching experiments supported by theoretical calculations indicate that plasmon-mediated substrate excitation to an electronic triplet state may play the key role in plasmon-assisted chemical transformation. Last but not least, we demonstrated the possible applicability of plasmon assistance to biological systems by AAC coupling of biotin to gold nanoparticles performed at -35 °C.

Protein-Ligand Interactions in the STING Binding Site Probed by Rationally Designed Single-Point Mutations: Experiment and Theory

STING protein (stimulator of interferon genes) plays an important role in the innate immune system. A number of potent compounds regulating its activity have been reported, mostly derivatives of cyclic dinucleotides (CDNs), natural STING agonists. Here, we aim to provide complementary information to large-scale "ligand-profiling" studies by probing the importance of STING-CDN protein-ligand interactions on the protein side. We examined in detail six typical CDNs each in complex with 13 rationally devised mutations in STING: S162A, S162T, Y167F, G230A, R232K, R232H, A233L, A233I, R238K, T263A, T263S, R293Q, and G230A/R293Q. The mutations switch on and off various types of protein-ligand interactions: π-π stacking, hydrogen bonding, ionic pairing, and nonpolar contacts. We correlated experimental data obtained by differential scanning fluorimetry, X-ray crystallography, and isothermal titration calorimetry with theoretical calculations. This enabled us to provide a mechanistic interpretation of the differences in the binding of representative CDNs to STING. We observed that the G230A mutation increased the thermal stability of the protein-ligand complex, indicating an increased level of ligand binding, whereas R238K and Y167F led to a complete loss of stabilization (ligand binding). The effects of the other mutations depended on the type of ligand (CDN) and varied, to some extent. A very good correlation (R² = 0.6) between the experimental binding affinities and interaction energies computed by quantum chemical methods enabled us to explain the effect of the studied mutations in detail and evaluate specific interactions quantitatively. Our work may inspire development of high-affinity ligands against the common STING haplotypes by targeting the key (sometimes non-intuitive) protein-ligand interactions.

Mapping Conformational Space of All 8000 Tripeptides by Quantum Chemical Methods: What Strain Is Affordable within Folded Protein Chains?

To gain more insight into the physicochemical aspects of a protein structure from the first principles, conformational space of all 8000 "capped" tripeptides (i.e., N-Ac-X₁X₂X₃-NH-CH₃, where X_i is one of the 20 natural amino acids) was investigated computationally. An enormous dataset (denoted P-CONF_1.6M and containing close to 1 600 000 conformers in total) has been obtained by employing a composite protocol combining density functional theory, semiempirical quantum mechanics (SQM), and state-of-the-art solvation methods with 1000 K molecular dynamics (MD) used to generate initial structures (200 snapshots for each tripeptide). This allowed us to present the first rigorous QM-based glimpse at the vast conformational space spanned by small protein fragments. The same computational procedure was repeated for tripeptide fragments taken from the SCOPe database of three-dimensional protein folds, by restraining them to their geometry in a protein. Such complementary data allowed us to compare the distribution of conformational strain energies of unrestrained tripeptidic fragments "in solvent" with those in existing protein chains. Besides providing a rigorous (ab initio) proof of a few well-known concepts and hypotheses concerning protein structures, such as the distribution of (φ, ψ) angles in Ramachandran plots, we have made several observations that came as a certain surprise: (1) distribution of conformational energies does not significantly differ between the "unbiased/unrestrained" conformers obtained from MD sampling in solvent and the biased conformers, i.e., those of a given tripeptide obtained from protein structures; (2) conformational (strain) energy window up to ∼20 to 25 kcal·mol^-1 is readily available to tripeptide fragments within the context of a protein chain; (3) overpopulation in certain regions of Ramachandran plot was observed for the unbiased conformers. Last but not least, the massive dataset of accurate (DFT-D3//COSMO-RS) conformational (free) energies of ∼1.6 M peptide conformers, P-CONF_1.6M, obtained throughout this work may serve as excellent dataset for calibrating and benchmarking of popular force fields.

Closed Shell Iron(IV) Oxo Complex with an Fe-O Triple Bond: Computational Design, Synthesis, and Reactivity

Iron(IV)-oxo intermediates in nature contain two unpaired electrons in the Fe-O antibonding orbitals, which are thought to contribute to their high reactivity. To challenge this hypothesis, we designed and synthesized closed-shell singlet iron(IV) oxo complex [(quinisox)Fe(O)]+ (1+ ; quinisox-H=(N-(2-(2-isoxazoline-3-yl)phenyl)quinoline-8-carboxamide). We identified the quinisox ligand by DFT computational screening out of over 450 candidates. After the ligand synthesis, we detected 1+ in the gas phase and confirmed its spin state by visible and infrared photodissociation spectroscopy (IRPD). The Fe-O stretching frequency in 1+ is 960.5 cm-1 , consistent with an Fe-O triple bond, which was also confirmed by multireference calculations. The unprecedented bond strength is accompanied by high gas-phase reactivity of 1+ in oxygen atom transfer (OAT) and in proton-coupled electron transfer reactions. This challenges the current view of the spin-state driven reactivity of the Fe-O complexes.

Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ9-Desaturase

A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe₂) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe₂ family of enzymes: soluble Δ⁹ desaturase (Δ⁹D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[Fe^IV]₂ unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ⁹D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe₂ chemistry, provided by an amino acid sequence analysis through the families of the NHFe₂ enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe₂ enzymes and may inspire further work in NHFe₍₂₎ biomimetic chemistry.

Interplay between Conformational Strain and Intramolecular Interaction in Protein Structures: Which of Them Is Evolutionarily Conserved?

By computing strain energies of peptide fragments within protein structures and their intramolecular interaction energies, we attempt to reveal general biophysical trends behind the secondary structure formation in the context of protein evolution. Our "protein basis set" consisted of 1143 representatives of different folds obtained from curated SCOPe database, and for each member of the set, the strain and intramolecular energy was calculated on the "rolling tripeptide" basis, employing the DFT-D3/COSMO-RS method for the former and the QM-calibrated force field method (MM) for the latter. The calculated data, strain and interactions, were correlated with the conservation of amino acid residues in secondary structure elements and also with the level of the residue burial within the protein three-dimensional structure. It allowed us to formulate several observations concerning fundamental differences between two main secondary structure motifs: α-helices and β-strands. We have shown that a strong interaction is one of the determining characteristics of the β-sheet formation, at least at the level of tripeptides (and likely penta- or heptapeptides, too), and that the β-strand is a prevailing secondary structure in the strongly-interacting regions of the protein folds conserved by evolution. On the other hand, low strain was neither proven to be an important physicochemical property conserved by evolution nor does it correlate with the propensity for the α-helix and β-strand. Finally, it has been demonstrated that the strong interaction has a certain level of connection with residue burial; however, we demonstrate that these two characteristics should be rather regarded as two complementary factors. These findings represent an important contribution to understanding protein folding from first principles, which is a complementary approach to ongoing efforts to solve the protein folding problem by knowledge-based approaches and machine-learning.

Dinucleoside polyphosphates act as 5'-RNA caps in bacteria

It has been more than 50 years since the discovery of dinucleoside polyphosphates (Np_nNs) and yet their roles and mechanisms of action remain unclear. Here, we show that both methylated and non-methylated Np_nNs serve as RNA caps in Escherichia coli. Np_nNs are excellent substrates for T7 and E. coli RNA polymerases (RNAPs) and efficiently initiate transcription. We demonstrate, that the E. coli enzymes RNA 5′-pyrophosphohydrolase (RppH) and bis(5′-nucleosyl)-tetraphosphatase (ApaH) are able to remove the Np_nN-caps from RNA. ApaH is able to cleave all Np_nN-caps, while RppH is unable to cleave the methylated forms suggesting that the methylation adds an additional layer to RNA stability regulation. Our work introduces a different perspective on the chemical structure of RNA in prokaryotes and on the role of RNA caps. We bring evidence that small molecules, such as Np_nNs are incorporated into RNA and may thus influence the cellular metabolism and RNA turnover.

Nicotinamide adenine dinucleotide and coenzyme A serve as a 5′-cap of prokaryotic RNA. Here the authors report that methylated and non-methylated dinucleoside polyphosphates (Np_nNs) exist as Escherichia coli RNA caps which can be cleaved by 5′-pyrophosphohydrolase (RppH) and bis(5′-nucleosyl)-tetraphosphatase (ApaH).

Enzymatic Preparation of 2'-5',3'-5'-Cyclic Dinucleotides, Their Binding Properties to Stimulator of Interferon Genes Adaptor Protein, and Structure/Activity Correlations

Cyclic dinucleotides are second messengers in the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, which plays an important role in recognizing tumor cells and viral or bacterial infections. They bind to the STING adaptor protein and trigger expression of cytokines via TANK binding kinase 1 (TBK1)/interferon regulatory factor 3 (IRF3) and inhibitor of nuclear factor-κB (IκB) kinase (IKK)/nuclear factor-κB (NFκB) signaling cascades. In this work, we describe an enzymatic preparation of 2'-5',3'-5'-cyclic dinucleotides (2'3'CDNs) with use of cyclic GMP-AMP synthases (cGAS) from human, mouse, and chicken. We profile substrate specificity of these enzymes by employing a small library of nucleotide-5'-triphosphate (NTP) analogues and use them to prepare 33 2'3'CDNs. We also determine affinity of these CDNs to five different STING haplotypes in cell-based and biochemical assays and describe properties needed for their optimal activity toward all STING haplotypes. Next, we study their effect on cytokine and chemokine induction by human peripheral blood mononuclear cells (PBMCs) and evaluate their cytotoxic effect on monocytes. Additionally, we report X-ray crystal structures of two new CDNs bound to STING protein and discuss structure-activity relationship by using quantum and molecular mechanical (QM/MM) computational modeling.

Factors Stabilizing β-Sheets in Protein Structures from a Quantum-Chemical Perspective

Protein folds are determined by the interplay between various (de)stabilizing forces, which can be broadly divided into a local strain of the protein chain and intramolecular interactions. In contrast to the α-helix, the β-sheet secondary protein structure is significantly stabilized by long-range interactions between the individual β-strands. It has been observed that quite diverse amino acid sequences can form a very similar small β-sheet fold, such as in the three-β-strand WW domain. Employing "calibrated" quantum-chemical methods, we show herein on two sequentially diverse examples of the WW domain that the internal strain energy is higher in the β-strands and lower in the loops, while the interaction energy has an opposite trend. Low strain energy computed for peptide sequences in the loop 1 correlates with its postulated early formation in the folding process. The relatively high strain energy within the β-strands (up to 8 kcal mol^-1 per amino acid residue) is compensated by even higher intramolecular interaction energy (up to 15 kcal mol^-1 per residue). It is shown in a quantitative way that the most conserved residues across the structural family of WW domains have the highest contributions to the intramolecular interaction energy. On the other hand, the residues in the regions with the lowest strain are not conserved. We conclude that the internal interaction energy is the physical quantity tuned by evolution to define the β-sheet protein fold.

Toward Ab Initio Protein Folding: Inherent Secondary Structure Propensity of Short Peptides from the Bioinformatics and Quantum-Chemical Perspective

By combining bioinformatics with quantum-chemical calculations, we attempt to address quantitatively some of the physical principles underlying protein folding. The former allowed us to identify tripeptide sequences in existing protein three-dimensional structures with a strong preference for either helical or extended structure. The selected representatives of pro-helical and pro-extended sequences were converted into "isolated" tripeptides-capped at N- and C-termini-and these were subjected to an extensive conformational sampling and geometry optimization (typically thousands to tens of thousands of conformers for each tripeptide). For each conformer, the QM(DFT-D3)/COSMO-RS free-energy value was then calculated, G_conf(solv). The Δ G_conf(solv) is expected to provide an objective, unbiased, and quantitatively accurate measure of the conformational preference of the particular tripeptide sequence. It has been shown that irrespective of the helical vs extended preferences of the selected tripeptide sequences in context of the protein, most of the low-energy conformers of isolated tripeptides prefer the R-helical structure. Nevertheless, pro-helical tripeptides show slightly stronger helix preference than their pro-extended counterparts. Furthermore, when the sampling is repeated in the presence of a partner tripeptide to mimic the situation in a β-sheet, pro-extended tripeptides (exemplified by the VIV) show a larger free-energy benefit than pro-helical tripeptides (exemplified by the EAM). This effect is even more pronounced in a hydrophobic solvent, which mimics the less polar parts of a protein. This is in line with our bioinformatic results showing that the majority of pro-extended tripeptides are hydrophobic. The preference for a specific secondary structure by the studied tripeptides is thus governed by the plasticity to adopt to its environment. In addition, we show that most of the "naturally occurring" conformations of tripeptide sequences, i.e., those found in existing three-dimensional protein structures, are within ∼10 kcal·mol^-1 from their global minima. In summary, our "ab initio" data suggest that complex protein structures may start to emerge already at the level of their small oligopeptidic units, which is in line with a hierarchical nature of protein folding.

Spectroscopic and Computational Evidence of Intramolecular AuI ⋅⋅⋅H+ -N Hydrogen Bonding

Despite substantial evidence of short Au⋅⋅⋅H−X contacts derived from a number of X‐ray structures of Au^I compounds, the nature of Au^I⋅⋅⋅H bonding in these systems has not been clearly understood. Herein, we present the first spectroscopic evidence for an intramolecular Au^I⋅⋅⋅H⁺−N hydrogen bond in a [Cl−Au−L]⁺ complex, where L is a protonated N‐heterocyclic carbene. The complex was isolated in the gas phase and characterized with helium‐tagging infrared photodissociation (IRPD) spectra, in which H⁺−N‐mode‐derived bands evidence the intramolecular Au^I⋅⋅⋅H⁺−N bond. Quantum chemical calculations reproduce the experimental IRPD spectra and allow to characterize the intramolecular Au⋅⋅⋅H⁺−N bonding with a short r_Au⋅⋅⋅H distance of 2.17 Å and an interaction energy of approximately −10 kcal mol⁻¹. Various theoretical descriptors of chemical bonding calculated for the Au⋅⋅⋅H⁺−N interaction provide strong evidence for a hydrogen bond of moderate strength.

Structural and computational basis for potent inhibition of glutamate carboxypeptidase II by carbamate-based inhibitors

A series of carbamate-based inhibitors of glutamate carboxypeptidase II (GCPII) were designed and synthesized using ZJ-43, N-[[[(1S)-1-carboxy-3-methylbutyl]amino]carbonyl]-l-glutamic acid, as a molecular template in order to better understand the impact of replacing one of the two nitrogen atoms in the urea-based GCPII inhibitor with an oxygen atom. Compound 7 containing a C-terminal 2-oxypentanedioic acid was more potent than compound 5 containing a C-terminal glutamic acid (2-aminopentanedioic acid) despite GCPII's preference for peptides containing an N-terminal glutamate as substrates. Subsequent crystallographic analysis revealed that ZJ-43 and its two carbamate analogs 5 and 7 with the same (S,S)-stereochemical configuration adopt a nearly identical binding mode while (R,S)-carbamate analog 8 containing a d-leucine forms a less extensive hydrogen bonding network. QM and QM/MM calculations have identified no specific interactions in the GCPII active site that would distinguish ZJ-43 from compounds 5 and 7 and attributed the higher potency of ZJ-43 and compound 7 to the free energy changes associated with the transfer of the ligand from bulk solvent to the protein active site as a result of the lower ligand strain energy and solvation/desolvation energy. Our findings underscore a broader range of factors that need to be taken into account in predicting ligand-protein binding affinity. These insights should be of particular importance in future efforts to design and develop GCPII inhibitors for optimal inhibitory potency.

Photodissociative Cross-Linking of Non-covalent Peptide-Peptide Ion Complexes in the Gas Phase

We report a gas-phase UV photodissociation study investigating non-covalent interactions between neutral hydrophobic pentapeptides and peptide ions incorporating a diazirine-tagged photoleucine residue. Phenylalanine (Phe) and proline (Pro) were chosen as the conformation-affecting residues that were incorporated into a small library of neutral pentapeptides. Gas-phase ion-molecule complexes of these peptides with photo-labeled pentapeptides were subjected to photodissociation. Selective photocleavage of the diazirine ring at 355 nm formed short-lived carbene intermediates that underwent cross-linking by insertion into H-X bonds of the target peptide. The cross-link positions were established from collision-induced dissociation tandem mass spectra (CID-MS³) providing sequence information on the covalent adducts. Effects of the amino acid residue (Pro or Phe) and its position in the target peptide sequence were evaluated. For proline-containing peptides, interactions resulting in covalent cross-links in these complexes became more prominent as proline was moved towards the C-terminus of the target peptide sequence. The photocross-linking yields of phenylalanine-containing peptides depended on the position of both phenylalanine and photoleucine. Density functional theory calculations were used to assign structures of low-energy conformers of the (GLPMG + GLL*LK + H)⁺ complex. Born-Oppenheimer molecular dynamics trajectory calculations were used to capture the thermal motion in the complexes within 100 ps and determine close contacts between the incipient carbene and the H-X bonds in the target peptide. This provided atomic-level resolution of potential cross-links that aided spectra interpretation and was in agreement with experimental data. Graphical Abstract ᅟ.