Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions

Unraveling the Structure-Spectrum Relationship of Yeast Phenylalanine Transfer RNA: Insights from Theoretical Modeling of Infrared Spectroscopy

Yeast phenylalanine tRNA (tRNAphe) is a paradigmatic model in structural biology. In this work, we combine molecular dynamics simulations and spectroscopy modeling to establish a direct link between its structure, conformational dynamics, and infrared (IR) spectra. Employing recently developed vibrational frequency maps and coupling models, we apply a mixed quantum/classical treatment of the line shape theory to simulate the IR spectra of tRNAphe in the 1600-1800 cm-1 region across its folded and unfolded conformations and under varying concentrations of Mg2+ ions. The predicted IR spectra of folded and unfolded tRNAphe are in good agreement with experimental measurements, validating our theoretical framework. We then elucidate how the characteristic L-shaped tertiary structure of the tRNA and its modulation in response to diverse chemical environments give rise to distinct IR absorption peaks and line shapes. These calculations effectively bridge IR spectroscopy experiments and atomistic molecular simulations, unraveling the molecular origins of the observed IR spectra of tRNAphe. This work presents a robust theoretical protocol for modeling the IR spectroscopy of nucleic acids, which will facilitate its application as a sensitive probe for detecting the fluctuating secondary and tertiary structures of these essential biological macromolecules.

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.

BAT2: an Open-Source Tool for Flexible, Automated, and Low Cost Absolute Binding Free Energy Calculations

Absolute binding free energy (ABFE) calculations with all-atom molecular dynamics (MD) have the potential to greatly reduce costs in the first stages of drug discovery. Here, we introduce BAT2, the new version of the Binding Affinity Tool (BAT.py), designed to combine full automation of ABFE calculations with high-performance MD simulations, making it a potential tool for virtual screening. We describe and test several changes and new features that were incorporated into the code, such as relative restraints between the protein and the ligand instead of using fixed dummy atoms, support for the OpenMM simulation engine, a merged approach to the application/release of restraints, support for cobinders and proteins with multiple chains, and many others. We also reduced the simulation times for each ABFE calculation, assessing the effect on the expected robustness and accuracy of the calculations.

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.

Unveiling the Mechanism of Phenamacril Resistance in F. graminearum: Computational and Experimental Insights into the C423A Mutation in FgMyoI

Phenamacril (PHA) is a highly selective fungicide for controlling fusarium head blight (FHB) mainly caused by F. graminearum and F. asiaticum. However, the C423A mutation in myosin I of F. graminearum (FgMyoI) leads to natural resistance to PHA. Here, based on the computational approaches and biochemical validation, we elucidate the atomic-level mechanism behind the natural resistance of F. graminearum to the fungicide PHA due to the C423A mutation in FgMyoI. The mutation leads to a rearrangement of pocket residues, resulting in increased size and flexibility of the binding pocket, which impairs the stable binding of PHA. MST experiments confirm that the mutant protein FgMyoIC423A exhibits significantly reduced affinity for PHA compared to wild-type FgMyoI and the nonresistant C423K mutant. This decreased binding affinity likely underlies the development of PHA resistance in F. graminearum. Conversely, the nonresistant C423K mutant retains sensitivity to PHA due to the introduction of a strong hydrogen bond donor, which facilitates stable binding of PHA in the pocket. These findings shed light on the molecular basis of PHA resistance and provide new directions for the creation of new myosin inhibitors.

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.

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.

Molecular Insights into the Variability in Infection and Immune Evasion Capabilities of SARS-CoV-2 Variants: A Sequence and Structural Investigation of the RBD Domain

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continuously emerge, an increasing number of mutations are accumulating in the Spike protein receptor-binding domain (RBD) region. Through sequence analysis of various Variants of Concern (VOC), we identified that they predominantly fall within the ο lineage although recent variants introduce any novel mutations in the RBD. Molecular dynamics simulations were employed to compute the binding free energy of these variants with human Angiotensin-converting enzyme 2 (ACE2). Structurally, the binding interface of the ο RBD displays a strong positive charge, complementing the negatively charged binding interface of ACE2, resulting in a significant enhancement in the electrostatic potential energy for the ο variants. Although the increased potential energy is partially offset by the rise in polar solvation free energy, enhanced electrostatic potential contributes to the long-range recognition between the ο variant's RBD and ACE2. We also conducted simulations of glycosylated ACE2-RBD proteins. The newly emerged ο (JN.1) variant has added a glycosylation site at N-354@RBD, which significantly weakened its binding affinity with ACE2. Further, our interaction studies with three monoclonal antibodies across multiple SARS-CoV-2 strains revealed a diminished neutralization efficacy against the ο variants, primarily attributed to the electrostatic repulsion between the antibodies and RBD interface. Considering the characteristics of the ο variant and the trajectory of emerging strains, we propose that newly developed antibodies against SARS-CoV-2 RBD should have surfaces rich in negative potential and, postbinding, exhibit strong van der Waals interactions. These findings provide invaluable guidance for the formulation of future therapeutic strategies.

Reversible pH-Gated MXene Membranes with Ultrahigh Mono-/Divalent-Ion Selectivity

