General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Prediction and biological significance of small changes in binding of leiomodin to tropomyosin

In cardiac muscle, regulation of actin polymerization at the thin filament pointed end is controlled by two structurally similar but functionally antagonistic proteins, leiomodin-2 and tropomodulin-1. Both proteins contain tropomyosin-binding site 1, which is essential for their recruitment to the pointed end. Using circular dichroism, we determined changes in melting temperatures (ΔTm) for complexes of tropomyosin and leiomodin-2 fragments containing several hypomorphic mutations, which moderately affect binding to tropomyosin. We ran molecular dynamics simulations for the complexes and calculated standard Gibbs free energies of binding, which we found to strongly correlate with the ΔTm. We found that the E34Q mutation in leiomodin-2 resulted in a decrease in the melting temperature of the complex of tropomyosin and leiomodin-2 fragments, indicating a decrease in the affinity of leiomodin-2 for tropomyosin. Although modest, this change in in vitro affinity made leiomodin-2 a weaker competitor for tropomyosin than tropomodulin-1 in cardiomyocytes. This mutation significantly reduced the ability of leiomodin-2 to displace tropomodulin-1 at thin filament pointed ends and affected the ability of leiomodin-2 to elongate thin filaments. Our results highlight the essential role of the tropomyosin-binding site in the dynamic equilibrium between tropomodulin-1 and leiomodin-2 at the pointed end of thin filaments. Our data also suggest the potential use of the correlation between ΔTm and the modeled standard Gibbs free energies of binding to predict changes in the stability of complexes between tropomyosin and leiomodin or tropomodulin isoforms.

Nanopores with an Engineered Selective Entropic Gate Detect Proteins at Nanomolar Concentration in Complex Biological Sample

Biological nanopores enable the electrical detection of biomolecules, making them ideal sensors for use in health-monitoring devices. Proteins are widely recognized as biomarkers for various diseases, but they present a unique challenge due to their vast diversity and concentration range in biological samples. Here, inspired by the nuclear pore complex, we incorporated a layer of disordered polypeptides into the biological nanopore YaxAB. This polypeptide mesh formed an entropic gate, significantly reducing the entry of proteins from a highly concentrated mixture, including blood. The introduction of a specific recognition element within the disordered polypeptides allowed targeted proteins to penetrate through the nanopores, where they were recognized by specific current signatures. This biosensing approach allowed for the recognition of nanomolar proteins directly from blood samples without prior sample preparation. This work paves the way for the next generation of nanopore sensors for the real-time detection of proteins in blood.

Unraveling the Tether: Exploring Representative Protein Linkers and Their Structural and Thermodynamical Properties

This study explores the thermodynamic and structural behaviors of linker peptides, short polypeptide segments that often bridge protein domains. We are focusing on three prototypical classes─glycine-serine (GS), glycine-glycine (GG), and alanine-proline (AP)─and exploring their conformational dynamics as isolated entities outside a multidomain protein context. Using extensive molecular dynamics (MD) simulations and free energy perturbation (FEP) analyses, we characterize the free energy landscapes, entropic properties, and solvation energetics of 20 representative linkers. Our results reveal a pronounced linear relationship between linker length and key thermodynamic contributions, including zero-point vibrational energy (ZPVE), potential energy, and entropy. Notably, vibrational entropy emerges as a dominant stabilizing term. We also found that AP linkers display more rigid, yet extended conformations compared to the highly flexible GS and moderately flexible GG linkers. These findings underscore the nuanced role of linker composition in contributing to multidomain protein architecture and dynamics, and highlight how thermodynamic forces shape linker conformational behavior. Collectively, our work enhances the mechanistic understanding of protein linkers, offering valuable insights for the rational design of peptide-based systems and informing future efforts to modulate interdomain flexibility and stability in multidomain proteins.

Structural elucidation of four fungal hydrophobins belonging to classes I and II: Results from Alphafold and accelerated molecular dynamics simulations

Hydrophobins (HFB) find application in various industries including biotechnology and medical devices; therefore, it is imperative to elucidate and learn more about their folded structures. Few fungal HFB protein structures are available in the Protein Data Bank (PDB), and fewer have been elucidated using homology modeling or short molecular dynamics (MD) simulations in the literature. Moreover, many homology modeling algorithms will only model the region with sequence identity. Therefore, we turned towards the state-of-the-art, artificial intelligence powered AlphaFold. It performed well in predicting the core β-barrel, a characteristic of HFBs, except for HFB9A which was unfolded with low confidence scores. These initial structures were then prepared for accelerated MD simulation in the hope of observing higher protein folding. With 500 ns long aMD simulations, we were able to obtain folded and energetically stable conformations for all the proteins except HFB9A, which exhibited much higher disorder, connected with higher atomic fluctuation, lowest hydrophobicity and overall compactness, and lesser secondary structure formation as visualized during the aMD simulation. The underlying intrinsic disorder in the HFBs was found to be the basis of harder-to-reach folding by AF2 which can be compensated by enhanced sampling MD simulations like the aMD technique. The characteristic two disordered loops of class I HFBs were obtained for SC3, while HFB9A showed that not all class I HFBs contain them. Class II HFBs were more stable, folded and compact with secondary structure motifs conserved throughout the trajectory which can be correlated with their comparatively much lower intrinsic disorder.

Intrinsically Disordered Proteins Can Behave as Different Polymers across Their Conformational Ensemble

Intrinsically disordered proteins (IDPs) are macromolecules, which in contrast to well-folded proteins explore a large number of conformationally heterogeneous states. In this work, we investigate the conformational space of the disordered protein β-casein using Hamiltonian replica exchange atomistic molecular dynamics (MD) simulations in explicit water. The energy landscape contains a global minimum along with two shallow funnels. Employing static polymeric scaling laws separately for individual funnels, we find that they cannot be described by the same polymeric scaling exponent. Around the global minimum, the conformations are globular, whereas in the vicinity of local minima, we recover coil-like scaling. To elucidate the implications of structural diversity on equilibrium dynamics, we initiated standard MD simulations in the NVT ensemble with representative conformations from each funnel. Global and internal motions for different classes of trajectories show heterogeneous dynamics with globule to coil-like signatures. Thus, IDPs can behave as entirely different polymers in different regions of the conformational space.

Studying the Signaling Mechanism of Neuropilin-1's Intracellular Disorder Region via Conformational Mining and Dynamic Interaction Characterization

