Structure and dynamics of the antifungal molecules Syringotoxin-B and Syringopeptin-25A from molecular dynamics simulation

The phytopathogen Pseudomonas syringae pv. syringae produces toxic cyclic lipodepsipeptides (CLPs): nona-peptides and syringopeptins. All CLPs inhibit the growth of many fungal species, including human pathogens, although different fungi display different degrees of sensitivity. The best studied CLPs are Syringomycin-E (SR-E), Syringotoxin-B (ST-B) and Syringopeptin-25A (SP-25A). Their biological activity is affected by membrane composition and their structural differences. We previously (Matyus et al. in Eur Biophys J 35:459–467, 2006) reported the molecular features and structural preferences of SR-E in water and octane environments. Here we investigate in atomic detail the molecular features of the two other main CLP components, ST-B and SP-25A, in water and octane by 200 ns molecular dynamics simulations (MD), using distance restraints derived from NMR NOE data (Ballio et al. in Eur J Biochem 234:747–758, 1995). We have obtained three-dimensional models of ST-B and SP-25A CLPs in different environments. These models can now be used as a basis to investigate the interactions of ST-B and SP-25A with lipid membranes an important further step towards a better understanding of the antifungal and antibacterial activity of these peptides. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

Molecular dynamics computations and solid state nuclear magnetic resonance of the gramicidin cation channel.

This paper reports on a coupled approach to determining the structure of the gramicidin A ion channel, utilizing solid state nuclear magnetic resonance (NMR) of isotopically labeled gramicidin channels aligned parallel to the magnetic field direction, and molecular dynamics (MD). MD computations using an idealized right-handed beta-helix as a starting point produce a refined molecular structure that is in excellent agreement with atomic resolution solid state NMR data. The data provided by NMR and MD are complementary to each other. When applied in a coordinated manner they provide a powerful approach to structure determination in molecular systems not readily amenable to x-ray diffraction.     

Probing the interaction between cHAVc3 peptide and the EC1 domain of E-cadherin using NMR and molecular dynamics simulations

The goal of this work is to probe the interaction between cyclic cHAVc3 peptide and the EC1 domain of human E-cadherin protein. Cyclic cHAVc3 peptide (cyclo(1,6)Ac-CSHAVC-NH₂) binds to the EC1 domain as shown by chemical shift perturbations in the 2D ¹H,-¹⁵N-HSQC NMR spectrum. The molecular dynamics (MD) simulations of the EC1 domain showed folding of the C-terminal tail region into the main head region of the EC1 domain. For cHAVc3 peptide, replica exchange molecular dynamics (REMD) simulations generated five structural clusters of cHAVc3 peptide. Representative structures of cHAVc3 and the EC1 structure from MD simulations were used in molecular docking experiments with NMR constraints to determine the binding site of the peptide on EC1. The results suggest that cHAVc3 binds to EC1 around residues Y36, S37, I38, I53, F77, S78, H79, and I94. The dissociation constants (K_d values) of cHAVc3 peptide to EC1 were estimated using the NMR chemical shifts data and the estimated K_ds are in the range of .5 × 10^?5–7.0 × 10^?5 M.     

Irregular structure of the HIV fusion peptide in membranes demonstrated by solid-state NMR and MD simulations

To better understand peptide-induced membrane fusion at a molecular level, we set out to determine the structure of the fusogenic peptide FP23 from the HIV-1 protein gp41 when bound to a lipid bilayer. An established solid-state ¹⁹F nuclear magnetic resonance (NMR) approach was used to collect local orientational constraints from a series of CF₃-phenylglycine-labeled peptide analogues in macroscopically aligned membranes. Fusion assays showed that these ¹⁹F-labels did not significantly affect peptide function. The NMR spectra were characteristic of well-behaved samples, without any signs of heterogeneity or peptide aggregation at 1:300 in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC). We can conclude from these NMR data that FP23 has a well-defined (time-averaged) conformation and undergoes lateral diffusion in the bilayer plane, presumably as a monomer or small oligomer. Attempts to evaluate its conformation in terms of various secondary structures, however, showed that FP23 does not form any type of regular helix or β-strand. Therefore, all-atom molecular dynamics (MD) simulations were carried out using the orientational NMR constraints as pseudo-forces to drive the peptide into a stable alignment and structure. The resulting picture suggests that FP23 can adopt multiple β-turns and insert obliquely into the membrane. Such irregular conformation explains why the structure of the fusion peptide could not be reliably determined by any biophysical method so far.     

