Using thioamides to site-specifically interrogate the dynamics of hydrogen bond formation in β-sheet folding

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

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

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

Thioimidate Solutions to Thioamide Problems during Thionopeptide Deprotection

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

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.

Thioamide Analogues of MHC I Antigen Peptides

Short, synthetic peptides that are displayed by major histocompatibility complex I (MHC I) can stimulate CD8 T cells in vivo to destroy virus-infected or cancer cells. The development of such peptides as vaccines that provide protective immunity, however, is limited by rapid proteolytic degradation. Introduction of unnatural amino acid residues can suppress MHC I antigen proteolysis, but the modified peptides typically display lower affinity for MHC I and/or diminished ability to activate CD8 T cells relative to native antigen. Here, we report a new strategy for modifying MHC I antigens to enhance resistance to proteolysis while preserving MHC I affinity and T cell activation properties. This approach, replacing backbone amide groups with thioamides, was evaluated in two well-characterized antigens presented by HLA-A2, a common human MHC I. For each antigen, singly modified thioamide analogues retained affinity for HLA-A2 and activated T cells specific for the native antigen, as measured via interferon-γ secretion. In each system, we identified a highly potent triply substituted thioamide antigen ("thio-antigen") that displayed substantial resistance to proteolytic cleavage. Collectively, our results suggest that thio-antigens may represent a general and readily accessible source of potent vaccine candidates that resist degradation.

Structural impact of thioamide incorporation into a β-hairpin

The thioamide is a naturally-occurring single atom substitution of the canonical amide bond. The exchange of oxygen to sulfur alters the amide's physical and chemical characteristics, thereby expanding its functionality. Incorporation of thioamides in prevalent secondary structures has demonstrated that they can either have stabilizing, destabilizing, or neutral effects. We performed a systematic investigation of the structural impact of thioamide incorporation in a β-hairpin scaffold with nuclear magnetic resonance (NMR). Thioamides as hydrogen bond donors did not increase the foldedness of the more stable "YKL" variant of this scaffold. In the less stable "HPT" variant of the scaffold, the thioamide could be stabilizing as a hydrogen bond donor and destabilizing as a hydrogen bond acceptor, but the extent of the perturbation depended upon the position of incorporation. To better understand these effects we performed structural modelling of the macrocyclic folded HPT variants. Finally, we compare the thioamide effects that we observe to previous studies of both side-chain and backbone perturbations to this β-hairpin scaffold to provide context for our observations.

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

Hydrogen Bond and Geometry Effects of Thioamide Backbone Modifications

Thioamide substitution of backbone peptide bonds can probe interactions along the main chain of proteins. Despite theoretical predictions of the enhanced hydrogen bonding propensities of thioamides, previous studies often do not consider the geometric constraints imposed by folded peptide secondary structure. This work addresses drawbacks in previous studies that ignored the geometry dependence and local dielectric properties of thioamide hydrogen bonding and identifies cases where thioamides may be either stronger or weaker hydrogen-bonding partners than amides.

Comparing C2=O and C2=S Barbiturates: Different Hydrogen-Bonding Patterns of Thiobarbiturates in Solution and the Solid State

Carbonyl-centered hydrogen bonds with various strength and geometries are often exploited in materials to embed dynamic and adaptive properties, with the use of thiocarbonyl groups as hydrogen-bonding acceptors remaining only scarcely investigated. We herein report a comparative study of C2=O and C2=S barbiturates in view of their differing hydrogen bonds, using the 5,5-disubstituted barbiturate B and the thiobarbiturate TB as model compounds. Owing to the different hydrogen-bonding strength and geometries of C2=O vs. C2=S, we postulate the formation of different hydrogen-bonding patterns in C2=S in comparison to the C2=O in conventional barbiturates. To study differences in their association in solution, we conducted concentration- and temperature-dependent NMR experiments to compare their association constants, Gibbs free energy of association ∆Gassn., and the coalescence behavior of the N-H‧‧‧S=C bonded assemblies. In Langmuir films, the introduction of C2=S suppressed 2D crystallization when comparing B and TB using Brewster angle microscopy, also revealing a significant deviation in morphology. When embedded into a hydrophobic polymer such as polyisobutylene, a largely different rheological behavior was observed for the barbiturate-bearing PB compared to the thiobarbiturate-bearing PTB polymers, indicative of a stronger hydrogen bonding in the thioanalogue PTB. We therefore prove that H-bonds, when affixed to a polymer, here the thiobarbiturate moieties in PTB, can reinforce the nonpolar PIB matrix even better, thus indicating the formation of stronger H-bonds among the thiobarbiturates in polymers in contrast to the effects observed in solution.

