Interplay of hydrogen bonds and n→π* interactions in proteins

An efficient method by combining different basis sets and SAPT levels

In symmetry-adapted perturbation theory (SAPT), accurate calculations on non-covalent interaction (NCI) for large complexes with more than 50 atoms are time-consuming using large basis sets. More efficient ones with smaller basis sets usually result in poor prediction in terms of dispersion and overall energies. In this study, we propose two composite methods with baseline calculated at SAPT2/aug-cc-pVDZ and SAPT2/aug-cc-pVTZ with dispersion term corrected at SAPT2+ level using bond functions and smaller basis set with δ MP2 corrections respectively. Benchmark results on representative NCI data sets, such as S22, S66, and so forth, show significant improvements on the accuracy compared to the original SAPT Silver standard and comparable to SAPT Gold standard in some cases with much less computational cost.

Differing Effects of Nonlinearity around the Proton Acceptor on CH··O and NH··O H-Bond Strength within Proteins

The effects of deviations from nonlinearity around the carbonyl proton acceptor of an amide group are assessed by DFT quantum chemical calculations for both CH··O and NH··O H-bonds. The proton donors are the imidazole functional group of His and the indole of Trp, which are paired respectively with N-methylacetamide and acetamide. The displacement of either CH or NH group toward the carbonyl O sp2 lone pairs stabilizes the system and strengthens the H-bond. But the two donor groups differ in their response to a shift out of the amide plane. While the NH··O H-bond is weakened by this displacement, a substantial strengthening is observed when the CH donor is moved out of this plane, in one direction versus the other. This pattern is explained on the basis of simple Coulombic considerations.

Direct α-Hydroxy Acid Loading onto a Bacterial Thiotemplate Assembly Line via a Multienzyme Gateway

Various nonribosomal peptide synthetases (NRPSs) create structural and functional diversity by incorporating α-hydroxy acids into peptide backbones. Trigonic acid, an unusual cyclopropanol-substituted hydroxy acid, is the source of the molecular warhead of malleicyprol, a critical virulence factor of human and animal pathogens of the Burkholderia pseudomallei (BP) group. The process of selecting and loading this building block remained enigmatic as the NRPS module designated for this task is incomplete. Using a combination of bioinformatics, mutational analyses, targeted metabolomics, and in vitro biochemical assays, we show that two trans-acting enzymes are required to load this central building block onto the modular assembly line. An adenylation-thiolation didomain enzyme (BurJ) activates trigonic acid, followed by the translocation of the enzyme-bound α-hydroxy acid thioester by an FkbH-like protein with a mutated phosphatase domain (BurH). This specialized gateway is the first reported direct loading of an α-hydroxy acid onto a bona fide NRPS module in bacteria and expands the synthetic biology toolbox for the site-specific incorporation of non-canonical building blocks. Moreover, insight into the biochemical basis of virulence factor biosynthesis can provide a foundation for developing enzyme inhibitors as anti-virulence therapeutics against BP pathogen infections.

Combining Flexible Region Design and Automatic Design to Enhance the Thermal Stability and Catalytic Efficiency of Leucine Dehydrogenase

Leucine dehydrogenase (LeuDH, EC 1.4.1.9) can reversibly catalyze the oxidative deamination of l-leucine and some other specific α-amino acids to form the corresponding α-ketoacids. This reaction has great significance in the field of food additives and the pharmaceutical industry. The LeuDH from Exiguobacterium sibiricum (EsLeuDH) has high catalytic efficiency but limited thermal stability, hindering its widespread industrial application. In this study, a mutant N5F/I12L/A352Y of EsLeuDH (referred to as M2) was developed with enhanced thermal stability and catalytic activity through rational modification. The M2 mutant exhibits a half-life at 60 °C (t1/2(60 °C)) of 975.7 min and a specific activity of 69.6 U mg-1, which is 5.4 and 2.1 times higher than those of EsLeuDH, respectively. This research may facilitate the utilization of EsLeuDH at elevated temperatures, enhancing its potential for industrial applications. The findings offer a practical and efficient approach for optimizing LeuDH and other industrial enzymes.

The First Genetic Characterization of the SPRN Gene in Pekin Ducks (Anas platyrhynchos domesticus)

