A prevalent intraresidue hydrogen bond stabilizes proteins

Helical twists and β-turns in structures at serine-proline sequences: Stabilization of cis-proline and type VI β-turns via C-H/O interactions

Structures at serine-proline sites in proteins were analyzed using a combination of peptide synthesis with structural methods and bioinformatics analysis of the PDB. Dipeptides were synthesized with the proline derivative (2S,4S)-(4-iodophenyl)hydroxyproline [hyp(4-I-Ph)]. The crystal structure of Boc-Ser-hyp(4-I-Ph)-OMe had two molecules in the unit cell. One molecule exhibited cis-proline and a type VIa2 β-turn (BcisD). The cis-proline conformation was stabilized by a C-H/O interaction between Pro C-Hα and the Ser side-chain oxygen. NMR data were consistent with stabilization of cis-proline by a C-H/O interaction in solution. The other crystallographically observed molecule had trans-Pro and both residues in the PPII conformation. Two conformations were observed in the crystal structure of Ac-Ser-hyp(4-I-Ph)-OMe, with Ser adopting PPII in one and the β conformation in the other, each with Pro in the δ conformation and trans-Pro. Structures at Ser-Pro sequences were further examined via bioinformatics analysis of the PDB and via DFT calculations. Ser-Pro versus Ala-Pro sequences were compared to identify bases for Ser stabilization of local structures. C-H/O interactions between the Ser side-chain Oγ and Pro C-Hα were observed in 45% of structures with Ser-cis-Pro in the PDB, with nearly all Ser-cis-Pro structures adopting a type VI β-turn. 53% of Ser-trans-Pro sequences exhibited main-chain COi•••HNi+3 or COi•••HNi+4 hydrogen bonds, with Ser as the i residue and Pro as the i + 1 residue. These structures were overwhelmingly either type I β-turns or N-terminal capping motifs on α-helices or 310-helices. These results indicate that Ser-Pro sequences are particularly potent in favoring these structures. In each, Ser is in either the PPII or β conformation, with the Ser Oγ capable of engaging in a hydrogen bond with the amide N-H of the i + 2 (type I β-turn or 310-helix; Ser χ1 t) or i + 3 (α-helix; Ser χ1 g+) residue. Non-proline cis amide bonds can also be stabilized by C-H/O interactions.

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.

Hydrogen Bonds Induced Ultralong Stability of Conductive π-d Conjugated FeCo3(DDA)2 with High OER Activity

Conductive π-d conjugated metal-organic frameworks (MOFs) have attracted wide concerns in electrocatalysis due to their intrinsic high conductivity. However, the poor electrocatalytic stability is still a major problem that hinders the practical application of MOFs. Herein, a novel approach to enhancing the stability of MOF-based electrocatalyst, namely, the introduction of hydrogen bonds (H-bonds), is reported. Impressively, the π-d conjugated MOF FeCo3(DDA)2 (DDA = 1,5-diamino-4,8-dihydroxy-9,10-anthraceneedione) exhibits ultrahigh oxygen evolution reaction (OER) stability (up to 2000 h). The experimental studies demonstrate that the presence of H-bonds in FeCo3(DDA)2 is responsible for its ultrahigh OER stability. Besides that, FeCo3(DDA)2 also displays a prominent OER activity (an overpotential of 260 mV vs reversible hydrogen electrode (RHE) at a current density of 10 mA cm-2 and a Tafel slope of 46.86 mV dec-1). Density functional theory (DFT) calculations further indicate that the synergistic effect of the Fe and Co sites in FeCo3(DDA)2 contributes to its prominent OER performance. This work provides a new avenue of boosting the electrocatalytic stability of conductive π-d conjugated MOFs.

Reconstructing dynamics correlation network to simultaneously improve activity and stability of 2,3-butanediol dehydrogenase by design of distal interchain disulfide bonds

The complete enzyme catalytic cycle includes substrate binding, chemical reaction and product release, in which different dynamic conformations are adopted. Due to the complex relationship among enzyme activity, stability and dynamics, the directed evolution of enzymes for improved activity or stability commonly leads to a trade-off in stability or activity. It hence remains a challenge to engineer an enzyme to have both enhanced activity and stability. Here, we have attempted to reconstruct the dynamics correlation network involved with active center to improve both activity and stability of a 2,3-butanediol dehydrogenase (2,3-BDH) by introducing inter-chain disulfide bonds. A computational strategy was first applied to evaluate the effect of introducing inter-chain disulfide bond on activity and stability of three 2,3-BDHs, and the N258C mutation of 2,3-BDH from Corynebacterium glutamicum (CgBDH) was proved to be effective in improving both activity and stability. In the results, CgBDH-N258C showed a different unfolding curve from the wild type, with two melting temperatures (Tm) of 68.3 °C and 50.8 °C, 19.7 °C and 2 °C higher than 48.6 °C of the wild type. Its half-life was also improved by 14.8-fold compared to the wild type. Catalytic efficiency (kcat/Km) of the mutant was increased by 7.9-fold toward native substrate diacetyl and 8.8-fold toward non-native substrate 2,5-hexanedione compared to the wild type. Molecular dynamics simulations revealed that an interaction network formed by Cys258, Arg162, Ala144 and the catalytic residues was reconstructed in the mutant and the dynamics change caused by the disulfide bond could be propagated through the interactions network. This improved the enzyme stability and activity by decreasing the flexibility and locking more "reactive" pose, respectively. Further construction of mutations including A144G showing a 44-fold improvement in catalytic efficiency toward meso-2,3-BD confirmed the role of modifying dynamics correlation network in tunning enzyme activity and selectivity. This study provided important insights into the relationship among dynamics, enzyme catalysis and stability, and will be useful in the designing new enzymes with co-evolution of stability, activity and selectivity.

