This review describes recent progress in the area of molecular simulations of peptide assemblies, including peptide-amphiphiles, and drug-amphiphiles. The ability to predict the structure and stability of peptide self-assemblies from the molecular level up is vital to the field of nanobiotechnology. Computational methods such as molecular dynamics offer the opportunity to characterize intermolecular forces between peptide-amphiphiles that are critical to the self-assembly process. Furthermore, these computational methods provide the ability to computationally probe the structure of these supramolecular assemblies at the molecular level, which is a challenge experimentally. Herein, we briefly highlight progress in the areas of all-atomistic and coarse-grained simulation studies investigating the self-assembly process of short peptides and peptide amphiphiles. We also discuss recent all-atomistic and coarse-grained simulations of the self-assembly of a drug-amphiphile into elongated filaments. Next, we discuss how these computational methods can provide further insight on the pathway of cylindrical nanofiber formation and predict their biocompatibility by studying the interaction of these peptide-amphiphile nanostructures with model cell membranes.

Peptide self-assembly has been used to design an array of nanostructures that possess functional biomedical applications. Experimental studies have reported nanofilament and nanotube formation from peptide-based drug amphiphiles (DAs). These DAs have shown to possess an inherently high drug loading with a tunable release mechanism. Herein, we report rational coarse-grained molecular dynamics simulations of the self-assembly process and the structure and stability of preassembled nanotubes at longer timescales (μs). We find that aggregation between these DAs at the submicrosecond timescale is driven by directional aromatic interactions between the drugs. The drugs form a large and high-density nucleus that is stable throughout microsecond timescales. Simulations of nanotubes characterize the drug−drug stacking and find correlations at nanometer length scales. These simulations can inform the rational molecular design of drug amphiphiles.

Accurate and efficient prediction of drug partitioning in model membranes is of significant interest to the pharmaceutical industry. Herein, we utilize advanced sampling methods, specifically, the adaptive biasing force methodology to calculate the potential of mean force for a model hydrophobic anticancer drug, camptothecin (CPT), across three model interfaces. We consider an octanol bilayer, a thick octanol/water interface, and a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/water interface. We characterize the enthalpic and entropic contributions of the drug to the potential of mean force. We show that the rotational entropy of the drug is inversely related to the probability of hydrogen bond formation of the drug with the POPC membrane. In addition, in long-time microsecond simulations of a high concentration of CPT above the POPC membrane, we show that strong drug–drug aromatic interactions shift the spatial orientation of the drug with the membrane. Stacks of hydrophobic drugs form, allowing penetration of the drug just under the POPC head groups. These results imply that inhomogeneous membrane models need to take into account the effect of drug aggregation on the membrane environment.

No abstract

The shape of drug delivery vehicles impacts both the circulation time and the effectiveness of the vehicle. Peptide-based drug amphiphiles (DAs) are promising new candidates as drug delivery vehicles that...

The pathway for supramolecular fiber formation is coupled with the underlying order of the self-assembling molecules. Here, we report on atomistic molecular dynamics simulations to characterize the initial stages of the self-assembly of a model drug amphiphile in an aqueous solution. We perform two-dimensional metadynamics calculations to characterize the assembly space of this model drug amphiphile�Tubustecan, TT1. TT1 is composed of the hydrophobic anticancer drug, Camptothecin (CPT), conjugated to a hydrophilic polyethylene glycol (PEG) chain. We find that the aromatic stacking of CPT drives the formation of a higherdensity liquid droplet. This droplet elongates and can form a higher-ordered supramolecular assembly upon reorganizing and forming an interface and additional aromatic stacking of the drugs. We show that novel reaction coordinates tailored to this class of molecules are essential in capturing the underlying degree of molecular order upon assembly. This approach can be refined and extended to characterize the supramolecular assembly pathway of other molecules containing aromatic compounds.

The nucleosome core particle (NCP) is the basic packaging unit of DNA. Recently reported structures of the NCP suggest that the histone octamer undergoes conformational changes during the process of DNA translocation around the histone octamer [1]. We perform long-time MD simulations of the nucleosome core particle in aqueous medium at varying salt concentrations to explore the unwrapping process [2]. Initial results show the formation of a loop/bulge of nucleosomal DNA after 2 microsecond simulation that continues to grow throughout our 5 microseconds simulation. Indeed, the loop structure agrees well with recently reported cryoEM studies of the distorted nucleosome core particle (NCP) [1]. Next, we explore with additional longer trajectories the process of nucleosome unwrapping at varying salt concentrations, comparing two nucleosomal DNA sequences, the alpha-satellite palindromic (ASP) nucleosomal DNA sequence and the Widom-601 sequence. We also probe the free energy profiles of nucleosome unwrapping that are associating with the ASP and Widom-601 sequence using Umbrella Sampling. Our results provide insight into the critical sequence-dependent barriers to nucleosome unwrapping.