Effect of conformation on the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH to its cyclic imide degradation product.

The objective of this study was to explain the increased propensity for the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH (1), a vitronectin-selective inhibitor, to its cyclic imide counterpart cyclo-(1,7)-Gly-Arg-Gly-Asu-Ser-Pro-Asp-Gly-OH (2). Therefore, we present the conformational analysis of peptides 1 and 2 by NMR and molecular dynamic simulations (MD). Several different NMR experiments, including COSY, COSY-Relay, HOHAHA, NOESY, ROESY, DQF-COSY and HMQC, were used to: (a) identify each proton in the peptides; (b) determine the sequential assignments; (c) determine the cis-trans isomerization of X-Pro peptide bond; and (d) measure the NH-HCalpha coupling constants. NOE- or ROE-constraints were used in the MD simulations and energy minimizations to determine the preferred conformations of cyclic peptides 1 and 2. Both cyclic peptides 1 and 2 have a stable solution conformation; MD simulations suggest that cyclic peptide 1 has a distorted type I beta-turn at Arg2-Gly3-Asp4-Ser5 and cyclic peptide 2 has a pseudo-type I beta-turn at Ser5-Pro6-Asp7-Gly1. A shift in position of the type I beta-turn at Arg2-Gly3-Asp4-Ser5 in peptide 1 to Ser5-Pro6-Asp7-Gly1 in peptide 2 occurs upon formation of the cyclic imide at the Asp4 residue. Although the secondary structure of cyclic peptide 1 is not conducive to succinimide formation, the reaction proceeds via neighbouring group catalysis by the Ser5 side chain. This mechanism is also supported by the intramolecular hydrogen bond network between the hydroxyl side chain and the backbone nitrogen of Ser5. Based on these results, the stability of Asp-containing peptides cannot be predicted by conformational analysis alone; the influence of anchimeric assistance by surrounding residues must also be considered.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Unfolding the fold of cyclic cysteine-rich peptides

We propose a method to extensively characterize the native state ensemble of cyclic cysteine-rich peptides. The method uses minimal information, namely, amino acid sequence and cyclization, as a topological feature that characterizes the native state. The method does not assume a specific disulfide bond pairing for cysteines and allows the possibility of unpaired cysteines. A detailed view of the conformational space relevant for the native state is obtained through a hierarchic multi-resolution exploration. A crucial feature of the exploration is a geometric approach that efficiently generates a large number of distinct cyclic conformations independently of one another. A spatial and energetic analysis of the generated conformations associates a free-energy landscape to the explored conformational space. Application to three long cyclic peptides of different folds shows that the conformational ensembles and cysteine arrangements associated with free energy minima are fully consistent with available experimental data. The results provide a detailed analysis of the native state features of cyclic peptides that can be further tested in experiment.     

NMR study of cyclic peptides with renin inhibitor activity

Several cyclic analogues of renin inhibitors, based on Glu-D-Phe-Lys motif have been investigated by NMR spectroscopy and molecular dynamics calculations (MD). The 15 membered macrocycle, resulting from Glu and Lys side-chain cyclization, exhibits conformational preference. The structural evidence from NMR shows the presence of hydrogen bond between Lys NH and Glu side-chain carbonyl, resulting in a 10 membered pseudo beta-turn-like structure. The structure of the cyclic moiety is similar in all the peptides, which takes at least two conformations around Calpha-Cbeta in Glu side chain. The restrained MD calculations further support such observations and show that the macrocycle is fairly rigid, with two conformations about the Glu Calpha-Cbeta bond. The linear peptide appendages, which are essential for activity in cyclic peptides, show an extended structure in the beta-region of Ramchandran plot. These calculations also demonstrate that for the most active peptide, two major conformers each exist about the Calpha-CO bond of the Lys, D-Trp and Leu residues. In this peptide, the cyclic moiety presents a negatively charged surface formed due to the carbonyl oxygens, which are thus available to form hydrogen bonds with the receptor. The linear fragment presents further binding sites with a surface which has the hydrophobic side chains of D-Trp, Leu and D-Met on one side and carbonyls on the other side.     

163 Enhanced sampling of peptides and proteins with a new biasing replica exchange method

