In silico study of amphiphilic nanotubes based on cyclic peptides in polar and non-polar solvent

Hydrogen-Bond-Driven Peptide Nanotube Formation: A DFT Study

DFT calculations were carried out to examine geometries and binding energies of H-bond-driven peptide nanotubes. A bolaamphiphile molecule, consisting of two N-α amido glycylglycine head groups linked by either one CH2 group or seven CH2 groups, is used as a building block for nanotube self-assembly. In addition to hydrogen bonds between adjacent carboxy or amide groups, nanotube formation is also driven by weak C-H· · ·O hydrogen bonds between a methylene group and the carboxy OH group, and between a methylene group and an amide O=C group. The intratubular O-H· · ·O=C hydrogen bonds account for approximately a third of the binding energies. Binding energies calculated with the wB97XD/DGDZVP method show that the hydrocarbon chains play a stabilizing role in nanotube self-assembly. The shortest nanotube has the length of a single monomer and a diameter than increases with the number of monomers. Lengthening of the tubular structure occurs through intertubular O-H· · ·O=C hydrogen bonds. The average intertubular O-H· · ·O=C hydrogen bond binding energy is estimated to change with the size of the nanotubes, decreasing slightly towards some plateau value near 15 kcal/mol according to the wB97XD/DGDZVP method.

Influence of salt and temperature on the self-assembly of cyclic peptides in water: a molecular dynamics study

It is found in the literature that cyclic peptides (CPs) are able to self-assemble in water to form cyclic peptide nanotubes (CPNTs) and are used extensively in the field of nanotechnology. Several factors influence the formation and stability of these nanotubes in water. However, an extensive study of the contribution of several important factors is still lacking. The purpose of this study is to explore the effect of temperature and salt (NaCl) on the association tendency of CPs. Furthermore, the self-association behavior of CPs in aqueous solutions at various temperatures is also thoroughly discussed. Cyclo-[(Asp-D-Leu-Lys-D-Leu)₂] is considered for this study and a series of classical molecular dynamics (MD) simulations at three different temperatures, viz. 280 K, 300 K, and 320 K, both in pure water and in NaCl solutions of different concentrations are carried out. The calculations of radial distribution functions, preferential interaction parameters, cluster formation and hydrogen bonding properties suggest a strong influence of NaCl concentration on the association propensity of CPs. Low NaCl concentration hinders CP association while high NaCl concentration facilitates the association of CPs. Besides this, the association of CPs is found to be enhanced at low temperature. Furthermore, the thermodynamics of CP association is predominantly found to be enthalpy driven in both the presence and absence of salt. No crossover between enthalpy and entropy in CP association is observed. In addition, the MM-GBSA method is used to investigate the binding free energies of the CP rings that self-assembled to form nanotube like structures at all three temperatures.

Molecular Dynamics Simulations of Transmembrane Cyclic Peptide Nanotubes Using Classical Force Fields, Hydrogen Mass Repartitioning, and Hydrogen Isotope Exchange Methods: A Critical Comparison

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5–7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Pulling Simulations and Hydrogen Sorption Modelling on Carbon Nanotube Bundles

Recent progress in molecular simulation technology has developed an interest in modernizing the usual computational methods and approaches. For instance, most of the theoretical work on hydrogen adsorption on carbon nanotubes was conducted a decade ago. It should be insightful to reinvestigate the field and take advantage of code improvements and features implemented in contemporary software. One example of such features is the pulling simulation modules now available in many molecular dynamics programs. We conduct pulling simulations on pairs of carbon nanotubes and measure the inter-tube distance before they dissociate in water. We use this distance to set the interval size between adjacent nanotubes as we arrange them in bundle configurations. We consider bundles with triangular, intermediate and honeycomb patterns, and armchair nanotubes with a chiral index from n = 5 to n = 10. Then, we simulate low pressure hydrogen adsorption isotherms at 77 K, using the grand canonical Monte Carlo method. The different bundle configurations adsorb great hydrogen amounts that may exceed 2% wt at ambient pressures. The computed hydrogen capacities are considered large for physisorption on carbon nanostructures and attributed to the ultra-microporous network and extraordinary high surface area of the configured models.