Nanostructure design using protein building blocks enhanced by conformationally constrained synthetic residues

Impacts of glycation and transglutaminase-catalyzed glycosylation with glucosamine on the conformational structure and allergenicity of bovine β-lactoglobulin

β-Lactoglobulin (β-LG) is recognized as the major milk allergen. In this study, the effects of transglutaminase (TGase) and glucosamine (GlcN)-catalyzed glycosylation and glycation on the conformational structure and allergenicity of β-LG were investigated. The formations of cross-linked peptides were demonstrated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and GlcN-conjugated modification was identified using matrix-assisted laser desorption ionization-time of flight-mass spectrometry (MALDI-TOF-MS). Structural analysis revealed that glycosylation and glycation of β-LG induced unfolding of the primary protein structure followed by a loss of the secondary structure. As revealed by circular dichroism (CD) spectroscopy, glycosylated β-LG exhibited the highest increase in the β-sheets from 32.6% to 40.4% (25 °C) and 44.2% (37 °C), and the percentage of α-helices decreased from 17.7% to 14.4% (25 °C) and 12.3% (37 °C), respectively. The tertiary and quaternary structures of β-LG also changed significantly during glycosylation and glycation, along with reduced free amino groups and variation in surface hydrophobicity. Immunoblotting and indirect enzyme-linked immuno sorbent assay (ELISA) analyses demonstrated that the lowest IgG- and IgE-binding capacities of β-LG were obtained following glycosylation at 37 °C, which were 52.7% and 56.3% lower than that of the native protein, respectively. The reduction in the antigenicity and potential allergenicity of glycosylated β-LG was more pronounced compared to TGase treated- and glycated β-LG, which correlated well with the structural changes. These results suggest that TGase-catalyzed glycosylation has more potential compared to glycation for mitigating the allergenic potential of milk products.

A Protocol for the Design of Protein and Peptide Nanostructure Self-Assemblies Exploiting Synthetic Amino Acids

In recent years there has been increasing interest in nanostructure design based on the self-assembly properties of proteins and polymers. Nanodesign requires the ability to predictably manipulate the properties of the self-assembly of autonomous building blocks, which can fold or aggregate into preferred conformational states. The design includes functional synthetic materials and biological macromolecules. Autonomous biological building blocks with available 3D structures provide an extremely rich and useful resource. Structural databases contain large libraries of protein molecules and their building blocks with a range of sizes, shapes, surfaces, and chemical properties. The introduction of engineered synthetic residues or short peptides into these building blocks can greatly expand the available chemical space and enhance the desired properties. Herein, we summarize a protocol for designing nanostructures consisting of self-assembling building blocks, based on our recent works. We focus on the principles of nanostructure design with naturally occurring proteins and synthetic amino acids, as well as hybrid materials made of amyloids and synthetic polymers.

The bactofilin cytoskeleton protein BacM of Myxococcus xanthus forms an extended β-sheet structure likely mediated by hydrophobic interactions

Bactofilins are novel cytoskeleton proteins that are widespread in Gram-negative bacteria. Myxococcus xanthus, an important predatory soil bacterium, possesses four bactofilins of which one, BacM (Mxan_7475) plays an important role in cell shape maintenance. Electron and fluorescence light microscopy, as well as studies using over-expressed, purified BacM, indicate that this protein polymerizes in vivo and in vitro into ~3 nm wide filaments that further associate into higher ordered fibers of about 10 nm. Here we use a multipronged approach combining secondary structure determination, molecular modeling, biochemistry, and genetics to identify and characterize critical molecular elements that enable BacM to polymerize. Our results indicate that the bactofilin-determining domain DUF583 folds into an extended β-sheet structure, and we hypothesize a left-handed β-helix with polymerization into 3 nm filaments primarily via patches of hydrophobic amino acid residues. These patches form the interface allowing head-to-tail polymerization during filament formation. Biochemical analyses of these processes show that folding and polymerization occur across a wide variety of conditions and even in the presence of chaotropic agents such as one molar urea. Together, these data suggest that bactofilins are comprised of a structure unique to cytoskeleton proteins, which enables robust polymerization.

Exploring the energy landscape of a molecular engineered analog of a tumor-homing peptide

