Folding Stability Research Articles

Statistical modeling to quantify the uncertainty of FoldX-predicted protein folding and binding stability

BackgroundComputational methods of predicting protein stability changes upon missense mutations are invaluable tools in high-throughput studies involving a large number of protein variants. However, they are limited by a wide variation in accuracy and difficulty of assessing prediction uncertainty. Using a popular computational tool, FoldX, we develop a statistical framework that quantifies the uncertainty of predicted changes in protein stability.ResultsWe show that multiple linear regression models can be used to quantify the uncertainty associated with FoldX prediction for individual mutations. Comparing the performance among models with varying degrees of complexity, we find that the model precision improves significantly when we utilize molecular dynamics simulation as part of the FoldX workflow. Based on the model that incorporates information from molecular dynamics, biochemical properties, as well as FoldX energy terms, we can generally expect upper bounds on the uncertainty of folding stability predictions of ± 2.9 kcal/mol and ± 3.5 kcal/mol for binding stability predictions. The uncertainty for individual mutations varies; our model estimates it using FoldX energy terms, biochemical properties of the mutated residue, as well as the variability among snapshots from molecular dynamics simulation.ConclusionsUsing a linear regression framework, we construct models to predict the uncertainty associated with FoldX prediction of stability changes upon mutation. This technique is straightforward and can be extended to other computational methods as well.

Structural characterization of PaaX, the main repressor of the phenylacetate degradation pathway in Escherichia coli W: A novel fold of transcription regulator proteins

PaaX is a transcriptional repressor of the phenylacetic acid (PAA) catabolic pathway, a central route for bacterial aerobic degradation of aromatic compounds. Induction of the route is achieved through the release of PaaX from its promoter sequences by the first compound of the pathway, phenylacetyl-coenzyme A (PA-CoA). We report the crystal structure of PaaX from Escherichia coli W. PaaX displays a novel type of fold for transcription regulators, showing a dimeric conformation where the monomers present a three-domain structure: an N-terminal winged helix-turn-helix domain, a dimerization domain similar to the Cas2 protein and a C-terminal domain without structural homologs. The domains are separated by a crevice amenable to harbour a PA-CoA molecule. The biophysical characterization of the protein in solution confirmed several hints predicted from the structure, i.e. its dimeric conformation, a modest importance of cysteines and a high dependence of solubility and thermostability on ionic strength. At a moderately acidic pH, the protein formed a stable folding intermediate with remaining α-helical structure, a disrupted tertiary structure and exposed hydrophobic patches. Our results provide valuable information to understand the stability and mechanism of PaaX and pave the way for further analysis of other regulators with similar structural configurations.

Anti-CD19/CD8 bispecific T cell engager for the potential treatment of B cell malignancies

The administration of blinatumomab was accompanied by several adverse effects, including activation of regulatory T-cells and cytokine storm. The objective of this study was to produce and evaluate a novel αCD8/CD19 BiTE (αCD8/CD19) with the potency to directly target CD8+T-cells. In-silico studies were utilized for determining proper folding, receptor binding, and structural stability of αCD8/CD19 protein. Western blotting and indirect surface staining were used to evaluate the size accuracy and binding potency of the purified protein. Functionality was assessed for granzyme B production, cytotoxicity, and proliferation. TheαCD8/CD19recombinant protein was produced in the CHO-K1 cell line with a final concentration of 1.94 mg/l. The αCD8/CD19 bound to CD8+and CD19+cell lines and induced significant granzyme B production, cytotoxic activity and proliferation potential in the presence of IL-2 and tumor target cells. The maximum CD8+T-cell biological activity was observed on the 10th day with 10:1 effector-to-target ratio.

Protein Folding Stability and Kinetics in Alginate Hydrogels.

