Re-Balancing Replica Exchange with Solute Tempering for Sampling Dynamic Protein Conformations

Likely Overstabilization of Charge-Charge Interactions in CHARMM36m(w): A Case for a99SB-disp Water

Recent years have witnessed drastic improvements in general-purpose explicit solvent protein force fields, partially driven by the need to study intrinsically disordered proteins (IDPs), and yet the state-of-the-art force fields such as CHARMM36m (c36m) and a99SB-disp still provide different performances in simulating disordered protein states, where c36m has a bias toward overcompaction for large IDPs. Here, we examine the performance of c36m and a99SB-disp in describing the stabilities of a set of 46 amino acid backbone and side chain pairs in various configurations. The free energy results show that c36m systematically predicts stronger interactions compared to a99SB-disp by an average of 0.2 kcal/mol for nonpolar pairs, 0.6 kcal/mol for polar pairs, and 0.8 kcal/mol for salt bridges. The most severe overstabilization in c36m is observed for charged pairs involving the Arg and Glu side chains by up to 2.9 kcal/mol. Importantly, the systematic overstabilization of c36m is only marginally alleviated by c36mw, an ad hoc patch to c36m that increases the dispersion interactions between TIP3P hydrogens and protein atoms. Guided by free energy decomposition, we evaluated if revising the charges alone could alleviate the severe overstabilization of salt bridges of c36m(w) vs a99SB-disp. The results suggested that the direct modification of protein-water interactions is also necessary. Toward this end, we proposed a tentative modification to c36m, referred to as c36mrb-disp, which combines modified Arg side chain charges, retuned backbone hydrogen bonding strength, and the a99SB-disp water model. The modified force field successfully reproduces the secondary structures of several intrinsically disordered peptides and proteins including (AAQAA)3, GB1p, and p53 transactivation domain, while maintaining the stability of a set of folded proteins. This work provides a set of useful systems for benchmarking and optimizing protein force fields and highlights the importance of balancing protein-protein and protein-water electrostatic interactions for accurately describing both folded and disordered proteins.

Validating Small-Molecule Force Fields for Macrocyclic Compounds Using NMR Data in Different Solvents

Macrocycles are a promising class of compounds as therapeutics for difficult drug targets due to a favorable combination of properties: They often exhibit improved binding affinity compared to their linear counterparts due to their reduced conformational flexibility, while still being able to adapt to environments of different polarity. To assist in the rational design of macrocyclic drugs, there is need for computational methods that can accurately predict conformational ensembles of macrocycles in different environments. Molecular dynamics (MD) simulations remain one of the most accurate methods to predict ensembles quantitatively, although the accuracy is governed by the underlying force field. In this work, we benchmark four different force fields for their application to macrocycles by performing replica exchange with solute tempering (REST2) simulations of 11 macrocyclic compounds and comparing the obtained conformational ensembles to nuclear Overhauser effect (NOE) upper distance bounds from NMR experiments. Especially, the modern force fields OpenFF 2.0 and XFF yield good results, outperforming force fields like GAFF2 and OPLS/AA. We conclude that REST2 in combination with modern force fields can often produce accurate ensembles of macrocyclic compounds. However, we also highlight examples for which all examined force fields fail to produce ensembles that fulfill the experimental constraints.

Illuminating Protein Allostery by Chemically Accurate Contact Response Analysis (ChACRA)

Decoding allostery at the atomic level is essential for understanding the relationship between a protein's sequence, structure, and dynamics. Recently, we have shown that decomposing temperature responses of inter-residue contacts can reveal allosteric couplings and provide useful insight into the functional dynamics of proteins. The details of this Chemically Accurate Contact Response Analysis (ChACRA) are presented here along with its application to two well-known allosteric proteins. The first protein, IGPS, is a model of ensemble allostery that lacks clear structural differences between the active and inactive states. We show that the application of ChACRA reveals the experimentally identified allosteric coupling between effector and active sites of IGPS. The second protein, ATCase, is a classic example of allostery with distinct active and inactive structural states. Using ChACRA, we directly identify the most significant residue level interactions underlying the enzyme's cooperative behavior. Both test cases demonstrate the utility of ChACRA's unsupervised machine learning approach for dissecting allostery at the residue level.

