Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent

Selective adsorption of divalent and trivalent cations in porous electrodes

The capacitive deionization technology uses the electrochemical adsorption of ions in porous electrodes to desalinate seawater or brackish water. Recently, capacitive deionization has gained significant attention as a technology for selective adsorption of ionic species from multicomponent aqueous electrolytes. To investigate the mechanism of selective adsorption at the molecular level, we performed molecular dynamics simulations of aqueous electrolytes and porous electrodes with different divalent or trivalent ions, electrode pore sizes, and applied voltages. We calculated the free energy barriers preventing ions from entering the pores of the electrode and the structure of the water molecules near the ions and the electrode surface under various conditions. Our results suggest that, when the pore and ion sizes are comparable, the steric and electrostatic interactions between the hydrated ions and electrode pores are comparable in magnitude. Moreover, the relative importance of the two interactions can be reversed by slight changes in the external conditions, such as the ion size, valence of the ions, electrode pore size, and applied voltage. Thus, by finely tuning the electrode pore size and the applied voltage, it may be possible to selectively adsorb a particular ionic species from a multicomponent electrolyte through capacitive deionization using a porous electrode.

BH3-mimetics or DNA-damaging agents in combination with RG7388 overcome p53 mutation-induced resistance to MDM2 inhibition

The development of drug resistance reduces the efficacy of cancer therapy. Tumor cells can acquire resistance to MDM2 inhibitors, which are currently under clinical evaluation. We generated RG7388-resistant neuroblastoma cells, which became more proliferative and metabolically active and were less sensitive to DNA-damaging agents in vitro and in vivo, compared with wild-type cells. The resistance was associated with a mutation of the p53 protein (His193Arg). This mutation abated its transcriptional activity via destabilization of the tetrameric p53-DNA complex and was observed in many cancer types. Finally, we found that Cisplatin and various BH3-mimetics could enhance RG7388-mediated apoptosis in RG7388-resistant neuroblastoma cells, thereby partially overcoming resistance to MDM2 inhibition.

The mechanisms underlying Li+/Mg2+ separation in ZIF-8 under an electric field from atomistic simulations

In this study, we explore the mass transfer and separation mechanism of Li+ and Mg2+ confined within the flexible nanoporous zeolite imidazolate framework ZIF-8 under the influence of an electric field, employing molecular dynamics simulation. Our results highlight that the electric field accelerates the dehydration process of ions and underscore the critical importance of ZIF-8 framework flexibility in determining the separation selectivity of the ZIF-8 membrane. The electric field is shown to diminish ion hydration in the confined space of ZIF-8, notably disrupting the orientation of water molecules in the first hydration shells of ions, leading to an asymmetrical ionic hydration structure characterized by the uniform alignment of water dipoles. Furthermore, despite the geometrical constraints imposed by the ZIF-8 framework, the electric field significantly enhances ionic mobility. Notably, the less stable hydration shell of Li+ facilitates its rapid, dehydration-induced transit through ZIF-8 nanopores, unlike Mg2+, whose stable hydration shell impedes dehydration. Further investigation into the structural characteristics of the six-ring windows traversed by Li+ and Mg2+ ions reveals distinct mechanisms of passage: for Mg2+ ions, significant window expansion is necessary, while for Li+ ions, the mechanism involves both window expansion and partial dehydration. These findings reveal the profound impact of the electric field and framework flexibility on the separation of Li+ and Mg2+, offering critical insights for the potential application of flexible nanoporous materials in the selective extraction of Li+ from salt-lake brine.

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

Machine-learned molecular mechanics force fields from large-scale quantum chemical data

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest.

Simulating Metal-Imidazole Complexes

One commonly observed binding motif in metalloproteins involves the interaction between a metal ion and histidine's imidazole side chains. Although previous imidazole-M(II) parameters established the flexibility and reliability of the 12-6-4 Lennard-Jones (LJ)-type nonbonded model by simply tuning the ligating atom's polarizability, they have not been applied to multiple-imidazole complexes. To fill this gap, we systematically simulate multiple-imidazole complexes (ranging from one to six) for five metal ions (Co(II), Cu(II), Mn(II), Ni(II), and Zn(II)) which commonly appear in metalloproteins. Using extensive (40 ns per PMF window) sampling to assemble free energy association profiles (using OPC water and standard HID imidazole charge models from AMBER) and comparing the equilibrium distances to DFT calculations, a new set of parameters was developed to focus on energetic and geometric features of multiple-imidazole complexes. The obtained free energy profiles agree with the experimental binding free energy and DFT calculated distances. To validate our model, we show that we can close the thermodynamic cycle for metal-imidazole complexes with up to six imidazole molecules in the first solvation shell. Given the success in closing the thermodynamic cycles, we then used the same extended sampling method for six other metal ions (Ag(I), Ca(II), Cd(II), Cu(I), Fe(II), and Mg(II)) to obtain new parameters. Since these new parameters can reproduce the one-imidazole geometry and energy accurately, we hypothesize that they will reasonably predict the binding free energy of higher-level coordination numbers. Hence, we did not extend the analysis of these ions up to six imidazole complexes. Overall, the results shed light on metal-protein interactions by emphasizing the importance of ligand-ligand interaction and metal-π-stacking within metalloproteins.

Electroosmotic Sensing of Uncharged Peptides and Differentiating Their Phosphorylated States Using Nanopores