Synthesis, derivatization, and conformational scanning of peptides containing N-Aminoglycine

N-alkylated glycine residues are the main constituent of peptoids and peptoid-peptide hybrids that are employed across the biomedical and materials sciences. While the impact of backbone N-alkylation on peptide conformation has been extensively studied, less is known about the effect of N-amination on the secondary structure propensity of glycine. Here, we describe a convenient protocol for the incorporation of N-aminoglycine into host peptides on solid support. Amide-to-hydrazide substitution also affords a nucleophilic handle for further derivatization of the backbone. To demonstrate the utility of late-stage hydrazide modification, we synthesized and evaluated the stability of polyproline II helix and β-hairpin model systems harboring N-aminoglycine derivatives. The described procedures provide facile entry into peptidomimetic libraries for conformational scanning.

Nucleic-acid-base photofunctional cocrystal for information security and antimicrobial applications

Cocrystal engineering is an efficient and simple strategy to construct functional materials, especially for the exploitation of novel and multifunctional materials. Herein, we report two kinds of nucleic-acid-base cocrystal systems that imitate the strong hydrogen bond interactions constructed in the form of complementary base pairing. The two cocrystals studied exhibit different colors of phosphorescence from their monomeric counterparts and show the feature of rare high-temperature phosphorescence. Mechanistic studies reveal that the strong hydrogen bond network stabilizes the triplet state and suppresses non-radiative transitions, resulting in phosphorescence even at 425 K. Moreover, the isolation effects of the hydrogen bond network regulate the interactions between the phosphor groups, realizing the manipulation from aggregation to single-molecule phosphorescence. Benefiting from the long-lived triplet state with a high quantum yield, the generation of reactive oxygen species by energy transfer is also available to utilize for some applications such as in photodynamic therapy and broad-spectrum microbicidal effects. In vitro experiments show that the cocrystals efficiently kill bacteria on a tooth surface and significantly help prevent dental caries. This work not only provides deep insight into the relationship of the structure-properties of cocrystal systems, but also facilitates the design of multifunctional cocrystal materials and enriches their potential applications.

Artificial Intelligence-based Amide-II Infrared Spectroscopy Simulation for Monitoring Protein Hydrogen Bonding Dynamics

The structurally sensitive amide II infrared (IR) bands of proteins provide valuable information about the hydrogen bonding of protein secondary structures, which is crucial for understanding protein dynamics and associated functions. However, deciphering protein structures from experimental amide II spectra relies on time-consuming quantum chemical calculations on tens of thousands of representative configurations in solvent water. Currently, the accurate simulation of amide II spectra for whole proteins remains a challenge. Here, we present a machine learning (ML)-based protocol designed to efficiently simulate the amide II IR spectra of various proteins with an accuracy comparable to experimental results. This protocol stands out as a cost-effective and efficient alternative for studying protein dynamics, including the identification of secondary structures and monitoring the dynamics of protein hydrogen bonding under different pH conditions and during protein folding process. Our method provides a valuable tool in the field of protein research, focusing on the study of dynamic properties of proteins, especially those related to hydrogen bonding, using amide II IR spectroscopy.

Fibril structures of TFG protein mutants validate the identification of TFG as a disease-related amyloid protein by the IMPAcT method

We previously presented a bioinformatic method for identifying diseases that arise from a mutation in a protein's low-complexity domain that drives the protein into pathogenic amyloid fibrils. One protein so identified was the tropomyosin-receptor kinase-fused gene protein (TRK-fused gene protein or TFG). Mutations in TFG are associated with degenerative neurological conditions. Here, we present experimental evidence that confirms our prediction that these conditions are amyloid-related. We find that the low-complexity domain of TFG containing the disease-related mutations G269V or P285L forms amyloid fibrils, and we determine their structures using cryo-electron microscopy (cryo-EM). These structures are unmistakably amyloid in nature and confirm the propensity of the mutant TFG low-complexity domain to form amyloid fibrils. Also, despite resulting from a pathogenic mutation, the fibril structures bear some similarities to other amyloid structures that are thought to be nonpathogenic and even functional, but there are other factors that support these structures' relevance to disease, including an increased propensity to form amyloid compared with the wild-type sequence, structure-stabilizing influence from the mutant residues themselves, and double-protofilament amyloid cores. Our findings elucidate two potentially disease-relevant structures of a previously unknown amyloid and also show how the structural features of pathogenic amyloid fibrils may not conform to the features commonly associated with pathogenicity.

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

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.

