An electronic effect on protein structure

Lone Pair-π Interactions in Organic Reactions

Noncovalent interactions between a lone pair of electrons and π systems can be categorized into two types based on the nature of π systems. Lone pair-π(C═O) interactions with π systems of unsaturated, polarized bonds are primarily attributed to orbital interactions, whereas lone pair-π(Ar) interactions with π systems of aromatic functional groups result from electrostatic attractions (for electron-deficient aryls) or dispersion attractions and Pauli repulsions (for electron-rich/neutral aryls). Unlike well-established noncovalent interactions, lone pair-π interactions have been comparatively underappreciated or less used to influence reaction outcomes. This review emphasizes experimental and computational studies aimed at integrating lone pair-π interactions into the design of catalytic systems and utilizing these interactions to regulate the reactivity and selectivity of chemical transformations. The role of lone pair-π interactions is highlighted in the stabilization or destabilization of transition states and ground-state binding. Examples influenced by lone pair-π interactions with both unsaturated, polarized bonds and aromatic rings as π systems are included. At variance with previous reviews, the present review is not structured according to the physical origin of particular classes of lone pair-π interactions but is divided into chapters according to ways in which lone pair-π interactions affect kinetics and/or selectivity of reactions.

Kinetic and Mechanistic Studies of Native Chemical Ligation with Phenyl α-Selenoester Peptides

Native chemical ligation (NCL) ligates two unprotected peptides in an aqueous buffer. One of the fragments features a C-terminal α-thioester functional group, and the second bears an N-terminal cysteine. The reaction mechanism depicts two steps: an intermolecular thiol-thioester exchange resulting in a transient thioester, followed by an intramolecular S-to-N acyl shift to yield the final native peptide bond. Although this mechanism is well established, the direct observation of the transient thioester has been elusive because the fast intramolecular rearrangement prevents its accumulation. Here, the use of α-selenoester peptides allows a faster first reaction and an early buildup of the intermediate, enabling its quantification and the kinetic monitoring of the first and second steps. The results show a correlation between the steric hindrance in the α-thioester residue and the rearrangement rate. In bulky residues, the S-to-N acyl shift has a significant contribution to the overall reaction rate. This is particularly notable for valine and likely for other similar β-branched amino acids.

Electronic Control of Polyproline II Helix Stability via the Identity of Acyl Capping Groups: the Pivaloyl Group Particularly Promotes PPII

The type II polyproline helix (PPII) is a fundamental secondary structure of proteins, important in globular proteins, in intrinsically disordered proteins, and at protein-protein interfaces. PPII is stabilized in part by n→π* interactions between consecutive carbonyls, via electron delocalization between an electron-donor carbonyl lone pair (n) and an electron-acceptor carbonyl (π*) on the subsequent residue. We previously demonstrated that changes to the electronic properties of the acyl donor can predictably modulate the strength of n→π* interactions, with data from model compounds, in solution in chloroform, in the solid state, and computationally. Herein, we examined whether the electronic properties of acyl capping groups could modulate the stability of PPII in peptides in water. In X-PPGY-NH2 peptides (X=10 acyl capping groups), the effect of acyl group identity on PPII was quantified by circular dichroism and NMR spectroscopy. Electron-rich acyl groups promoted PPII relative to the standard acetyl (Ac-) group, with the pivaloyl and iso-butyryl groups most significantly increasing PPII. In contrast, acyl derivatives with electron-withdrawing substituents and the formyl group relatively disfavored PPII. Similar results, though lesser in magnitude, were also observed in X-APPGY-NH2 peptides, indicating that the capping group can impact PPII conformation at both proline and non-proline residues. The pivaloyl group was particularly favorable in promoting PPII. The effects of acyl capping groups were further analyzed in X-DfpPGY-NH2 and X-ADfpPGY-NH2 peptides, Dfp=4,4-difluoroproline. Data on these peptides indicated that acyl groups induced order Piv- > Ac- > For-. These results suggest that greater consideration should be given to the identity of acyl capping groups in inducing structure in peptides.

Proline Analogues

Within the canonical repertoire of the amino acid involved in protein biogenesis, proline plays a unique role as an amino acid presenting a modified backbone rather than a side-chain. Chemical structures that mimic proline but introduce changes into its specific molecular features are defined as proline analogues. This review article summarizes the existing chemical, physicochemical, and biochemical knowledge about this peculiar family of structures. We group proline analogues from the following compounds: substituted prolines, unsaturated and fused structures, ring size homologues, heterocyclic, e.g., pseudoproline, and bridged proline-resembling structures. We overview (1) the occurrence of proline analogues in nature and their chemical synthesis, (2) physicochemical properties including ring conformation and cis/trans amide isomerization, (3) use in commercial drugs such as nirmatrelvir recently approved against COVID-19, (4) peptide and protein synthesis involving proline analogues, (5) specific opportunities created in peptide engineering, and (6) cases of protein engineering with the analogues. The review aims to provide a summary to anyone interested in using proline analogues in systems ranging from specific biochemical setups to complex biological systems.