Prion diseases are fatal neurodegenerative disorders characterized by an accumulation of misfolded prion protein (PrPSc) in brain tissues. The shadow of prion protein (Sho) encoded by the shadow of prion protein gene (SPRN) is involved in prion disease progress. The interaction between Sho and PrP accelerates the PrPSc conversion rate while the SPRN gene polymorphisms have been associated with prion disease susceptibility in several species. Until now, the SPRN gene has not been investigated in ducks. We identified the duck SPRN gene sequence and investigated the genetic polymorphisms of 184 Pekin ducks. We compared the duck SPRN nucleotide sequence and the duck Sho protein amino acid sequence with those of several other species. Finally, we predicted the duck Sho protein structure and the effects of non-synonymous single nucleotide polymorphisms (SNPs) using computational programs. We were the first to report the Pekin duck SPRN gene sequence. The duck Sho protein sequence showed 100% identity compared with the chicken Sho protein sequence. We found 27 novel SNPs in the duck SPRN gene. Four amino acid substitutions were predicted to affect the hydrogen bond distribution in the duck Sho protein structure. Although MutPred2 and SNPs&GO predicted that all non-synonymous polymorphisms were neutral or benign, SIFT predicted that four variants, A22T, G49D, A68T, and M105I, were deleterious. To the best of our knowledge, this is the first report about the genetic and structural characteristics of the duck SPRN gene.

Thioimidates provide general access to thioamide, amidine, and imidazolone peptide-bond isosteres

Thioamides, amidines, and heterocycles are three classes of modifications that can act as peptide-bond isosteres to alter the peptide backbone. Thioimidate protecting groups can address many of the problematic synthetic issues surrounding installation of these groups. Historically, amidines have received little attention in peptides due to limitations in methods to access them. The first robust and general procedure for the introduction of amidines into peptide backbones exploits the utility of thioimidate protecting groups as a means to side-step reactivity that ultimately renders existing methods unsuitable for the installation of amidines along the main-chain of peptides. Further, amidines formed on-resin can be reacted to form (4H)-imidazolone heteorcycles which have recently been shown to act as cis-amide isosteres. General methods for heterocyclic installation capable of geometrically restricting peptide conformation are also under-developed. This work is significant because it describes a generally applicable and divergent approach to access unexplored peptide designs and architectures.

Thioimidate Solutions to Thioamide Problems during Thionopeptide Deprotection

Thioamides have structural and chemical similarity to peptide bonds, offering valuable insights when probing peptide backbone interactions, but are prone to side reactions during solid-phase peptide synthesis (SPPS). Thioimidates have been demonstrated to be effective protecting groups for thioamides during peptide elongation. We further demonstrate how thioimidates can assist thioamides through the most yield-crippling step of thionopeptide deprotection, allowing for the first isolation of an important benchmark α-helical peptide that had previously eluded synthesis and isolation.

Membrane-Bound Redox Enzyme Cytochrome bd-I Promotes Carbon Monoxide-Resistant Escherichia coli Growth and Respiration

The terminal oxidases of bacterial aerobic respiratory chains are redox-active electrogenic enzymes that catalyze the four-electron reduction of O2 to 2H2O taking out electrons from quinol or cytochrome c. Living bacteria often deal with carbon monoxide (CO) which can act as both a signaling molecule and a poison. Bacterial terminal oxidases contain hemes; therefore, they are potential targets for CO. However, our knowledge of this issue is limited and contradictory. Here, we investigated the effect of CO on the cell growth and aerobic respiration of three different Escherichia coli mutants, each expressing only one terminal quinol oxidase: cytochrome bd-I, cytochrome bd-II, or cytochrome bo3. We found that following the addition of CO to bd-I-only cells, a minimal effect on growth was observed, whereas the growth of both bd-II-only and bo3-only strains was severely impaired. Consistently, the degree of resistance of aerobic respiration of bd-I-only cells to CO is high, as opposed to high CO sensitivity displayed by bd-II-only and bo3-only cells consuming O2. Such a difference between the oxidases in sensitivity to CO was also observed with isolated membranes of the mutants. Accordingly, O2 consumption of wild-type cells showed relatively low CO sensitivity under conditions favoring the expression of a bd-type oxidase.

n → π* Interaction Enabling Transient Inversion of Chirality

This study reports the observation and characterization of two isomers of the acrolein dimer by using high-resolution rotational spectroscopy in pulsed jets. The first isomer is stabilized by two hydrogen bonds, adopting a planar configuration, and is energetically favored over the second isomer, which exhibits a dominant n → π* interaction in a nearly orthogonal arrangement. Surprisingly, the n → π* interaction was revealed to enable a concerted tunneling motion of two moieties along the carbonyl group. This motion leads to the inversion of transient chirality associated with the exchange of donor-acceptor roles, as revealed by the spectral feature of quadruplets. Inversion of transient chirality is a fundamental phenomenon in quantum mechanics and commonly observed for only inversional motions of protons. It is the first discovery, to the best of our knowledge, that such heavy moieties can also undergo chirality inversion.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally Determined Enthalpic and Entropic Contributions to Interactions of Amide sp2O, N, C and sp3C Unified Atoms with Naphthalene sp2C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp2O with amide sp2N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp2C (lone pair π, n-π*) are particularly favorable. Sp3C-sp3C (hydrophobic), sp3C-sp2C (hydrophobic, CH-π), sp2C-sp2C (hydrophobic, π-π), and sp3C-sp2N interactions are favorable, sp2C-sp2N interactions are neutral, while sp2O-sp2O and sp2N-sp2N self-interactions and sp2O-sp3C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp2O, sp2N, sp2C, and sp3C unified atoms with aromatic sp2C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp2O-aromatic sp2C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp3C- and sp2C-aromatic sp2C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp2 atoms in protein processes.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally-determined Enthalpic and Entropic Contributions to Interactions of Amide sp 2 O, N, C and sp 3 C Unified Atoms with Naphthalene sp 2 C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics and other groups occur in protein self-assembly processes including folding, oligomerization and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C=O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair-π, n-π * ) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π) and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of fourteen amides on naphthalene solubility at 10, 25 and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g. lone pair-π) while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