Convergence-Adaptive Roundtrip Method Enables Rapid and Accurate FEP Calculations

The free energy perturbation (FEP) method is a powerful technique for accurate binding free energy calculations, which is crucial for identifying potent ligands with a high affinity in drug discovery. However, the widespread application of FEP is limited by the high computational cost required to achieve equilibrium sampling and the challenges in obtaining converged predictions. In this study, we present the convergence-adaptive roundtrip (CAR) method, which is an enhanced adaptive sampling approach, to address the key challenges in FEP calculations, including the precision-efficiency tradeoff, sampling efficiency, and convergence assessment. By employing on-the-fly convergence analysis to automatically adjust simulation times, enabling efficient traversal of the important phase space through rapid propagation of conformations between different states and eliminating the need for multiple parallel simulations, the CAR method increases convergence and minimizes computational overhead while maintaining calculation accuracy. The performance of the CAR method was evaluated through relative binding free energy (RBFE) calculations on benchmarks comprising four diverse protein-ligand systems. The results demonstrated a significant speedup of over 8-fold compared to conventional FEP methods while maintaining high accuracy. The overall R2 values of 0.65 and 0.56 were obtained using the combined-structure FEP approach and the single-step FEP approach, respectively, in conjunction with the CAR method. In-depth case studies further highlighted the superior performance of the CAR method in terms of convergence acceleration, improved predicted correlations, and reduced computational costs. The advancement of the CAR method makes it a highly effective approach, enhancing the applicability of FEP in drug discovery.

Biomolecular dynamics in the 21st century

The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.

Molecular simulations integrated with experiments for probing the interaction dynamics and binding mechanisms of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) exploit their plasticity to deploy a rich panoply of soft interactions and binding phenomena. Advances in tailoring molecular simulations for IDPs combined with experimental cross-validation offer an atomistic view of the mechanisms that control IDP binding, function, and dysfunction. The emerging theme is that unbound IDPs autonomously form transient local structures and self-interactions that determine their binding behavior. Recent results have shed light on whether and how IDPs fold, stay disordered or drive condensation upon binding; how they achieve binding specificity and select among competing partners. The disorder-binding paradigm is now being proactively used by researchers to target IDPs for rational drug design and engineer molecular responsive elements for biosensing applications.

Improved Prediction of Molecular Response to Pulling by Combining Force Tempering with Replica Exchange Methods

Small mechanical forces play important functional roles in many crucial cellular processes, including in the dynamic behavior of the cytoskeleton and in the regulation of osmotic pressure through membrane-bound proteins. Molecular simulations offer the promise of being able to design the behavior of proteins that sense and respond to these forces. However, it is difficult to predict and identify the effect of the relevant piconewton (pN) scale forces due to their small magnitude. Previously, we introduced the Infinite Switch Simulated Tempering in Force (FISST) method, which allows one to estimate the effect of a range of applied forces from a single molecular dynamics simulation, and also demonstrated that FISST additionally accelerates sampling of a molecule's conformational landscape. For some problems, we find that this acceleration is not sufficient to capture all relevant conformational fluctuations, and hence, here we demonstrate that FISST can be combined with either temperature replica exchange or solute tempering approaches to produce a hybrid method that enables more robust prediction of the effect of small forces on molecular systems.

