New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions

Structural basis for malate-driven, pore lipid-regulated activation of the Arabidopsis vacuolar anion channel ALMT9

In plant cells, ALMTs are key plasma and vacuolar membrane-localized anion channels regulating plant responses to the environment. Vacuolar ALMTs control anion accumulation in plant cells and, in guard cells, they regulate stomata aperture. The activation of vacuolar ALMTs depends on voltage and cytosolic malate, but the underlying molecular mechanisms remain elusive. Here we report the cryo-EM structures of ALMT9 from Arabidopsis thaliana (AtALMT9), a malate-activated vacuolar anion channel, in plugged and unplugged lipid-bound states. In all these states, membrane lipids interact with the ion conduction pathway of AtALMT9. We identify two unplugged states presenting two distinct pore width profiles. Combining structural and functional analysis we identified conserved residues involved in ion conduction and in the pore lipid interaction. Molecular dynamics simulations revealed a peculiar anion conduction mechanism in AtALMT9. We propose a voltage-dependent activation mechanism based on the competition between pore lipids and malate at the cytosolic entrance of the channel.

Role of Divalent Ions in Membrane Models of Polymyxin-Sensitive and Resistant Gram-Negative Bacteria

Polymyxins, critical last-resort antibiotics, impact the distribution of membrane-bound divalent cations in the outer membrane of Gram-negative bacteria. We employed atomistic molecular dynamics simulations to model the effect of displacing these ions. Two polymyxin-sensitive and two polymyxin-resistant models of the outer membrane of Salmonella enterica were investigated. First, we found that the removal of all calcium ions induces global stress on the model membranes, leading to substantial membrane restructuring. Next, we used enhanced sampling methods to explore the effects of localized stress by displacing membrane-bound ions. Our findings indicate that creating defects in the membrane-bound ion network facilitates polymyxin permeation. Additionally, our study of polymyxin-resistant mutations revealed that divalent ions in resistant model membranes are less likely to be displaced, potentially contributing to the increased resistance associated with these mutations. Lastly, we compared results from all-atom molecular dynamics simulations with coarse-grained simulations, demonstrating that the choice of force field significantly influences the behavior of membrane-bound ions under stress.

Targeting the G-quadruplex as a novel strategy for developing antibiotics against hypervirulent drug-resistant Staphylococcus aureus

BACKGROUND

The rapid emergence of multiple drug-resistant (MDR) bacterial pathogens and the lack of a novel antibiotic pipeline pose a serious threat to global healthcare. The limited number of established targets further restricts the identification of novel antibiotics to treat life-threatening MDR infections caused by Staphylococcus aureus strains. Therefore, novel targets for developing antibiotics are urgently required. In this study, we hypothesized that the G-quadruplex (G4)-binding ligands can be used as novel antibiotics as their binding can possibly downregulate/block the expression of vital genes.

METHODS

To test this, first we screened the antibiotic properties of representative G4-binding ligands against hypervirulent and MDR S. aureus USA300 and determined the in vitro and in vivo antibacterial activity; and proposed the mechanism of action by applying various microbiological, infection, microscopic, and biophysicochemical techniques.

RESULTS

Herein, among screened G4-binding ligands, N-methyl mesoporphyrin IX (NMM) showed the highest antibacterial activity against S. aureus USA300. NMM exhibited a minimum inhibitory concentration (MIC) of 5 μM against S. aureus USA300, impacting cell division and the cell wall by repressing the expressions of genes in the division cell wall (dcw) gene cluster. Genome-wide bioinformatics analysis of G4 motifs and their mapping on S. aureus genome, identified the presence of G4-motif in the promoter of mraZ, a conserved master regulator of the dcw cluster regulating the coordinated cell division and cell wall synthesis. Physicochemical assessments using UV-visible, circular dichroism, and nuclear magnetic resonance spectroscopy confirmed that the G4-motif present in the mraZ promoter formed an intramolecular parallel G4 structure, interacting with NMM. In vivo reporter followed by coupled in vitro transcription/translation (IVT) assays confirmed the role of mraZ G4 as a target interacting NMM to impose extreme antibacterial activity against both the gram-positive and -negative bacteria. In-cell and in vivo validation of NMM using RAW264.7 cells and Galleria mellonella; respectively, demonstrated that NMM exhibited superior antibiotic activity compared to well-established antibiotics, with no observed cytotoxicity.

CONCLUSIONS

In summary, the current study identified NMM as a broad-spectrum potent antibacterial agent and elucidated its plausible mechanism of action primarily by targeting G4-motif in the mraZ promoter of the dcw gene cluster.

Molecular Interactions of the Pioneer Transcription Factor GATA3 With DNA

The GATA3 transcription factor is a pioneer transcription factor that is critical in the development, proliferation, and maintenance of several immune cell types. Identifying the detailed conformational dynamics and interactions of this transcription factor, as well as its clinically important population variants will allow us to unravel its mode of action. In this study, we analyze the molecular interactions of the GATA3 transcription factor bound to dsDNA as well as three clinically important population variants by atomistic molecular dynamics simulations. We identify the effect of the variants on the DNA conformational dynamics and delineate the differences compared to the wildtype transcription factor that could be related to impaired function. We highlight the structural plasticity in the binding of the GATA3 transcription factor and identify important DNA-protein contacts. Although the DNA-protein contacts are persistent and appear to be stable, they exhibit nanosecond timescale fluctuations and several binding/unbinding events. Further, we identify differential DNA binding in the three variants and show that the N-terminal binding is reduced in two of the variants. Our results indicate that reduced minor groove width and DNA diameter are important hallmarks for the binding of GATA3. Our work is an important step towards understanding the functional dynamics of the GATA3 protein and its clinically significant population variants.