N···C═O n → π* Interaction: Gas-Phase Electronic and Vibrational Spectroscopy Combined with Quantum Chemistry Calculations

Herein, we have used gas-phase electronic and vibrational spectroscopic techniques for the first time to study the N···C═O n → π* interaction in ethyl 2-(2-(dimethylamino) phenyl) acetate (NMe2-Ph-EA). We have measured the electronic spectra of NMe2-Ph-EA in the mass channels of its two distinct fragments of m/z = 15 and 192 using a resonant two-photon ionization technique as there was extensive photofragmentation of NMe2-Ph-EA. Identical electronic spectra obtained in the mass channels of both fragments confirm the dissociation of NMe2-Ph-EA in the ionic state, and hence, the electronic spectrum of the fragment represents that of NMe2-Ph-EA only. UV-UV hole-burning spectroscopy proved the presence of a single conformer of NMe2-Ph-EA in the experiment. Detailed quantum chemistry calculations reveal the existence of a N···C═O n → π* interaction in all six low-energy conformers of NMe2-Ph-EA. A comparison of the IR spectrum of NMe2-Ph-EA acquired from the gas-phase experiment with those obtained from theoretical calculations indicates that the experimentally observed conformer has a N···C═O n → π* interaction. The present finding might be further valuable in drug design and their recognition based on the N···C═O n → π* interaction.

Probing Non-Covalent Interactions through Molecular Balances: A REG-IQA Study

The interaction energies of two series of molecular balances (1-X with X = H, Me, OMe, NMe2 and 2-Y with Y = H, CN, NO2, OMe, NMe2) designed to probe carbonyl…carbonyl interactions were analysed at the B3LYP/6-311++G(d,p)-D3 level of theory using the energy partitioning method of Interacting Quantum Atoms/Fragments (IQA/IQF). The partitioned energies are analysed by the Relative Energy Gradient (REG) method, which calculates the correlation between these energies and the total energy of a system, thereby explaining the role atoms have in the energetic behaviour of the total system. The traditional "back-of-the-envelope" open and closed conformations of molecular balances do not correspond to those of the lowest energy. Hence, more care needs to be taken when considering which geometries to use for comparison with the experiment. The REG-IQA method shows that the 1-H and 1-OMe balances behave differently to the 1-Me and 1-NMe2 balances because the latter show more prominent electrostatics between carbonyl groups and undergoes a larger dihedral rotation due to the bulkiness of the functional groups. For the 2-Y balance, REG-IQA shows the same behaviour across the series as the 1-H and 1-OMe balances. From an atomistic point of view, the formation of the closed conformer is favoured by polarisation and charge-transfer effects on the amide bond across all balances and is counterbalanced by a de-pyramidalisation of the amide nitrogen. Moreover, focusing on the oxygen of the amide carbonyl and the α-carbon of the remaining carbonyl group, electrostatics have a major role in the formation of the closed conformer, which goes against the well-known n-π* interaction orbital overlap concept. However, REG-IQF shows that exchange-correlation energies overtake electrostatics for all the 2-Y balances when working with fragments around the carbonyl groups, while they act on par with electrostatics for the 1-OMe and 1-NMe2. REG-IQF also shows that exchange-correlation energies in the 2-Y balance are correlated to the inductive electron-donating and -withdrawing trends on aromatic groups. We demonstrate that methods such as REG-IQA/IQF can help with the fine-tuning of molecular balances prior to the experiment and that the energies that govern the probed interactions are highly dependent on the atoms and functional groups involved.

Contemporary Applications of Thioamides and Methods for Their Synthesis

Thioamides are naturally occurring isosteres of amide bonds in which the chalcogen atom of the carbonyl is changed from oxygen to sulfur. This substitution gives rise to altered nucleophilicity and hydrogen bonding properties with importance for both chemical reactivity and non-covalent interactions. As such, thioamides have been introduced into biologically active compounds to achieve improved target affinity and/or stability towards hydrolytic enzymes but have also been applied as probes of protein and peptide folding and dynamics. Recently, a series of new methods have been developed for the synthesis of thioamides as well as their utilization in peptide chemistry. Further, novel strategies for the incorporation of thioamides into proteins have been developed, enabling both structural and functional studies to be performed. In this Review, we highlight the recent developments in the preparation of thioamides and their applications for peptide modification and study of protein function.

