Salt bridges are ionic interactions that are of great importance in protein recognition. However, their structural description using X-ray crystallography or NMR may be inconclusive. Classical molecular dynamics (MD) used for the interpretation neglects electronic polarization, which results in artifactual overbinding. Here, we resolve the problem via charge scaling, which accounts for electronic polarization in a mean-field way. We study three salt bridges in insulin analogue. New NMR ensembles are generated via NOE-restrained MD using ff19SB and CHARMM36m force fields and the scaled-charge prosECCo75. Tens of μs of unrestrained MD show in a statistically converged manner that ff19SB induces a non-native salt bridge. This behavior is quantified via umbrella sampling of salt bridge dissociation, which indicates a rather high strength of up to 4 and 5 kcal mol-1 for CHARMM36m and ff19SB, respectively. In contrast, prosECCo75 gives a biologically reasonable dissociation barrier of 1 kcal mol-1. Our results indicate that a physically justified description of charge-charge interactions within a nonpolarizable MD framework reliably describes aqueous biomolecular systems.

• Publikační typ
• časopisecké články MeSH

Additives in vaping products, such as flavors, preservatives, or thickening agents, are commonly used to enhance user experience. Among these, Vitamin E acetate (VEA) was initially thought to be harmless but has been implicated as the primary cause of e-cigarette or vaping product use-associated lung injury, a serious lung disease. In our study, VEA serves as a proxy for other e-cigarette additives. To explore its harmful effects, we developed an exposure system to subject a pulmonary surfactant (PSurf) model to VEA-rich vapor. Through detailed analysis and atomic-level simulations, we found that VEA tends to cluster into aggregates on the PSurf surface, inducing deformations and weakening its essential elastic properties, critical for respiratory cycle function. Apart from VEA, our experiments also indicate that propylene glycol and vegetable glycerin, widely used in e-liquid mixtures, or their thermal decomposition products, alter surfactant properties. This research provides molecular-level insights into the detrimental impacts of vaping product additives on lung health.

• Klíčová slova
• EVALI, Lung surfactant, Molecular dynamics simulation, Pulmonary surfactant, Vaping-associated pulmonary injury,
• MeSH
• acetáty analýza   chemie   MeSH
• biologické modely MeSH
• lidé MeSH
• plicní surfaktanty * chemie   MeSH
• propylenglykol chemie   MeSH
• systémy dodávající nikotin elektronicky * MeSH
• vaping * škodlivé účinky   MeSH
• vitamin E * MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• acetáty MeSH
• plicní surfaktanty * MeSH
• propylenglykol MeSH
• vitamin E * MeSH

prosECCo75 is an optimized force field effectively incorporating electronic polarization via charge scaling. It aims to enhance the accuracy of nominally nonpolarizable molecular dynamics simulations for interactions in biologically relevant systems involving water, ions, proteins, lipids, and saccharides. Recognizing the inherent limitations of nonpolarizable force fields in precisely modeling electrostatic interactions essential for various biological processes, we mitigate these shortcomings by accounting for electronic polarizability in a physically rigorous mean-field way that does not add to computational costs. With this scaling of (both integer and partial) charges within the CHARMM36 framework, prosECCo75 addresses overbinding artifacts. This improves agreement with experimental ion binding data across a broad spectrum of systems─lipid membranes, proteins (including peptides and amino acids), and saccharides─without compromising their biomolecular structures. prosECCo75 thus emerges as a computationally efficient tool providing enhanced accuracy and broader applicability in simulating the complex interplay of interactions between ions and biomolecules, pivotal for improving our understanding of many biological processes.

• MeSH
• peptidy chemie   MeSH
• proteiny chemie   MeSH
• simulace molekulární dynamiky * MeSH
• statická elektřina * MeSH
• voda chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• peptidy MeSH
• proteiny MeSH
• voda MeSH

Investigating the electrical double layer (EDL) structure has been a long-standing challenge and has seen the emergence of several sophisticated techniques able to probe selectively the few molecular layers of a solid/water interface. While a qualitative estimation of the thickness of the EDL can be obtained using simple theoretical models, following experimentally its evolution is not straightforward and can be even more complicated in nano- or microscale systems, particularly when changing the ionic concentration by several orders of magnitude. Here, we bring insight into the structure of the EDL of SiO2 nanoparticle suspensions and its evolution with increasing ionic concentration using angle-resolved second harmonic scattering (AR-SHS). Below millimolar salt concentrations, we can successively characterize inner-sphere adsorption, diffuse layer formation, and outer-sphere adsorption. Moreover, we show for the first time that, by appropriately selecting the nanoparticle size, it is possible to retrieve information also in the millimolar range. There, we observe a decrease in the magnitude of the surface potential corresponding to a compression in the EDL thickness, which agrees with the results of several other electroanalytical and optical techniques. Molecular dynamics simulations suggest that the EDL compression mainly results from the diffuse layer compression rather than outer-sphere ions (Stern plane) moving closer to the surface.

• Publikační typ
• časopisecké články MeSH

The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

• Publikační typ
• časopisecké články MeSH

We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data.

• Publikační typ
• časopisecké články MeSH

