Molecular dynamics simulation of adrenocorticotropin (1-10) peptide in a solvated dodecylphosphocholine micelle

Adrenocorticotropin (ACTH) (1-10), an adrenocorticotropin hormone fragment, has been studied by molecular dynamics (MD) simulation in an NPT ensemble in an explicit dodecylphosphocholine (DPC) micelle. Two starting configurations of the peptide/micelle system, corresponding to the insertion and surface-binding modes, were used. A common equilibrated configuration, in which the peptide lies parallel to the micellar surface, was reached from both simulations. In the initial part of the simulations, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated before the peptide reached an equilibrium configuration with respect to the micelle. Analyses of the trajectories from the subsequent free (unrestrained) MD simulation showed that ACTH (1-10) does not conform strictly to a helical structure. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide protons of the peptide and the micellar head groups. However, the extent of the latter interaction is less pronounced than in the negatively charged SDS micelle. The final structure enhances the amphipathic nature of the peptide, facilitating better interactions at the water-hydrophobic interface. The primary hydrophobic interactions with the micelle came from the side chains of Met4, Phe7, and Trp9. All peptide bonds were either hydrated or were involved in intramolecular hydrogen bonding. The interactions with the DPC micelle, the conformation of the bound peptide, and the dynamics of the peptide, as revealed by the time correlation functions of the N-H bonds, were compared with those of the ACTH (1-10)/SDS system studied previously by MD simulations.     

Molecular dynamics simulations of adrenocorticotropin (1-24) peptide in a solvated dodecylphosphocholine (DPC) micelle and in a dimyistoylphosphatidylcholine (DMPC) bilayer

The structure and interactions of the 1-24 fragment of the adrenocorticotropin hormone, ACTH (1-24), with membrane have been studied by molecular dynamics (MD) simulation in an NPT ensembles in two explicit membrane mimics, a dodecylphosphocholine (DPC) micelle and a dimyristoylphosphatidylcholine (DMPC) bilayer. The starting configuration of the peptide/lipid systems had the 1-10 segment of the peptide lying on the surface of the model membrane, the same as the equilibrated structure (by MD) of ACTH (1-10) in a DPC micelle. The simulations showed that the peptide adopts the surface-binding mode and essentially the same structure in both systems. Thus the results of this work lend support to the assumption that micelles are reasonable mimics for biological membranes for the study of peptide binding. The 1-10 segment is slightly tilted from the parallel orientation to the interface and interacts strongly with the membrane surface while the more polar 11-24 segment shows little tendency to interact with the membrane surface, preferring to reside primarily in the aqueous phase. Furthermore, the 1-10 segment of the peptide binds to the DPC micelle in essentially the same way as ACTH (1-10). Thus the MD results are in excellent agreement with the model of interaction of ACTH (1-24) with membrane derived from NMR experiments. The secondary structure and the hydration of the peptide and the interactions of specific residues with the lipid head groups have also been analyzed.     

Molecular dynamics study of substance P peptides partitioned in a sodium dodecylsulfate micelle

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) during the equilibration stage of the simulation. It was shown that when SP-Y8 was initially placed in an insertion (perpendicular) configuration, the peptide equilibrated to a surface-bound (parallel) configuration in approximately 450 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained from the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbones of all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions come from the hydrophobic interaction between the side chains of Lys-3, Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, from Lys-3, Phe-7, Leu-10, and Met-11 in SP-Y8, and the micellar interior. Significant interactions, electrostatic and hydrogen bonding, between the N-terminal residues, Arg-Pro-Lys, and the micellar headgroups were observed. These latter interactions served to affect both the structure and, especially, the flexibility, of the N-terminus. The results from simulation of the same peptides in a water/CCl4 biphasic cell were compared with the results of the present study, and the validity of using the biphasic system as an approximation for peptide-micelle or peptide-bilayer systems is discussed.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Molecular dynamics study of substance P peptides in a biphasic membrane mimic

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.     

Simulation of the N-terminus of HIV-1 glycoprotein 41000 fusion peptide in micelles.