The correct characterization and identification of different kinds of proteins is crucial for the survival and development of living organisms, and proteomics research promotes the analysis and understanding of future genome functions. Nanopore technique has been proved to accurately identify individual nucleotides. However, accurate and rapid protein sequencing is difficult due to the variability of protein structures that contains more than 20 amino acids, and it remains very challenging especially for uncharged peptides as they can not be electrophoretically driven through the nanopore. Graphene nanopores have the advantages of high accuracy, sensitivity and low cost in identifying protein phosphorylation modifications. Here, by using all-atom molecular dynamics simulations, charged graphene nanopores are employed to electroosmotically capture and sense uncharged peptides. By further mimicking AFM manipulation of single molecules, it is also found that the uncharged peptides and their phosphorylated states could also be differentiated by both the ionic current and pulling force signals during their pulling processes through the nanopore with a slow and constant velocity. The results shows ability of using nanopores to detect and discriminate single amino acid and its phosphorylation, which is essential for the future low-cost and high-throughput sequencing of protein residues and their post-translational modifications.

Cas12a domain flexibility guides R-loop formation and forces RuvC resetting

The specific nature of CRISPR-Cas12a makes it a desirable RNA-guided endonuclease for biotechnology and therapeutic applications. To understand how R-loop formation within the compact Cas12a enables target recognition and nuclease activation, we used cryo-electron microscopy to capture wild-type Acidaminococcus sp. Cas12a R-loop intermediates and DNA delivery into the RuvC active site. Stages of Cas12a R-loop formation-starting from a 5-bp seed-are marked by distinct REC domain arrangements. Dramatic domain flexibility limits contacts until nearly complete R-loop formation, when the non-target strand is pulled across the RuvC nuclease and coordinated domain docking promotes efficient cleavage. Next, substantial domain movements enable target strand repositioning into the RuvC active site. Between cleavage events, the RuvC lid conformationally resets to occlude the active site, requiring re-activation. These snapshots build a structural model depicting Cas12a DNA targeting that rationalizes observed specificity and highlights mechanistic comparisons to other class 2 effectors.

Reconstruction of zinc-metal battery solvation structures operating from -50 ~ +100 °C

Serious solvation effect of zinc ions has been considered as the cause of the severe side reactions (hydrogen evolution, passivation, dendrites, and etc.) of aqueous zinc metal batteries. Even though the regulation of cationic solvation structure has been widely studied, effects of the anionic solvation structures on the zinc metal were rarely examined. Herein, co-reconstruction of anionic and cationic solvation structures was realized through constructing a new multi-component electrolyte (Zn(BF4)2-glycerol-boric acid-chitosan-polyacrylamide, simplified as ZGBCP), which incorporates double crosslinking network via the esterification, protonation and polymerization reactions, thereby combining multiple advantages of 'liquid-like' high conductivity, 'gel-like' robust interface, and 'solid-like' high Zn2+ transfer number. Based on the ZGBCP electrolyte, the Zn anodes achieve record-low polarization and stable cycling. Furthermore, the ZGBCP electrolyte renders the AZMBs ultrawide working temperature (-50 °C ~ +100 °C) and ultralong cycle life (30000 cycles), which further validates the feasibility of the dual solvation structure strategy and provides a innovative perspective for the development of high-performance AZMBs.

Molecular Dynamics Simulation on the Suppression Mechanism of Phosphorylation to Ser222 by Allosteric Inhibitors Targeting MEK1/2 Kinase

Allosteric inhibitors of mitogen-activated protein kinase 1 (MEK1) reveal distinct interactions with MEK1 activation loop residues. The structural analyses will determine whether, and how, distinct inhibitors suppress the phosphorylation of MEK1 and may guide future therapeutic development. In this study, we explored the suppression mechanism of the phosphorylation process in the presence of MEK allosteric inhibitors, such as selumetinib, trametinib, cobimetinib, and CH5126766, by employing molecular dynamics simulations accompanied by principal component analysis. The simulations of wildtype MEK1 show that Ser222 can come close to γ-phosphate but not Ser218. We have found the conformation where Ser222 is within 5 Å of distance, which makes Ser222 accessible for γ-phosphate. The conformation analysis from the simulations of MEK1 in the presence of allosteric inhibitors reveals that the inhibitor restricts the flexibility of Ser222 through strong interactions with the activation loop, Lys97, and water mediates interactions with amino acids in the vicinity. The results reveal that all the inhibitors act as screeners between the activation loop and Mg-ATP and restricting the flexibility of the activation loop through strong interaction causes the suppression of the phosphorylation process of MEK1. The results conclude that a strong interaction of allosteric inhibitors with the activation loop restricts the movement of Ser222 toward Mg-ATP, which could be the dominant factor for the suppression of phosphorylation in MEK1. This research will provide novel insights to design effective anticancer therapeutics for targeting MEK1 in the future.

Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Molecular design of hydroxamic acid-based derivatives as urease inhibitors of Helicobacter pylori

Helicobacter pylori is the main causative agent of gastric cancer, especially non-cardiac gastric cancers. This bacterium relies on urease producing much ammonia to colonize the host. Herein, the study provides valuable insights into structural patterns driving urease inhibition for high-activity molecules designed via exploring known inhibitors. Firstly, an ensemble model was devised to predict the inhibitory activity of novel compounds in an automated workflow (R2 = 0.761) that combines four machine learning approaches. The dataset was characterized in terms of chemical space, including molecular scaffolds, clustering analysis, distribution for physicochemical properties, and activity cliffs. Through these analyses, the hydroxamic acid group and the benzene ring responsible for distinct activity were highlighted. Activity cliff pairs uncovered substituents of the benzene ring on hydroxamic acid derivatives are key structures for substantial activity enhancement. Moreover, 11 hydroxamic acid derivatives were designed, named mol1-11. Results of molecular dynamic simulations showed that the mol9 exhibited stabilization of the active site flap's closed conformation and are expected to be promising drug candidates for Helicobacter pylori infection and further in vitro, in vivo, and clinical trials to demonstrate in future.

Tailoring Electric Double Layer by Cation Specific Adsorption for High-Voltage Quasi-Solid-State Lithium Metal Batteries