Many single-pass membrane proteins contain an intrinsically disordered region (IDR) within their intracellular domain, playing a key role in regulating cellular signaling. However, understanding the functional mechanisms of these disordered regions has remained a challenge. In this study, we focus on the cytoplasmic IDR of neuropilin-1 (NRP-1 IDR) and employ a combination of experimental and computational methods to investigate its dynamics and function. We compare several enhanced sampling molecular simulations, structural statistics-based methods, and AI-driven conformation mining techniques, emphasizing the strengths and limitations of each with respect to sampling diversity and energy landscape exploration. Subsequently, we investigate the broad array of potential binding partners for the NRP-1 IDR and employ AlphaFold3 for complex structure prediction, highlighting the promiscuous binding behavior of the NRP-1 IDR. Finally, we focus on high-confidence binding partners, GIPC-1 and SNX-5, validating the interaction of the NRP-1 IDR with these proteins and investigating the effects of membrane context and phosphorylation on these interactions. Our findings provide critical insights into how a flexible cytoplasmic region in signal-transmembrane proteins can modulate transmembrane signaling.

Optimal scaling of protein-water interactions coupled with targeted torsional refinements yields balanced force fields suitable for simulations of single-chain folded proteins, disordered polypeptides, and protein-protein complexes

All-atom molecular dynamics (MD) simulations based on physics-based force fields, serve as an essential complement to experiments for investigating protein structure, dynamics, and interactions. Despite significant advances in force field development, achieving a consistent balance of molecular interactions that stabilize folded proteins and protein-protein complexes while simultaneously capturing the conformational dynamics of intrinsically disordered polypeptides (IDPs), remains challenging. In this work, we systematically evaluated two current state-of-the-art force fields (i) AMBER ff03ws, and (ii) AMBER ff99SBws, by comprehensively assessing their performance on both folded domains and IDPs. By selectively scaling side chain-water interactions for uncharged residues, the refined AMBER ff03w-sc force field demonstrated improved conformational stability of folded proteins while maintaining accurate representations of IDPs. However, AMBER ff03w-sc failed to correct the discrepancies in NMR-derived ps-ns timescale backbone dynamics associated with flexible loops. Interestingly, AMBER ff99SBws retained its structural stability despite the application of upscaled interactions with water for both sidechain and backbone atoms and displayed robust agreement with NMR-derived backbone dynamics. Further, a targeted refinement of glutamine backbone torsion parameters, yielded AMBER ff99SBws-STQ', which effectively resolved discrepancies associated with glutamine α-helicity predictions. Extensive validation against small angle X-ray scattering (SAXS) and NMR chemical shifts, revealed that both refined force fields accurately reproduced chain dimensions and secondary structure propensities of disordered peptides and prion-like domains. Importantly, both force fields reliably maintained the stability of protein-protein complexes over microsecond timescales. Our systematic refinement strategies provide improved accuracy and transferability for simulating diverse protein systems, from folded domains to IDPs and protein complexes.

Small-Molecule Inhibitors of Amyloid Beta: Insights from Molecular Dynamics-Part A: Endogenous Compounds and Repurposed Drugs

The amyloid hypothesis is the predominant model of Alzheimer's disease (AD) pathogenesis, suggesting that amyloid beta (Aβ) peptide is the primary driver of neurotoxicity and a cascade of pathological events in the central nervous system. Aβ aggregation into oligomers and deposits triggers various processes, such as vascular damage, inflammation-induced astrocyte and microglia activation, disrupted neuronal ionic homeostasis, oxidative stress, abnormal kinase and phosphatase activity, tau phosphorylation, neurofibrillary tangle formation, cognitive dysfunction, synaptic loss, cell death, and, ultimately, dementia. Molecular dynamics (MD) is a powerful structure-based drug design (SBDD) approach that aids in understanding the properties, functions, and mechanisms of action or inhibition of biomolecules. As the only method capable of simulating atomic-level internal motions, MD provides unique insights that cannot be obtained through other techniques. Integrating experimental data with MD simulations allows for a more comprehensive understanding of biological processes and molecular interactions. This review summarizes and evaluates MD studies from the past decade on small molecules, including endogenous compounds and repurposed drugs, that inhibit amyloid beta. Furthermore, it outlines key considerations for future MD simulations of amyloid inhibitors, offering a potential framework for studies aimed at elucidating the mechanisms of amyloid beta inhibition by small molecules.

Grappa - a machine learned molecular mechanics force field

Simulating large molecular systems over long timescales requires force fields that are both accurate and efficient. In recent years, E(3) equivariant neural networks have lifted the tension between computational efficiency and accuracy of force fields, but they are still several orders of magnitude more expensive than established molecular mechanics (MM) force fields. Here, we propose Grappa, a machine learning framework to predict MM parameters from the molecular graph, employing a graph attentional neural network and a transformer with symmetry-preserving positional encoding. The resulting Grappa force field outperforms tabulated and machine-learned MM force fields in terms of accuracy at the same computational efficiency and can be used in existing Molecular Dynamics (MD) engines like GROMACS and OpenMM. It predicts energies and forces of small molecules, peptides, and RNA at state-of-the-art MM accuracy, while also reproducing experimentally measured values for J-couplings. With its simple input features and high data-efficiency, Grappa is well suited for extensions to uncharted regions of chemical space, which we show on the example of peptide radicals. We demonstrate Grappa's transferability to macromolecules in MD simulations from a small fast-folding protein up to a whole virus particle. Our force field sets the stage for biomolecular simulations closer to chemical accuracy, but with the same computational cost as established protein force fields.

The evolution of the Amber additive protein force field: History, current status, and future

Molecular dynamics simulations are pivotal in elucidating the intricate properties of biological molecules. Nonetheless, the reliability of their outcomes hinges on the precision of the molecular force field utilized. In this perspective, we present a comprehensive review of the developmental trajectory of the Amber additive protein force field, delving into researchers' persistent quest for higher precision force fields and the prevailing challenges. We detail the parameterization process of the Amber protein force fields, emphasizing the specific improvements and retained features in each version compared to their predecessors. Furthermore, we discuss the challenges that current force fields encounter in balancing the interactions of protein-protein, protein-water, and water-water in molecular dynamics simulations, as well as potential solutions to overcome these issues.

Effects of All-Atom and Coarse-Grained Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of a Tau K18 Monomer

