Exploring helix structures of γ-peptides based on 2-(aminomethyl)cyclopentanecarboxylic acid.

Conformational search and density functional theory calculations were performed to explore the preferences of helical structures for chiro-specific oligo-γ-peptides of 2-(aminomethyl)cyclopentanecarboxylic acid (γAmc5) with a cyclopentyl constraint on the Cα-Cß bond in solution. The dimer and tetramer of γAmc5 (1) with homochiral (1S, 2S) configurations exhibited a strong preference for the 9-membered helix foldamer in solution, except for the tetramer in water. However, the oligomers of γAmc5 (1) longer than tetramer preferentially adopted a right-handed (P)-2.614-helix (H1-14) as the peptide sequence becomes longer and as solvent polarity increases. The high stabilities for H1-14 foldamers of γAmc5 (1) in solution were ascribed to the favored solvation free energies. The calculated mean backbone torsion angles for H1-14 helix foldamers of γAmc5 (1) were similar to those calculated for oligomers of other γ-residues with cyclopentane or cyclohexane rings. However, the substitution of cyclopentane constraints on the Cα-Cß bond of the γAmc5 (1) residue resulted in different conformational preferences and/or handedness of helix foldamers. In particular, the pyrrolidine-substituted analogs of the H1-14 foldamers of γAmc5 (1) with adjacent amine diads substituted at a proximal distance are expected to be potential catalysts for the crossed aldol condensation in nonpolar and polar solvents.

Prospect of Thiazole-based γ-Peptide Foldamers in Enamine Catalysis: Exploration of the Nitro-Michael Addition.

As three-dimensional folding is prerequisite to biopolymer activity, complex functions may also be achieved through foldamer science. Because of the diversity of sizes, shapes and folding available with synthetic monomers, foldamer frameworks enable a numerous opportunities for designing new generations of catalysts. We herein demonstrate that heterocyclic γ-peptide scaffolds represent a versatile platform for enamine catalysis. One central feature was to determine how the catalytic activity and the transfer of chiral information might be under the control of the conformational behaviours of the oligomer.

C9/12 Ribbon-Like Structures in Hybrid Peptides Alternating α- and Thiazole-Based γ-Amino Acids.

According to their restricted conformational freedom, heterocyclic γ-amino acids are usually considered to be related to Z-vinylogous γ-amino acids. In this context, oligomers alternating α-amino acids and thiazole-based γ-amino acids (ATCs) were expected to fold into a canonical 12-helical shape as described for α/γ-hybrid peptides composed of cis-α/ß-unsaturated γ-amino acids. However, through a combination of X-ray crystallography, NMR spectroscopy, FTIR experiments, and DFT calculations, it was determined that the folding behavior of ATC-containing hybrid peptides is much more complex. The homochiral α/(S)-ATC sequences were unable to adopt a stable conformation, whereas the heterochiral α/(R)-ATC peptides displayed novel ribbon structures stabilized by unusual C9/12 -bifurcated hydrogen bonds. These ribbon structures could be considered as a succession of pre-organized γ/α dipeptides and may provide the basis for designing original α-helix mimics.

12/10-Helical ß-Peptide with Dynamic Folding Propensity: Coexistence of Right- and Left-Handed Helices in an Enantiomeric Foldamer.

We present the first examples of atomic-resolution crystal data for the ß-peptide 12/10-helix from oligomers of cis-2-aminocyclohexane carboxylic acid (cis-ACHC) with alternating chirality. The local conformations of two enantiomeric cis-ACHC dimer units suggested that a chiral ß-peptide may adopt both right-handed and left-handed helical conformations in solution. To probe the conformational behavior of 12/10-helical ß-peptides, the two reference helices with a single handedness were synthesized with a more rigidified cis-ACHC derivative. Comparison with these reference helices at low temperature revealed that a chiral cis-ACHC oligomer with alternating chirality indeed displays 12/10-helical conformations with both handedness that equilibrate rapidly in solution. This is a very rare example of chiral oligomers with helix inversion ability. The 12/10-helical backbone should be a valuable addition to potential scaffolds for applications involving helices with dynamic folding propensity.

Propensities of peptides containing the Asn-Gly segment to form ß-turn and ß-hairpin structures.

