Dimerization Free Energy Research Articles

Thermodynamics of the gas-phase dimerization of formic acid: Fully anharmonic finite temperature calculations at the CCSD(T) and many DFT levels.

A proof-of-concept study is undertaken to demonstrate the utility of the machine learning combined with the thermodynamic perturbation theory (MLPT) to test the accuracy of electronic structure methods in finite-temperature thermodynamic calculations. As a test example, formic acid dimer is chosen, which is one of the systems included in the popular benchmark set S22 [Jurečka et al., Phys. Chem. Chem. Phys. 8, 1985-1993 (2006)]. Starting from the explicit molecular dynamics and thermodynamic integration performed at the PBE + D2 level, the MLPT is used to obtain fully anharmonic dimerization free and internal energies at the reference quality CCSD(T) level and 19 different density functional approximations, including GGA, meta-GGA, non-local, and hybrid functionals with and without dispersion corrections. Our finite-temperature results are shown to be both qualitatively and quantitatively different from those obtained using the conventional benchmarking strategy based on fixed structures. The hybrid functional HSE06 is identified as the best performing approximate method tested, with the errors in free and internal energies of dimerization being 36 and 41meV, respectively.

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes.

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Improved Protein Model in SPICA Force Field.

The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.

A thermodynamic analysis of CLC transporter dimerization in lipid bilayers

The CLC-ec1 chloride/proton antiporter is a membrane-embedded homodimer with subunits that can dissociate and associate, but the thermodynamic driving forces favor the assembled dimer at biological densities. Yet, the physical reasons for this stability are confounding as dimerization occurs via the burial of hydrophobic interfaces away from the lipid solvent. For binding of nonpolar surfaces in aqueous solution, the driving force is often attributed to the hydrophobic effect, but this should not apply in the membrane since there is very little water. To investigate this further, we quantified the thermodynamic changes associated with CLC dimerization in membranes by carrying out a van 't Hoff analysis of the temperature dependency of the free energy of dimerization, ΔG°. To ensure that the reaction reached equilibrium at different temperatures, we utilized a Förster resonance energy transfer assay to report on relaxation kinetics of subunit exchange as a function of temperature. Equilibration times were then applied to measure CLC-ec1 dimerization isotherms at different temperatures using the single-molecule subunit-capture photobleaching analysis approach. The results demonstrate that the dimerization free energy of CLC in Escherichia coli-like membranes exhibits a nonlinear temperature dependency corresponding to a large, negative change in heat capacity, a signature of solvent ordering effects such as the hydrophobic effect. Consolidating this with our previous molecular analyses suggests that the nonbilayer defect required to solvate the monomeric state is one source of the observed change in heat capacity and indicates the existence of a generalizable driving force for protein association in membranes.

Decoupling of Interactions between Model-Charged Peptides Reveals Key Factors Responsible for Liquid-Liquid Phase Separation.

Liquid-liquid phase separation (LLPS) by disordered proteins has been shown to govern biological processes and cause numerous diseases. Therefore, a deeper understanding of the interactions and their variation with external factors is key to modulating the LLPS behavior of different systems and protecting proteins from pathological aggregation. In this context, we have looked at interactions between similarly charged peptides to understand the molecular features that may drive or prevent condensate formation under various conditions. We have studied dimer formation for model peptides where charged and noncharged amino acids have been placed alternatively. Using arginine and glutamic acid as the charged residues and varying the other residues with glycine, alanine, and proline to alter hydrophobicity, we have obtained the free-energy surface (FES) for the dimer formation for these systems under high salt concentration at two different temperatures using all-atom molecular dynamics simulations and the well-tempered metadynamics method. Our results indicate that a combination of effects such as hydrophobicity, arginine-arginine interactions, or water release from the solvation shell makes dimerization free energy more favorable for the positively charged peptides with lower flexibility. For the negatively charged peptides, the crucial role of water has been found in governing the FES. Systems having charged residues and phenylalanine in the peptide sequence also have been studied at high salt concentrations using unbiased simulations. In this case, only the positively charged peptides were found to aggregate through temperature-dependent hydrophobic and cation-π interactions. Overall, our study indicates that the negatively charged peptides are more likely to remain in the dilute phase under various conditions compared to the positively charged systems. The findings from our study would be helpful in designing and controlling systems to obtain LLPS behavior for therapeutic usage.