Table of Contents Graphic

Synergistic n → π* and nN → π*Ar interactions in C-terminal modified prolines: effect on Xaa-Pro cis/trans equilibrium

Carbonyl-carbonyl (CO⋯CO) n → π* interaction often coexists with a hydrogen bond (HB) or another n → π* interaction. Although the interplay between HB and n → π* interaction was previously studied, there is no systematic investigation that shows a synergistic relationship of n → π* with another noncovalent interaction. Herein, we have studied a set of proline-diacylhydrazine (Pro-DAH) molecules and observed that increase in the strength of the n_N → π*_Ar interaction on their DAH side strengthened the n → π* interaction on the Pro side, which was experimentally determined by measuring the K_trans/cis of the Xaa-Pro amide bond. Overall, we describe a simple C-terminal modification strategy to stabilize the trans-Pro geometry that could be useful to stabilize PPII helices and collagen triple helices that require Pro to adopt the trans amide geometry.

Effect of Substituents on the Intramolecular n→π* Interaction in 3-[2-(Dimethylamino) phenyl] propanal: A Computational Study

n→π* non-covalent interaction (NCI) and hydrogen bond have similarity in terms of delocalization of the electron density between the two orbitals involved in the interaction. Hydrogen bond (X-H···Y) involves delocalization of the lone pair electrons (n) on the Y atom into the σ* orbital of the X-H bond. In contrast, the n→π* interaction deals with delocalizing the lone pair electrons (n) on the N, O, or S atom into the π* orbital of a C═O group or aromatic ring. Herein, we have shown a resemblance of this weak n→π* interaction with the relatively stronger hydrogen bond in terms of folding the side chains in flexible molecules. This work reports the study of folding of the flexible side-chain in 3-[2-(dimethylamino) phenyl] propanal (DMAPhP) through a N···C═O n→π* interaction using various computational approaches such as NBO, QTAIM, and NCI analyses. The folding of the molecule by the n→π* interaction observed in this study is found to be similar to that present in the secondary structures of peptides or proteins through hydrogen bonding interactions. Interestingly, the stabilization of the global minimum conformer of DMAPhP by the n→π* interaction demonstrates the importance of this NCI in providing conformational preferences in molecular systems. Another important finding of this study is that the theoretical redshift obtained in the C═O stretching frequency of the most stable conformer of the DMAPhP is contributed mostly by the n→π* interaction as the C═O group is not involved in hyperconjugation with any neighboring heteroatom, which is a common phenomenon in any ester or amide. We have also demonstrated here that the strength of the intramolecular n→π* interaction can be modulated by varying the electronic substituents at the para position of the donor group involved in the interaction.

Importance of Noncovalent Interactions Involving Sulfur Atoms in Thiopeptide Antibiotics─Glycothiohexide α and Nocathiacin I

Noncovalent interactions involving sulfur atoms play essential roles in protein structure and function by significantly contributing to protein stability, folding, and biological activity. Sulfur is a highly polarizable atom that can participate in many types of noncovalent interactions including hydrogen bonding, sulfur-π interactions, and S-lone pair interactions, but the impact of these sulfur-based interactions on molecular recognition and drug design is still often underappreciated. Here, we examine, using quantum chemical calculations, the roles of sulfur-based noncovalent interactions in complex naturally occurring molecules representative of thiopeptide antibiotics: glycothiohexide α and its close structural analogue nocathiacin I. While donor-acceptor orbital interactions make only very small contributions, electrostatic and dispersion contributions are predicted to be significant in many cases. In pursuit of understanding the magnitudes and nature of these noncovalent interactions, we made potential structural modifications that could significantly expand the chemical space of effective thiopeptide antibiotics.

Unexpected luminescence of non-conjugated biomass-based polymers: new approach in photothermal imaging

