Different Drug Mobilities in Hydrophobic Cavities of Host-Guest Complexes between ß-Cyclodextrin and 5-Fluorouracil at Different Stoichiometries: A Molecular Dynamics Study in Water.

Cyclodextrins (CDs) are cyclic oligosaccharides able to form noncovalent water-soluble complexes useful in many different applications for the solubilization, delivery, and greater bioavailability of hydrophobic drugs. The complexation of 5-fluorouracil (5-FU) with natural or synthetic cyclodextrins permits the solubilization of this poorly soluble anticancer drug. In this theoretical work, the complexes between ß-CD and 5-FU are investigated using molecular mechanics (MM) and molecular dynamics (MD) simulations in water. The inclusion complexes are formed thanks to the favorable intermolecular interactions between ß-CD and 5-FU. Both 1:1 and 1:2 ß-CD/5-FU stoichiometries are investigated, providing insight into their interaction geometries and stability over time in water. In the 1:2 ß-CD/5-FU complexes, the intermolecular interactions affect the drug's mobility, suggesting a two-step release mechanism: a fast release for the more exposed and hydrated drug molecule, with greater freedom of movement near the ß-CD rims, and a slow one for the less-hydrated and well-encapsulated and confined drug. MD simulations study the intermolecular interactions between drugs and specific carriers at the atomistic level, suggesting a possible release mechanism and highlighting the role of the impact of the drug concentration on the kinetics process in water. A comparison with experimental data in the literature provides further insights.

Adsorption Affinities of Small Volatile Organic Molecules on Graphene Surfaces for Novel Nanofiller Design: A DFT Study.

The adsorption of organic molecules on graphene surfaces is a crucial process in many different research areas. Nano-sized carbon allotropes, such as graphene and carbon nanotubes, have shown promise as fillers due to their exceptional properties, including their large surface area, thermal and electrical conductivity, and potential for weight reduction. Surface modification methods, such as the "pyrrole methodology", have been explored to tailor the properties of carbon allotropes. In this theoretical work, an ab initio study based on Density Functional Theory is performed to investigate the adsorption process of small volatile organic molecules (such as pyrrole derivatives) on graphene surface. The effects of substituents, and different molecular species are examined to determine the influence of the aromatic ring or the substituent of pyrrole's aromatic ring on the adsorption energy. The number of atoms and presence of π electrons significantly influence the corresponding adsorption energy. Interestingly, pyrroles and cyclopentadienes are 10 kJ mol-1 more stable than the corresponding unsaturated ones. Pyrrole oxidized derivatives display more favorable supramolecular interactions with graphene surface. Intermolecular interactions affect the first step of the adsorption process and are important to better understand possible surface modifications for carbon allotropes and to design novel nanofillers in polymer composites.

Adsorption and Self-Aggregation of Chiral [5]-Aza[6]helicenes on DNA Architecture: A Molecular Dynamics Study.

Helicenes are an extremely interesting class of conjugated molecules without asymmetric carbon atoms but with intrinsic chirality. These molecules can interact with double-stranded chiral B-DNA architecture, modifying after their adsorption the hydrophilicity exposed by DNA to the biological environment. They also form ordered structures due to self-aggregation processes with possible different light emissions. Following initial studies based on molecular mechanics (MM) and molecular dynamics (MD) simulations regarding the adsorption and self-aggregation process of 5-aza[5]helicenes on double-stranded B-DNA, this theoretical work investigates the interaction between (M)- and (P)-5-aza[6]helicenes with double-helix DNA. Initially, the interaction of the pure single enantiomer with DNA is studied. Possible preferential absorption in minor or major grooves can occur. Afterward, the interaction of enantiopure compounds (M)- and (P)-5-aza[6]helicenes, potentially occurring in a racemic mixture at different concentrations, was investigated, taking into consideration both competitive adsorption on DNA and the possible helicenes' self-aggregation process. The structural selectivity of DNA binding and the role of helicene concentration in adsorption and the self-aggregation process are interesting. In addition, the ability to form ordered structures on DNA that follow its chiral architecture, thanks to favorable van der Waals intermolecular interactions, is curious.