Folding and Unfolding of the Tryptophan Zipper in the Presence of Two Thioamide Substitutions

We studied the stability and folding and unfolding kinetics of the tryptophan zipper, containing different double thioamide subsitutions. Conformation change was triggered by photoisomerization of an integrated AMPP photoswitch in the turn region of the hairpin, and transient spectra were recorded in the deep UV and the mid-IR, covering the time window of the (un)folding transition from picoseconds to tens of microseconds. Thio-substitution of inward-pointing backbone carbonyls was found to strongly destabilize the β-hairpin structures, whereas molecules with two outward pointing thio-carbonyls showed similar or enhanced stability with respect to the unsubstituted sequence, which we attribute to stronger interstrand hydrogen bonding. Thiolation of the two Trp residues closest to the turn can even prevent the opening of the hairpin after cis-trans isomerization of the switch. The circular dichroism due to the two thioamide ππ* transitions is spectrally well-separated from the aromatic tryptophan signal. It changes upon photoswitching, reflecting a local change in coupling and geometry.

Opportunities and challenges in the synthesis of thioamidated peptides

Chemical modifications of peptides hold great promise for modulating their pharmacological properties. In the last few decades amide to thioamide substitution has been widely explored to modulate the conformation, non-covalent interactions, and proteolytic stability of peptides. Despite widespread utilization, there are some potential limitations including epimerization and degradation under basic and acidic conditions, respectively. In this chapter, we present the synthetic method to build thio-precursors, their site-specific incorporation onto a growing peptide chain, and troubleshooting during the elongation of thioamidated peptides. This highly efficient, rapid, and robust method can be used for positional scanning of the thioamide bond.

Enantioselective, Catalytic Multicomponent Synthesis of Homoallylic Amines Enabled by Hydrogen-Bonding and Dispersive Interactions

We report a one-step catalytic, enantioselective method for the preparation of homoallylic N-Boc amines directly from acetals. Reactive iminium ion intermediates are generated in situ through the combination of an acetal, a chiral thiourea catalyst, trialkylsilyl triflate, and N-Boc carbamate and are subsequently trapped by a variety of allylsilane nucleophiles. No homoallylic ether byproducts are detected, consistent with allylation of the iminium intermediate being highly favored over allylation of the intermediate oxocarbenium ion. Acetals derived from aromatic aldehydes possessing a variety of functional groups and substitution patterns yield homoallylic amines with excellent levels of enantiomeric enrichment. Experimental and computational data are consistent with an anchoring hydrogen-bond interaction between the protioiminium ion and the amide of the catalyst in the enantiodetermining transition state, and with stereodifferentiation achieved through specific noncovalent interactions (NCIs) with the catalyst pyrenyl moiety. Evidence is provided that the key NCI in the major pathway is a π-stacking interaction, contrasting with the cation-π interactions invoked in previously studied reactions promoted by the same family of aryl-pyrrolidino-H-bond-donor catalysts.

Nickel-catalyzed C-alkylation of thioamide, amides and esters by primary alcohols through a hydrogen autotransfer strategy

A simple catalyst of Ni(OAc)2 and P(t-Bu)3 enables selective C-alkylation of thioacetamides and primary acetamides with alcohols for the first time. Monoalkylation of thioamides, amides and t-butyl esters occurs in excellent yields (>95%). Mechanistic studies reveal that the reaction proceeds via a hydrogen autotransfer pathway.

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.

Convenient synthesis of thioamidated peptides and proteins

