DNA origami protection and molecular interfacing through engineered sequence-defined peptoids.

DNA nanotechnology has established approaches for designing programmable and precisely controlled nanoscale architectures through specific Watson-Crick base-pairing, molecular plasticity, and intermolecular connectivity. In particular, superior control over DNA origami structures could be beneficial for biomedical applications, including biosensing, in vivo imaging, and drug and gene delivery. However, protecting DNA origami structures in complex biological fluids while preserving their structural characteristics remains a major challenge for enabling these applications. Here, we developed a class of structurally well-defined peptoids to protect DNA origamis in ionic and bioactive conditions and systematically explored the effects of peptoid architecture and sequence dependency on DNA origami stability. The applicability of this approach for drug delivery, bioimaging, and cell targeting was also demonstrated. A series of peptoids (PE1-9) with two types of architectures, termed as "brush" and "block," were built from positively charged monomers and neutral oligo-ethyleneoxy monomers, where certain designs were found to greatly enhance the stability of DNA origami. Through experimental and molecular dynamics studies, we demonstrated the role of sequence-dependent electrostatic interactions of peptoids with the DNA backbone. We showed that octahedral DNA origamis coated with peptoid (PE2) can be used as carriers for anticancer drug and protein, where the peptoid modulated the rate of drug release and prolonged protein stability against proteolytic hydrolysis. Finally, we synthesized two alkyne-modified peptoids (PE8 and PE9), conjugated with fluorophore and antibody, to make stable DNA origamis with imaging and cell-targeting capabilities. Our results demonstrate an approach toward functional and physiologically stable DNA origami for biomedical applications.

Compact Peptoid Molecular Brushes for Nanoparticle Stabilization.

Controlling the interfaces and interactions of colloidal nanoparticles (NPs) via tethered molecular moieties is crucial for NP applications in engineered nanomaterials, optics, catalysis, and nanomedicine. Despite a broad range of molecular types explored, there is a need for a flexible approach to rationally vary the chemistry and structure of these interfacial molecules for controlling NP stability in diverse environments, while maintaining a small size of the NP molecular shell. Here, we demonstrate that low-molecular-weight, bifunctional comb-shaped, and sequence-defined peptoids can effectively stabilize gold NPs (AuNPs). The generality of this robust functionalization strategy was also demonstrated by coating of silver, platinum, and iron oxide NPs with designed peptoids. Each peptoid (PE) is designed with varied arrangements of a multivalent AuNP-binding domain and a solvation domain consisting of oligo-ethylene glycol (EG) branches. Among designs, a peptoid (PE5) with a diblock structure is demonstrated to provide a superior nanocolloidal stability in diverse aqueous solutions while forming a compact shell (∼1.5 nm) on the AuNP surface. We demonstrate by experiments and molecular dynamics simulations that PE5-coated AuNPs (PE5/AuNPs) are stable in select organic solvents owing to the strong PE5 (amine)-Au binding and solubility of the oligo-EG motifs. At the vapor-aqueous interface, we show that PE5/AuNPs remain stable and can self-assemble into ordered 2D lattices. The NP films exhibit strong near-field plasmonic coupling when transferred to solid substrates.

Electromagnetic bioeffects: a multiscale molecular simulation perspective.

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

Surface Dynamics and Ligand-Core Interactions of Quantum Sized Photoluminescent Gold Nanoclusters.

Quantum-sized metallic clusters protected by biological ligands represent a new class of luminescent materials; yet the understanding of structural information and photoluminescence origin of these ultrasmall clusters remains a challenge. Herein we systematically study the surface ligand dynamics and ligand-metal core interactions of peptide-protected gold nanoclusters (AuNCs) with combined experimental characterizations and theoretical molecular simulations. We show that the peptide sequence plays an important role in determining the surface peptide structuring, interfacial water dynamics and ligand-Au core interaction, which can be tailored by controlling peptide acetylation, constituent amino acid electron donating/withdrawing capacity, aromaticity/hydrophobicity and by adjusting environmental pH. Specifically, emission enhancement is achieved through increasing the electron density of surface ligands in proximity to the Au core, discouraging photoinduced quenching, and by reducing the amount of surface-bound water molecules. These findings provide key design principles for understanding the surface dynamics of peptide-protected nanoparticles and maximizing the photoluminescence of metallic clusters through the exploitation of biologically relevant ligand properties.