Proteins are commonly encapsulated in alginate gels for drug delivery and tissue-engineering applications. However, there is limited knowledge of how encapsulation impacts intrinsic protein properties such as folding stability or unfolding kinetics. Here, we use fast relaxation imaging (FReI) to image protein unfolding in situ in alginate hydrogels after applying a temperature jump. Based on changes in the Förster resonance energy transfer (FRET) response of FRET-labeled phosphoglycerate kinase (PGK), we report the quantitative impact of multiple alginate hydrogel concentrations on protein stability and folding dynamics. The gels stabilize PGK by increasing its melting temperature up to 18.4 °C, and the stabilization follows a nonmonotonic dependence on the alginate density. In situ kinetic measurements also reveal that PGK deviates more from two-state folding behavior in denser gels and that the gel decreases the unfolding rate and accelerates the folding rate of PGK, compared to buffer. Phi-value analysis suggests that the folding transition state of an encapsulated protein is structurally similar to that of folded protein. This work reveals both beneficial and negative impacts of gel encapsulation on protein folding, as well as potential mechanisms contributing to altered stability.

Retracing the evolution of a modern periplasmic binding protein.

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Small-Incisional Double Eyelid Blepharoplasty: A Retrospective Study of Our Minimally Invasive Technique with Three Mini Incisions.

Small-incisional double eyelid surgery has been increasingly performed these years and achieved good aesthetic results, but the techniques vary greatly between literatures. The authors reviewed the cases of three-small-incisional double eyelid surgery performed in the past three years and introduced their surgical technique in detail. A total of 87 patients receiving bilateral three-small-incisional double eyelid surgery were included in this retrospective study. The pretarsal folds were designed meticulously, along which three evenly distributed 2mm-long incisions were made. A minimal amount of orbicularis oculi muscle and pretarsal soft tissue were removed to expose the pretarsal fascia for further fixation. An appropriate amount of orbital septal fat was removed through the lateral incision if required. The superficial orbicularis oculi muscle and dermis on the lower margin of the incision were fixed onto the pretarsal fascia with some underlying tarsus on the upper margin of the tarsus. The skin was closed by one stitch for each incision. The average follow-up was 9.9 ± 5.2 months (range: 6-27 months). All the patients were satisfied with the result. None of them experienced loss of the pretarsal fold, bilateral asymmetry, scar hyperplasia, or persistent swelling after surgery. Our three-small-incisional technique with minor soft tissue debulking offers a simple, safe, and reproducible approach to double eyelids. It can create a stable and natural-looking pretarsal fold with a short recovery period. This journal requires that authors assign a level of evidence to each article. For a full description of these Evidence-Based Medicine ratings, please refer to the Table of Contents or the online Instructions to Authors www.springer.com/00266 .

The lyophilized chloroplasts store synthetic DARPin G3 as bioactive encapsulated organelles

BackgroundThe high cost of fermentation, purification, cold storage and transportation, short shelf life, and sterile delivery methods of biopharmaceuticals, is a matter for producers and consumers as well. Since the FDA has now approved plant cells for large-scale, cost-effective biopharmaceutical production, the isolation and lyophilization of transplastomic chloroplasts can cover concerns about limitations. DARPins are engineered small single-domain proteins that have been selected to bind to HER2 with high affinity and specificity. HER2 is an oncogene involved in abnormal cell growth in some cancers and the target molecule for cancer immunotherapy.ResultsIn this study, we reported the prolonged stability and functionality of DARPin G3 in lyophilized transplastomic tobacco leaves and chloroplasts. Western blot analysis of lyophilized leaves and chloroplasts stored at room temperature for up to nine months showed that the DARPin G3 protein was stable and preserved proper folding. Lyophilization of leaves and isolated chloroplasts increased DARPin G3 protein concentrations by 16 and 32-fold, respectively. The HER2-binding assay demonstrated that the chloroplast-made DARPin G3 can maintain its stability and binding activity without any affinity drop in lyophilized leaf materials throughout this study for more than nine months at room temperature.ConclusionLyophilization of chloroplasts expressing DARPin G3 would further reduce costs and simplify downstream processing, purification, and storage. Compressed packages of lyophilized chloroplasts were much more effective than lyophilized transplastomic leaves considering occupied space and downstream extraction and purification of DARPin G3 after nine months. These methods facilitate any relevant formulation practices for these compounds to meet any demand-oriented needs.

Important Features for Protein Foldings in Two Acyl Carrier Proteins from Enterococcus faecalis.