Artificial ion channel membranes hold high promise in water treatment, nanofluidics, and energy conversion, but it remains a great challenge to construct such smart membranes with both reversible ion-gating capability and desirable ion selectivity. Herein, we constructed a smart MXene-based membrane via p-phenylenediamine functionalization (MLM-PPD) with highly stable and aligned two-dimensional subnanochannels, which exhibits reversible ion-gating capability and ultrahigh metal ion selectivity similar to biological ion channels. The pH-sensitive groups within the MLM-PPD channel confers excellent reversible Mg2+-gating capability with a pH-switching ratio of up to 100. The mono/divalent metal-ion selectivity up to 1243.8 and 400.9 for K+/Mg2+ and Li+/Mg2+, respectively, outperforms other reported membranes. Theoretical calculations combined with experimental results reveal that the steric hindrance and stronger PPD-ion interactions substantially enhance the energy barrier for divalent metal ions passing through the MLM-PPD, and thus leading to ultrahigh mono/divalent metal-ion selectivity. This work provides a new strategy for developing artificial-ion channel membranes with both reversible ion-gating functionality and high-ion selectivity for various applications.

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods

RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.

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.

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.

Thermodynamics of Metal-Acetate Interactions

Metal ions play crucial roles in protein- and ligand-mediated interactions. They not only act as catalysts to facilitate biological processes but are also important as protein structural elements. Accurately predicting metal ion interactions in computational studies has always been a challenge, and various methods have been suggested to improve these interactions. One such method is the 12-6-4 Lennard-Jones (LJ)-type nonbonded model. Using this model, it has been possible to successfully reproduce the experimental properties of metal ions in aqueous solution. The model includes induced dipole interactions typically ignored in the standard 12-6 LJ nonbonded model. In this we expand the applicability of this model to metal ion-carboxylate interactions. Using 12-6-4 parameters that reproduce the solvation free energies of the metal ions leads to an overestimation of metal ion-acetate interactions, thus, prompting us to fine-tune the model to specifically handle the latter. We also show that the standard 12-6 LJ model significantly falls short in reproducing the experimental binding free energy between acetate and 11 metal ions (Ni(II), Mg(II), Cu(II), Zn(II), Co(II), Cu(I), Fe(II), Mn(II), Cd(II), Ca(II), and Ag(I)). In this study, we describe optimized C4 parameters for the 12-6-4 LJ nonbonded model to be used with three widely employed water models (Transferable Intermolecular Potential with 3 Points (TIP3P), Simple Point Charge Extended (SPC/E), and Optimal Point Charge (OPC) water models). These parameters can accurately match the experimental binding free energy between 11 metal ions and acetate. These parameters can be applied to the study of metalloproteins and transition metal ion channels and transporters, as acetate serves as a representative of the negatively charged amino acid side chains from aspartate and glutamate.

Characterizing the Molecular Mechanism of the Lethal C423D Mutation in FgMyoI: A Molecular Perspective

The lethal mutation C423D in Fusarium graminearum myosin I (FgMyoI) occurs close to the binding pocket of the allosteric inhibitor phenamacril and causes severe inhibition on mycelial growth of F. graminearum strain PH-1. Here, based on extensive Gaussian accelerated molecular dynamics simulations and wet experiments, we elucidate the underlying molecular mechanism of the abnormal functioning of the FgMyoIC423D mutant at the atomistic level. Our results suggest that the damaging mutation C423D exhibits a synergistic allosteric inhibition mechanism similar to but more robust than that of phenamacril, including effects on the active site and actin binding. Unlike phenamacril-induced closure of Switch2, the mutation results in unfolding of the N-terminal relay helix with a partially opened Switch2 and blocks the structural rearrangement of the relay/SH1 helices, impairing the proper initiation of the recovery stroke. Due to the significant influence of C423D mutation on the function of FgMyoI, designing covalent inhibitors targeting this site holds tremendous potential.

DP/MM: A Hybrid Model for Zinc-Protein Interactions in Molecular Dynamics

Zinc-containing proteins are vital for many biological processes, yet accurately modeling them using classical force fields is hindered by complicated polarization and charge transfer effects. This study introduces DP/MM, a hybrid force field scheme that utilizes a deep potential model to correct the atomic forces of zinc ions and their coordinated atoms, elevating them from MM to QM levels of accuracy. Trained on the difference between MM and QM atomic forces across diverse zinc coordination groups, the DP/MM model faithfully reproduces structural characteristics of zinc coordination during simulations, such as the tetrahedral coordination of Cys4 and Cys3His1 groups. Furthermore, DP/MM allows water exchange in the zinc coordination environment. With its unique blend of accuracy, efficiency, flexibility, and transferability, DP/MM serves as a valuable tool for studying structures and dynamics of zinc-containing proteins and also represents a pioneering approach in the evolving landscape of machine learning potentials for molecular modeling.

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.

Unambiguous discrimination of all 20 proteinogenic amino acids and their modifications by nanopore

Natural proteins are composed of 20 proteinogenic amino acids and their post-translational modifications (PTMs). However, due to the lack of a suitable nanopore sensor that can simultaneously discriminate between all 20 amino acids and their PTMs, direct sequencing of protein with nanopores has not yet been realized. Here, we present an engineered hetero-octameric Mycobacterium smegmatis porin A (MspA) nanopore containing a sole Ni2+ modification. It enables full discrimination of all 20 proteinogenic amino acids and 4 representative modified amino acids, Nω,N'ω-dimethyl-arginine (Me-R), O-acetyl-threonine (Ac-T), N4-(β-N-acetyl-D-glucosaminyl)-asparagine (GlcNAc-N) and O-phosphoserine (P-S). Assisted by machine learning, an accuracy of 98.6% was achieved. Amino acid supplement tablets and peptidase-digested amino acids from peptides were also analyzed using this strategy. This capacity for simultaneous discrimination of all 20 proteinogenic amino acids and their PTMs suggests the potential to achieve protein sequencing using this nanopore-based strategy.

Structure and properties variations of regenerated cellulose fibers induced by metal ion impurity