Intra- and Intersubunit Ion-Pair Interactions Determine the Ability of Apolipoprotein C-II Mutants To Form Hybrid Amyloid Fibrils.

The apolipoprotein family is structurally defined by amphipathic α-helical regions that interact with lipid surfaces. In the absence of lipid, human apolipoprotein (apo) C-II also forms well-defined amyloid fibrils with cross-ß structure. Formation of this ß-structure is accompanied by the burial of two charged residues, K30 and D69, that form an ion-pair within the amyloid fibril core. Molecular dynamics (MD) simulations indicate these buried residues form both intra- and intersubunit ion-pair interactions that stabilize the fibril. Mutations of the ion-pair (either K30D or D69K) reduce fibril stability and prevent fibril formation by K30D apoC-II under standard conditions. We investigated whether mixing K30D apoC-II with other mutants would overcome this loss of fibril forming ability. Co-incubation of equimolar mixtures of K30D apoC-II with wild-type, D69K, or double-mutant (K30D/D69K) apoC-II promoted the incorporation of K30D apoC-II into hybrid fibrils with increased stability. MD simulations showed an increase in the number of intersubunit ion-pair interactions accompanied the increased stability of the hybrid fibrils. These results demonstrate the important role of both intra- and intersubunit charge interactions in stabilizing apoC-II amyloid fibrils, a process that may be a key factor in determining the general ability of proteins to form amyloid fibrils.

Solution Conditions Affect the Ability of the K30D Mutation To Prevent Amyloid Fibril Formation by Apolipoprotein C-II: Insights from Experiments and Theoretical Simulations.

Apolipoproteins form amphipathic helical structures that bind lipid surfaces. Paradoxically, lipid-free apolipoproteins display a strong propensity to form cross-ß structure and self-associate into disease-related amyloid fibrils. Studies of apolipoprotein C-II (apoC-II) amyloid fibrils suggest that a K30-D69 ion pair accounts for the dual abilities to form helix and cross-ß structure. Consistent with this is the observation that a K30D mutation prevents fibril formation under standard fibril forming conditions. However, we found that fibril formation by K30D apoC-II proceeded readily at low pH and a higher salt or protein concentration. Structural analysis demonstrated that K30D apoC-II fibrils at pH 7 have a structure similar to that of the wild-type fibrils but are less stable. Molecular dynamics simulations of the wild-type apoC-II fibril model at pH 7 and 3 showed that the loss of charge on D69 at pH 3 leads to greater separation between residues K30 and D69 within the fibril with a corresponding reduction in ß-strand content around residue 30. In contrast, in simulations of the K30D mutant model at pH 7 and 3, residues D30 and D69 moved closer at pH 3, accompanied by an increase in ß-strand content around residue 30. The simulations also demonstrated a strong dominance of inter- over intramolecular contacts between ionic residues of apoC-II and suggested a cooperative mechanism for forming favorable interactions between the individual strands under different conditions. These observations demonstrate the important role of the buried K30-D69 ion pair in the stability and solution properties of apoC-II amyloid fibrils.

Gold Nanoparticles: Understanding and Designing the Gold-Bio Interface: Insights from Simulations (Small 18/2016).

Computational modeling is a virtual microscope providing molecular insight into the complex interactions occurring at nano-bio interfaces. The image shows the atomistic landscape of a peptide-adlayer on gold nanoparticles, providing the physicochemical information needed to tailor the design of efficient colorimetric biosensors. The successes and challenges of theoretical simulations in facilitating gold nanomaterial design for biomedicine are reviewed on page 2395 by I. Yarovsky and co-workers.