The unique physicochemical properties of a thioamide bond, which is an ideal isostere of an amide bond, have not been fully exploited because of the tedious synthesis of thionated amino acid building blocks. Here, we report a purification-free and highly efficient synthesis of thiobenzotriazolides of Fmoc-protected and orthogonally protected 20 naturally occurring amino acids including asparagine, glutamine, and histidine. The near-quantitative conversion to the respective thioamidated peptides on solid support demonstrates the robustness of the synthetic route. Furthermore, the unaltered incorporation efficiency of thiobenzotriazolides from their stock solution till 48 h suggests their compatibility toward automated peptide synthesis. Finally, utilizing an optimized cocktail of 2% DBU + 5% piperazine for fast Fmoc-deprotection, we report the synthesis of a thioamidated Pin1 WW domain and thioamidated GB1 directly on solid support.

Broad-Band Ultraviolet CD Spectroscopy of Ultrafast Peptide Backbone Conformational Dynamics

The far-UV spectral window widely used for the conformational analysis of biomolecules is not easily covered with broad-band lasers. This has made it difficult to use circular dichroism (CD) spectroscopy to directly follow fast structure changes. By combining transient CD spectroscopy in the deep-UV with thioamide substitution, we demonstrate a method to overcome this difficulty. We investigated a dipeptide whose two carbonyl oxygen atoms were replaced by sulfur, red-shifting the strong lowest-lying ππ* transitions into the more accessible 250-370 nm spectral window. Coupling of the two thioamide units cannot be resolved by achiral 2D-UV spectroscopy, but it gives rise to a pronounced bisignate CD spectrum. The transient CD spectra reveal weakening of this coupling in the electronically excited state, where conformational constraints are released. Our results show that direct local probing of fast backbone conformational change via CD spectroscopy is possible in combination with site-selective thio substitution in peptides and proteins.

Biosynthesis and Chemical Applications of Thioamides

Thioamidation as a posttranslational modification is exceptionally rare, with only a few reported natural products and exactly one known protein example (methyl-coenzyme M reductase from methane-metabolizing archaea). Recently, there has been significant progress in elucidating the biosynthesis and function of several thioamide-containing natural compounds. Separate developments in the chemical installation of thioamides into peptides and proteins have enabled cell biology and biophysical studies to advance the current understanding of natural thioamides. This review highlights the various strategies used by Nature to install thioamides in peptidic scaffolds and the potential functions of this rare but important modification. We also discuss synthetic methods used for the site-selective incorporation of thioamides into polypeptides with a brief discussion of the physicochemical implications. This account will serve as a foundation for the further study of thioamides in natural products and their various applications.

Protein stabilization by tuning the steric restraint at the reverse turn

The incorporation of pseudoallylic strain by N-methylation at the solvent exposed loop in proteins leads to a stark increase in their thermodynamic stability that can be tuned by altering the amino acid composition.

Reverse turns are solvent-exposed motifs in proteins that are crucial in nucleating β-sheets and drive the protein folding. The solvent-exposed nature makes reverse turns more amenable to chemical modifications than α-helices or β-sheets towards modulating the stability of re-engineered proteins. Here, we utilize van der Waals repulsive forces in tuning the steric restraint at the reverse turn. The steric restraint induced upon N-methylation of the i+1–i+2 amide bond at the reverse turn results in well-folded and stable β-sheets in aqueous solution at room temperature. The developed superactive turn inducing motif is tolerant to a wide variety of functional groups present on coded amino acids making the designed turn fully compatible with bioactive loops in proteins. We demonstrate that the steric restraint and the functional groups at the reverse turn act in synergy to modulate the folding of re-engineered β-sheets. Introduction of the turn motifs onto a three-stranded β-sheet protein, Pin 1 WW domain, resulted in various analogs showing a cooperative two-state transition with thermal stability (T_M) ranging from 62 °C to 82 °C. Despite modulating the stability of Pin 1 variants by ∼2.8 kcal mol^–1 (ΔΔG_f), the native fold in all the protein variants was found to be unperturbed. This structural stability is brought about by conformational preorganization at the engineered reverse turn that results in strong intramolecular hydrogen bonds along the three dimensional structure of the protein. Thus, this simple loop engineering strategy via two amino acid substitution provides us a “toolkit” to modulate the stability of β-sheet containing peptides and proteins in aqueous solution that will greatly expand the scope of de novo protein and foldamer design.