Design, synthesis, and characterization of type I collagen mimetic peptides

The current wound-healing collagen mimetic peptides (CMPs) have limitations such as poor membrane permeability and protease susceptibility. Herein, the solid-phase peptide synthesis of CMPs containing the integrin binding motif GFOGER is reported. The peptide sequences also consist of lipophilic moieties (adamantane and palmitic acid) for improved membrane permeability and different collagen-inducing tripeptides, namely, Thr-Thr-Lys (TTK), Gly-His-Lys (GHK), Gln-Pro-Arg (QPR), and Glu-Glu-Met (EEM). The synthesized peptides were successfully characterized and purified using liquid chromatography-mass spectrometry and preparative high-performance liquid chromatography techniques, respectively. The palmitic acid moiety increased the hydrophobic nature of the peptides, and they were retained longer on the stationary material of the reverse phase C-18 column. The three-dimensional parallel-strand helical structure of peptide DGD-GG-GFOGER-GG-TTK-palmitate was obtained using nuclear magnetic resonance spectroscopy and circular dichroism. The synthesized peptides have the desired helical structure, which can promote integrin binding.

Molecular orbital and topological electron density study of n → π* interactions: amides and thioamides cases

The n → π* interactions were studied in amides and thioamides systems models, through the analysis of the electron density topology along with the Natural Bonding Orbital (NBO) approach. The effect of the dispersion terms was assessed using different DFT functionals. The NBO, independent gradient model (IGM), and the analysis of the reduced density gradient outcomes show that dispersion forces play a significant role in the strength of n → π* interactions. The IGM results indicate that δg height values for n → π* interactions do not extend beyond 0.025. All the methods used in this work predict that n → π* interaction between pairs of thioamides is stronger than those between amides. However, the electron density topology-based methods were not able to replicate the trends in the relative force of this interaction found in the experimental and NBO results.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

An Inherent Difference between Serine and Threonine Phosphorylation: Phosphothreonine Strongly Prefers a Highly Ordered, Compact, Cyclic Conformation

Phosphorylation and dephosphorylation of proteins by kinases and phosphatases are central to cellular responses and function. The structural effects of serine and threonine phosphorylation were examined in peptides and in proteins, by circular dichroism, NMR spectroscopy, bioinformatics analysis of the PDB, small-molecule X-ray crystallography, and computational investigations. Phosphorylation of both serine and threonine residues induces substantial conformational restriction in their physiologically more important dianionic forms. Threonine exhibits a particularly strong disorder-to-order transition upon phosphorylation, with dianionic phosphothreonine preferentially adopting a cyclic conformation with restricted ϕ (ϕ ∼ -60°) stabilized by three noncovalent interactions: a strong intraresidue phosphate-amide hydrogen bond, an n → π* interaction between consecutive carbonyls, and an n → σ* interaction between the phosphate Oγ lone pair and the antibonding orbital of C-Hβ that restricts the χ2 side-chain conformation. Proline is unique among the canonical amino acids for its covalent cyclization on the backbone. Phosphothreonine can mimic proline's backbone cyclization via noncovalent interactions. The preferred torsions of dianionic phosphothreonine are ϕ,ψ = polyproline II helix > α-helix (ϕ ∼ -60°); χ1 = g-; χ2 ∼ +115° (eclipsed C-H/O-P bonds). This structural signature is observed in diverse proteins, including in the activation loops of protein kinases and in protein-protein interactions. In total, these results suggest a structural basis for the differential use and evolution of threonine versus serine phosphorylation sites in proteins, with serine phosphorylation typically inducing smaller, rheostat-like changes, versus threonine phosphorylation promoting larger, step function-like switches, in proteins.

The Yin and Yang of How N-Terminal Acyl Caps Affect Collagen Triple Helices

N-terminal acylation is a common tool for the installation of functional moieties (e.g., sensors or bioactive molecules) on collagen model peptides (CMPs). The N-acyl group and its length are generally assumed to have little or no influence on the properties of the collagen triple helix formed by the CMP. Here, we show that the length of short (C1-C4) acyl capping groups has different effects on the thermal stability of collagen triple helices in POG, OGP, and GPO frames. While the effect of different capping groups on the stability of triple helices in the GPO frame is negligible, longer acyl chains stabilize OGP triple helices but destabilize POG analogues. The observed trends arise from a combination of steric repulsion, the hydrophobic effect, and n → π* interactions. Our study provides a basis for the design of N-terminally functionalized CMPs with predictable effects on triple helix stability.

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.

A fluoro-alkene mimic of Gly-trans-Pro produces a stable collagen triple helix