Understanding and Designing the Gold-Bio Interface: Insights from Simulations.

Gold nanoparticles (AuNPs) are an integral part of many exciting and novel biomedical applications, sparking the urgent need for a thorough understanding of the physicochemical interactions occurring between these inorganic materials, their functional layers, and the biological species they interact with. Computational approaches are instrumental in providing the necessary molecular insight into the structural and dynamic behavior of the Au-bio interface with spatial and temporal resolutions not yet achievable in the laboratory, and are able to facilitate a rational approach to AuNP design for specific applications. A perspective of the current successes and challenges associated with the multiscale computational treatment of Au-bio interfacial systems, from electronic structure calculations to force field methods, is provided to illustrate the links between different approaches and their relationship to experiment and applications.

Electromagnetic-field effects on structure and dynamics of amyloidogenic peptides.

Electromagnetic fields (EMFs) are ever-present, and so is the need to better understand their influence on human health and biological matter in general. The interaction between a molecular system and external EMF can alter the structure, and dynamical behaviour, and, hence, biological function of proteins with uncertain health consequences. This urges a detailed investigation of EMF-induced effects on basic protein biophysics. Here, we used all-atom non-equilibrium molecular dynamics simulations to understand and quantify the response mechanisms of the amyloidogenic apoC-II(60-70) peptides to non-ionising radiation by modelling their behaviour under external electromagnetic and electric fields of different strengths. Our simulations show high strength fields (>0.04 V/nm) cause structural changes in apoC-II(60-70) due to the peptide dipole alignment along the applied field direction, which disrupts the inherent ß-hairpin conformation known to be the intermediate state for fibril formation. The intermediate field-strength range (0.04-0.004 V/nm) causes a significant acceleration in peptide dynamics, which leads to the increased population of structures with fibril-inhibiting characteristics, such as the separated N- and C-termini and colocation of the aromatic residues at the same peptide face. In contrast, lower field strengths (<0.004 V/nm) promote the formation of the amyloid-prone hairpin structures relative to the ambient conditions. These findings suggest that intermediate-strength electromagnetic fields could be considered for designing alternative treatments of amyloid diseases, while the very high and low field strengths could be employed for engineering well-ordered fibrillar aggregates for non-medicinal applications.

Hydrogen/Deuterium Exchange and Molecular Dynamics Analysis of Amyloid Fibrils Formed by a D69K Charge-Pair Mutant of Human Apolipoprotein C-II.

Plasma apolipoproteins form amphipathic α helices in lipid environments but in the lipid-free state show a high propensity to form ß structure and self-associate into amyloid fibrils. The widespread occurrence of apolipoproteins in amyloid plaques suggests disease-related roles, specifically in atherosclerosis. To reconcile the dual abilities of apolipoproteins to form either α helices or cross-ß sheet structures, we examined fibrils formed by human apolipoprotein C-II (apoC-II). A structural model for apoC-II fibrils shows a cross-ß core with parallel ß strands, including a buried K30-D69 charge pair. We investigated the effect of abolishing this charge pair in mutant D69K apoC-II. Fluorescence studies indicated more rapid fibril formation and less solvent accessibility of tryptophan (W26) in D69K apoC-II fibrils than in wild-type (WT) fibrils. X-ray diffraction data of aligned D69K apoC-II fibrils yielded a typical cross-ß structure with increased ß sheet spacing compared to that of WT fibrils. Hydrogen/deuterium (H/D) exchange patterns were similar for D69K apoC-II fibrils compared to WT fibrils, albeit with an overall reduction in the level of slow H/D exchange, particularly around residues 29-32. Molecular dynamics simulations indicated reduced ß strand content for a model D69K apoC-II tetramer compared to the WT tetramer and confirmed an expansion of the cross-ß spacing that contributed to the formation of a stable charge pair between K69 and E27. The results highlight the importance of charge-pair interactions within the apoC-II fibril core, which together with numerous salt bridges in the flexible connecting loop play a major role in the ability of lipid-free apoC-II to form stable cross-ß fibrils.