Molecular modeling and phylogenetic analyses highlight the role of amino acid 347 of the N1 subtype neuraminidase in influenza virus host range and interspecies adaptation.

The N1 neuraminidases (NAs) of avian and pandemic human influenza viruses contain tyrosine and asparagine, respectively, at position 347 on the rim of the catalytic site; the biological significance of this difference is not clear. Here, we used molecular dynamics simulation to model the effects of amino acid 347 on N1 NA interactions with sialyllacto-N-tetraoses 6'SLN-LC and 3'SLN-LC, which represent NA substrates in humans and birds, respectively. Our analysis predicted that Y347 plays an important role in the NA preference for the avian-type substrates. The Y347N substitution facilitates hydrolysis of human-type substrates by resolving steric conflicts of the Neu5Ac2-6Gal moiety with the bulky side chain of Y347, decreasing the free energy of substrate binding, and increasing the solvation of the Neu5Ac2-6Gal bond. Y347 was conserved in all N1 NA sequences of avian influenza viruses in the GISAID EpiFlu database with two exceptions. First, the Y347F substitution was present in the NA of a specific H6N1 poultry virus lineage and was associated with the substitutions G228S and/or E190V/L in the receptor-binding site (RBS) of the hemagglutinin (HA). Second, the highly pathogenic avian H5N1 viruses of the Gs/Gd lineage contained sporadic variants with the NA substitutions Y347H/D, which were frequently associated with substitutions in the HA RBS. The Y347N substitution occurred following the introductions of avian precursors into humans and pigs with N/D347 conserved during virus circulation in these hosts. Comparative evolutionary analysis of site 347 revealed episodic positive selection across the entire tree and negative selection within most host-specific groups of viruses, suggesting that substitutions at NA position 347 occurred during host switches and remained under pervasive purifying selection thereafter. Our results elucidate the role of amino acid 347 in NA recognition of sialoglycan substrates and emphasize the significance of substitutions at position 347 as a marker of host range and adaptive evolution of influenza viruses.

L-Arginine-Derived Polyamidoamine Oligomers Bearing at Both Ends ß-Cyclodextrin Units as pH-Sensitive Curcumin Carriers.

The aza-Michael polyaddition of L-arginine and N,N'-methylene-bis-acrylamide gives the biocompatible and easily cell-internalized polyamidoamine ARGO7. By controlled synthesis, two ARGO7 oligomers, namely a trimer and a pentamer, bearing acrylamide terminal units, were obtained as precursors of the ß-cyclodextrin-end-terminated oligomers P3 and P5, which have been shown to encapsulate curcumin at both pH 7.4 and 4.5. After lyophilization, P3- and P5-curcumin complexes gave stable water solutions. The apparent solubility of encapsulated curcumin was in the range 20-51 µg mL-1, that is, three orders of magnitude higher than the water solubility of free curcumin (0.011 µg mL-1). The drug release profiles showed induction periods both at pH levels 4.5 and 7.4, suggesting a diffusive release mechanism, as confirmed by kinetic studies. The release rate of curcumin was higher at pH 7.4 than at pH 4.5 and, in both cases, it was higher for the P5 complex. Encapsulated curcumin was more photostable than the free drug. Molecular mechanics and molecular dynamics simulations explain at atomistic level the formation of aggregates due to favorable van der Waals interactions. The drug molecules interact with the external surface of carriers or form inclusion complexes with the ß-cyclodextrin cavities. The aggregate stability is higher at pH 4.5.

Surface Interactions between Ketoprofen and Silica-Based Biomaterials as Drug Delivery System Synthesized via Sol-Gel: A Molecular Dynamics Study.