Thioamide substitution to probe the hydroxyproline recognition of VHL ligands

Thioamide substitution influences hydrogen bond and n → π^∗ interactions involved in the conformational stability of protein secondary structures and oligopeptides. Hydroxyproline is the key recognition element of small molecules targeting the von Hippel-Lindau (VHL) E3 ligase, which are of interest as probes of hypoxia signaling and ligands for PROTAC conjugation. We hypothesized that VHL ligands could be a privileged model system to evaluate the contribution of these interactions to protein:ligand complex formation. Herein we report the synthesis of VHL ligands bearing thioamide substitutions at the central hydroxyproline moiety, and characterize their binding by fluorescence polarization, isothermal titration calorimetry, X-ray crystallography and molecular modeling. In spite of a conserved binding mode, the substitution pattern had a pronounced impact on the ligand affinities. Together the results underscore the role of hydrogen bond and n → π^∗ interactions in fine tuning hydroxyproline recognition by VHL.

Dithioamide substitutions in proteins: effects on thermostability, peptide binding, and fluorescence quenching in calmodulin

Thioamide substitutions in the backbones of proteins can modulate their structure and thermostability, or serve as spectroscopic probes in fluorescence quenching experiments. Using native chemical ligation, we have produced the first examples of a protein (calmodulin) containing two thioamides. Dithioamide variants were made to explore the effects of combining stabilizing, neutral, and destabilizing single thioamide substitutions. One of the dithioamide calmodulin variants exhibited stabilization greater than any monothioamide variant, although the effect could not easily be anticipated from the results of single substitutions. Each of the calmodulin variants retained the ability to bind a target peptide, and the dithioamide proteins exhibited an increase in fluorescence quenching of tryptophan relative to their single thioamide counterparts. These results show that multiply thioamidated proteins can be synthesized, and that properly placed thioamides can be used to increase protein thermostability or enhance fluorecsence quenching in peptide binding experiments.

Electronic interactions of i, i + 1 dithioamides: increased fluorescence quenching and evidence for n-to-π* interactions

Thioamide residues can be effective, minimally-perturbing fluorescence quenching probes for studying protein folding and proteolysis. In order to increase the level of quenching, we have here explored the use of adjacent dithioamides. We have found that they are more effective fluorescence quenchers, as expected, but we have also observed unexpected changes in the thioamide absorption spectra that may arise from n-to-π* interactions of the thiocarbonyls. We have made use of the increased quenching to improve the fluorescence turn-on of thioamide protease sensors.

Increasing the bioactive space of peptide macrocycles by thioamide substitution

We show that substituting a single atom, O to S (amide to thioamide), in a peptide bond results in global restriction of the conformational flexibility in peptide macrocycles with minimal perturbation of the parent conformation. The van der Waals interactions between the C 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 S group and the surrounding atoms are the major driving force in inducing the conformational restriction, resulting in well-defined structures of these cyclic peptides with static 3-D presentation of the pharmacophores. Utilizing this property of thioamides, we report the development of a superactive antagonist of pro-angiogenic αvβ3, αvβ5 and α5β1 integrins, which are responsible for cancer cell proliferation and survival. Using simple thio-scanning and spatial screening of a non-efficacious and conformationally flexible cyclic peptide, we could achieve a more than 10⁵ fold enhancement in its efficacy in cellulo via a single O to S substitution. The developed peptide shows better efficacy in inhibiting the pro-angiogenic integrins than the drug candidate cilengitide, with a significantly enhanced serum half-life of 36 h compared to that of cilengitide (12 h). The long shelf-life, absence of non-specific toxicity and resistance to degradation of the thioamidated macrocyclic peptides in human serum suggest the promise of thioamides in markedly improving the affinity, efficacy and pharmacology of peptide macrocycles.

