Conformational stability of helical peptides containing a thioamide linkage

[structure: see text] Thioxo peptide analogues of the alpha-helical peptide GCN4-p1 were synthesized and evaluated for helicity and oligomeric state. Sedimentation equilibrium and CD measurements indicate that the thioxo peptides fold into parallel alpha-helical coiled coil structures essentially identical to the native structure. This work marks the first incorporation of a thioamide linkage into the backbone of an alpha-helix and demonstrates that a thioamide linkage is compatible with positions within the helix as well as near the C-terminus.     

Theoretical conformational analysis of opiate peptides Leu-Enkephalin (H-Tyr-Gly-Gly-Phe-Leu-OH) and its two thioamide analogs (H-Tyr-Glyψ[CSNH]Gly-Phe-Leu-OH) and (H-Tyr-Gly-Glyψ[CSNH]Phe-Leu-OH)

In order to find informations on the native structure of the Leu-Enkephalin opiate peptide, the parent peptide and its two thioamide analogs (Thio-Gly2)Leu-Enkephalin and (Thio-Gly3)Leu-Enkephalin were studied by the theoretical method PEPSEA. This comparative conformational analysis showed that the active conformation is a β turn structure centered on Gly3 and Phe4. Moreover, this study showed also that the more active analog (Thio-Gly2)Leu-Enk has a lower tendency to adopt this structure. Consequently, its high activity can only be explained by its long lifetime due to its resistance to enzymatic hydrolysis, following the substitution of the amide linkage by the thioamide one. The weakly active analog (Thio-Gly3)Leu-Enk does not adopt this structure and prefers instead a β turn structure centered on Gly2 and Gly3. This study also confirmed the importance of the distances between the Tyr and Phe residues at positions 1 and 4, and that of the terminal Tyrosine N-H group which must be free of any intramolecular hydrogen bond in order to be available in the molecular recognition process.     

Rational design of thioamide peptides as selective inhibitors of cysteine protease cathepsin L

Aberrant levels of cathepsin L (Cts L), a ubiquitously expressed endosomal cysteine protease, have been implicated in many diseases such as cancer and diabetes. Significantly, Cts L has been identified as a potential target for the treatment of COVID-19 due to its recently unveiled critical role in SARS-CoV-2 entry into the host cells. However, there are currently no clinically approved specific inhibitors of Cts L, as it is often challenging to obtain specificity against the many highly homologous cathepsin family cysteine proteases. Peptide-based agents are often promising protease inhibitors as they offer high selectivity and potency, but unfortunately are subject to degradation in vivo. Thioamide substitution, a single-atom O-to-S modification in the peptide backbone, has been shown to improve the proteolytic stability of peptides addressing this issue. Utilizing this approach, we demonstrate herein that good peptidyl substrates can be converted into sub-micromolar inhibitors of Cts L by a single thioamide substitution in the peptide backbone. We have designed and scanned several thioamide stabilized peptide scaffolds, in which one peptide, R^S_1A, was stabilized against proteolysis by all five cathepsins (Cts L, Cts V, Cts K, Cts S, and Cts B) while inhibiting Cts L with >25-fold specificity against the other cathepsins. We further showed that this stabilized R^S_1A peptide could inhibit Cts L in human liver carcinoma lysates (IC₅₀ = 19 μM). Our study demonstrates that one can rationally design a stabilized, specific peptidyl protease inhibitor by strategic placement of a thioamide and reaffirms the place of this single-atom modification in the toolbox of peptide-based rational drug design.

Information on the effects of sidechain and backbone modification on the activity of cathepsin (Cts) L, V, K, S, and B was used to design a thioamide peptide that is inert to all Cts and selectively inhibits Cts L.     

Filaggrin Peptides with β‐Hairpin Structure Bind Rheumatoid Arthritis Antibodies

In the early detection of rheumatoid arthritis (RA) synthetic filaggrin peptides serve as antigens for rheumatoid‐specific autoantibodies (anti‐citrullinated peptide antibody, ACPA) in ELISA tests. In this work we present a peptide that exhibits the binding epitope of ACPA in the form of a stable folding β‐hairpin. The homogeneity of the peptide folding was confirmed by NMR spectroscopy and might lead to the first proposed structure of the antibody‐bound conformation of the epitope.     

Minimalist probes for studying protein dynamics: thioamide quenching of selectively excitable fluorescent amino acids