Recently a new non-coded amino acid was designed as a replacement for Arg, to protect the tumor-homing pentapeptide CREKA (Cys-Arg-Glu-Lys-Ala) from proteases. This constrained Arg analog, denoted c(5)Arg, was engineered to also promote the stability of the CREKA bioactive conformation. The conformational profile of the CREKA analog obtained by replacing Arg by c(5)Arg has been extensively investigated in this work. Two molecular dynamics simulations-based strategies have been employed: a modified simulated annealing and replica exchange. Results obtained using both techniques show that the conformational features of the new analog fulfill the purpose of its design. The new CREKA analog not only preserves the main structural attributes found for the bioactive conformation of the parent peptide but also shows lower flexibility. Moreover, the conformational profile of the mutated peptide narrows towards the most stable structures previously observed for the parent CREKA peptide.

LpxA: a natural nanotube

UDP-N-acetylglucosamine 3-O-acyltransferase is a protein with a left-handed parallel beta-helix, which is a natural nanotube. They are associated with unusual high stability. To identify the reason behind the structural stability of beta-helical nanotubular structure, we have performed a total of 4 mus molecular dynamics simulations of the protein in implicit solvent at four different temperatures and monitored the unfolding pathway. The correlation in movement between different regions of the nanotubular structure has been identified from the dynamical cross-correlation map and contribution of some specific residues towards unfolding transition has been identified by principal component analysis. Difference in stability of the three loop regions has also been characterized. Construction of the unfolding conformational energy landscape identifies the probable intermediates that can appear in the unfolding pathway of the protein.

Molecular modeling of the misfolded insulin subunit and amyloid fibril

Insulin, a small hormone protein comprising 51 residues in two disulfide-linked polypeptide chains, adopts a predominantly alpha-helical conformation in its native state. It readily undergoes protein misfolding and aggregates into amyloid fibrils under a variety of conditions. Insulin is a unique model system in which to study protein fibrillization, since its three disulfide bridges are retained in the fibrillar state and thus limit the conformational space available to the polypeptide chains during misfolding and fibrillization. Taking into account this unique conformational restriction, we modeled possible monomeric subunits of the insulin amyloid fibrils using beta-solenoid folds, namely, the beta-helix and beta-roll. Both models agreed with currently available biophysical data. We performed molecular dynamics simulations, which allowed some limited insights into the relative structural stability, suggesting that the beta-roll subunit model may be more stable than the beta-helix subunit model. We also constructed beta-solenoid-based insulin fibril models and conducted fiber diffraction simulation to identify plausible fibril architectures of insulin amyloid. A comparison of simulated fiber diffraction patterns of the fibril models to the experimental insulin x-ray fiber diffraction data suggests that the model fibers composed of six twisted beta-roll protofilaments provide the most reasonable fit to available experimental diffraction patterns and previous biophysical studies.

Modeling an electronic conductor based on natural peptide sequences

In this work we used quantum mechanical calculations, molecular dynamics simulations, and QM/MM methods to examine the formation of a pi-stacking ladder and the existence of charge transfer processes in a tubular nanostructure constructed using protein building blocks. Initially, the conformational properties of beta-3-thienylalanine, a synthetic amino acid with an aromatic side group, were studied using quantum mechanical calculations. Next, this amino acid was used to substitute the natural residues at selected positions of the nanotube. Molecular dynamics simulations showed that a rational design of the targeted replacements enhances the stability of the tubular nanostructure and forms a pi-stacking ladder in the inner core of the tube. Finally, QM/MM calculations on the oxidized nanotube evidenced the delocalization of the pi-electron deficiency between the thienyl side groups of the different synthetic residues. The overall results reflected that the system under study is a potential candidate to be used as a nanowire.

Stimuli-responsive properties of peptide-based copolymers studied via directional growth of self-assembled patterns on solid substrate

We studied the self-assembly of peptide-based ABA and CBC triblock-copolymers (obtained by bacterial expression) containing random coiled hydrophilic central B blocks flanked with helical A or C blocks. The A and C blocks were of different compositions with respect to the fraction of lysine residues which provided a higher pH sensitivity of the copolymer solutions. The interchain interactions of the copolymers driven by external stimuli (pH and temperature) were explored in the process of macromolecular self-assembling in the thin films of the copolymer solutions deposited on the solid substrate. The interactions involved in the macromolecular association affected the morphology of the developed patterns. The polypeptide of the B block was not involved in the formation of the secondary structures, while the A and C blocks demonstrated helical folding responsible for the intermolecular association. The mechanism of the responsive behavior of the copolymers is based on the reversible assembling of the helices into coiled-coil structures upon the change of pH or temperature. It was found that at low pH values, when electrostatic repulsion was strong and the A/C blocks unfolded, assembling yielded fractal dendrites. Increasing the pH resulted in the recovery of the helical conformation of the A/C blocks and caused a transition from the fractal to compact structures. An elevation of temperature resulted in the disruption of the dendritic structures. The reported here approach to the evaluation of the intermolecular interactions, based on the analysis of the dendritic patterns, provides a rapid and simple method for the characterization of complex processes of self-assembling biomacromolecules.