In this paper, the N-terminus of glycoprotein-41, the HIV-1 fusion peptide, was studied by molecular dynamics simulations in an explicit sodium dodecyl sulfate micelle. The simulation provides a detailed picture of the equilibrium structure and peptide stability as it interacts with the micelle. The equilibrium location of the peptide shows the peptide at the surface of the micelle with hydrophobic residues interacting with the micelle's core. At equilibrium, the peptide adopts an alpha-helical structure from residues 5-16 and a type-1 beta-turn from 17-20 with the other residues exhibiting more flexible conformations. The primary hydrophobic interactions with the micelle are from the leucine and phenylalanine residues (Leu-7, Phe-8, Leu-9, Phe-11, Leu-12) while the alanine and glycine residues (Ala-1, Gly-3, Gly-5, Ala-6, Gly-10, Gly-13, Ala-14, Ala-15, Gly-16, Gly-10, Ala-21) interact favorably with water molecules. The results suggest that Phe-8, part of the highly conserved FLG motif of the fusion peptide, plays a key role in the interaction of the peptide with membranes. Our simulations corroborate experimental investigations of the fusion peptide in SDS micelles, providing a high-resolution picture that explains the experimental findings.     

Helical conformations of three crystalline pentapeptide fragments of suzukacillin,a membrane channel forming polypeptide

The crystal structures of three pentapeptide fragments of suzukacillin-A have been determined. Boc-Aib-Pro-Val-Aib-Val-OMe (peptide 1–5) adopts a distorted helical conformation, stabilized by three intramolecular hydrogen bonds (two 5→1, one 4→1). Boc-Ala-Aib-Ala-Aib-Aib-OMe (peptide 6–10) and Boc-Leu-Aib-Pro-Val-Aib-OMe (peptide 16–20) adopt 3₁₀ helical structures stabilized by three and two 4→1 intramolecular hydrogen bonds, respectively. These structures provide substantial support for a largely helical conformation for the suzukacillin membrane channel.     

Evidence for a folded structure of Met-enkephalin in membrane mimetic systems: 1H-nmr studies in sodiumdodecylsulfate, lyso-phosphatidylcholine, and mixed lyso-phosphatidylcholine/sulfatide micelles

¹H-nmr spectra of Met-enkephalin dissolved in aqueous solution of sodiumdodecylsulfate (SDS) micelles are reported as a function of pH and temperature. The temperature behavior of the amide protons is compared with that observed for the same peptide dissolved in aqueous solution of lyso-phosphatidylcholine (LPC) and lyso-phosphatidylcholine-sulfatide (LPC-SH) micelles. The temperature coefficients are affected by the micelle polarity, which suggests that the peptide backbone is not remote from the micelle surface. pH titration performed in the presence of SDS micelles gives a number of intrinsic and extrinsic pK_a values, indicative of a folded structure of the opioid molecule. This conformation is characterized by the existence of an intramolecular hydrogen bond involving the Met-5 amide proton and an interaction of the N-terminal residue with the aliphatic side chains of both Phe-4 and Met-5.     

Molecular Dynamics Simulations of Adrenocorticotropin (1–24) Peptide in a Solvated Dodecylphosphocholine (DPC) Micelle and in a Dimyistoylphosphatidylcholine (DMPC) Bilayer

Abstract

The structure and interactions of the 1–24 fragment of the adrenocorticotropin hormone, ACTH (1–24), with membrane have been studied by molecular dynamics (MD) simulation in an NPT ensembles in two explicit membrane mimics, a dodecylphosphocholine (DPC) micelle and a dimyristoylphosphatidylcholine (DMPC) bilayer. The starting configuration of the peptide/lipid systems had the 1–10 segment of the peptide lying on the surface of the model membrane, the same as the equilibrated structure (by MD) of ACTH (1–10) in a DPC micelle. The simulations showed that the peptide adopts the surface-binding mode and essentially the same structure in both systems. Thus the results of this work lend support to the assumption that micelles are reasonable mimics for biological membranes for the study of peptide binding. The 1–10 segment is slightly tilted from the parallel orientation to the interface and interacts strongly with the membrane surface while the more polar 11–24 segment shows little tendency to interact with the membrane surface, preferring to reside primarily in the aqueous phase. Furthermore, the 1–10 segment of the peptide binds to the DPC micelle in essentially the same way as ACTH (1–10). Thus the MD results are in excellent agreement with the model of interaction of ACTH (1–24) with membrane derived from NMR experiments. The secondary structure and the hydration of the peptide and the interactions of specific residues with the lipid head groups have also been analyzed.     

Molecular dynamics investigation of the influence of anionic and zwitterionic interfaces on antimicrobial peptides' structure: implications for peptide toxicity and activity