Effects of chirality and side chain length in Cα,α-dialkylated residues on β-hairpin peptide folded structure and stability

Strategic incorporation of achiral Cα,α-dialkylated amino acids with bulky substituents into peptides can be used to promote extended strand conformations and inhibit protein-protein interactions associated with amyloid formation. In this work, we evaluate the thermodynamic impact of chiral Cα,α monomers on folding preferences in such systems through introduction of a series of Cα-methylated and Cα-ethylated residues into a β-hairpin host sequence. Depending on stereochemical configuration of the artificial monomer and potential for additional hydrophobic packing, a Cα-ethyl-Cα-propyl glycine residue can provide similar or enhanced folded stability relative to an achiral Cα,α-diethyl analogue.

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

The "Beacon" Structural Model of Protein Folding: Application for Trp-Cage in Water

Protein folding is a process in which a polypeptide must undergo folding process to obtain its three-dimensional structure. Thermodynamically, it is a process of enthalpy to overcome the loss of conformational entropy in folding. Folding is primarily related to hydrophobic interactions and intramolecular hydrogen bondings. During folding, hydrophobic interactions are regarded to be the driving forces, especially in the initial structural collapse of a protein. Additionally, folding is guided by the strong interactions within proteins, such as intramolecular hydrogen bondings related to the α-helices and β-sheets of proteins. Therefore, a protein is divided into the folding key (FK) regions related to intramolecular hydrogen bondings and the non-folding key (non-FK) regions. Various conformations are expected for FK and non-FK regions. Different from non-FK regions, it is necessary for FK regions to form the specific conformations in folding, which are regarded as the necessary folding pathways (or "beacons"). Additionally, sequential folding is expected for the FK regions, and the intermediate state is found during folding. They are reflected on the local basins in the free energy landscape (FEL) of folding. To demonstrate the structural model, molecular dynamics (MD) simulations are conducted on the folding pathway of the TRP-cage in water.

Length-Dependent Transition from Extended to Folded Shapes in Short Oligomers of an Azetidine-Based α-Amino Acid: The Critical Role of NH···N H-Bonds

Hydrogen bonds (H-bonds) are ubiquitous in peptides and proteins and are central to the stabilization of their structures. Inter-residue H-bonds between non-adjacent backbone amide NH and C=O motifs lead to the well-known secondary structures of helices, turns and sheets, but it is recognized that other H-bonding modes may be significant, including the weak intra-residue H-bond (called a C5 H-bond) that implicates the NH and C=O motifs of the same amino acid residue. Peptide model compounds that adopt stable C5 H-bonds are not readily available and the so-called 2.0₅-helix, formed by successive C5 H-bonds, is an elusive secondary structure. Using a combination of theoretical chemistry and spectroscopic studies in both the gas phase and solution phase, we have demonstrated that derivatives of 3-amino-1-methylazetidine-3-carboxylic acid, Aatc(Me) can form sidechain-backbone N-H···N C6γ H-bonds that accompany-and thereby stabilize-C5 H-bonds. In the capped trimer of Aatc(Me), extended C5/C6γ motifs are sufficiently robust to challenge classical 3₁₀-helix formation in solution and the fully-extended 2.0₅-helix conformer has been characterized in the gas phase. Concurrent H-bonding support for successive C5 motifs is a new axiom for stabilizing the extended backbone secondary structure in short peptides.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Disassembly of Amyloid Fibril with Infrared Free Electron Laser

Amyloid fibril causes serious amyloidosis such as neurodegenerative diseases. The structure is composed of rigid β-sheet stacking conformation which makes it hard to disassemble the fibril state without denaturants. Infrared free electron laser (IR-FEL) is an intense picosecond pulsed laser that is oscillated through a linear accelerator, and the oscillation wavelengths are tunable from 3 μm to 100 μm. Many biological and organic compounds can be structurally altered by the mode-selective vibrational excitations due to the wavelength variability and the high-power oscillation energy (10-50 mJ/cm²). We have found that several different kinds of amyloid fibrils in amino acid sequences were commonly disassembled by the irradiation tuned to amide I (6.1-6.2 μm) where the abundance of β-sheet decreased while that of α-helix increased by the vibrational excitation of amide bonds. In this review, we would like to introduce the IR-FEL oscillation system briefly and describe combination studies of experiments and molecular dynamics simulations on disassembling amyloid fibrils of a short peptide (GNNQQNY) from yeast prion and 11-residue peptide (NFLNCYVSGFH) from β2-microglobulin as representative models. Finally, possible applications of IR-FEL for amyloid research can be proposed as a future outlook.

Non-covalent interactions reveal the protein chain δ conformation in a flexible single-residue model

