Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations.

With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

Activation Pathways and Free Energy Landscapes of the SARS-CoV-2 Spike Protein.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses a spike protein (S-protein) to recognize the receptor protein ACE2 of human cells and initiate infection, during which the conformational transition of the S-protein from inactive (down) state to active (up) state is one of the key molecular events determining the infectivity but the underlying mechanism remains poorly understood. In this work, we investigated the activation pathways and free energy landscape of the S-protein of SARS-CoV-2 and compared with those of the closely related counterpart SARS-CoV using molecular dynamics simulations. Our results revealed a large difference between the activation pathways of the two S-proteins. The transition from inactive to an active state for the S-protein of SARS-CoV-2 is more cooperative, involving simultaneous disruptions of several key interfacial hydrogen bonds, and the transition encounters a much higher free energy barrier. In addition, the conformational equilibrium of the SARS-CoV-2 S-protein is more biased to the inactive state compared to that of the SARS-CoV S-protein, suggesting the transient feature of the active state before binding to the receptor protein of the host cell. The key interactions contributing to the difference of the activation pathways and free energy landscapes were discussed. The results provide insights into the molecular mechanism involved in the initial stage of the SARS-CoV-2 infection.

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding.

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.

Chirality-Dependent Adsorption between Amphipathic Peptide and POPC Membrane.

The interactions between chiral molecules and cell membranes have attracted more and more attention in recent decades, due to their importance in molecular science and medical applications. It is observed that some peptides composed of different chiral amino acids may have distinct interactions with a membrane. How does the membrane exhibit a selective behavior related to the chirality of the peptides? Microscopically, the interactions between the peptides and the membrane are poorly understood. In this work, we study the interactions between an amphipathic peptide (C6) and POPC membrane with simulations. The kinetics and thermodynamics of peptide enantiomers during the adsorption to the membrane are characterized with direct simulations and umbrella sampling. It is observed that there are slow kinetics for the peptide composed of D-type amino acids. Along the observed pathways, the free energy landscapes are determined with umbrella sampling techniques. A free-energy barrier for the peptide composed of D-amino acids is observed, which is consistent with the kinetic observations. The results indicate the concurrent adsorption and rotation of the peptide helix. The local interactions between the peptides and the membrane are examined in detail, including the contact interactions between the peptides and the membrane, and the distributions of the lipids around the peptide. There are observable differences of the local interactions for the cases related to different peptide enantiomers. These results further demonstrate the importance of the rotation of peptide helix during the adsorption. More interestingly, all these kinetic differences between peptide enantiomers can be explained based on the conformations of the residue Trp and interactions between Trp and lipid molecules. These results give us a molecular understanding of the mechanism of the chirality-dependent peptide-membrane interactions, and may provide clues to designing systems which are sensitive to the chirality of membranes.

Surface-assisted assembly of a histidine-rich lipidated peptide for simultaneous exfoliation of graphite and functionalization of graphene nanosheets.

Biological molecules have promising potential to exfoliate graphite and produce biocompatible graphene nano-materials for biomedical applications. Here, a systematic design of a histidine-rich lipidated peptide sequence is presented that simultaneously exfoliates graphite flakes and functionalizes the resulting graphene nanosheets (∼150 nm lateral size) with long-term dispersion stability in aqueous solution (>8 months). The details of peptide/peptide and peptide/graphite interactions are probed using various microscopy, spectroscopy and molecular dynamics simulation methods. The results show that histidine and stearic acid interact with the graphite surface through π-π stacking and hydrophobic forces, respectively. Surface-assisted assembly of peptide molecules is then initiated via hydrogen bonds between deprotonated histidine segments, and a textured peptide nano-structure is formed. The work of adhesion between the peptide and graphite is found to be high enough to promote exfoliation of graphite flakes through layer-by-layer peeling of graphene nanosheets. The positively charged arginine in the peptide is exposed outward, and is responsible for the stable dispersion. The peptide molecules are sufficiently small, presenting the possibility to insert into and increase the spacing between the graphitic layers for enhanced exfoliation. The peptide-functionalized graphene nanosheets not only show great biocompatibility with cells in vitro, but also enhance cancer drug uptake by the cells.

Functional Control of Peptide Amphiphile Assemblies via Modulation of Internal Cohesion and Surface Chemistry Switch.

Understanding the impacts of the internal cohesion and surface chemistry of supramolecular systems on the collective behaviors in the contacts between the systems and biomolecules can greatly expand the functional diversity and adaptivity of supramolecular nanostructures. Here we show how the tuned molecular interactions modulate the morphologies and internal cohesion of peptide amphiphile (PA) self-assemblies and their resultant functions. Circular dichroism spectroscopy, fluorescence probing, atomic force and electron microscopy, along with molecular dynamics simulations, revealed that the PA self-assembly formed compact long fibers when surface charge repulsion was screened, but formed loose short fibers or micelle-like assemblies when hydrogen bonding was disrupted or hydrophobic core was enhanced. More importantly, depending on the strength of the phospholipid affinity for the cationic segment of the PA, the same internal cohesion of PA nanostructures can lead to either cell death or cell survival, providing unique insights into the design of supramolecular materials.