In the ionic liquids (ILs) method for processing regenerated cellulose fiber (RCF), which is a high-performance ecologically benign product, metal ion impurities (such as Fe3+ and Cu2+) of cellulose might inevitably remain in the recycled ILs and coagulation bath. The presence of metal ions might negatively impact the properties of the manufactured RCFs and obstruct their applications, which are urgent to be made clear. For this research, the solvent for dissolving wood pulp cellulose (WPC) was 1-ethyl-3-methylimidazolium diethyl phosphate ([Emim]DEP) with various metal ion concentrations. The effect of metal ions in IL on the dissolution of cellulose was investigated by Molecular Dynamics simulations. Rheological analysis and degree of polymerization (DP) analysis were adopted to evaluate the influence on fiber spinnability of different spinning solution metal ion concentrations and various dissolving times. Further, the morphology and mechanical performances of the RCFs variation regulation were also thoroughly researched. The findings showed that the presence of metal ions in the spinning solution affected the DP, crystallinity, and orientation factor of RCFs, which will influence their stress more sensitively than the strain. These findings can serve as a practical guide for the commercial manufacture of emerging fiber.

Molecular binding of different classes of organophosphates to methyl parathion hydrolase from Ochrobactrum species

Methyl parathion hydrolase (MPH) is an enzyme of the metallo-β-lactamase superfamily, which hydrolyses a wide range of organophosphates (OPs). Recently, MPH has attracted attention as a promising enzymatic bioremediator. The crystal structure of MPH enzyme shows a dimeric form, with each subunit containing a binuclear metal ion center. MPH also demonstrates metal ion-dependent selectivity patterns. The origins of these patterns remain unclear but are linked to open questions about the more general role of metal ions in functional evolution and divergence within enzyme superfamilies. We aimed to investigate and compare the binding of different OP pesticides to MPH with cobalt(II) metal ions. In this study, MPH was modeled from Ochrobactrum sp. with different OP pesticides bound, including methyl paraoxon and dichlorvos and profenofos. The docked structures for each substrate optimized by DFT calculation were selected and subjected to atomistic molecular dynamics simulations for 500 ns. It was found that alpha metal ions did not coordinate with all the pesticides. Rather, the pesticides coordinated with less buried beta metal ions. It was also observed that the coordination of beta metal ions was perturbed to accommodate the pesticides. The binding free energy calculations and structure-based pharmacophore model revealed that all the three substrates could bind well at the active site. However, profenofos exhibit a stronger binding affinity to MPH in comparison to the other two substrates. Therefore, our findings provide molecular insight on the binding of different OP pesticides which could help us design the enzyme for OP pesticides degradation.

Be2+ Causes Hypersensitivity but Mg2+ and Ca2+ Do Not─Favorable Metal Coordination Is the Key for Differential Allosteric Modulation and Binding Affinities

Although the ion selectivity of metalloproteins has been well established, selective metal antigen recognition by immunoproteins remains elusive. One such case is the recognition of the Be2+ ion against its heavier congeners, Mg2+ and Ca2+, by the human leukocyte antigen immunoprotein (HLA-DP2), leading to immunotoxicity. Integrating with our previous mechanistic study on Be2+ toxicity, herein, we have explored the basis of characteristic nontoxicity of Mg2+ and Ca2+ ions despite their in vivo abundance. The ion binding cleft of the HLA-DP2-peptide complex is composed of four acidic residues, p4D and p7E from the peptide and β26E and β69E from the protein. While the tetrahedral coordination site of the smaller Be2+ ion is located deep inside the cavity, hexa- to octa-coordination sites of Mg2+ and Ca2+ ions are located closer to the protein surface. The intrinsic high coordination number of Mg2+/Ca2+ ions induces allosteric modifications on the HLA-DP2_M2 surface, which are atypical for TCR recognition. Furthermore, the lower binding energy of larger Mg2+ and Ca2+ ions with the cavity residues can be correlated to the lower charge density and reduced covalent bonding nature as compared to those of the smaller Be2+ ion. In short, weak binding of Mg2+ and Ca2+ ions and the unfavorable allosteric surface modifications are probably the major determinants for the absence of Mg2+/Ca2+ ion-mediated hypersensitivity in humans.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches

BACKGROUND

Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear.

RESULTS

This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin.

CONCLUSION

Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

Exploring Major Flavonoid Phytochemicals from Nelumbo nucifera Gaertn. as Potential Skin Anti-Aging Agents: In Silico and In Vitro Evaluations

Nelumbo nucifera Gaertn., an aquatic medicinal plant (Nelumbonaceae family), has a history of use in traditional medicine across various regions. Our previous study demonstrated the skin anti-aging potential of its stamen ethanolic extract by effectively inhibiting collagenase and tyrosinase enzymes. While the major constituents of this extract are well documented, there is a lack of research on the individual compounds' abilities to inhibit skin aging enzymes. Therefore, this study aimed to evaluate the anti-aging potential of the primary flavonoids found in N. nucifera using both in silico and in vitro approaches. Our initial step involved molecular docking to identify compounds with the potential to inhibit collagenase, elastase, and tyrosinase. Among the seven flavonoids studied, kaempferol-3-O-robinobioside (Kae-3-Rob) emerged as the most promising candidate, exhibiting the highest docking scores for three skin aging-related enzymes. Subsequent enzyme-based inhibition assays confirmed that Kae-3-Rob displayed robust inhibitory activity against collagenase (58.24 ± 8.27%), elastase (26.29 ± 7.16%), and tyrosinase (69.84 ± 6.07%). Furthermore, we conducted extensive 200-ns molecular dynamics (MD) simulations, revealing the stability of the complexes formed between Kae-3-Rob and each enzyme along the MD simulation time. MM/PBSA-based binding free energy calculations indicated the considerably stronger binding affinity of Kae-3-Rob for collagenase and tyrosinase compared to elastase, which was related to the greater percentage of hydrogen bond occupations. These computational findings were consistent with the relatively high inhibitory activity of Kae-3-Rob against collagenase and tyrosinase observed in our in vitro experiment. In conclusion, the results obtained from this comprehensive study suggest that Kae-3-Rob, a key flavonoid from N. nucifera, holds significant potential as a source of bioactive compounds for anti-aging cosmeceutical and other phytopharmaceutical application.