Fluorescent probe pairs that can be selectively excited in the presence of Trp and Tyr are of great utility in studying conformational changes in proteins. However, the size of these probe pairs can restrict their incorporation to small portions of a protein sequence where their effects on secondary and tertiary structure can be tolerated. Our findings show that a thioamide bond-a single atom substitution of the peptide backbone-can quench fluorophores that are red-shifted from intrinsic protein fluorescence, such as acridone. Using steady-state and fluorescence lifetime measurements, we further demonstrate that this quenching occurs through a dynamic electron-transfer mechanism. In a proof-of-principle experiment, we apply this technique to monitor unfolding in a model peptide system, the villin headpiece HP35 fragment. Thioamide analogues of the natural amino acids can be placed in a variety of locations in a protein sequence, allowing one to make a large number of measurements to model protein folding.     

The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation.

The structural and thermodynamic properties of a 6-residue beta-peptide that was designed to form a hairpin conformation have been studied by NMR spectroscopy and MD simulation in methanol solution. The predicted hairpin would be characterized by a 10-membered hydrogen-bonded turn involving residues 3 and 4, and two extended antiparallel strands. The interproton distances and backbone torsional dihedral angles derived from the NMR experiments at room temperature are in general terms compatible with the hairpin conformation. Two trajectories of system configurations from 100-ns molecular-dynamics simulations of the peptide in solution at 298 and 340 K have been analyzed. In both simulations reversible folding to the hairpin conformation is observed. Interestingly, there is a significant conformational overlap between the unfolded state of the peptide at each of the temperatures. As already observed in previous studies of peptide folding, the unfolded state is composed of a (relatively) small number of predominant conformers and in this case lacks any type of secondary-structure element. The trajectories provide an excellent ground for the interpretation of the NMR-derived data in terms of ensemble averages and distributions as opposed to single-conformation interpretations. From this perspective, a relative population of the hairpin conformation of 20% to 30% would suffice to explain the NMR-derived data. Surprisingly, however, the ensemble of structures from the simulation at 340 K reproduces more accurately the NMR-derived data than the ensemble from the simulation at 298 K, a question that needs further investigation.     

@-Tide-stabilized beta-hairpins

As minimalist versions of beta-structure, two-stranded beta-hairpins are commonly employed as platforms for assessing the interactions that stabilize beta-sheets in proteins. We have found that the presence of a 1,6-dihydro-3(2H)-pyridinone moiety (the "@-unit") as an amino acid replacement at the i - 1 or i + 4 positions relative to a beta-turn strongly stabilizes the hairpin conformation. Hybrids of this type bridge the gap between natural beta-hairpins and unnatural beta-sheets because the @-unit only replaces one residue in a peptide while stabilizing the hairpin conformation to a greater extent than a normal amino acid. In this report, we describe the synthesis of a variety of @-tide-templated hairpins and the NMR and CD characterization of their conformations in both polar and nonpolar solvents.     

Reaction and characterization of thioamide dianions derived from N-benzyl thioamides

Thioamide dianions were generated by the highly efficient reaction of N-benzyl thioamides with 2 equiv of BuLi. Alkylation, allylation, and silylation took place selectively at the carbon atom adjacent to the nitrogen atom of the thioamide dianions. Oxiranes and an aldehyde were also used as electrophiles in the reaction of thioamide dianions to form N-thioacyl 1,3- or 1,2-amino alcohols. The insertion reaction of elemental sulfur to a thioamide dianion and subsequent ethylation afforded a N-thioacyl hemithioaminal. NMR studies on the thioamide mono- and dianions derived from N-benzyl 2-methoxythiobenzamide showed a linear relationship between the chemical shifts of all carbon atoms of thioamide mono- and dianions. The results also suggested that the negative charge at the benzylic carbon atom of the dianion is not fully delocalized. The charge distribution patterns of the dianion are consistent with those of pi polarization.     

Probing the structural and electronic effects to stabilize nonplanar forms of thioamide derivatives: A computational study