The interfacial instability of high-nickel layered oxides severely plagues practical application of high-energy quasi-solid-state lithium metal batteries (LMBs). Herein, a uniform and highly oxidation-resistant polymer layer within inner Helmholtz plane is engineered by in situ polymerizing 1-vinyl-3-ethylimidazolium (VEIM) cations preferentially adsorbed on LiNi0.83Co0.11Mn0.06O2 (NCM83) surface, inducing the formation of anion-derived cathode electrolyte interphase with fast interfacial kinetics. Meanwhile, the copolymerization of [VEIM][BF4] and vinyl ethylene carbonate (VEC) endows P(VEC-IL) copolymer with the positively-charged imidazolium moieties, providing positive electric fields to facilitate Li+ transport and desolvation process. Consequently, the Li||NCM83 cells with a cut-off voltage up to 4.5 V exhibit excellent reversible capacity of 130 mAh g-1 after 1000 cycles at 25 °C and considerable discharge capacity of 134 mAh g-1 without capacity decay after 100 cycles at -20 °C. This work provides deep understanding on tailoring electric double layer by cation specific adsorption for high-voltage quasi-solid-state LMBs.

Elucidating the transport of water and ions in the nanochannel of covalent organic frameworks by molecular dynamics

Inspired by biological channels, achieving precise separation of ion/water and ion/ion requires finely tuned pore sizes at molecular dimensions and deliberate exposure of charged groups. Covalent organic frameworks (COFs), a class of porous crystalline materials, offer well-defined nanoscale pores and diverse structures, making them excellent candidates for nanofluidic channels that facilitate ion and water transport. In this study, we perform molecular simulations to investigate the structure and kinetics of water and ions confined within the typical COFs with varied exposure of charged groups. The COFs exhibit vertically arrayed nanochannels, enabling diffusion coefficients of water molecules within COFs to remain within the same order of magnitude as in the bulk. The motion of water molecules manifests in two distinct modes, creating a mobile hydration layer around acid groups. The ion diffusion within COFs displays a notable disparity between monovalent (M+) and divalent (M2+) cations. As a result, the selectivity of M+/M2+ can exceed 100, while differentiation among M+ is less pronounced. In addition, our simulations indicate a high rejection (R > 98%) in COFs, indicating their potential as ideal materials for desalination. The chemical flexibility of COFs indicates that would hold significant promise as candidates for advanced artificial ion channels and separation membranes.

Improving Aqueous Zinc Ion Batteries with Alkali Metal Ions

Aqueous zinc (Zn) ion batteries have received broad attention recently. However, their practical application is limited by severe Zn dendrite growth and the hydrogen evolution reaction. In this study, three alkali metal ions (Li+, Na+, and K+) are added in ZnSO4 electrolytes, which are subjected to electrochemical measurements and molecular dynamics simulations. The studies show that since K+ has the highest mobility and self-diffusion coefficient among the four ions (Li+, Na+, K+, and Zn2+), it enables K+ to preferentially approach a zinc dendrite at an earlier time, driven by a negative electric field during a cathodic process. The electric double layer, with K+ around the negatively charged Zn dendrite, inhibits dendrite growth and mitigates the hydrogen evolution reaction on the Zn anode. Under this kinetic effect, the Zn-Zn symmetric cell with K+ exhibits a long cycling stability of 1000 h at 1 mA·cm-2 of 1 mAh·cm-2 and 190 h at 30 mA·cm-2 of 2 mAh·cm-2. Such a kinetic effect is also observed with additives Na+ and Li+, though less profound than that of K+.

Inactivation of polyphenol oxidase by low intensity DC field: Experiment and mechanism analysis via molecular dynamics simulation and molecular docking

In this study, inactivation of mushroom polyphenol oxidase (PPO) by low intensity direct current (DC) electric field and its molecular mechanism were investigated. In the experiments under 3 V/cm, 5 V/cm, 7 V/cm and 9 V/cm electric fields, PPOs were all completely inactivated after different exposure times. Under 1 V/cm, a residual activity of 11.88 % remained. The inactivation kinetics confirms to Weibull model. Under 1-7 V/cm, n value closes to a constant about 1.3. The structural analysis of PPO under 3 V/cm and 5 V/cm by fluorescence emission spectroscopy and molecular dynamics (MD) simulation showed that the tertiary structure was slightly changed with increased radius of gyration, higher potential energy and rate of C-alpha fluctuation. After exposure to the electric field, most of the hydrophobic tryptophan (TRP) residues turned to the hydrophilic surface, resulting the fluorescence red-shifted and quenched. Molecular docking indicated that the receptor binding domain of catechol in PPO was changed. PPO under electric field was MD simulated the first time, revealing the changing mechanism of the electric field itself on PPO, a binuclear copper enzyme, which has a metallic center. All these suggest that the low intensity DC electric field would be a promising option for enzymatic browning inhibition or even enzyme activity inactivation.

Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure. The infrastructure provides access to QM/MM models that may serve as the foundation for QM/MM-ΔMLP potentials, which supplement the semiempirical QM/MM model with a MLP correction trained to reproduce ab initio QM/MM energies and forces. Efficient optimization of minimum free energy pathways is enabled through a new surface-accelerated finite-temperature string method implemented in the FE-ToolKit package. Furthermore, we interfaced Sander with the i-PI software by implementing the socket communication protocol used in the i-PI client-server model. The new interface with i-PI allows for the treatment of nuclear quantum effects with semiempirical QM/MM-ΔMLP models. The modular interoperable software is demonstrated on proton transfer reactions in guanine-thymine mispairs in a B-form deoxyribonucleic acid helix. The current work represents a considerable advance in the development of modular software for performing free energy simulations of chemical reactions that are important in a wide range of applications.

Cation-Binding of Glutamate in Aqueous Solution