Thioamide Substitution Selectively Modulates Proteolysis and Receptor Activity of Therapeutic Peptide Hormones

Peptide hormones are attractive as injectable therapeutics and imaging agents, but they often require extensive modification by mutagenesis and/or chemical synthesis to prevent rapid in vivo degradation. Alternatively, the single-atom, O-to-S modification of peptide backbone thioamidation has the potential to selectively perturb interactions with proteases while preserving interactions with other proteins, such as target receptors. Here, we use the validated diabetes therapeutic, glucagon-like peptide-1 (GLP-1), and the target of clinical investigation, gastric inhibitory polypeptide (GIP), as proof-of-principle peptides to demonstrate the value of thioamide substitution. In GLP-1 and GIP, a single thioamide near the scissile bond renders these peptides up to 750-fold more stable than the corresponding oxopeptides toward cleavage by dipeptidyl peptidase 4, the principal regulator of their in vivo stability. These stabilized analogues are nearly equipotent with their parent peptide in cyclic AMP activation assays, but the GLP-1 thiopeptides have much lower β-arrestin potency, making them novel agonists with altered signaling bias. Initial tests show that a thioamide GLP-1 analogue is biologically active in rats, with an in vivo potency for glycemic control surpassing that of native GLP-1. Taken together, these experiments demonstrate the potential for thioamides to modulate specific protein interactions to increase proteolytic stability or tune activation of different signaling pathways.

Thioamide, a Hydrogen Bond Acceptor in Proteins and Nucleic Acids

Thioamides are used as potential surrogates of amides to study the structure and dynamics of proteins and nucleic acids. However, incorporation of thioamides in biomolecules leads to changes in their structures and conformations mostly attributed to the strength of the amide-N-H···S═C hydrogen bond. In most cases, it is considered weak owing to the small electronegativity of sulfur, and in some cases, it is as strong as conventional H-bonds. Herein, adopting PDB structure analysis, NMR spectroscopy, and quantum chemistry calculations, we have shown that thioamides in a geometrical and structural constraint-free environment are capable of forming strong H-bonds like their amide counterparts. These studies also enabled us to determine the amide-N-H···S═C H-bond enthalpy (ΔH) very precisely. The estimated ΔH for the amide-N-H···S═C H-bond is ∼-30 kJ/mol, which suggests that the amide-N-H···S═C H-bond is a strong H-bond and merits its inclusion in computational force fields for biomolecular structure simulations to explore the role of amide-N-H···S═C H-bonds in nucleobase pairing and protein folding.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.

The effects of thioamide backbone substitution on protein stability: a study in α-helical, β-sheet, and polyproline II helical contexts

Thioamides are single atom substitutions of the peptide bond that serve as versatile probes of protein structure. Effective use of thioamides requires a robust understanding of the impact that the substitution has on a protein of interest. However, the thermodynamic effects of thioamide incorporation have only been studied in small structural motifs, and their influence on secondary structure in the context of full-length proteins is not known. Here we describe a comprehensive survey of thioamide substitutions in three benchmark protein systems (calmodulin, the B1 domain of protein G, and collagen) featuring the most prevalent secondary structure motifs: α-helix, β-sheet, and polyproline type II helix. We find that in most cases, effects on thermostability can be understood in terms of the positioning and local environment of the thioamide relative to proximal structural elements and hydrogen bonding networks. These observations set the stage for the rational design of thioamide substituted proteins with predictable stabilities.

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

Effects of isosteric substitutions on the conformational preference and cis–trans isomerization of proline-containing peptides

Isosteric substitutions of the peptide CO group by CS and CSe groups increased thetranspopulation and rotational barrier to the prolylcis–transisomerization of proline-containing peptides.

Molecular dynamics simulations elucidate conformational selection and induced fit mechanisms in the binding of PD-1 and PD-L1

PD-L1 binds PD-1 through a complex mechanism including both the conformational selection and induced fit pathways.

Thioxylated cyclosporin A for studying protein-drug interactions