Population growth, depletion of world resources and persistent toxic chemical production underline the need to seek new smart materials from inexpensive, biodegradable, and renewable feedstocks. Hence, "metal-free" ring-opening copolymerization to convert biomass carvone-based monomers into non-conventional luminescent biopolymers is considered a sustainable approach to achieve these goals. The non-conventional emission was studied in terms of steady-state and time-resolved spectroscopy in order to unravel the structure-properties for different carvone-based copolymers. The results highlighted the importance of the final copolymer folding structure as well as its environment in luminescent behavior (cluster-triggered emission). In all cases, their luminescent behavior is sensitive to small temperature fluctuations (where the minimum detected temperature is T_m ∼ 2 °C and relative sensitivity is S_r ∼ 6% °C) even at the microscopic scale, which endows these materials a great potential as thermosensitive smart polymers for photothermal imaging.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

Evaluating the impact of X-ray damage on conformational heterogeneity in room-temperature (277 K) and cryo-cooled protein crystals

Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02–1.52 Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277 K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to ∼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9–1.16 Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.

Structure–Activity Relationship Analysis of Flavonoids and Its Inhibitory Activity Against BACE1 Enzyme Toward a Better Therapy for Alzheimer’s Disease

Drug development in Alzheimer’s disease (AD) suffers from a high attrition rate. In 2021, 117 agents tested in phases I and II and 36 agents tested in phase III were discontinued. Natural product compounds may be good lead compounds for AD as they contain functional groups that are important for binding against key AD targets such as β-secretase enzyme (BACE1). Hence, in this study, 64 flavonoids collected from rigorous literature search and screening that have been tested from 2010 to 2022 against BACE1, which interferes in the formation of amyloid plaque, were analyzed. The 64 unique flavonoids can be further classified into five core fragments. The flavonoids were subjected to clustering analysis based on its structure, and each representative of the clusters was subjected to molecular docking. There were 12 clusters formed, where only 1 cluster contained compounds from two different core fragments. Several observations can be made where 1) flavanones with sugar moieties showed higher inhibitory activity compared to flavanones without sugar moieties. The number of sugar moieties and position of glycosidic linkage may also affect the inhibitory activity. 2) Non-piperazine-substituted chalcones when substituted with functional groups with decreasing electronegativity at the para position of both rings result in a decrease in inhibitory activity. Molecular docking indicates that ring A is involved in hydrogen bond, whereas ring B is involved in van der Waals interaction with BACE1. 3) Hydrogen bond is an important interaction with the catalytic sites of BACE1, which are Asp32 and Asp228. As flavonoids contain favorable structures and properties, this makes them an interesting lead compound for BACE1. However, to date, no flavonoids have made it through clinical trials. Hence, these findings may aid in the design of highly potent and specific BACE1 inhibitors, which could delay the progression of AD.

Solvation stabilizes intercarbonyl n→π* interactions and polyproline II helix

n→π* interactions between consecutive carbonyls stabilize the α-helix and polyproline II helix (PPII) conformations in proteins. n→π* interactions have been suggested to provide significant conformational biases to the disordered states of proteins. To understand the roles of solvation on the strength of n→π* interactions, computational investigations were conducted on a model n→π* interaction, the twisted-parallel-offset formaldehyde dimer, as a function of explicit solvation of the donor and acceptor carbonyls, using water and HF. In addition, the effects of urea, thiourea, guanidinium, and monovalent cations on n→π* interaction strength were examined. Solvation of the acceptor carbonyl significantly strengthens the n→π* interaction, while solvation of the donor carbonyl only modestly weakens the n→π* interaction. The n→π* interaction strength was maximized with two solvent molecules on the acceptor carbonyl. Urea stabilized the n→π* interaction via simultaneous engagement of both oxygen lone pairs on the acceptor carbonyl. Solvent effects were further investigated in the model peptides Ac-Pro-NMe₂, Ac-Ala-NMe₂, and Ac-Pro₂-NMe₂. Solvent effects in peptides were similar to those in the formaldehyde dimer, with solvation of the acceptor carbonyl increasing n→π* interaction strength and resulting in more compact conformations, in both the proline endo and exo ring puckers, as well as a reduction in the energy difference between these ring puckers. Carbonyl solvation leads to an energetic preference for PPII over both the α-helix and β/extended conformations, consistent with experimental data that protic solvents and protein denaturants both promote PPII. Solvation of the acceptor carbonyl weakens the intraresidue C5 hydrogen bond that stabilizes the β conformation.

SM-COLSARSPROT: Highly Immunogenic Supramutational Synthetic Peptides Covering the World's Population