Biomaterial-based drug delivery systems for a controlled drug release are drawing increasing attention thanks to their possible pharmaceutical and biomedical applications. It is important to control the local administration of drugs, especially when the drug exhibits problems diffusing across biological barriers. Thus, in an appropriate concentration, it would be released in situ, reducing side effects due to interactions with the biological environment after implantation. A theoretical study based on Molecular Mechanics and Molecular Dynamics methods is performed to investigate possible surface interactions between the amorphous SiO2 surface and the ketoprofen molecules, an anti-inflammatory drug, considering the role of drug concentration. These theoretical results are compared with experimental data obtained by analyzing, through Fourier transform infrared spectroscopy (FT-IR), the interaction between the SiO2 amorphous surface and two percentages of the ketoprofen drug entrapped in a silica matrix obtained via the sol-gel method and dried materials. The loaded drug in these amorphous bioactive material forms hydrogen bonds with the silica surface, as found in this theoretical study. The surface interactions are essential to have a new generation of biomaterials not only important for biocompatibility, with specific structural and functional properties, but also able to incorporate anti-inflammatory agents for release into the human body.

Adsorption of Chiral [5]-Aza[5]helicenes on DNA Can Modify Its Hydrophilicity and Affect Its Chiral Architecture: A Molecular Dynamics Study.

Helicenes are interesting chiral molecules without asymmetric carbon atoms but with intrinsic chirality. Functionalized 5-Aza[5]helicenes can form non-covalent complexes with anticancer drugs and therefore be potential carriers. The paper highlights the different structural selectivity for DNA binding for two enantiopure compounds and the influence of concentration on their adsorption and self-aggregation process. In this theoretical study based on atomistic molecular dynamics simulations the interaction between (M)- and (P)-5-Aza[5]helicenes with double helix B-DNA is investigated. At first the interaction of single pure enantiomer with DNA is studied, in order to find the preferred site of interaction at the major or minor groove. Afterwards, the interaction of the enantiomers at different concentrations was investigated considering both competitive adsorption on DNA and possible helicenes self-aggregation. Therefore, racemic mixtures were studied. The helicenes studied are able to bind DNA modulating or locally modifying its hydrophilic surface into hydrophobic after adsorption of the first helicene layer partially covering the negative charge of DNA at high concentration. The (P)-enantiomer shows a preferential binding affinity of DNA helical structure even during competitive adsorption in the racemic mixtures. These DNA/helicenes non-covalent complexes exhibit a more hydrophobic exposed surface and after self-aggregation a partially hidden DNA chiral architecture to the biological environment.

Understanding Surface Interaction and Inclusion Complexes between Piroxicam and Native or Crosslinked ß-Cyclodextrins: The Role of Drug Concentration.

Drug concentration plays an important role in the interaction with drug carriers affecting the kinetics of release process and toxicology effects. Cyclodextrins (CDs) can solubilize hydrophobic drugs in water enhancing their bioavailability. In this theoretical study based on molecular mechanics and molecular dynamics methods, the interactions between ß-cyclodextrin and piroxicam, an important nonsteroidal anti-inflammatory drug, were investigated. At first, both host-guest complexes with native ß-CD in the 1:1 and in 2:1 stoichiometry were considered without assuming any initial a priori inclusion: the resulting inclusion complexes were in good agreement with literature NMR data. The interaction between piroxicam and a ß-CD nanosponge (NS) was then modeled at different concentrations. Two inclusion mechanisms were found. Moreover, piroxicam can interact with the external NS surface or with its crosslinkers, also forming one nanopore. At larger concentration, a nucleation process of drug aggregation induced by the first layer of adsorbed piroxicam molecules is observed. The flexibility of crosslinked ß-CDs, which may be swollen or quite compact, changing the surface area accessible to drug molecules, and the dimension of the aggregate nucleated on the NS surface are important factors possibly affecting the kinetics of release, which shall be theoretically studied in more detail at specific concentrations.

Dendrimer Dynamics: A Review of Analytical Theories and Molecular Simulation Methods.

The theoretical study of dendrimers is reviewed, considering both analytical approaches and molecular simulation methods. We discuss the effect of molecular symmetry on the degeneracy of the relaxation times, and then the calculation of observable quantities, in particular the intrinsic viscosity, and then the viscoelastic complex modulus and the dynamic structure factor, in comparison with the available experimental data. In particular, the maximum intrinsic viscosity with increasing molar mass is analyzed in some detail. The approximations and/or assumptions of the adopted methods are also described in connection with analogous results for polymer of a different topology, in particular linear and star polymers.