Strong Replaces Weak: Design of H-Bond Interactions Enables Cryogenic Aqueous Zn Metal Batteries

Despite the numerous advantages of aqueous Zn batteries, their practical application under cryogenic conditions is hindered by the freezing of the electrolyte because the abundance of hydrogen bonds (H-bonds) between H2O molecules drives the aqueous system to transform to an orderly frozen structure. Here, a design of H-bond interactions based on the guiding ideology of "strong replaces weak" is proposed. The strong H-bonds formed between introduced eutectic components and water molecules break down the weak H-bonds in the original water molecule network, which contributes to an ultralow freezing point and a high ionic conductivity of 1.7 mS cm-1 at -40 °C. Based on multiperspective theoretical simulations and tailor-made in situ cooling Raman characterizations, it has been demonstrated that substituting weak H-bonds with strong H-bonds facilitates the structural reshaping of Zn2+ solvation and remodeling of the H-bond network in the electrolyte. Endowed with this advantage, reversible and stable Zn plating/stripping behaviors could be realized at -40 °C, and the full cells display a high discharge capacity (200 mA h g-1) at -40 °C with ∼75% capacity retention after 1000 cycles. This study will expand the design philosophy of antifreezing aqueous electrolytes and provide a perspective to promote the adoption of Zn metal batteries for cryogenic environment large-scale energy storage.

Mechanism of Fluoride Ion Encapsulation by Magnesium Ions in a Bacterial Riboswitch

Riboswitches sense various ions in bacteria and activate gene expression to synthesize proteins that help maintain ion homeostasis. The crystal structure of the aptamer domain (AD) of the fluoride riboswitch shows that the F- ion is encapsulated by three Mg2+ ions bound to the ligand-binding domain (LBD) located at the core of the AD. The assembly mechanism of this intricate structure is unknown. To this end, we performed computer simulations using coarse-grained and all-atom RNA models to bridge multiple time scales involved in riboswitch folding and ion binding. We show that F- encapsulation by the Mg2+ ions bound to the riboswitch involves multiple sequential steps. Broadly, two Mg2+ ions initially interact with the phosphate groups of the LBD using water-mediated outer-shell coordination and transition to a direct inner-shell interaction through dehydration to strengthen their interaction with the LBD. We propose that the efficient binding mode of the third Mg2+ and F- is that they form a water-mediated ion pair and bind to the LBD simultaneously to minimize the electrostatic repulsion between three Mg2+ bound to the LBD. The tertiary stacking interactions among the LBD nucleobases alone are insufficient to stabilize the alignment of the phosphate groups to facilitate Mg2+ binding. We show that the stability of the whole assembly is an intricate balance of the interactions among the five phosphate groups, three Mg2+, and the encapsulated F- ion aided by the Mg2+ solvated water. These insights are helpful in the rational design of RNA-based ion sensors and fast-switching logic gates.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.

Translocation Processes of Pt(II)-Based Drugs through Human Breast Cancer Cell Membrane: In Silico Experiments

Breast cancer is one of the most frequent modalities of cancer worldwide, with notable mortality. The medication based on platinum drugs (cisplatin (cddp), carboplatin (cpx), and oxaliplatin (oxa)) is a conventional chemotherapy despite severe side effects and the development of drug resistance. In order to provide a deeper molecular description of the influx and efflux processes of platinum drugs through breast cancer tissues, this study focuses on molecular dynamics (MD) simulations of the passive translocation process through a realistic plasma membrane prototype of human breast cancer cell (c_memb). The results showed that the permeation events were mainly mediated by neutral lipids (DOPC, DOPE, and cholesterol), producing a low and temporary membrane deformation. The drug insertion in the region of polar heads was the most favorable stage of the translocation mechanism, especially for cddp and oxa with potential wells of -8.6 and -9.8 kcal mol-1, respectively. However, the potentials of mean force (PMF) revealed unfavorable kinetics for the permeation of these drugs through lipid tails, with energy barriers of 28.3 (cddp), 32.2 (cpx), and 30.4 kcal mol-1 (oxa). The low permeability coefficients (P) of cpx and oxa, which were 3 and 1 orders of magnitude inferior than for cddp, resulted from the high energy barriers for their traslocation processes through the membrane. The obtained results provide a more accurate picture of the permeation of Pt(II)-based drugs through breast cancer cells, which may be relevant for the design and evaluation of new platinum complexes.

Insights into substrate coordination and glycosyl transfer of poplar cellulose synthase-8