Fifty ~20-amino acid (aa)-long peptides were selected from functionally relevant SARS-CoV-2 S, M, and E proteins for trial B-21 and another 53 common ones, plus some new ones derived from the virus' main genetic variants for complementary trial C-21. Peptide selection was based on tremendous SARS-CoV-2 genetic variability for analysing them concerning vast human immunogenetic polymorphism for developing the first supramutational, Colombian SARS-protection (SM-COLSARSPROT), peptide mixture. Specific physicochemical rules were followed, i.e., aa predilection for polyproline type II left-handed (PPIIL) formation, replacing β-branched, aromatic aa, short-chain backbone H-bond-forming residues, π-π interactions (n→π* and π-CH), aa interaction with π systems, and molecular fragments able to interact with them, disrupting PPIIL propensity formation. All these modified structures had PPIIL formation propensity to enable target peptide interaction with human leukocyte antigen-DRβ1* (HLA-DRβ1*) molecules to mediate antigen presentation and induce an appropriate immune response. Such modified peptides were designed for human use; however, they induced high antibody titres against S, M, and E parental mutant peptides and neutralising antibodies when suitably modified and chemically synthesised for immunising 61 major histocompatibility complex class II (MHCII) DNA genotyped Aotus monkeys (matched with their corresponding HLA-DRβ1* molecules), predicted to cover 77.5% to 83.1% of the world's population. Such chemically synthesised peptide mixture represents an extremely pure, stable, reliable, and cheap vaccine for COVID-19 pandemic control, providing a new approach for a logical, rational, and soundly established methodology for other vaccine development.

Sp2- and sp3-C⋯O tetrel bonds in the 3-oxetanone homodimer

The structures and non-covalent interactions at play in the 3-oxetanone homodimer have been investigated using a pulsed jet Fourier transform microwave spectrometer supplemented with quantum chemical calculations. Two isomers were identified in the pulsed jet. With the analyses of non-covalent intermolecular interactions including the quantum theory of atoms, Johnson's non-covalent interactions and natural bond orbital, the observed global minimum is stabilized by a combination of one sp²-C⋯O tetrel bond and a network of multiple C-H⋯O weak hydrogen bonds. The second isomer is characterized by carbonyl-carbonyl interactions, with the formation of one sp²- and one sp³-C⋯O tetrel bond. The conformational population of the two observed isomers in the supersonic expansion was estimated to be N_CE1/N_CC1 ≈ 7/5.

Ethyl 5-Formyl-1-(pyridin-3-yl)-1H-1,2,3-triazole-4-carboxylate: Synthesis, Crystal Structure, Hirshfeld Surface Analysis, and DFT Calculation

For the first time, 5-formyl-1-(pyridin-3-yl)-1H-1,2,3-triazole-4-carboxylate was synthesized via a two-step scheme. The molecular structure of the compound was determined by a single-crystal X-ray diffraction analysis. The Hirshfeld surface analysis was used to study various intermolecular interactions. The crystalline structure is marked by the presence of three types of π-interactions (n→π*, lp···π, and π···π) between the -C(H)=O group and triazole rings. The compound is a versatile polyfunctional building block for construction of annulated 1,2,3-triazoles.

Analysis of the Orbital and Electrostatic Contributions to the Lone Pair-Aromatic Interaction Using Molecular Rotors

The attractive interaction between carbonyl oxygens and the π-face of aromatic surfaces was studied using N-phenylimide molecular rotors. The C═O···Ar interactions could stabilize the transition states but were half the strength of comparable C═O···C═O interactions. The C═O···Ar interaction had a significant electrostatic component but only a small orbital delocalization component.

Tuning Supramolecular Polymer Assembly through Stereoelectronic Interactions

The supramolecular polymerization of 2,11-dithia[3.3]paracyclophanes through self-complementary intermolecular and transannular amide hydrogen bonding is presented. An n → π* interaction between the amide hydrogen bonding units and the central bridging atom results from the single-point exchange of a carbon atom for a sulfur atom. This orbital donor-acceptor interaction can be strengthened by oxidizing the sulfide to a sulfone which acts to shorten the donor···acceptor distance and increase orbital overlap. Experimental signatures of the increased n → π* interaction include larger isodesmic polymerization elongation constants in solution, changes in characteristic bond stretching frequencies, and geometric/structural changes evaluated by X-ray crystallography. The experimental data are supported by extensive computational investigations of both assembling and nonassembling 2,11-dithia[3.3]paracyclophanes as well as a rationally designed model system to confirm the role of stereoelectronic effects on supramolecular polymer assembly.

Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation

Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π^∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.

Dynamic Community Composition Unravels Allosteric Communication in PDZ3

