Capping β-hairpin with N-terminal d-amino acid stabilizes peptide scaffold

D-Amino acid substituted peptides as potential alternatives of homochiral L-configurations

On the primitive Earth, both L- and D-amino acids would have been present. However, only L-amino acids are essential blocks to construct proteins in modern life. To study the relative stability of D-amino acid substituted peptides, a variety of computational methods were applied. Ten prebiotic amino acids (Gly, Ala, Asp, Glu, Ile, Leu, Pro, Ser, Thr, and Val) were previously determined by multiple meteorite, spark discharge, and hydrothermal vent studies. Some previously reported early Earth polypeptide analogs were focused on in this study. Tripeptides composed of only Asp, Ser, and Val exemplified that different positions (i.e., N-terminus, C-terminus, and middle) made a difference in the minimal folding energy of peptides, while the chemical classification of amino acid (hydrophobic, acidic, or hydroxylic) did not show a significant difference. Hierarchical cluster analysis for dipeptides with all possible combinations of the proposed ten prebiotic amino acids and their D-amino acid substituted derivatives generated five clusters. Primordial simple polypeptides were modeled to test the significance of molecular fluctuations, secondary structure occupancies, and folding energy differences based on these clusters. We found peptides with α-helices, long β-sheets, and long loops are usually less sensitive to D-amino acid replacements in comparison to short β-sheets. Intriguingly, amongst 129 D-amino acid residues, mutation sensitivity profiles presented that the ratio of more to less stable residues was about 1. In conclusion, some combinations of a mixture of L- and D-amino acids can potentially act as essential building blocks of life.

Computational study on the polymerization reaction of d-aminopeptidase for the synthesis of d-peptides

Almost all natural proteins are composed exclusively of l-amino acids, and this chirality influences their properties, functions, and selectivity. Proteases can recognize proteins composed of l-amino acids but display lower selectivity for their stereoisomers, d-amino acids. Taking this as an advantage, d-amino acids can be used to develop polypeptides or biobased materials with higher biostability. Chemoenzymatic peptide synthesis is a technique that uses proteases as biocatalysts to synthesize polypeptides, and d-stereospecific proteases can be used to synthesize polypeptides incorporating d-amino acids. However, engineered proteases with modified catalytic activities are required to allow the incorporation of d-amino acids with increased efficiency. To understand the stereospecificity presented by proteases and their involvement in polymerization reactions, we studied d-aminopeptidase. This enzyme displays the ability to efficiently synthesize poly d-alanine-based peptides under mild conditions. To elucidate the mechanisms involved in the unique specificity of d-aminopeptidase, we performed quantum mechanics/molecular mechanics simulations of its polymerization reaction and determined the energy barriers presented by the chiral substrates. The enzyme faces higher activation barriers for the acylation and aminolysis reactions with the l-stereoisomer than with the d-substrate (10.7 and 17.7 kcal mol⁻¹ higher, respectively). The simulation results suggest that changes in the interaction of the substrate with Asn155 influence the stereospecificity of the polymerization reaction.

We studied the molecular mechanism of d-aminopeptidase for the synthesis of polypeptides incorporating d-amino acids.

Molecular dynamics-derived rotamer libraries for d-amino acids within homochiral and heterochiral polypeptides

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d-amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d-amino acids are usually constructed as the mirror images of l-amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host-guest peptide systems containing d-amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d-amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d-amino acid rotamer libraries derived from MD simulations of alanine-based host-guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d-amino acids as guest residues within the chiral l-Ala and d-Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues.

Assessing the Ability of Spectroscopic Methods to Determine the Difference in the Folding Propensities of Highly Similar β-Hairpins

We have evaluated the ability of nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopies to describe the difference in the folding propensities of two structurally highly similar cyclic β-hairpins, comparing the outcome to that of molecular dynamics simulations. NAMFIS-type NMR ensemble analysis and CD spectroscopy were observed to accurately describe the consequence of altering a single interaction site, whereas a single-site ¹³C NMR chemical shift melting curve-based technique was not.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.