Surface presentation of functional peptides in solution determines cell internalization efficiency of TAT conjugated nanoparticles.

Functionalizing nanoparticles with cell-penetrating peptides is a popular choice for cellular delivery. We investigated the effects of TAT peptide concentration and arrangement in solution on functionalized nanoparticles' efficacy for membrane permeation. We found that cell internalization correlates with the positive charge distribution achieved prior to nanoparticle encountering interactions with membrane. We identified a combination of solution based properties required to maximize the internalization efficacy of TAT-functionalized nanoparticles.

Dimensionality of carbon nanomaterials determines the binding and dynamics of amyloidogenic peptides: multiscale theoretical simulations.

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer's, Parkinson's and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth.

Graphitic nanoflakes modulate the structure and binding of human amylin.

Human amylin is an inherently disordered protein whose ability to form amyloid fibrils is linked to the onset of type II diabetes. Graphitic nanomaterials have potential in managing amyloid diseases as they can disrupt protein aggregation processes in biological settings, but optimising these materials to prevent fibrillation is challenging. Here, we employ bias-exchange molecular dynamics simulations to systematically study the structure and adsorption preferences of amylin on graphitic nanoflakes that vary in their physical dimensions and surface functionalisation. Our findings reveal that nanoflake size and surface oxidation both influence the structure and adsorption preferences of amylin. The purely hydrophobic substrate of pristine graphene (PG) nanoflakes encourages non-specific protein adsorption, leading to unrestricted lateral mobility once amylin adheres to the surface. Particularly on larger PG nanoflakes, this induces structural changes in amylin that may promote fibril formation, such as the loss of native helical content and an increase in ß-sheet character. In contrast, oxidised graphene nanoflakes form hydrogen bonds between surface oxygen sites and amylin, and as such restricting protein mobility. Reduced graphene oxide (rGO) flakes, featuring lower amounts of surface oxidation, are amphiphilic and exhibit substantial regions of bare carbon which promote protein binding and reduced conformational flexibility, leading to conservation of the native structure of amylin. In comparison, graphene oxide (GO) nanoflakes, which are predominantly hydrophilic and have a high degree of surface oxidation, facilitate considerable protein structural variability, resulting in substantial contact area between the protein and GO, and subsequent protein unfolding. Our results indicate that tailoring the size, oxygen concentration and surface patterning of graphitic nanoflakes can lead to specific and robust protein binding, ultimately influencing the likelihood of fibril formation. These atomistic insights provide key design considerations for the development of graphitic nanoflakes that can modulate protein aggregation by sequestering protein monomers in the biological environment and inhibit conformational changes linked to amyloid fibril formation.

Atomic Scale Structure of Self-Assembled Lipidated Peptide Nanomaterials.

ß-Peptides have great potential as novel biomaterials and therapeutic agents, due to their unique ability to self-assemble into low dimensional nanostructures, and their resistance to enzymatic degradation in vivo. However, the self-assembly mechanisms of ß-peptides, which possess increased flexibility due to the extra backbone methylene groups present within the constituent ß-amino acids, are not well understood due to inherent difficulties of observing their bottom-up growth pathway experimentally. A computational approach is presented for the bottom-up modelling of the self-assembled lipidated ß3-peptides, from monomers, to oligomers, to supramolecular low-dimensional nanostructures, in all-atom detail. The approach is applied to elucidate the self-assembly mechanisms of recently discovered, distinct structural morphologies of low dimensional nanomaterials, assembled from lipidated ß3-peptide monomers. The resultant structures of the nanobelts and the twisted fibrils are stable throughout subsequent unrestrained all-atom molecular dynamics simulations, and these assemblies display good agreement with the structural features obtained from X-ray fiber diffraction and atomic force microscopy data. This is the first reported, fully-atomistic model of a lipidated ß3-peptide-based nanomaterial, and the computational approach developed here, in combination with experimental fiber diffraction analysis and atomic force microscopy, will be useful in elucidating the atomic scale structure of self-assembled peptide-based and other supramolecular nanomaterials.

