Browse
Submit Data
Databases
API
Help

Dataset Information

21 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Using molecular dynamics simulations to prioritize and understand AI-generated cell penetrating peptides.

ABSTRACT: Cell-penetrating peptides have important therapeutic applications in drug delivery, but the variety of known cell-penetrating peptides is still limited. With a promise to accelerate peptide development, artificial intelligence (AI) techniques including deep generative models are currently in spotlight. Scientists, however, are often overwhelmed by an excessive number of unannotated sequences generated by AI and find it difficult to obtain insights to prioritize them for experimental validation. To avoid this pitfall, we leverage molecular dynamics (MD) simulations to obtain mechanistic information to prioritize and understand AI-generated peptides. A mechanistic score of permeability is computed from five steered MD simulations starting from different initial structures predicted by homology modelling. To compensate for variability of predicted structures, the score is computed with sample variance penalization so that a peptide with consistent behaviour is highly evaluated. Our computational pipeline involving deep learning, homology modelling, MD simulations and synthesizability assessment generated 24 novel peptide sequences. The top-scoring peptide showed a consistent pattern of conformational change in all simulations regardless of initial structures. As a result of wet-lab-experiments, our peptide showed better permeability and weaker toxicity in comparison to a clinically used peptide, TAT. Our result demonstrates how MD simulations can support de novo peptide design by providing mechanistic information supplementing statistical inference.

SUBMITTER: Tran DP

PROVIDER: S-EPMC8137933 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Characterizing the conformational landscape of MDM2-binding p53 peptides using Molecular Dynamics simulations.

Project description:The conformational landscapes of p53 peptide variants and phage derived peptide (12/1) variants, all known to bind to MDM2, are studied using hamiltonian replica exchange molecular dynamics simulations. Complementing earlier observations, the current study suggests that the p53 peptides largely follow the 'conformational selection' paradigm in their recognition of and complexation by MDM2 while the 12/1 peptides likely undergo some element of conformational selection but are mostly driven by 'binding induced folding'. This hypothesis is further supported by pulling simulations that pull the peptides away from their bound states with MDM2. This data extends the earlier mechanisms proposed to rationalize the entropically driven binding of the p53 set and the enthalpically driven binding of the 12/1 set. Using our hypothesis, we suggest mutations to the 12/1 peptide that increase its helicity in simulations and may, in turn, shift the binding towards conformational selection. In summary, understanding the conformational landscapes of the MDM2-binding peptides may suggest new peptide designs with bespoke binding mechanisms.

| S-EPMC5688104 | biostudies-literature

Molecular dynamics simulations of Ibuprofen binding to Abeta peptides.

Project description:Using replica exchange molecular dynamics simulations and the implicit solvent model we probed binding of ibuprofen to Abeta(10-40) monomers and amyloid fibrils. We found that the concave (CV) fibril edge has significantly higher binding affinity for ibuprofen than the convex edge. Furthermore, binding of ibuprofen to Abeta monomers, as compared to fibrils, results in a smaller free energy gain. The difference in binding free energies is likely to be related to the presence of the groove on the CV fibril edge, in which ibuprofen tends to accumulate. The confinement effect of the groove promotes the formation of large low-energy ibuprofen clusters, which rarely occur on the surface of Abeta monomers. These observations led us to suggest that the ibuprofen binding mechanism for Abeta fibrils is different from that for monomers. In general, ibuprofen shows a preference to bind to those regions of Abeta monomers (amino terminal) and fibrils (the CV edge) that are also the primary aggregation interfaces. Based on our findings and on available experimental data, we propose a rationale for the ibuprofen antiaggregation effect.

| S-EPMC2756381 | biostudies-literature

Deep Boosted Molecular Dynamics: Accelerating Molecular Simulations with Gaussian Boost Potentials Generated Using Probabilistic Bayesian Deep Neural Network.

Project description:We have developed a new deep boosted molecular dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30 ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1 μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300 ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on a deep learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

| S-EPMC10355842 | biostudies-literature

Deep Boosted Molecular Dynamics (DBMD): Accelerating molecular simulations with Gaussian boost potentials generated using probabilistic Bayesian deep neural network.

Project description:We have developed a new Deep Boosted Molecular Dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on Deep Learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

| S-EPMC10081221 | biostudies-literature

Structural Analysis and Design of Chionodracine-Derived Peptides Using Circular Dichroism and Molecular Dynamics Simulations.