Toward Accurate Simulation of Coupling between Protein Secondary Structure and Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate phase separation that underlies the formation of a biomolecular condensate. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding the sequence-specific phase separation of IDPs. However, the widely used Cα-only models are limited in capturing the peptide nature of IDPs, particularly backbone-mediated interactions and effects of secondary structures, in phase separation. Here, we describe a hybrid resolution (HyRes) protein model toward a more accurate description of the backbone and transient secondary structures in phase separation. With an atomistic backbone and coarse-grained side chains, HyRes can semiquantitatively capture the residue helical propensity and overall chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for the direct simulation of spontaneous phase separation and, at the same time, appears accurate enough to resolve the effects of single His to Lys mutations. HyRes simulations also successfully predict increased β-structure formation in the condensate, consistent with available experimental CD data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate the phase separation propensity as measured by the saturation concentration. The simulations successfully recapitulate the effect of these mutants on the helicity and phase separation propensity of TDP-43 CR. Analyses reveal that the balance between backbone and side chain-mediated interactions, but not helicity itself, actually determines phase separation propensity. These results support that HyRes represents an effective protein model for molecular simulation of IDP phase separation and will help to elucidate the coupling between transient secondary structures and phase separation.

Replica Exchange with Hybrid Tempering Efficiently Samples PGLa Peptide Binding to Anionic Bilayer

We evaluated the utility of a variant of the replica exchange method, a replica exchange with hybrid tempering (REHT), for all-atom explicit water biomolecular simulations and compared it with a more traditional replica exchange with the solute tempering (REST) algorithm. As a test system, we selected a 21-mer antimicrobial peptide PGLa binding to an anionic DMPC/DMPG lipid bilayer. Application of REHT revealed the following binding mechanism. Due to the strong hydrophobic moment, the bound PGLa adopts an extensive helical structure. The binding free energy landscape identifies two major bound states, a metastable surface bound state and a dominant inserted state. In both states, positively charged PGLa amino acids maintain electrostatic interactions with anionic phosphate groups by rotating the PGLa helix around its axis. PGLa binding causes an influx of anionic DMPG and an efflux of zwitterionic DMPC lipids from the peptide proximity. PGLa thins the bilayer and disorders the adjacent fatty acid tails. Deep invasion of water wires into the bilayer hydrophobic core is detected in the inserted peptide state. The analysis of charge density distributions indicated that peptide positive charges are nearly compensated for by lipid negative charges and water dipole ordering, whereas ions play no role in peptide binding. Thus, electrostatic interactions are the key energetic factor in binding cationic PGLa to an anionic DMPC/DMPG bilayer. Comparison of REHT and REST shows that due to exclusion of lipids from tempered partition, REST lags behind REHT in peptide equilibration, particularly, with respect to peptide insertion and helix acquisition. As a result, REST struggles to provide accurate details of PGLa binding, although it still qualitatively maps the bimodal binding mechanism. Importantly, REHT not only equilibrates PGLa in the bilayer faster than REST, but also with less computational effort. We conclude that REHT is a preferable choice for studying interfacial biomolecular systems.

Accurate Simulation of Coupling between Protein Secondary Structure and Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate liquid-liquid phase separation (LLPS) that underlies the formation of membraneless organelles. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding sequence-specific phase separation of IDPs. However, the widely-used Cα-only models are severely limited in capturing the peptide nature of IDPs, including backbone-mediated interactions and effects of secondary structures, in LLPS. Here, we describe a hybrid resolution (HyRes) protein model for accurate description of the backbone and transient secondary structures in LLPS. With an atomistic backbone and coarse-grained side chains, HyRes accurately predicts the residue helical propensity and chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for direct simulation of spontaneous phase separation, and at the same time accurate enough to resolve the effects of single mutations. HyRes simulations also successfully predict increased beta-sheet formation in the condensate, consistent with available experimental data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate LLPS propensity. The simulations successfully recapitulate the effect of these mutants on the helicity and LLPS propensity of TDP-43 CR. Analyses reveal that the balance between backbone and sidechain-mediated interactions, but not helicity itself, actually determines LLPS propensity. We believe that the HyRes model represents an important advance in the molecular simulation of LLPS and will help elucidate the coupling between IDP transient secondary structures and phase separation.