Interactions of the cations Li+, Na+, Mg2+, and Ca2+ with L-glutamate (Glu-) in aqueous solution were studied at room temperature with dielectric relaxation spectroscopy in the gigahertz region. Spectra of ∼0.4 M NaGlu with added LiCl, NaCl, MgCl2, or CaCl2 (c(MCln) ≤ 1.5 M) were evaluated and experiments supplemented by density functional theory and 3D reference interaction site model (3D-RISM) calculations. In addition to the modes found for aqueous NaGlu, namely, the reorientation of free Glu- ions (peaking at ∼1.6 GHz), of moderately retarded H2O molecules hydrating the carboxylate moieties of Glu- (∼8.4 GHz), of the cooperative resettling of the H-bond network of bulk water (∼20 GHz), and its preceding fast H-bond flip (∼400 GHz), an additional low-frequency relaxation at ∼0.4 GHz was detected upon the addition of the four salts. In the case of NaGlu + MgCl2(aq) and NaGlu + CaCl2(aq), this mode could be unequivocally assigned to an ion pair formed by the cation and the side-chain carboxylate moiety of Glu-. For NaGlu + LiCl(aq), either this species or a backbone-[Li+-H2O-Cl--Glu-] triple ion is formed. Binding constants increase in the order Li+ Collapse

Key Words Collapse

MESH Headings Collapse

Grants

• Fonds der Chemischen Industrie
• G.A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences

Collapse

Affiliation(s)

Sergej Friesen Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany
Sergey E. Kruchinin G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya st. 1, Ivanovo 153045, Russian Federation
Marina V. Fedotova G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya st. 1, Ivanovo 153045, Russian Federation
Richard Buchner Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany

Collapse

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Dynamic molecular portraits of ion-conducting pores characterize functional states of TRPV channels

Structural biology is solving an ever-increasing number of snapshots of ion channel conformational ensembles. Deciphering ion channel mechanisms, however, requires understanding the ensemble dynamics beyond the static structures. Here, we present a molecular modeling-based approach characterizing the ion channel structural intermediates, or their "dynamic molecular portraits", by assessing water and ion conductivity along with the detailed evaluation of pore hydrophobicity and residue packing. We illustrate the power of this approach by analyzing structures of few vanilloid-subfamily transient receptor potential (TRPV) channels. Based on the pore architecture, there are three major states that are common for TRPVs, which we call α-closed, π-closed, and π-open. We show that the pore hydrophobicity and residue packing for the open state is most favorable for the pore conductance. On the contrary, the α-closed state is the most hydrophobic and always non-conducting. Our approach can also be used for structural and functional classification of ion channels.

Binding of Phosphate Species to Ca2+ and Mg2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4-, HPO42-, PO43- coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal-phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation-anion interactions.

Secondary Conformational Checkpoint in CRISPR-Cas9

A specific checkpoint between target DNA binding and cleavage primarily governs the precision of Cas9 gene editing. Although various CRISPR-Cas9 variants have been developed to improve DNA cleavage accuracy, we still lack a comprehensive understanding of how they work at the molecular level. Herein, we have focused on studying the late-stage conformational transitions of Cas9 and an evolved Cas9 mutant (evoCas9) that start from the precleavage state. Our submilliseconds of dynamic simulations reveal that the presence of base mismatches leads the HNH nuclease domain of Cas9 to alter its principal functional modes of motion, thereby impairing its conformational activation. This observation suggests the existence of a secondary conformational checkpoint that fine-tunes the final DNA cleavage activation. Remarkably, evoCas9 is prone to deviating from the normal activation pathway with base mismatches. This is characterized by a noticeable shift in the positioning of the HNH domain and a significantly perturbed allosteric communication network within the enzyme. Therefore, the mutations evolved in evoCas9 also reinforce the secondary checkpoint in addition to the previously identified primary checkpoint, collectively ensuring this variant's high gene-editing accuracy. This mechanism should also apply to other Cas9-guide RNA variants with enhanced fidelity.

Quantum sensing of microRNAs with nitrogen-vacancy centers in diamond

Label-free detection of nucleic acids such as microRNAs holds great potential for early diagnostics of various types of cancers. Measuring intrinsic biomolecular charge using methods based on field effect has been a promising way to accomplish label-free detection. However, the charges of biomolecules are screened by counter ions in solutions over a short distance (Debye length), thereby limiting the sensitivity of these methods. Here, we measure the intrinsic magnetic noise of paramagnetic counter ions, such as Mn2+, interacting with microRNAs using nitrogen-vacancy (NV) centers in diamond. All-atom molecular dynamics simulations show that microRNA interacts with the diamond surface resulting in excess accumulation of Mn ions and stronger magnetic noise. We confirm this prediction by observing an increase in spin relaxation contrast of the NV centers, indicating higher Mn2+ local concentration. This opens new possibilities for next-generation quantum sensing of charged biomolecules, overcoming limitations due to the Debye screening.

Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe3+ Ion

The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.

RNA targeting and cleavage by the type III-Dv CRISPR effector complex

CRISPR-Cas are adaptive immune systems in bacteria and archaea that utilize CRISPR RNA-guided surveillance complexes to target complementary RNA or DNA for destruction1-5. Target RNA cleavage at regular intervals is characteristic of type III effector complexes6-8. Here, we determine the structures of the Synechocystis type III-Dv complex, an apparent evolutionary intermediate from multi-protein to single-protein type III effectors9,10, in pre- and post-cleavage states. The structures show how multi-subunit fusion proteins in the effector are tethered together in an unusual arrangement to assemble into an active and programmable RNA endonuclease and how the effector utilizes a distinct mechanism for target RNA seeding from other type III effectors. Using structural, biochemical, and quantum/classical molecular dynamics simulation, we study the structure and dynamics of the three catalytic sites, where a 2'-OH of the ribose on the target RNA acts as a nucleophile for in line self-cleavage of the upstream scissile phosphate. Strikingly, the arrangement at the catalytic residues of most type III complexes resembles the active site of ribozymes, including the hammerhead, pistol, and Varkud satellite ribozymes. Our work provides detailed molecular insight into the mechanisms of RNA targeting and cleavage by an important intermediate in the evolution of type III effector complexes.