Project description:Antimicrobial peptides have been identified as one of the alternatives to the extensive use of common antibiotics as they show a broad spectrum of activity against human pathogens. Among these is Chionodracine (Cnd), a host-defense peptide isolated from the Antarctic icefish Chionodraco hamatus, which belongs to the family of Piscidins. Previously, we demonstrated that Cnd and its analogs display high antimicrobial activity against ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa and Enterobacter species). Herein, we investigate the interactions with lipid membranes of Cnd and two analogs, Cnd-m3 and Cnd-m3a, showing enhanced potency. Using a combination of Circular Dichroism, fluorescence spectroscopy, and all-atom Molecular Dynamics (MD) simulations, we determined the structural basis for the different activity among these peptides. We show that all peptides are predominantly unstructured in water and fold, preferentially as ?-helices, in the presence of lipid vesicles of various compositions. Through a series of MD simulations of 400 ns time scale, we show the effect of mutations on the structure and lipid interactions of Cnd and its analogs. By explaining the structural basis for the activity of these analogs, our findings provide structural templates to design minimalistic peptides for therapeutics.

| S-EPMC7073106 | biostudies-literature

Molecular dynamics simulations of spontaneous fibril formation by random-coil peptides.

Project description:Assembly of normally soluble proteins into amyloid fibrils is a cause or associated symptom of numerous human disorders, including Alzheimer's and the prion diseases. We report molecular-level simulation of spontaneous fibril formation. Systems containing 12-96 model polyalanine peptides form fibrils at temperatures greater than a critical temperature that decreases with peptide concentration and exceeds the peptide's folding temperature, consistent with experimental findings. Formation of small amorphous aggregates precedes ordered nucleus formation and subsequent rapid fibril growth through addition of beta-sheets laterally and monomeric peptides at fibril ends. The fibril's structure is similar to that observed experimentally.

| S-EPMC526199 | biostudies-literature

Molecular dynamics simulations of peptides and proteins with amplified collective motions.

Project description:We present a novel method that uses the collective modes obtained with a coarse-grained model/anisotropic network model to guide the atomic-level simulations. Based on this model, local collective modes can be calculated according to a single configuration in the conformational space of the protein. In the molecular dynamics simulations, the motions along the slowest few modes are coupled to a higher temperature by the weak coupling method to amplify the collective motions. This amplified-collective-motion (ACM) method is applied to two test systems. One is an S-peptide analog. We realized the refolding of the denatured peptide in eight simulations out of 10 using the method. The other system is bacteriophage T4 lysozyme. Much more extensive domain motions between the N-terminal and C-terminal domain of T4 lysozyme are observed in the ACM simulation compared to a conventional simulation. The ACM method allows for extensive sampling in conformational space while still restricting the sampled configurations within low energy areas. The method can be applied in both explicit and implicit solvent simulations, and may be further applied to important biological problems, such as long timescale functional motions, protein folding/unfolding, and structure prediction.

| S-EPMC1302944 | biostudies-literature

Molecular dynamics simulations to understand glycosaminoglycan interactions in the free- and protein-bound states.

Project description:Natural glycosaminoglycans (GAGs) are informational molecules with astounding structural diversity. Understanding the behavior of GAGs in the free and protein-bound states is critical for harnessing this diversity. Molecular dynamics (MD) offers atomistic insight into principles governing GAG recognition by proteins. Here, we discuss how MD can be used to understand local and global properties of GAGs in free solution, including torsions, puckering, hydrogen bonding, flexibility, and energetics. We discuss MD studies on GAG-protein complexes, which help elucidate the strength of interacting residues, role of water, energetics, and so on. The MD results accumulated so far suggest that GAG recognition of proteins is a continuum from the highly selective on one end to the fully non-selective on the other with intermediate levels of selectivity, including moderately selective and plastic. The advancements in MD technology, such as coarse-grained MD, coupled with really long simulations will help understand macroscale molecular movements in the future.

| S-EPMC9189024 | biostudies-literature

Molecular dynamics simulations of glycoproteins using CHARMM.

Project description:Molecular dynamics simulations are an effective tool to study the structure, dynamics, and thermodynamics of carbohydrates and proteins. However, the simulations of heterogeneous glycoprotein systems have been limited due to the lack of appropriate molecular force field parameters describing the linkage between the carbohydrate and the protein regions as well as the tools to prepare these systems for modeling studies. In this work we outline the recent developments in the CHARMM carbohydrate force field to treat glycoproteins and describe in detail the step-by-step procedures involved in building glycoprotein geometries using CHARMM-GUI Glycan Reader.

| S-EPMC4537648 | biostudies-literature

Fighting COVID-19 Using Molecular Dynamics Simulations.

Project description: Not available

| S-EPMC7571292 | biostudies-literature