Understanding Molecular Aggregation of Ligand-Protected Atomically-Precise Metal Nanoclusters.

Monolayer-protected atomically precise nanoclusters (MPCs) are an important class of molecules due to their unique structural features and diverse applications, including bioimaging, sensors, and drug carriers. Understanding the atomistic and dynamical details of their self-assembly process is crucial for designing system-specific applications. Here, we applied molecular dynamics and on-the-fly probability-based enhanced sampling simulations to study the aggregation of Au25(pMBA)18 MPCs in aqueous and methanol solutions. The MPCs interact via both hydrogen bonds and π-stacks between the aromatic ligands to form stable dimers, oligomers, and crystals. The dimerization free energy profiles reveal a pivotal role of the ligand charged state and solvent mediating the molecular aggregation. Furthermore, MPCs' ligands exhibit suppressed conformational flexibility in the solid phase due to facile intercluster hydrogen bonds and π-stacks. Our work provides unprecedented molecular-level dynamical details of the aggregation process and conformational dynamics of MPCs ligands in solution and crystalline phases.

Counteractive Effects of Choline Geranate (CAGE) ILs and Ethanol on Insulin’s Stability—A Leap Forward towards Oral Insulin Formulation

Choline geranate (CAGE) ionic liquids (ILs) stabilize insulin, thereby aiding its oral delivery, whereas ethanol (EtOH) affects its stability by disrupting the hydrophobic interactions. In this study, cognizance of the stabilization mechanism of insulin dimer in the presence of both CAGE ILs and EtOH mixtures is achieved through biased and unbiased molecular dynamics (MD) simulations. Here, two order parameters are employed to study the insulin dimer dissociation using well-tempered metadynamics (WT-MetaD). The stability of insulin is found to be strongly maintained until a 0.20 mole fraction of EtOH. Besides, higher concentrations of EtOH marginally affect the insulin stability. Moreover, geranate anions form a higher number of H-bonding interactions with water molecules, which aids insulin stabilization. Conversely, the addition of EtOH minimizes the water-mediated H-bonding interactions of geranate. Additionally, geranate traps the EtOH molecules, thereby preventing the interactions between insulin and EtOH. Furthermore, the free energy landscape (FEL) reveals the absence of dimer dissociation along with noticeable deviations in the distances R and the number of contacts Q. The dimerization free energy of insulin was calculated to be −16.1 kcal/mol at a 0.20 mole fraction of EtOH. Moreover, increments in mole fractions of EtOH effectuate a decrease in the insulin stability. Thus, the present study represents CAGE ILs as efficient insulin dimer stabilizes at low concentrations of EtOH.

On Computing Equilibrium Binding Constants for Protein-Protein Association in Membranes.

Protein association in lipid membranes is fundamental to membrane protein function and of great biomedical relevance. All-atom and coarse-grained models have been extensively used to understand the protein-protein interactions in the membrane and to compute equilibrium association constants. However, slow translational and rotational diffusion of protein in membrane presents challenges to the effective sampling of conformations defining the ensembles of free and bound states contributing to the association equilibrium and the free energy of dimerization. We revisit the homodimerization equilibrium of the TM region of glycophorin A. Conformational sampling is performed using umbrella sampling along previously proposed one-dimensional collective variables and compared with sampling over a two-dimensional collective variable space using the MARTINI v2.2 force field. We demonstrate that the one-dimensional collective variables suffer from restricted sampling of the native homodimer conformations leading to a biased free energy landscape. Conversely, simulations along the two-dimensional collective variable effectively characterize the thermodynamically relevant native and non-native interactions contributing to the association equilibrium. These results demonstrate the challenges associated with accurately characterizing binding equilibria when multiple poses contribute to the bound state ensemble.