An Active Site Inhibitor Induces Conformational Penalties for ACE2 Recognition by the Spike Protein of SARS-CoV-2.

The novel RNA virus, severe acute respiratory syndrome coronavirus II (SARS-CoV-2), is currently the leading cause of mortality in 2020, having led to over 1.6 million deaths and infecting over 75 million people worldwide by December 2020. While vaccination has started and several clinical trials for a number of vaccines are currently underway, there is a pressing need for a cure for those already infected with the virus. Of particular interest in the design of anti-SARS-CoV-2 therapeutics is the human protein angiotensin converting enzyme II (ACE2) to which this virus adheres before entry into the host cell. The SARS-CoV-2 virion binds to cell-surface bound ACE2 via interactions of the spike protein (s-protein) on the viral surface with ACE2. In this paper, we use all-atom molecular dynamics simulations and binding enthalpy calculations to determine the effect that a bound ACE2 active site inhibitor (MLN-4760) would have on the binding affinity of SARS-CoV-2 s-protein with ACE2. Our analysis indicates that the binding enthalpy could be reduced for s-protein adherence to the active site inhibitor-bound ACE2 protein by as much as 1.48-fold as an upper limit. This weakening of binding strength was observed to be due to the destabilization of the interactions between ACE2 residues Glu-35, Glu-37, Tyr-83, Lys-353, and Arg-393 and the SARS-CoV-2 s-protein receptor binding domain (RBD). The conformational changes were shown to lead to weakening of ACE2 interactions with SARS-CoV-2 s-protein, therefore reducing s-protein binding strength. Further, we observed increased conformational lability of the N-terminal helix and a conformational shift of a significant portion of the ACE2 motifs involved in s-protein binding, which may affect the kinetics of the s-protein binding when the small molecule inhibitor is bound to the ACE2 active site. These observations suggest potential new ways for interfering with the SARS-CoV-2 adhesion by modulating ACE2 conformation through distal active site inhibitor binding.

Translocation of silica nanospheres through giant unilamellar vesicles (GUVs) induced by a high frequency electromagnetic field.

Membrane model systems capable of mimicking live cell membranes were used for the first time in studying the effects arising from electromagnetic fields (EMFs) of 18 GHz where membrane permeability was observed following exposure. A present lack of understanding of the mechanisms that drive such a rapid change in membrane permeabilization as well as any structural or dynamic changes imparted on biomolecules affected by high-frequency electromagnetic irradiation limits the use of 18 GHz EMFs in biomedical applications. A phospholipid, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) labelled with a fluorescent marker 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(lissamine rhodamine B sulfonyl) (rhodamine-DOPE) was used in constructing the giant unilamellar vesicles (GUVs). After three cycles of exposure, enhanced membrane permeability was observed by the internalisation of hydrophilic silica nanospheres of 23.5 nm and their clusters. All-atom molecular dynamics simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes exposed to high frequency electric fields of different field strengths showed that within the simulation timeframe only extremely high strength fields were able to cause an increase in the interfacial water dynamics characterized by water dipole realignments. However, a lower strength, high frequency EMF induced changes of the water hydrogen bond network, which may contribute to the mechanisms that facilitate membrane permeabilization in a longer timeframe.

Effects of mutation on the amyloidogenic propensity of apolipoprotein C-II(60-70) peptide.