The propensities of peptides that contain the Asn-Gly segment to form ß-turn and ß-hairpin structures were explored using the density functional methods and the implicit solvation model in CH2 Cl2 and water. The populations of preferred ß-turn structures varied depending on the sequence and solvent polarity. In solution, ß-hairpin structures with ßI' turn motifs were most preferred for the heptapeptides containing the Asn-Gly segment regardless of the sequence of the strands. These preferences in solution are consistent with the corresponding X-ray structures. The sequence, H-bond strengths, solvent polarity, and conformational flexibility appeared to interact to determine the preferred ß-hairpin structure of each heptapeptide, although the ß-turn segments played a role in promoting the formation of ß-hairpin structures and the ß-hairpin propensity varied. In the heptapeptides containing the Asn-Gly segment, the ß-hairpin formation was enthalpically favored and entropically disfavored at 25°C in water. The calculated results for ß-turns and ß-hairpins containing the Asn-Gly segment imply that these structural preferences may be useful for the design of bioactive macrocyclic peptides containing ß-hairpin mimics and the design of binding epitopes for protein-protein and protein-nucleic acid recognitions. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 653-664, 2016.

Remote Stereoinductive Intramolecular Nitrile Oxide Cycloaddition: Asymmetric Total Synthesis and Structure Revision of (-)-11ß-Hydroxycurvularin.

The first total synthesis and structure revision of (-)-11ß-hydroxycurvularin (1b), a macrolide possessing a ß-hydroxyketone moiety, were accomplished. The ß-hydroxyketone moiety in this natural product was introduced by cleavage of the N-O bond in an isoxazoline ring that was formed diastereoselectively in a 1,5-remote stereocontrolled fashion by employing intramolecular nitrile oxide cycloaddition.

Hairpin formation promoted by the heterochiral dinipecotic acid segment: A DFT study.

Conformational preferences for the turn and ß-hairpin structures of Ala-based peptides [Ac-Ala(n)-(R)-Nip-(S)-Nip-Ala(n)-X (n = 0-2; X = NHMe or NMe2)] containing nipecotic acid (Nip) residues were carried out using the density functional M06-2X and the implicit solvation model SMD in CH2Cl2 and/or water. The turn structure of the (R)-Nip-(S)-Nip segment with a C10 H-bond between two terminal groups was found to be most preferred (populated at 98.9%) in CH2Cl2; this structure is consistent with IR and (1)H NMR results. The stabilities of the ß-hairpins containing the (R)-Nip-(S)-Nip segment as a turn motif relative to the extended structures increased with peptide sequence length. The relative strengths of the H-bonds between the carbonyl oxygen and the amide hydrogen appeared to be responsible for stabilizing the turn and ß-hairpin structures in CH2Cl2. In addition, the (R)-Nip-(S)-Nip segment exhibited the capability to be incorporated into one of the two ß-turn motifs of gramicidin S (GS). The structure of this GS derivative (GS-Nip2 ) was generally similar to the native peptide but was less hydrophobic and it is therefore expected to exhibit lower hemolytic activity; however, further experiments are needed to evaluate its antimicrobial activity. The structure of GS-Nip2 was somewhat more flexible than GS in solvents of higher polarity. Thus, our calculated results regarding the turn and ß-hairpin motifs of the (R)-Nip-(S)-Nip segment indicate that this structure might be useful for the design of bioactive macrocyclic peptides containing ß-hairpin mimics as well as binding epitopes in protein-protein and protein-nucleic acid recognitions.

Thiazole-based γ-building blocks as reverse-turn mimetic to design a gramicidin S analogue: conformational and biological evaluation.

This paper describes the ability of a new class of heterocyclic γ-amino acids named ATCs (4-amino(methyl)-1,3-thiazole-5-carboxylic acids) to induce turns when included in a tetrapeptide template. Both hybrid Ac-Val-(R or S)-ATC-Ile-Ala-NH2 sequences were synthesized and their conformations were studied by circular dichroism, NMR spectroscopy, MD simulations, and DFT calculations. It was demonstrated that the ATCs induced highly stable C9 pseudocycles in both compounds promoting a twist turn and a reverse turn conformation depending on their absolute configurations. As a proof of concept, a bioactive analogue of gramicidin S was successfully designed using an ATC building block as a turn inducer. The NMR solution structure of the analogue adopted an antiparallel ß-pleated sheet conformation similar to that of the natural compound. The hybrid α,γ-cyclopeptide exhibited significant reduced haemotoxicity compared to gramicidin S, while maintaining strong antibacterial activity.

Conformational preferences of helix foldamers of γ-peptides based on 2-(aminomethyl)cyclohexanecarboxylic acid.