Single atom substitution of cyclosporin A (CsA) through thioxylation has been used to study the structure-activity relationship of the immunosuppressive complex, involving the CsA receptor protein cyclophilin 18 (Cyp18) and the immunological target protein phosphatase calcineurin (CaN), illustrating the contributions of peptide backbone in protein-drug interaction. Moreover, the subtle difference between thioxylation positions in CsA has led to a remarkable change in the quenching effect on Cyp18 intrinsic fluorescence. Using the thioxylated compound Cs7 as an isosteric derivative of CsA in competition assay, the experiment has led to the determination of k_off value in solution. Whereas the conformational heterogeneity of CsA has been found to be associated with its two-phase binding kinetics to Cyp18, the dissociation rate of CsA from complex is independent from the initial ligand structure.

A prevalent intraresidue hydrogen bond stabilizes proteins

Current limitations in de novo protein structure prediction and design suggest an incomplete understanding of the interactions that govern protein folding. Here we demonstrate that previously unappreciated hydrogen bonds occur within proteins between the amide proton and carbonyl oxygen of the same residue. Quantum calculations, infrared spectroscopy, and nuclear magnetic resonance spectroscopy show that these interactions share hallmark features of canonical hydrogen bonds. Biophysical analyses demonstrate that selective attenuation or enhancement of these C5 hydrogen bonds affects the stability of synthetic β-sheets. These interactions are common, affecting approximately 5% of all residues and 94% of proteins, and their cumulative impact provides several kcal/mol of conformational stability to a typical protein. C5 hydrogen bonds stabilize, especially, the flat β-sheets of the amyloid state, which is linked with Alzheimer’s disease and other neurodegenerative disorders. Inclusion of these interactions in computational force fields would improve models of protein folding, function, and dysfunction.

Minimalist IR and fluorescence probes of protein function

Spectroscopic studies of small proteins and peptides, especially those requiring fine spatial and/or temporal resolution, demand synthetic probes that confer the minimal possible steric and functional change on the native properties. Here we review the recent progress in development of minimally disruptive probes for fluorescence and infrared spectroscopies, as well as the methods to efficiently incorporate them into proteins. Advances in spectroscopy on the one hand result in high specialization of synthetic probes for a particular purpose, but on the other hand allow for the same probes be used for different techniques to gather complementary biochemical information.

Engineering β-sheets employing N-methylated heterochiral amino acids

There is a lack of functional group diversity in the reverse turn motifs nucleating a β-sheet conformation in designed peptides, proteins and foldamers. The majority of these sequences consist of d-Pro-l-Pro, d-Pro-Gly or Asn-Gly as the turn inducing motif restricting their biological application and physicochemical modulation. In this report, for the first time we elucidate that N-methylation of heterochiral amino acids in linear peptides nucleates β-sheet conformation without the necessity of having a ring or covalent constraint at the reverse turn. Our results show that d-Pro can be conveniently substituted by any other N-methylated d-amino acid followed by an N-methylated l-amino acid or sarcosine to adopt a βII' turn inducing the β-sheet folding. Furthermore, we reveal that a single amino acid either at the i + 1 or i + 2 site of the reverse turn can modulate the right-handed twist, which eventually dictates the extent of the foldedness of the β-hairpin.

Chemoselective modifications for the traceless ligation of thioamide-containing peptides and proteins

Thioamides are single-atom substitutions of canonical amide bonds, and have been proven to be versatile and minimally perturbing probes in protein folding studies. Previously, our group showed that thioamides can be incorporated into proteins by native chemical ligation (NCL) with Cys as a ligation handle. In this study, we report the expansion of this strategy into non-Cys ligation sites, utilizing radical initiated desulfurization to "erase" the side chain thiol after ligation. The reaction exhibited high chemoselectivity against thioamides, which can be further enhanced with thioacetamide as a sacrificial scavenger. As a proof-of-concept example, we demonstrated the incorporation of a thioamide probe into a 56 amino acid protein, the B1 domain of Protein G (GB1). Finally, we showed that the method can be extended to β-thiol amino acid analogs and selenocysteine.