In silico molecular engineering for a targeted replacement in a tumor-homing peptide

A new amino acid has been designed as a replacement for arginine (Arg, R) to protect the tumor-homing pentapeptide CREKA (Cys-Arg-Glu-Lys-Ala) from proteases. This amino acid, denoted (Pro)hArg, is characterized by a proline skeleton bearing a specifically oriented guanidinium side chain. This residue combines the ability of Pro to induce turn-like conformations with the Arg side-chain functionality. The conformational profile of the CREKA analogue incorporating this Arg substitute has been investigated by a combination of simulated annealing and molecular dynamics. Comparison of the results with those previously obtained for the natural CREKA shows that (Pro)hArg significantly reduces the conformational flexibility of the peptide. Although some changes are observed in the backbone...backbone and side-chain...side-chain interactions, the modified peptide exhibits a strong tendency to accommodate turn conformations centered at the (Pro)hArg residue and the overall shape of the molecule in the lowest energy conformations characterized for the natural and the modified peptides exhibit a high degree of similarity. In particular, the turn orients the backbone such that the Arg, Glu, and Lys side chains face the same side of the molecule, which is considered important for bioactivity. These results suggest that replacement of Arg by (Pro)hArg in CREKA may be useful in providing resistance against proteolytic enzymes while retaining conformational features which are essential for tumor-homing activity.

Mutational analysis and allosteric effects in the HIV-1 capsid protein carboxyl-terminal dimerization domain

The carboxyl-terminal domain (CTD, residues 146-231) of the HIV-1 capsid (CA) protein plays an important role in the CA-CA dimerization and viral assembly of the human immunodeficiency virus type 1. Disrupting the native conformation of the CA is essential for blocking viral capsid formation and viral replication. Thus, it is important to identify the exact nature of the structural changes and driving forces of the CTD dimerization that take place in mutant forms. Here, we compare the structural stability, conformational dynamics, and association force of the CTD dimers for both wild-type and mutated sequences using all-atom explicit-solvent molecular dynamics (MD). The simulations show that Q155N and E159D at the major homology region (MHR) and W184A and M185A at the helix 2 region are energetically less favorable than the wild-type, imposing profound negative effects on intermolecular CA-CA dimerization. Detailed structural analysis shows that three mutants (Q155N, E159D, and W184A) display much more flexible local structures and weaker CA-CA association than the wildtype, primarily due to the loss of interactions (hydrogen bonds, side chain hydrophobic contacts, and pi-stacking) with their neighboring residues. Most interestingly, the MHR that is far from the interacting dimeric interface is more sensitive to the mutations than the helix 2 region that is located at the CA-CA dimeric interface, indicating that structural changes in the distinct motif of the CA could similarly allosterically prevent the CA capsid formation. In addition, the structural and free energy comparison of the five residues shorter CA (151-231) dimer with the CA (146-231) dimer further indicates that hydrophobic interactions, side chain packing, and hydrogen bonds are the major, dominant driving forces in stabilizing the CA interface.

Analysis of the sequence and structural features of the left-handed beta-helical fold

The left-handed parallel beta-helix (LbetaH) is a structurally repetitive, highly regular, and symmetrical fold formed by coiling of elongated beta-sheets into helical "rungs." This canonical fold has recently received interest as a possible solution to the fibril structure of amyloid and as a building block of self-assembled nanotubular structures. In light of this interest, we aimed to understand the structural requirements of the LbetaH fold. We first sought to determine the sequence characteristics of the repeats by analyzing known structures to identify positional preferences of specific residues types. We then used molecular dynamics simulations to demonstrate the stabilizing effect of successive rungs and the hydrophobic core of the LbetaH. We show that a two-rung structure is the minimally stable LbetaH structure. In addition, we defined the structure-based sequence preference of the LbetaH and undertook a genome-wide sequence search to determine the prevalence of this unique protein fold. This profile-based LbetaH search algorithm predicted a large fraction of LbetaH proteins from microbial origins. However, the relative number of predicted LbetaH proteins per specie was approximately equal across the genomes from prokaryotes to eukaryotes.