The δ conformation is a local secondary structure in proteins that implicates a π_amide N-H⋯N interaction between a backbone N atom and the NH of the following residue. Small-molecule models thereof have been limited so far to rigid proline-type compounds. We show here that in derivatives of a cyclic amino acid with a sulphur atom in the γ-position, specific side-chain/backbone N-H⋯S interactions stabilize the δ conformation sufficiently to allow it to compete with classical C5 and C7 H-bonded conformers.

Influence of Proline Chirality on Neighbouring Azaproline Residue Stereodynamic Nitrogen Preorganization

We report herein the first systematic crystal structural investigation of azaproline incorporated in homo- and heterochiral diprolyl peptides. The X-ray crystallography data of peptides 1-5 illustrates that stereodynamic nitrogen in azaproline adopted the stereochemistry of neighbouring proline residue without depending on its position in the peptide sequence. Natural bond orbital analysis of crystal structures indicates O_azPro -C'_Pro of peptides 4 and 5 participating in n→π* interaction with stabilization energy about 1.21-1.33 kcal/mol. Density functional theory calculations suggested that the endo-proline ring puckering favoured over exo-conformation by 6.72-7.64 kcal/mol. NBO and DFT data reveals that the n→π* interactions and proline ring puckering stabilize azaproline chirality with the neighbouring proline stereochemistry. The CD, solvent titration, variable-temperature and 2D NMR experimental results further supported the crystal structures conformation.

Detecting the Hydrogen Bond Cooperativity in a Protein β-Sheet by H/D Exchange

The hydrogen bond (H-bond) cooperativity in the β-sheet of GB3 is investigated by a NMR hydrogen/deuterium (H/D) exchange method. It is shown that the weakening of one backbone N-H…O=C H-bond between two β-strands, β1 and β2, due to the exchange of NH to ND of the H-bond donor in β1, perturbs the chemical shift of ¹³C_α, ¹³C_β, ¹H_α, ¹H_N, and ¹⁵N of the H-bond acceptor and its following residue in β2. Quantum mechanical calculations suggest that the -H-bond chemical shift isotope effect is caused by the structural reorganization in response to the H-bond weakening. This structural reorganization perturbs four neighboring H-bonds, with three being weaker and one being stronger, indicating that three H-bonds are cooperative and one is anticooperative with the perturbed H-bond. The sign of the cooperativity depends on the relative position of the H-bonds. This H-bond cooperativity, which contributes to β-sheet stability overall, can be important for conformational coupling across the β-sheet.

Aza-PNA: Engineering E-Rotamer Selectivity Directed by Intramolecular H-bonding

The replacement of α(CH₂) by NH in monomers of standard aeg PNA and its homologue β-ala PNA leads to respective aza-PNA monomers (1 and 2) in which the N_αH can form either an 8-membered H-bonded ring with folding of the backbone (DMSO and water) or a 5-membered N_αH─αCO (water) to stabilize E-type rotamers. Such aza-PNA oligomers with exclusive E rotamers and intraresidue backbone H-bonding can modulate its DNA/RNA binding and assembling properties.

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.

Right- and left-handed PPI helices in cyclic dodecapeptoids

Enantiomorphic right- and left-handed polyproline type I helices in four cyclic dodecapeptoids with methoxyethyl and propargyl side chains are observed for the first time by single crystal X-ray diffraction. The peculiar absence of NH⋯OC hydrogen bonds in peptoids unveils the role of intramolecular backbone-to-backbone CO⋯CO interactions and CH⋯OC hydrogen bonds in the stabilization of the macrocycle conformation. Moreover, intramolecular backbone-side chain C5 CH⋯OC hydrogen bonds emerge as a stabilizing factor.

Extended β-Strands Contribute to Reversible Amyloid Formation

The assembly of proteins into fibrillar amyloid structures was once considered to be pathologic and essentially irreversible. Recent studies reveal amyloid-like structures that form reversibly, derived from protein low-complexity domains which function in cellular metabolism. Here, by comparing atomic-level structures of reversible and irreversible amyloid fibrils, we find that the β-sheets of reversible fibrils are enriched in flattened (as opposed to pleated) β-sheets formed by stacking of extended β-strands. Quantum mechanical calculations show that glycine residues favor extended β-strands which may be stabilized by intraresidue interactions between the amide proton and the carbonyl oxygen, known as C5 hydrogen-bonds. Larger residue side chains favor shorter strands and pleated sheets. These findings highlight a structural element that may regulate reversible amyloid assembly.

Newly identified C–H⋯O hydrogen bond in histidine

Histidine C–H bonds observed in protein structures include (clockwise from top left): myoglobin, β-lactamase, and photoactive yellow protein; calculations indicate that tautomeric/protonation state influences H-bonding ability (bottom left).

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.

Bispidine as a β-strand nucleator: from a β-arch to self-assembled cages and vesicles