An aberrant fused in sarcoma liquid droplet of amyotrophic lateral sclerosis pathological variant, R495X, accelerates liquid-solid phase transition

Intracellular aggregation of fused in sarcoma (FUS) is associated with the pathogenesis of familial amyotrophic lateral sclerosis (ALS). Under stress, FUS forms liquid droplets via liquid-liquid phase separation (LLPS). Two types of wild-type FUS LLPS exist in equilibrium: low-pressure LLPS (LP-LLPS) and high-pressure LLPS (HP-LLPS); the former dominates below 2 kbar and the latter over 2 kbar. Although several disease-type FUS variants have been identified, the molecular mechanism underlying accelerated cytoplasmic granule formation in ALS patients remains poorly understood. Herein, we report the reversible formation of the two LLPS states and the irreversible liquid-solid transition, namely droplet aging, of the ALS patient-type FUS variant R495X using fluorescence microscopy and ultraviolet-visible absorption spectroscopy combined with perturbations in pressure and temperature. Liquid-to-solid phase transition was accelerated in the HP-LLPS of R495X than in the wild-type variant; arginine slowed the aging of droplets at atmospheric conditions by inhibiting the formation of HP-LLPS more selectively compared to that of LP-LLPS. Our findings provide new insight into the mechanism by which R495X readily forms cytoplasmic aggregates. Targeting the aberrantly formed liquid droplets (the HP-LLPS state) of proteins with minimal impact on physiological functions could be a novel therapeutic strategy for LLPS-mediated protein diseases.

Exploring Product Release from Yeast Cytosine Deaminase with Metadynamics

The yeast cytosine deaminase (yCD) enzyme/5-fluorocytosine prodrug system is a promising candidate for targeted chemotherapeutics. After conversion of the prodrug into the toxic chemotherapeutic drug, 5-fluorouracil (5-FU), the slow product release from the enzyme limits the overall catalytic efficiency of the enzyme/prodrug system. Here, we present a computational study of the product release of the anticancer drug, 5-FU, from yCD using metadynamics. We present a comparison of the 5-FU drug to the natural enzyme product, uracil. We use volume-based metadynamics to compute the free energy landscape for product release and show a modest binding affinity for the product to the enzyme, consistent with experiments. Next, we use infrequent metadynamics to estimate the unbiased release rate from Kramers time-dependent rate theory and find a favorable comparison to experiment with a slower rate of product release for the 5-FU system. Our work demonstrates how adaptive sampling methods can be used to study the protein-ligand unbinding process for engineering enzyme/prodrug systems and gives insights into the molecular mechanism of product release for the yCD/5-FU system.

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

Anti-leishmanial activity of Eleutherine plicata Herb. and predictions of isoeleutherin and its analogues

Introduction: Leishmaniasis is caused by protozoa of the genus Leishmania, classified as tegumentary and visceral. The disease treatment is still a serious problem, due to the toxic effects of available drugs, the costly treatment and reports of parasitic resistance, making the search for therapeutic alternatives urgent. This study assessed the in vitro anti-leishmanial potential of the extract, fractions, and isoeleutherin from Eleutherine plicata, as well as the in silico interactions of isoeleutherin and its analogs with Trypanothione Reductase (TR), in addition to predicting pharmacokinetic parameters. Methods: From the ethanolic extract of E. plicata (EEEp) the dichloromethane fraction (FDEp) was obtained, and isoeleutherin isolated. All samples were tested against promastigotes, and parasite viability was evaluated. Isoeleutherin analogues were selected based on similarity in databases (ZINK and eMolecules) to verify the impact on structural change. Results and Discussion: The extract and its fractions were not active against the promastigote form (IC50 > 200 μg/mL), while isoeleutherin was active (IC50 = 25 μg/mL). All analogues have high intestinal absorption (HIA), cell permeability was moderate in Caco2 and low to moderate in MDCK. Structural changes interfered with plasma protein binding and blood-brain barrier permeability. Regarding metabolism, all molecules appear to be CYP3A4 metabolized and inhibited 2-3 CYPs. Molecular docking and molecular dynamics assessed the interactions between the most stable configurations of isoeleutherin, analogue compound 17, and quinacrine (control drug). Molecular dynamics simulations demonstrated stability and favorable interactions with TR. In summary, fractionation contributed to antileishmanial activity and isoleutherin seems to be promising. Structural alterations did not contribute to improve pharmacokinetic aspects and analogue 17 proved to be more promising than isoeleutherin, presenting better stabilization in TR.

Critical Solvation Structures Arrested Active Molecules for Reversible Zn Electrochemistry

Aqueous Zn-ion batteries (AZIBs) have attracted increasing attention in next-generation energy storage systems due to their high safety and economic. Unfortunately, the side reactions, dendrites and hydrogen evolution effects at the zinc anode interface in aqueous electrolytes seriously hinder the application of aqueous zinc-ion batteries. Here, we report a critical solvation strategy to achieve reversible zinc electrochemistry by introducing a small polar molecule acetonitrile to form a "catcher" to arrest active molecules (bound water molecules). The stable solvation structure of [Zn(H2O)6]2+ is capable of maintaining and completely inhibiting free water molecules. When [Zn(H2O)6]2+ is partially desolvated in the Helmholtz outer layer, the separated active molecules will be arrested by the "catcher" formed by the strong hydrogen bond N-H bond, ensuring the stable desolvation of Zn2+. The Zn||Zn symmetric battery can stably cycle for 2250 h at 1 mAh cm-2, Zn||V6O13 full battery achieved a capacity retention rate of 99.2% after 10,000 cycles at 10 A g-1. This paper proposes a novel critical solvation strategy that paves the route for the construction of high-performance AZIBs.