Density functional calculations have been performed to examine the stability of nonplanar conformations of thioamide derivatives. Electrostatic, orbital, and ring strain effects were invoked to stabilize the nonplanar conformations of thioamide systems 2 – 7 . Electrostatic interactions helped to achieve the twisted forms of thioamide derivatives; however, pyramidal forms predicted to be the global minimum. Negative hyperconjugative type interactions enhanced the stability of the twisted form 4b when compared with the planar form 4a . The influence of ring strain effect to achieve the twisted form of thioamide was observed with azirine ring. The predictions made with B3LYP/cc‐pVDZ+ level of theory was found to be in good agreement with more accurate CBS‐QB3 method. The solvent calculations performed with polarized continuum solvation model suggest that the relative stabilities of the nonplanar forms of thioamide derivatives are in general similar to the gas phase results. The importance of hydrogen bonding interactions between the solvent molecules and thioamide derivatives was observed toward the enhanced stability of twisted forms using a combination of explicit solvent molecules and continuum model. The natural bond orbital analysis confirmed the participation of n_N → π*_C?S delocalizations in the planar forms and corroborated the earlier reports on larger delocalizations in thioamide systems. Furthermore, the influence of electrostatic and ring strain effects on the amide, natural amides, and selenoamide has also been studied. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011     

Versatile bridging ability of secondary thioamide group for constructing metal cluster based on pincer complex

The crystallization of the SCS pincer platinum complex with secondary thioamide groups, [4-methyl-2,6-bis(anilinothiocarbonyl)-κ²S,S’-phenyl-κC¹]chloroplatinum(II), spontaneously forms a trinuclear cluster complex accompanied by liberating HCl in DMSO/EtOH mixture. The crystal structure of the SNS pincer copper complex with secondary thioamide groups, [2,5-bis(benzylaminothiocarbonyl)-κ²S,S’-pyrrolyl-κN¹]chlorocopper(II), exhibits a dimer structure bearing the bridging thioamide group. The thioamide group acts as the bridging ligand of the multimetallic complexes in different anionic and neutral manners.     

Engineering a β‐Helical d,l‐Peptide for Folding in Polar Media

β Helices—helices formed by alternating d,l ‐peptides and stabilized by β‐sheet hydrogen bonding—are found naturally in only a handful of highly hydrophobic peptides. This paper explores the scope of β‐helical structure by presenting the first design and biophysical characterization of a hydrophilic d,l ‐peptide, 1 , that forms a β helix in methanol. The design of 1 is based on the β‐hairpin/β helix—a new supersecondary that had been characterized previously only for hydrophobic peptides in nonpolar solvents. Incorporating polar residues in 1 provided solubility in methanol, in which the peptide adopts the expected β‐hairpin/β‐helical structure, as evidenced by CD, analytical ultracentrifugation (AUC), NMR spectroscopy, and NMR‐based structure calculations. Upon titration with water (at constant peptide concentration), the structure in methanol ( 1 m ) transitions cooperatively to an extended conformation ( 1 w ) resembling a cyclic β‐hairpin; observation of an isodichroic point in the solvent‐dependent CD spectra indicates that this transition is a two‐state process. In contrast, neither 1 m nor 1 w show cooperative thermal melting; instead, their structures appear intact at temperatures as high as 65 °C; this observation suggests that steric constraint is dominant in stabilizing these structures. Finally, the ¹H NMR CαH spectroscopic resonances of 1 m are downfield‐shifted with respect to random‐coil values, a hitherto unreported property for β helices that appears to be a general feature of these structures. These results show for the first time that an appropriately designed β‐helical peptide can fold stably in a polar solvent; furthermore, the structural and spectroscopic data reported should prove useful in the future design and characterization of water‐soluble β helices.     

Restricting the conformational heterogeneity of RNA by specific incorporation of 8-bromoguanosine

In an effort to reduce the conformational heterogeneity of RNA, the modified nucleobase 8-bromoguanosine (8BrG) was introduced into oligonucleotides having the hairpin tetraloop motif YNMG (Y = U or C and M = C or A). Purine nucleobases with bromine at position eight are known to preferentially adopt the syn conformation as nucleosides. The hairpin tetraloop motif YNMG was chosen as a model system because it has a syn guanosine at position four of the loop that is essential for thermodynamic stability. Thermodynamic and structural characterization of modified oligonucleotides with the hairpin sequences UUCG, CGCG, and CGAG by UV-melting and NMR spectroscopy revealed that 8BrG substitution has a small effect upon the hairpin conformation, while the duplex conformation is strongly destabilized (DeltaDeltaG degrees 37 approximately +4.7 kcal mol-1), thus inhibiting dimerization. These results support a model in which 8BrG substitution shifts the hairpin-duplex equilibrium constant toward the hairpin conformation by destabilizing the duplex. This methodology should be useful for limiting conformational heterogeneity in large RNAs, with potential applications in structural biology and enzymology.     