SPICA Force Field for Proteins and Peptides.

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin.

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a relaxed state at pH 4-8 (with the typical alignment of the crystallographic structure) and a compact state at pH 3-4 (with a tighter association and rotated alignment). This work shows that an approach based on constant-pH MD simulations can produce rich detailed pictures of pH-dependent protein associations, as illustrated for BLG dimerization.

Van’t Hoff analysis of CLC-ec1 dimerization in Escherichia coli lipid membranes

What are the thermodynamic forces that drive hydrophobic membrane proteins to form stable complexes with other hydrophobic proteins in the greasy lipid solvent? In order to investigate this question, we must have a way of dissecting the changes in thermodynamic properties, i.e. enthalpy and entropy, associated with membrane protein reactions. The CLC-ec1 chloride/proton antiporter forms a homodimer that participates in an association equilibrium with a dimerization free energy of −13 kcal/mole in E.

Structural and energetic features of the dimerization of the main proteinase of SARS-CoV-2 using molecular dynamic simulations.

The COVID-19 pandemic caused by SARS-CoV-2 has been declared a global health crisis. The development of anti-SARS-CoV-2 drugs heavily depends on the systematic study of the critical biological processes of key proteins of coronavirus among which the main proteinase (Mpro) dimerization is a key step for virus maturation. Because inhibiting the Mpro dimerization can efficiently suppress virus maturation, the key residues that mediate dimerization can be treated as targets of drug and antibody developments. In this work, the structure and energy features of the Mpro dimer of SARS-CoV-2 and SARS-CoV were studied using molecular dynamics (MD) simulations. The free energy calculations using the Generalized Born (GB) model showed that the dimerization free energy of the SARS-CoV-2 Mpro dimer (-107.5 ± 10.89 kcal mol-1) is larger than that of the SARS-CoV Mpro dimer (-92.83 ± 9.81 kcal mol-1), indicating a more stable and possibly a quicker formation of the Mpro dimer of SARS-CoV-2. In addition, the energy decomposition of each residue revealed 11 key attractive residues. Furthermore, Thr285Ala weakens the steric hindrance between the two protomers of SARS-CoV-2 that can form more intimate interactions. It is interesting to find 11 repulsive residues which effectively inhibit the dimerization process. At the interface of the Mpro dimer, we detected three regions that are rich in interfacial water which stabilize the SARS-CoV-2 Mpro dimer by forming hydrogen bonds with two protomers. The key residues and rich water regions provide important targets for the future design of anti-SARS-CoV-2 drugs through inhibiting Mpro dimerization.

Designing a Solution-Stable Distannene: The Decisive Role of London Dispersion Effects in the Structure and Properties of {Sn(C6H2-2,4,6-Cy3)2}2 (Cy = Cyclohexyl).

The reaction of 1 equiv of the dimeric lithium salt of a new London dispersion effect donor ligand {Li(C6H2-2,4,6-Cy3)·OEt2}2 (Cy = cyclohexyl) with SnCl2 afforded the distannene {Sn(C6H2-2,4,6-Cy3)2}2 (1). The distannene remains dimeric in solution, as indicated by its room-temperature 119Sn NMR signal (δ = 361.3 ppm) and its electronic spectrum, which is invariant over the temperature range of -10 to 100 °C. The formation of the distannene, which has a short Sn-Sn distance of 2.7005(7) Å and greatly enhanced stability in solution compared to that of other distannenes, is due to increased interligand London dispersion (LD) attraction arising from multiple close approaches of ligand C-H moieties across the Sn-Sn bond. DFT-D4 calculations revealed a dispersion stabilization of dimer 1 of 38 kcal mol-1 and a dimerization free energy of ΔGdimer = -6 kcal mol-1. In contrast, the reaction of 2 equiv of the similarly shaped but less bulky, less hydrogen-rich Li(C6H2-2,4,6-Ph3)·(OEt2)2 with SnCl2 yielded the monomeric stannylene Sn(C6H2-2,4,6-Ph3)2 (2), which is unstable in solution at ambient temperature.

