Mesoscopic Simulation of Phospholipid Membranes, Peptides, and Proteins with Molecular Fragment Dynamics

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

Advancements have been made to dissipative particle dynamics (DPD), a robust coarse-grained (CG) simulation method, to study the folded structures of four miniproteins (1L2Y, 1WN8, 1YRF, and 2I9M) in explicit and implicit solvents. In this endeavor, we aim to establish model parametrization and enhance computational efficiency. Unlike traditional CG models that use empirical force parameters, ex-force parameters (r0(ex), a ~ , δd, δp) of DPD particles constructed for specific research purposes can be obtained from atomistic molecular dynamics simulations. On the other hand, im-force parameters (r0(im), c, σ) can be derived from ex-DPD simulations, according to the underlying thermodynamic theory. Based on a mapping scheme proposed for the modeling of amino acids, all-atom proteins can be converted into a CG model. Both ex-/im-DPDs are then carried out to investigate the folding pathways of the four mini-proteins. Structural analysis of the RMSDs shows that the im-simulated proteins have greater structural similarity to native proteins than the ex-simulated ones. The constructed CG models achieve a resolution of Angstrom (Å), a level normally associated with atomic models. Additionally, speed tests reveal that im-DPD accelerates the simulation process and significantly improves simulation efficiency.

A Structure Independent Molecular Fragment Interfuse Model for Mesoscale Dissipative Particle Dynamics Simulation of Peptides

There is a need to develop robust computational models for mesoscale simulation of the structure of peptides over large length scales toward the discovery of novel peptides for medical applications to address the issues of peptide aggregation, enzymatic degradation, and short half-life. The primary objective was to predict the structure and conformation of peptides whose native structures are not known. This work presents a new model for computation of interaction parameters between the beads in coarse-grained dissipative particle dynamics (DPD) simulation that is properly calibrated for amino acids, supports compressibility requirement of water molecules, and accounts for subtle differences in the structure of amino acids and the charge in the side chain of charged amino acids. This new model is referred to as Structure Independent Molecular Fragment Interfuse Model, abbreviated as SIMFIM, because it accounts for specific interactions between different beads, which represent molecular fragments of the amino acids, in calculating nonbonded interaction parameters in the absence of knowing the actual peptide structure. The electrostatic interactions are incorporated in this model by using a normal distribution of charges around the center of the beads to prevent the collapse of oppositely charged soft beads. The uniquely parameterized DPD force field in the SIMFIM model is optimized for a given peptide with respect to the degree of coarse-grained graining for simulating the peptide over long times and length scales. The SIMFIM model was tested in this work using four peptides, namely, TrpZip2, Rubrivinodin, Lihuanodin, and IC3-CB1/Gai peptides, whose structures were sourced from the Protein Data Bank. The SIMFIM model predicted radius of gyration (Rg) values for the peptides closer to the actual structures as compared to the conventional model, and there was less deviation between the predicted and actual structures of the peptides.

Notes on the Treatment of Charged Particles for Studying Cyclotide/Membrane Interactions with Dissipative Particle Dynamics

Different charge treatment approaches are examined for cyclotide-induced plasma membrane disruption by lipid extraction studied with dissipative particle dynamics. A pure Coulomb approach with truncated forces tuned to avoid individual strong ion pairing still reveals hidden statistical pairing effects that may lead to artificial membrane stabilization or distortion of cyclotide activity depending on the cyclotide's charge state. While qualitative behavior is not affected in an apparent manner, more sensitive quantitative evaluations can be systematically biased. The findings suggest a charge smearing of point charges by an adequate charge distribution. For large mesoscopic simulation boxes, approximations for the Ewald sum to account for mirror charges due to periodic boundary conditions are of negligible influence.

A generic force field for simulating native protein structures using dissipative particle dynamics