Structure–function relationships in ShKT domain peptides: ShKT-Ts1 from the sea anemone Telmatactis stephensoni

Diverse structural scaffolds have been described in peptides from sea anemones, with the ShKT domain being a common scaffold first identified in ShK toxin from Stichodactyla helianthus. ShK is a potent blocker of voltage-gated potassium channels (K_V1.x), and an analog, ShK-186 (dalazatide), has completed Phase 1 clinical trials in plaque psoriasis. The ShKT domain has been found in numerous other species, but only a tiny fraction of ShKT domains has been characterized functionally. Despite adopting the canonical ShK fold, some ShKT peptides from sea anemones inhibit K_V1.x, while others do not. Mutagenesis studies have shown that a Lys–Tyr (KY) dyad plays a key role in K_V1.x blockade, although a cationic residue followed by a hydrophobic residue may also suffice. Nevertheless, ShKT peptides displaying an ShK-like fold and containing a KY dyad do not necessarily block potassium channels, so additional criteria are needed to determine whether new ShKT peptides might show activity against potassium channels. In this study, we used a combination of NMR and molecular dynamics (MD) simulations to assess the potential activity of a new ShKT peptide. We determined the structure of ShKT-Ts1, from the sea anemone Telmatactis stephensoni, examined its tissue localization, and investigated its activity against a range of ion channels. As ShKT-Ts1 showed no activity against K_V1.x channels, we used MD simulations to investigate whether solvent exposure of the dyad residues may be informative in rationalizing and potentially predicting the ability of ShKT peptides to block K_V1.x channels. We show that either a buried dyad that does not become exposed during MD simulations, or a partially exposed dyad that becomes buried during MD simulations, correlates with weak or absent activity against K_V1.x channels. Therefore, structure determination coupled with MD simulations, may be used to predict whether new sequences belonging to the ShKT family may act as potassium channel blockers.     

'Wave-type' structure of a synthetic hexaglycosylated decapeptide: A part of the extracellular domain of human glycophorin A

The three-dimensional structure of a glycopeptide, His-Thr*-Ser*-Thr*-Ser*-Ser*-Ser*-Val-Thr-Lys, with 2-acetamido-2-deoxy--D-galactose (GalNAc) residues linked to six adjacent amino acids from Thr-10 to Ser-15, was studied by NMR spectroscopy and molecular dynamics (MD) simulations. The hexaglycosylated decapeptide is part of the extracellular domain of human glycophorin A and shows an extended structure of the peptide backbone due to O-glycosylation. Furthermore, each GalNAc residue exhibits one and only one NOE contact from the NHAc proton to the backbone amide proton of the amino acid that the sugar is directly bound to. This indicates a strong preference for the orientation of all GalNAc residues towards the N-terminus. NOE build-up curves were used to determine 42 inter-proton distances that, in connection with angles of the peptide backbone obtained from 3J-coupling constants, resulted in constraints for a MD simulation in water. The NMR data and the MD simulations show a preference for an extended backbone structure. The GalNAc residues are located alternatingly on opposite sides of the backbone and reduce the flexibility of the peptide backbone. The conformation of the molecule is relatively rigid and shows a 'wave-type' 3D structure of the peptide backbone within the glycosylation cluster. This new structural element is also supported by the unusual CD spectrum of the glycopeptide.     

Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.     

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using “realistic” potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.     

Mapping interactions of gastric inhibitory polypeptide with GIPR N‐terminus using NMR and molecular dynamics simulations