Split Membrane: A New Model to Accelerate All-Atom MD Simulation of Phospholipid Bilayers

All-atom molecular dynamics simulations are powerful tools for studying cell membranes and their interactions with proteins and other molecules. However, these processes occur on time scales determined by the diffusion rate of phospholipids, which are challenging to achieve in all-atom models. Here, we present a new all-atom model that accelerates lipid diffusion by splitting phospholipid molecules into head and tail groups. The bilayer structure is maintained by using external lateral potentials, which compensate for the lipid split. This split model enhances lateral lipid diffusion more than ten times, allowing faster and cheaper equilibration of large systems with different phospholipid types. The current model has been tested on membranes containing PSM, POPC, POPS, POPE, POPA, and cholesterol. We have also evaluated the interaction of the split model membranes with the Disheveled DEP domain and amphiphilic helix motif of the transcriptional repressor Opi1 as representative of peripheral proteins as well as the dimeric fragment of the epidermal growth factor receptor transmembrane domain and the Human A2A Adenosine of G protein-coupled receptors as representative of transmembrane proteins. The split model can predict the interaction sites of proteins and their preferred phospholipid type. Thus, the model could be used to identify lipid binding sites and equilibrate large membranes at an affordable computational cost.

How Binding Site Flexibility Promotes RNA Scanning by TbRGG2 RRM: A Molecular Dynamics Simulation Study

RNA recognition motifs (RRMs) are a key class of proteins that primarily bind single-stranded RNAs. In this study, we applied standard atomistic molecular dynamics simulations to obtain insights into the intricate binding dynamics between uridine-rich RNAs and TbRGG2 RRM using the recently developed OL3-Stafix AMBER force field, which improves the description of single-stranded RNA molecules. Complementing structural experiments that unveil a primary binding mode with a single uridine bound, our simulations uncover two supplementary binding modes in which adjacent nucleotides encroach upon the binding pocket. This leads to a unique molecular mechanism through which the TbRGG2 RRM is capable of rapidly transitioning the U-rich sequence. In contrast, the presence of non-native cytidines induces stalling and destabilization of the complex. By leveraging extensive equilibrium dynamics and a large variety of binding states, TbRGG2 RRM effectively expedites diffusion along the RNA substrate while ensuring robust selectivity for U-rich sequences despite featuring a solitary binding pocket. We further substantiate our description of the complex dynamics by simulating the fully spontaneous association process of U-rich sequences to the TbRGG2 RRM. Our study highlights the critical role of dynamics and auxiliary binding states in interface dynamics employed by RNA-binding proteins, which is not readily apparent in traditional structural studies but could represent a general type of binding strategy employed by many RNA-binding proteins.

Understanding Folding of bFGF and Potential Cellular Protective Mechanisms of Neural Cells

Traumatic brain injury (TBI) is a serious health condition that affects an increasing number of people, especially veterans and athletes. TBI causes serious consequences because of its long-lasting impact on the brain and its alarming frequency of occurrence. Although the brain has some natural protective mechanisms, the processes that trigger them are poorly understood. Fibroblast growth factor (FGF) proteins interact with receptor proteins to protect cells. Gaps in the literature include how basic-FGF (bFGF) is activated by heparin, can heparin-like molecules induce neural protection, and the effect of allosteric binding on bFGF activity. To fill the gap in our understanding, we applied temperature replica exchange to study the influence of heparin binding to bFGF and how mutations in bFGF influence stability. A new favorable binding site was identified by comparing free energies computed from the potential of mean force (PMF). Although the varied sugars studied resulted in different interactions with bFGF compared to heparin, they each produced structural effects similar to those of bFGF that likely facilitate receptor binding and signaling. Our results also demonstrate how point mutations can trigger the same conformational change that is believed to promote favorable interactions with the receptor. A deeper atomic-level understanding of how chemicals are released during TBI is needed to improve the development of new treatments for TBI and could contribute to a better understanding of other diseases.

Molecular Dynamics Study of Polyacrylamide and Polysaccharide-Derived Flocculants Adsorption on Mg(OH)2 Surfaces at pH 11

Brucite (Mg(OH)2) is a typical precipitate in the mining industry that adversely affects processes such as flotation and thickening. Gaining insights into the physicochemical properties of this mineral is critical for developing strategies to mitigate these challenges and improve operational efficiency. Additionally, incorporating natural-origin polymers aligns with the shift toward more sustainable mining practices. In this study, molecular dynamics simulations were employed to investigate the interaction of brucite with polysaccharides such as cellulose, guar gum, and alginate and to compare these with conventional polymers, including polyacrylamide, hydrolyzed polyacrylamide, and polyacrylic acid, under conditions of pH 11 in low-salinity water. The methodology enhanced adsorption sampling by incorporating additional temporary interactions between the polymer and the brucite surface. The results reveal that neutral polymers exhibit stronger and more stable interactions with brucite compared to charged polymers, which is consistent with the neutral nature of brucite under the studied conditions. Van der Waals forces predominantly govern the adsorption of polysaccharides, while Coulombic forces primarily drive interactions involving polyacrylamides. These findings provide valuable insights into the molecular mechanisms of polymer-brucite interactions, facilitating the development of more effective and sustainable mining additives.

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of the UUCG RNA Tetraloop and Other Systems

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (ff). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ffs using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ffs did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ffs, OL3CP-gHBfix21, DES-Amber, and OL3R2.7, successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ffs, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG°fold) varied significantly among these three ffs, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3CP-gHBfix21, DES-Amber, and OL3R2.7, respectively. The ΔG°fold value predicted by the OL3CP-gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher ΔG°fold values obtained using DES-Amber and OL3R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3R2.7 ffs on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ffs and highlights key areas for future ff development.