Adsorption of arginine-rich positively charged peptides onto neutral zwitterionic phosphocholine (PC) bilayers is a key step in the translocation of those potent cell-penetrating peptides into the cell interior. In the past, we have shown both theoretically and experimentally that polyarginines adsorb to the neutral PC-supported lipid bilayers in contrast to polylysines. However, comparing our results with previous studies showed that the results often do not match even at the qualitative level. The adsorption of arginine-rich peptides onto 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) may qualitatively depend on the actual experimental conditions where binding experiments have been performed. In this work, we systematically studied the adsorption of R9 and K9 peptides onto the POPC bilayer, aided by molecular dynamics (MD) simulations and fluorescence cross-correlation spectroscopy (FCCS) experiments. Using MD simulations, we tested a series of increasing peptide concentrations, in parallel with increasing Na+ and Ca2+ salt concentrations, showing that the apparent strength of adsorption of R9 decreases upon the increase of peptide or salt concentration in the system. The key result from the simulations is that the salt concentrations used experimentally can alter the picture of peptide adsorption qualitatively. Using FCCS experiments with fluorescently labeled R9 and K9, we first demonstrated that the binding of R9 to POPC is tighter by almost 2 orders of magnitude compared to that of K9. Finally, upon the addition of an excess of either Na+ or Ca2+ ions with R9, the total fluorescence correlation signal is lost, which implies the unbinding of R9 from the PC bilayer, in agreement with our predictions from MD simulations.

• MeSH
• adsorpce MeSH
• arginin MeSH
• fosfatidylcholiny chemie   MeSH
• fosforylcholin MeSH
• lecitiny MeSH
• lipidové dvojvrstvy * chemie   MeSH
• osmolární koncentrace MeSH
• penetrační peptidy * chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• arginin MeSH
• fosfatidylcholiny MeSH
• fosforylcholin MeSH
• lecitiny MeSH
• lipidové dvojvrstvy * MeSH
• penetrační peptidy * MeSH

Electrostatic interactions have a determining role in the conformational and dynamic behavior of polyelectrolyte molecules. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in a water solution with the most commonly used K+ or Na+ counterions, were investigated using atomistic molecular dynamics (MD) simulations. We performed a comparison of seven popular force fields, namely AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and using two common corrections for overbinding of ions, the non-bonded fix (NBFIX), and electronic continuum corrections (ECC). These corrections were originally introduced to correct for the often-reported problem concerning the overbinding of ions to the charged groups of polyelectrolytes. In this work, a comparison of the simulation results with existing experimental data revealed several differences between the investigated force fields. The data from these simulations and comparisons with previous experimental data were then used to determine the limitations and strengths of these force fields in the context of the structural and dynamic properties of anionic polyamino acids. Physical properties, such as molecular sizes, local structure, and dynamics, were studied using two types of common counterions, namely potassium and sodium. The results show that, in some cases, both the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of the ions and charged side chain groups. The local structures and dynamics are more sensitive to dihedral angle parameterization, resulting in a preference for defined monomer conformations and the type of correction used. We also provide recommendations based on the results.

• Klíčová slova
• carboxyls, counterions, force fields, ions, molecular dynamics, peptides and proteins,
• Publikační typ
• časopisecké články MeSH

Ion-specific effects play a crucial role in controlling the stability of colloidal systems and regulating interfacial processes. Although mechanistic pictures have been developed to explain the electrostatic structure of solid/water colloidal interfaces, ion-specific effects remain poorly understood. Here we quantify the average interfacial water orientation and the electrostatic surface potential around 100 nm SiO2 and TiO2 colloidal particles in the presence of NaCl, RbCl, and CaCl2 using polarimetric angle-resolved second harmonic scattering. We show that these two parameters can be used to establish the ion adsorption mechanism in a low ionic strength regime (<1 mM added salt). The relative differences between salts as a function of the ionic strength demonstrate cation- and surface-specific preferences for inner- vs outer-sphere adsorption. Compared to monovalent Rb+ and Na+, Ca2+ is found to be preferentially adsorbed as outer-sphere on SiO2 surfaces, while a dominant inner-sphere adsorption is observed for Ca2+ on TiO2. Molecular dynamics simulations performed on crystalline SiO2 and TiO2 surfaces support the experimental conclusions. This work contributes to the understanding of the electrostatic environment around colloidal nanoparticles on a molecular level by providing insight into ion-specific effects with micromolar sensitivity.

• Publikační typ
• časopisecké články MeSH

The Ca2+ sensor STIM1 and the Ca2+ channel Orai1 that form the store-operated Ca2+ (SOC) channel complex are key targets for drug development. Selective SOC inhibitors are currently undergoing clinical evaluation for the treatment of auto-immune and inflammatory responses and are also deemed promising anti-neoplastic agents since SOC channels are linked with enhanced cancer cell progression. Here, we describe an investigation of the site of binding of the selective inhibitor Synta66 to the SOC channel Orai1 using docking and molecular dynamics simulations, and live cell recordings. Synta66 binding was localized to the extracellular site close to the transmembrane (TM)1 and TM3 helices and the extracellular loop segments, which, importantly, are adjacent to the Orai1-selectivity filter. Synta66-sensitivity of the Orai1 pore was, in fact, diminished by both Orai1 mutations affecting Ca2+ selectivity and permeation of Na+ in the absence of Ca2+. Synta66 also efficiently blocked SOC in three glioblastoma cell lines but failed to interfere with cell viability, division and migration. These experiments provide new structural and functional insights into selective drug inhibition of the Orai1 Ca2+ channel by a high-affinity pore blocker.

• Klíčová slova
• Ca2+, SOCE, GBM, Orai, STIM, Synta66, binding, docking, glioblastoma multiforme, pocket, pore,
• Publikační typ
• časopisecké články MeSH