Stereo‐dependent dimerization, boiling points, diffusion coefficients, and dielectric constants of E/Z‐HFO‐1234ze

AbstractStereo‐dependent properties of E‐ and Z‐HFO‐1234ze were investigated by density functional theory, extrapolated complete basis set limit, explicitly correlation, and symmetry‐adapted perturbation theory together with molecular dynamics simulations. The enthalpy of isomerization of E‐HFO‐1234ze to the Z isomer is 9.7 kJ/mol, showing the reverse cis effect. The HFO‐1234ze molecules interact predominantly via π‐π stacking of the double bonds and hydrogen bonding of CH⋯F. Although the E‐dimers involve less steric exchange‐repulsion, their binding energies are generally smaller by 2–6 kJ/mol than the Z‐dimers, in which the Coulombic multipole–multipole‐type electrostatic interactions and the long‐range dispersions are more profound. Local ordering in liquid occurs primarily via π stacking of the CHF‐end. The significant difference in boiling points of E‐ and Z‐HFO‐1234ze is clarified in terms of the dimerization free energy. Self‐diffusion coefficients and dielectric constants of HFO‐1234ze depend strongly on the E‐ and Z‐conformations.

Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields.

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly.

The Role of Counterions in Intermolecular Radical Coupling of Ru-bda Catalysts

Intermolecular radical coupling (also interaction of two metal centers I2M) is one of the main mechanisms for O–O bond formation in water oxidation catalysts. For Ru(bda)L2 (H2bda = 2,2′-bipyridine-6,6′-dicarboxylate, L = pyridine or similar nitrogen containing heterocyclic ligands) catalysts a significant driving force in water solution is the hydrophobic effects driven by the solvent. The same catalyst has been successfully employed to generate N2 from ammonia, also via I2M, but here the solvent was acetonitrile where hydrophobic effects are absent. We used a classical force field for the key intermediate [RuVIN(bda)(py)2]+ to simulate the dimerization free energy by calculation of the potential mean force, in both water and acetonitrile to understand the differences and similarities. In both solvents the complex dimerizes with similar free energy profiles. In water the complexes are essentially free cations with limited ion paring, while in acetonitrile the ion-pairing is much more significant. This ion-pairing leads to significant screening of the charges, making dimerization possible despite lower solvent polarity that could lead to repulsion between the charged complexes. In water the lower ion pairing is compensated by the hydrophobic effect leading to favorable dimerization despite repulsion of the charges. A hypothetical doubly charged [RuVIIN(bda)py2]2+ was also studied for deeper understanding of the charge effect. Despite the double charge the complexes only dimerized favorably in the lower dielectric solvent acetonitrile, while in water the separated state is more stable. In the doubly charged catalyst the effect of ion-pairing is even more pronounced in acetonitrile where it is fully paired similar to the 1+ complex, while in water the separation of the ions leads to greater repulsion between the two catalysts, which prevents dimerization.Graphic

Addressing the Excessive Aggregation of Membrane Proteins in the MARTINI Model.

The MARTINI model is a widely used coarse-grained force field popular for its capacity to represent a diverse array of complex biomolecules. However, efforts to simulate increasingly realistic models of membranes, involving complex lipid mixtures and multiple proteins, suggest that membrane protein aggregates are overstabilized by the MARTINI v2.2 force field. In this study, we address this shortcoming of the MARTINI model. We determined the free energy of dimerization of four transmembrane protein systems using the nonpolarizable MARTINI model. Comparison with experimental FRET-based estimates of the dimerization free energy was used to quantify the significant overstabilization of each protein homodimer studied. To improve the agreement between simulation and experiment, a single uniform scaling factor, α, was used to enhance the protein-lipid Lennard-Jones interaction. A value of α = 1.04-1.045 was found to provide the best fit to the dimerization free energies for the proteins studied while maintaining the specificity of contacts at the dimer interface. To further validate the modified force field, we performed a multiprotein simulation using both MARTINI v2.2 and the reparameterized MARTINI model. While the original MARTINI model predicts oligomerization of protein into a single aggregate, the reparameterized MARTINI model maintains a dynamic equilibrium between monomers and dimers as predicted by experimental studies. The proposed reparameterization is an alternative to the standard MARTINI model for use in simulations of realistic models of a biological membrane containing diverse lipids and proteins.