To propose new mechanism-based therapeutics for Alzheimer's disease (AD), it is crucial to study the kinetics and oligomerization/aggregation mechanisms of the hallmark tau proteins, which have various isoforms and are intrinsically disordered. In this study, multiple all-atom (AA) and coarse-grained (CG) force fields (FFs) have been benchmarked on molecular dynamics (MD) simulations of K18 tau (M243-E372), which is a truncated form (130 residues) of full-length tau (441 residues). FF19SB is first excluded because the dynamics are too slow, and the conformations are too stable. All other benchmarked AAFFs (Charmm36m, FF14SB, Gromos54A7, and OPLS-AA) and CGFFs (Martini3 and Sirah2.0) exhibit a trend of shrinking K18 tau into compact structures with the radius of gyration (ROG) around 2.0 nm, which is much smaller than the experimental value of 3.8 nm, within 200 ns of AA-MD or 2000 ns of CG-MD. Gromos54A7, OPLS-AA, and Martini3 shrink much faster than the other FFs. To perform meaningful postanalysis of various properties, we propose a strategy of selecting snapshots with 2.5 < ROG < 4.5 nm, instead of using all sampled snapshots. The calculated chemical shifts of all C, CA, and CB atoms have very good and close root-mean-square error (RMSE) values, while Charmm36m and Sirah2.0 exhibit better chemical shifts of N than other FFs. Comparing the calculated distributions of the distance between the CA atoms of CYS291 and CYS322 with the results of the FRET experiment demonstrates that Charmm36m is a perfect match with the experiment while other FFs exhibit limitations. In summary, Charmm36m is recommended as the best AAFF, and Sirah2.0 is recommended as an excellent CGFF for simulating tau K18.

Survey of the Aβ-peptide structural diversity: molecular dynamics approaches

The review deals with the application of Molecular Dynamics (MD) to the structure modeling of beta-amyloids (Aβ), currently classified as intrinsically disordered proteins (IDPs). In this review, we strive to relate the main advances in this area but specifically focus on the approaches and methodology. All relevant papers on the Aβ modeling are cited in the Tables in Supplementary Data, including a concise description of the applied approaches, sorted according to the types of the studied systems: modeling of the monomeric Aβ and Aβ aggregates. Similar sections focused according to the type of modeled object are present in the review. In the final part of the review, novel methods of general IDP modeling not confined to Aβ are described.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01253-y.

Impact of water models on the structure and dynamics of enzyme tunnels

Protein hydration plays a vital role in many biological functions, and molecular dynamics simulations are frequently used to study it. However, the accuracy of these simulations is often sensitive to the water model used, a phenomenon particularly evident in intrinsically disordered proteins. Here, we investigated the extent to which the choice of water model alters the behavior of complex networks of tunnels within proteins. Tunnels are essential because they allow the exchange of substrates and products between buried enzyme active sites and the bulk solvent, directly affecting enzyme efficiency and selectivity, making the study of tunnels crucial for a holistic understanding of enzyme function at the molecular level. By performing simulations of haloalkane dehalogenase LinB and its two variants with engineered tunnels using TIP3P and OPC models, we investigated their effects on the overall tunnel topology. We also analyzed the properties of the primary tunnels, including their conformation, bottleneck dimensions, sampling efficiency, and the duration of tunnel openings. Our data demonstrate that all three proteins exhibited similar conformational behavior in both models but differed in the geometrical characteristics of their auxiliary tunnels, consistent with experimental observations. Interestingly, the results indicate that the stability of the open tunnels might be sensitive to the water model used. Because TIP3P can provide comparable data on the overall tunnel network, it is a valid choice when computational resources are limited or compatibility issues impede the use of OPC. However, OPC seems preferable for calculations requiring an accurate description of transport kinetics.

Dimerization of Full-Length Aβ-42 Peptide: A Comparison of Different Force Fields and Water Models

Among the two isoforms of amyloid-β i. e., Aβ-40 and Aβ-42, Aβ-42 is more toxic due to its increased aggregation propensity. The oligomerization pathways of amyloid-β may be investigated by studying its dimerization process at an atomic level. Intrinsically disordered proteins (IDPs) lack well-defined structures and are associated with numerous neurodegenerative disorders. Molecular dynamics simulations of these proteins are often limited by the choice of parameters due to inconsistencies in the empirically developed protein force fields and water models. To evaluate the accuracy of recently developed force fields for IDPs, we study the dimerization of full-length Aβ-42 in aqueous solution with three different combinations of AMBER force field parameters and water models such as ff14SB/TIP3P, ff19SB/OPC, and ff19SB/TIP3P using classical MD and Umbrella Sampling method. This work may be used as a benchmark to compare the performance of different force fields for the simulations of IDPs.

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.

Protamine-Mediated Tangles Produce Extreme Deoxyribonucleic Acid Compaction

In sperm cells, protamine replaces histones to compact DNA 10-20 times more than in somatic cells. To characterize the extreme compaction, we employed confocal microscopy and optical tweezers to determine the conformations and stability of protamine-bound λ-DNA. Confocal images show increasing compaction of λ-DNA at increasing protamine concentration. In the presence of protamine, single λ-DNA molecules form tangles that withstand forces strong enough (∼55 pN) for strand separation and shorten the contour length by up to 40% even at high forces, as well as bends and loops that rupture at 10-40 pN forces. Strand separation nucleates tangles, implicating protamine interactions with DNA bases. Molecular dynamics simulations show that Arg sidechains of protamine each form hydrogen bonds with multiple bases, frequently in the form of a wedge between the two strands of DNA. Protamine may participate in both local and higher-order chromatin organization, leading to extreme compaction and global transcription silencing.

Sequence and structural determinants of RNAPII CTD phase-separation and phosphorylation by CDK7

The intrinsically disordered carboxy-terminal domain (CTD) of the largest subunit of RNA Polymerase II (RNAPII) consists of multiple tandem repeats of the consensus heptapeptide Y1-S2-P3-T4-S5-P6-S7. The CTD promotes liquid-liquid phase-separation (LLPS) of RNAPII in vivo. However, understanding the role of the conserved heptad residues in LLPS is hampered by the lack of direct biochemical characterization of the CTD. Here, we generated a systematic array of CTD variants to unravel the sequence-encoded molecular grammar underlying the LLPS of the human CTD. Using in vitro experiments and molecular dynamics simulations, we report that the aromaticity of tyrosine and cis-trans isomerization of prolines govern CTD phase-separation. The cis conformation of prolines and β-turns in the SPXX motif contribute to a more compact CTD ensemble, enhancing interactions among CTD residues. We further demonstrate that prolines and tyrosine in the CTD consensus sequence are required for phosphorylation by Cyclin-dependent kinase 7 (CDK7). Under phase-separation conditions, CDK7 associates with the surface of the CTD droplets, drastically accelerating phosphorylation and promoting the release of hyperphosphorylated CTD from the droplets. Our results highlight the importance of conformationally restricted local structures within spacer regions, separating uniformly spaced tyrosine stickers of the CTD heptads, which are required for CTD phase-separation.

Harnessing conformational dynamics in enzyme catalysis to achieve nature-like catalytic efficiencies: the shortest path map tool for computational enzyme redesign