The emergence of multi-drug resistant Enterococcus faecalis raises a serious threat to global public health. E. faecalis is a gram-positive intestinal commensal bacterium found in humans. E. faecalis can endure extreme environments such as high temperature, pressure, and high salt, which facilitates them to cause infection in hospitals. E. faecalis has two acyl carrier proteins, AcpA (EfAcpA) in de novo fatty acid synthesis (FAS) and AcpB (EfAcpB) which utilizes exogenous fatty acids. Previously, we determined the tertiary structures of these two ACPs and investigated their structure-function relationships. Solution structures revealed that overall folding of these two ACPs is similar to those of other bacterial ACPs. However, circular dichroism (CD) experiments showed that the melting temperature of EfAcpA is 76.3°C and that of EfAcpB is 79.2°C, which are much higher than those of other bacterial ACPs. In this study, to understand the origin of their structural stabilities, we verified the important residues for stable folding of these two ACPs by monitoring thermal and chemical denaturation. Hydrogen/deuterium exchange and chemical denaturation experiments on wild-type and mutant proteins revealed that Ile10 of EfAcpA and Ile14 of EfAcpB mediate compact intramolecular packing and promote high thermostability and stable folding. E. faecalis may maximize efficiency of FAS and increase adaptability to the environmental stress by having two thermostable ACPs. This study may provide insight into bacterial adaptability and development of antibiotics against multi-drug-resistant E. faecalis.

Impact of Deep Eutectic Solvents on Kinetics and Folding Stability of Formate Dehydrogenase

Specifically designed co-solvent mixtures are an efficient way to enhance the kinetics of enzyme-catalyzed reactions without compromising enzyme stability; among them, several deep eutectic solvents have emerged as exciting co-solvent mixtures for biocatalytic reactions. DESs nature allows one to tailor the enzyme-co-solvent interactions by using DESs constituents of diverse functional groups. In this work, the influence of co-solvents (betaine, glycerol, and sorbitol) and two DESs (betaine:glycerol and betaine:sorbitol) on the kinetics of candida boidinii Formate dehydrogenase was evaluated. The results showed a 30% increase in catalytic efficiency by adding 15 wt.-% of betaine to the buffered aqueous reaction media. Further, cbFDH folded-state stability was evaluated using differential scanning fluorimetry to finally obtain the binding affinity, unfolding curves, and thermodynamic unfolding parameters. The addition of glycerol, sorbitol, and DESs increased cbFDH protection against thermal stress, and this effect could be improved by increasing co-solvent concentrations. Moreover, DESs showed the ability to reduce the irreversibility of the unfolding process. Betaine was the only co-solvent that had a negative stability effect, which was offset by using betaine-based DESs. The latter was a result of the additivity of certain individual co-solvent effects on thermal stability. Non-monotonous stability effects were obtained by adding sorbitol to the buffer solutions, probably because hydrogen bond dynamics between cbFDH/co-solvent/water change dramatically with the amount of water present. Finally, DESs improved NAD+ binding affinity with cbFDH interestingly without direct correlation with the results obtained for kinetics.

Optimization of fracture reduction robot controller based on improved sparrow algorithm

The accuracy of a fracture reduction robot (FRR) is critical for ensuring the safety of surgery. Improving the repositioning accuracy of a FRR, reducing the error, and realizing a safer and more stable folding motion is critical. To achieve this, a sparrow search algorithm (SSA) based on the Levy flight operator was proposed in this study for self-tuning the robot controller parameters. An inverse kinematic analysis of the FRR was also performed. The robot dynamics model was established using Simulink, and the inverse dynamics controller for the fracture reduction mechanism was designed using the computed torque control method. Both simulation and physical experiments were also performed. The actual motion trajectory of the actuator drive rod and its error with a desired trajectory was obtained through simulation. An optimized Levy-sparrow search algorithm (Levy-SSA) crack reduction robot controller demonstrated an overall reduction of two orders of magnitude in the reduction error, with an average error reduction of 98.74% compared with the traditional unoptimized controller. The Levy-SSA increased the convergence of the crack reduction robot control system to the optimal solution, improved the accuracy of the motion trajectory, and exhibited important implications for robot controller optimization.

Removal of N-Terminal Peptide Impacts Structural Aspects of an IgE-Reactive Recombinant Der p 5