The development of synthetic scaffolds that nucleate well-folded secondary structures is highly challenging. Herein, we designed and synthesized a series of core-modified peptides (F1, F2, F3, and F4) that fold into β-strand structures. These bispidine-scaffolded peptides were studied by CD, IR, NMR, single crystal XRD, and Molecular Dynamics (MD) simulations to investigate their conformational preferences. Solid-state and solution studies revealed that bispidine is a versatile scaffold that could be placed either at the terminal or at the middle of the peptide strand for nucleating the β-strand structure. Scaffolds that nucleate an isolated β-strand conformation are rare. Bispidine placed at the C-terminus of the peptide chain could nucleate a β-strand conformation, while bispidine placed at the middle resulted in a β-arch conformation. This nucleation activity stems from the ability to restrict the psi torsion angle (ψ) through intramolecular C5 hydrogen bonding between the equatorial hydrogen(s) of bispidine and the carbonyl oxygen(s) of the amino acid close to the scaffold. Furthermore, the bispidine peptidomimetic with a super secondary structure, namely β-arch, assembled into single-hole submicron cages and spherical vesicles as evident from microscopic studies. The design logic defined here will be a significant strategy for the development of β-strand mimetics and super secondary structures.

N-H⋯X interactions stabilize intra-residue C5 hydrogen bonded conformations in heterocyclic α-amino acid derivatives

Nature makes extensive and elaborate use of hydrogen bonding to assemble and stabilize biomolecular structures. The shapes of peptides and proteins rely significantly on N–H⋯O Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> C interactions, which are the linchpins of turns, sheets and helices. The C5 H-bond, in which a single residue provides both donor and acceptor, is generally considered too weak to force the backbone to adopt extended structures. Exploiting the synergy between gas phase (experimental and quantum chemistry) and solution spectroscopies to decipher IR spectroscopic data, this work demonstrates that the extended C5-based conformation in 4-membered ring heterocyclic α-amino acid derivatives is significantly stabilized by the formation of an N–H⋯X H-bond. In this synergic system the strength of the C5 interaction remains constant while the N–H⋯X H-bond strength, and thereby the support provided by it, varies with the heteroatom.

In 4-membered ring heterocyclic α-amino acid derivatives, extended conformations based on intraresidue C5 H-bonds can be stabilized by N–H⋯X H-bonds, making the combined C5–C6γ structures prominent in both gas phase and in weakly polar solutions.

SARS-CoV-2 Variants, RBD Mutations, Binding Affinity, and Antibody Escape

Since 2020, the receptor-binding domain (RBD) of the spike protein of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been constantly mutating, producing most of the notable missense mutations in the context of "variants of concern", probably in response to the vaccine-driven alteration of immune profiles of the human population. The Delta variant, in particular, has become the most prevalent variant of the epidemic, and it is spreading in countries with the highest vaccination rates, causing the world to face the risk of a new wave of the contagion. Understanding the physical mechanism responsible for the mutation-induced changes in the RBD's binding affinity, its transmissibility, and its capacity to escape vaccine-induced immunity is the "urgent challenge" in the development of preventive measures, vaccines, and therapeutic antibodies against the coronavirus disease 2019 (COVID-19) pandemic. In this study, entropy-enthalpy compensation and the Gibbs free energy change were used to analyze the impact of the RBD mutations on the binding affinity of SARS-CoV-2 variants with the receptor angiotensin converting enzyme 2 (ACE2) and existing antibodies. Through the analysis, we found that the existing mutations have already covered almost all possible detrimental mutations that could result in an increase of transmissibility, and that a possible mutation in amino-acid position 498 of the RBD can potentially enhance its binding affinity. A new calculation method for the binding energies of protein-protein complexes is proposed based on the entropy-enthalpy compensation rule. All known structures of RBD-antibody complexes and the RBD-ACE2 complex comply with the entropy-enthalpy compensation rule in providing the driving force behind the spontaneous protein-protein docking. The variant-induced risk of breakthrough infections in vaccinated people is attributed to the L452R mutation's reduction of the binding affinity of many antibodies. Mutations reversing the hydrophobic or hydrophilic performance of residues in the spike RBD potentially cause breakthrough infections of coronaviruses due to the changes in geometric complementarity in the entropy-enthalpy compensations between antibodies and the virus at the binding sites.

Entropy-Enthalpy Compensations Fold Proteins in Precise Ways

Exploring the protein-folding problem has been a longstanding challenge in molecular biology and biophysics. Intramolecular hydrogen (H)-bonds play an extremely important role in stabilizing protein structures. To form these intramolecular H-bonds, nascent unfolded polypeptide chains need to escape from hydrogen bonding with surrounding polar water molecules under the solution conditions that require entropy-enthalpy compensations, according to the Gibbs free energy equation and the change in enthalpy. Here, by analyzing the spatial layout of the side-chains of amino acid residues in experimentally determined protein structures, we reveal a protein-folding mechanism based on the entropy-enthalpy compensations that initially driven by laterally hydrophobic collapse among the side-chains of adjacent residues in the sequences of unfolded protein chains. This hydrophobic collapse promotes the formation of the H-bonds within the polypeptide backbone structures through the entropy-enthalpy compensation mechanism, enabling secondary structures and tertiary structures to fold reproducibly following explicit physical folding codes and forces. The temperature dependence of protein folding is thus attributed to the environment dependence of the conformational Gibbs free energy equation. The folding codes and forces in the amino acid sequence that dictate the formation of β-strands and α-helices can be deciphered with great accuracy through evaluation of the hydrophobic interactions among neighboring side-chains of an unfolded polypeptide from a β-strand-like thermodynamic metastable state. The folding of protein quaternary structures is found to be guided by the entropy-enthalpy compensations in between the docking sites of protein subunits according to the Gibbs free energy equation that is verified by bioinformatics analyses of a dozen structures of dimers. Protein folding is therefore guided by multistage entropy-enthalpy compensations of the system of polypeptide chains and water molecules under the solution conditions.

