Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations.

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Molecular Dynamics Simulations of Polyelectrolyte Complexes.

Polyelectrolyte complexes (PECs) are currently of great interest due to their applications toward developing new adaptive materials and their relevance in membraneless organelles. These complexes emerge during phase separation when oppositely charged polymers are mixed in aqueous media. Peptide-based PECs are particularly useful toward developing new drug delivery methods due to their inherent biocompatibility. The underlying peptide sequence can be tuned to optimize specific material properties of the complex, such as interfacial tension and viscosity. Given their applicability, it would be advantageous to understand the underlying sequence-dependent phase behavior of oppositely charged peptides. Here, we report microsecond molecular dynamic simulations to characterize the effect of hydrophobicity on the sequence-dependent peptide conformation for model polypeptide sequences that were previously reported by Tabandeh et al. These sequences are designed with alternating chirality of the peptide backbone. We present microsecond simulations of six oppositely charged peptide pairs, characterizing the sequence-dependent effect on peptide size, degree of hydrogen bonding, secondary structure, and conformation. This analysis recapitulates sensible trends in peptide conformation and degree of hydrogen bonding, consistent with experimentally reported results. Ramachandran plots reveal that backbone conformation at the single amino acid level is highly influenced by the neighboring sequence in the chain. These results give insight into how subtle changes in hydrophobic side chain size and chirality influence the strength of hydrogen bonding between the chains and, ultimately, the secondary structure. Furthermore, principal component analysis reveals that the minimum energy structures may be subtly modulated by the underlying sequence.

Exploration of the Nucleation Pathway for Supramolecular Fibers.

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 higher-density 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.

PEO-b-PBD Diblock Copolymers Induce Packing Defects in Lipid/Hybrid Membranes and Improve Insertion Rates of Natively Folded Peptides.

Hybrid membranes assembled from biological lipids and synthetic polymers are a promising scaffold for the reconstitution and utilization of membrane proteins. Recent observations indicate that inclusion of small fractions of polymer in lipid membranes can improve protein folding and function, but the exact structural and physical changes a given polymer sequence imparts on a membrane often remain unclear. Here, we use all-atom molecular dynamics simulations to study the structure of hybrid membranes assembled from DOPC phospholipids and PEO-b-PBD diblock copolymers. We verified our computational model using new and existing experimental data and obtained a detailed picture of the polymer conformations in the lipid membrane that we can relate to changes in membrane elastic properties. We find that inclusion of low polymer fractions induces transient packing defects into the membrane. These packing defects act as insertion sites for two model peptides, and in this way, small amounts of polymer content in lipid membranes can lead to large increases in peptide insertion rates. Additionally, we report the peptide conformational space in both pure lipid and hybrid membranes. Both membranes support similar alpha helical peptide structures, exemplifying the biocompatibility of hybrid membranes.

The Interaction of Supramolecular Anticancer Drug Amphiphiles with Phospholipid Membranes.

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 can self-assemble into shapes such as nanofilament and nanotube (diameter ~ 6-10 nm). The number of conjugated drugs affects the IC50 of these DAs, which is correlated to the effective cellular uptake. Characterizing and optimizing the interaction of these DAs and their assemblies with the cellular membrane is experimentally challenging. Long-time molecular dynamics can determine if the DA molecular structure affects the translocation across and interaction with the cellular membrane. Here, we report long-time atomistic simulation on Anton 2 (up to 25 µs) of these DAs with model cellular membranes. Results indicate that the interaction of these DAs with model cellular membranes is dependent on the number of conjugated drugs. We find that, with increased drug loading, the hydrophobic drug (camptothecin) builds up in the outer hydrophobic core of the membrane, pulling in positively charged peptide groups. Next, we computationally probe the interaction of differing shapes of these model drug delivery vehicles-nanofilament and nanotube-with the same model membranes, finding that the interaction of these nanostructures with the membrane is strongly repulsive. Results suggest that the hydrogen bond density between the nanostructure and the membrane may play a key role in modulating the interaction between the nanostructure and the membrane. Taken together, these results offer important insights for the rational design of peptide-based drug delivery vehicles.

Interaction of Camptothecin with Model Cellular Membranes.

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.

Rational Coarse-Grained Molecular Dynamics Simulations of Supramolecular Anticancer Nanotubes.

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.

Molecular Dynamics Simulations of Supramolecular Anticancer Nanotubes.

We report here on long-time all-atomistic molecular dynamics simulations of functional supramolecular nanotubes composed by the self-assembly of peptide-drug amphiphiles (DAs). These DAs have been shown to possess an inherently high drug loading of the hydrophobic anticancer drug camptothecin. We probe the self-assembly mechanism from random with ∼0.4 µs molecular dynamics simulations. Furthermore, we also computationally characterize the interfacial structure, directionality of π-π stacking, and water dynamics within several peptide-drug nanotubes with diameters consistent with the reported experimental nanotube diameter. Insight gained should inform the future design of these novel anticancer drug delivery systems.

Isomeric control of the mechanical properties of supramolecular filament hydrogels.