We report the first experimental evidence for a fluoro-alkene amide isostere participating in n→π* donation, which stabilizes the collagen triple helix. Of the three amide positions in canonical collagen-like peptides, Gly-Pro, Pro-Hyp, and Hyp-Gly, triple helix stability stands to benefit from substitution of only the isomerizable 3° Gly-Pro amide bond with a trans-locked fluoro-alkene. A (Z)-fluoro-alkene isostere of Gly-trans-Pro was synthesized, and its effect on the thermostability of a collagen-like peptide triple helix was measured. The mixture of enantiomers, Boc-Gly-Ψ[(Z)CFC]-L/D-Pro-OH, was synthesized in 8 steps with 27% overall yield, and the Fmoc-Gly-Ψ[(Z)CFC]-L/D-Pro-Hyp-OBn diastereomers were separated. The Gly-Ψ[(Z)CFC]-Pro isostere installed in a collagen-like peptide forms a stable triple helix. By CD, the thermal melting (T_m) value of the fluoro-alkene peptide was +42.2 ± 0.4 °C, and the T_m value of the control peptide was +48.4 ± 0.5 °C, a difference in stability of ΔT_m -6.2 °C. Deshielding of the fluorine nucleus in the ¹⁹F NMR spectra is evidence of a stabilizing n→π* electronic interaction.

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.

Inviting C5-Trifluoromethylated Pseudoprolines into Collagen Mimetic Peptides

Numerous collagen mimetic peptides (CMPs) have been engineered using proline derivatives substituted at their C(3) and/or C(4) position in order to stabilize or functionalize collagen triple-helix mimics. However, no example has been reported so far with C(5) substitutions. Here, we introduce a fluorinated CMP incorporating trifluoromethyl groups at the C(5) position of pseudoproline residues. In tripeptide models, our CD, NMR, and molecular dynamics (MD) studies have shown that, when properly arranged, these residues meet the structural requirements for a triple-helix assembly. Two host-guest CMPs were synthesized and analyzed by CD spectroscopy. The NMR analysis in solution of the most stable confirmed the presence of structured homotrimers that we interpret as triple helices. MD calculations showed that the triple-helix model remained stable throughout the simulation with all six trifluoromethyl groups pointing outward from the triple helix. Pseudoprolines substituted at the C(5) positions appeared as valuable tools for the design of new fluorinated collagen mimetic peptides.

A Unique and Stable Polyproline I Helix Sorted out from Conformational Equilibrium by Solvent Polarity

Polyproline I helical structures are often considered as the hidden face of their most famous geminal sibling, Polyproline II, as PPI is generally spotted only within a conformational equilibrium. We designed and synthesized a stable Polyproline I structure exploiting the striking tendency of (S)-indoline-2-carboxylic acid to drive the peptide bond conformation toward the cis amide isomer, when dissolved in polar solvents. The cooperative effect of only four amino acidic units is sufficient to form a preferential structure in solution. We shed light on this rare secondary structure with a thorough analysis of the spectroscopic and chiroptical properties of the tetramer, supported by X-ray crystallography and computational studies.

C. elegans as a model for inter-individual variation in metabolism

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors1-4. It is difficult to connect differences in metabolism to genomic variation and derive underlying molecular mechanisms in humans, owing to differences in diet and lifestyle, among others. Here we use the nematode Caenorhabditis elegans as a model to study inter-individual variation in metabolism. By comparing three wild strains and the commonly used N2 laboratory strain, we find differences in the abundances of both known metabolites and those that have not to our knowledge been previously described. The latter metabolites include conjugates between 3-hydroxypropionate (3HP) and several amino acids (3HP-AAs), which are much higher in abundance in one of the wild strains. 3HP is an intermediate in the propionate shunt pathway, which is activated when flux through the canonical, vitamin-B12-dependent propionate breakdown pathway is perturbed5. We show that increased accumulation of 3HP-AAs is caused by genetic variation in HPHD-1, for which 3HP is a substrate. Our results suggest that the production of 3HP-AAs represents a 'shunt-within-a-shunt' pathway to accommodate a reduction-of-function allele in hphd-1. This study provides a step towards the development of metabolic network models that capture individual-specific differences of metabolism and more closely represent the diversity that is found in entire species.

Characterizing the n→π* interaction of pyridine with small ketones: a rotational study of pyridine⋯acetone and pyridine⋯2-butanone