Using experimental and computational methods we identified the effects of mutation on the structure and dynamics of the amyloidogenic peptide apoC-II(60-70), in monomeric and oligomeric states. Methionine (Met60) substitutions to hydrophilic Gln, hydrophobic Val, and methionine sulfoxide residues were investigated and the results compared with observations of fibril formation by the wild-type, Met60Gln, Met60Val, and oxidised Met60 (oxi-Met) apoC-II(60-70) peptides. ThT fluorescence measurements showed fibril formation by all peptides, however with different kinetics. The wild-type and Met60Val peptides formed fibrils fastest, while oxi-Met and Met60Gln peptides exhibited significantly longer lag phases. Molecular dynamics simulations showed that the mutated monomers exhibited structural features consistent with fibril-forming propensity, such as ß-hairpin conformation and a hydrophobic core. However, important differences to the wild-type were also noted, such as increased structural flexibility (oxi-Met and Met60Gln systems) and a broader distribution of the aromatic angle orientation, which could contribute to the different fibrillation kinetics observed in these peptides. Our results also showed that the critical nucleus size for fibril formation by apoC-II(60-70) may not be very large, since tetrameric oligomers in anti-parallel configuration were very stable within the 100 ns of simulations. The single-point mutations Met60Val and Met60Gln had no significant effect on the structural stability of the tetramer. The rate of fibril formation by apoC-II(60-70) peptides was generally much faster compared to longer apoC-II(56-76) peptides. Also, the effects of amino acid modifications on the kinetics of peptide fibril formation differ from the effects observed for apoC-II(56-76) and full-length apoC-II, suggesting that additional mechanisms are involved in fibril formation by mature apoC-II.

Surface enhanced Raman scattering artificial nose for high dimensionality fingerprinting.

Label-free surface-enhanced Raman spectroscopy (SERS) can interrogate systems by directly fingerprinting their components' unique physicochemical properties. In complex biological systems however, this can yield highly overlapping spectra that hinder sample identification. Here, we present an artificial-nose inspired SERS fingerprinting approach where spectral data is obtained as a function of sensor surface chemical functionality. Supported by molecular dynamics modeling, we show that mildly selective self-assembled monolayers can influence the strength and configuration in which analytes interact with plasmonic surfaces, diversifying the resulting SERS fingerprints. Since each sensor generates a modulated signature, the implicit value of increasing the dimensionality of datasets is shown using cell lysates for all possible combinations of up to 9 fingerprints. Reliable improvements in mean discriminatory accuracy towards 100% are achieved with each additional surface functionality. This arrayed label-free platform illustrates the wide-ranging potential of high-dimensionality artificial-nose based sensing systems for more reliable assessment of complex biological matrices.

Exploring the folding free energy landscape of insulin using bias exchange metadynamics.

The bias exchange metadynamics (BE-META) technique was applied to investigate the folding mechanism of insulin, one of the most studied and biologically important proteins. The BE-META simulations were performed starting from an extended conformation of chain B of insulin, using only eight replicas and seven reaction coordinates. The folded state, together with the intermediate states along the folding pathway were identified and their free energy was determined. Three main basins were found separated from one another by a large free energy barrier. The characteristic native fold of chain B was observed in one basin, while the other two most populated basins contained "molten-globule" conformations stabilized by electrostatic and hydrophobic interactions, respectively. Transitions between the three basins occur on the microsecond time scale. The implications and relevance of this finding to the folding mechanisms of insulin were investigated.

Systematic comparison of empirical forcefields for molecular dynamic simulation of insulin.

The use of atomistic simulation methodologies based on empirical forcefields has enhanced our understanding of many physical processes governing protein structure and dynamics. However, the forcefields used in classical modeling studies are often designed for a particular class of proteins and rely on continuous improvement and validation by comparison of simulations with experimental data. We present a comprehensive comparison of five popular forcefields for simulating insulin. The effect of each forcefield on the conformational evolution and structural properties of the peptide is analyzed in detail and compared with available experimental results. In this study we observed that different forcefields favor different structural trends. However, the all-atom forcefield CHARMM27 and the united-atom forcefield GROMOS 43A1 delivered the best representation of the experimentally observed dynamic behavior of chain B of insulin.