Browse
Submit Data
Databases
API
Help

Dataset Information

20 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

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 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 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

Confronting pitfalls of AI-augmented molecular dynamics using statistical physics.

Project description:Artificial intelligence (AI)-based approaches have had indubitable impact across the sciences through the ability to extract relevant information from raw data. Recently, AI has also found use in enhancing the efficiency of molecular simulations, wherein AI derived slow modes are used to accelerate the simulation in targeted ways. However, while typical fields where AI is used are characterized by a plethora of data, molecular simulations, per construction, suffer from limited sampling and thus limited data. As such, the use of AI in molecular simulations can suffer from a dangerous situation where the AI-optimization could get stuck in spurious regimes, leading to incorrect characterization of the reaction coordinate (RC) for the problem at hand. When such an incorrect RC is then used to perform additional simulations, one could start to deviate progressively from the ground truth. To deal with this problem of spurious AI-solutions, here, we report a novel and automated algorithm using ideas from statistical mechanics. It is based on the notion that a more reliable AI-solution will be one that maximizes the timescale separation between slow and fast processes. To learn this timescale separation even from limited data, we use a maximum caliber-based framework. We show the applicability of this automatic protocol for three classic benchmark problems, namely, the conformational dynamics of a model peptide, ligand-unbinding from a protein, and folding/unfolding energy landscape of the C-terminal domain of protein G. We believe that our work will lead to increased and robust use of trustworthy AI in molecular simulations of complex systems.

| S-EPMC7863682 | 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

Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange.

Project description:Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.

| S-EPMC1877756 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data