What can we learn from highly connected beta-rich structures for structural interface design?

Most hubs' binding sites are able to transiently interact with numerous proteins. We focus on beta-rich hubs with the goal of inferring features toward design. Since they are able to interact with many partners and association of beta-conformations may lead to amyloid fibrils, we ask whether there is some property that distinguishes them from low-connectivity beta-rich proteins, which may be more interaction specific. Identification of such features should be useful as they can be incorporated in interface design while avoiding polymerization into fibrils. We classify the proteins in the yeast interaction map according to the types of their secondary structures. The small number of the obtained beta-rich protein structures in the Protein Data Bank likely reflects their low occurrence in the proteome. Analysis of the obtained structures indicates that highly connected beta-rich proteins tend to have clusters of conserved residues in their cores, unlike beta-rich structures with low connectivity, suggesting that the highly packed conserved cores are important to the stability of proteins, which have residue composition and sequence prone to beta-structure and amyloid formation. The enhanced stability may hinder partial unfolding, which, depending on the conditions, is more likely to lead to polymerization of these sequences.

Conformational preferences of beta- and gamma-aminated proline analogues

Quantum mechanical calculations have been used to investigate how the incorporation of an amino group to the Cbeta- or Cgamma-positions of the pyrrolidine ring affects the intrinsic conformational properties of the proline. Specifically, a conformational study of the N-acetyl-N'-methylamide derivatives of four isomers of aminoproline, which differ not only in the beta- or gamma-position of the substituent but also in its cis or trans relative disposition, has been performed. To further understand the role of the intramolecular hydrogen bonds between the backbone carbonyl groups and the amino side group, a conformational study was also performed on the corresponding four analogues of (dimethylamino)proline. In addition, the effects of solvation on aminoproline and (dimethylamino)proline dipeptides have been evaluated using a self-consistent reaction field model, and considering four different solvents (carbon tetrachloride, chloroform, methanol and water). Results indicate that the incorporation of the amino substituent into the pyrrolidine ring affects the conformational properties, with backbone...side chain intramolecular hydrogen bonds detected when it is incorporated in a cis relative disposition. In general, the incorporation of the amino side group tends to stabilize those structures where the peptide bond involving the pyrrolidine nitrogen is arranged in cis. The aminoproline isomer with the substituent attached to the Cgamma-position with a cis relative disposition is the most stable in the gas phase and in chloroform, methanol and water solutions. Replacement of the amino side group by the dimethylamino substituent produces significant changes in the potential energy surfaces of the four investigated (dimethylamino)proline-containing dipeptides. Thus, these changes affect not only the number of minima, which increases considerably, but also the backbone and pseudorotational preferences. In spite of these effects, comparison of the conformational preferences, i.e., the more favored conformers, calculated for different isomers of aminoproline and (dimethylamino)proline dipeptides showed a high degree of consistency for the two families of compounds.

Structural analysis of a beta-helical protein motif stabilized by targeted replacements with conformationally constrained amino acids

Here we study conformational stabilization induced in a beta-helical nanostructure by position-specific mutations. The nanostructure is constructed through the self-assembly of the beta-helical building block excised from E. coli galactoside acetyltransferase (PDB code 1krr , chain A; residues 131-165). The mutations involve substitutions by cyclic, conformationally constrained amino acids. Specifically, a complete structural analysis of the Pro-Xaa-Val sequence [with Xaa being Gly, Ac 3c (1-aminocyclopropane-1-carboxylic acid) and Ac 5c (1-aminocyclopentane-1-carboxylic acid)], corresponding to the 148-150 loop region in the wild-type (Gly) and mutated (Ac 3c and Ac 5c) 1krr , has been performed using Molecular Dynamics simulations and X-ray crystallography. Simulations have been performed for the wild-type and mutants of three different systems, namely the building block, the nanoconstruct and the isolated Pro-Xaa-Val tripeptide. Furthermore, the crystalline structures of five peptides of Pro-Xaa-Val or Xaa-Val sequences have been solved by X-ray diffraction analysis and compared with theoretical predictions. Both the theoretical and crystallographic studies indicate that the Pro-Ac n c-Val sequences exhibit a high propensity to adopt turn-like conformations, and this propensity is little affected by the chemical environment. Overall, the results indicate that replacement of Gly149 by Ac 3c or Ac 5c significantly reduce the conformational flexibility of the target site enhancing the structural specificity of the building block and the nanoconstruct derived from the 1krr beta-helical motif.

The energy landscape of a selective tumor-homing pentapeptide

