Improved Parameters for the Martini Coarse-Grained Protein Force Field.

The formation of the Cluster of Differentiation 47 (CD47, PDB code: 2JJT)/signal regulatory protein α (SIRPα) complex is very important as it protects healthy cells from immune clearance while promoting macrophage phagocytosis for tumour elimination. Although several antibodies have been developed for cancer therapy, new function-blocking aptamers are still under development. This study aims to design the aptamer AptCD47, which can block the formation of the CD47/SIRPα complex. This study employs the MARTINI coarse-grained (CG) force field and the stochastic tunnelling-basin hopping-discrete molecular dynamics (STUN-BH-DMD) method to identify the most stable AptCD47/CD47 complexes. Coarse-grained molecular dynamics (CGMD) simulations were used to obtain root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analyses. The results demonstrate that the formation of AptCD47/CD47 complexes renders the CD47 structure more stable than the single CD47 molecule in a water environment. The minimum energy pathway (MEP) obtained by the nudged elastic band (NEB) method indicates that the binding processes of 5′-ATTCAATTCC-3′ and 5′-AGTGCAATCT-3′ to CD47 are barrierless, which is much lower than the binding barrier of SIRPα to CD47 of about 14.23 kcal/mol. Therefore, these two AptCD47/CD47 complexes can create a high spatial binding barrier for SIRPα, preventing the formation of a stable CD47/SIRPα complex. The proposed numerical process with the MARTINI CG force field can be used to design CD47 aptamers that efficiently block SIRPα from binding to CD47. Communicated by Ramaswamy H. Sarma

In this article we study the lateral molecular diffusion in homogeneous lipid 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers by molecular dynamics simulations with the all-atom Optimised Potentials for Liquid Simulations and the coarse-grained MARTINI force field. On the statistically relevant time scales, the centre-of-mass mean-square displacement exhibits in both cases the subdiffusive asymptotic form , with and , respectively. The diffusive dynamics obtained by the MARTINI force field is, however, faster by a factor of about 3. The subdiffusive characteristics of the diffusion process is confirmed by comparing the integral of the centre-of-mass velocity autocorrelation function with its analytical long-time tail. The agreement is particularly good for the MARTINI force field, which permits to extend the simulation length and the system size considerably. Our results are in agreement with experimental observations of subdiffusion in lipid bilayers on longer time scales and do not support the finding of some authors that the latter should be considered as a transient phenomenon.

A Structurally Flexible Protein Backbone for the MARTINI Coarse Grained Force Field

We simulate poly(ethylene glycol) (PEG) oligomers and model Tween 80 (polyoxyethylene sorbitan monooleate) molecules at water/alkane interfaces. Using the weighted histogram analysis method (WHAM), including an extension of WHAM to two reaction coordinates to remove hysteresis, we calculate interfacial potentials of mean force (PMFs) for PEG and Tween 80 using three force fields: the atomistic GROMOS 53a6OXY+D and two coarse-grained (CG) MARTINI force fields. Because the force fields have not yet been validated for PEO adsorption to hydrophobic interfaces, we calculate PMFs for alcohol ethoxylates C12E2 and C12E8 and find that they agree with semiempirical results of Mulqueen and Blankschtein [Langmuir 2002, 18 (2), 365-376] for the GROMOS 53a6OXY+D force field, whereas for both MARTINI force fields, PEO adsorbs too weakly to a clean hydrophobic interface. One MARTINI force field incorrectly shows depletion rather than adsorption to a clean hydrophobic interface. We find that the adsorption free energy for PEG oligomers at a clean, planar water/alkane interface is around 1.3 kBT per monomer for the atomistic force field but is less than half of this for the two CG force fields. With the newly validated GROMOS 53a6OXY+D force field, we bracket the dilute adsorption free energy for a model Tween 80 molecule at the clean water/squalane interface. We also calculate the pressure-area isotherm. We exploit these data with the Nikas-Mulqueen-Blankschtein (NMB) theory and a simple transport model to demonstrate a transition from irreversible to reversible adsorption with increasing surface coverage, consistent with experimental results of Reichert and Walker [Langmuir 2013, 29 (6), 1857-1867].