Protein S-palmitoylation is markedly inhibited by 4″-alkyl ether lipophilic derivatives of EGCG, the major green tea polyphenol: In vitro and in silico studies

S-palmitoylation is a dynamic lipid-based protein post-translational modification facilitated by a family of protein acyltransferases (PATs) commonly known as DHHC-PATs or DHHCs. It is the only lipid modification that is reversible, and this very fact uniquely qualifies it for therapeutic interventions through the development of DHHC inhibitors. Herein, we report that 4″-alkyl ether lipophilic derivatives of EGCG can effectively inhibit protein S-palmitoylation in vitro. With the help of metabolic labeling followed by copper(I)-catalyzed azide-alkyne cycloaddition Click reaction, we demonstrate that 4″-C14 EGCG and 4″-C16 EGCG markedly inhibited S-palmitoylation in various mammalian cells including HEK 293T, HeLa, and MCF-7 using both in gel fluorescence as well as confocal microscopy. Further, these EGCG derivatives were able to attenuate the S-palmitoylation to the basal level in DHHC3-overexpressed cells, suggesting that they are plausibly targeting DHHCs. Confocal microscopy data qualitatively reflected spatial and temporal distribution of S-palmitoylated proteins in different sub-cellular compartments and the inhibitory effects of 4″-C14 EGCG and 4″-C16 EGCG were clearly observed in the native cellular environment. Our findings were further substantiated by in silico analysis which revealed promising binding affinity and interactions of 4″-C14 EGCG and 4″-C16 EGCG with key amino acid residues present in the hydrophobic cleft of the DHHC20 enzyme. We also demonstrated the successful inhibition of S-palmitoylation of GAPDH by 4″-C16 EGCG. Taken together, our in vitro and in silico data strongly suggest that 4″-C14 EGCG and 4″-C16 EGCG can act as potent inhibitors for S-palmitoylation and can be employed as a complementary tool to investigate S-palmitoylation.

The Deinococcus protease PprI senses DNA damage by directly interacting with single-stranded DNA

Bacteria have evolved various response systems to adapt to environmental stress. A protease-based derepression mechanism in response to DNA damage was characterized in Deinococcus, which is controlled by the specific cleavage of repressor DdrO by metallopeptidase PprI (also called IrrE). Despite the efforts to document the biochemical, physiological, and downstream regulation of PprI-DdrO, the upstream regulatory signal activating this system remains unclear. Here, we show that single-stranded DNA physically interacts with PprI protease, which enhances the PprI-DdrO interactions as well as the DdrO cleavage in a length-dependent manner both in vivo and in vitro. Structures of PprI, in its apo and complexed forms with single-stranded DNA, reveal two DNA-binding interfaces shaping the cleavage site. Moreover, we show that the dynamic monomer-dimer equilibrium of PprI is also important for its cleavage activity. Our data provide evidence that single-stranded DNA could serve as the signal for DNA damage sensing in the metalloprotease/repressor system in bacteria. These results also shed light on the survival and acquired drug resistance of certain bacteria under antimicrobial stress through a SOS-independent pathway.

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

An alpha-helical lid guides the target DNA toward catalysis in CRISPR-Cas12a

CRISPR-Cas12a is a powerful RNA-guided genome-editing system that generates double-strand DNA breaks using its single RuvC nuclease domain by a sequential mechanism in which initial cleavage of the non-target strand is followed by target strand cleavage. How the spatially distant DNA target strand traverses toward the RuvC catalytic core is presently not understood. Here, continuous tens of microsecond-long molecular dynamics and free-energy simulations reveal that an α-helical lid, located within the RuvC domain, plays a pivotal role in the traversal of the DNA target strand by anchoring the crRNA:target strand duplex and guiding the target strand toward the RuvC core, as also corroborated by DNA cleavage experiments. In this mechanism, the REC2 domain pushes the crRNA:target strand duplex toward the core of the enzyme, while the Nuc domain aids the bending and accommodation of the target strand within the RuvC core by bending inward. Understanding of this critical process underlying Cas12a activity will enrich fundamental knowledge and facilitate further engineering strategies for genome editing.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Unraveling potential mechanism of different metal ions effect on anammox through big data analysis, molecular docking and molecular dynamics simulation

Heavy metals (HMs) have been widely reported to pose an adverse effect on anaerobic ammonia oxidation (anammox) bacteria, yet the underlying mechanisms remain unclear. This study provides new insights into the potential mechanisms of interaction between HMs and functional enzymes through big date analysis, molecular docking and molecular dynamics simulation. The statistical analysis indicated that 10 mg/L Cu(II) and Cd(II) reduced nitrogen removal rate (NRR) by 85% and 43%, while 5 mg/L Fe(II) enhanced NRR by 29%. Additionally, the results of molecular simulations provided a microscopic interpretation for these macroscopic data. Molecular docking revealed that Hg(II) formed a distinctive binding site on ferritin, while other HMs resided at iron oxidation sites. Furthermore, HMs exhibited distinct binding sites on hydrazine dehydrogenase. Concurrently, the molecular dynamics simulation results further substantiated their capacity to form complexes. Cu(II) displayed the strongest binding affinity with ferritin for -1576 ± 79 kJ/mol in binding free energy calculation. Moreover, Cd(II) bound to ferritin and HDH for -1052.67 ± 58.49 kJ/mol, -290.02 ± 49.68 kJ/mol, respectively. This research addressed a crucial knowledge gap, shedding light on potential applications for remediating heavy metal-laden industrial wastewater.

Development of a Fluctuating Charge Model for Zinc-Containing Metalloproteins