Cellulose is an abundant cell wall component of land plants. It is synthesized from UDP-activated glucose molecules by cellulose synthase, a membrane-integrated processive glycosyltransferase. Cellulose synthase couples the elongation of the cellulose polymer with its translocation across the plasma membrane. Here, we present substrate- and product-bound cryogenic electron microscopy structures of the homotrimeric cellulose synthase isoform-8 (CesA8) from hybrid aspen (poplar). UDP-glucose binds to a conserved catalytic pocket adjacent to the entrance to a transmembrane channel. The substrate's glucosyl unit is coordinated by conserved residues of the glycosyltransferase domain and amphipathic interface helices. Site-directed mutagenesis of a conserved gating loop capping the active site reveals its critical function for catalytic activity. Molecular dynamics simulations reveal prolonged interactions of the gating loop with the substrate molecule, particularly across its central conserved region. These transient interactions likely facilitate the proper positioning of the substrate molecule for glycosyl transfer and cellulose translocation.

Enhancement by pyrazolones of colistin efficacy against mcr-1-expressing E. coli: an in silico and in vitro investigation

Owing to the emergence of antibiotic resistance, the polymyxin colistin has been recently revived to treat acute, multidrug-resistant Gram-negative bacterial infections. Positively charged colistin binds to negatively charged lipids and damages the outer membrane of Gram-negative bacteria. However, the MCR-1 protein, encoded by the mobile colistin resistance (mcr) gene, is involved in bacterial colistin resistance by catalysing phosphoethanolamine (PEA) transfer onto lipid A, neutralising its negative charge, and thereby reducing its interaction with colistin. Our preliminary results showed that treatment with a reference pyrazolone compound significantly reduced colistin minimal inhibitory concentrations in Escherichia coli expressing mcr-1 mediated colistin resistance (Hanpaibool et al. in ACS Omega, 2023). A docking-MD combination was used in an ensemble-based docking approach to identify further pyrazolone compounds as candidate MCR-1 inhibitors. Docking simulations revealed that 13/28 of the pyrazolone compounds tested are predicted to have lower binding free energies than the reference compound. Four of these were chosen for in vitro testing, with the results demonstrating that all the compounds tested could lower colistin MICs in an E. coli strain carrying the mcr-1 gene. Docking of pyrazolones into the MCR-1 active site reveals residues that are implicated in ligand-protein interactions, particularly E246, T285, H395, H466, and H478, which are located in the MCR-1 active site and which participate in interactions with MCR-1 in ≥ 8/10 of the lowest energy complexes. This study establishes pyrazolone-induced colistin susceptibility in E. coli carrying the mcr-1 gene, providing a method for the development of novel treatments against colistin-resistant bacteria.

PyConSolv: A Python Package for Conformer Generation of (Metal-Containing) Systems in Explicit Solvent

We introduce PyConSolv, a freely available Python package that automates the generation of conformers of metal- and nonmetal-containing complexes in explicit solvent, through classical molecular dynamics simulations. Using a streamlined workflow and interfacing with widely used computational chemistry software, PyConSolv is an all-in-one tool for the generation of conformers in any solvent. Input requirements are minimal; only the geometry of the structure and the desired solvent in xyz (XMOL) format are needed. The package can also account for charged systems, by including arbitrary counterions in the simulation. A bonded model parametrization is performed automatically, utilizing AmberTools, ORCA, and Multiwfn software packages. PyConSolv provides a selection of preparametrized solvents and counterions for use in classical molecular dynamics simulations. We show the applicability of our package on a number of (transition-metal-containing) systems. The software is provided open source and free of charge.

Molecular Dynamics Simulations of the Human Ecto-5'-Nucleotidase (h-ecto-5'-NT, CD73): Insights into Protein Flexibility and Binding Site Dynamics

Human ecto-5'-nucleotidase (h-ecto-5'-NT, CD73) is a homodimeric Zn2+-binding metallophosphoesterase that hydrolyzes adenosine 5'-monophosphate (5'-AMP) to adenosine and phosphate. h-Ecto-5'-NT is a key enzyme in purinergic signaling pathways and has been recognized as a promising biological target for several diseases, including cancer and inflammatory, infectious, and autoimmune diseases. Despite its importance as a biological target, little is known about h-ecto-5'-NT dynamics, which poses a considerable challenge to the design of inhibitors of this target enzyme. Here, to explore h-ecto-5'-NT flexibility, all-atom unbiased molecular dynamics (MD) simulations were performed. Remarkable differences in the dynamics of the open (catalytically inactive) and closed (catalytically active) conformations of the apo-h-ecto-5'-NT were observed during the simulations, and the nucleotide analogue inhibitor AMPCP was shown to stabilize the protein structure in the closed conformation. Our results suggest that the large and complex domain motion that enables the h-ecto-5'-NT open/closed conformational switch is slow, and therefore, it could not be completely captured within the time scale of our simulations. Nonetheless, we were able to explore the faster dynamics of the h-ecto-5'-NT substrate binding site, which is mainly located at the C-terminal domain and well conserved among the protein's open and closed conformations. Using the TRAPP ("Transient Pockets in Proteins") approach, we identified transient subpockets close to the substrate binding site. Finally, conformational states of the substrate binding site with higher druggability scores than the crystal structure were identified. In summary, our study provides valuable insights into h-ecto-5'-NT structural flexibility, which can guide the structure-based design of novel h-ecto-5'-NT inhibitors.

Repurposing of World-Approved Drugs for Potential Inhibition against Human Carbonic Anhydrase I: A Computational Study

