Correction to "Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models"

Molecular Mechanisms Underlying Medium-Chain Free Fatty Acid-Regulated Activity of the Phospholipase PlaF from Pseudomonas aeruginosa

PlaF is a membrane-bound phospholipase A1 from Pseudomonas aeruginosa that is involved in remodeling membrane glycerophospholipids (GPLs) and modulating virulence-associated signaling and metabolic pathways. Previously, we identified the role of medium-chain free fatty acids (FFAs) in inhibiting PlaF activity and promoting homodimerization, yet the underlying molecular mechanism remained elusive. Here, we used unbiased and biased molecular dynamics simulations and free energy computations to assess how PlaF interacts with FFAs localized in the water milieu surrounding the bilayer or within the bilayer and how these interactions regulate PlaF activity. Medium-chain FFAs localized in the upper bilayer leaflet can stabilize inactive dimeric PlaF, likely through interactions with charged surface residues, as has been experimentally validated. Potential of mean force (PMF) computations indicate that membrane-bound FFAs may facilitate the activation of monomeric PlaF by lowering the activation barrier for changing into a tilted, active configuration. We estimated that the coupled equilibria of PlaF monomerization-dimerization and tilting at the physiological concentration of PlaF lead to the majority of PlaF forming inactive dimers when in a cell membrane loaded with decanoic acid (C10). This is in agreement with a suggested in vivo product feedback loop and gas chromatography-mass spectrometry profiling results, indicating that PlaF catalyzes the release of C10 from P. aeruginosa membranes. Additionally, we found that C10 in the water milieu can access the catalytic site of active monomeric PlaF, contributing to the competitive component of C10-mediated PlaF inhibition. Our study provides mechanistic insights into how medium-chain FFAs may regulate the activity of PlaF, a potential bacterial drug target.

Design and Characterization of Compact, Programmable, Multistranded Nonimmunostimulatory Nucleic Acid Nanoparticles Suitable for Biomedical Applications

We report a thorough investigation of the role of single-stranded thymidine (ssT) linkers in the stability and flexibility of minimal, multistranded DNA nanostructures. We systematically explore the impact of varying the number of ssTs in three-way junction motifs (3WJs) on their formation and properties. Through various UV melting experiments and molecular dynamics simulations, we demonstrate that while the number of ssTs minimally affects thermodynamic stability, the increasing ssT regions significantly enhance the structural flexibility of 3WJs. Utilizing this knowledge, we design triangular DNA nanoparticles with varying ssTs, all showing exceptional assembly efficiency except for the 0T triangle. All triangles demonstrate enhanced stability in blood serum and are nonimmunostimulatory and nontoxic in mammalian cell lines. The 4T 3WJ is chosen as the building block for constructing other polygons due to its enhanced flexibility and favorable physicochemical characteristics, making it a versatile choice for creating cost-effective, stable, and functional DNA nanostructures that can be stored in the dehydrated forms while retaining their structures. Our study provides valuable insights into the design and application of nucleic acid nanostructures, emphasizing the importance of understanding stability and flexibility in the realm of nucleic acid nanotechnology. Our findings suggest the intricate connection between these ssTs and the structural adaptability of DNA 3WJs, paving the way for more precise design and engineering of nucleic acid nanosystems suitable for broad biomedical applications.

DECTIN-1: A modifier protein in CTLA-4 haploinsufficiency

Autosomal dominant loss-of-function (LoF) variants in cytotoxic T-lymphocyte associated protein 4 (CTLA4) cause immune dysregulation with autoimmunity, immunodeficiency and lymphoproliferation (IDAIL). Incomplete penetrance and variable expressivity are characteristic of IDAIL caused by CTLA-4 haploinsufficiency (CTLA-4h), pointing to a role for genetic modifiers. Here, we describe an IDAIL proband carrying a maternally inherited pathogenic CTLA4 variant and a paternally inherited rare LoF missense variant in CLEC7A, which encodes for the β-glucan pattern recognition receptor DECTIN-1. The CLEC7A variant led to a loss of DECTIN-1 dimerization and surface expression. Notably, DECTIN-1 stimulation promoted human and mouse regulatory T cell (Treg) differentiation from naïve αβ and γδ T cells, even in the absence of transforming growth factor-β. Consistent with DECTIN-1's Treg-boosting ability, partial DECTIN-1 deficiency exacerbated the Treg defect conferred by CTL4-4h. DECTIN-1/CLEC7A emerges as a modifier gene in CTLA-4h, increasing expressivity of CTLA4 variants and acting in functional epistasis with CTLA-4 to maintain immune homeostasis and tolerance.

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

HYPOTHESIS

The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). We showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers.

EXPERIMENTS

First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances.

FINDINGS

Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential

Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/E_b, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C₄ term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C₄ term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C₄ term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.