Metalloproteins widely exist in biology and play important roles in various processes. To accurately simulate metalloprotein systems, modeling polarization and charge transfer effects is vital. The fluctuating charge (FQ) model can efficiently generate atomic charges and simulate the charge transfer effect; it has been developed for a wide range of applications, but few models have been specifically tailored for metalloproteins. In this study, we present a fluctuating charge model specifically for zinc-containing metalloproteins based on the extended charge equilibration (EQeq) scheme. Our model was parametrized to reproduce CM5 charges instead of RESP/CHELPG charges because the former is less dependent on the conformation or basis set, does not suffer from unphysical charges for buried atoms, and is still being able to well reproduce the molecular dipoles. During our study, we found that adding the Pauling-bond-order-like term (referred to as the "+C term" in a previous study) between the zinc ion and ligating atoms significantly improves the model's performance. Although our model was trained for four-coordinated zinc sites only, our results indicated it can well describe the atomic charges in diverse zinc sites. Morever, our model was used to generate partial charges for the metal sites in three different zinc-containing metalloproteins (with four-, five-, and six-coordinated metal sites, respectively). These charges exhibited performance comparable to that of the RESP charges in molecular dynamics (MD) simulations. Additional tests indicated our model could also well reproduce the CM5 charges when geometric changes were involved. Those results indicate that our model can efficiently calculate the atomic charges for metal sites and well simulate the charge transfer effect, which marks an important step toward developing versatile polarizable force fields for metalloproteins.

Reconciling ASPP-p53 binding mode discrepancies through an ensemble binding framework that bridges crystallography and NMR data

ASPP2 and iASPP bind to p53 through their conserved ANK-SH3 domains to respectively promote and inhibit p53-dependent cell apoptosis. While crystallography has indicated that these two proteins employ distinct surfaces of their ANK-SH3 domains to bind to p53, solution NMR data has suggested similar surfaces. In this study, we employed multi-scale molecular dynamics (MD) simulations combined with free energy calculations to reconcile the discrepancy in the binding modes. We demonstrated that the binding mode based solely on a single crystal structure does not enable iASPP's RT loop to engage with p53's C-terminal linker-a verified interaction. Instead, an ensemble of simulated iASPP-p53 complexes facilitates this interaction. We showed that the ensemble-average inter-protein contacting residues and NMR-detected interfacial residues qualitatively overlap on ASPP proteins, and the ensemble-average binding free energies better match experimental KD values compared to single crystallgarphy-determined binding mode. For iASPP, the sampled ensemble complexes can be grouped into two classes, resembling the binding modes determined by crystallography and solution NMR. We thus propose that crystal packing shifts the equilibrium of binding modes towards the crystallography-determined one. Lastly, we showed that the ensemble binding complexes are sensitive to p53's intrinsically disordered regions (IDRs), attesting to experimental observations that these IDRs contribute to biological functions. Our results provide a dynamic and ensemble perspective for scrutinizing these important cancer-related protein-protein interactions (PPIs).

Characterization of Various Noncovalent Polyphenol-Starch Complexes and Their Prebiotic Activities during In Vitro Digestion and Fermentation

This study explores the structural characterization of six noncovalent polyphenol-starch complexes and their prebiotic activities during in vitro digestion and fermentation. Ferulic acid, caffeic acid, gallic acid, isoquercetin, astragalin, and hyperin were complexed with sweet potato starch (SPS). The polyphenols exhibited high binding capacity (>70%) with SPS. A partial release of flavonoids from the complexes was observed via in vitro digestion, while the phenolic acids remained tightly bound. Molecular dynamics (MD) simulation revealed that polyphenols altered the spatial configuration of polysaccharides and intramolecular hydrogen bonds formed. Additionally, polyphenol-SPS complexes exerted inhibitory effects on starch digestion compared to gelatinized SPS, owing to the increase in resistant starch fraction. It revealed that the different complexes stimulated the growth of Lactobacillus rhamnosus and Bifidobacterium bifidum, while inhibiting the growth of Escherichia coli. Moreover, in vitro fermentation experiments revealed that complexes were utilized by the gut microbiota, resulting in the production of short-chain fatty acids and a decrease in pH. In addition, the polyphenol-SPS complexes altered the composition of gut microbiota by promoting the growth of beneficial bacteria and decreasing pathogenic bacteria. Polyphenol-SPS complexes exhibit great potential for use as a prebiotic and exert dual beneficial effects on gut microbiota.

Solvation energies of the ferrous ion in water and in ammonia at various temperatures

CONTEXT

The solvation of metal ions is crucial to understanding relevant properties in physics, chemistry, or biology. Therefore, we present solvation enthalpies and solvation free energies of the ferrous ion in water and ammonia. Our results agree well with the experimental reports for the hydration free energy and hydration enthalpy. We obtained [Formula: see text] kJ mol[Formula: see text] for the hydration free energy and [Formula: see text] kJ mol[Formula: see text] for the hydration enthalpy of ferrous ion in water at room temperature. At ambient temperature, we obtained [Formula: see text] kJ mol[Formula: see text] as the [Formula: see text] ammoniation free energy and [Formula: see text] kJ mol[Formula: see text] for the ammoniation enthalpy. In addition, the free energy of solvation is deeply affected when the temperature increases. This pattern can be attributed to the rise of entropy when the temperature rises. Besides, the temperature does not affect the ammoniation enthalpies and the hydration enthalpy of the [Formula: see text] ion.

METHOD

All the geometry optimizations are performed at the MP2 methods associated with the 6-31++g(d,p) basis set of Pople. solvated phase structures of [Formula: see text] ion in water or in ammonia are performed using the PCM model. The [Formula: see text] program suite was used to perform all the calculations. The program TEMPO was also used to evaluate the temperature sensitivity of the different obtained geometries.

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.

Role of ATP Hydrolysis and Product Release in the Translocation Mechanism of SARS-CoV-2 NSP13