Glucose‐dependent insulinotropic polypeptide (gastric inhibitory polypeptide, or GIP), a 42‐amino acid incretin hormone, modulates insulin secretion in a glucose‐concentration‐dependent manner. Its insulinotropic action is highly dependent on glucose concentration that surmounts the hypoglycemia side effects associated with current therapy. In order to develop a GIP‐based anti‐diabetic therapy, it is essential to establish the 3D structure of the peptide and study its interaction with the GIP receptor (GIPR) in detail. This will give an insight into the GIP‐mediated insulin release process. In this article, we report the solution structure of GIP(1–42, human)NH₂ deduced by NMR and the interaction of the peptide with the N‐terminus of GIPR using molecular modelling methods. The structure of GIP(1–42, human)NH₂ in H₂O has been investigated using 2D‐NMR (DQF‐COSY, TOCSY, NOESY, ¹H‐¹³C HSQC) experiments, and its conformation was built by constrained MD simulations with the NMR data as constraints. The peptide in H₂O exhibits an α‐helical structure between residues Ser8 and Asn39 with some discontinuity at residues Gln29 to Asp35; the helix is bent at Gln29. This bent gives the peptide an ‘L’ shape that becomes more pronounced upon binding to the receptor. The interaction of GIP with the N‐terminus of GIPR was modelled by allowing GIP to interact with the N‐terminus of GIPR under a series of decreasing constraints in a molecular dynamics simulation, culminating with energy minimization without application of any constraints on the system. The canonical ensemble obtained from the simulation was subjected to a detailed energy analysis to identify the peptide–protein interaction patterns at the individual residue level. These interaction energies shed some light on the binding of GIP with the GIPR N‐terminus in a quantitative manner. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Voltage-dependent insertion of alamethicin at phospholipid/water and octane/water interfaces

Understanding the binding and insertion of peptides in lipid bilayers is a prerequisite for understanding phenomena such as antimicrobial activity and membrane-protein folding. We describe molecular dynamics simulations of the antimicrobial peptide alamethicin in lipid/water and octane/water environments, taking into account an external electric field to mimic the membrane potential. At cis-positive potentials, alamethicin does not insert into a phospholipid bilayer in 10 ns of simulation, due to the slow dynamics of the peptide and lipids. However, in octane N-terminal insertion occurs at field strengths from 0.33 V/nm and higher, in simulations of up to 100 ns duration. Insertion of alamethicin occurs in two steps, corresponding to desolvation of the Gln7 side chain, and the backbone of Aib10 and Gly11. The proline induced helix kink angle does not change significantly during insertion. Polyalanine and alamethicin form stable helices both when inserted in octane and at the water/octane interface, where they partition in the same location. In water, both polyalanine and alamethicin partially unfold in multiple simulations. We present a detailed analysis of the insertion of alamethicin into the octane slab and the influence of the external field on the peptide structure. Our findings give new insight into the mechanism of channel formation by alamethicin and the structure and dynamics of membrane-associated helices.     

Molecular dynamics simulations of peptides containing an unnatural amino acid: dimerization, folding, and protein binding

We have performed molecular dynamics (MD) simulations to study the dimerization, folding, and binding to a protein of peptides containing an unnatural amino acid. NMR studies have shown that the substitution of one residue in a tripeptide beta-strand by the unnatural amino acid Hao (5-HO2CCONH-2-MeO-C6H3-CO-NHNH2) modifies the conformational flexibility of the beta-strand and the hydrogen-bonding properties of its two edges: The number of hydrogen-bond donors and acceptors increases at one edge, whereas at the other, they are sterically hindered. In simulations in chloroform, the Hao-containing peptide 9 (i-PrCO-Phe-Hao-Val-NHBu) forms a beta-sheet-like hydrogen-bonded dimer, in good agreement with the available experimental data. Addition of methanol to the solution induces instability of this beta-sheet, as confirmed by the experiments. MD simulations also reproduce the folding of the synthetic peptide 1a (i-PrCO-Hao-Ut-Phe-Ile-Leu-NHMe) into a beta-hairpin-like structure in chloroform. Finally, the Hao-containing peptide, Ac-Ala-Hao-Ala-NHMe, is shown to form a stable complex with the Ras analogue, Rap1 A, in water at room temperature. Together with the available experimental data, these simulation studies indicate that Hao-containing peptides may serve as inhibitors of beta-sheet interactions between proteins.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Molecular statistics of cytochrome <Emphasis Type="Italic">c</Emphasis>: structural plasticity and molecular environment