Single-molecule dynamic structural biology with vertically arranged DNA on a fluorescence microscope

The intricate interplay between DNA and proteins is key for biological functions such as DNA replication, transcription and repair. Dynamic nanoscale observations of DNA structural features are necessary for understanding these interactions. Here we introduce graphene energy transfer with vertical nucleic acids (GETvNA), a method to investigate DNA-protein interactions that exploits the vertical orientation adopted by double-stranded DNA on graphene. This approach enables the dynamic study of DNA conformational changes via energy transfer from a probe dye to graphene, achieving spatial resolution down to the Ångström scale at subsecond temporal resolution. We measured DNA bending induced by adenine tracts, bulges, abasic sites and the binding of endonuclease IV. In addition, we observed the translocation of the O6-alkylguanine DNA alkyltransferase on DNA, reaching single base-pair resolution and detecting preferential binding to adenine tracts. This method promises widespread use for dynamical studies of nucleic acids and nucleic acid-protein interactions with resolution so far reserved for traditional structural biology techniques.

Spatiotemporal characterization of water diffusion anomalies in saline solutions using machine learning force field

Understanding water behavior in salt solutions remains a notable challenge in computational chemistry. Conventional force fields have shown limitations in accurately representing water's properties across different salt types (chaotropes and kosmotropes) and concentrations, demonstrating the need for better methods. Machine learning force field applications in computational chemistry, especially through deep potential molecular dynamics (DPMD), offer a promising alternative that closely aligns with the accuracy of first-principles methods. Our research used DPMD to study how salts affect water by comparing its results with ab initio molecular dynamics, SPC/Fw, AMOEBA, and MB-Pol models. We studied water's behavior in salt solutions by examining its spatiotemporally correlated movement. Our findings showed that each model's accuracy in depicting water's behavior in salt solutions is strongly connected to spatiotemporal correlation. This study demonstrates both DPMD's advanced abilities in studying water-salt interactions and contributes to our understanding of the basic mechanisms that control these interactions.

The Influence of Ionic Environment on Nucleosome-Mica Interactions Revealed via Molecular Dynamics Simulations

Nucleosomes are the fundamental units of DNA compaction, playing a key role in modulating gene expression. As such, they are widely studied through both experimental and computational methods. While atomic force microscopy (AFM) is a powerful tool for visualizing and characterizing both canonical and modified nucleosomes, it relies on nucleosome interactions with mica surfaces. These interactions occur through cations adsorbed on the negatively charged mica, but the specific influences of monovalent and divalent cations on nucleosome adsorption remain unclear. In this study, we used molecular dynamics simulations to investigate how monovalent potassium ions and divalent magnesium ions affect nucleosome binding to mica surfaces. We also explored the impact of pretreated mica surfaces on nucleosome binding and structure. Our findings reveal that nucleosome-mica interactions depend on the type of cations present, which leads to distinct effects on nucleosome structure. Notably, nucleosomes bind effectively to mica surfaces in the presence of potassium ions with minimal structural perturbations.

Compared Inhibitory Activities of Tamoxifen and Avenanthramide B on Liver Esterase and Correlation Based on the Superimposed Structure Between Porcine and Human Liver Esterase

Exposure to tamoxifen can exert effects on the human liver, and esterases process prodrugs such as antibiotics and convert them to less toxic metabolites. In this study, the porcine liver esterase (PLE)-inhibitory activity of tamoxifen has been investigated. PLE showed inhibition of a PLE isoenzyme (PLE5). In addition, avenanthramides, which have a similar structure to that of tamoxifen, have been used to determine the PLE-inhibitory effect. Among the avenanthramide derivatives, avenanthramide B has been shown to inhibit PLE. Avenanthramide B interacts with Lys284 of PLE, whereas avenanthramide A and C counteract with Lys284. Avenanthramide B has shown a similar inhibitory effect to that of tamoxifen. Given that avenanthramide B can modulate the action of PLE, it can be used in pharmaceutical and industrial applications for modulating the effects of PLE. Based on superimposed structures between PLE and human liver esterase, the impact of tamoxifen use in humans is discussed. In addition, this study can serve as a fundamental basis for future investigations regarding the potential risk of tamoxifen and other drugs. Thus, this study presents an insight into the comparison of structurally similar tamoxifen and avenanthramides on liver esterases, which can have implications for the pharmaceutical and agricultural industries.

Cysteine hyperoxidation rewires communication pathways in the nucleosome and destabilizes the dyad

Gene activity is tightly controlled by reversible chemical modifications called epigenetic marks, which are of various types and modulate gene accessibility without affecting the DNA sequence. Despite an increasing body of evidence demonstrating the role of oxidative-type modifications of histones in gene expression regulation, there remains a complete absence of structural data at the atomistic level to understand the molecular mechanisms behind their regulatory action. Owing to μs time-scale MD simulations and protein communication networks analysis, we describe the impact of histone H3 hyperoxidation (i.e., S-sulfonylation) on the nucleosome core particle dynamics. Our results reveal the atomic-scale details of the intrinsic structural networks within the canonical histone core and their perturbation by hyperoxidation of the histone H3 C110. We show that this modification involves local rearrangements of the communication networks and destabilizes the dyad, and that one modification is enough to induce a maximal structural signature. Our results suggest that cysteine hyperoxidation in the nucleosome core particle might favor its disassembly.

Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics

All-atom constant pH molecular dynamics simulations offer a powerful tool for understanding pH-mediated and proton-coupled biological processes. As the protonation equilibria of protein side chains are shifted by electrostatic interactions and desolvation energies, pKa values calculated from the constant pH simulations may be sensitive to the underlying protein force field and water model. Here we investigated the force field dependence of the all-atom particle mesh Ewald (PME) continuous constant pH (PME-CpHMD) simulations of a mini-protein BBL. The replica-exchange titration simulations based on the Amber ff19sb and ff14sb force fields with the respective water models showed significantly overestimated pKa downshifts for a buried histidine (His166) and for two glutamic acids (Glu141 and Glu161) that are involved in salt-bridge interactions. These errors (due to undersolvation of neutral histidines and overstabilization of salt bridges) are consistent with the previously reported pKa's based on the CHARMM c22/CMAP force field, albeit in larger magnitudes. The pKa calculations also demonstrated that ff19sb with OPC water is significantly more accurate than ff14sb with TIP3P water, and the salt-bridge related pKa downshifts can be partially alleviated by the atom-pair specific Lennard-Jones corrections (NBFIX). Together, these data suggest that the accuracies of the protonation equilibria of proteins from constant pH simulations can significantly benefit from improvements of force fields.

Bending DNA increases its helical repeat

In all biological systems, DNA is under high mechanical stress from bending and twisting. For example, DNA is tightly bent in nucleosome complexes, virus capsids, bacterial chromosomes, or complexes with transcription factors that regulate gene expression. A structurally and mechanically accurate model of DNA is therefore necessary to understand some of the most fundamental molecular mechanisms in biology including DNA packaging, replication, transcription and gene regulation. An iconic feature of DNA is its double helical nature with an average repeath 0 of ~10.45 base pairs per turn, which is commonly believed to be independent of curvature. We developed a ligation assay on nicked DNA circles of variable curvature that reveals a strong unwinding of DNA to over 11 bp/turn for radii around 3-4 nm. Our work constitutes a major modification of the standard mechanical model of DNA and requires reassessing the molecular mechanisms and energetics of all processes involving tightly bent DNA.

On the specificity of the recognition of m6A-RNA by YTH reader domains

Most processes of life are the result of polyvalent interactions between macromolecules, often of heterogeneous types and sizes. Frequently, the times associated with these interactions are prohibitively long for interrogation using atomistic simulations. Here, we study the recognition of N6-methylated adenine (m6A) in RNA by the reader domain YTHDC1, a prototypical, cognate pair that challenges simulations through its composition and required timescales. Simulations of RNA pentanucleotides in water reveal that the unbound state can impact (un)binding kinetics in a manner that is both model- and sequence-dependent. This is important because there are two contributions to the specificity of the recognition of the Gm6AC motif: from the sequence adjacent to the central adenine and from its methylation. Next, we establish a reductionist model consisting of an RNA trinucleotide binding to the isolated reader domain in high salt. An adaptive sampling protocol allows us to quantitatively study the dissociation of this complex. Through joint analysis of a data set including both the cognate and control sequences (GAC, Am6AA, and AAA), we derive that both contributions to specificity, sequence, and methylation, are significant and in good agreement with experimental numbers. Analysis of the kinetics suggests that flexibility in both the RNA and the YTHDC1 recognition loop leads to many low-populated unbinding pathways. This multiple-pathway mechanism might be dominant for the binding of unstructured polymers, including RNA and peptides, to proteins when their association is driven by polyvalent, electrostatic interactions.

Molecular Dynamics Simulations of Nucleosomes Containing Histone Variant H2A.J

Histone proteins form the building blocks of chromatin-nucleosomes. Incorporation of alternative histone variants instead of the major (canonical) histones into nucleosomes is a key mechanism enabling epigenetic regulation of genome functioning. In humans, H2A.J is a constitutively expressed histone variant whose accumulation is associated with cell senescence, inflammatory gene expression, and certain cancers. It is sequence-wise very similar to the canonical H2A histones, and its effects on the nucleosome structure and dynamics remain elusive. This study employed all-atom molecular dynamics simulations to reveal atomistic mechanisms of structural and dynamical effects conferred by the incorporation of H2A.J into nucleosomes. We showed that the H2A.J C-terminal tail and its phosphorylated form have unique dynamics and interaction patterns with the DNA, which should affect DNA unwrapping and the availability of nucleosomes for interactions with other chromatin effectors. The dynamics of the L1-loop and the hydrogen bonding patterns inside the histone octamer were shown to be sensitive to single amino acid substitutions, potentially explaining the higher thermal stability of H2A.J nucleosomes. Taken together, our study demonstrated unique dynamical features of H2A.J-containing nucleosomes, which contribute to further understanding of the molecular mechanisms employed by H2A.J in regulating genome functioning.

Single-Molecule-Based, Label-Free Monitoring of Molecular Glue Efficacies for Promoting Protein-Protein Interactions Using YaxAB Nanopores

Modulating protein-protein interactions (PPIs) is an attractive strategy in drug discovery. Molecular glues, bifunctional small-molecule drugs that promote PPIs, offer an approach to targeting traditionally undruggable targets. However, the efficient discovery of molecular glues has been hampered by the current limitations of conventional ensemble-averaging-based methods. In this study, we present a YaxAB nanopore for probing the efficacy of molecular glues in inducing PPIs. Using YaxAB nanopores, we demonstrate single-molecule-based, label-free monitoring of protein complex formation between mammalian target of rapamycin (mTOR) and FK506-binding proteins (FKBPs) triggered by the molecular glue, rapamycin. Owing to its wide entrance and adjustable pore size, in combination with potent electro-osmotic flow (EOF), a single funnel-shaped YaxAB nanopore enables the simultaneous detection and single-molecule-level quantification of multiprotein states, including single proteins, binary complexes, and ternary complexes induced by rapamycin. Notably, YaxAB nanopores could sensitively discriminate between the binary complexes or ternary complexes induced by rapamycin and its analogues, despite the subtle size differences of ∼122 or ∼116 Da, respectively. Taken together, our results provide proof-of-concept for single-molecule-based, label-free, and ultrasensitive screening and structure-activity relationship (SAR) analysis of molecular glues, which will contribute to low-cost, highly efficient discovery, and rational design of bifunctional modality of drugs, such as molecular glues.

Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics

All-atom constant pH molecular dynamics simulations offer a powerful tool for understanding pH-mediated and proton-coupled biological processes. As the protonation equilibria of protein sidechains are shifted by electrostatic interactions and desolvation energies, pK a values calculated from the constant pH simulations may be sensitive to the underlying protein force field and water model. Here we investigated the force field dependence of the all-atom particle mesh Ewald (PME) continuous constant pH (PME-CpHMD) simulations of a mini-protein BBL. The replica-exchange titration simulations based on the Amber ff19sb and ff14sb force fields with the respective water models showed significantly overestimated pK a downshifts for a buried histidine (His166) and for two glutamic acids (Glu141 and Glu161) that are involved in salt-bridge interactions. These errors (due to undersolvation of neutral histidines and overstabilization of salt bridges) are consistent with the previously reported pK a's based on the CHARMM c22/CMAP force field, albeit in larger magnitudes. The pK a calculations also demonstrated that ff19sb with OPC water is significantly more accurate than ff14sb with TIP3P water, and the salt-bridge related pK a downshifts can be partially alleviated by the atom-pair specific Lennard-Jones corrections (NBFIX). Together, these data suggest that the accuracies of the protonation equilibria of proteins from constant pH simulations can significantly benefit from improvements of force fields.

Computer-Aided Identification and Design of Ligands for Multi-Targeting Inhibition of a Molecular Acute Myeloid Leukemia Network

BACKGROUND/OBJECTIVES

Acute myeloid leukemia (AML) is characterized by therapeutic failure and long-term risk for disease relapses. As several therapeutic targets participate in networks, they can rewire to eventually evade single-target drugs. Hence, multi-targeting approaches are considered on the expectation that interference with many different components could synergistically hinder activation of alternative pathways and demolish the network one-off, leading to complete disease remission.

METHODS

Herein, we established a network-based, computer-aided approach for the rational design of drug combinations and de novo agents that interact with many AML network components simultaneously.

RESULTS

A reconstructed AML network guided the selection of suitable protein hubs and corresponding multi-targeting strategies. For proteins responsive to existing drugs, a greedy algorithm identified the minimum amount of compounds targeting the maximum number of hubs. We predicted permissible combinations of amiodarone, artenimol, fostamatinib, ponatinib, procaine, and vismodegib that interfere with 3-8 hubs, and we elucidated the pharmacological mode of action of procaine on DNMT3A. For proteins that do not respond to any approved drugs, namely cyclins A1, D2, and E1, we used structure-based de novo drug design to generate a novel triple-targeting compound of the chemical formula C15H15NO5, with favorable pharmacological and drug-like properties.

CONCLUSIONS

Overall, by integrating network and structural pharmacology with molecular modeling, we determined two complementary strategies with the potential to annihilate the AML network, one in the form of repurposable drug combinations and the other as a de novo synthesized triple-targeting agent. These target-drug interactions could be prioritized for preclinical and clinical testing toward precision medicine for AML.

Charge Scaling in Classical Force Fields for Lithium Ions in Polymers

Polymer electrolytes are of interest for applications in energy storage. Molecular simulations of ion transport in polymer electrolytes have been widely used to study the conductivity in these materials. Such simulations have generally relied on classical force fields. A peculiar feature of such force fields has been that in the particular case of lithium ions (Li+), their charge must be scaled down by approximately 20% to achieve agreement with experimental measurements of ion diffusivity. In this work, we present first-principles calculations that serve to justify the charge-scaling factor and van der Waals interaction parameters for Li+ diffusion in poly(ethylene glycol) (PEO) with bistriflimide (TFSI-) counterions. Our results indicate that a scaling factor of 0.79 provides good agreement with DFT calculations over a relatively wide range of Li+ concentrations and temperatures, consistent with past reports where that factor was adjusted by trial and error. We also show that such a scaling factor leads to diffusivities that are in quantitative agreement with experimental measurements.

In silico study of DNA mononucleotide self-assembly

Recent experiments have demonstrated the self-assembly and long-range ordering of concentrated aqueous solutions of DNA and RNA mononucleotides. These are found to form Watson-Crick pairs that stack into columns that become spatially organized into a columnar liquid-crystalline phase. In this work, we numerically investigate this phase behavior by adopting an extremely coarse-grained model in which nucleotides are represented as semi-disk-like polyhedra decorated with attractive (patchy) sites that mimic the stacking and pairing interactions. We carry out Monte Carlo simulations of these patchy polyhedra by adapting algorithms borrowed from computer graphics. This model reproduces the features of the experimental phase behavior, which essentially depends on the combination of pairing and stacking interactions.

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters via Ab Initio Molecular Dynamics Simulations To Accurately Model Glycosaminoglycan Electrostatic Interactions

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

Complexes of HMO1 with DNA: Structure and Affinity