Complexes formed by pyridine and small ketones such as acetone and 2-butanone have been generated in a supersonic jet and characterized by broadband Fourier transform microwave spectroscopy combined with high-level theoretical computations. The spectra of the complexes show a quadrupole coupling hyperfine structure due to the presence of a nitrogen atom and the splittings owing to the low barriers of the internal rotation of the methyl groups bonded to the carbonyl group. The corresponding barriers have been determined from the analysis of the spectra. We show in both complexes that pyridine closes a cycle with a ketone carbonyl group through an N⋯CO n→π* tetrel interaction and a C-H⋯O contact. The n→π* tetrel bond involves the pyridine N atom lone pair and the ketone carbonyl group with a geometry approaching the Bürgi-Dunitz trajectory for the nucleophilic attack to a carbonyl group.

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.

Prolyl Isomerization-Mediated Conformational Changes Define ATR Subcellular Compartment-Specific Functions

ATR is a PI3K-like kinase protein, regulating checkpoint responses to DNA damage and replication stress. Apart from its checkpoint function in the nucleus, ATR actively engages in an antiapoptotic role at mitochondria following DNA damage. The different functions of ATR in the nucleus and cytoplasm are carried out by two prolyl isomeric forms of ATR: trans- and cis-ATR, respectively. The isomerization occurs at the Pin1 Ser428-Pro429 motif of ATR. Here, we investigated the structural basis of the subcellular location-specific functions of human ATR. Using a mass spectrometry-based footprinting approach, the surface accessibility of ATR lysine residues to sulfo-NHS-LC-biotin modification was monitored and compared between the cis- and the trans-isomers. We have identified two biotin-modified lysine residues, K459 and K469, within the BH3-like domain of cis-ATR that were not accessible in trans-ATR, indicating a conformational change around the BH3 domain between cis- and trans-ATR. The conformational alteration also involved the N-terminal domain and the middle HEAT domain. Moreover, experimental results from an array of complementary assays show that cis-ATR with the accessible BH3 domain was able to bind to tBid while trans-ATR could not. In addition, both cis- and trans-ATR can directly form homodimers via their C-terminal domains without ATRIP, while nuclear (trans-ATR) in the presence of ATRIP forms dimer-dimer complexes involving both N- and C-termini of ATR and ATRIP after UV. Structural characteristics around the Ser428-Pro429 motif and the BH3 domain region are also analyzed by molecular modeling and dynamics simulation. In support, cis conformation was found to be significantly more energetically favorable than trans at the Ser428-Pro429 bond in a 20-aa wild-type ATR peptide. Taken together, our results suggest that the isomerization-induced structural changes of ATR define both its subcellular location and compartment-specific functions and play an essential role in promoting cell survival and DNA damage responses.

Fluorine-Containing Prolines: Synthetic Strategies, Applications, and Opportunities

Fluorinated prolines play an important role in peptide studies, protein engineering, medicinal chemistry, drug discovery, and agrochemistry. Since the first synthesis of 4-fluoroprolines by Gottlieb and Witkop in 1965, their popularity started to grow exponentially. For example, during the past two decades, all isomeric trifluoromethyl-substituted prolines have been synthesized. In this Perspective, chemical properties and applications of fluorinated prolines are discussed. Synthetic approaches to all known fluorine-containing prolines are also discussed and analyzed. This analysis unexpectedly revealed an unsolved problem: in strict contrast to fluoro- and trifluoromethyl-substituted prolines, the corresponding analogues with fluoromethyl and difluoromethyl groups are mostly unknown. At the end of the paper, structures of several interesting, yet unknown, fluorinated prolines are disclosed─a good opportunity for chemists to make them.

Leveraging the n→π* Interaction in Alkene Isomerization by Selective Energy Transfer Catalysis

Examples of geometric alkene isomerization in nature are often limited to the net exergonic direction (ΔG°<0), with the antipodal net endergonic processes (ΔG°>0) comparatively under-represented. Inspired by the expansiveness of the maleate to fumarate (Z→E) isomerization in biochemistry, we investigated the inverse E→Z variant to validate nO →πC=O * interactions as a driving force for contra-thermodynamic isomerization. A general protocol involving selective energy transfer catalysis with inexpensive thioxanthone as a sensitizer (λmax =402 nm) is disclosed. Whilst in the enzymatic process nO →πC=O * interactions commonly manifest themselves in the substrate, these same interactions are shown to underpin directionality in the antipodal reaction by shortening the product alkene chromophore. The process was validated with diverse fumarate derivatives (>30 examples, up to Z:E>99:1), including the first examples of tetrasubstituted alkenes, and the involvement of nO →πC=O * interactions was confirmed by X-ray crystallography.

Facile synthesis of eight-membered cyclic(ester-amide)s and their organocatalytic ring-opening polymerizations

Facile modular synthesis of eight-membered cyclic(ester-amide)s based on phthalic anhydride and β-amino alcohols and organocatalitic ROP of the monomers to afford degradable semi-aromatic poly(ester-amides)s with tunable thermal properties.