Why Should Metformin Not Be Given in Advanced Kidney Disease? Potential Leads from Computer Simulations

Metformin is considered as the go-to drug in the treatment of diabetes. However, it is either prescribed in lower doses or not prescribed at all to patients with kidney problems. To find a potential explanation for this practice, we employed atomistic-level computer simulations to simulate the transport of metformin through multidrug and toxin extrusion 1 (MATE1), a protein known to play a key role in the expulsion of metformin into urine. Herein, we examine the hydrogen bonding between MATE1 and one or more metformin molecules. The simulation results indicate that metformin continuously forms and breaks off hydrogen bonds with MATE1 residues. However, the mean hydrogen bond lifetimes increase for an order of magnitude when three metformin molecules are inserted instead of one. This new insight into the metformin transport process may provide the molecular foundation behind the clinical practice of not prescribing metformin to kidney disease patients.

Stabilization of Azapeptides by Namide···H-Namide Hydrogen Bonds

An unusual N_amide···H-N_amide hydrogen bond (HB) was previously proposed to stabilize the azapeptide β-turns. Herein we provide experimental evidence for the N_amide···H-N_amide HB and show that this HB endows a stabilization of 1-3 kcal·mol^-1 and enforces the trans-cis-trans (t-c-t) and cis-cis-trans (c-c-t) amide bond conformations in azapeptides and N-methyl-azapeptides, respectively. Our results indicate that these N_amide···H-N_amide HBs can have stabilizing contributions even in short azapeptides that cannot fold to form β-turns.

Synthesis and conformation of backbone N-aminated peptides

The chemical modification of peptides is a promising approach for the design of protein-protein interaction inhibitors and peptide-based drug candidates. Among several peptidomimetic strategies, substitution of the amide backbone maintains side-chain functionality that may be important for engagement of biological targets. Backbone amide substitution has been largely limited to N-alkylation, which can promote cis amide geometry and disrupt important H-bonding interactions. In contrast, N-amination of peptides induces distinct backbone geometries and maintains H-bond donor capacity. In this chapter we discuss the conformational characteristics of designed N-amino peptides and present a detailed protocol for their synthesis on solid support. The described methods allow for backbone N-amino scanning of biologically active parent sequences.

Studies of the highly potent lantibiotic peptide nisin Z in aqueous solutions of salts and biological buffer components

The lantibiotic nisin, usually used as a 2.5%w/w in NaCl and milk solids, has activity against a wide range of Gram-positive bacteria, especially food-borne pathogens, and has been used as a food preservative for decades without the development of significant resistance. It has been reported that the high purity (>95%) nisin Z form has activity against the Gram-negative speciesE. coli, which is significantly reduced in the presence of NaCl. This current study examined, by¹H NMR spectroscopy, the effects of NaCl, and a range of other salts, on the observed aqueous solution¹H NMR spectra of nisin Z in the pH 3-4 range, where nisin Z has its maximum stability. Nisin's mechanism of action involves binding to the polyoxygenated pyrophosphate moiety of lipid II, and in acidic solution the positively charged C-terminus region is reported to interact with the negative sulfate groups of SDS micelles, so the study was extended to include a number of polyoxygenated anions commonly used as buffers in many biological assays. In general, the biggest changes found were in the chemical shifts of protons in the hydrophobic N-terminus region, rather than the more polar C-terminus region. The effects seen on the addition of the salts (cations and anions) were not just an overall non-specific ionic strength effect, as different salts caused different effects, in an unpredictive manner. Similarly, the polyoxygenated anions behaved differently and not predictably, and neither the cations/anions, or polyoxygenated anions, constitute a Hofmeister or inverse Hofmeister series.

Flat, Cα,β -Didehydroalanine Foldamers with Ferrocene Pendants: Assessing the Role of α-Peptide Dipolar Moments

The foldamer field is continuously expanding as it allows to produce molecules endowed with 3D-structures and functions never observed in nature. We synthesized flat foldamers based on the natural, but non-coded, Cα,β -didehydroalanine α-amino acid, and covalently linked to them two ferrocene (Fc) moieties, as redox probes. These conjugates retain the flat and extended conformation of the 2.05 -helix, both in solution and in the crystal state (X-ray diffraction). Cyclic voltammetry measurements agree with the adoption of the 2.05 -helix, characterized by a negligible dipole moment. Thus, elongated α-peptide stretches of this type are insulators rather than charge conductors, the latter being constituted by peptide α-helices. Also, our homo-tetrapeptide has a N-to-C length of about 18.2 Å, almost double than that (9.7 Å) of an α-helical α-tetrapeptide.