Saccharomyces cerevisiae HMO1 is an architectural nuclear DNA-binding protein that stimulates the activity of some remodelers and regulates the transcription of ribosomal protein genes, often binding to a DNA motif called IFHL. However, the molecular mechanism dictating this sequence specificity is unclear. Our circular dichroism spectroscopy studies show that the HMO1:DNA complex forms without noticeable changes in the structure of DNA and HMO1. Molecular modeling/molecular dynamics studies of the DNA complex with HMO1 Box B reveal two extended sites at the N-termini of helices I and II of Box B that are involved in the formation of the complex and stabilize the DNA bend induced by intercalation of the F114 side chain between base pairs. A comparison of the affinities of HMO1 for 24 bp DNA fragments containing either randomized or IFHL sequences reveals a twofold increase in the stability of the complex in the latter case, which may explain the selectivity in the recognition of the IFHL-containing promoter regions.

Cation-induced intramolecular coil-to-globule transition in poly(ADP-ribose)

Poly(ADP-ribose) (PAR), a non-canonical nucleic acid, is essential for DNA/RNA metabolism and protein condensation, and its dysregulation is linked to cancer and neurodegeneration. However, key structural insights into PAR's functions remain largely uncharacterized, hindered by the challenges in synthesizing and characterizing PAR, which are attributed to its length heterogeneity. A central issue is how PAR, comprised solely of ADP-ribose units, attains specificity in its binding and condensing proteins based on chain length. Here, we integrate molecular dynamics simulations with small-angle X-ray scattering to analyze PAR structures. We identify diverse structural ensembles of PAR that fall into distinct subclasses and reveal distinct compaction of two different lengths of PAR upon the addition of small amounts of Mg2+ ions. Unlike PAR15, PAR22 forms ADP-ribose bundles via local intramolecular coil-to-globule transitions. Understanding these length-dependent structural changes could be central to deciphering the specific biological functions of PAR.

Ubiquitinated histone H2B as gatekeeper of the nucleosome acidic patch

Monoubiquitination of histones H2B-K120 (H2BK120ub) and H2A-K119 (H2AK119ub) play opposing roles in regulating transcription and chromatin compaction. H2BK120ub is a hallmark of actively transcribed euchromatin, while H2AK119ub is highly enriched in transcriptionally repressed heterochromatin. Whereas H2BK120ub is known to stimulate the binding or activity of various chromatin-modifying enzymes, this post-translational modification (PTM) also interferes with the binding of several proteins to the nucleosome H2A/H2B acidic patch via an unknown mechanism. Here, we report cryoEM structures of an H2BK120ub nucleosome showing that ubiquitin adopts discrete positions that occlude the acidic patch. Molecular dynamics simulations show that ubiquitin remains stably positioned over this nucleosome region. By contrast, our cryoEM structures of H2AK119ub nucleosomes show ubiquitin adopting discrete positions that minimally occlude the acidic patch. Consistent with these observations, H2BK120ub, but not H2AK119ub, abrogates nucleosome interactions with acidic patch-binding proteins RCC1 and LANA, and single-domain antibodies specific to this region. Our results suggest a mechanism by which H2BK120ub serves as a gatekeeper to the acidic patch and point to distinct roles for histone H2AK119 and H2BK120 ubiquitination in regulating protein binding to nucleosomes.

nP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat

Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.

Structure of the human dopamine transporter and mechanisms of inhibition

The neurotransmitter dopamine has central roles in mood, appetite, arousal and movement1. Despite its importance in brain physiology and function, and as a target for illicit and therapeutic drugs, the human dopamine transporter (hDAT) and mechanisms by which it is inhibited by small molecules and Zn2+ are without a high-resolution structural context. Here we determine the structure of hDAT in a tripartite complex with the competitive inhibitor and cocaine analogue, (-)-2-β-carbomethoxy-3-β-(4-fluorophenyl)tropane2 (β-CFT), the non-competitive inhibitor MRS72923 and Zn2+ (ref. 4). We show how β-CFT occupies the central site, approximately halfway across the membrane, stabilizing the transporter in an outward-open conformation. MRS7292 binds to a structurally uncharacterized allosteric site, adjacent to the extracellular vestibule, sequestered underneath the extracellular loop 4 (EL4) and adjacent to transmembrane helix 1b (TM1b), acting as a wedge, precluding movement of TM1b and closure of the extracellular gate. A Zn2+ ion further stabilizes the outward-facing conformation by coupling EL4 to EL2, TM7 and TM8, thus providing specific insights into how Zn2+ restrains the movement of EL4 relative to EL2 and inhibits transport activity.

Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Tailoring Phage Nanosomes for Enhanced Theranostic Properties of Near Infrared Dyes

Near-infrared (NIR) phototherapies offer noninvasive, cost-effective solutions for treating tumors and microbial infections. However, organic NIR dyes commonly used suffer from solubility and stability issues requiring frequent dosing. We address this challenge by exploring the bacteriophage-mediated enhancement of NIR dye properties. Upon encapsulation within phage nanosomes, IR780 and Indocyanine green (ICG), with similar optical properties but distinct water solubility and exhibit enhanced UV-vis absorbance and photothermal transduction efficacy compared to liposomes. Experimental characterization corroborated with all-atom molecular dynamics simulations imprints the nanoscale structure, solubility, dynamics, and binding of these NIR dye molecules to the membrane and protein molecules present in Phage capsid. These NIR dye-loaded phage nanosomes, coencapsulated with mitoxantrone, demonstrate enhanced anticancer activity, and when combined with amphotericin B, these dye molecules exhibit superior photothermal effects against fungal infections. Our findings present a simple and efficient approach for tuning the photothermal performance of existing NIR dyes through a rational design for enhanced therapeutic outcomes.

Pathways to a Shiny Future: Building the Foundation for Computational Physical Chemistry and Biophysics in 2050