The third domain of PSD-95 (PDZ3) is a model for investigating allosteric communication in protein and ligand interactions. While motifs contributing to its binding specificity have been scrutinized, a conformational dynamical basis is yet to be established. Despite the miniscule structural changes due to point mutants, the observed significant binding affinity differences have previously been assessed with a focus on two α-helices located at the binding groove (α₂) and the C-terminus (α₃). Here, we employ a new computational approach to develop a generalized view on the molecular basis of PDZ3 binding selectivity and interaction communication for a set of point mutants of the protein (G330T, H372A, G330T-H372A) and its ligand (CRIPT, named L₁, and its T-2F variant, L₂) along with the wild type (WT). To analyze the dynamical aspects hidden in the conformations that are produced by molecular dynamics simulations, we utilize variations in community composition calculated based on the betweenness centrality measure from graph theory. We find that the highly charged N-terminus, which is located far from the ligand, has the propensity to share the same community with the ligand in the biologically functional complexes, indicating a distal segment might mediate the binding dynamics. N- and C-termini of PDZ3 share communities, and α₃ acts as a hub for the whole protein by sustaining the communication with all structural segments, albeit being a trait not unique to the functional complexes. Moreover, α₂ which lines the binding cavity frequently parts communities with the ligand and is not a controller of the binding but is rather a slave to the overall dynamics coordinated by the N-terminus. Thus, ligand binding fate in PDZ3 is traced to the population of community compositions extracted from dynamics despite the lack of significant conformational changes.

Spectroscopic evidence of n → π* interactions involving carbonyl groups

n → π* has emerged as an important noncovalent interaction that can affect the conformations of both small- and macromolecules including peptides and proteins. Carbonyl-carbonyl (COCO) n → π* interactions involving CO groups are well studied. Recent studies have shown that the COCO n → π* interactions are the most abundant secondary interactions in proteins with a frequency of 33 interactions per 100 residues and, among the various secondary interactions, n → π* interactions are expected to provide the highest enthalpic contributions to the conformational stability of globular proteins. However, n → π* interactions are relatively weak and provide an average stabilization of about 0.25 kcal mol-1 per interaction in proteins. The strongest n → π* interaction could be as strong as a moderate hydrogen bond. Therefore, it is challenging to detect and quantify these weak interactions, especially in solution in the presence of perturbation from other intermolecular interactions. Accordingly, spectroscopic investigations that can provide direct evidence of n → π* interaction are limited, and the majority of papers published in this area have relied on X-ray crystallography and/or theoretical calculations to establish the presence of this interaction. The aim of this perspective is to discuss the studies where a spectroscopic signature in support of n → π* interaction was observed. As the "n → π* interaction" is a relatively new terminology, there remains the possibility of there being earlier studies where spectroscopic evidence for n → π* interactions was obtained but it was not discussed in light of the n → π* terminology. We noticed several such studies and, as can be expected, these studies were often overlooked in the discussion of n → π* interactions in the recent literature. In this perspective, we have also discussed these studies and provided computational support for the presence of n → π* interaction.

Formation of n → π+ interaction facilitating dissociative electron transfer in isolated tyrosine-containing molecular peptide radical cations

Long-range electron transfer in proteins can be rationalized as a sequential short-distance electron-hopping processes via amino acid residues having low ionization energy as relay stations. Tyrosine residues can serve as such redox-active intermediates through one-electron oxidation to form a π-radical cation at its phenol side chain. An electron transfer from a vicinal functional group to this π-electron hole completes an elementary step of charge migration. However, transient oxidized/reduced intermediates formed at those relay stations during electron transfer processes have not been observed. In this study, formation of analog reactive intermediates via electron donor-acceptor coupling is observed by using IRMPD action spectroscopy. An elementary charge migration at the molecular level in model tyrosine-containing peptide radical cations [M]˙+ in the gas phase is revealed with its unusual Cα-Cβ bond cleavage at the side chain of the N-terminal residue. This reaction is induced by the radical character of the N-terminal amino group (-NH2˙+) resulting from an n → π+ interaction between the nonbonding electron pair of NH2 (n) and the π-electron hole at the Tyr side chain (π+). The formation of -NH2˙+ is supported by the IRMPD spectrum showing a characteristic NH2 scissor vibration coupled with Tyr side-chain stretches at 1577 cm-1. This n → π+ interaction facilitates a dissociative electron transfer with NH2 as the relay station. The occurrence of this side-chain cleavage may be an indicator of the formation of reactive conformers featuring the n → π+ interaction.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Fluorescence in "Nonfluorescent" Polymers