Hydrogen Bonding in a l-Glutamine-Based Polyamidoamino Acid and its pH-Dependent Self-Ordered Coil Conformation.

This paper reports on synthesis, acid-base properties, and self-structuring in water of a chiral polyamidoamino acid, M-l-Gln, obtained from the polyaddition of N,N'-methylenebisacrylamide with l-glutamine, with the potential of establishing hydrogen bonds through its prim-amide pendants. The M-l-Gln showed pH-responsive circular dichroism spectra, revealing ordered conformations. Structuring was nearly insensitive to ionic strength but sensitive to denaturing agents. The NMR diffusion studies were consistent with a population of unimolecular nanoparticles thus excluding aggregation. The M-l-Gln had the highest molecular weight and hydrodynamic radius among all polyamidoamino acids described. Possibly, transient hydrogen bonds between l-glutamine molecules and M-l-Gln growing chains facilitated the polyaddition reaction. Theoretical modeling showed that M-l-Gln assumed pH-dependent self-ordered coil conformations with main chain transoid arrangements reminiscent of the protein hairpin motif owing to intramolecular dipole moments and hydrogen bonds. The latter were most numerous at the isoelectric point (pH 4.5), where they mainly involved even topologically distant main chain amide N-H and side chain amide C=O brought to proximity by structuring. Hydrogen bonds at pH 4.5 were also suggested by variable temperature NMR. The 2D NOESY experiments at pH 4.5 confirmed the formation of compact structures through the analysis of the main chain/side chain hydrogen contacts, in line with MD simulations.

A Molecular Dynamics Study of a Photodynamic Sensitizer for Cancer Cells: Inclusion Complexes of γ-Cyclodextrins with C70.

Photodynamic therapy is an emerging treatment of tumor diseases. The complexes with γ-cyclodextrins (γ-CD) and fullerenes or their derivatives can be used as photosensitizers by direct injection into cancer cells. Using molecular mechanics and molecular dynamics methods, the stability and the geometry of the 2:1 complexes [(γ-CD)2/C70] are investigated analyzing the differences with the analogous C60 complexes, studied in a previous theoretical work and experimentally found to be much less efficient in cancer therapy. The inclusion complex of γ-CD and C70 has a 2:1 stoichiometry, the same as C60, but is significantly less stable and displays an unlike arrangement. In vacuo, mimicking an apolar solvent, the complex is compact, whereas in water the two γ-CDs encapsulate C70 forming a relatively stable complex by interacting through their primary rims, however exposing part of C70 to the solvent. Other higher-energy complexes with the γ-CDs facing different rims can form in water, but in all cases part of the hydrophobic C70 surface remains exposed to water. The stability and arrangement of these peculiar amphiphilic inclusion complexes having non-covalent interactions in water can be an important key for cancer therapy to enhance both the solubilization and the fullerene insertion into liposomes or cell membranes.

Self-Structuring in Water of Polyamidoamino Acids with Hydrophobic Side Chains Deriving from Natural α-Amino Acids.

This paper reports on synthesis, acid-base properties and self-structuring in water of chiral polyamidoamino acids (PAACs) obtained by polyaddition of N,N'-methylenebisacrylamide with l-alanine, l-valine and l-leucine (M-l-Ala, M-l-Val, M-l-Leu) with potential for selective interactions with biomolecules. The polymers maintained the acid-base properties of amino acids. In water, the circular dichroism spectra of PAACs revealed pH-dependent structuring in the range 3⁻11 and in the wavelength interval 200⁻280 nm. Taking as reference the values at pH 3, the differential molar ellipticities were plotted in the pH interval 3⁻11. Sigmoidal curves were obtained presenting inflection points at pH 8.1, 6.8 and 7.3 for M-l-Ala, M-l-Val and M-l-Leu, respectively, corresponding to the amine half-ionization. Theoretical modeling showed that PAACs assumed stable folded conformations. Intramolecular interactions led to transoid arrangements of the main chain reminiscent of protein hairpin motif. Oligomers with ten repeat units had simulated gyration radii consistent with the hydrodynamic radii obtained by dynamic light scattering.

