Explicit Representation of Cation-π Interactions in Force Fields with 1/r4 Nonbonded Terms

Cation-π Interactions Greatly Influence Ion Transportability of the Light-Driven Sodium Pump KR2: A Molecular Dynamics Study

Krokinobacter eikastus rhodopsin 2 (KR2) is a typical light-driven sodium pump. Although wild-type KR2 exhibits high Na+ selectivity, mutagenesis performed on the residues constituting the entrance enables permeation of K+ and Cs+, while the underlying mechanism remains elusive. This study presents a comprehensive molecular dynamics investigation, including force field optimization, metadynamics, and alchemical free energy methods, to explore the N61L/G263F mutant of KR2, which exhibits transportability for K+ and Cs+. The introduced Phe263 residue can directly promote ion binding at the entrance through cation-π interactions, while the N61L mutation can enhance ion binding at Phe46 by relieving steric hindrance. These results suggest that cation-π interactions may significantly influence the ion transportability and selectivity of KR2, which can provide important insights for protein engineering and the design of artificial ion transporters.

OPLS/2020 Force Field for Unsaturated Hydrocarbons, Alcohols, and Ethers

The OPLS all-atom force field was updated and applied to modeling unsaturated hydrocarbons, alcohols, and ethers. Testing has included gas-phase conformational energetics, properties of pure liquids, and free energies of hydration. Monte Carlo statistical mechanics (MC) calculations were used to model 60 liquids. In addition, a robust, automated procedure was devised to compute the free energies of hydration with high precision via free-energy perturbation (FEP) calculations using double annihilation. Testing has included larger molecules than in the past, and parameters are reported for the first time for some less common groups including alkynes, allenes, dienes, and acetals. The average errors in comparison with experimental data for the computed properties of the pure liquids were improved with the modified force field (OPLS/2020). For liquid densities and heats of vaporization, the average unsigned errors are 0.01 g/cm3 and 0.2 kcal/mol. The average error and signed error for free energies of hydration are both 1.2 kcal/mol. As noted before, this reflects a systematic overestimate of the hydrophobicity of organic molecules when the parametrization is done to minimize the errors for properties of pure liquids. Implications for the modeling of biomolecular systems with standard force fields are considered.

Linear and ring polypeptides complexed with oppositely charged surfactants: the cohesion of the complexes as revealed in atomistic simulations

The use of linear supercharged unfolded polypeptides (SUPs) and oppositely charged surfactants in aqueous solution has demonstrated impressive adhesive properties. These substances possess biocompatibility, biodegradability and other necessary properties for practical application as a biomedical glue in wound repair. The success of these substances, coupled with limited knowledge about such systems, provides hope for enhancing the performance of the final product. One potential approach involves altering the topology of the polypeptide chain. In this article, we conduct a comparative analysis to examine the behavior of the ring and linear chains of a polypeptide in aqueous solution. This analysis utilizes full-atomic computer modeling to monitor the properties of the chains. We investigate the temperature dependence of the shape and size of individual polypeptides in the solution, as well as the formation of complexes via mixing the polypeptide chains with oppositely charged sodium dodecylbenzene sulfonate (SDBS) surfactant molecules in a stoichiometric ratio. Additionally, we explore the cohesive properties of the resulting complex through power experiments involving the extraction of single polypeptide chains out of the SUP-SDBS complexes.

Progress on the Pharmacological Targeting of Janus Pseudokinases

The Janus kinases (JAKs) are key components of the JAK-STAT signaling pathway and are involved in myriad physiological processes. Though they are the molecular targets of many FDA-approved drugs, these drugs manifest adverse effects due in part to their inhibition of the requisite JAK kinase activity. However, the JAKs uniquely possess an integrated pseudokinase domain (JH2) that regulates the adjacent kinase domain (JH1). The therapeutic targeting of JH2 domains has been less thoroughly explored and may present an avenue to modulate the JAKs without the adverse effects associated with targeting the adjacent JH1 domain. The potential of this strategy was recently demonstrated with the FDA approval of the TYK2 JH2 ligand deucravacitinib for treating plaque psoriasis. In this light, the structure and targetability of the JAK pseudokinases are discussed, in conjunction with the state of development of ligands that bind to these domains.

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides

To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.

Information-Rich, Dual-Function 13C/2H-Isotopic Crosstalk NMR Assay for Human Serine Racemase (hSR) Provides a PLP-Enzyme "Partitioning Fingerprint" and Reveals Disparate Chemotypes for hSR Inhibition

