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

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

• MeSH
• glykosaminoglykany * chemie   MeSH
• kyseliny sulfonové chemie   MeSH
• simulace molekulární dynamiky * MeSH
• sírany * chemie   MeSH
• statická elektřina * MeSH
• vápník chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• glykosaminoglykany * MeSH
• kyseliny sulfonové MeSH
• sírany * MeSH
• sulfamic acid MeSH Prohlížeč
• vápník 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

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

Colloidal nanoparticles exhibit unique size-dependent properties differing from their bulk counterpart, which can be particularly relevant for catalytic applications. To optimize surface-mediated chemical reactions, the understanding of the microscopic structure of the nanoparticle-liquid interface is of paramount importance. Here we use polarimetric angle-resolved second harmonic scattering (AR-SHS) to determine surface potential values as well as interfacial water orientation of ∼100 nm diameter amorphous TiO2 nanoparticles dispersed in aqueous solutions, without any initial assumption on the distribution of interfacial charges. We find three regions of different behavior with increasing NaCl concentration. At very low ionic strengths (0-10 μM), the Na+ ions are preferentially adsorbed at the TiO2 surface as inner-sphere complexes. At low ionic strengths (10-100 μM), a distribution of counterions equivalent to a diffuse layer is observed, while at higher ionic strengths (>100 μM), an additional layer of hydrated condensed ions is formed. We find a similar behavior for TiO2 nanoparticles in solutions of different basic pH. Compared to identically sized SiO2 nanoparticles, the TiO2 interface has a higher affinity for Na+ ions, which we further confirm with molecular dynamics simulations. With its ability to monitor ion adsorption at the surface with micromolar sensitivity and changes in the surface potential, AR-SHS is a powerful tool to investigate interfacial properties in a variety of catalytic and photocatalytic applications.

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

The microscopic description of the interface of colloidal particles in solution is essential to understand and predict the stability of these systems, as well as their chemical and electrochemical reactivity. However, this description often relies on the use of simplified electrostatic mean field models for the structure of the interface, which give only theoretical estimates of surface potential values and do not provide properties related to the local microscopic structure, such as the orientation of interfacial water molecules. Here we apply polarimetric angle-resolved second harmonic scattering (AR-SHS) to 300 nm diameter SiO2 colloidal suspensions to experimentally determine both surface potential and interfacial water orientation as a function of pH and NaCl concentration. The surface potential values and interfacial water orientation change significantly over the studied pH and salt concentration range, whereas zeta-potential (ζ) values remain constant. By comparing the surface and ζ-potentials, we find a layer of hydrated condensed ions present in the high pH case, and for NaCl concentrations ≥1 mM. For milder pH ranges (pH < 11), as well as for salt concentrations <1 mM, no charge condensation layer is observed. These findings are used to compute the surface charge densities using the Gouy-Chapman and Gouy-Chapman-Stern models. Furthermore, by using the AR-SHS data, we are able to determine the preferred water orientation in the layer directly in contact with the silica interface. Molecular dynamics simulations confirm the experimental trends and allow deciphering of the contributions of water layers to the total response.

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