Deciphering the Backbone Noncovalent Interactions that Stabilize Polyproline II Conformation and Reduce cis Proline Abundance in Polyproline Tracts

Proline (Pro) has a higher propensity to adopt cis amide geometry than the other natural amino acids, and a poly-Pro (poly-P) tract can adopt either a polyproline I (PPI, all cis amide) or a polyproline II (PPII, all trans amide) helical conformation. Recent studies have revealed a reduced abundance of cis amide geometry among the inner Pro residues of a poly-P tract. However, the forces that stabilize the polyproline helices and the reason for the higher trans amide propensity of the inner Pro residues of a poly-P tract are poorly understood. Herein, we have studied both Pro and non-Pro PPII helical sequences and identified the backbone noncovalent interactions that are crucial to the higher stability of the trans Pro-amide geometry and the preference for a PPII helical conformation. We show the presence of reciprocal CO···CO interactions that extend over the whole PPII helical region. Interestingly, the CO···CO interactions strengthen with the increase in the PPII helical chain length and the inner CO groups possess stronger CO···CO interactions, which could explain the reduced cis abundance of the inner Pro residues of a poly-P tract. We also identified a much stronger (∼0.9 kcal·mol^-1) n_O → σ*_Cα-Cβ interaction between the N-terminal CO oxygen lone pair and the antibonding orbital (σ*) of their C_α-C_β bonds. As the n_O → σ*_Cα-Cβ interaction is possible only in the trans isomers of Pro, this interaction should be crucial for the stabilization of a PPII helix. Finally, an unusual n_N(amide) → σ*_C-N interaction (∼0.3 kcal·mol^-1) was observed between the peptidic nitrogen lone pair (n_N) and the antibonding orbital (σ*_C-N) of the subsequent C-terminal peptide C-N bond. We propose a cumulative effect of these interactions in the stabilization of a PPII helix.

A Proline Mimetic for the Design of New Stable Secondary Structures: Solvent-Dependent Amide Bond Isomerization of (S)-Indoline-2-carboxylic Acid Derivatives

A thorough experimental and computational study on the conformational properties of (S)-indoline-2-carboxylic acid derivatives has been conducted. Methyl (S)-1-acetylindoline-2-carboxylate, both a mimetic of proline and phenylalanine, shows a remarkable tendency toward the cis amide isomer when dissolved in polar solvents. This behavior is opposite to the general preference of proline for the trans isomer, making indoline-2-carboxylic acid a good candidate for the design of different secondary structures and new materials.

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.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

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.

A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach

The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.

Intrinsic Conformational Dynamics of Alanine in Water/Ethanol Mixtures: An Experiment-Driven Molecular Dynamics Study

In vitro, cationic glycylalanylglycine (GAG) forms a hydrogel in binary mixtures of water and ethanol. In water, alanine residue is known for its high polyproline II (pPII) content. Spectroscopic data, including three J-coupling constants and amide I' profiles, indicate that addition of 42% ethanol to water significantly reduces the pPII content of alanine residue in GAG. Here, experiment-based Gaussian Ramachandran distributions of alanine in GAG at different ethanol fractions are examined and three MD force fields are evaluated with respect to their ability to capture these ethanol-induced conformational changes. MD simulations on monomeric GAG in eight different water/ethanol mixtures within Amber ff14SB, OPLS-AA/M, and CHARMM36m reveal that only Amber ff14SB partially captures the ethanol-induced conformational changes of alanine residue in monomeric GAG when 42% ethanol is added to water. MD simulations of 200 mM GAG ensembles in pure water and in the aqueous solution with 42% ethanol showcase the ability of CHARMM36m to capture the effect of ethanol on the average pPII content of alanine in GAG and provide a plausible explanation for this effect, which may stem from an increased propensity of GAG to form oligomers in the presence of ethanol.

Reconciling Electrostatic and n→π* Orbital Contributions in Carbonyl Interactions

Interactions between carbonyl groups are prevalent in protein structures. Earlier investigations identified dominant electrostatic dipolar interactions, while others implicated lone pair n→π* orbital delocalisation. Here these observations are reconciled. A combined experimental and computational approach confirmed the dominance of electrostatic interactions in a new series of synthetic molecular balances, while also highlighting the distance-dependent observation of inductive polarisation manifested by n→π* orbital delocalisation. Computational fiSAPT energy decomposition and natural bonding orbital analyses correlated with experimental data to reveal the contexts in which short-range inductive polarisation augment electrostatic dipolar interactions. Thus, we provide a framework for reconciling the context dependency of the dominance of electrostatic interactions and the occurrence of n→π* orbital delocalisation in C=O⋅⋅⋅C=O interactions.