Human carbonic anhydrases (hCAs) have enzymatic activities for reversible hydration of CO2 and are acknowledged as promising targets for the treatment of various diseases. Using molecular docking and molecular dynamics simulation approaches, we hit three compounds of methyl 4-chloranyl-2-(phenylsulfonyl)-5-sulfamoyl-benzoate (84Z for short), cyclothiazide, and 2,3,5,6-tetrafluoro-4-piperidin-1-ylbenzenesulfonamide (3UG for short) from the existing hCA I inhibitors and word-approved drugs. As a Zn2+-dependent metallo-enzyme, the influence of Zn2+ ion models on the stability of metal-binding sites during MD simulations was addressed as well. MM-PBSA analysis predicted a strong binding affinity of -18, -16, and -14 kcal/mol, respectively, for these compounds, and identified key protein residues for binding. The sulfonamide moiety bound to the Zn2+ ion appeared as an essential component of hCA I inhibitors. Vina software predicted a relatively large (unreasonable) Zn2+-sulfonamide distance, although the relative binding strength was reproduced with good accuracy. The selected compounds displayed potent inhibition against other hCA isoforms of II, XIII, and XIV. This work is valuable for molecular modeling of hCAs and further design of potent inhibitors.

A transferable double exponential potential for condensed phase simulations of small molecules

The Lennard-Jones potential is the most widely-used function for the description of non-bonded interactions in transferable force fields for the condensed phase. This is not because it has an optimal functional form, but rather it is a legacy resulting from when computational expense was a major consideration and this potential was particularly convenient numerically. At present, it persists because the effort that would be required to re-write molecular modelling software and train new force fields has, until now, been prohibitive. Here, we present Smirnoff-plugins as a flexible framework to extend the Open Force Field software stack to allow custom force field functional forms. We deploy Smirnoff-plugins with the automated Open Force Field infrastructure to train a transferable, small molecule force field based on the recently-proposed double exponential functional form, on over 1000 experimental condensed phase properties. Extensive testing of the resulting force field shows improvements in transfer free energies, with acceptable conformational energetics, run times and convergence properties compared to state-of-the-art Lennard-Jones based force fields.

Ca2+ ions facilitate the organization of the Annexin A2/S100A10 heterotetramer

Annexin A2 (A2) is a member of the Annexin family, which contains Ca2+ -regulated phospholipid-binding proteins. Annexins associate with S100 proteins to form heterotetramers. The A2/S100A10 heterotetramer (A2t) is the most extensively studied of these heterotetramers. It induces membrane microdomain formation, causes membrane budding, and facilitates proliferation of some cancers. In this work, the first molecular dynamics (MD) study on the complete A2t of 868 amino acids was performed. MD trajectories of more than 600 ns each were generated for complete A2t complexes with and without Ca2+ ions. The outward extension of membrane-binding residues A2-K279 and A2-K281 was shown to be inhibited in the absence of Ca2+ as they were captured by Ca2+ -binding residue D322. F-actin binding residue A2-D339 was observed to occupy either an exposed or buried state in the absence of Ca2+ , while it only occupied the buried state in the presence of Ca2+ . The observed motions of the A2t subunits are highly organized with a strongly correlated central region which is negatively correlated with the periphery of the complex. The central region contains the S100A10 (p11) dimer, A2-N, and A2-I, while the periphery contains A2-II, A2-III, and A2-IV. Novel interactions between A2 and p11 were identified. A2 residues outside of A2-N (K80, R77, E82, and R145) had strong interactions with p11. Residue R145 of A2 may have a significant effect on the dynamics of the system, with its interaction resulting in asymmetric motions of A2. The presented results provide novel insights to inform future experimental studies.

The Modulation of Sucrose Nonfermenting 1-Related Protein Kinase 2.6 State by Persulfidation and Phosphorylation: Insights from Molecular Dynamics Simulations

SnRK2.6 (SUCROSE NONFERMENTING 1-RELATED PROTEIN KINASE2.6) has been characterized as a molecular switch for the intracellular abscisic acid (ABA) signal-transduction pathway. Normally, SnRK2.6 is kept in an "off" state, forming a binary complex with protein phosphatase type 2Cs (PP2Cs). Upon stressful conditions, SnRK2.6 turns into an "on" state by its release from PP2Cs and then phosphorylation at Ser175. However, how the "on" and "off" states for SnRK2.6 are fine-tuned, thereby controlling the initiation and braking processes of ABA signaling, is still largely unclear. SnRK2.6 activity was tightly regulated through protein post-translational modifications (PTM), such as persulfidation and phosphorylation. Taking advantage of molecular dynamics simulations, our results showed that Cys131/137 persulfidation on SnRK2.6 induces destabilized binding and weakened interactions between SnRK2.6 and HAB1 (HYPERSENSITIVE TO ABA1), an important PP2C family protein. This unfavorable effect on the association of the SnRK2.6-HAB1 complex suggests that persulfidation functions are a positive regulator of ABA signaling initiation. In addition, Ser267 phosphorylation in persulfidated SnRK2.6 renders a stable physical association between SnRK2.6 and HAB1, a key characterization for SnRK2.6 inhibition. Rather than Ser175, HAB1 cannot dephosphorylate Ser267 in SnRK2.6, which implies that the retained phosphorylation status of Ser267 could ensure that the activated SnRK2.6 reforms the binary complex to cease ABA signaling. Taken together, our findings expand current knowledge concerning the regulation of persulfidation and phosphorylation on the state transition of SnRK2.6 and provide insights into the fine-tuned mechanism of ABA signaling.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7, TPM1, and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes such as MYH2, TPM2, and TNNI2 that encode parts of the skeletal muscle sarcomere cause muscle diseases affecting skeletal muscle, such as distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7) encoding sarcomeric proteins in which the same pathogenic variant affects skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain pathogenic variants that also cause cardiac abnormalities. We report five families with DA because of heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 (ACTC1). ACTC1 encodes a highly conserved actin that binds to myosin in cardiac and skeletal muscle. Pathogenic variants in ACTC1 have been found previously to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition because of variants in ACTC1 and suggests that some functions of ACTC1 are shared in cardiac and skeletal muscle.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Molecular Dynamics Study of the Effect of Charge and Glycosyl on Superoxide Anion Distribution near Lipid Membrane