Coarse-grained (CG) simulation methods are now widely used to model the structure and dynamics of large biomolecular systems. One important issue for using such methods - especially with regard to using them to model, for example, intracellular environments - is to demonstrate that they can reproduce experimental data on the thermodynamics of protein-protein interactions in aqueous solutions. To examine this issue, we describe here simulations performed using the popular coarse-grained MARTINI force field, aimed at computing the thermodynamics of lysozyme and chymotrypsinogen self-interactions in aqueous solution. Using molecular dynamics simulations to compute potentials of mean force between a pair of protein molecules, we show that the original parameterization of the MARTINI force field is likely to significantly overestimate the strength of protein-protein interactions to the extent that the computed osmotic second virial coefficients are orders of magnitude more negative than experimental estimates. We then show that a simple down-scaling of the van der Waals parameters that describe the interactions between protein pseudo-atoms can bring the simulated thermodynamics into much closer agreement with experiment. Overall, the work shows that it is feasible to test explicit-solvent CG force fields directly against thermodynamic data for proteins in aqueous solutions, and highlights the potential usefulness of osmotic second virial coefficient measurements for fully parameterizing such force fields.

The MARTINI coarse-grained force field [Monticelli, L. et al. J. Chem. Theory Comput.2008, 4, 819-834] is examined for use in molecular dynamics simulations of the dissociation of gas-phase protein complexes. Coarse-grained force fields allow longer time scales and larger systems to be treated compared with all-atom force fields. In this work, results for the dissociation of the cytochrome c' dimer using MARTINI are compared with published studies using the OPLS-AA/L all-atom force field. Several structural parameters such as the minimum distance between monomers, radius of gyration, and root-mean-square deviation as well as potential energy contributions (Lennard-Jones and Coulomb) are calculated as a function of the center of mass distance. The MARTINI force field semiquantitatively reproduces the results of previous all-atom studies but appears to be somewhat too attractive.

Covering All the Bases: A Martini Coarse-Grained Force Field for DNA

Hybrid molecular dynamics simulations of atomistic (AA) solutes embedded in coarse-grained (CG) environment can substantially reduce the computational cost with respect to fully atomistic simulations. However, interfacing both levels of resolution is a major challenge that includes a balanced description of the relevant interactions. This is especially the case for polar solvents such as water, which screen the electrostatic interactions and thus require explicit electrostatic coupling between AA and CG subsystems. Here, we present and critically test computationally efficient hybrid AA/CG models. We combined the Gromos atomistic force field with the MARTINI coarse-grained force field. To enact electrostatic coupling, two recently developed CG water models with explicit electrostatic interactions were used: the polarizable MARTINI water model and the BMW model. The hybrid model was found to be sensitive to the strength of the AA-CG electrostatic coupling, which was adjusted through the relative dielectric permittivity εr(AA-CG). Potentials of mean force (PMFs) between pairs of amino acid side chain analogues in water and partitioning free enthalpies of uncharged amino acid side chain analogues between apolar solvent and water show significant differences between the hybrid simulations and the fully AA or CG simulations, in particular for charged and polar molecules. For apolar molecules, the results obtained with the hybrid AA/CG models are in better agreement with the fully atomistic results. The structures of atomistic ubiquitin solvated in CG water and of a single atomistic transmembrane α-helix and the transmembrane portion of an atomistic mechanosensitive channel in CG lipid bilayers were largely maintained during 50-100 ns of AA/CG simulations, partly due to an overstabilization of intramolecular interactions. This work highlights some key challenges on the way toward hybrid AA/CG models that are both computationally efficient and sufficiently accurate for biomolecular simulations.

The investigation of the compounding mechanisms of mixed collectors in low-rank coal flotation has gradually progressed from a macro to a micro level, and the accuracy of research outcomes significantly depends on the development of micro models. While All-Atom Molecular Dynamics (AAMD) simulations have been used for this purpose, the employment of an excessive amount of collector in the system leads to a failure in replicating the actual flotation process, which in turn, limits the practical applicability of research findings. To overcome this drawback, this paper suggests using the MARTINI Coarse-Grained force field (CGFF), developed by the S.J.Marrink group, to construct a simulation system that accurately mimics the amount of collector employed in practical flotation processes. In light of this, to make the MARTINI force field more applicable for research in the low-rank coal field, we employ density functional theory (DFT) and coarse-grained molecular dynamics (CGMD) methods to investigate water cluster interactions, identify magic number clusters, and fit the corresponding CGFF parameters to develop a CG water model suitable for exploring the compounding mechanism of mixed collectors. This research lays the groundwork for future investigations into the aggregation behavior of mixed collectors under actual usage conditions, as well as their synergistic adsorption at the coal/water interface.

