This work deals with an analysis of the covalent functionalization of a carbon nanotube using polyethylene glycol chains terminated by folic acid fragments. The analysis is focused on theoretical predictions, using molecular dynamics simulations, of the properties of such constructs as pH controlled carriers of the anticancer drug doxorubicin. The analyzed systems are expected to hold the doxorubicin in the inner cavity of the carbon nanotube at neutral pH and unload the drug at slightly acidic pH. This property comes from incorporation into the nanotube of some dye molecules (p-phenylenediamine or neutral red) which undergo protonation at slightly acidic pH. We found that both dyes lead to the formation of a stable, co-absorbed phase of a doxorubicin-dye mixture inside the nanotube at physiological pH. At acidic pH we observed a spontaneous release of dyes from the nanotube, leading finally to the state with only doxorubicin encapsulated in the nanotube interior. Thus, the analyzed constructs can be considered as carriers of doxorubicin that are selective to tumor microenvironments (which exhibit reduced pH due to hypoxia and overexpression of folate receptors). However, we also found that the release of doxorubicin from the nanotube at acidic pH is kinetically blocked, at least in the case of the system sizes studied here. Thus, we also discussed some possible ways of reducing the activation barriers against doxorubicin release at acidic pH.

Molecular dynamics simulations were employed to study the properties of G-quadruplex and i-motif secondary DNA structures formed within the canonical telomere fragment of the Watson–Crick duplex. These secondary structures were built symmetrically in the same place of the duplex and were subjected to the analysis in standard unbiased simulations and using metadynamics scheme for the determination of potential of mean force associated with the enforced unfolding of the i-motif parts of the systems. Also, enforced formation of i-motif structures, starting from partially unfolded duplex, were studied in order to find whether formation of i-motif facilitates spontaneous formation of G-quadruplex. We found that i-motif formed from single stranded DNA is unstable at neutral pH and room temperature. On the other hand, the i-motif is strongly stabilized by the presence of complementary G-quadruplex, which should be the most likely configuration when these secondary structures form from double stranded DNA. The stabilization is observed either in neutral or in acidic pH though in the neutral case the i-motif can also reveal considerable stability in the hairpin configuration. We did not observe spontaneous folding of the guanine-rich strand into the G-quadruplex when the cytosine rich strand was dragged to i-motif configuration. This observation suggests that both folding and unfolding transitions are kinetically blocked.

We studied the permeability of graphene sheets with designed nanopores using the classical molecular dynamics. To characterize the energy profile for transmission we calculated the potential of the mean force. A high selectivity for methane + butane mixture with the hydrogen-passivated pore diameter of 0.32 nm was found where the volume exclusion mechanism governs the separation process. In the case of a slightly larger pore diameter of 0.64 nm the same alkane mixture separates completely unexpectedly: a larger butane molecule permeates much faster than a small methane one. The blocking effect of the permeation path by a larger mixture component when it worked like a cork was also observed. This is a promising perspective for using graphene to design intelligent membranes which can maintain a constant composition of mixtures in the permeable area.