Search

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. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)

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., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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.

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.