Enzymes exhibit diverse conformations, as represented in the free energy landscape (FEL). Such conformational diversity provides enzymes with the ability to evolve towards novel functions. The challenge lies in identifying mutations that enhance specific conformational changes, especially if located in distal sites from the active site cavity. The shortest path map (SPM) method, which we developed to address this challenge, constructs a graph based on the distances and correlated motions of residues observed in nanosecond timescale molecular dynamics (MD) simulations. We recently introduced a template based AlphaFold2 (tAF2) approach coupled with 10 nanosecond MD simulations to quickly estimate the conformational landscape of enzymes and assess how the FEL is shifted after mutation. In this study, we evaluate the potential of SPM when coupled with tAF2-MD in estimating conformational heterogeneity and identifying key conformationally-relevant positions. The selected model system is the beta subunit of tryptophan synthase (TrpB). We compare how the SPM pathways differ when integrating tAF2 with different MD simulation lengths from as short as 10 ns until 50 ns and considering two distinct Amber forcefield and water models (ff14SB/TIP3P versus ff19SB/OPC). The new methodology can more effectively capture the distal mutations found in laboratory evolution, thus showcasing the efficacy of tAF2-MD-SPM in rapidly estimating enzyme dynamics and identifying the key conformationally relevant hotspots for computational enzyme engineering.

Computational studies on the catalytic potential of the double active site for enzyme engineering

Proteins possessing double active sites have the potential to revolutionise enzyme design strategies. This study extensively explored an enzyme that contains both a natural active site (NAS) and an engineered active site (EAS), focusing on understanding its structural and functional properties. Metadynamics simulations were employed to investigate how substrates interacted with their respective active sites. The results revealed that both the NAS and EAS exhibited similar minimum energy states, indicating comparable binding affinities. However, it became apparent that the EAS had a weaker binding site for the substrate due to its smaller pocket and constrained conformation. Interestingly, the EAS also displayed dynamic behaviour, with the substrate observed to move outside the pocket, suggesting the possibility of substrate translocation. To gain further insights, steered molecular dynamics (SMD) simulations were conducted to study the conformational changes of the substrate and its interactions with catalytic residues. Notably, the substrate adopted distinct conformations, including near-attack conformations, in both the EAS and NAS. Nevertheless, the NAS demonstrated superior binding minima for the substrate compared to the EAS, reinforcing the observation that the engineered active site was less favourable for substrate binding due to its limitations. The QM/MM (Quantum mechanics and molecular mechanics) analyses highlight the energy disparity between NAS and EAS. Specifically, EAS exhibited elevated energy levels due to its engineered active site being located on the surface. This positioning exposes the substrate to solvents and water molecules, adding to the energy challenge. Consequently, the engineered enzyme did not provide a significant advantage in substrate binding over the single active site protein. Further, the investigation of internal channels and tunnels within the protein shed light on the pathways facilitating transport between the two active sites. By unravelling the complex dynamics and functional characteristics of this double-active site protein, this study offers valuable insights into novel strategies of enzyme engineering. These findings establish a solid foundation for future research endeavours aimed at harnessing the potential of double-active site proteins in diverse biotechnological applications.

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

Advancements have been made to dissipative particle dynamics (DPD), a robust coarse-grained (CG) simulation method, to study the folded structures of four miniproteins (1L2Y, 1WN8, 1YRF, and 2I9M) in explicit and implicit solvents. In this endeavor, we aim to establish model parametrization and enhance computational efficiency. Unlike traditional CG models that use empirical force parameters, ex-force parameters (r0(ex), a ~ , δd, δp) of DPD particles constructed for specific research purposes can be obtained from atomistic molecular dynamics simulations. On the other hand, im-force parameters (r0(im), c, σ) can be derived from ex-DPD simulations, according to the underlying thermodynamic theory. Based on a mapping scheme proposed for the modeling of amino acids, all-atom proteins can be converted into a CG model. Both ex-/im-DPDs are then carried out to investigate the folding pathways of the four mini-proteins. Structural analysis of the RMSDs shows that the im-simulated proteins have greater structural similarity to native proteins than the ex-simulated ones. The constructed CG models achieve a resolution of Angstrom (Å), a level normally associated with atomic models. Additionally, speed tests reveal that im-DPD accelerates the simulation process and significantly improves simulation efficiency.

Design and Discovery of New Collagen V-Derived FGF2-Blocking Natural Peptides Inhibiting Lung Squamous Cell Carcinoma In Vitro and In Vivo

Aberrant FGF2/FGFR signaling is implicated in lung squamous cell carcinoma (LSCC), posing treatment challenges due to the lack of targeted therapeutic options. Designing drugs that block FGF2 signaling presents a promising strategy different from traditional kinase inhibitors. We previously reported a ColVα1-derived fragment, HEPV (127AA), that inhibits FGF2-induced angiogenesis. However, its large size may limit therapeutic application. This study combines rational peptide design, molecular dynamics simulations, knowledge-based prediction, and GUV and FRET assays to identify smaller peptides with FGF2-blocking properties. We synthesized two novel peptides, HBS-P1 (45AA) and HBS-P2 (66AA), that retained the heparin-binding site. Both peptides demonstrated anti-LSCC and antiangiogenesis properties in cell viability and microvessel network induction assays. In two LSCC subcutaneous models, HBS-P1, with its affinity for FGF2 and enhanced penetration ability, demonstrated substantial therapeutic potential without apparent toxicities. Our study provides the first evidence supporting the development of collagen V-derived natural peptides as FGF2-blocking agents for LSCC treatment.

Transient Non-local Interactions Dominate the Dynamics of Measles Virus NTAIL

The RNA genome of measles virus is encapsidated by the nucleoprotein within a helical nucleocapsid that serves as template for both transcription and replication. The intrinsically disordered domain of the nucleoprotein (NTAIL), partly protruding outward from the nucleocapsid, is essential for binding the polymerase complex responsible for viral transcription and replication. As for many IDPs, binding of NTAIL occurs through a short molecular recognition element (MoRE) that folds upon binding, with the majority of NTAIL remaining disordered. Though NTAIL regions far from the MoRE influence the binding affinity, interactions between them and the MoRE have not been investigated in depth. Using an integrated approach, relying on photo-induced electron transfer (PET) experiments between tryptophan and cysteine pairs placed at different positions in the protein under varying salt and pH conditions, combined with simulations and analytical models, we identified transient interactions between two disordered regions distant in sequence, which dominate NTAIL dynamics, and regulate the conformational preferences of both the MoRE and the entire NTAIL domain. Co-evolutionary analysis corroborates our findings, and suggests an important functional role for the same intramolecular interactions. We propose mechanisms by which these non-local interactions may regulate binding to the phosphoprotein, polymerase recruitment, and ultimately viral transcription and replication. Our findings may be extended to other IDPs, where non-local intra-protein interactions affect the conformational preferences of intermolecular binding sites.

Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data