Inclusion complexes of ß-cyclodextrin with tricyclic drugs: an X-ray diffraction, NMR and molecular dynamics study.

Tricyclic fused-ring cyclobenzaprine (1) and amitriptyline (2) form 1:1 inclusion complexes with ß-cyclodextrin (ß-CD) in the solid state and in water solution. Rotating frame NOE experiments (ROESY) showed the same geometry of inclusion for both 1/ß-CD and 2/ß-CD complexes, with the aromatic ring system entering the cavity from the large rim of the cyclodextrin and the alkylammonium chain protruding out of the cavity and facing the secondary OH rim. These features matched those found in the molecular dynamics (MD) simulations in solution and in the solid state from single-crystal X-ray diffraction of 1/ß-CD and 2/ß-CD complexes. The latter complex was found in a single conformation in the solid state, whilst the MD simulations in explicit water reproduced the conformational transitions observed experimentally for the free molecule.

Self-Ordering Secondary Structure of d- and l-Arginine-Derived Polyamidoamino Acids.

This paper reports on synthesis, acid-base properties and pH-dependent structuring in water of d-, l- and d,l-ARGO7, bioinspired polymers obtained by polyaddition of the corresponding arginine stereoisomers with N,N'-methylenebis(acrylamide). The circular dichroism spectra of d- and l-ARGO7 showed a peak at 228 nm and quickly and reversibly responded to pH changes, but were nearly unaffected by temperature, ionic strength, and denaturating agents. Theoretical modeling studies of L-ARGO7 showed that it assumed a folded structure. Intramolecular interactions led to transoid arrangements of the main chain reminiscent of the protein hairpin motif. Torsion angles showed a quite similar distribution at pH 6 and 14 consistent with the similarity of the CD spectra from pH 6 upward.

Aggregation behavior of amphiphilic cyclodextrins in a nonpolar solvent: evidence of large-scale structures by atomistic molecular dynamics simulations and solution studies.

Chemically modified cyclodextrins carrying both hydrophobic and hydrophilic substituents may form supramolecular aggregates or nanostructures of great interest. These systems have been usually investigated and characterized in water for their potential use as nanocarriers for drug delivery, but they can also aggregate in apolar solvents, as shown in the present paper through atomistic molecular dynamics simulations and dynamic light scattering measurements. The simulations, carried out with a large number of molecules in vacuo adopting an unbiased bottom-up approach, suggest the formation of bidimensional structures with characteristic length scales of the order of 10 nm, although some of these sizes are possibly affected by the assumed periodicity of the simulation cell, in particular at longer lengths. In any case, these nanostructures are stable at least from the kinetic viewpoint for relatively long times thanks to the large number of intermolecular interactions of dipolar and dispersive nature. The dynamic light scattering experiments indicate the presence of aggregates with a hydrodynamic radius of the order of 80 nm and a relatively modest polydispersity, even though smaller nanometer-sized aggregates cannot be fully ruled out. Taken together, these simulation and experimental results indicate that amphiphilically modified cyclodextrins do also form large-scale nanoaggregates even in apolar solvents.

Atomistic simulation of hydrophobin HFBII conformation in aqueous and fluorous media and at the water/vacuum interface.

Hydrophobins are proteins of interest for numerous applications thanks to their unique conformational and surface properties and their ability to self-assemble at interfaces. Here we report fully atomistic molecular mechanics and molecular dynamics results together with circular dichroism experimental data, aimed to study the conformational properties of the hydrophobin HFBII in a fluorinated solvent in comparison with a water solution and/or at an aqueous/vacuum interface. Both the atomistic simulations and the circular dichroism data show the remarkable structural stability of HFBII at all scales in all these environments, with no significant structural change, although a small cavity is formed in the fluorinated solvent. The combination of theoretical calculations and circular dichroism data can describe in detail the protein conformation and flexibility in different solvents and/or at an interface, and constitutes a first step towards the study of their self-assembly.