The first dual-function assay for human serine racemase (hSR), the only bona fide racemase in human biology, is reported. The hSR racemization function is essential for neuronal signaling, as the product, d-serine (d-Ser), is a potent N-methyl d-aspartate (NMDA) coagonist, important for learning and memory, with dysfunctional d-Ser-signaling being observed in some neuronal disorders. The second hSR function is β-elimination and gives pyruvate; this activity is elevated in colorectal cancer. This new NMR-based assay allows one to monitor both α-proton-exchange chemistry and β-elimination using only the native l-Ser substrate and hSR and is the most sensitive such assay. The assay judiciously employs segregated dual 13C-labeling and 13C/2H crosstalk, exploiting both the splitting and shielding effects of deuterium. The assay is deployed to screen a 1020-compound library and identifies an indolo-chroman-2,4-dione inhibitor family that displays allosteric site binding behavior (noncompetitive inhibition vs l-Ser substrate; competitive inhibition vs adenosine 5'-triphosphate (ATP)). This assay also reveals important mechanistic information for hSR; namely, that H/D exchange is ∼13-fold faster than racemization, implying that K56 protonates the carbanionic intermediate on the si-face much faster than does S84 on the re-face. Moreover, the 13C NMR peak pattern seen is suggestive of internal return, pointing to K56 as the likely enamine-protonating residue for β-elimination. The 13C/2H-isotopic crosstalk assay has also been applied to the enzyme tryptophan synthase and reveals a dramatically different partition ratio in this active site (β-replacement: si-face protonation ∼6:1 vs β-elimination: si-face protonation ∼1:3.6 for hSR), highlighting the value of this approach for fingerprinting the pyridoxal phosphate (PLP) enzyme mechanism.

Refinement of the Optimized Potentials for Liquid Simulations Force Field for Thermodynamics and Dynamics of Liquid Alkanes

Torsion and Lennard-Jones parameters of the optimized potentials for liquid simulations (OPLS) all-atom force field have been refined for describing thermodynamics and dynamics of a wide range of liquid alkanes. Monte Carlo statistical mechanics (MC) and molecular dynamics (MD) simulations were carried out. For thermodynamics properties, MC simulations with truncated electrostatic interactions performed very closely to MD simulations with a Verlet neighbor list and the particle mesh Ewald algorithm. The average errors in comparison with experimental data for computed properties were improved with the modified force field (OPLS/2020), especially for long-chain alkanes. For liquid densities, heats of vaporization, and free energies of hydration, the average errors are 0.01 g/cm³, 0.2 kcal/mol, and ca. 0.5 kcal/mol, respectively; significant gains were made for relative heats of vaporization of isomeric series. Results for self-diffusion coefficients also reproduce experimental data well for linear alkane liquids up to hexadecane. The new force field is suitable for use in improved modeling of myriad systems of importance in chemistry, biology, and materials science.

Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model

In this work, a set of 12-6-4 force fields (FFs) parameters were developed for the actinyl molecular cations, AnO2n+ (n = 1, 2), from uranium to plutonium for classical molecular dynamics (MD) for four water models: TIP3P, SPC/E, OPC3, and TIP4Pew. Such a non-bonded potential model taking into account the induced dipole between the metallic center and the surrounding molecules has shown better performances for various cations than the classic 12-6 non-bonded potentials. The parametrization method proposed elsewhere for metallic cations has been extended to these molecular cations. In contrast to the actinyl 12-6 FFs from the literature, the new models reproduce correctly both solvation and thermodynamic properties, thanks to the inclusion of the induced dipole term (C4). The transferability of such force fields was assessed by performing MD simulations of carbonato actinyl species, which are highly implicated in actinide migration or actinide extraction from seawater. A highly satisfying agreement was found when comparing the EXAFS signals computed from our MD simulation to the experimental ones. The set of FFs developed here opens new possibilities for the study of actinide chemistry.

Conversion of a False Virtual Screen Hit into Selective JAK2 JH2 Domain Binders Using Convergent Design Strategies

The Janus kinase 2 (JAK2) pseudokinase domain (JH2) is an ATP-binding domain that regulates the activity of the catalytic tyrosine kinase domain (JH1). Dysregulation of JAK2 JH1 signaling caused by the V617F mutation in JH2 is implicated in various myeloproliferative neoplasms. To explore if JAK2 activity can be modulated by a small molecule binding to the ATP site in JH2, we have developed several ligand series aimed at selectively targeting the JAK2 JH2 domain. We report here the evolution of a false virtual screen hit into a new JAK2 JH2 series. Optimization guided by computational modeling has yielded analogues with nanomolar affinity for the JAK2 JH2 domain and >100-fold selectivity for the JH2 domain over the JH1 domain. A crystal structure for one of the potent compounds bound to JAK2 JH2 clarifies the origins of the strong binding and selectivity. The compounds expand the platform for seeking molecules to regulate JAK2 signaling, including V617F JAK2 hyperactivation.

Complexes of Alkaline and Ammonium Cations with Dopamine and Eumelanin Precursors: Dissecting the Role of Noncovalent Cation-π and Cation-Lone Pair (σ-Type) Interactions

Cation-π interactions and their possible competition with other noncovalent interactions (NCI) might play a key role in both dopamine- and eumelanin-based bioinspired materials. In this contribution, to unravel the delicate interplay between cation-π interactions and other possible competing forces, the configurational space of noncovalent complexes formed by dopamine or eumelanin precursors (o-benzoquinone, DHI and a semiquinone dimer) and three different cations (Na⁺, K⁺, and NH₄⁺) is sampled by means of accurate ab initio calculations. To this end, we resort to the mp2^mod method, recently validated by us for benzene-, phenol-, and catechol-cation complexes, whose computational convenience allows for an extensive exploration of the cation-molecule interaction energy surface, by sampling a total of more than 10⁴ arrangements. The mp2^mod interaction energy landscapes reveal that, besides the expected cation-π driven arrangements, for all considered molecule-cation pairs the most stable complexes are found when the cation lies within the plane containing the six-membered ring, thus maximizing the σ-type interaction with the oxygen's lone pairs. Due to the loss of aromaticity, the σ-type/cation-π strength ratio is remarkably large in o-benzoquinone, where cation-π complexes seem unlikely to be formed. The above features are shared among all considered cations but are significantly larger when considering the smaller Na⁺. Besides delivering a deeper insight onto the NCI network established by the considered precursors in the presence of ions, the present results can serve as a reference database to validate or refine lower level methods, as, for instance, the force fields employed in classical simulations.

Molecular interactions contributing to FUS SYGQ LC-RGG phase separation and co-partitioning with RNA polymerase II heptads

The RNA-binding protein FUS (Fused in Sarcoma) mediates phase separation in biomolecular condensates and functions in transcription by clustering with RNA polymerase II. Specific contact residues and interaction modes formed by FUS and the C-terminal heptad repeats of RNA polymerase II (CTD) have been suggested but not probed directly. Here we show how RGG domains contribute to phase separation with the FUS N-terminal low-complexity domain (SYGQ LC) and RNA polymerase II CTD. Using NMR spectroscopy and molecular simulations, we demonstrate that many residue types, not solely arginine-tyrosine pairs, form condensed-phase contacts via several interaction modes including, but not only sp2-π and cation-π interactions. In phases also containing RNA polymerase II CTD, many residue types form contacts, including both cation-π and hydrogen-bonding interactions formed by the conserved human CTD lysines. Hence, our data suggest a surprisingly broad array of residue types and modes explain co-phase separation of FUS and RNA polymerase II.

Physisorption of bio oil nitrogen compounds onto montmorillonite

We assess computationally the adsorption of a series of nitrogen containing heterocycles and fatty acid amides from bio-oil on a model clay surface, Na-montmorillonite. The adsorption energies and conformations predicted by atomistic detail molecular dynamics (MD) simulations are compared against density functional theory (DFT) based molecular electrostatic potentials (MEP) and Hirshfeld, AIM, Merz-Singh-Kollman, and ChelpG charges. MD predicts systematically adsorption via cation bridging with adsorption strength of the heterocycles following purine > pyridine > imidazole > pyrrole > indole > quinoline. The fatty acid amides adsorption strength follows the steric availability and bulkiness of the head group. A comparison against the DFT calculations shows that MEP predicts adsorption geometries and the MD simulations reproduce the conformations for single adsorption site species. However, the DFT derived charge distibutions show that MD force-fields with non-polarizable fixed partial charge representations parametrized for aqueous environments cannot be used in apolar solvent environments without careful accuracy considerations. The overall trends in adsorption energies are reproduced by the Charmm GenFF employed in the MD simulations but the adsorption energies are systematically overestimated in this apolar solvent environment. The work has significance both for revealing nitrogen compound adsorption trends in technologically relevant bio oil environments but also as a methodological assessment revealing the limits of state of the art biomolecular force-fields and simulation protocols in apolar bioenvironments.

Cation-π Interactions and their Functional Roles in Membrane Proteins

Cation-π interactions arise as a result of strong attractive forces between positively charged entities and the π-electron cloud of aromatic groups. The physicochemical characteristics of cation-π interactions are particularly well-suited to the dual hydrophobic/hydrophilic environment of membrane proteins. As high-resolution structural data of membrane proteins bring molecular features into increasingly sharper view, cation-π interactions are gaining traction as essential contributors to membrane protein chemistry, function, and pharmacology. Here we review the physicochemical properties of cation-π interactions and present several prominent examples which demonstrate significant roles for this specialized biological chemistry.

Indoloxytriazines as binding molecules for the JAK2 JH2 pseudokinase domain and its V617F variant

Small molecules that selectively bind to the pseudokinase JH2 domain over the JH1 kinase domain of JAK2 kinase are sought. Virtual screening led to the purchase of 17 compounds among which 9 were found to bind to V617F JAK2 JH2 with affinities of 40 - 300 μM in a fluorogenic assay. Ten analogues were then purchased yielding 9 additional active compounds. Aminoanilinyltriazine 22 was particularly notable as it shows no detectable binding to JAK2 JH1, and it has a 65-μM dissociation constant K _d with V617F JAK2 JH2. A crystal structure for 22 in complex with wild-type JAK2 JH2 was obtained to elucidate the binding mode. Additional de novo design led to the synthesis of 19 analogues of 22 with the most potent being 33n with K _d values of 2-3 μM for WT and V617F JAK2 JH2, and with 16-fold selectivity relative to binding with WT JAK2 JH1.

Accuracy of Alternate Nonpolarizable Force Fields for the Determination of Protein-Ligand Binding Affinities Dominated by Cation-π Interactions

Modifying pair-specific Lennard-Jones parameters through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field has proven cost-effective for improving the description of cation-π interactions in biological objects by means of pairwise additive potential energy functions. Here, two sets of newly optimized CHARMM36 force-field parameters including NBFIX corrections, coined CHARMM36m-NBF and CHARMM36-WYF, and the original force fields, namely CHARMM36m and Amber ff14SB, are used to determine the standard binding free energies of seven protein-ligand complexes containing cation-π interactions. Compared with precise experimental measurements, our results indicate that the uncorrected, original force fields significantly underestimate the binding free energies, with a mean error of 5.3 kcal/mol, while the mean errors of CHARMM36m-NBF and CHARMM36-WYF amount to 0.8 and 2.1 kcal/mol, respectively. The present study cogently demonstrates that the use of modified parameters jointly with NBFIX corrections dramatically increases the accuracy of the standard binding free energy of protein-ligand complexes dominated by cation-π interactions, most notably with CHARMM36m-NBF.

Potency- and Selectivity-Enhancing Mutations of Conotoxins for Nicotinic Acetylcholine Receptors Can Be Predicted Using Accurate Free-Energy Calculations

Nicotinic acetylcholine receptor (nAChR) subtypes are key drug targets, but it is challenging to pharmacologically differentiate between them because of their highly similar sequence identities. Furthermore, α-conotoxins (α-CTXs) are naturally selective and competitive antagonists for nAChRs and hold great potential for treating nAChR disorders. Identifying selectivity-enhancing mutations is the chief aim of most α-CTX mutagenesis studies, although doing so with traditional docking methods is difficult due to the lack of α-CTX/nAChR crystal structures. Here, we use homology modeling to predict the structures of α-CTXs bound to two nearly identical nAChR subtypes, α3β2 and α3β4, and use free-energy perturbation (FEP) to re-predict the relative potency and selectivity of α-CTX mutants at these subtypes. First, we use three available crystal structures of the nAChR homologue, acetylcholine-binding protein (AChBP), and re-predict the relative affinities of twenty point mutations made to the α-CTXs LvIA, LsIA, and GIC, with an overall root mean square error (RMSE) of 1.08 ± 0.15 kcal/mol and an R2 of 0.62, equivalent to experimental uncertainty. We then use AChBP as a template for α3β2 and α3β4 nAChR homology models bound to the α-CTX LvIA and re-predict the potencies of eleven point mutations at both subtypes, with an overall RMSE of 0.85 ± 0.08 kcal/mol and an R2 of 0.49. This is significantly better than the widely used molecular mechanics-generalized born/surface area (MM-GB/SA) method, which gives an RMSE of 1.96 ± 0.24 kcal/mol and an R2 of 0.06 on the same test set. Next, we demonstrate that FEP accurately classifies α3β2 nAChR selective LvIA mutants while MM-GB/SA does not. Finally, we use FEP to perform an exhaustive amino acid mutational scan of LvIA and predict fifty-two mutations of LvIA to have greater than 100X selectivity for the α3β2 nAChR. Our results demonstrate the FEP is well-suited to accurately predict potency- and selectivity-enhancing mutations of α-CTXs for nAChRs and to identify alternative strategies for developing selective α-CTXs.

Benchmarking London dispersion corrected density functional theory for noncovalent ion-π interactions

The strongly attractive noncovalent interactions of charged atoms or molecules with π-systems are important binding motifs in many chemical and biological systems. These so-called ion-π interactions play a major role in enzymes, molecular recognition, and for the structure of proteins. In this work, a molecular test set termed IONPI19 is compiled for inter- and intramolecular ion-π interactions, which is well balanced between anionic and cationic systems. The IONPI19 set includes interaction energies of significantly larger molecules (up to 133 atoms) than in other ion-π test sets and covers a broad range of binding motifs. Accurate (local) coupled cluster values are provided as reference. Overall, 19 density functional approximations, including seven (meta-)GGAs, eight hybrid functionals, and four double-hybrid functionals combined with three different London dispersion corrections, are benchmarked for interaction energies. DFT results are further compared to wave function based methods such as MP2 and dispersion corrected Hartree-Fock. Also, the performance of semiempirical QM methods such as the GFNn-xTB and PMx family of methods is tested. It is shown that dispersion-uncorrected DFT underestimates ion-π interactions significantly, even though electrostatic interactions dominate the overall binding. Accordingly, the new charge dependent D4 dispersion model is found to be consistently better than the standard D3 correction. Furthermore, the functional performance trend along Jacob's ladder is generally obeyed and the reduction of the self-interaction error leads to an improvement of (double) hybrid functionals over (meta-)GGAs, even though the effect of the SIE is smaller than expected. Overall, the double-hybrids PWPB95-D4/QZ and revDSD-PBEP86-D4/QZ turned out to be the most reliable among all assessed methods for the description of ion-π interactions, which opens up new perspectives for systems where coupled cluster calculations are no longer computationally feasible.

Biomolecular Condensates: Sequence Determinants of Phase Separation, Microstructural Organization, Enzymatic Activity, and Material Properties

This perspective article highlights recent progress and emerging challenges in understanding the formation and function of membraneless organelles (MLOs). A long-term goal in the MLO field is to identify the sequence-encoded rules that dictate the formation of compositionally controlled biomolecular condensates, which cells utilize to perform a wide variety of functions. The molecular organization of the different components within a condensate can vary significantly, ranging from a homogeneous mixture to core-shell droplet structures. We provide many examples to highlight the richness of the observed behavior and potential research directions for improving our mechanistic understanding. The tunable environment within condensates can, in principle, alter enzymatic activity significantly. We examine recent examples where this was demonstrated, including applications in synthetic biology. An important question about MLOs is the role of liquid-like material properties in biological function. We discuss the need for improved quantitative characterization tools and the development of sequence-structure-dynamics relationships.

The phenoxyl group-modulated interplay of cation-π and σ-type interactions in the alkali metal series

The interaction potential energy surfaces (IPESs) of four alkaline metal cations (Na+, K+, Rb+ and Cs+) complexed with phenol and catechol were explored by accurate ab initio calculations to investigate the interplay of different noncovalent interactions and their behavior along the alkali metal series and upon -OH substitution. Selected one-dimensional interaction energy curves revealed two different minimum energy configurations for all phenol- and catechol-metal complexes, characterized either by cation-π or σ-type interactions. For each investigated complex several two-dimensional IPES maps were also computed, exploiting the computational advantages of the MP2mod approach. The size of the alkali cation was found to play a similar role in modulating both kinds of complexes, as the interaction strength always decreases along the metal series, from Na+ to Cs+. Conversely, the number of hydroxyl substituents markedly affected cation-π complexes vs. σ-type ones. As a most relevant finding, in catechol-metal complexes the strength of cation-π interactions is around half that of the σ-type ones. It is argued that the combined effect of cation dimensions and hydroxyl substitution in catechol-Na+ complexes makes σ-type configurations remarkably more stable and easily accessible than cation-π ones. Besides shedding new light on the origin of biological phenomena connected with underwater adhesion, the quantum mechanical interaction energy database provided herein may offer a useful reference for tuning accurate force fields, suitable for molecular dynamics simulations, where environmental effects might be also taken into account.