Molecular dynamics simulations of three related helical antimicrobial peptides have been carried out in zwitterionic diphosphocholine (DPC) micelles and anionic sodiumdodecylsulfate (SDS) micelles. These systems can be considered as model mammalian and bacterial membrane interfaces, respectively. The goal of this study is to dissect the differences in peptide composition which make the mutant peptides (novispirin-G10 and novispirin-T7) less toxic than the parent peptide ovispirin (OVIS), although all three peptides have highly antibacterial properties. Compared to G10 and T7, OVIS inserts deepest into the DPC micelle. This correlates well with the lesser toxicity of G10 and T7. There is strong evidence which suggests that synergistic binding of hydrophobic residues drives binding of OVIS to the micelle. The helical content of G10 and T7 is reduced in the presence of DPC, and this leads to less amphipathic peptide structures, which bind weakly to the micelle. Simulations in SDS were carried out to compare the influence of membrane electrostatics on peptide structure. All three peptides bound strongly to SDS, and retained helical form. This corresponds well with their equally potent antibacterial properties. Based on the simulations, we argue that secondary structure stability often leads to toxic properties. We also propose that G10 and T7 operate by the carpet mechanism of cell lysis. Toxicity of peptides operating by the carpet mechanism can be attenuated by reducing the peptide helical content. The simulations successfully capture experimental binding states, and the different depths of binding of the three peptides to the two micelles correlate with their antibacterial and toxic properties.     

Solution conformation of a cyclic neurokinin antagonist: A nmr and molecular dynamics study

The solution structure of a hexapeptide, cyclo(Gln-Trp-Phe-Gly-Leu-Met), which is a selective NK-2 antagonist, has been studied by a combination of two-dimensional nmr and molecular dynamics (MD) techniques. The simulation based on nmr and MD data resulted in the convergence to a family of structures. Free molecular dynamics for 50 ps in the presence of DMSO solvent molecules shows that the structure is energetically stable. One intramolecular hydrogen bond between the amide proton of Gin and the carbonyl oxygen of Gly was revealed. This result is consistent with the results from the measurement of the temperature coefficient of the amide protons. The extent of intermolecular hydrogen bonding between the amide protons of the peptide and DMSO was also revealed by the free MD simulation. The resulting structure of the cyclic peptide contains a variation type I′ β-turn in the Gly-Leu-Met-Gln segment. Comparison of the structure of this peptide with that of other NK-2 antagonist cyclic hexapeptides was made, and the activity of cyclic antagonists appears to be inversely related to the conformational rigidity of the cyclic peptides. © 1994 John Wiley & Sons, Inc.     

X-Pro peptides: Synthesis and solution conformation of benzyloxycarbonyl-(Aib-Pro)4methyl ester. Evidence for a novel helical structure

The synthesis of the octapeptide, benzyloxycarbonyl-(α-aminoisobutyryl-L-prolyl)₄-methyl ester [Z-(Aib-Pro)₄-OMe] and an analysis of its solution conformation is reported. The octapeptide is shown to possess three strong intramolecular hydrogen bonds on the basis of studies of the solvent and temperature dependence of NH chemical shifts and rates of hydrogen–deuterium exchange. ¹³C studies are consistent with a structure involving only trans Aib-Pro bonds, while ir experiments support a hydrogen-bonded conformation. The Aib 3, 5, and 7 NH groups are shown to participate in hydrogen bonding. A 3₁₀ helical conformation compatible with the spectroscopic data is suggested. The proposed conformation consists of three type III β-turns with Aib and Pro at the corners and stabilized by 4 → 1 intramolecular hydrogen bonds.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Peptide folding driven by Van der Waals interactions

Contrary to the widespread view that hydrogen bonding and its entropy effect play a dominant role in protein folding, folding into helical and hairpin-like structures is observed in molecular dynamics (MD) simulations without hydrogen bonding in the peptide-solvent system. In the widely used point charge model, hydrogen bonding is calculated as part of the interaction between atomic partial charges. It is removed from these simulations by setting atomic charges of the peptide and water to zero. Because of the structural difference between the peptide and water, van der Waals (VDW) interactions favor peptide intramolecular interactions and are a major contributing factor to the structural compactness. These compact structures are amino acid sequence dependent and closely resemble standard secondary structures, as a consequence of VDW interactions and covalent bonding constraints. Hydrogen bonding is a short range interaction and it locks the approximate structure into the specific secondary structure when it is included in the simulation. In contrast to standard molecular simulations where the total energy is dominated by charge-charge interactions, these simulation results will give us a new view of the folding mechanism.     

Antimicrobial peptide RP-1 structure and interactions with anionic versus zwitterionic micelles