Separation of chiral nanotubes with an opposite handedness by chiral oligopeptide adsorption: A molecular dynamics study.

The separation of enantiomeric chiral nanotubes that can form non-covalent complexes with an unlike stability upon adsorption of chiral molecules is a process of potential interest in different fields and applications. Using fully atomistic molecular dynamics simulations, we report in this paper a theoretical study of the adsorption and denaturation of an oligopeptide formed by 16 chiral amino acids having a helical structure in the native state on both the inner and the outer surface of the chiral (10, 20) and (20, 10) single-walled carbon nanotubes having an opposite handedness, and of the armchair (16, 16) nanotube with a similar diameter for comparison. In the final adsorbed state, the oligopeptide loses in all cases its native helical conformation, assuming elongated geometries that maximize its contact with the surface through all the 16 amino acids. We find that the complexes formed by the two chiral nanotubes and the chosen oligopeptide have a strongly unlike stability both when adsorption takes place on the outer convex surface of the nanotube, and when it occurs on the inner concave surface. Thus, our molecular simulations indicate that separation of chiral, enantiomeric carbon nanotubes for instance by chromatographic methods can indeed be carried out using oligopeptides of a sufficient length.

Aggregation behaviour of amphiphilic cyclodextrins: the nucleation stage by atomistic molecular dynamics simulations.

Amphiphilically modified cyclodextrins may form various supramolecular aggregates. Here we report a theoretical study of the aggregation of a few amphiphilic cyclodextrins carrying hydrophobic thioalkyl groups and hydrophilic ethylene glycol moieties at opposite rims, focusing on the initial nucleation stage in an apolar solvent and in water. The study is based on atomistic molecular dynamics methods with a "bottom up" approach that can provide important information about the initial aggregates of few molecules. The focus is on the interaction pattern of amphiphilic cyclodextrin (aCD), which may interact by mutual inclusion of the substituent groups in the hydrophobic cavity of neighbouring molecules or by dispersion interactions at their lateral surface. We suggest that these aggregates can also form the nucleation stage of larger systems as well as the building blocks of micelles, vesicle, membranes, or generally nanoparticles thus opening new perspectives in the design of aggregates correlating their structures with the pharmaceutical properties.

Understanding the topography effects on competitive adsorption on a nanosized anatase crystal: a molecular dynamics study.

Using Molecular Mechanics and Molecular Dynamics methods, we investigated at the atomistic level the topography effects both on physisorption on different crystalline planes of TiO2 anatase and on the competitive adsorption when three different crystallographic faces were simultaneously present in an idealized nanosized crystal interacting with a simple heteroaromatic molecule experimentally used in sunlight-induced photosynthetic reaction.

Surface topography effects in protein adsorption on nanostructured carbon allotropes.

We report a molecular dynamics (MD) simulation study of protein adsorption on the surface of nanosized carbon allotropes, namely single-walled carbon nanotubes (SWNT) considering both the convex outer surface and the concave inner surface, together with a graphene sheet for comparison. These systems are chosen to investigate the effect of the surface curvature on protein adsorption at the same surface chemistry, given by sp(2) carbon atoms in all cases. The simulations show that proteins do favorably interact with these hydrophobic surfaces, as previously found on graphite which has the same chemical nature. However, the main finding of the present study is that the adsorption strength does depend on the surface topography: in particular, it is slightly weaker on the outer convex surfaces of SWNT and is conversely enhanced on the inner concave SWNT surface, being therefore intermediate for flat graphene. We additionally find that oligopeptides may enter the cavity of common SWNT, provided their size is small enough and the tube diameter is large enough for both entropic and energetic reasons. Therefore, we suggest that proteins can effectively be used to solubilize in water single-walled (and by analogy also multiwalled) carbon nanotubes through adsorption on the outer surface, as indeed experimentally found, and to functionalize them after insertion of oligopeptides within the cavity of nanotubes of appropriate size.