Classical MD simulations (cMD) are limited by the sampling of relevant states of the peptides. Replica exchange (REMD) methods aim to search the conformational space of proteins more efficiently (reviewed in Ostermeir & Zacharias, 2013). We have developed a Hamiltonian REMD method that takes advantage of an intrinsic property of proteins, the specific Φ ? dihedral angle combinations along the polymer backbone. By employing a coupled two-dimensional biasing potential the energy barriers along the polymer backbone are reduced more effectively than by a previous approach based on a one-D biasing potential (Kannan & Zacharias, 2007). Thus, adjacent amino acids along the polymers backbone can easily switch between favourable regions in the Ramachandran plot. Additionally, energy barriers of rotameric states of amino acid side chains of proteins are also biased in the replica runs. The method improves the sampling of conformational substates of proteins at a modest number of replicas (nine replicas in the standard set-up with one replica running without biasing potential) compared to much larger numbers necessary in the case of standard temperature (T)-REMD simulations. A further improvement is achieved by a dynamical adjustment of the penalty potential levels in the replicas such that high exchange rates and improved mixing of conformations between different replicas are guaranteed. The biasing potential (BP)-REMD method turns out to be suitable to speed up both the folding of spaghetti-like test peptides and the refinement of loop decoy structures. Starting from extended structures, an α-helical oligo-alanine and β-hairpin chignolin and the Trp-cage protein fold more rapidly in near-native structures than in cMD simulations. The BP-REMD simulations not only accelerate the folding process of test proteins but also enlarge the variety of sampled configurations in conformational space. Since flexible parts of the protein can be penalized selectively, this method provides a precise tool to investigate regions of interest of the protein.     

Systematic conformational investigations of peptoids and peptoid-peptide chimeras

Peptoids are originally defined as N-substituted oligoglycine derivatives, and in a broader definition as N-substituted peptides (peptoid-peptide chimeras). Both types were systematically investigated by force field calculations. The Merck MMFF and YASARA2 force fields were shown to be, among others, the most suitable ones for conformational investigations of peptoids with no missing parameterizations, in contrast to AMBER or CHARMM. Ramachandran-like plots were calculated for dipeptoids and chimeras using energy calculations and grid searches by varying the dihedral angels PHI and PSI in steps of 10 degrees for s-cis- and s-trans amide bonds. Barriers as well as low energy conformations are compared to peptide Ramachandran plots, showing that peptoids have both, more barriers due to additional steric interactions as well as access to minimum conformations not accessible by peptides. Low energy conformations of dimers were used as starting conformations of higher oligomers of the peptoids for extensive molecular dynamics simulations over 10 or 20 ns with the YASARA2 force field and an explicit water solvent box to evaluate their potential to form secondary structural elements. Especially peptoids with aminoisobutyric acid-like monomer units were found to form left-handed or polyproline-like helices also known from less common natural peptides. Furthermore, new secondary structures appear feasible based on stable conformations outside the allowed areas of the Ramachandran plot for peptides, but allowed for peptoids.     

Modeling loop reorganization free energies of acetylcholinesterase: a comparison of explicit and implicit solvent models

The treatment of hydration effects in protein dynamics simulations varies in model complexity and spans the range from the computationally intensive microscopic evaluation to simple dielectric screening of charge-charge interactions. This paper compares different solvent models applied to the problem of estimating the free-energy difference between two loop conformations in acetylcholinesterase. Molecular dynamics (MD) simulations were used to sample potential energy surfaces of the two basins with solvent treated by means of explicit and implicit methods. Implicit solvent methods studied include the generalized Born (GB) model, atomic solvation potential (ASP), and the distance-dependent dieletric constant. By using the linear response approximation (LRA), the explicit solvent calculations determined a free-energy difference that is in excellent agreement with the experimental estimate, while rescoring the protein conformations with GB or the Poisson equation showed inconsistent and inferior results. While the approach of rescoring conformations from explicit water simulations with implicit solvent models is popular among many applications, it perturbs the energy landscape by changing the solvent contribution to microstates without conformational relaxation, thus leading to non-optimal solvation free energies. Calculations applying MD with a GB solvent model produced results of comparable accuracy as observed with LRA, yet the electrostatic free-energy terms were significantly different due to optimization on a potential energy surface favored by an implicit solvent reaction field. The simpler methods of ASP and the distance-dependent scaling of the dielectric constant both produced considerable distortions in the protein internal free-energy terms and are consequently unreliable.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Multiple loop conformations of peptides predicted by molecular dynamics simulations are compatible with nuclear magnetic resonance