Recently, a great deal of research has been started on generating fairly strong photoluminescence from organic molecules without having any conjugated π-system or fluorophore. Discrete chromophores or auxochromophores termed as "subfluorophores" may undergo "space conjugation" via co-operative intramolecular conformation followed by intermolecular aggregation to generate fluorescence or sometimes phosphorescence emission. Polymeric materials are important in this regard as nonconjugated polymers self-assemble/aggregate in a moderately concentrated solution and also in the solid state, producing membranes, films, and so forth with good physical and mechanical properties. Therefore, promoting fluorescence in these commodity polymers is very much useful for sensing, organic light emitting diodes (OLED), and biological applications. In this perspective, we have discussed the aggregation-induced emission from four different types of architectures, for example, (i) dendrimers or hyperbranched polymers, (ii) entrapped polymeric micellar self-assembly, (iii) cluster formation, and (iv) stretching-induced aggregation, begining with the genesis of fluorescence from aggregation of propeller-shaped small organic molecules. The mechanism of induced fluorescence of polymers with subfluorophoric groups is also discussed from the theoretical calculations of the energy bands in the aggregated state. Also, an attempt has been made to highlight some useful applications in the sensing of surfactants, bacteria, cell imaging, drug delivery, gene delivery, OLED, and so forth.

Experimentally determined strengths of favorable and unfavorable interactions of amide atoms involved in protein self-assembly in water

Folding and other protein self-assembly processes are driven by favorable interactions between O, N, and C unified atoms of the polypeptide backbone and side chains. These processes are perturbed by solutes that interact with these atoms differently than water does. Amide NH···O=C hydrogen bonding and various π-system interactions have been better characterized structurally or by simulations than experimentally in water, and unfavorable interactions are relatively uncharacterized. To address this situation, we previously quantified interactions of alkyl ureas with amide and aromatic compounds, relative to interactions with water. Analysis yielded strengths of interaction of each alkylurea with unit areas of different hybridization states of unified O, N, and C atoms of amide and aromatic compounds. Here, by osmometry, we quantify interactions of 10 pairs of amides selected to complete this dataset. An analysis yields intrinsic strengths of six favorable and four unfavorable atom-atom interactions, expressed per unit area of each atom and relative to interactions with water. The most favorable interactions are sp2O-sp2C (lone pair-π, presumably n-π*), sp2C-sp2C (π-π and/or hydrophobic), sp2O-sp2N (hydrogen bonding) and sp3C-sp2C (CH-π and/or hydrophobic). Interactions of sp3C with itself (hydrophobic) and with sp2N are modestly favorable, while sp2N interactions with sp2N and with amide/aromatic sp2C are modestly unfavorable. Amide sp2O-sp2O interactions and sp2O-sp3C interactions are more unfavorable, indicating the preference of amide sp2O to interact with water. These intrinsic interaction strengths are used to predict interactions of amides with proteins and chemical effects of amides (including urea, N-ethylpyrrolidone [NEP], and polyvinylpyrrolidone [PVP]) on protein stability.

n→π* Interactions Modulate the Disulfide Reduction Potential of Epidithiodiketopiperazines

Epithiodiketopiperazines (ETPs) are a structurally complex class of fungal natural products with potent anticancer activity. In ETPs, the diketopiperazine ring is spanned by a disulfide bond that is constrained in a high-energy eclipsed conformation. We employed computational, synthetic, and spectroscopic methods to investigate the physicochemical attributes of this atypical disulfide bond. We find that the disulfide bond is stabilized by two n→π* interactions, each with large energies (3-5 kcal/mol). The n→π* interactions in ETPs make disulfide reduction much more difficult, endowing stability in physiological environments in a manner that could impact their biological activity. These data reveal a previously unappreciated means to stabilize a disulfide bond and highlight the utility of the n→π* interaction in molecular design.

Increasing protein stability by engineering the n → π* interaction at the β-turn

Abundant n → π* interactions between adjacent backbone carbonyl groups, identified by statistical analysis of protein structures, are predicted to play an important role in dictating the structure of proteins. However, experimentally testing the prediction in proteins has been challenging due to the weak nature of this interaction. By amplifying the strength of the n → π* interaction via amino acid substitution and thioamide incorporation at a solvent exposed β-turn within the GB1 proteins and Pin 1 WW domain, we demonstrate that an n → π* interaction increases the structural stability of proteins by restricting the ϕ torsion angle. Our results also suggest that amino acid side-chain identity and its rotameric conformation play an important and decisive role in dictating the strength of an n → π* interaction.

Amino acid residues adopt a right-handed α-helical conformation with increasing strength of the n → π* interaction. We also demonstrate a direct consequence of n → π* interactions on enhancing the structural stability of proteins.

Intra- and intermolecular interactions in a series of chlorido-tricarbonyl-diazabutadienerhenium(I) complexes: structural and theoretical studies