Infrared and Fluorescence Assessment of Protein Dynamics: From Folding to Function

While folding or performing functions, a protein can sample a rich set of conformational space. However, experimentally capturing all of the important motions with sufficient detail to allow a mechanistic description of their dynamics is nontrivial since such conformational events often occur over a wide range of time and length scales. Therefore, many methods have been employed to assess protein conformational dynamics, and depending on the nature of the conformational transition in question, some may be more advantageous than others. Herein, we describe our recent efforts, and also those of others, wherever appropriate, to use infrared- and fluorescence-based techniques to interrogate protein folding and functional dynamics. Specifically, we focus on discussing how to use extrinsic spectroscopic probes to enhance the structural resolution of these techniques and how to exploit various cross-linking strategies to acquire dynamic and mechanistic information that was previously difficult to attain.

Semi-synthesis of thioamide containing proteins

Our laboratory has shown that the thioamide, a single atom O-to-S substitution, can be a versatile fluorescence quenching probe that is minimally-perturbing when placed at many locations in a protein sequence. In order to make these and other thioamide experiments applicable to full-sized proteins, we have developed methods for incorporating thioamides by generating thiopeptide fragments through solid phase synthesis and ligating them to protein fragments expressed in E. coli. To install donor fluorophores, we have adapted unnatural amino acid mutagenesis methods, including the generation of new tRNA synthetases for the incorporation of small, intrinsically fluorescent amino acids. We have used a combination of these two methods, as well as chemoenzymatic protein modification, to efficiently install sidechain and backbone modifications to generate proteins labeled with fluorophore/thioamide pairs.

Modeling Vibrational Spectra of Ester Carbonyl Stretch in Water and DMSO Based on Molecular Dynamics Simulation

On the basis of molecular dynamics simulation, we model the ester carbonyl stretch FTIR signals of methyl acetate in D2O and DMSO. An ab initio map is constructed at the B3LYP/6-311++G** level to relate the carbonyl stretch frequency to the external electric field. Using this map, fluctuating Hamiltonian of the carbonyl stretch is constructed from the MD simulation trajectory. The IR spectra calculated based on this Hamiltonian are found to be in good agreement with the experiment. For methyl acetate in D2O, hydrogen bonding on alkoxy oxygen causes a blue shift of frequency, while that on carbonyl oxygen causes a red shift. Two peaks observed in FTIR signals originate from the balance of these two effects. Furthermore, in both D2O and DMSO solutions, correlations are found between the instantaneous electric field on C═O and the frequencies. Broader line width of the signal in D2O suggests a more inhomogeneous electric field distribution due to the complicated hydrogen-bonding environment.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs

Residue-level unfolding of two helix-turn-helix proteins--one naturally occurring and one de novo designed--is reconstructed from multiple sets of site-specific (13)C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa-Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako-Saitô-Muñoz-Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and (13)C-amide I' bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for "experimental" reaction coordinates--namely, the degree of local folding as sensed by site-specific (13)C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.

Environmentally responsive histidine-carboxylate zipper formation between proteins and nanoparticles

Interfacing synthetic materials with biomacromolecules provides new systems for biological applications. We report the creation of a reversible multivalent supramolecular "zipper" recognition motif between gold nanoparticles and proteins. In this assembly, carboxylate-functionalized nanoparticles interact strongly with oligohistidine tags. This interaction can be tuned through His-tag length, and offers unique binding profiles based on the pH and electrolyte concentration of the medium.

On the use of thioamides as fluorescence quenching probes for tracking protein folding and stability

Our laboratory has developed thioamide analogs of the natural amino acids as minimally-perturbing fluorescence quenching probes that can be placed at many locations in a protein sequence. We have shown that the mechanism of quenching can be either Förster resonance energy transfer (FRET) or photoinduced electron transfer (PET), depending on the identity of the donor fluorophore. Furthermore, we have shown that one can use a combination of semi-synthetic methods to label full-sized proteins with fluorophore-thioamide pairs. These probes can be used to study protein-protein interactions, protein folding or misfolding, and proteolysis.