In response to the emergence of COVID-19, caused by SARS-CoV-2, there has been a growing interest in understanding the functional mechanisms of the viral proteins to aid in the development of new therapeutics. Nonstructural protein 13 (nsp13) helicase is an attractive target for antivirals because it is essential for viral replication and has a low mutation rate, yet the structural mechanisms by which this enzyme binds and hydrolyzes ATP to cause unidirectional RNA translocation remain elusive. Using Gaussian accelerated molecular dynamics (GaMD), we generated comprehensive conformational ensembles of all substrate states along the ATP-dependent cycle. Shape-GMM clustering of the protein yields four protein conformations that describe an opening and closing of both the ATP pocket and the RNA cleft that is achieved through a combination of conformational selection and induction along the ATP hydrolysis cycle. Furthermore, three protein-RNA conformations are observed that implicate motifs Ia, IV, and V as playing a pivotal role in an ATP-dependent inchworm translocation mechanism. Finally, based on a linear discriminant analysis of protein conformations, we identify L405 as a pivotal residue for the opening and closing mechanism and propose a L405D mutation as a way to disrupt translocation. This research enhances our understanding of nsp13's role in viral replication and could contribute to the development of antiviral strategies.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Structural and dynamical changes of the Streptococcus gordonii metalloregulatory ScaR protein induced by Mn2+ ion binding

Divalent metal ions are essential micronutrients for many intercellular reactions. Maintaining their homeostasis is necessary for the survival of bacteria. In Streptococcus gordonii, one of the primary colonizers of the tooth surface, the cellular concentration of manganese ions (Mn2+) is regulated by the manganese-sensing transcriptional factor ScaR which controls the expression of proteins involved in manganese homeostasis. To resolve the molecular mechanism through which the binding of Mn2+ ions increases the binding affinity of ScaR to DNA, a variety of computational (QM and MD) and experimental (ITC, DSC, EMSA, EPR, and CD) methods were applied. The computational results showed that Mn2+ binding induces a conformational change in ScaR that primarily affects the position of the DNA binding domains and, consequently, the DNA binding affinity of the protein. In addition, experimental results revealed a 1:4 binding stoichiometry between ScaR dimer and Mn2+ ions, while the computational results showed that the binding of Mn2+ ions in the primary binding sites is sufficient to induce the observed conformational change of ScaR.

Artificial Sodium Channels for Enhanced Osmotic Energy Harvesting

Inspired by the ion channels of electric eels, we can use biomimetic nanofluidic materials to harvest the osmotic power released by mixing seawater and river water. While biological ion channels have both cation/anion and inter-cation selectivity, previous nanofluidic materials neglected the latter. As a result, NaCl solutions were generally used to simulate river water, ignoring the fact that the dominating cation in river water is typically Ca2+. In this work, we show that the different ionic compositions of seawater and river water can be exploited to improve osmotic power density by employing biomimetic sodium selective materials. Inspired by a range of properties of biological sodium channels, we constructed artificial sodium channels with zeolitic imidazolate framework-65 crystals, which selectively transport Na+ but almost completely block Ca2+. Resultantly, the effective concentration gradient of seawater/river water is dramatically increased by preventing the major cations in the river water from participating in the ion diffusion. As a result, the osmotic power density can be increased by more than 1 order of magnitude. These results should open new avenues to develop high-performance osmotic generators and may advance other applications based on biomimetic ion channels such as neuromorphic information processing.

Monovalent metal ion binding promotes the first transesterification reaction in the spliceosome

Cleavage and formation of phosphodiester bonds in nucleic acids is accomplished by large cellular machineries composed of both protein and RNA. Long thought to rely on a two-metal-ion mechanism for catalysis, structure comparisons revealed many contain highly spatially conserved second-shell monovalent cations, whose precise function remains elusive. A recent high-resolution structure of the spliceosome, essential for pre-mRNA splicing in eukaryotes, revealed a potassium ion in the active site. Here, we employ biased quantum mechanics/ molecular mechanics molecular dynamics to elucidate the function of this monovalent ion in splicing. We discover that the K+ ion regulates the kinetics and thermodynamics of the first splicing step by rigidifying the active site and stabilizing the substrate in the pre- and post-catalytic state via formation of key hydrogen bonds. Our work supports a direct role for the K+ ion during catalysis and provides a mechanistic hypothesis likely shared by other nucleic acid processing enzymes.

Hydrated Eutectic Electrolyte Induced Bilayer Interphase for High-Performance Aqueous Zn-Ion Batteries with 100 °C Wide-Temperature Range

The practical implementation of aqueous zinc-ion batteries (AZIBs) encounters challenges such as dendrite growth, parasitic reactions, and severe decay in battery performance under harsh environments. Here, a novel hydrated eutectic electrolyte (HEE) composed of Zn(ClO4 )2 ·6H2 O, ethylene glycol (EG), and InCl3 solution is introduced to effectively extend the lifespan of AZIBs over a wide temperature range from -50 to 50 °C. Molecular dynamics simulations and spectroscopy analysis demonstrate that the H2 O molecules are confined within the liquid eutectic network through dual-interaction, involving coordination with Zn2+ and hydrogen bonding with EG, thus weakening the activity of free water and extending the electrochemical window. Importantly, cryo-transmission electron microscopy and spectroscopy techniques reveal that HEE in situ forms a zincophobic/zincophilic bilayer interphase by the dissociation-reduction of eutectic molecules. Specifically, the zincophilic interphase reduces the energy barrier for Zn nucleation, promoting uniform Zn deposition, while the zincophobic interphase prevents active water from contacting the Zn surface, thus inhibiting the side reactions. Furthermore, the relationships between the structural evolution of the liquid eutectic network and interfacial chemistry at electrode/electrolyte interphase are further discussed in this work. The scalability of this design strategy can bring benefits to AZIBs operating over a wide temperature range.