Increasing the Functional Group Diversity in Helical β-Peptoids: Achievement of Solvent- and pH-Dependent Folding

We report the synthesis of a series of bis-functionalized β-peptoid oligomers of the hexamer length. This was achieved by synthesizing and incorporating protected amino- or azido-functionalized chiral building blocks into precursor oligomers by a trimer segment coupling strategy. The resulting hexamers were readily elaborated to provide target compounds displaying amino groups, carboxy groups, hydroxy groups, or triazolo-pyridines, which should enable metal ion binding. Analysis of the novel hexamers by circular dichroism (CD) spectroscopy and ¹H-¹³C heteronuclear single quantum coherence nuclear magnetic resonance (HSQC NMR) spectroscopy revealed robust helical folding propensity in acetonitrile. CD analysis showed a solvent-dependent degree of helical content in the structural ensembles when adding different ratios of protic solvents including an aqueous buffer. These studies were enabled by a substantial increase in solubility compared to previously analyzed β-peptoid oligomers. This also allowed for the investigation of the effect of pH on the folding propensity of the amino- and carboxy-functionalized oligomers, respectively. Interestingly, we could show a reversible effect of sequentially adding acid and base, resulting in a switching between compositions of folded ensembles with varying helical content. We envision that the present discoveries can form the basis for the development of functional peptidomimetic materials responsive to external stimuli.

Turn Conformation of β-Amino Acid-Based Short Peptides Promoted by an Amidothiourea Moiety at C-Terminus

A C-terminal amidothiourea motif is shown to promote a β-turn-like folded conformation in a series of β-amino acid-based short peptides in both the solid state and solution phase by an intramolecular 11-membered ring hydrogen bond.

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.

Phosphorylation-Dependent Pin1 Isomerization of ATR: Its Role in Regulating ATR's Anti-apoptotic Function at Mitochondria, and the Implications in Cancer

Peptidyl-prolyl isomerization is an important post-translational modification of protein because proline is the only amino acid that can stably exist as cis and trans, while other amino acids are in the trans conformation in protein backbones. This makes prolyl isomerization a unique mechanism for cells to control many cellular processes. Isomerization is a rate-limiting process that requires a peptidyl-prolyl cis/trans isomerase (PPIase) to overcome the energy barrier between cis and trans isomeric forms. Pin1, a key PPIase in the cell, recognizes a phosphorylated Ser/Thr-Pro motif to catalyze peptidyl-prolyl isomerization in proteins. The significance of the phosphorylation-dependent Pin1 activity was recently highlighted for isomerization of ATR (ataxia telangiectasia- and Rad3-related). ATR, a PIKK protein kinase, plays a crucial role in DNA damage responses (DDR) by phosphorylating hundreds of proteins. ATR can form cis or trans isomers in the cytoplasm depending on Pin1 which isomerizes cis-ATR to trans-ATR. Trans-ATR functions primarily in the nucleus. The cis-ATR, containing an exposed BH3 domain, is anti-apoptotic at mitochondria by binding to tBid, preventing activation of pro-apoptotic Bax. Given the roles of apoptosis in many human diseases, particularly cancer, we propose that cytoplasmic cis-ATR enables cells to evade apoptosis, thus addicting cancer cells to cis-ATR formation for survival. But in normal DDR, a predominance of trans-ATR in the nucleus coordinates with a minimal level of cytoplasmic cis-ATR to promote DNA repair while preventing cell death; however, cells can die when DNA repair fails. Therefore, a delicate balance/equilibrium of the levels of cis- and trans-ATR is required to ensure the cellular homeostasis. In this review, we make a case that this anti-apoptotic role of cis-ATR supports oncogenesis, while Pin1 that drives the formation of trans-ATR suppresses tumor growth. We offer a potential, novel target that can be specifically targeted in cancer cells, without killing normal cells, to significantly reduce the adverse effects usually seen in cancer treatment. We also raise important issues regarding the roles of phosphorylation-dependent Pin1 isomerization of ATR in diseases and propose areas of future studies that would shed more understanding on this important cellular mechanism.

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.

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.

Analysis of Density Functional Tight Binding with Natural Bonding Orbitals