The affinity and selectivity of protein-protein interactions can be fine-tuned by varying the size, flexibility, and amino acid composition of involved surface loops. As a model for such surface loops, we study the conformational landscape of an octapeptide, whose flexibility is chemically steered by a covalent ring closure integrating an azobenzene dye into and by a disulfide bridge additionally constraining the peptide backbone. Because the covalently integrated azobenzene dyes can be switched by light between a bent cis state and an elongated trans state, six cyclic peptide models of strongly different flexibilities are obtained. The conformational states of these peptide models are sampled by NMR and by unconstrained molecular dynamics (MD) simulations. Prototypical conformations and the free-energy landscapes in the high-dimensional space spanned by the phi/psi angles at the peptide backbone are obtained by clustering techniques from the MD trajectories. Multiple open-loop conformations are shown to be predicted by MD particularly in the very flexible cases and are shown to comply with the NMR data despite the fact that such open-loop conformations are missing in the refined NMR structures.     

Chain reversals in model peptides: studies of cystine-containing cyclic peptides. II. Effects of valyl residues and possible i-to-(i + 3) attractive ionic interactions on cyclization of [Cys1], [Cys6] hexapeptides

The synthesis of four N-acetyl N'-methylamide cystine-containing hexapeptides, CVPGVC, CGVVGC, CKPGEC, and CEPGKC, is described. These were used in disulfide-exchange reactions with the peptide CVPGGC as the formal oxidant. The relative propensities for peptide cyclization were thus deduced, and the tendency toward the formation of chain-reversal conformations was established quantitatively. An additional peptide, CVVVVC, was prepared but was never obtained as the cyclic monomer, demonstrating that the formation of chain-reversals in this peptide was of very low probability. Incorporation of pairs of valyl residues decreased the ease of cyclization, but it appeared that conformational flexibility in the cystine-containing hexapeptides may have compensated for substitutions which would have been expected to hinder the adoption of certain beta-turn conformations. The peptides containing ionic residues were cyclized more readily than expected, and this process was relatively insensitive to salt concentration. This observation is discussed with regard to the stabilization of beta-turns by i-to-(i + 3) ionic interactions in peptides and proteins. A method for blocking thiols was introduced as an improvement in the analysis of the equilibrium mixtures.     

Conformational interconversion in compstatin probed with molecular dynamics simulations

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR).     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

NMR and molecular modeling characterization of RGD containing peptides.

The tripeptide sequence arginine-glycine-aspartic acid (RGD) has been shown to be the key recognition segment in numerous cell adhesion proteins. The solution conformation and dynamics in DMSO-d6 of the cyclic pentapeptides, [formula: see text], a potent fibrinogen receptor antagonist, and [formula: see text], a weak fibrinogen receptor antagonist, have been characterized by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. 1H-1H distance constraints derived from two-dimensional NOE spectroscopy and torsional angle constraints obtained from 3JNH-H alpha coupling constants, combined with computer-assisted modeling using conformational searching algorithms and energy minimization have allowed several low energy conformations of the peptides to be determined. Low temperature studies in combination with molecular dynamics simulations suggest that each peptide does not exist in a single, well-defined conformation, but as an equilibrating mixture of conformers in fast exchange on the NMR timescale. The experimental results can be fit by considering pairs of low energy conformers. Despite this inherent flexibility, distinct conformational preferences were found which may be related to the biological activity of the peptides.     

Cyclic peptides — Small and big and their conformational aspects

Cyclic peptides form an interesting class of compounds for study by conformational analysis, by virtue of their unique conformational features and biological properties. The small cyclic peptides having 3-6 peptide units in their ring, show a variety of conformational characteristics such as occurrence ofcis peptide units, flexibility of peptide dimension and variety in hydrogen bonding. The different possible conformations of cyclic tri- and hexa-peptides are given and certain specific conformational features are discussed for cyclic tetra and pentapeptides. For higher cyclic peptides, the hydrogen bonding requirement for stability of the backbone of the ring, is seen to be kept to a minimum. These various features and their significance are examined and discussed in the light of energy minimization studies and analysis of available experimental data.     

Combined use of molecular dynamics simulations and NMR to explore peptide bond isomerization and multiple intramolecular hydrogen-bonding possibilities in a cyclic pentapeptide, cyclo(Gly-Pro-D-Phe-Gly-Val).