The conformational preferences of helix foldamers having different sizes of the H-bonded pseudocycles have been studied for di- to octa-γ(2,3)-peptides based on 2-(aminomethyl)cyclohexanecarboxylic acid (γAmc6) with a cyclohexyl constraint on the C(α)-C(ß) bond using density functional methods. The helical structures of the γAmc6 oligopeptides with homochiral configurations are known to be much stable than those with heterochiral configurations in the gas phase and in solution (chloroform and water). In particular, it is found that the (P/M)-2.5(14)-helices are most preferred in the gas phase and in chloroform, whereas the (P/M)-2.3(12)-helices become most populated in water due to the larger helix dipole moments. As the peptide sequence becomes longer, the helix propensities of 14- and 12-helices are found to increase both in the gas phase and in solution. The γAmc6 peptides longer than octapeptide are expected to exist as a mixture of 12- and 14-helices with the similar populations in water. The mean backbone torsion angles and helical parameters of the 14-helix foldamers of γAmc6 oligopeptides are quite similar to those of 2-aminocyclohexylacetic acid oligopeptides and γ(2,3,4)-aminobutyric acid tetrapeptide in the solid state, despite the different substituents on the backbone.

A mechanistic study supports a two-step mechanism for peptide bond formation on the ribosome.

We report the feasible pathways of the quaternary model system for the ribosome-catalyzed PT reaction obtained by density functional calculations. Our results indicate that the step from the reactant complex to the first six-membered TS involving a proton shuttle via the 2'-OH of the P-site A76 in the stepwise pathway is the most favored rate-limiting step in solution. It is found that the C-O3' bond-breaking of A76 is not significant but the C-N bond formation with a tetrahedral intermediate occurs in the rate-limiting step and that the fast breakdown of the C-O3' bond is followed in the second transition state. These are consistent with recent kinetic experiments.

Impact of aza-substitutions on the preference of helix handedness for ß-peptide oligomers: a DFT study.

We studied the helix preference of the heterochiral pentamers of cis-2-aminocyclohexanecarboxylic acid (c-ACHC) and cis-2-aminocyclopentanecarboxylic acid (c-ACPC) with alternating backbone configurations by replacing Cß-aza- or Cα-aza-peptide residues using DFT methods in solution. The helix-handedness preferences of two pentamers were strongly affected by the replacement positions (i.e., chiralities) but not depending on the solvent polarity.

Exploring Conformational Preferences of Leu-enkephalin Using the Conformational Search and Double-Hybrid DFT Energy Calculations.

The conformational preferences of Leu-enkephalin (Leu-Enk) were explored by the conformational search and density functional theory (DFT) calculations. By a combination of low-energy conformers of each residue, the initial structures of the neutral Leu-Enk were generated and optimized using the ECEPP3 force field in the gas phase. These structures were reoptimized at the HF/3-21G(d) and M06-2X levels of theory with 6-31G(d) and 6-31+G(d) basis functions. We finally located the 139 structures with the relative energy <10 kcal mol-1 in the gas phase, from which the structures of the corresponding zwitterionic Leu-Enk were generated and reoptimized at the M06-2X/6-31+G(d) level of theory using the implicit solvation model based on density (SMD) in water. The conformational preferences of Leu-Enk were analyzed using Gibbs free energies corrected by single-point energies calculated at the double-hybrid DSD-PBEP86-D3BJ/def2-TZVP level of theory in the gas phase and in water. The neutral Leu-Enk dominantly adopted a folded structure in the gas phase stabilized by three H-bonds with a ßII'-bend-like motif at the Gly3-Phe4 sequence and a close contact between the side chains of Phe4 and Leu5. The zwitterionic Leu-Enk exhibited a folded structure in water stabilized by three H-bonds with double ß-bends such as a ßII' bend at the Gly2-Gly3 sequence and a ßI bend at the Gly3-Phe4 sequence. The calculated ensemble-averaged distance between CGly2 α and CLeu5 α of the zwitterionic Leu-Enk in water is consistent with the value estimated from the simulated annealing using the distance constraints derived from nuclear Overhauser effect spectroscopy (NOESY) spectra in water. Interestingly, the preferred conformations of the neutral and zwitterionic Leu-Enk are new folded structures not predicted by earlier computational studies. According to the refined model of the zwitterionic Leu-Enk bound to δ-opioid receptor (δOR), there were favorable interactions of the terminal charged groups of Leu-Enk with the side chains of charged residues of δOR as well as a favorable CAryl···H interaction of the Phe4 residue of Leu-Enk with Trp284 of δOR. Hence, these favorable interactions would induce the folded structure of the zwitterionic Leu-Enk with double ß-bends isolated in water into the "bioactive conformation" like an extended structure when binding to δOR.