The description of chemical bonding by the density functional tight binding (DFTB) model is analyzed using natural bonding orbitals (NBOs) and compared to results from density functional theory (B3LYP/aug-cc-pVTZ) calculations. Several molecular systems have been chosen to represent fairly diverse bonding scenarios that include standard covalent bonds, hypervalent interactions, multicenter bonds, metal-ligand interactions (with and without the pseudo-Jahn-Teller effect), and through-space donor-acceptor interactions. Overall, the results suggest that DFTB3/3OB provides physically sound descriptions for the different bonding scenarios analyzed here, as reflected by the general agreement between DFTB3 and B3LYP NBO properties, such as the nature of the NBOs, the magnitudes of natural charges and bond orders, and the dominant donor-acceptor interactions. The degree of ligand-to-metal charge transfer and the ionic nature of pentavalent phosphate are overestimated, likely reflecting the minimal-basis nature of DFTB3/3OB. Moreover, certain orbital interactions, such as geminal interactions, are observed to be grossly overestimated by DFTB3 for hypervalent phosphate and several transition metal compounds that involve copper and nickel. The study indicates that results from NBO analysis can be instructive for identifying electronic structure descriptions at the approximate quantum-mechanical level that require improvement and thus for guiding the systematic improvement of these methods.

Electronic and Steric Control of n→π* Interactions: Stabilization of the α-Helix Conformation without a Hydrogen Bond

The preferred conformations of peptides and proteins are dependent on local interactions that bias the conformational ensemble. The n→π* interaction between consecutive carbonyls promotes compact conformations, including the α-helix and polyproline II helix. In order to further understand the n→π* interaction and to develop methods to promote defined conformational preferences through acyl N-capping motifs, a series of peptides was synthesized in which the electronic and steric properties of the acyl group were modified. Using NMR spectroscopy, van't Hoff analysis of enthalpies, X-ray crystallography, and computational investigations, we observed that more electron-rich donor carbonyls (pivaloyl, iso-butyryl, propionyl) promote stronger n→π* interactions and more compact conformations than acetyl or less electron-rich donor carbonyls (methoxyacetyl, fluoroacetyl, formyl). X-ray crystallography indicates a strong, electronically tunable preference for the α-helix conformation, as observed directly on the φ and ψ torsion angles. Electron-donating acyl groups promote the α-helical conformation, even in the absence of the hydrogen bonding that stabilizes the α-helix. In contrast, electron-withdrawing acyl groups led to more extended conformations. More sterically demanding groups can promote trans amide bonds independent of the electronic effect on n→π* interactions. Chloroacetyl groups additionally promote n→π* interactions through the interaction of the chlorine lone pair with the proximal carbonyl π*. These data provide additional support for an important role of n→π* interactions in the conformational ensemble of disordered or unfolded proteins. Moreover, this work suggests that readily incorporated acyl N-capping motifs that modulate n→π* interactions may be employed rationally to promote conformational biases in peptides, with potential applications in molecular design and medicinal chemistry.

Strategy and Effects of Polyproline Peptide Stapling by Copper(I)-Catalyzed Alkyne-Azide Cycloaddition Reaction

Polyproline is a unique type of peptide that has a stable, robust, and well-defined helical structure in an aqueous environment. These features have allowed polyproline to be used as a nanosized scaffold for applications in chemical biology and related fields. To understand its structural properties and to expand the applications, this secondary structure was tested systematically by stapling the peptide at different locations with staples of various lengths. Using the efficient copper(I)-catalyzed alkyne-azide cycloaddition (CuAAC), we successfully prepared stapled polyproline and investigated the impact of this peptide macrocyclization through circular dichroism analysis. Whereas the stapling seems to have no significant effect on polyproline helix II (PPII) conformation in water, the location and the length of the staple affect the transformation of conformation in n-propanol. These results provide valuable information for future research using peptide stapling to manipulate polyproline conformation for various applications.

Conformational Analysis of n→π* Interactions in Collagen Triple Helix Models

Ab initio calculations of three models of collagen at positions Pro-Pro-Gly (1), Pro-Gly-Pro (2), and Gly-Pro-Pro (3) were performed to assess the conformational variation of n→π* contributions to the stability of the collagen triple helix. Full conformational analyses by relaxed potential-energy scans of the Ψ dihedral angle of the central residue in models 1, 2, and 3 revealed the presence of several n→π* interactions. In model 2, with Gly as the central residue, both the Φ and Ψ dihedral angles of Gly were scanned. Most minima of each model contained one or two n→π* interactions, with pyramidalization at the π* carbon. We also observed pyramidalization at the n→π* donor amide nitrogens. Minima with hydrogen-bond or non-native n→π* interactions compete with the collagen stabilizing n→π* interactions. The collagen-like n→ re-π* conformation was found as the global minimum only in model 3. The global minimum of 1 had a 5-membered ring hydrogen bond with an additional weak n→ si-π* interaction. The global minimum of 2 was in the extended conformation. We predict that the n→π* interactions found in native collagen, while individually small, cumulatively contribute to the stability of the triple helix conformation of collagen.