The conformational behavior of a model cyclic pentapeptide--cyclo(Gly-L-Pro-D-Phe-Gly-L-Val)--has been explored through the combined use of in vacuo molecular dynamics simulations and a range of nmr experiments (preceding paper). The molecular dynamics analysis suggests that, despite the conformational constraints imposed by formation of the pentapeptide cycle, this pentapeptide undergoes conformational transitions between various hydrogen-bonded conformations, characterized by low energy barriers. An inverse gamma turn with Pro in position i + 1 and a gamma turn with D-Phe in position i + 1 are two alternatives occurring frequently. Like other DLDDL cyclic pentapeptides, cyclo(Gly-Pro-D-Phe-Gly-Val) is also stabilized by an inverse gamma-turn structure with the beta-branched Val residue in position i + 1, and this hydrogen bond is retained in the different conformational families. The gamma-turn around D-Phe3 and the inverse gamma turn around Val5 are consistent with the nmr observations. 3JNH-CH alpha coupling constants of the all-trans forms were calculated from one of the molecular dynamics trajectories and are comparable to nmr experimental data, suggesting that the conformational states visited during the simulation are representative of the conformational distribution in solution. In addition to the equilibrium among various hydrogen-bonded all-trans conformers, the observation in nmr spectra of two sets of resonances for all peptide protons indicated a slow conformational interconversion of the Gly-Pro peptide bond between trans and cis isomers. The activation energy between these two conformers was determined experimentally by magnetization transfer and was calculated by high temperature constrained molecular dynamics simulation. Both methods yield a free energy of activation of ca. 20 kcal/mol. Furthermore, the free energy of activation is dependent on the direction of rotation of the Gly-Pro peptide bond.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

A cyclic peptide analogue of the loop III region of platelet-derived growth factor-BB is a synthetic antigen for the native protein.

We report the synthesis and characterization of a cyclic peptide analogue of the loop III region of platelet-derived growth factor (PDGF) B-chain sequence, cyclo(73Arg-Lys-Ile-Glu-Ile-Val-Arg-Lys-Lys81-Cys), incorporating a C-terminus cysteine residue for the conjugation to a carrier protein. The synthesis involved solid-phase chemistry, utilizing Fmoc-tBu chemistry and acid labile side-chain protecting groups, followed by 'head-to-tail' cyclization using the allyl-protected glutamic acid anchored on its side chain to the solid support with HATU/HOAt as the coupling agent. Conformational differences between the cyclic and its linear counterpart PDGF peptides were determined by circular dichroism measurements in aqueous media. High titre antisera were raised to both cyclic and linear peptide immunogens. Antisera raised to the cyclic peptide cross-reacted with PDGF-BB in both Western blot and ELISA, whereas antisera raised to the linear peptide had no reactivity with PDGF-BB. The cyclic peptide (conformational design analogue) produces an immunogen which is able to antigenically mimic the secondary structure of loop III of PDGF-BB and forms a basis from which further small molecular mimetics of PDGF may be designed for use as both immunogens and also potential agonists/antagonists of PDGF. Similarly constructed immunogens may also be useful in the design of vaccines which direct responses to loop regions in other target proteins.     

Computational studies on the backbone-dependent side-chain orientation induced by the (S,S)-CXC motif

Disulfide cyclization is a well-known procedure to impose conformational restriction to peptides undergoing backbone flexibility. Rigid conformations are induced only for small rings with a specific combination of amino acids. In this work, we present a computational search of the backbone and backbone-dependent side-chain orientation of two series of linear and cyclic peptide analogs. The -C[XY]C- scaffold (where X,Y is arginine, aspartic acid or alanine residue) in its open and (S,S) cyclic form was used for the design of the studied analogs. Thirty-six compounds, resulting from the extension with one residue at either the N- or the C-terminus were studied with classical MD. The local backbone conformation and the relative orientation of the X and Y side chains induced by either cyclization and/or the presence of the charged residues are discussed. From the present study it is concluded that cyclization has a great impact on the synplanar orientation of the X and Y side chains in the (S,S)Ac-XCYC-NH2 series of compounds while charge-charge interaction has only a weak synergic effect. On the contrary, the antiplanar orientation is favored in the case of (S,S)Ac-CXCY-NH2.