Experimental and computational diagnosis of the fluxional nature of a benzene-1,3,5-tricarboxamide-based hydrogen-bonded dimer

Precise characterization of the hydrogen bond network present in discrete self-assemblies of benzene-1,3,5-tricarboxamide monomers derived from amino-esters (ester BTAs) is crucial for the construction of elaborated functional co-assemblies. For all ester BTA dimeric structures previously reported, ester carbonyls in the side chain acted as hydrogen bond acceptors, yielding well-defined dimers stabilized by six hydrogen bonds. The ester BTA monomer derived from glycine (BTA Gly) shows a markedly different self-assembly behaviour. We report herein a combined experimental and computational investigation aimed at determining the nature of the dimeric species formed by BTA Gly. Two distinct dimeric structures were characterized by single-crystal X-ray diffraction measurements. Likewise, a range of spectroscopic and scattering techniques as well as molecular modelling were employed to diagnose the nature of dynamic dimeric structures in toluene. Our results unambiguously establish that both ester and amide carbonyls are involved in the hydrogen bond network of the discrete dimeric species formed by BTA Gly. The participation of roughly 4.5 ester carbonyls and 1.5 amide carbonyls per dimer as determined by FT-IR spectroscopy implies that several conformations coexist in solution. Moreover, NMR analysis and modelling data reveal rapid interconversion between these different conformers leading to a symmetric structure on the NMR timescale. Rapid hydrogen bond shuffling between conformers having three (three), two (four), one (five) and zero (six) amide carbonyl groups (ester carbonyl groups, respectively) as hydrogen bond acceptors is proposed to explain the magnetic equivalence of the amide N-H on the NMR timescale. When compared to other ester BTA derivatives in which only ester carbonyls act as hydrogen bond acceptors, the fluxional behaviour of the hydrogen-bonded dimers of BTA Gly likely originates from a larger range of energetically favorable conformations accessible through rotation of the BTA side chains.

The first non-helical Aib-containing hexapeptide: The crystal structure of Z-Gly-Aib-Gly-Aib-Gly-Aib-OtBu

The synthetic peptide Z-Gly-Aib-Gly-Aib-Gly-Aib-OtBu was crystallized from a mixture of ethyl acetate and n-hexane. The crystals belong to the centrosymmetric space group Pbca. There are three molecules in the asymmetric unit. The three molecules differ mainly in the Z-group conformation. The first Gly residue adopts a fully extended conformation, residues 2 and 3 lie in the left-handed helical region, residues 4 and 5 in the right-handed helical region, and residue 6 again in the left-handed helical region of the Ramachandran plot. There are only two of four possible intramolecular hydrogen bonds formed, namely, between Aib4 and Gly1 forming a β-turn of type III' and between Aib6 and Gly3 forming a β-turn of type I. The inverted molecules (by space group symmetry) lie in the regions with opposite handedness and form β-turns of type III and I'. In contrast to all known long synthetic and naturally occurring Aib-containing peptides that fold as 3₁₀ - or α-helix, Z-(Gly-Aib)₃ -OtBu folds in a quite flat structure from which only the protecting groups bulge out.

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.

Effects of Excipients on the Structure and Dynamics of Filgrastim Monitored by Thermal Unfolding Studies by CD and NMR Spectroscopy

Product excipients are used to confer a number of desirable properties on the drug substance to maintain or improve stability and facilitate drug delivery. This is especially important for products where the active pharmaceutical ingredient (API) is a recombinant protein. In this study, we aimed to determine if excipients and formulation conditions affect the structure and/or modulate the dynamics of the protein API of filgrastim products. Samples of uniformly labeled ¹⁵N-Met-granulocyte-colony stimulating factor (GCSF) were prepared at 100 μM (near formulation concentration) with various concentrations of individual components (polysorbate-20 and -80, sorbitol) and three pH values. Nuclear magnetic resonance (NMR) spectroscopy techniques were applied to measure chemical shift perturbation (CSP) to detect structural changes, and relaxation parameters (T ₁, T ₂, and heteronuclear Overhauser effect) were measured to probe the effects on protein backbone motions. In parallel, the same solution conditions were subjected to protein thermal unfolding studies monitored by circular dichroism spectropolarimetry (CD). Detergents (polysorbate-20 and 80) do not induce any observable changes on the protein structure and do not modify its dynamics at formulation concentration. Lowering pH to 4.0, a condition known to stabilize the conformation of filgrastim, as well as the addition of sorbitol produced changes of the fast motion dynamics in the nanosecond and picosecond timescale. NMR-derived order parameters, which measure the local conformational entropy of the protein backbone, show that lowering pH leads to a compaction of the four-helix bundle while the addition of sorbitol relaxes helices B and C, thereby reducing the mobility of loop CD. CSPs and measurements of protein dynamics via NMR-derived order parameters provide a description in structural and motional terms at an atomic resolution on how formulation components contribute to the stabilization of filgrastim products.