To examine the effects of membrane charge, the electrolyte species and glycosyl on the distribution of negatively charged radical of superoxide anion (·O₂^-) around the cell membrane, different phospholipid bilayer systems containing ·O₂^- radicals, different electrolytes and phospholipid bilayers were constructed through Charmm-GUI and Amber16. These systems were equilibrated with molecular dynamics by using Gromacs 5.0.2 to analyze the statistical behaviors of ·O₂^- near the lipid membrane under different conditions. It was found that in the presence of potassium rather than sodium, the negative charge of the phospholipid membrane is more likely to rarefy the superoxide anion distribution near the membrane surface. Further, the presence of glycosyl significantly reduced the density of ·O₂^- near the phospholipid bilayer by 78.3% compared with that of the neutral lipid membrane, which may have a significant contribution to reducing the lipid peroxidation from decreasing the ·O₂^- density near the membrane.

Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

Monovalent Ionic Atmosphere Modulates the Selection of Suboptimal RNA Sequences by Splicing Factors' RNA Recognition Motifs

The U2AF2 splicing factor is involved in the RNA recognition of the pre-mRNA poly-pyrimidine signaling sequence. This protein contains two RRM domains connected by a flexible linker, which ensure the preferential selection of a poly-uridine sequence over a poly-cytosine one. In this work, all-atom simulations provide insights into the U2AF2 recognition mechanism and on the features underlying its selectivity. Our outcomes show that U2AF2's RNA recognition is driven by cooperative events modulated by RNA-protein and RNA-ion interactions. Stunningly, monovalent ions contribute to mediating the binding of the weakly binding polyC strand, thus contributing to the selection of suboptimal poly-pyrimidine tracts. This finding broadens our understanding of the diverse traits tuning splicing factors' selectivity and adaptability to precisely handle and process diverse pre-mRNA sequences.

Uncovering the Mechanism of the Proton-Coupled Fluoride Transport in the CLCF Antiporter

Fluoride is a natural antibiotic abundantly present in the environment and, in micromolar concentrations, is able to inhibit enzymes necessary for bacteria to survive. However, as is the case with many antibiotics, bacteria have evolved resistance methods, including through the use of recently discovered membrane proteins. One such protein is the CLC^F F^-/H⁺ antiporter protein, a member of the CLC superfamily of anion-transport proteins. Though previous studies have examined this F^- transporter, many questions are still left unanswered. To reveal details of the transport mechanism used by CLC^F, we have employed molecular dynamics simulations and umbrella sampling calculations. Our results have led to several discoveries, including the mechanism of proton import and how it is able to aid in the fluoride export. Additionally, we have determined the role of the previously identified residues Glu118, Glu318, Met79, and Tyr396. This work is among the first studies of the CLC^F F^-/H⁺ antiporter and is the first computational investigation to model the full transport process, proposing a mechanism which couples the F^- export with the H⁺ import.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Comparison of Empirical Zn2+ Models in Protein–DNA Complexes

Zinc ions are the second most abundant ions found in humans. Their role in proteins can be merely structural but also catalytic, owing to their transition metal character. Modelling their geometric–coordination versatility by empirical force fields is, thus, a challenging task. In this work, we evaluated three popular models, specifically designed to represent zinc ions with regard to their capability of preserving structural integrity. To this end, we performed molecular dynamics simulations of two zinc-containing protein–DNA complexes, which differed in their zinc coordination, i.e., four cysteines or two cysteines and two histidines. The most flexible non-bonded 12-6-4 Lennard–Jones-type model shows a preference for six-fold coordination of the Zn2+-ions in contradiction to the crystal structure. The cationic dummy atom model favours tetrahedral geometry, whereas the bonded extended zinc AMBER force field model, by construction, best preserves the initial geometry of a regular or slightly distorted tetrahedron. Our data renders the extended zinc AMBER force field the best model for structural zinc ions in a given geometry. In more complicated cases, though, more flexible models may be advantageous.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7 , TPM1 , and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes MYH2 , TPM2 , and TNNI2 , that encode parts of the skeletal muscle sarcomere, cause muscle diseases affecting skeletal muscle, such as the distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7 ) encoding sarcomeric proteins in which the same pathogenic variant affects both skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain mutations that also cause cardiac abnormalities. We report five families with DA due to heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 ( ACTC1 ). ACTC1 encodes a highly conserved actin that binds to myosin in both cardiac and skeletal muscle. Mutations in ACTC1 have previously been found to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition due to mutations in ACTC1 and suggests that some functions of actin, alpha, cardiac muscle 1 are shared in cardiac and skeletal muscle.

Mechanism of Chiral-Selective Aminoacylation of an RNA Minihelix Explored by QM/MM Free-Energy Simulations

Aminoacylation of a primordial RNA minihelix composed of D-ribose shows L-amino acid preference over D-amino acid without any ribozymes or enzymes. This preference in the amino acylation reaction likely plays an important role in the establishment of homochirality in L-amino acid in modern proteins. However, molecular mechanisms of the chiral selective reaction remain unsolved mainly because of difficulty in direct observation of the reaction at the molecular scale by experiments. For seeking a possible mechanism of the chiral selectivity, quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations of the aminoacylation reactions in a modeled RNA were performed to investigate differences in their free-energy profiles along the reactions for L- and D-alanine and its physicochemical origin. The reaction is initiated by approaching a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. The QM/MM umbrella sampling MD calculations showed that the height of the free-energy barrier for L-alanine aminoacylation reaction was 17 kcal/mol, which was 9 kcal/mol lower than that for the D-alanine system. At the transition state, the distance between the negatively charged 3′-oxygen and the positively charged amino group of L-alanine was shorter than that of D-alanine, which was caused by the chirality difference of the amino acid. These results indicate that the transition state for L-alanine is more electrostatically stabilized than that for D-alanine, which would be a plausible mechanism previously unexplained for chiral selectivity in the RNA minihelix aminoacylation.

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