Dimeric states of transmembrane domains of insulin and IGF-1R receptors: Structures and possible role in activation

Despite the biological significance of insulin signaling, the molecular mechanisms of activation of the insulin receptor (IR) and other proteins from its family remain elusive. Current hypothesis on signal transduction suggests ligand-triggered structural changes in the extracellular domain followed by transmembrane (TM) domains closure and dimerization leading to trans-autophosphorylation and kinase activity in intracellular segments of the receptor. Using NMR spectroscopy, we detected dimerization of isolated TM segments of IR in different membrane-mimicking environments and observed multiple signals of NH groups of protein backbone possibly corresponding to several dimer conformations. Taking available experimental data as constraints, several atomistic models of dimeric TM domains of IR and insulin-like growth factor 1 (IGF-1R) receptors were elaborated. Molecular dynamics simulations of IR ectodomain revealed noticeable collective movements potentially responsible for closure of the C-termini of FnIII-3 domains and spatial approaching of TM helices upon insulin-induced receptor activation. In addition, we demonstrated that the intracellular part of the receptor does not impose restrictions on the positioning of TM helices in the membrane. Finally, we used two independent structure prediction methods to generate a series of dimer conformations followed by their cluster analysis and dimerization free energy estimation to select the best dimer models. Biological relevance of the later was further tested via comparison of the hydrophobic organization of TM helices for both wild-type receptors and their mutants. Based on these data, the ability of several segments from other proteins to functionally replace IR and/or IGF-1R TM domains was explained.

Why do G-quadruplexes dimerize through the 5’-ends? Driving forces for G4 DNA dimerization examined in atomic detail

G-quadruplexes (G4) are secondary structures formed by guanine-rich nucleic acid sequences and shown to exist in living cells where they participate in regulation of gene expression and chromosome maintenance. G-quadruplexes with solvent-exposed guanine tetrads show the tendency to associate together through cofacial stacking, which may be important for packaging of G4-forming sequences and allows for the design of higher-order G4 DNA structures. To understand the molecular driving forces for G4 association, here, we study the binding interaction between two parallel-stranded G-quadruplexes using all-atom molecular dynamics simulations. The predicted dimerization free energies show that direct binding through the 5'-G-tetrads is the most preferred of all possible end-to-end stacking orientations, consistently with all available experimental data. Decomposition of dimerization enthalpies in combination with simulations at varying ionic strength further indicate that the observed orientational preferences arise from a fine balance between the electrostatic repulsion of the sugar-phosphate backbones and favorable counterion binding at the dimeric interface. We also demonstrate how these molecular-scale findings can be used to devise means of controlling G4 dimerization equilibrium, e.g., by altering salt concentration and using G4-targeted ligands.

A Martini coarse-grained model of the calcein fluorescent dye

Calcein leakage assays are a standard experimental set-up for probing the extent of damage induced by external agents on synthetic lipid vesicles. The fluorescence signal associated with calcein release from liposomes is the signature of vesicle disruption, transient pore formation or vesicle fusion. This type of assay is widely used to test the membrane disruptive effect of biological macromolecules, such as proteins, antimicrobial peptides and RNA and is also used on synthetic nanoparticles with a polymer, metal or oxide core. Little is known about the effect that calcein and other fluorescent dyes may have on the properties of lipid bilayers, potentially altering their structure and permeability. Here we develop a coarse-grained model of calcein that is compatible with the Martini force field for lipids. We validate the model by comparing its dimerization free energy, aggregation behavior at different concentrations and interaction with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane to those obtained at atomistic resolution. Our coarse-grained description of calcein makes it suitable for the simulation of large calcein-filled liposomes and of their interactions with external agents, allowing for a direct comparison between simulations and experimental liposome leakage assays.