A Geometric Definition of Short to Medium Range Hydrogen-Mediated Interactions in Proteins

We present a method to rapidly identify hydrogen-mediated interactions in proteins (e.g., hydrogen bonds, hydrogen bonds, water-mediated hydrogen bonds, salt bridges, and aromatic π-hydrogen interactions) through heavy atom geometry alone, that is, without needing to explicitly determine hydrogen atom positions using either experimental or theoretical methods. By including specific real (or virtual) partner atoms as defined by the atom type of both the donor and acceptor heavy atoms, a set of unique angles can be rapidly calculated. By comparing the distance between the donor and the acceptor and these unique angles to the statistical preferences observed in the Protein Data Bank (PDB), we were able to identify a set of conserved geometries (15 for donor atoms and 7 for acceptor atoms) for hydrogen-mediated interactions in proteins. This set of identified interactions includes every polar atom type present in the Protein Data Bank except OE1 (glutamate/glutamine sidechain) and a clear geometric preference for the methionine sulfur atom (SD) to act as a hydrogen bond acceptor. This method could be readily applied to protein design efforts.

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.

Application of a Distance-Dependent Sigmoidal Dielectric Constant to the REMC/SAAP3D Simulations of Chignolin, Trp-Cage, and the G10q Mutant

The replica-exchange Monte Carlo method based on the single amino acid potential (SAAP) force field, i.e., REMC/SAAP3D, was recently developed by our group for the molecular simulation of short peptides. In this study, the method has been improved by applying a distance-dependent dielectric (DDD) constant and extended to the peptides containing D-amino acid (AA) residues. For chignolin (10 AAs), a sigmoidal DDD model reasonably allocated the native-like β-hairpin structure with all-atom root mean square deviation (RMSD) = 2.0 Å as a global energy minimum. The optimal DDD condition was subsequently applied for Trp-cage (20 AAs) and its G10q mutant. The native-like α-rich folded structures with main-chain RMSD = 3.7 and 3.8 Å were obtained as global energy minima for Trp-cage and G10q, respectively. The results suggested that the REMC/SAAP3D method with the sigmoidal DDD model is useful for structural prediction for the short peptides comprised of up to 20 AAs. In addition, the relative contributions of SAAP to the total energy (%SAAP) were evaluated by energetic component analysis. The ratios of %SAAP were about 40 and 20% for chignolin and Trp-cage (or G10q), respectively. It was proposed that SAAP is more important for the secondary structure formation than for the assembly to a higher-order folded structure, in which the attractive van der Waals interaction may play a more important role.

Polyol and sugar osmolytes can shorten protein hydrogen bonds to modulate function

Polyol and sugar osmolytes are commonly used in therapeutic protein formulations. How they may affect protein structure and function is an important question. In this work, through NMR measurements, we show that glycerol and sorbitol (polyols), as well as glucose (sugar), can shorten protein backbone hydrogen bonds. The hydrogen bond shortening is also captured by molecular dynamics simulations, which suggest a hydrogen bond competition mechanism. Specifically, osmolytes weaken hydrogen bonds between the protein and solvent to strengthen those within the protein. Although the hydrogen bond change is small, with the average experimental cross hydrogen bond 3hJNC' coupling of two proteins GB3 and TTHA increased by ~ 0.01 Hz by the three osmolytes (160 g/L), its effect on protein function should not be overlooked. This is exemplified by the PDZ3-peptide binding where several intermolecular hydrogen bonds are formed and osmolytes shift the equilibrium towards the bound state.

Conformation control through concurrent N-H⋯S and N-H⋯O[double bond, length as m-dash]C hydrogen bonding and hyperconjugation effects

In addition to the classical N-H⋯O[double bond, length as m-dash]C non-covalent interaction, less conventional types of hydrogen bonding, such as N-H⋯S, may play a key role in determining the molecular structure. In this work, using theoretical calculations in combination with spectroscopic analysis in both gas phase and solution phase, we demonstrate that both these H-bonding modes exist simultaneously in low-energy conformers of capped derivatives of Attc, a thietane α-amino acid. 6-Membered ring inter-residue N-H⋯S interactions (C6^γ), assisted by hyperconjugation between the thietane ring and the backbone, combine with 5-membered ring intra-residue backbone N-H⋯O[double bond, length as m-dash]C interactions (C5) to provide a C5-C6^γ feature that stabilizes a planar geometry in the monomer unit. Two contiguous C5-C6^γ features in the planar dimer implicate an unprecedented three-centre H-bond of the type C[double bond, length as m-dash]O⋯H(N)⋯SR₂, while the trimer adopts two C5-C6^γ features separated by a Ramachandran α-type backbone configuration. These low-energy conformers are fully characterized in the gas phase and support is presented for their existence in solution state.