Solution structure, dimerization, and dynamics of a lipophilic alpha/3(10)-helical, C alpha-methylated peptide. Implications for folding of membrane proteins.

The solution structure and the dimerization behavior of the lipophilic, highly C(alpha)-methylated model peptide, mBrBz-Iva(1)-Val(2)-Iva(3)-(alphaMe)Val(4)-(alphaMe)Phe(5)-(alphaMe)Val(6)-Iva(7)-NHMe, was studied by NMR spectroscopy and molecular dynamics simulations. The conformational analysis resulted in a right-handed 3(10)/alpha-helical equilibrium fast on the NMR time scale with a slight preference for the alpha-helical conformation. The NOESY spectrum showed intermolecular NOEs due to an aggregation of the heptapeptide. In addition, temperature-dependent diffusion measurements were performed to calculate the hydrodynamic radius. All these findings are consistent with an antiparallel side-by-side dimerization. The structure of the dimeric peptide was calculated with a simulated annealing strategy. The lipophilic dimer is held together by favorable van der Waals interactions in the sense of a bulge fitting into a groove. The flexibility of the helical conformations concerning an alpha/3(10)-helical equilibrium is shown in a 3 ns molecular dynamics simulation of the resulting dimeric structure. Both overall helical structures of each monomer and the antiparallel mode of dimerization are stable. However, transitions were seen of several residues from a 3(10)-helical into an alpha-helical conformation and vice versa. Hence, this peptide represents a good model in which two often-discussed aspects of hierarchical transmembrane protein folding are present: i <-- i + 3 and i <-- i + 4 local H-bonding interactions cause a specific molecular shape which is then recognized as attractive by other surrounding structures.     

Effects of turn stability on the kinetics of refolding of a hairpin in a beta-sheet

As part of our continuing study of the effects of the turn sequence on the conformational stability as well as the mechanism of folding of a beta-sheet structure, we have undertaken a parallel investigation of the solution structure, conformational stability, and kinetics of refolding of the beta-sheet VFIVDGOTYTEV(D)PGOKILQ. The latter peptide is an analogue of the original Gellman beta-sheet VFITS(D)PGKTYTEV(D)PGOKILQ, wherein the TS(D)PGK turn sequence in the first hairpin has been replaced by VDGO. Thermodynamics studies revealed comparable conformational stability of the two peptides. However, unlike the Gellman peptide, which showed extremely rapid refolding of the first hairpin, early kinetic events associated with the refolding of the corresponding hairpin in the VDGO mutant were found to be significantly slower. A detailed study of the conformation of the modified peptide suggested that hydrophobic interactions might be contributing to its stability. Accordingly, we surmise that the early kinetic events are sensitive to whether the formation of the hairpin is nucleated at the turn or by sequestering of the hydrophobic residues across the strand, before structural rearrangements to produce the nativelike topology. Nucleation of the hairpin at the turn is expected to be intrinsically rapid for a strong turn. However, if the process must involve collapse of hydrophobic side chains, the nucleation should be slower as solvent molecules must be displaced to sequester the hydrophobic residues. These findings reflect the contribution of different forces toward nucleation of hairpins in the mechanism of folding of beta-sheets.     

Metal ion complexation by acyclic polyethers with lipophilic amide, thioamide, and amine end groups

A series of acyclic polyethers with lipophilic amide, thioamide, and amine end groups was synthesized. Metal ion transport across bulk liquid membranes and measurement of thermodynamic parameters for ligand-metal ion complexation by titration calorimetry show strong selectivity for complexation of lead ion over other metal ion species for the diamide ligand. Lead ion complexation by the acyclic polyether diamide involves the amide oxygens and silver ion coordination by a dithioamide analog involves the thioamide sulfurs. With a proper length of the ethereal linkage, the ligand wraps around the metal ion in a pseudocyclic fashion.     

Kinetic resolution of racemic alcohols using thioamide modified 1-methyl-histidine methyl ester

The low-molecular-weight and easily prepared N-thiobenzoyl 1-methyl-histidine methyl ester 3k was utilized to efficiently catalyze the kinetic resolution of racemic secondary alcohols. Comparison of the conformations of amide catalyst 3c and thioamide catalyst 3k was made to understand the origin of the improvement of the enantioselectivity by thioamide modification.     