The conversion of coarse-grained to atomistic models is an important step in obtaining insight about atomistic scale processes from coarse-grained simulations. For this process, called backmapping or reverse transformation, several tools are available, but these commonly require libraries of molecule fragments or they are linked to a specific software package. In addition, the methods are usually restricted to specific molecules and to a specific force field. Here, we present an alternative method, consisting of geometric projection and subsequent force-field based relaxation. This method is designed to be simple and flexible, and offers a generic solution for resolution transformation. For simple systems, the conversion only requires a list of particle correspondences on the two levels of resolution. For special cases, such as nondefault protonation states of amino acids and virtual sites, a target particle list can be specified. The mapping uses simple building blocks, which list the particles on the different levels of resolution. For conversion to higher resolution, the initial model is relaxed with several short cycles of energy minimization and position-restrained MD. The reconstruction of an atomistic backbone from a coarse-grained model is done using a new dedicated algorithm. The method is generic and can be used to map between any two particle based representations, provided that a mapping can be written. The focus of this work is on the coarse-grained MARTINI force field, for which mapping definitions are written to allow conversion to and from the higher-resolution force fields GROMOS, CHARMM, and AMBER, and to and from a simplified three-bead lipid model. Together, these offer the possibility to simulate mesoscopic membrane structures, to be transformed to MARTINI and subsequently to an atomistic model for investigation of detailed interactions. The method was tested on a set of systems ranging from a simple, single-component bilayer to a large protein-membrane-solvent complex. The results demonstrate the efficiency and the efficacy of the new approach.

The coarse-grained MARTINI force field, initially developed for membranes, has proven to be an exceptional tool for investigating supramolecular peptide assemblies. Over the years, the force field underwent refinements to enhance accuracy, enabling, for example, the reproduction of protein-ligand interactions and constant pH behavior. However, these protein-focused improvements seem to have compromised its ability to model short peptide self-assembly. In this study, we assess the performance of MARTINI 3 in reproducing peptide self-assembly using the well-established diphenylalanine (FF) as our test case. Unlike its success in version 2.1, FF does not even exhibit aggregation in version 3. By systematically exploring parameters for the aromatic side chains and charged backbone beads, we established a parameter set that effectively reproduces tube formation. Remarkably, these parameter adjustments also replicate the self-assembly of other di- and tripeptides and coassemblies. Furthermore, our analysis uncovers pivotal insights for enhancing the performance of MARTINI in modeling short peptide self-assembly. Specifically, we identify issues stemming from overestimated hydrophilicity arising from charged termini and disruptions in π-stacking interactions due to insufficient planarity in aromatic groups and a discrepancy in intermolecular distances between this and backbone-backbone interactions. This investigation demonstrates that strategic modifications can harness the advancements offered by MARTINI 3 for the realm of short peptide self-assembly.

The interactions between amino acid side chains govern protein secondary, tertiary, and quaternary structure formation. For molecular modeling approaches to be able to realistically describe these phenomena, the underlying force fields have to represent these interactions as accurately as possible. Here, we compare the side chain-side chain interactions for a number of commonly used force fields, namely the all-atom OPLS, the united-atom GROMOS, and the coarse-grain MARTINI force field. We do so by calculating the dimerization free energies between selected pairs of side chains and structural characterization of their binding modes. To mimic both polar and nonpolar environments, the simulations are performed in water, n-octanol, and decane. In general, reasonable correlations are found between all three force fields, with deviations on the order of 1 kT in aqueous solvent. In apolar solvent, however, significantly larger differences are found, especially for charged amino acid pairs between the OPLS and GROMOS force fields, and for polar interactions in the MARTINI force field in comparison to the higher resolution models. Interestingly, even in cases where the dimerization free energies are similar, the binding mode may differ substantially between the force fields. This was found to be especially the case for aromatic residues. In addition to the inter-force-field comparison, we compared the various force fields to a knowledge-based potential. The two independent approaches show good correlation in aqueous solvent with an exception of aromatic residues for which the interaction strength is lower in the knowledge-based potentials.

Aggregation kinetics of short peptides: All-atom and coarse-grained molecular dynamics study

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models.

The MARTINI coarse-grained (CG) force field is used to test the ability of CG models to simulate ionic transport through protein nanopores. The ionic conductivity of CG ions in solution was computed and compared with experimental results. Next, we studied the electrostatic behavior of a solvated CG lipid bilayer in salt solution under an external electric field. We showed this approach correctly describes the experimental conditions under a potential bias. Finally, we performed CG molecular dynamics simulations of the ionic transport through a protein nanopore (α-hemolysin) inserted in a lipid bilayer, under different electric fields, for 2–3 microseconds. The resulting I − V curve is qualitatively consistent with experiments, although the computed current is one order of magnitude smaller. Current saturation was observed for potential biases over ±350 mV. We also discuss the time to reach a stationary regime and the role of the protein flexibility in our CG simulations.