Topologically, platelet factor-4 kinocidins consist of distinct N-terminal extended, C-terminal helical, and interposing gamma-core structural domains. The C-terminal alpha-helices autonomously confer direct microbicidal activity, and the synthetic antimicrobial peptide RP-1 is modeled upon these domains. In this study, the structure of RP-1 was assessed using several complementary techniques. The high-resolution structure of RP-1 was determined by NMR in anionic sodium dodecyl sulfate (SDS) and zwitterionic dodecylphosphocholine (DPC) micelles, which approximate prokaryotic and eukaryotic membranes, respectively. NMR data indicate the peptide assumes an amphipathic alpha-helical backbone conformation in both micelle environments. However, small differences were observed in the side-chain orientations of lysine, tyrosine, and phenylalanine residues in SDS versus DPC environments. NMR experiments with a paramagnetic probe indicated differences in positioning of the peptide within the two micelle types. Molecular dynamics (MD) simulations of the peptide in both micelle types were also performed to add insight into the peptide/micelle interactions and to assess the validity of this technique to predict the structure of peptides in complex with micelles. MD independently predicted RP-1 to interact only peripherally with the DPC micelle, leaving its spherical shape intact. In contrast, RP-1 entered deeply into and significantly distorted the SDS micelle. Overall, the experimental and MD results support a preferential specificity of RP-1 for anionic membranes over zwitterionic membranes. This specificity likely derives from differences in RP-1 interaction with distinct lipid systems, including subtle differences in side chain orientations, rather than gross changes in RP-1 structure in the two lipid environments.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

Dielectric screening effect of electronic polarization and intramolecular hydrogen bonding

Recent site‐resolved hydrogen exchange measurements have uncovered significant discrepancies between simulations and experimental data during protein folding, including the excessive intramolecular hydrogen bonds in simulations. This finding indicates a possibility that intramolecular charge–charge interactions have not included sufficient dielectric screening effect of the electronic polarization. Scaling down peptide atomic charges according to the optical dielectric constant is tested in this study. As a result, the number of intramolecular hydrogen bonds is lower than using unscaled atomic charges while reaching the same levels of helical contents or β‐hairpin backbone hydrogen bonds, because van der Waals interactions contribute substantially to peptide folding in water. Reducing intramolecular charge–charge interactions and hydrogen bonding increases conformational search efficiency. In particular, it reduces the equilibrium helical content in simulations using AMBER force field and the energy barrier in folding simulations using CHARMM force field.     

Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.     

Conformational analysis of GT1B ganglioside and its interaction with botulinum neurotoxin type B: a study by molecular modeling and molecular dynamics

The conformational property of oligosaccharide GT1B in aqueous environment was studied by molecular dynamics (MD) simulation using all-atom model. Based on the trajectory analysis, three prominent conformational models were proposed for GT1B. Direct and water-mediated hydrogen bonding interactions stabilize these structures. The molecular modeling and 15 ns MD simulation of the Botulinum Neuro Toxin/B (BoNT/B) - GT1B complex revealed that BoNT/B can accommodate the GT1B in the single binding mode. Least mobility was seen for oligo-GT1B in the binding pocket. The bound conformation of GT1B obtained from the MD simulation of the BoNT/B-GT1B complex bear a close conformational similarity with the crystal structure of BoNT/A-GT1B complex. The mobility noticed for Arg 1268 in the dynamics was accounted for its favorable interaction with terminal NeuNAc. The internal NeuNAc1 tends to form 10 hydrogen bonds with BoNT/B, hence specifying this particular site as a crucial space for the therapeutic design that can restrict the pathogenic activity of BoNT/B.     

Study of the Alzheimer's Aβ40 peptide in SDS micelles using molecular dynamics simulations

The interaction of the Alzheimer's amyloid beta peptide, Aβ40, with sodium dodecyl sulfate (SDS) micelles, together with the self-assembly of SDS molecules around the peptide from an initial random distribution were studied using atomistic and coarse-grained (CG) molecular dynamics simulations. In atomistic simulations, the peptide structure in the micelle was characterized by two helical regions connected through a short hinge. The initial structure of the system was shown to affect the simulation results. The atomistic self-assembly of SDS molecules resulted in a 38-molecule micelle around the peptide, along with some globules and individual molecules. Coarse-grained simulation results, however, did not show such a difference, and at the end of all CG simulations, a complete 60-molecule micelle was obtained, with the peptide located at the interface of the micelle with water. The obtained CG radial density profiles and SDS micelle size and shape properties were identical for all CG simulations.