Recently, a potentially powerful strategy based on phage-display libraries has been presented to target tumors via homing peptides attached to nanoparticles. The Cys-Arg-Glu-Lys-Ala (CREKA) peptide sequence has been identified as a tumor-homing peptide that binds to clotted plasmas proteins present in tumor vessels and interstitium. The aim of this work consists of mapping the conformational profile of CREKA to identify the bioactive conformation. For this purpose, a conformational search procedure based on modified simulated annealing combined with molecular dynamics was applied to three systems that mimic the experimentally used conditions: (i) the free peptide; (ii) the peptide attached to a nanoparticle; and (iii) the peptide inserted in a phage display protein. In addition, the free peptide was simulated in an ionized aqueous solution environment, which mimics the ionic strength of the physiological medium. Accessible minima of all simulated systems reveal a multiple interaction pattern involving the ionized side chains of Arg, Glu, and Lys, which induces a beta-turn motif in the backbone observed in all simulated CREKA systems.

Application of 1-aminocyclohexane carboxylic acid to protein nanostructure computer design

Conformationally restricted amino acids are promising candidates to serve as basic pieces in redesigned protein motifs which constitute the basic modules in synthetic nanoconstructs. Here we study the ability of constrained cyclic amino acid 1-aminocyclohexane-1-carboxylic acid (Ac6c) to stabilize highly regular beta-helical motifs excised from naturally occurring proteins. Calculations indicate that the conformational flexibility observed in both the ring and the main chain is significantly higher than that detected for other 1-aminocycloalkane-1-carboxylic acids (Acnc, where n refers to the size of the ring) with smaller cycles. Incorporation of Ac6c into the flexible loops of beta-helical motifs indicates that the stability of such excised building blocks as well as the nanoassemblies derived from them is significantly enhanced. Thus, the intrinsic Ac6c tendency to adopt folded conformations combined with the low structural strain of the cyclohexane ring confers the ability to both self-adapt to the beta-helix motif and to stabilize the overall structure by absorbing part of its conformational fluctuations. Comparison with other Acnc residues indicates that the ability to adapt to the targeted position improves considerably with the ring size, i.e., when the rigidity introduced by the strain of the ring decreases.

Testing β-helix terminal coils stability by targeted substitutions with non-proteogenic amino acids: A molecular dynamics study

The search for new building block templates useful for nanostructures design, targets protein motifs with a wide range of structures. Stabilizing these building blocks when extracted from their natural environment becomes a fundamental goal in order to successfully control their assembly. Targeted replacements of natural residues by conformationally constrained amino acids were shown to be a successful strategy to achieve such stabilization. In this work, the effect of replacing natural amino acids by non-proteogenic residues in a beta-helix building block has been evaluated using extensive molecular dynamics simulations. Here, we focus on systematic substitutions of valine residues present in beta-sheet segments of a beta-helical building block excised from Escherichia coli galactoside acetyltransferase, residues 131-165. Four different types of non-proteogenic amino acids have been considered for substitution: (i) one dehydroamino acid, (ii) two d-amino acids, (iii) one beta-amino acid and (iv) two alpha,alpha-dialkylamino acids. Our results indicate that the ability of non-proteogenic amino acids to stabilize small building block motifs is site-dependent. We conclude that if the replacement does not alter the energy balance between attractive non-covalent interactions and steric hindrance, synthetic residues are suitable candidates to nucleate beta-helix formation.

RNAJunction: a database of RNA junctions and kissing loops for three-dimensional structural analysis and nanodesign

We developed a database called RNAJunction that contains structure and sequence information for RNA structural elements such as helical junctions, internal loops, bulges and loop-loop interactions. Our database provides a user-friendly way of searching structural elements by PDB code, structural classification, sequence, keyword or inter-helix angles. In addition, the structural data was subjected to energy minimization. This database is useful for analyzing RNA structures as well as for designing novel RNA structures on a nanoscale. The database can be accessed at: http://rnajunction.abcc.ncifcrf.gov/

Metal ion-dependent, reversible, protein filament formation by designed beta-roll polypeptides

BACKGROUND

A right-handed, calcium-dependent beta-roll structure found in secreted proteases and repeat-in-toxin proteins was used as a template for the design of minimal, soluble, monomeric polypeptides that would fold in the presence of Ca2+. Two polypeptides were synthesised to contain two and four metal-binding sites, respectively, and exploit stacked tryptophan pairs to stabilise the fold and report on the conformational state of the polypeptide.