Molecular dynamics (MD) is a great tool for elucidating conformational dynamics of proteins and peptides in water at the atomistic level that often surpasses the level of detail available experimentally. Structure predictions, however, are limited by the accuracy of the underlying MD force field. This limitation is particularly stark in the case of intrinsically disordered peptides and proteins, which are characterized by solvent-accessible and disordered peptide regions and domains. Recent studies show that most additive MD force fields, including CHARMM36m, do not reproduce the intrinsic conformational distributions of guest amino acid residues x in cationic GxG peptides in water in line with experimental data. Positing that a lack of polarizability in additive MD force fields may be the culprit for the reported discrepancies, we here examine the conformational dynamics of guest glycine and alanine residues in cationic GxG peptides in water using two polarizable MD force fields, CHARMM Drude and AMOEBA. Our results indicate that while AMOEBA captures the experimental data better than CHARMM Drude, neither of the two polarizable force fields offers an improvement of the Ramachandran distributions of glycine and alanine residues in cationic GGG and GAG peptides, respectively, over CHARMM36m.

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-β (Aβ42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the β-content of Aβ42; however, it increases the tendency of the monomeric Aβ42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aβ42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aβ42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

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.

Differentiable simulation to develop molecular dynamics force fields for disordered proteins

Implicit solvent force fields are computationally efficient but can be unsuitable for running molecular dynamics on disordered proteins. Here I improve the a99SB-disp force field and the GBNeck2 implicit solvent model to better describe disordered proteins. Differentiable molecular simulations with 5 ns trajectories are used to jointly optimise 108 parameters to better match explicit solvent trajectories. Simulations with the improved force field better reproduce the radius of gyration and secondary structure content seen in experiments, whilst showing slightly degraded performance on folded proteins and protein complexes. The force field, called GB99dms, reproduces the results of a small molecule binding study and improves agreement with experiment for the aggregation of amyloid peptides. GB99dms, which can be used in OpenMM, is available at https://github.com/greener-group/GB99dms. This work is the first to show that gradients can be obtained directly from nanosecond-length differentiable simulations of biomolecules and highlights the effectiveness of this approach to training whole force fields to match desired properties.

Probing the hierarchical dynamics of DNA-sperm nuclear transition protein complexes through fuzzy interaction and mesoscale condensation

Nuclear transition protein TNP1 is a crucial player mediating histone-protamine exchange in condensing spermatids. A unique combination of intrinsic disorder and multivalent properties turns TNP1 into an ideal agent for orchestrating the formation of versatile TNP-DNA assemblies. Despite its significance, the physicochemical property and the molecular mechanism followed by TNP1 for histone replacement and DNA condensation are still poorly understood. This study reports the first-time in vitro expression and purification of human TNP1 and investigates the hierarchical dynamics of TNP1-DNA interaction using a combination of computational simulations, biochemical assays, fluorescence imaging, and atomic force microscopy. We explored three crucial facets of TNP1-DNA interactions. Initially, we delve into the molecular binding process that entails fuzzy interactions between TNP1 and DNA at the atomistic scale. Subsequently, we analyze how TNP1 binding affects the electrostatic and mechanical characteristics of DNA and influences its morphology. Finally, we study the biomolecular condensation of TNP1-DNA when subjected to high concentrations. The findings of our study set the foundation for comprehending the potential involvement of TNP1 in histone replacement and DNA condensation in spermatogenesis.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

Structure Comparison of Beta Amyloid Peptide Aβ 1-42 Isoforms. Molecular Dynamics Modeling

Beta amyloid peptide Aβ 1-42 (Aβ42) has a unique dual role in the human organism, as both the peptide with an important physiological function and one of the most toxic biological compounds provoking Alzheimer's disease (AD). There are several known Aβ42 isoforms that we discuss here that are highly neurotoxic and lead to the early onset of AD. Aβ42 is an intrinsically disordered protein with no experimentally solved structure under physiological conditions. The objective of this research was to establish the appropriate molecular dynamics (MD) methodology and model a uniform set of structures for the Aβ42 isoforms that form the core of this study. For that purpose, force field selection and verification including convergence testing for MD simulations was made. Replica exchange MD and conventional MD modeling of several Aβ42 and Aβ16 isoforms that have neurotoxic and amyloidogenic effects impacting the severity of Alzheimer's disease were carried out with the optimal force field and solvent parameters. A standardized ensemble of structures for the Aβ42 and Aβ16 isoforms covering 30-50% of the conformational ensembles extracted from the free energy minima was calculated from MD trajectories. The resulting data set of modeled structures includes Aβ42 wild type, isoD7, pS8, D7H, and H6R-Aβ42 and Aβ16 wild type, isoD7, pS8, D7H, and H6R-Aβ16. The representative structures are given in the Supporting Information; they are open for public access. In the study, we also evaluated the differences between the structures of Aβ42 isoforms and speculate on their possible relevance to the known functions. Utilizing several representative structures for a single disordered protein for docking, with their subsequent averaging by conformations, would markedly increase the reliability of docking results.

Unraveling the Behavior of Intrinsically Disordered Protein c-Myc: A Study Utilizing Gaussian-Accelerated Molecular Dynamics

The protein c-Myc is a transcription factor that remains largely intrinsically disordered and is known to be involved in various biological processes and is overexpressed in various cancers, making it an attractive drug target. However, intrinsically disordered proteins such as c-Myc do not show funnel-like basins in their free-energy landscapes; this makes their druggability a challenge. For the first time, we propose a heterodimer model of c-Myc/Max in full length in this work. We used Gaussian-accelerated molecular dynamics (GaMD) simulations to explore the behavior of c-Myc and its various regions, including the transactivation domain (TAD) and the basic helix-loop-helix-leucine-zipper (bHLH-Zipper) motif in three different conformational states: (a) monomeric c-Myc, (b) c-Myc when bound to its partner protein, Max, and (c) when Max was removed after binding. We analyzed the GaMD trajectories using root-mean-square deviation (RMSD), radius of gyration, root-mean-square fluctuation, and free-energy landscape (FEL) calculations to elaborate the behaviors of these regions. The results showed that the monomeric c-Myc structure showed a higher RMSD fluctuation as compared with the c-Myc/Max heterodimer in the bHLH-Zipper motif. This indicated that the bHLH-Zipper motif of c-Myc is more stable when it is bound to Max. The TAD region in both monomeric and Max-bound states showed similar plasticity in terms of RMSD. We also conducted residue decomposition calculations and showed that the c-Myc and Max interaction could be driven mainly by electrostatic interactions and the residues Arg299, Ile403, and Leu420 seemed to play important roles in the interaction. Our work provides insights into the behavior of c-Myc and its regions that could support the development of drugs that target c-Myc and other intrinsically disordered proteins.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

Spidroins under the Influence of Alcohol: Effect of Ethanol on Secondary Structure and Molecular Level Solvation of Silk-Like Proteins

Future sustainable materials based on designer biomolecules require control of the solution assembly, but also interfacial interactions. Alcohol treatments of protein materials are an accessible means to this, making understanding of the process at the molecular level of seminal importance. We focus here on the influence of ethanol on spidroins, the main proteins of silk. By large-scale atomistically detailed molecular dynamics (MD) simulations and interconnected experiments, we characterize the protein aggregation, secondary structure changes, molecular level origins of them, and solvation environment changes for the proteins, as induced by ethanol as a solvation additive. The MD and circular dichoroism (CD) findings jointly show that ethanol promotes ordered structure in the protein molecules, leading to an increase of helix content and turns but also increased aggregation, as revealed by dynamic light scattering (DLS) and light microscopy. The structural changes correlate at the molecular level with increased intramolecular hydrogen bonding. The simulations reveal that polar amino acids, such as glutamine and serine, are most influenced by ethanol, whereas glycine residues are most prone to be involved in the ethanol-induced secondary structure changes. Furthermore, ethanol engages in interactions with the hydrophobic alanine-rich regions of the spidroin, significantly decreasing the hydrophobic interactions of the protein with itself and its surroundings. The protein solutes also change the microstructure of water/ethanol mixtures, essentially decreasing the level of larger local clustering. Overall, the work presents a systematic characterization of ethanol effects on a widely used, common protein type, spidroins, and generalizes the findings to other intrinsically disordered proteins by pinpointing the general features of the response. The results can aid in designing effective alcohol treatments for proteins, but also enable design and tuning of protein material properties by a relatively controllable solvation handle, the addition of ethanol.

Conformational and Interaction Landscape of Histone H4 Tails in Nucleosomes Probed by Paramagnetic NMR Spectroscopy

The fundamental repeat unit of chromatin, the nucleosome, consists of approximately 147 base pairs of double-stranded DNA and a histone protein octamer containing two copies each of histones H2A, H2B, H3, and H4. Each histone possesses a dynamically disordered N-terminal tail domain, and it is well-established that the tails of histones H3 and H4 play key roles in chromatin compaction and regulation. Here we investigate the conformational ensemble and interactions of the H4 tail in nucleosomes by means of solution NMR measurements of paramagnetic relaxation enhancements (PREs) in recombinant samples reconstituted with 15N-enriched H4 and nitroxide spin-label tagged H3. The experimental PREs, which report on the proximities of individual H4 tail residues to the different H3 spin-label sites, are interpreted by using microsecond time-scale molecular dynamics simulations of the nucleosome core particle. Collectively, these data enable improved localization of histone H4 tails in nucleosomes and support the notion that H4 tails engage in a fuzzy complex interaction with nucleosomal DNA.

Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering

Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.

Intrinsic disorder and conformational coexistence in auxin coreceptors

AUXIN/INDOLE 3-ACETIC ACID (Aux/IAA) transcriptional repressor proteins and the TRANSPORT INHIBITOR RESISTANT 1/AUXIN SIGNALING F-BOX (TIR1/AFB) proteins to which they bind act as auxin coreceptors. While the structure of TIR1 has been solved, structural characterization of the regions of the Aux/IAA protein responsible for auxin perception has been complicated by their predicted disorder. Here, we use NMR, CD and molecular dynamics simulation to investigate the N-terminal domains of the Aux/IAA protein IAA17/AXR3. We show that despite the conformational flexibility of the region, a critical W-P bond in the core of the Aux/IAA degron motif occurs at a strikingly high (1:1) ratio of cis to trans isomers, consistent with the requirement of the cis conformer for the formation of the fully-docked receptor complex. We show that the N-terminal half of AXR3 is a mixture of multiple transiently structured conformations with a propensity for two predominant and distinct conformational subpopulations within the overall ensemble. These two states were modeled together with the C-terminal PB1 domain to provide the first complete simulation of an Aux/IAA. Using MD to recreate the assembly of each complex in the presence of auxin, both structural arrangements were shown to engage with the TIR1 receptor, and contact maps from the simulations match closely observations of NMR signal-decreases. Together, our results and approach provide a platform for exploring the functional significance of variation in the Aux/IAA coreceptor family and for understanding the role of intrinsic disorder in auxin signal transduction and other signaling systems.

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau

A known imbalance between intra-protein and protein-water interactions in many empirical force fields results in collapsed conformational ensembles of intrinsically disordered proteins in explicit solvent simulations that disagree with experiments. Multiple strategies have been introduced in the literature to modify protein-water interactions, which improve agreement between experiments and simulations. In this work, we combine simulations with standard and modified force fields with a spatially resolved analysis of solvation free energy contributions and compare the consequences of each strategy. We find that enhanced Lennard-Jones (LJ) interactions between protein atoms and water oxygens primarily improve the solvation of nonpolar functional groups of the protein. In contrast, modified electrostatics in the water model or strengthened LJ interactions between the protein and water hydrogens mainly affect the hydration of polar functional groups. Modified electrostatics further impact the average orientation of water molecules in the hydration shell. As a result, protein-water interactions with the first hydration layers are strengthened, while interactions with water molecules in higher hydration shells are weakened. Hence, distinct strategies to balance intra-protein and protein-water interactions in simulations have qualitatively different effects on protein solvation. These differences are not necessarily captured by comparisons to experiments that report on global parameters describing protein conformational ensembles, e.g., the radius of gyration, but will influence the tendency of a protein to form aggregates or phase-separated droplets.

Tuning the ATP-ATP and ATP-disordered protein interactions in high ATP concentration by altering water models

The adenosine triphosphate (ATP)-protein interactions have been of great interest since the recent experimental finding of ATP's role as a hydrotrope. The interaction between ATP and disordered proteins is fundamental to the dissolution of protein aggregates and the regulation of liquid-liquid phase separation by ATP. Molecular dynamics simulation is a powerful tool for analyzing these interactions in molecular detail but often suffers from inaccuracies in describing disordered proteins and ATPs in high concentrations. Recently, several water models have been proposed to improve the description of the protein-disordered states, yet how these models work with ATP has not been explored. To this end, here, we study how water models affect ATP and alter the ATP-ATP and ATP-protein interactions for the intrinsically disordered protein, α-Synuclein. Three water models, TIP4P-D, OPC, and TIP3P, are compared, while the protein force field is fixed to ff99SBildn. The results show that ATP over-aggregates into a single cluster in TIP3P water, but monomers and smaller clusters are found in TIP4P-D and OPC waters. ATP-protein interaction is also over-stabilized in TIP3P, whereas repeated binding/unbinding of ATP to α-Synuclein is observed in OPC and TIP4P-D waters, which is in line with the recent nuclear magnetic resonance experiment. The adenine ring-mediated interaction is found to play a major role in ATP-ATP and ATP-protein contacts. Interestingly, changing Mg2+ into Na+ strengthened the electrostatic interaction and promoted ATP oligomerization and ATP-α-Synuclein binding. Overall, this study shows that changing the water model can be an effective approach to improve the properties of ATP in high concentration.

Modeling the effects of phosphorylation on phase separation of the FUS low-complexity domain

Aggregation of the RNA-binding protein fused in sarcoma (FUS) is a hallmark of neurodegenerative diseases. Phosphorylation of Ser/Thr in the FUS low-complexity domain (FUS-LC) may regulate phase separation of FUS and prevent pathological aggregation in cells. However, many details of this process remain elusive to date. In this work, we systematically investigated the phosphorylation of FUS-LC and the underlying molecular mechanism by molecular dynamics (MD) simulations and free energy calculations. The results clearly show that phosphorylation can destroy the fibril core structure of FUS-LC by breaking interchain interactions, particularly contacts involving residues like Tyr, Ser, and Gln. Among the six phosphorylation sites, Ser61 and Ser84 may have more important effects on the stability of the fibril core. Our study reveals structural and dynamic details of FUS-LC phase separation modulated by phosphorylation.

Benchmarking Molecular Dynamics Force Fields for All-Atom Simulations of Biological Condensates

Proteins containing intrinsically disordered regions are integral parts of the cellular signaling pathways and common components of biological condensates. Point mutations in the protein sequence, genetic at birth or acquired through aging, can alter the properties of the condensates, marking the onset of neurodegenerative diseases such as ALS and dementia. While the all-atom molecular dynamics method can, in principle, elucidate the conformational changes that arise from point mutations, the applications of this method to protein condensate systems is conditioned upon the availability of molecular force fields that can accurately describe both structured and disordered regions of such proteins. Using the special-purpose Anton 2 supercomputer, we benchmarked the efficacy of nine presently available molecular force fields in describing the structure and dynamics of a Fused in sarcoma (FUS) protein. Five-microsecond simulations of the full-length FUS protein characterized the effect of the force field on the global conformation of the protein, self-interactions among its side chains, solvent accessible surface area, and the diffusion constant. Using the results of dynamic light scattering as a benchmark for the FUS radius of gyration, we identified several force fields that produced FUS conformations within the experimental range. Next, we used these force fields to perform ten-microsecond simulations of two structured RNA binding domains of FUS bound to their respective RNA targets, finding the choice of the force field to affect stability of the RNA-FUS complex. Taken together, our data suggest that a combination of protein and RNA force fields sharing a common four-point water model provides an optimal description of proteins containing both disordered and structured regions and RNA-protein interactions. To make simulations of such systems available beyond the Anton 2 machines, we describe and validate implementation of the best performing force fields in a publicly available molecular dynamics program NAMD. Our NAMD implementation enables simulations of large (tens of millions of atoms) biological condensate systems and makes such simulations accessible to a broader scientific community.

Accelerated Molecular Dynamics for Peptide Folding: Benchmarking Different Combinations of Force Fields and Explicit Solvent Models

Accelerated molecular dynamics (aMD) protocols were assessed on predicting the secondary structure of eight peptides, of which two are helical, three are β-hairpins, and three are disordered. Protocols consisted of combinations of three force fields (ff99SB, ff14SB, ff19SB) and two explicit solvation models (TIP3P and OPC), and were evaluated in two independent aMD simulations, one starting from an extended conformation, the other starting from a misfolded conformation. The results of these analyses indicate that all three combinations performed well on helical peptides. As for β-hairpins, ff19SB performed well with both solvation methods, with a slight preference for the TIP3P solvation model, even though performance was dependent on both peptide sequence and initial conformation. The ff19SB/OPC combination had the best performance on intrinsically disordered peptides. In general, ff14SB/TIP3P suffered the strongest helical bias.

Critical Assessment of Self-Consistency Checks in the All-Atom Molecular Dynamics Simulation of Intrinsically Disordered Proteins

All atom simulations can be used to quantify conformational properties of Intrinsically Disordered Proteins (IDP). However, simulations must satisfy convergence checks to ensure observables computed from simulation are reliable and reproducible. While absolute convergence is purely a theoretical concept requiring infinitely long simulation, a more practical, yet rigorous, approach is to impose Self Consistency Checks (SCCs) to gain confidence in the simulated data. Currently there is no study of SCCs in IDPs, unlike their folded counterparts. In this paper, we introduce different criteria for self-consistency checks for IDPs. Next, we impose these SCCs to critically assess the performance of different simulation protocols using the N terminal domain of HIV Integrase and the linker region of SARS-CoV-2 Nucleoprotein as two model IDPs. All simulation protocols begin with all-atom implicit solvent Monte Carlo (MC) simulation and subsequent clustering of MC generated conformations to create the representative structures of the IDPs. These representative structures serve as the initial structure for subsequent molecular dynamics (MD) runs with explicit solvent. We conclude that generating multiple short (∼3 μs) MD simulation trajectories─all starting from the most representative MC generated conformation─and merging them is the protocol of choice due to (i) its ability to satisfy multiple SCCs, (ii) consistently reproducing experimental data, and (iii) the efficiency of running independent trajectories in parallel by harnessing multiple cores available in modern GPU clusters. Running one long trajectory (greater than 20 μs) can also satisfy the first two criteria but is less desirable due to prohibitive computation time. These findings help resolve the challenge of identifying a usable starting configuration, provide an objective measure of SCC, and establish rigorous criteria to determine the minimum length (for one long simulation) or number of trajectories needed in all-atom simulation of IDPs.

Benchmarking Molecular Dynamics Force Fields for All-Atom Simulations of Biological Condensates

Proteins containing intrinsically disordered regions are integral components of the cellular signaling pathways and common components of biological condensates. Point mutations in the protein sequence, genetic at birth or acquired through aging, can alter the properties of the condensates, marking the onset of neurodegenerative diseases such as ALS and dementia. While the all-atom molecular dynamics method can, in principle, elucidate the conformational changes that arise from point mutations, the applications of this method to protein condensate systems is conditioned upon the availability of molecular force fields that can accurately describe both structured and disordered regions of such proteins. Using the special-purpose Anton 2 supercomputer, we benchmarked the efficacy of nine presently available molecular force fields in describing the structure and dynamics of a Fused in sarcoma (FUS) protein. Five-microsecond simulations of the full-length FUS protein characterized the effect of the force field on the global conformation of the protein, self-interactions among its side chains, solvent accessible surface area and the diffusion constant. Using the results of dynamic light scattering as a benchmark for the FUS radius of gyration, we identified several force fields that produced FUS conformations within the experimental range. Next, we used these force fields to perform ten-microsecond simulations of two structured RNA binding domains of FUS bound to their respective RNA targets, finding the choice of the force field to affect stability of the RNA-FUS complex. Taken together, our data suggest that a combination of protein and RNA force fields sharing a common four-point water model provides an optimal description of proteins containing both disordered and structured regions and RNA-protein interactions. To make simulations of such systems available beyond the Anton 2 machines, we describe and validate implementation of the best performing force fields in a publicly available molecular dynamics program NAMD. Our NAMD implementation enables simulations of large (tens of millions of atoms) biological condensate systems and makes such simulations accessible to a broader scientific community.

Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL

The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Incorporation of D2O-Induced Fluorine Chemical Shift Perturbations into Ensemble-Structure Characterization of the ERalpha Disordered Region

Structural characterization of intrinsically disordered proteins (IDPs) requires a concerted effort between experiments and computations by accounting for their conformational heterogeneity. Given the diversity of experimental tools providing local and global structural information, constructing an experimental restraint-satisfying structural ensemble remains challenging. Here, we use the disordered N-terminal domain (NTD) of the estrogen receptor alpha (ERalpha) as a model system to combine existing small-angle X-ray scattering (SAXS) and hydroxyl radical protein footprinting (HRPF) data and newly acquired solvent accessibility data via D2O-induced fluorine chemical shifting (DFCS) measurements. A new set of DFCS data for the solvent exposure of a set of 12 amino acid positions were added to complement previously acquired HRPF measurements for the solvent exposure of the other 16 nonoverlapping amino acids, thereby improving the NTD ensemble characterization considerably. We also found that while choosing an initial ensemble of structures generated from a different atomic-level force field or sampling/modeling method can lead to distinct contact maps even when the same sets of experimental measurements were used for ensemble-fitting, comparative analyses from these initial ensembles reveal commonly recurring structural features in their ensemble-averaged contact map. Specifically, nonlocal or long-range transient interactions were found consistently between the N-terminal segments and the central region, sufficient to mediate the conformational ensemble and regulate how the NTD interacts with its coactivator proteins.

Targeting disordered-structured domain interactions in Galectin-3 based on NMR and enhanced MD

Intrinsically disordered regions (IDRs) are common and important functional domains in many proteins. However, IDRs are difficult to target for drug development due to the lack of defined structures that would facilitate the identification of possible drug-binding pockets. Galectin-3 is a carbohydrate-binding protein of which overexpression has been implicated in a wide variety of disorders, including cancer and inflammation. Apart from its carbohydrate-recognition/binding domain (CRD), Galectin-3 also contains a functionally important disordered N-terminal domain (NTD) that contacts the C-terminal domain (CTD) and could be a target for drug development. To overcome challenges involved in inhibitor design due to lack of structure and the highly dynamic nature of the NTD, we used a protocol combining nuclear magnetic resonance data from recombinant Galectin-3 with accelerated molecular dynamics (MD) simulations. This approach identified a pocket in the CTD with which the NTD makes frequent contact. In accordance with this model, mutation of residues L131 and L203 in this pocket caused loss of Galectin-3 agglutination ability, signifying the functional relevance of the cavity. In silico screening was used to design candidate inhibitory peptides targeting the newly discovered cavity, and experimental testing of only three of these yielded one peptide that inhibits the agglutination promoted by wild-type Galectin-3. NMR experiments further confirmed that this peptide indeed binds to a cavity in the CTD, not within the actual CRD. Our results show that it is possible to apply a combination of MD simulations and NMR experiments to precisely predict the binding interface of a disordered domain with a structured domain, and furthermore use this predicted interface for designing inhibitors. This procedure can potentially be extended to many other targets in which similar IDR interactions play a vital functional role.

Fast Polarizable Water Model for Atomistic Simulations

Simulating water accurately has been a major challenge in atomistic simulations for decades. Inclusion of electronic polarizability effects holds considerable promise, yet existing approaches suffer from significant computational overheads compared to the widely used nonpolarizable water models. We have developed a globally optimal polarizable water model, OPC3-pol, that explicitly accounts for electronic polarizability with minimal impact on the computational efficiency. OPC3-pol reproduces five key bulk water properties at room temperature with an average relative error of 0.6%. In atomistic simulations, OPC3-pol's computational efficiency is in between that of 3- and 4-point nonpolarizable models; the model supports increased (4 fs) integration time step. OPC3-pol is tested in simulations of globular protein ubiquitin and a B-DNA dodecamer with several AMBER force fields, ff99SB, ff14SB, ff19SB, and OL15, demonstrating structure stability close to reference on multi-microsecond time scale. Simulation of an intrinsically disordered amyloid β-peptide yields an ensemble with the radius of gyration of a random coil. The proposed water model can be trivially adopted by any package that supports standard nonpolarizable force fields and water models; its intended use is in long classical atomistic simulations where water polarization effects are expected to be important.

The performance of OPC water model in prediction of the phase equilibria of methane hydrate

Molecular dynamics (MD) simulations were performed to determine the three-phase coexistence line of sI methane hydrates. The MD simulations were carried out at four different pressures (4, 10, 40 and 100 MPa) by using direct phase coexistence method. In current simulations, water was described by either TIP4P/Ice or OPC models and methane was described as a simple Lennard-Jones (LJ) interaction site. Lorentz-Berthelot combining rules were used to calculate the parameters of the cross interactions. For OPC model, positive deviations from the energetic Lorentz-Berthelot rule were also considered based on the solubility of methane in water. For TIP4P/Ice water model, the obtained three phase coexistence temperatures showed good agreement with experiment data at higher pressures, which is consistent with previous predictions. For OPC water model, simulations using the classic and the modified LB parameters both showed negative deviations to the experimental values. Our results also indicated that the deviation of the T3 prediction by OPC model not much correlated with the predicted melting point of ice. At 4 MPa, the modified OPC model showed outstanding prediction of hydrate equilibrium temperature, even better than the prediction by TIP4P/Ice. The relative higher accuracy in biomolecular MD of OPC model suggests that this model may have a better performance in hydrate MD simulations of biomolecule-based additives.