Nuclear magnetic resonance experiments performed on yeast mitochondrial cytochrome c (Cytc), a paradigmatic electron transfer protein, reveal that the two oxidation states have similar structures, but different mobility: despite the few structural differences compared with the reduced form, the oxidized form displays a larger unfolding propensity. Molecular dynamics simulations performed on both NMR reduced and NMR oxidized forms show that the reduced form has a larger solvent-accessible surface area (SASA). Starting from this observation, a molecular statistical approach was then applied in order to correlate the molecular surface to molecular mobility. Simulations started from biased initial conditions corresponding to different molecular sizes were combined with the maximal constrained entropy method. The NMR structure of oxidized Cytc is more suited to expose a smaller SASA than the NMR structure of the reduced form, but the accessible conformational landscape at 300 K around the NMR oxidized structure is flatter than for the NMR reduced structure. Protein configurations of smaller SASA and size display larger plasticity when they resemble the NMR oxidized structure, whereas they are more rigid when they resemble the NMR reduced structure. Implications of the results for the protein properties during its functional process are discussed. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

3D-Structure of the interior fusion peptide of HGV/GBV-C by 1H NMR, CD and molecular dynamics studies

In this work, we present a structural characterization of the putative fusion peptide E2(279-298) corresponding to the E2 envelope protein of the HGV/GBV-C virus by (1)H NMR, CD and MD studies performed in H(2)O/TFE and in lipid model membranes. The peptide is largely unstructured in water, whereas in H(2)O/TFE and in model membranes it adopts an helical structure (approximately 65-70%). The partitioning free energy DeltaG ranges from -6 to -7.5 kcal mol(-1). OCD measurements on peptide-containing hydrated and oriented lipid multilayers showed that the peptide adopts a predominantly surface orientation. The (1)H NMR data (observed NOEs, deuterium exchange rates, Halpha chemical shift index and vicinal coupling constants) and the molecular dynamics calculations support the conclusions that the peptide adopts a stable helix in the C-terminal 9-18 residues slightly inserted into the lipid bilayer and a major mobility in the amino terminus of the sequence (1-8 residues).     

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.     

The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study

The role of peptides Abeta40 and Abeta42 in the early pathogenesis of Alzheimer's disease (AD) is frequently emphasized in the literature. It is known that Abeta42 is more prone to aggregation than Abeta40, even though they differ in only two (IA) amino acid residues at the C-terminal end. A direct comparison of the ensembles of conformations adopted by the monomers in solution has been limited by the inherent flexibility of the unfolded peptides. Here, we characterize the conformations of Abeta40 and Abeta42 in water by using a combination of molecular dynamics (MD) and measured scalar (3)J(HNHalpha) data from NMR experiments. We perform replica exchange MD (REMD) simulations and find that classical forcefields reproduce the NMR data quantitatively when the sampling is extended to the microseconds time-scale. Using the quantitative agreement of the NMR data as a validation of the model, we proceed to compare the conformational ensembles of the Abeta40 and Abeta42 peptide monomers. Our analysis confirms the existence of structured regions within the otherwise flexible Abeta peptides. We find that the C terminus of Abeta42 is more structured than that of Abeta40. The formation of a beta-hairpin in the sequence (31)IIGLMVGGVVIA involving short strands at residues 31-34 and 38-41 (in bold) reduces the C-terminal flexibility of the Abeta42 peptide and may be responsible for the higher propensity of this peptide to form amyloids.     

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.     

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.