(1) Background: Modification of the structural elements of allergens is widely used in the field of allergies. The goal of the present research was to express, purify, and characterize the shortened recombinant group 5 allergen of Dermatophagoides pteronyssinus (rDer p 5). (2) Methods: rDer p 5 storage stability and aggregation capacity were explored through in silico analysis, dynamic light scattering (DLS), and SDS-PAGE. Serum IgE reactivity and cytokine amount were investigated in sera or cell culture supernatants through ELISA, MULTIPLEX®, and Western blot analysis using sera from sensitized humans from Brazil, Colombia, and Ecuador. (3) Results: Dimeric rDer p 5 was detected through native PAGE, and this result was confirmed by data from DLS. The protein was thermically stable, as it did not degrade at 4 °C for 21 days. The shortened rDer p 5 was classified as a major IgE allergen in Brazil and Colombia, but minor in Ecuador. IL-13, IL-10, IL-1β, and IL-6 were significantly elevated in the sera of rDer p 5-reactive patients. The same cytokines plus IL-5 were more secreted by human cells upon rDer p 5 stimulation. (4) Conclusions: N-terminal peptide deletion led to a higher rDer p 5 folding stability, which, even though dimeric, was an IgE-reactive protein. Therefore, rDer p 5 could be used for molecular diagnostic applications or as backbone for hypoallergen design.

Simulations of a protein fold switch reveal crowding-induced population shifts driven by disordered regions

Macromolecular crowding effects on globular proteins, which usually adopt a single stable fold, have been widely studied. However, little is known about crowding effects on fold-switching proteins, which reversibly switch between distinct folds. Here we study the mutationally driven switch between the folds of GA and GB, the two 56-amino acid binding domains of protein G, using a structure-based dual-basin model. We show that, in the absence of crowders, the fold populations PA and PB can be controlled by the strengths of contacts in the two folds, κA and κB. A population balance, PA ≈ PB, is obtained for κB/κA = 0.92. The resulting model protein is subject to crowding at different packing fractions, ϕc. We find that crowding increases the GB population and reduces the GA population, reaching PB/PA ≈ 4 at ϕc = 0.44. We analyze the ϕc-dependence of the crowding-induced GA-to-GB switch using scaled particle theory, which provides a qualitative, but not quantitative, fit of our data, suggesting effects beyond a spherical description of the folds. We show that the terminal regions of the protein chain, which are intrinsically disordered only in GA, play a dominant role in the response of the fold switch to crowding effects.

The emerging role of ATP as a cosolute for biomolecular processes.

ATP is an important small molecule that appears at outstandingly high concentration within the cellular medium. Apart from its use as a source of energy and a metabolite, there is increasing evidence for important functions as a cosolute for biomolecular processes. Owned to its solubilizing kosmotropic triphosphate and hydrophobic adenine moieties, ATP is a versatile cosolute that can interact with biomolecules in various ways. We here use three models to categorize these interactions and apply them to review recent studies. We focus on the impact of ATP on biomolecular solubility, folding stability and phase transitions. This leads us to possible implications and therapeutic interventions in neurodegenerative diseases.

EnzyHTP Computational Directed Evolution with Adaptive Resource Allocation.

Directed evolution facilitates enzyme engineering via iterative rounds of mutagenesis. Despite the wide applications of high-throughput screening, building "smart libraries" to effectively identify beneficial variants remains a major challenge in the community. Here, we developed a new computational directed evolution protocol based on EnzyHTP, a software that we have previously reported to automate enzyme modeling. To enhance the throughput efficiency, we implemented an adaptive resource allocation strategy that dynamically allocates different types of computing resources (e.g., GPU/CPU) based on the specific need of an enzyme modeling subtask in the workflow. We implemented the strategy as a Python library and tested the library using fluoroacetate dehalogenase as a model enzyme. The results show that compared to fixed resource allocation where both CPU and GPU are on-call for use during the entire workflow, applying adaptive resource allocation can save 87% CPU hours and 14% GPU hours. Furthermore, we constructed a computational directed evolution protocol under the framework of adaptive resource allocation. The workflow was tested against two rounds of mutational screening in the directed evolution experiments of Kemp eliminase (KE07) with a total of 184 mutants. Using folding stability and electrostatic stabilization energy as computational readout, we identified all four experimentally observed target variants. Enabled by the workflow, the entire computation task (i.e., 18.4 μs MD and 18,400 QM single-point calculations) completes in 3 days of wall-clock time using ∼30 GPUs and ∼1000 CPUs.