In the last quarter-century, the field of molecular dynamics (MD) has undergone a remarkable transformation, propelled by substantial enhancements in software, hardware, and underlying methodologies. In this Perspective, we contemplate the future trajectory of MD simulations and their possible look at the year 2050. We spotlight the pivotal role of artificial intelligence (AI) in shaping the future of MD and the broader field of computational physical chemistry. We outline critical strategies and initiatives that are essential for the seamless integration of such technologies. Our discussion delves into topics like multiscale modeling, adept management of ever-increasing data deluge, the establishment of centralized simulation databases, and the autonomous refinement, cross-validation, and self-expansion of these repositories. The successful implementation of these advancements requires scientific transparency, a cautiously optimistic approach to interpreting AI-driven simulations and their analysis, and a mindset that prioritizes knowledge-motivated research alongside AI-enhanced big data exploration. While history reminds us that the trajectory of technological progress can be unpredictable, this Perspective offers guidance on preparedness and proactive measures, aiming to steer future advancements in the most beneficial and successful direction.

Insights into mechanisms of MALT1 allostery from NMR and AlphaFold dynamic analyses

Mucosa-associated lymphoid tissue lymphoma-translocation protein 1 (MALT1) is an attractive target for the development of modulatory compounds in the treatment of lymphoma and other cancers. While the three-dimensional structure of MALT1 has been previously determined through X-ray analysis, its dynamic behaviour in solution has remained unexplored. We present here dynamic analyses of the apo MALT1 form along with the E549A mutation. This investigation used NMR 15N relaxation and NOE measurements between side-chain methyl groups. Our findings confirm that MALT1 exists as a monomer in solution, and demonstrate that the domains display semi-independent movements in relation to each other. Our dynamic study, covering multiple time scales, along with the assessment of conformational populations by Molecular Dynamic simulations, Alpha Fold modelling and PCA analysis, put the side chain of residue W580 in an inward position, shedding light at potential mechanisms underlying the allosteric regulation of this enzyme.

Asymmetric ion transport through "Janus" MoSSe sub-nanometer pores

We conduct all-atom molecular dynamics simulations to systematically investigate the underlying mechanisms governing ion transport through a sub-nanometer pore decorated with negative charges in a "Janus" MoSSe membrane. The charge imbalance between S and Se atoms on each side of the membrane induces different types of ion adsorption processes depending on the pore inner charge configuration, and the polarity of external biases, which leads to asymmetry in ionic I-V characteristics. Statistical analysis of the total translocation times including adsorption-desorption processes, and ion dwell times indicates that potassium ions predominantly remain adsorbed during their interaction with the membrane before undertaking a quick translocation through the pore. High applied biases suppress cation adsorption, which results in fast translocation with the current flow boosted by negative inner charges around the pore. We also show that in a membrane consisting of several "Janus" layers, the applied bias necessary to overcome the sub-nm pore barrier increases with the number of layers, providing control over the ionic current.

Electro-osmotic Flow Generation via a Sticky Ion Action

Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.

How arginine inhibits substrate-binding domain 2 elucidated using molecular dynamics simulations

The substrate-binding domain 2 (SBD2) is an important part of the bacterial glutamine (GLN) transporter and mediates binding and delivery of GLN to the transporter translocation subunit. The SBD2 consists of two domains, D1 and D2, that bind GLN in the space between domains in a closed structure. In the absence of ligand, the SBD2 adopts an open conformation with larger space between domains. The GLN binding and closing are essential for the subsequent transport into the cell. Arginine (ARG) can also bind to SBD2 but does not induce closing and inhibits GLN transport. We use atomistic molecular dynamics (MD) simulations in explicit solvent to study ARG binding in the presence of the open SBD2 structure and observed reversible binding to the native GLN binding site with similar contacts but no transition to a closed SBD2 state. Absolute binding free energy simulations predict a considerable binding affinity of ARG and GLN to the binding site on the D1 domain. Free energy simulations to induce subsequent closing revealed a strong free energy penalty in case of ARG binding in contrast to GLN binding that favors the closed SBD2 state but still retains a free energy barrier for closing. The simulations allowed the identification of the molecular origin of the closing penalty in case of bound ARG and suggested a mutation of lysine at position 373 to alanine that strongly reduced the penalty and allowed closing even in the presence of bound ARG. The study offers an explanation of the molecular mechanism of how ARG competitively inhibits GLN from binding to SBD2 and from triggering the transition to a closed conformation. The proposed Lys373Ala mutation shows promise as a potential tool to validate whether a conformational mismatch between open SBD2 and the translocator is responsible for preventing ARG uptake to the cell.

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

The Influence of Ionic Environment on Nucleosome-Mica Interactions Revealed via Molecular Dynamics Simulations

Mica serves as a crucial substrate in Atomic Force Microscopy (AFM) studies for visualizing and characterizing nucleosomes. Nucleosomes interact with the negatively charged mica surface via adsorbed cations. However, the specific influences of monovalent and divalent cations on nucleosome adsorption to the mica surface remain unclear. In this study, we investigated the binding of nucleosomes to the mica surface in the presence of monovalent potassium ions and divalent magnesium ions using molecular dynamics simulations. We also explored the impact of pre-treated mica surfaces on nucleosome binding and structure. Our findings reveal that nucleosome-mica interactions vary depending on the cations present, resulting in distinct effects on nucleosome structure. Notably, nucleosomes bind effectively to a mica surface in the presence of potassium ions with minimal structural perturbations.

The effect of functional groups on the glass transition temperature of atmospheric organic compounds: a molecular dynamics study