Exploring Helical Folding in Oligomers of Cyclopentane-Based ϵ-Amino Acids: A Computational Study.

The conformational preferences of oligopeptides of an ϵ-amino acid (2-((1R,3S)-3-(aminomethyl)cyclopentyl)acetic acid, Amc5 a) with a cyclopentane substituent in the Cß -Cγ -Cδ sequence of the backbone were investigated using DFT methods in chloroform and water. The most preferred conformation of Amc5 a oligomers (dimer to hexamer) was the H16  helical structure both in chloroform and water. Four residues were found to be sufficient to induce a substantial H16 helix population in solution. The Amc5 a hexamer adopted a stable left-handed (M)-2.316 helical conformation with a rise of 4.8 Šper turn. The hexamer of Ampa (an analogue of Amc5 a with replacing cyclopentane by pyrrolidine) adopted the right-handed mixed (P)-2.918/16 helical conformation in chloroform and the (M)-2.416 helical conformation in water. Therefore, hexamers of ϵ-amino acid residues exhibited different preferences of helical structures depending on the substituents in peptide backbone and the solvent polarity as well as the chain length.

Exploring Helical Folding in Oligomers of Cyclopentane-Based ϵ-Amino Acids: A Computational Study.

Invited for this month's cover picture is the group of Young Kee Kang at Chungbuk National University (Republic of Korea). The cover picture shows the preferred conformation of the hexamer of ϵ-amino acid Amc5 a with a cyclopentane substituent in the backbone investigated using DFT methods in chloroform and water. The Amc5 a hexamer adopted a stable left-handed conformation with a rise of 4.8 Šper turn both in chloroform and water. However, the hexamer of Ampa (an analogue of Amc5 a with replacing cyclopentane by pyrrolidine) adopted different conformations in chloroform and in water. Read the full text of their Research Article at 10.1002/open.202100253.

Conformational Preferences of Cyclopentane-Based Oligo-δ-peptides in the Gas Phase and in Solution.

The conformational preferences of oligomers of δ-amino acid (δAc5 a) with a cyclopentyl constraint in the Cß -Cγ bond of the backbone were investigated by using DFT methods in the gas phase and in solution. The folded structures with C10 H-bonded pseudocycles were most preferred for dimer and tetramer of δAc5 a residues both in chloroform and water. However, for the hexameric Ac-(δAc5 a)6 -NHMe, the mixed H16/14 helical structure was found to be most preferred in chloroform (populated at 68 %), whereas the H14 helical structure was the most dominant conformation in water (populated at 60 %). The stability of the former was ascribed to the intrinsic conformational energy, whereas the solvation free energy was crucial to stabilize the latter. Pyrrolidine-substituted analogues of the hexameric Ac-(δAc5 a)6 -NHMe, with adjacent amine diads that are almost exactly one turn apart with two nitrogen atoms separated by ca. 5.5 Å, adopted helical structures. They are potential catalysts in nonpolar and polar solvents as they have similar structures to a helical 1 : 2 α:ß-heptapeptide that exhibited good catalytic performance in the crossed aldol condensation.

Assessment of density functionals with long-range and/or empirical dispersion corrections for conformational energy calculations of peptides.

Density functionals with long-range and/or empirical dispersion corrections, including LC-ωPBE, B97-D, ωB97X-D, M06-2X, B2PLYP-D, and mPW2PLYP-D functionals, are assessed for their ability to describe the conformational preferences of Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) in the gas phase and in water, which have been used as prototypes for amino acid residues of peptides. For both dipeptides, the mean absolute deviation (MAD) is estimated to be 0.22-0.40 kcal/mol in conformational energy and 2.0-3.2° in torsion angles φ and ψ using these functionals with the 6-311++G(d,p) basis set against the reference values calculated at the MP2/aug-cc-pVTZ//MP2/aug-cc-pVDZ level of theory in the gas phase. The overall performance is obtained in the order B2PLYP-D ≈ mPW2PLYP-D > ωB97X-D ≈ M06-2X > MP2 > LC-ωPBE > B3LYP with the 6-311++G(d,p) basis set. The SMD model at the M06-2X/6-31+G(d) level of theory well reproduced experimental hydration free energies of the model compounds for backbone and side chains of peptides with MADs of 0.47 and 4.3 kcal/mol for 20 neutral and 5 charged molecules, respectively. The B2PLYP-D/6-311++G(d,p)//SMD M06-2X/6-31+G(d) level of theory provides the populations of backbone and/or prolyl peptide bond for the alanine and proline dipeptides in water that are consistent with the observed values.