Supramolecular filament hydrogels are an emerging class of biomaterials that hold great promise for regenerative medicine, tissue engineering, and drug delivery. However, fine-tuning of their bulk mechanical properties at the molecular level without altering their network structures remains a significant challenge. Here we report an isomeric strategy to construct amphiphilic peptides through the conjugation of isomeric hydrocarbons to influence the local viscoelastic properties of their resulting supramolecular hydrogels. In this case, the packing requirements of the chosen isomeric hydrocarbons within the supramolecular filaments are dictated by their atomic arrangements at the molecular and intermolecular levels. Atomistic molecular dynamics simulations suggest that this design strategy can subtly alter the molecular packing at the interface between the peptide domain and the hydrophobic core of the supramolecular assemblies, without changing both the filament width and morphology. Our results from wide-angle X-ray scattering and molecular simulations further confirm that alterations to the intermolecular packing at the interface impact the strength and degree of hydrogen bonding within the peptide domains. This subtle difference in the isomeric hydrocarbon design and their consequent packing difference led to variations in the persistence length of the individual supramolecular filaments. Microrheological analysis reveals that this difference in filament stiffness enables the fine-tuning of the mechanical properties of the hydrogel at the macroscopic scale. We believe that this isomeric platform provides an innovative method to tune the local viscoelastic properties of supramolecular polymeric hydrogels without necessarily altering their network structures.

Coarse-grained molecular dynamics studies of the structure and stability of peptide-based drug amphiphile filaments.

Peptide-based supramolecular filaments, in particular filaments self-assembled by drug amphiphiles (DAs), possess great potential in the field of drug delivery. These filaments possess one hundred percent drug loading, with a release mechanism that can be tuned based on the dissociation of the supramolecular filaments and the degradation of the DAs [Cheetham et al., J. Am. Chem. Soc., 2013, 135(8), 2907]. Recently, much attention has been drawn to the competing intermolecular interactions that drive the self-assembly of peptide-based amphiphiles into supramolecular filaments. Recently, we reported on long-time atomistic molecular dynamics simulations to characterize the structure and growth of chiral filaments by the self-assembly of a DA containing the aromatic anti-cancer drug camptothecin [Kang et al., Macromolecules, 2016, 49(3), 994]. We found that the π-π stacking of the aromatic drug governs the early stages of the self-assembly process, while also contributing towards the chirality of the self-assembled filament. Based on these all-atomistic simulations, we now build a chemically accurate coarse-grained model that can capture the structure and stability of these supramolecular filaments at long time-scales (microseconds). These coarse-grained models successfully recapitulate the growth of the molecular clusters (and their elongation trends) compared with previously reported atomistic simulations. Furthermore, the interfacial structure and the helicity of the filaments are conserved. Next, we focus on characterization of the disassembly process of a 0.675 µm DA filament at microsecond time-scales. These results provide very useful tools for the rational design of functional supramolecular filaments, in particular supramolecular filaments for drug delivery applications.

Molecular simulations of peptide amphiphiles.

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 into the pathway of cylindrical nanofiber formation and predict their biocompatibility by studying the interaction of these peptide-amphiphile nanostructures with model cell membranes.

TCR triggering by pMHC ligands tethered on surfaces via poly(ethylene glycol) depends on polymer length.

Antigen recognition by T cells relies on the interaction between T cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) at the interface between the T cell and the antigen presenting cell (APC). The pMHC-TCR interaction is two-dimensional (2D), in that both the ligand and receptor are membrane-anchored and their movement is limited to 2D diffusion. The 2D nature of the interaction is critical for the ability of pMHC ligands to trigger TCR. The exact properties of the 2D pMHC-TCR interaction that enable TCR triggering, however, are not fully understood. Here, we altered the 2D pMHC-TCR interaction by tethering pMHC ligands to a rigid plastic surface with flexible poly(ethylene glycol) (PEG) polymers of different lengths, thereby gradually increasing the ligands' range of motion in the third dimension. We found that pMHC ligands tethered by PEG linkers with long contour length were capable of activating T cells. Shorter PEG linkers, however, triggered TCR more efficiently. Molecular dynamics simulation suggested that shorter PEGs exhibit faster TCR binding on-rates and off-rates. Our findings indicate that TCR signaling can be triggered by surface-tethered pMHC ligands within a defined 3D range of motion, and that fast binding rates lead to higher TCR triggering efficiency. These observations are consistent with a model of TCR triggering that incorporates the dynamic interaction between T cell and antigen-presenting cell.

Molecular simulation of the concentration-dependent interaction of hydrophobic drugs with model cellular membranes.

We report here the interactions between a hydrophobic drug and a model cellular membrane at the molecular level using all-atom molecular dynamics simulations of paclitaxel, a hydrophobic cancer drug. The calculated free energy of a single drug across the bilayer interface displays a minimum in the outer hydrophobic zone of the membrane. The transfer free energy shows excellent agreement with reported experimental data. In two sets of long-time simulations of high concentrations of drug in the membrane (12 and 11 mol %), the drugs display substantial clustering and rotation with significant directional preference in their diffusion. The main taxane ring partitions in the outer hydrophobic zone, while the three phenyl rings prefer to be closer to the hydrophobic core of the membrane. The clustering of the drug molecules, order parameters of the lipid tails, and water penetration suggest that the fluidity and permeability of the membrane are affected by the concentration of drugs that it contains. Furthermore, at the high-concentration limit, the free energy minimum shifts closer to the hydrophobic core and the central barrier to cross the membrane decreases. Moreover, the transfer free energy change substantially increases, suggesting that increasing concentration facilitates drug partitioning into the membrane.

Filomicelles in nanomedicine - from flexible, fragmentable, and ligand-targetable drug carrier designs to combination therapy for brain tumors.

Nanoparticles that are made by self-assembly into non-spherical shapes are promising as drug delivery vehicles. This review focuses on flexible and fragmentable filamentous micelles referred to as filomicelles made of degradable block copolymer amphiphiles. They are inspired by filoviruses and also by tubular proplatelets that break up into smaller platelets in blood flow. The synthesis and assembly of the constituent block copolymers are described together with ligand targeting and fragmentation as well as drug release in therapeutic applications to model tumors and most recently brain tumors.