Organic compounds constitute a substantial part of atmospheric particulate matter not only in terms of mass concentration but also in terms of distinct functional groups. The glass transition temperature provides an indirect way to investigate the phase state of the organic compounds, playing a crucial role in understanding their behavior and influence on aerosol processes. Molecular dynamics (MD) simulations were implemented here to predict the glass transition temperature (Tg) of atmospherically relevant organic compounds as well as the influence of their functional groups and length of their carbon chain. The cooling step used in the simulations was chosen to be neither too low (to supress crystallization) nor too high (to avoid Tg overprediction). According to the MD simulations, the predicted Tg is sensitive to the functional groups as follows: carboxylic acid (-COOH) > hydroxyl (-OH) and (-COOH) > carbonyls (-CO). Increasing the number of carbon atoms leads to higher Tg for the linearly structured compounds. Linear compounds with lower molecular weight were found to exhibit a lower Tg. No clear correlation between O : C and Tg was observed. The architecture of the carbon chain (linear, or branched, or ring) was also found to impact the glass transition temperature. Compounds containing a non-aromatic carbon ring are characterized by a higher Tg compared to linear and branched ones with the same number of carbon atoms.

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

Binding to nucleosome poises human SIRT6 for histone H3 deacetylation

Sirtuin 6 (SIRT6) is an NAD+-dependent histone H3 deacetylase that is prominently found associated with chromatin, attenuates transcriptionally active promoters and regulates DNA repair, metabolic homeostasis and lifespan. Unlike other sirtuins, it has low affinity to free histone tails but demonstrates strong binding to nucleosomes. It is poorly understood how SIRT6 docking on nucleosomes stimulates its histone deacetylation activity. Here, we present the structure of human SIRT6 bound to a nucleosome determined by cryogenic electron microscopy. The zinc finger domain of SIRT6 associates tightly with the acidic patch of the nucleosome through multiple arginine anchors. The Rossmann fold domain binds to the terminus of the looser DNA half of the nucleosome, detaching two turns of the DNA from the histone octamer and placing the NAD+ binding pocket close to the DNA exit site. This domain shows flexibility with respect to the fixed zinc finger and moves with, but also relative to, the unwrapped DNA terminus. We apply molecular dynamics simulations of the histone tails in the nucleosome to show that in this mode of interaction, the active site of SIRT6 is perfectly poised to catalyze deacetylation of the H3 histone tail and that the partial unwrapping of the DNA allows even lysines close to the H3 core to reach the enzyme.

Structure and dynamics of double-stranded DNA rotaxanes

A DNA rotaxane, with its unique mechanically interlocked architecture consisting of a circular DNA molecule threaded onto a linear DNA axle, holds promise as a fundamental component for nanoscale functional devices. Nevertheless, its structural and dynamic behaviors, essential for advancing molecular machinery, remain largely unexplored. Using extensive all-atom molecular dynamics simulations, we investigated the behaviors of double-stranded DNA (dsDNA) rotaxanes, concentrating on the effects of shape distortion induced by torsional stress in small circular dsDNA containing 70-90 base pairs. We analyzed structural characteristics, including shape, intermolecular distances, and tilt angles, while also exploring dynamic properties such as translational diffusion and toroidal rotation. Our results indicate that shape distortion brings the circular and linear dsDNA components into closer proximity and causes a slight increase in translational diffusion yet a minor decrease in toroidal rotation. Nevertheless, there is no apparent evidence of coupling between translation and rotation. Overall, the insights from this study indicate that such shape distortion does not significantly alter their structure and dynamics. This finding provides flexibility for the design of DNA rotaxanes in nanoscale applications.

Molecular Determinants of Current Blockade Produced by Peptide Transport Through a Nanopore

The nanopore sensing method holds the promise of delivering a single molecule technology for identification of biological proteins, direct detection of post-translational modifications, and perhaps de novo determination of a protein's amino acid sequence. The key quantity measured in such nanopore sensing experiments is the magnitude of the ionic current passing through a nanopore blocked by a polypeptide chain. Establishing a relationship between the amino acid sequence of a peptide fragment confined within a nanopore and the blockade current flowing through the nanopore remains a major challenge for realizing the nanopore protein sequencing. Using the results of all-atom molecular dynamics simulations, here we compare nanopore sequencing of DNA with nanopore sequencing of proteins. We then delineate the factors affecting the blockade current modulation by the peptide sequence, showing that the current can be determined by (i) the steric footprint of an amino acid, (ii) its interactions with the pore wall, (iii) the local stretching of a polypeptide chain, and (iv) the local enhancement of the ion concentration at the nanopore constriction. We conclude with a brief discussion of the prospects for purely computational prediction of the blockade currents.

Hyaluronan-arginine enhanced and dynamic interaction emerges from distinctive molecular signature due to electrostatics and side-chain specificity

Hyaluronan is a natural carbohydrate polymer with a negative charge that fosters gel-like conditions crucial for its cellular functions and industrial applications. As a recognized ligand for proteins, understanding their mutual interactions provides solid ground to tune hyaluronan's gel properties using biocompatible peptides. This work employs NMR and molecular dynamics simulations to identify molecular motifs relevant to hyaluronan-peptide interactions using arginine, lysine, and glycine oligopeptides. Arginine-rich peptides exhibit the strongest binding to hyaluronan according to chemical shift perturbation measurements, followed distantly by the similarly charged lysine. This difference highlights the significance of electrostatics and the peculiarities of the guanidinium side chain in arginine, capable of non-polar interactions that further stabilize the binding. Additional nuclear Overhauser effect measurements do not show stable interaction partners, precluding strong and well-defined complexes. Finally, molecular simulations support our findings and show an extended but significant interaction region, especially for arginine, responsible for the observed enhanced binding, which can also promote cross-linking of hyaluronan polymers. Our findings pave the way for optimizing biocompatible peptides to alter hyaluronan gels' properties efficiently and also explain why hyaluronan-protein interaction typically involves positively charged arginine-rich regions also capable of forming hydrogen bonds and non-polar interactions.

van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD

Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.