A histidine substitution confers metal binding affinity to a Schistosoma japonicum Glutathione S-transferase.

Glutathione S-transferases (GSTs) are multifunctional enzymes that are used as fusion tags on recombinant proteins in mammalian and Escherichia coli expression systems. We recently found that the Schistosoma japonicum GST (SjGST) displays weak Ni(2+) ion binding affinity. Glu26 and His79 were assumed to be its Ni(2+) binding sites based on the structure of the 26-kDa Clonorchis sinensis GST. To enhance SjGST Ni(2+) binding affinity, Glu26 was mutated to His. SjGST-E26H was expressed and purified at a high concentration of imidazole to a higher purity than wild type SjGST. In addition, human biotin protein ligase fused to SjGST-E26H was purified with a immobilized Ni affinity column.

Conformational preferences of cationic ß-peptide in water studied by CCSD(T), MP2, and DFT methods.

The conformational preferences of the cationic nylon-3 ßNM [(3R,4)-diaminobutanoic acid, dAba] dipeptide in water were explored as the first step to understand the mode of action of polymers of ßNM against phylogenetically diverse and intrinsically drug-resistant pathogenic fungi. The CCSD(T), MP2, M06-2X, ωB97X-D, B2PLYP-D3BJ, and DSD-PBEP86-D3BJ levels of theory with various basis sets were assessed for relative energies of the 45 local minima of the cationic Ac-dAba-NHMe located at the SMD M06-2X/6-31+G(d) level of theory in water against the benchmark CCSD(T)/CBS-limit energies in water. The best performance was obtained at the double-hybrid DSD-PBEP86-D3BJ/def2-QZVP level of theory with RMSD = 0.12 kcal/mol in water. The M06-2X/def2-QZVP level of theory predicted reasonably the conformational preference with RMSD = 0.38 kcal/mol in water and may be an alternative level of theory with marginal deviations for the calculation of conformational energies of relatively longer cationic peptides in water. In particular, the H 14-helical structures appeared to be the most feasible conformations for the cationic Ac-dAba-NHMe populated at 48-64% by relative free energies in water. The hexamer built from the H 14-structure of the cationic Ac-dAba-NHMe adopted a left-handed 314-helix, which has a slightly narrower radius and a longer rise than the regular 314-helix of ß-peptides. Hence, the 314-helices of oligomers or polymers of the cationic dAba residues are expected to be the active conformation to exhibit the ability to bridge between charged lipid head groups that might cause a local depression or invagination of the membrane of fungi.

Tailoring the Physicochemical Properties of Antimicrobial Peptides onto a Thiazole-Based γ-Peptide Foldamer.

Antimicrobial peptides (AMPs) are amphipathic molecules displaying broad-spectrum bactericidal activity, providing opportunities to develop a new generation of antibiotics. However, their use is limited either by poor metabolic stability or by high hemolytic activity. We herein addressed the potential of thiazole-based γ-peptide oligomers named ATCs as tunable scaffolds to design polycationic AMP mimetics. Knowing the side chain distribution along the backbone, we rationally designed facially amphiphilic sequences with bactericidal effect in the micromolar range. Since no hemolytic activity was detected up to 100 µM, this class of compounds has shown the potential for therapeutic development.

PMFF: Development of a Physics-Based Molecular Force Field for Protein Simulation and Ligand Docking.

The physics-based molecular force field (PMFF) was developed by integrating a set of potential energy functions in which each term in an intermolecular potential energy function is derived based on experimental values, such as the dipole moments, lattice energy, proton transfer energy, and X-ray crystal structures. The term "physics-based" is used to emphasize the idea that the experimental observables that are considered to be the most relevant to each term are used for the parameterization rather than parameterizing all observables together against the target value. PMFF uses MM3 intramolecular potential energy terms to describe intramolecular interactions and includes an implicit solvation model specifically developed for the PMFF. We evaluated the PMFF in three ways. We concluded that the PMFF provides reliable information based on the structure in a biological system and interprets the biological phenomena accurately by providing more accurate evidence of the biological phenomena.