Intracellular spatially-targeted chemical chaperones increase native state stability of mutant SOD1 barrel.

Amyotrophic lateral sclerosis (ALS) is a progressive neurological disorder with currently no cure. Central to the cellular dysfunction associated with this fatal proteinopathy is the accumulation of unfolded/misfolded superoxide dismutase 1 (SOD1) in various subcellular locations. The molecular mechanism driving the formation of SOD1 aggregates is not fully understood but numerous studies suggest that aberrant aggregation escalates with folding instability of mutant apoSOD1. Recent advances on combining organelle-targeting therapies with the anti-aggregation capacity of chemical chaperones have successfully reduce the subcellular load of misfolded/aggregated SOD1 as well as their downstream anomalous cellular processes at low concentrations (micromolar range). Nevertheless, if such local aggregate reduction directly correlates with increased folding stability remains to be explored. To fill this gap, we synthesized and tested here the effect of 9 ER-, mitochondria- and lysosome-targeted chemical chaperones on the folding stability of truncated monomeric SOD1 (SOD1bar) mutants directed to those organelles. We found that compound ER-15 specifically increased the native state stability of ER-SOD1bar-A4V, while scaffold compound FDA-approved 4-phenylbutyric acid (PBA) decreased it. Furthermore, our results suggested that ER15 mechanism of action is distinct from that of PBA, opening new therapeutic perspectives of this novel chemical chaperone on ALS treatment.

Folding stabilities of ribosome-bound nascent polypeptides probed by mass spectrometry

The folding of most proteins occurs during the course of their translation while their tRNA-bound C termini are embedded in the ribosome. How the close proximity of nascent proteins to the ribosome influences their folding thermodynamics remains poorly understood. Here, we have developed a mass spectrometry-based approach for determining the stabilities of nascent polypeptide chains using methionine oxidation as a folding probe. This approach enables quantitative measurement subglobal folding stabilities of ribosome nascent chains within complex protein mixtures and extracts. To validate the methodology, we analyzed the folding thermodynamics of three model proteins (dihydrofolate reductase, chemotaxis protein Y, and DNA polymerase IV) in soluble and ribosome-bound states. The data indicate that the ribosome can significantly alter the stability of nascent polypeptides. Ribosome-induced stability modulations were highly variable among different folding domains and were dependent on localized charge distributions within nascent polypeptides. The results implicated electrostatic interactions between the ribosome surface and nascent polypeptides as the cause of ribosome-induced stability modulations. The study establishes a robust proteomic methodology for analyzing localized stabilities within ribosome-bound nascent polypeptides and sheds light on how the ribosome influences the thermodynamics of protein folding.

Dynamics of the Second Extracellular Loop Control Transducer Coupling of Peptide-Activated GPCRs.

Many peptide-activated rhodopsin-like GPCRs share a β-hairpin folding motif in the extracellular loop 2 (ECL2), which interacts with the peptide ligand while at the same time being connected to transmembrane helix 3 (TM3) via a highly conserved disulfide bond. Currently, it remains unknown whether the coupling of the specifically shaped ECL2 to TM3 influences the activation of peptide-activated GPCRs. We investigated this possibility in a selection of peptide GPCRs with known structures. Most of the receptors with cysteine to alanine mutations folded like the respective wild-type and resided in the cell membrane, challenging pure folding stabilization by the disulfide bridge. G-protein signaling of the disulfide mutants was retained to a greater extent in secretin-like GPCRs than in rhodopsin-like GPCRs, while recruitment of arrestin was completely abolished in both groups, which may be linked to alterations in ligand residence time. We found a correlation between receptor activity of the neuropeptide Y2 receptor and alterations in ECL2 dynamics using engineered disulfide bridges or site-directed spin labeling and EPR spectroscopy. These data highlight the functional importance of the TM3-ECL2 link for the activation of specific signaling pathways in peptide-activated GPCRs, which might have implications for future drug discovery.

Molecular evolution steered structural adaptations in the DNA polymerase III α subunit of halophilic bacterium Salinibacter ruber.