Atomistic picture for the folding pathway of a hybrid-1 type human telomeric DNA G-quadruplex.

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and K(+) ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.

Molecular binding of self-assembling peptide EAK16-II with anticancer agent EPT and its implication in cancer cell inhibition.

The current drug delivery techniques involve encapsulation, targeting and controlled release of the drug with various molecules or nanoparticles, but rarely has the drug molecular state or form been investigated. It is necessary to deliver a drug with a prescribed molecular state in order to maximize drug therapeutic effects. Here we present two facile methods to characterize molecular states of the anticancer drug ellipticine (EPT) encapsulated in the self-assembling peptide EAK, and relate the different molecular states of EPT to their respective cancer inhibition efficacies. The first method is UV-based, where drug loading capacity of a particular molecular state was determined. The experimental data corroborated a molecular binding model, where peptide-drug interaction was assumed to be electrostatic in nature. The developed model could elucidate a unique pH effect on protonated EPT loading capacity. The second method is based on fluorescence characteristics of EPT, which could differentiate the two molecular states: protonated and crystalline of EPT in situ. The inner filter effect was, however, found with this method, presenting an ineluctable obstacle in quantitative analysis of fluorescence data. A correction method for the inner filter effect was thus developed. With this approach, concentrations of EPT at different molecular states in its peptide complex solutions were determined. In vitro cytotoxicity assay was applied to evaluate the efficacy of the two molecular states of EPT, showing that protonated EPT was more efficient at killing cancer cells than crystalline EPT. The molecular binding model and two characterization methods for EAK-EPT complexation could be extended to other carrier-drug systems.

Interaction of an ionic complementary peptide with a hydrophobic graphite surface.

Protein adsorption on a surface plays an important role in biomaterial science and medicine. It is strongly related to the interaction between the protein residues and the surface. Here we report all-atom molecular dynamics simulations of the adsorption of an ionic complementary peptide, EAK16-II, to the hydrophobic highly ordered pyrolytic graphite surface. We find that, the hydrophobic interaction is the main force to govern the adsorption, and the peptide interchain electrostatic interaction affects the adsorption rate. Under neutral pH condition, the interchain electrostatic attraction facilitates the adsorption, whereas under acidic and basic conditions, because of the protonation and deprotonation of glutamic acid and lysine residues, respectively, the resulting electrostatic repulsion slows down the adsorption. We also found that under basic condition, during the adsorption peptide Chain II will be up against a choice to adsorb to the surface through the hydrophobic interaction or to form a temporary hydrophobic core with the deposited peptide Chain I. These results provide a basis for understanding some of the fundamental interactions governing peptide adsorption on the surface, which can shed new light on novel applications, such as the design of implant devices and drug delivery materials.

Tension at the surface: which phase is more important, liquid or vapor?

Tension at the surface is a most fundamental physicochemical property of a liquid surface. The concept of surface tension has widespread implications in numerous natural, engineering and biomedical processes. Research to date has been largely focused on the liquid side; little attention has been paid to the vapor--the other side of the surface, despite over 100 years of study. However, the question remains as to whether the vapor plays any role, and to what extent it affects the surface tension of the liquid. Here we show a systematic study of the effect of vapor on the surface tension and in particular, a surprising observation that the vapor, not the liquid, plays a dominant role in determining the surface tension of a range of common volatile organic solutions. This is in stark contrast to results of common surfactants where the concentration in the liquid plays the major role. We further confirmed our results with a modified adsorption isotherm and molecular dynamics simulations, where highly structured, hydrogen bonded networks, and in particular a solute depletion layer just beneath the Gibbs dividing surface, were revealed.

Comparative all-atomic study of unfolding pathways for proteins chymotrypsin inhibitor 2 and barnase.

The features of transition states and intermediates are important in the study on protein folding. However, transition states and intermediates could not be obviously identified from trajectories obtained by dynamic simulations. In this work, a different method to identify and characterize the transition states and intermediates by combining the root mean square deviation of C(alpha) atoms and the similarity factor Q to the native state is proposed. The unfolding processes based on all-atomic simulations for proteins chymotrypsin inhibitor 2 and barnase are studied, and the related transition states and intermediates are identified by observing an unfolding factor U = 1-F. Comparisons between the conformational cluster analysis and experimental results are also made. The various analyses on the unfolding behaviors indicate that our method can well define the transition states and intermediates, and the factor U (or F) can be used as a reaction coordinate of the folding and unfolding process. It is also found that three-state folding proteins might experience more complicated pathways and have more rugged energy landscapes than two-state folding proteins.