A series of new chlorido-tricarbonylrhenium(I) complexes bearing alkyl-substituted diazabutadiene (DAB) ligands, namely N,N'-bis(2,4-dimethylbenzene)-1,4-diazabutadiene (L1), N,N'-bis(2,4-dimethylbenzene)-2,3-dimethyl-1,4-diazabutadiene (L2), N,N'-bis(2,4,6-trimethylbenzene)-2,3-dimethyl-1,4-diazabutadiene (L3) and N,N'-bis(2,6-diisopropylbenzene)-1,4-diazabutadiene (L4), were synthesized and investigated. The crystal structures have been fully characterized by X-ray diffraction and spectroscopic methods. Density functional theory, natural bond orbital and non-covalent interaction index methods have been used to study the optimized geometry in the gas phase and intra- and intermolecular interactions in the complexes, respectively. The most important studied interactions in these metal carbonyl complexes are n→π*, n→σ* and π→π*. Among complexes 1-4, only 2 shows interesting intermolecular n→π* interactions due to lp(C[triple-bond]O)...π* and lp(Cl)...π* (lp = lone pair) contacts.

Cation-π Interactions in Graphene-Containing Systems for Water Treatment and Beyond

Cation-π interactions are common in nature, especially in organisms. Their profound influences in chemistry, physics, and biology have been continuously investigated since they were discovered in 1981. However, the importance of cation-π interactions in materials science, regarding carbonaceous nanomaterials, has just been realized. The interplay between cations and delocalized polarizable π electrons of graphene would bring about significant changes to the intrinsic characteristics of graphene and greatly affect the device performance based on graphene and its derivatives. Here, the cation-π interactions in graphene containing systems for water treatment applications (e.g., separation membranes, adsorbents) are highlighted. The cross-linking effects caused by cation-π interactions contribute to membrane stability and selectivity and enhanced adsorption. Their roles in dominating the performance of graphene-based structures for other specific applications are also discussed. Relevant theoretical modeling and calculations are summarized to offer an in-depth understanding of the underlying mechanisms which can help in designing more functional materials and structures. Perspectives on the potential directions that deserve effort are also presented.

Structural, Non-Covalent Interaction, and Natural Bond Orbital Studies on Bromido-Tricarbonyl Rhenium(I) Complexes Bearing Alkyl-Substituted 1,4-Diazabutadiene (DAB) Ligands

The synthesis, characterization, structural and computational studies of Re(I) tricarbonyl bromo complexes bearing alkyl-substituted 1,4-diazabutadiene ligands, [Re(CO)3(1,4-DAB)Br], where 1,4-DAB = N,N-bis(2,4-dimethylbenzene)-1,4-diazabutadiene, 2,4-Me2DAB (1); N,N-bis(2,4-dimethylbenzene)-2,3-dimethyl-1,4-diazabutadiene, 2,4-Me2DABMe (2); N,N-bis(2,4,6-trimethylbenzene)-1,4-diazabutadiene, 2,4,6-Me3DAB (3); and N,N-bis(2,6-diisopropylbenzene)-1,4-diazabutadiene, 2,6-ipr2DAB (4) are reported. The complexes were characterized by different spectroscopic methods such as FT-IR, 1H-NMR, 13C-NMR, and elemental analyses and their solid-state structures were confirmed by X-ray diffraction. In each complex, the Re(I) centre shows a distorted octahedral shape with a facial geometry of carbonyl groups. The gas phase geometry of the complexes was identified by density functional theory. Interesting intermolecular n…π* interactions of complexes 1 and 3 were investigated by non-covalent interaction index (NCI), and natural bond orbital (NBO) analyses. The intramolecular n…σ*, σ…π*, π…σ* interactions were also studied in complexes 3 and 4.

Supramolecular assembly through intermolecular n → π* interactions through a coordinated perrhenate formed via superoxidation of Re(i) to Re(vii) in the formation of substituted Re(CO)3 complexes bearing Diimine ligands

Two Re(i) tricarbonyl complexes bearing Ph₄TAP and dafone ligands and having a coordinated perrhenate group obtained via in situ superoxidation of Re(i) to Re(vii) have been characterized.

Conformational control of N-methyl-N,N′-diacylhydrazines by noncovalent carbon bonding in solution

Unusual N(amide)⋯C–X noncovalent carbon bonding interactions stabilize the trans–cis (t–c) amide bond rotamers of N-methyl-N,N′-diacylhydrazines over the expected trans–trans (t–t) rotamers in solution.

Non-covalent sulfoxide⋯(nitrosyl group) interactions involving coordinated nitrosyl in a Ru(ii) nitrosyl complex with an α-diimine ligand: structural and computational studies

Inspection of the X-ray structure of the newly prepared Ru–nitrosyl complex bearing an α-diimine ligand revealed for the first time the π-hole interaction involving the coordinated nitrosyl group with DMSO.

Measuring multiple 17O–13C J-couplings in naphthalaldehydic acid: a combined solid state NMR and density functional theory approach

A combined multinuclear solid-state NMR and a density functional theory computational approach, with SIMPSON simulations, is evaluated to determine the four heteronuclear ¹J(¹³C,¹⁷O) couplings in naphthalaldehydic acid.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.