A significant portion of the earth has a salty environment, and the literature on bacterial survival mechanisms in salty environments is limited. During molecular evolution, halophiles increase acidic amino acid residues on their protein surfaces which leads to a negatively charged surface potential that helps them to maintain the protein integrity and protect them from denaturation by competing with salt ions. Through protein family analysis, we have investigated the molecular-level adaptive features of DNA polymerase III's catalytic subunit (alpha) and its structure-function relationship. This study throws light on the novel understanding of halophilic bacterial replication and the molecular basis of salt adaptation. Comparisons of the amino acid contents and electronegativity of halophilic and mesophilic bacterial proteins revealed adaptations that allow halophilic bacteria to thrive in high salt concentrations. A significantly lower isoelectric point of halophilic bacterial proteins indicates the acidic nature. Also, an abundance of disordered regions in halophiles suggests the requirement of the salt ions that play a crucial role in their stable protein folding. Despite having similar topology, mesophilic and halophilic proteins, a set of very prominent molecular modifications was observed in the alpha subunit of halophiles.

Mega-scale experimental analysis of protein folding stability in biology and design.

Advances in DNA sequencing and machine learning are providing insights into protein sequences and structures on an enormous scale1. However, the energetics driving folding are invisible in these structures and remain largely unknown2. The hidden thermodynamics of folding can drive disease3,4, shape protein evolution5-7 and guide protein engineering8-10, and new approaches are needed to reveal these thermodynamics for every sequence and structure. Here we present cDNA display proteolysis, a method for measuring thermodynamic folding stability for up to 900,000 protein domains in a one-week experiment. From 1.8 million measurements in total, we curated a set of around 776,000 high-quality folding stabilities covering all single amino acid variants and selected double mutants of 331 natural and 148 de novo designed protein domains 40-72 amino acids in length. Using this extensive dataset, we quantified (1) environmental factors influencing amino acid fitness, (2) thermodynamic couplings (including unexpected interactions) between protein sites, and (3) the global divergence between evolutionary amino acid usage and protein folding stability. We also examined how our approach could identify stability determinants in designed proteins and evaluate design methods. The cDNA display proteolysis method is fast, accurate and uniquely scalable, and promises to reveal the quantitative rules for how amino acid sequences encode folding stability.

Computational Investigation of Conformational Properties of Short Azapeptides: Insights from DFT Study and NBO Analysis.

Azapeptides have gained much attention due to their ability to enhance the stability and bioavailability of peptide drugs. Their structural preferences, essential to understanding their function and potential application in the peptide drug design, remain largely unknown. In this work, we systematically investigated the conformational preferences of three azaamino acid residues in tripeptide models, Ac-azaXaa-Pro-NHMe [Xaa = Asn (4), Asp (5), Ala (6)], using the popular DFT functionals, B3LYP and B3LYP-D3. A solvation model density (SMD) was used to mimic the solvation effect on the conformational behaviors of azapeptides in water. During the calculation, we considered the impact of the amide bond in the azapeptide models on the conformational preferences of models 4-6. We analyzed the effect of the HB between the side-chain main chain and main-chain main-chain on the conformational behaviors of azapeptides 4-6. We found that the predicted lowest energy conformation for the three models differs depending on the calculation methods. In the gas phase, B3LYP functional indicates that the conformers tttANP-1 and tttADP-1 of azapeptides 4 and 5 correspond to the type I of β-turn, the lowest energy conformation with all-trans amide bonds. Considering the dispersion correction, B3LYP-D3 functional predicts the conformers tctANP-2 and tctADP-3 of azapeptide 4 and 5, which contain the cis amide bond preceding the Pro residue, as the lowest energy conformation in the gas phase. The results imply that azaAsx and Pro residues may involve cis-trans isomerization in the gas phase. In water, the predicted lowest energy conformer of azapeptides 4 and 5 differs from the gas phase results and depends on the calculational method. For azapeptide 6, regardless of calculation methods and phases, tttAAP-1 (β-I turn) is predicted as the lowest energy conformer. The results imply that the effect of the side chain that can form HBs on the conformational preferences of azapeptides 4 and 5 may not be negligible. We compared the theoretical results of azaXaa-Pro models with those of Pro-azaXaa models, showing that incorporating azaamino acid residue in peptides at different positions can significantly impact the folding patterns and stability of azapeptides.