A 310-helix single turn enforced by crosslinking of lysines with 1,1′-ferrocenedicarboxylic acid

Prior work has shown that covalently linking the side chains of amino acids in the i and i+3 and i and i+4 positions in a peptide will enforce a helical conformation. In this work the ability of an organometallic entity to enforce a helical conformation in a peptide was explored. The tetrapeptide Boc-Lys-Ala-Val-Lys-NHCH₃ was prepared, then reacted with 1,1′-ferrocenedicarboxylic acid chloride. Reaction of the lysine side chain amines with the diacid chloride resulted in a metallacyclicpeptide (1) in which the two lysines are crosslinked via the ferrocene. The solution conformation of the metallacyclicpeptide (1) was studied using CD and NMR spectroscopy. The NMR methods employed were Karplus analysis of coupling constants, chemical shift changes of NH protons and ROESY data. The results show that the metallacyclicpeptide (1) adopts a single turn of the 3₁₀-helix conformation.     

New binuclear Pd(II) thioamide complexes for the Heck reaction of aryl bromides

The reaction of binucleating thioamide ligands (L₁-L₃) with [PdCl₂(PPh₃)₂] in 1:2 molar ratio in methanol medium afforded a series of binuclear palladium(II) complexes. The synthesized complexes were characterized by elemental analysis and ¹H NMR. The molecular structure of one of the complexes was established by X-ray diffraction method. The binuclear palladium(II) thioamide complex has been shown to be an active catalyst for the Heck reaction of aryl bromides with alkenes.     

Probing the complete folding trajectory of a DNA hairpin using dual beam fluorescence fluctuation spectroscopy

The conformational fluctuations of dye-quencher labeled DNA hairpin molecules in aqueous solution were investigated using dual probe beam fluorescence fluctuation spectroscopy. The measurements revealed the flow and diffusion times of the DNA molecules through two spatially offset optical probe regions, the absolute and relative concentrations of each conformational substate of the DNA, and the kinetics of the DNA hairpin folding and unfolding reactions in the 1 micros to 10 ms time range. A DNA hairpin containing a 21-nucleotide polythymine loop and a 4-base pair stem exhibited double exponential relaxation kinetics, with time constants of 84 and 393 micros. This confirms that folding and melting of the DNA hairpin structure is not a two state process but proceeds by way of metastable intermediate states. The fast time constant corresponds to formation and unfolding of an intermediate, and the slow time constant is due to formation and disruption of the fully base-paired stem. This is consistent with a previous study of a similar DNA hairpin with a 5-base pair stem, in which the fast reaction was attributed to the fluctuations of an intermediate DNA conformation [J. Am. Chem. Soc. 2006, 128, 1240-1249]. In that case, reactions involving the native conformation could not be observed directly due to the limited observation time range of the fluorescence correlation spectroscopy experiment. The intermediate states of the DNA hairpins are suggested to be due to a collapsed ensemble of folded hairpins containing various partially folded or misfolded conformations.     

RNA architecture dictates the conformations of a bound peptide.

BACKGROUND: The biological function of several viral and bacteriophage proteins, and their arginine-rich subdomains, involves RNA-mediated interactions. It has been shown recently that bound peptides adopt either beta-hairpin or alpha-helical conformations in viral and phage peptide-RNA complexes. We have compared the structures of the arginine-rich peptide domain of HIV-1 Rev bound to two RNA aptamers to determine whether RNA architecture can dictate the conformations of a bound peptide. RESULTS: The core-binding segment of the HIV-1 Rev peptide class II RNA aptamer complex spans the two-base bulge and hairpin loop of the bound RNA and the carboxy-terminal segment of the bound peptide. The bound peptide is anchored in place by backbone and sidechain intermolecular hydrogen bonding and van der Waals stacking interactions. One of the bulge bases participates in U*(A*U) base triple formation, whereas the other is looped out and flaps over the bound peptide in the complex. The seven-residue hairpin loop is closed by a sheared G*A mismatch pair with several pyrimidines looped out of the hairpin fold. CONCLUSIONS: Our structural studies establish that RNA architecture dictates whether the same HIV-1 Rev peptide folds into an extended or alpha-helical conformation on complex formation. Arginine-rich peptides can therefore adapt distinct secondary folds to complement the tertiary folds of their RNA targets. This contrasts with protein-RNA complexes in which elements of RNA secondary structure adapt to fit within the tertiary folds of their protein targets.