We report QDπ-v1.0 for modeling the internal energy of drug molecules containing H, C, N, and O atoms. The QDπ model is in the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP) that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model that is corrected to a quantitatively high-level of accuracy through a deep-learning potential (DeepPot-SE). The model has the advantage that it is able to properly treat electrostatic interactions and handle changes in charge/protonation states. The model is trained against reference data computed at the ωB97X/6-31G* level (as in the ANI-1x data set) and compared to several other approximate semiempirical and machine learning potentials (ANI-1x, ANI-2x, DFTB3, MNDO/d, AM1, PM6, GFN1-xTB, and GFN2-xTB). The QDπ model is demonstrated to be accurate for a wide range of intra- and intermolecular interactions (despite its intended use as an internal energy model) and has shown to perform exceptionally well for relative protonation/deprotonation energies and tautomers. An example application to model reactions involved in RNA strand cleavage catalyzed by protein and nucleic acid enzymes illustrates QDπ has average errors less than 0.5 kcal/mol, whereas the other models compared have errors over an order of magnitude greater. Taken together, this makes QDπ highly attractive as a potential force field model for drug discovery.

Drug repurposing and molecular mechanisms of the antihypertensive drug candesartan as a TMEM16A channel inhibitor

TMEM16A, a Ca²⁺-activated chloride channel (CaCC), and its pharmacological inhibitors can inhibit the growth of lung adenocarcinoma cells. However,the poor efficacy, safety, and stability of TMEM16A inhibitors limit the development of these agents. Therefore, finding new therapeutic directions from already marketed drugs is a feasible strategy to obtain safe and effective therapeutic drugs. Here, we screened a library contain more than 2400 FDA, EMA, and NMPA-approved drugs through virtual screening. We identified a drug candidate, candesartan (CDST), which showed strong inhibitory effect on the TMEM16A in a concentration-dependent manner with an IC₅₀ of 24.40 ± 3.21 μM. In addition, CDST inhibited proliferation, migration and induced apoptosis of LA795 cells targeting TMEM16A, and significantly inhibited lung adenocarcinoma tumor growth in vivo. The molecular mechanism of CDST inhibiting TMEM16A channel indicated it bound to R515/R535/E623/E624 in the drug pocket, thereby blocked the pore. In conclusion, we identified a novel TMEM16A channel inhibitor, CDST, which exhibited excellent inhibitory activity against lung adenocarcinoma. Considering that CDST has been used in clinical treatment of hypertension, it may play an important role in the combined treatment of hypertension and lung adenocarcinoma as a multi-target drug in the future.

A modified bonded model approach for molecular dynamics simulations of New Delhi Metallo-β-lactamase

Modelling metalloproteins using the classical force fields is challenging. Several methods have been devised to model metalloproteins in force fields. Of these methods, the bonded model, combined with Restrained Electrostatic Potential (RESP) charge fitting, proved its superiority. The latter method was facilitated by the development of the python-based Metal Centre Parameter Builder (MCPB.py) AmberTool. However, the standard bonded model method offered by the MCPB.py tool may not be appropriate for validating and refining the binding modes predicted by docking when crystal structures are lacking. That is because the representation of coordination interactions between any bound ligand and metal ions by covalent bonds can hinder the flexibility of the ligand. Therefore, a new modification to the standard bonded model approach is proposed here. Molecular dynamics (MD) simulations based on the new modified bonded model (MBM) approach avoid the bias caused by coordination bonds and, unlike hybrid QM/MM MD, allow for sufficient sampling of the binding mode given the currently available computational power. The MBM MD approach reproduced the studied crystal structure conformations of New Delhi Metallo-β-lactamase 1 (NDM-1). Furthermore, the MBM approach described the binding interactions of intact β-lactams with NDM-1 reasonably, and predicted a non-productive binding mode for the poor NDM-1 substrate aztreonam whilst predicting productive binding modes for known good substrates. This study presents a useful MD method for metallo-β-lactamases and provides better understanding of β-lactam substrates recognition by NDM-1. The proposed MBM approach might also be useful in the investigation of other metal-containing protein targets.

A Combined in Silico and Structural Study Opens New Perspectives on Aliphatic Sulfonamides, a Still Poorly Investigated Class of CA Inhibitors

Aliphatic sulfonamides are an interesting class of carbonic anhydrase inhibitors (CAIs) proven to be effective for several carbonic anhydrase (CA) isoforms involved in pathologic states. Here we report the crystallographic structures of hCA II in complex with two aliphatic sulfonamides incorporating coumarin rings, which showed a good inhibition and selectivity for this isoform. Although these two molecules have a very similar chemical structure, differing only in the substitution of the two aliphatic hydrogen atoms with two fluorine atoms, they adopt a significantly different binding mode within the enzyme active site. Theoretical binding free energy calculations, performed to rationalize these data, showed that a delicate balance of electrostatic and steric effects modulate the protein-ligand interactions. Data presented here can be fruitfully used for the rational design of novel and effective isozyme-specific inhibitor molecules.