RESULTS

Initial analysis of the two polypeptides in the presence of calcium suggested the polypeptides were disordered. The addition of lanthanum to these peptides caused aggregation. Upon further study by right angle light scattering and electron microscopy, the aggregates were identified as ordered protein filaments that required lanthanum to polymerize. These filaments could be disassembled by the addition of a chelating agent. A simple head-to-tail model is proposed for filament formation that explains the metal ion-dependency. The model is supported by the capping of one of the polypeptides with biotin, which disrupts filament formation and provides the ability to control the average length of the filaments.

CONCLUSION

Metal ion-dependent, reversible protein filament formation is demonstrated for two designed polypeptides. The polypeptides form filaments that are approximately 3 nm in diameter and several hundred nm in length. They are not amyloid-like in nature as demonstrated by their behaviour in the presence of congo red and thioflavin T. A capping strategy allows for the control of filament length and for potential applications including the "decoration" of a protein filament with various functional moieties.

Stability of tubular structures based on beta-helical proteins: self-assembled versus polymerized nanoconstructs and wild-type versus mutated sequences

In this work we used atomistic molecular dynamics simulations to examine different aspects of tubular nanostructures constructed using protein building blocks with a beta-helical conformation. Initially, we considered two different natural protein building blocks, which were extracted from the protein data base, to compare the relative stabilities of the nanotubes obtained made of self-assembled and covalently linked repeats. Results show nanotubes constructed by linking building blocks through covalent bonds are very stable suggesting that the basic principles of polymer physics are valid when the repeating units are made of large fragments of proteins. In contrast, the stability of self-assembled nanostructures strongly depends on the attractive nonbonding interactions associated to building blocks aligned in a complementary manner. On the other hand, we investigated the ability of a conformationally constrained synthetic amino acid to enhance the stability of both self-assembled and polymerized nanotubes when it is used to substitute natural residues. Specifically, we considered 1-aminocyclopentane-1-caboxylic acid, which involves strong stereochemical constraints produced by the cyclopentane side chain. We found that the incorporation of this amino acid within the more flexible regions of the beta-helical building blocks is an excellent strategy to enhance the stability of the nanotubes. Thus, when a single mutation is performed in the loop region of the beta-helix, the bend architecture of the whole loop is stabilized since the conformational mobility is reduced not only at the mutated position but also at the adjacent positions.

Coarse-grained representation of beta-helical protein building blocks

A general strategy to develop coarse-grained models of beta-helical protein fragments is presented. The procedure has been applied to a building block formed by a two-turn repeat motif from E. coli galactoside acetyltransferase, which is able to provide a very stable self-assembled tubular nanoconstruct upon stacking of its replicas. For this purpose, first, we have developed a computational scheme to sample very efficiently the configurational space of the building block. This method, which is inspired by a strategy recently designed to study amorphous polymers and by an advanced Monte Carlo algorithm, provides a large ensemble of uncorrelated configurations at a very reasonable computational cost. The atomistic configurations provided by this method have been used to obtain a coarse-grained model that describes the amino acids with fewer particles than those required for full atomistic detail, i.e., two, three, or four depending on the chemical nature of the amino acid. Coarse-grained potentials have been developed considering the following types of interactions: (i) electrostatic and van der Waals interactions between residues i and i + n with n >/= 2; (ii) interactions between residues i and i + 1; and (c) intra-residue interactions. The reliability of the proposed model has been tested by comparing the atomistic and coarse-grained energies calculated for a large number of independent configurations of the beta-helical building block.

Conformations of proline analogues having double bonds in the ring

The intrinsic conformational preferences of proline analogues having double bonds between carbon atoms in their rings have been investigated using quantum mechanical calculations at the B3LYP/6-31+G(d,p) level. For this purpose, the potential energy surface of the N-acety-N'-methylamide derivatives of three dehydroprolines (proline analogues unsaturated at alpha,beta; beta,gamma; and gamma,delta) and pyrrole (proline analogue with unsaturations at both alpha,beta and gamma,delta) have been explored, and the results are compared with those obtained for the derivative of the nonmodified proline. We found that the double bonds affect the ring puckering and the geometric internal parameters, even though the backbone conformation was influenced the most. Results indicate that the formation of double bonds between carbon atoms in the pyrrolidine ring should be considered as an effective procedure to restrict the conformational flexibility of prolines. Interestingly, we also found that the N-acetyl-N'-methylamide derivative of pyrrole shows a high probability of having a cis peptide bond preceding the proline analogue.