A coarse-grained force field for molecular dynamics simulations of the native structures of proteins in a dissipative particle dynamics (DPD) framework is developed. The parameters for bonded interactions are derived by mapping the bonds and angles for 20 amino acids onto target distributions obtained from fully atomistic simulations in explicit solvent. A dual-basin potential is introduced for stabilizing backbone angles, to cover a wide spectrum of protein secondary structures. The backbone dihedral potential enables folding of the protein from an unfolded initial state to the folded native structure. The proposed force field is validated by evaluating the structural properties of several model peptides and proteins including the SARS-CoV-2 fusion peptide, consisting of α-helices, β-sheets, loops and turns. Detailed comparisons with fully atomistic simulations are carried out to assess the ability of the proposed force field to stabilize the different secondary structures present in proteins. The compact conformations of the native states were evident from the radius of gyration and the high intensity peaks of the root mean square deviation histograms, which were found to be within 0.4 nm. The Ramachandran-like energy landscape on the phase space of backbone angles (θ) and dihedrals (ϕ) effectively captured the conformational phase space of α-helices at ∼(ϕ = 50°,θ = 90°) and β-strands at ∼(ϕ = ±180°,θ = 90-120°). Furthermore, the residue-residue native contacts were also well reproduced by the proposed DPD model. The applicability of the model to multidomain complexes was assessed using lysozyme and a large α-helical bacterial pore-forming toxin, cytolysin A. Our study illustrates that the proposed force field is generic, and can potentially be extended for efficient in silico investigations of membrane bound polypeptides and proteins using DPD simulations.

Quantitative Estimation of Cyclotide-Induced Bilayer Membrane Disruption by Lipid Extraction with Mesoscopic Simulation

Cyclotide-induced membrane disruption is studied at the microsecond timescale by dissipative particle dynamics to quantitatively estimate a kinetic rate constant for membrane lipid extraction with a ″sandwich″ interaction model where two bilayer membranes enclose a cyclotide/water compartment. The obtained bioactivity trends for cyclotides Kalata B1, Cycloviolacin O2, and selected mutants with different membrane types are in agreement with experimental findings: For all membranes investigated, Cycloviolacin O2 shows a higher lipid extraction activity than Kalata B1. The presence of cholesterol leads to a decreased cyclotide activity compared to cholesterol-free membranes. Phosphoethanolamine-rich membranes exhibit an increased membrane disruption. A cyclotide's ″hydrophobic patch″ surface area is important for its bioactivity. A replacement of or with charged amino acid residues may lead to super-mutants with above-native activity but without simple charge-activity patterns. Cyclotide mixtures show linearly additive bioactivities without significant sub- or over-additive effects. The proposed method can be applied as a fast and easy-to-use tool for exploring structure-activity relationships of cyclotide/membrane systems: With the open software provided, the rate constant of a single cyclotide/membrane system can be determined in about 1 day by a scientific end-user without programming skills.

MFsim-an open Java all-in-one rich-client simulation environment for mesoscopic simulation

MFsim is an open Java all-in-one rich-client computing environment for mesoscopic simulation with Jdpd as its default simulation kernel for Molecular Fragment (Dissipative Particle) Dynamics. The new environment comprises the complete preparation-simulation-evaluation triad of a mesoscopic simulation task and especially enables biomolecular simulation tasks with peptides and proteins. Productive highlights are a SPICES molecular structure editor, a PDB-to-SPICES parser for particle-based peptide/protein representations, a support of polymer definitions, a compartment editor for complex simulation box start configurations, interactive and flexible simulation box views including analytics, simulation movie generation or animated diagrams. As an open project, MFsim allows for customized extensions for different fields of research.

SPICES: a particle-based molecular structure line notation and support library for mesoscopic simulation

Simplified Particle Input ConnEction Specification (SPICES) is a particle-based molecular structure representation derived from straightforward simplifications of the atom-based SMILES line notation. It aims at supporting tedious and error-prone molecular structure definitions for particle-based mesoscopic simulation techniques like Dissipative Particle Dynamics by allowing for an interplay of different molecular encoding levels that range from topological line notations and corresponding particle-graph visualizations to 3D structures with support of their spatial mapping into a simulation box. An open Java library for SPICES structure handling and mesoscopic simulation support in combination with an open Java Graphical User Interface viewer application for visual topological inspection of SPICES definitions are provided.

Jdpd: an open java simulation kernel for molecular fragment dissipative particle dynamics

Jdpd is an open Java simulation kernel for Molecular Fragment Dissipative Particle Dynamics with parallelizable force calculation, efficient caching options and fast property calculations. It is characterized by an interface and factory-pattern driven design for simple code changes and may help to avoid problems of polyglot programming. Detailed input/output communication, parallelization and process control as well as internal logging capabilities for debugging purposes are supported. The new kernel may be utilized in different simulation environments ranging from flexible scripting solutions up to fully integrated “all-in-one” simulation systems.

Multiscale simulations for understanding the evolution and mechanism of hierarchical peptide self-assembly

Multiscale molecular simulations that combine and systematically link several hierarchies can provide insights into the evolution and dynamics of hierarchical peptide self-assembly from the molecular level to the mesoscale.