Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

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.

Energy, water, and protein folding: A molecular dynamics-based quantitative inventory of molecular interactions and forces that make proteins stable

Protein folding energetics can be determined experimentally on a case-by-case basis but it is not understood in sufficient detail to provide deep control in protein design. The fundamentals of protein stability have been outlined by calorimetry, protein engineering, and biophysical modeling, but these approaches still face great difficulty in elucidating the specific contributions of the intervening molecules and physical interactions. Recently, we have shown that the enthalpy and heat capacity changes associated to the protein folding reaction can be calculated within experimental error using molecular dynamics simulations of native protein structures and their corresponding unfolded ensembles. Analyzing in depth molecular dynamics simulations of four model proteins (CI2, barnase, SNase, and apoflavodoxin), we dissect here the energy contributions to ΔH (a key component of protein stability) made by the molecular players (polypeptide and solvent molecules) and physical interactions (electrostatic, van der Waals, and bonded) involved. Although the proteins analyzed differ in length, isoelectric point and fold class, their folding energetics is governed by the same quantitative pattern. Relative to the unfolded ensemble, the native conformations are enthalpically stabilized by comparable contributions from protein-protein and solvent-solvent interactions, and almost equally destabilized by interactions between protein and solvent molecules. The native protein surface seems to interact better with water than the unfolded one, but this is outweighed by the unfolded surface being larger. From the perspective of physical interactions, the native conformations are stabilized by van de Waals and Coulomb interactions and destabilized by conformational strain arising from bonded interactions. Also common to the four proteins, the sign of the heat capacity change is set by interactions between protein and solvent molecules or, from the alternative perspective, by Coulomb interactions.

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.

Limitations of Atomistic Molecular Dynamics to Reveal Ejection of Proteins from Charged Nanodroplets

Atomistic molecular dynamics (MD) is frequently used to unravel the mechanisms of macroion release from electrosprayed droplets. However, atomistic MD is currently feasible for only the smallest window of droplet sizes appearing at the end steps of a droplet's lifetime. The relevance of the observations made to the actual droplet evolution, which is much longer than the simulated sizes, has not been addressed yet in the literature. Here, we perform a systematic study of the desolvation mechanisms of poly(ethylene glycol) (PEG), protonated peptides of different compositions, and proteins, to (a) obtain insight into the charging mechanism of macromolecules in larger droplets than those that are currently amenable to atomistic MD and (b) examine whether currently used atomistic MD modeling can establish the extrusion mechanism of proteins from droplets. To mimic larger droplets that are not amenable to MD modeling, we scale down the systems, by simulating a large droplet size relative to the macromolecule. MD of PEG charging reveals that, above a critical droplet size, ions are available near the backbone of the macromolecule, but charging occurs only transiently by transfer of ions from the solvent to the macroion, while below the critical size, the capture of the ion from PEG has a lifetime sufficiently long for the extrusion of a charged PEG from the aqueous droplet. This is the first report of the role of droplet curvature in the relation between macroion conformation and charging. Simulations of protonated peptides with a high degree of hydrophobicity show that partial extrusion of a peptide from the droplet surface is rare relative to desolvation by drying-out. Different from what has been presented in the literature, we argue that atomistic MD simulations have not sufficiently established the extrusion mechanism of proteins from droplets and their charging mechanism. We also argue that release of highly charged proteins can occur at an earlier stage of a droplet's lifetime than predicted by atomistic MD. In this earlier stage, we emphasize the key role of jets emanating from a droplet at the point of charge-induced instability in the release of proteins.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as ₂₁₆PTPP₂₁₉, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

Fuzzy Drug Targets: Disordered Proteins in the Drug-Discovery Realm

Intrinsically disordered proteins (IDPs) and regions (IDRs) form a large part of the eukaryotic proteome. Contrary to the structure-function paradigm, the disordered proteins perform a myriad of functions in vivo. Consequently, they are involved in various disease pathways and are plausible drug targets. Unlike folded proteins, that have a defined structure and well carved out drug-binding pockets that can guide lead molecule selection, the disordered proteins require alternative drug-development methodologies that are based on an acceptable picture of their conformational ensemble. In this review, we discuss various experimental and computational techniques that contribute toward understanding IDP "structure" and describe representative pursuances toward IDP-targeting drug development. We also discuss ideas on developing rational drug design protocols targeting IDPs.

Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach

This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.

Convergent views on disordered protein dynamics from NMR and computational approaches

Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.

Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM

We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the α_R region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.

Quantitative prediction of ensemble dynamics, shapes and contact propensities of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are highly dynamic systems that play an important role in cell signaling processes and their misfunction often causes human disease. Proper understanding of IDP function not only requires the realistic characterization of their three-dimensional conformational ensembles at atomic-level resolution but also of the time scales of interconversion between their conformational substates. Large sets of experimental data are often used in combination with molecular modeling to restrain or bias models to improve agreement with experiment. It is shown here for the N-terminal transactivation domain of p53 (p53TAD) and Pup, which are two IDPs that fold upon binding to their targets, how the latest advancements in molecular dynamics (MD) simulations methodology produces native conformational ensembles by combining replica exchange with series of microsecond MD simulations. They closely reproduce experimental data at the global conformational ensemble level, in terms of the distribution properties of the radius of gyration tensor, and at the local level, in terms of NMR properties including 15N spin relaxation, without the need for reweighting. Further inspection revealed that 10-20% of the individual MD trajectories display the formation of secondary structures not observed in the experimental NMR data. The IDP ensembles were analyzed by graph theory to identify dominant inter-residue contact clusters and characteristic amino-acid contact propensities. These findings indicate that modern MD force fields with residue-specific backbone potentials can produce highly realistic IDP ensembles sampling a hierarchy of nano- and picosecond time scales providing new insights into their biological function.

Challenges in describing the conformation and dynamics of proteins with ambiguous behavior

Traditionally, our understanding of how proteins operate and how evolution shapes them is based on two main data sources: the overall protein fold and the protein amino acid sequence. However, a significant part of the proteome shows highly dynamic and/or structurally ambiguous behavior, which cannot be correctly represented by the traditional fixed set of static coordinates. Representing such protein behaviors remains challenging and necessarily involves a complex interpretation of conformational states, including probabilistic descriptions. Relating protein dynamics and multiple conformations to their function as well as their physiological context (e.g., post-translational modifications and subcellular localization), therefore, remains elusive for much of the proteome, with studies to investigate the effect of protein dynamics relying heavily on computational models. We here investigate the possibility of delineating three classes of protein conformational behavior: order, disorder, and ambiguity. These definitions are explored based on three different datasets, using interpretable machine learning from a set of features, from AlphaFold2 to sequence-based predictions, to understand the overlap and differences between these datasets. This forms the basis for a discussion on the current limitations in describing the behavior of dynamic and ambiguous proteins.

Energetics of a protein disorder-order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered ‘lid’ region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Molecular simulations and biophysical measurements elucidate why the ligand AM-7209 orders a disordered region of the protein MDM2 on binding. This work expands strategies available to medicinal chemists for targeting disordered proteins.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

A unifying framework for amyloid-mediated membrane damage: The lipid-chaperone hypothesis

Over the past thirty years, researchers have highlighted the role played by a class of proteins or polypeptides that forms pathogenic amyloid aggregates in vivo, including i) the amyloid Aβ peptide, which is known to form senile plaques in Alzheimer's disease; ii) α-synuclein, responsible for Lewy body formation in Parkinson's disease and iii) IAPP, which is the protein component of type 2 diabetes-associated islet amyloids. These proteins, known as intrinsically disordered proteins (IDPs), are present as highly dynamic conformational ensembles. IDPs can partially (mis) fold into (dys) functional conformations and accumulate as amyloid aggregates upon interaction with other cytosolic partners such as proteins or lipid membranes. In addition, an increasing number of reports link the toxicity of amyloid proteins to their harmful effects on membrane integrity. Still, the molecular mechanism underlying the amyloidogenic proteins transfer from the aqueous environment to the hydrocarbon core of the membrane is poorly understood. This review starts with a historical overview of the toxicity models of amyloidogenic proteins to contextualize the more recent lipid-chaperone hypothesis. Then, we report the early molecular-level events in the aggregation and ion-channel pore formation of Aβ, IAPP, and α-synuclein interacting with model membranes, emphasizing the complexity of these processes due to their different spatial-temporal resolutions. Next, we underline the need for a combined experimental and computational approach, focusing on the strengths and weaknesses of the most commonly used techniques. Finally, the last two chapters highlight the crucial role of lipid-protein complexes as molecular switches among ion-channel-like formation, detergent-like, and fibril formation mechanisms and their implication in fighting amyloidogenic diseases.

Investigating the Disordered and Membrane-Active Peptide A-Cage-C Using Conformational Ensembles

The driving forces and conformational pathways leading to amphitropic protein-membrane binding and in some cases also to protein misfolding and aggregation is the subject of intensive research. In this study, a chimeric polypeptide, A-Cage-C, derived from α-Lactalbumin is investigated with the aim of elucidating conformational changes promoting interaction with bilayers. From previous studies, it is known that A-Cage-C causes membrane leakages associated with the sporadic formation of amorphous aggregates on solid-supported bilayers. Here we express and purify double-labelled A-Cage-C and prepare partially deuterated bicelles as a membrane mimicking system. We investigate A-Cage-C in the presence and absence of these bicelles at non-binding (pH 7.0) and binding (pH 4.5) conditions. Using in silico analyses, NMR, conformational clustering, and Molecular Dynamics, we provide tentative insights into the conformations of bound and unbound A-Cage-C. The conformation of each state is dynamic and samples a large amount of overlapping conformational space. We identify one of the clusters as likely representing the binding conformation and conclude tentatively that the unfolding around the central W23 segment and its reorientation may be necessary for full intercalation at binding conditions (pH 4.5). We also see evidence for an overall elongation of A-Cage-C in the presence of model bilayers.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Applications of KBFF20

Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native β-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20

A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.

Histone H4 Tails in Nucleosomes: a Fuzzy Interaction with DNA

The interaction of positively charged N-terminal histone tails with nucleosomal DNA plays an important role in chromatin assembly and regulation, modulating their susceptibility to post-translational modifications and recognition by chromatin-binding proteins. Here, we report residue-specific ¹⁵ N NMR relaxation rates for histone H4 tails in reconstituted nucleosomes. These data indicate that H4 tails are strongly dynamically disordered, albeit with reduced conformational flexibility compared to a free peptide with the same sequence. Remarkably, the NMR observables were successfully reproduced in a 2-μs MD trajectory of the nucleosome. This is an important step toward resolving an apparent inconsistency where prior simulations were generally at odds with experimental evidence on conformational dynamics of histone tails. Our findings indicate that histone H4 tails engage in a fuzzy interaction with nucleosomal DNA, underpinned by a variable pattern of short-lived salt bridges and hydrogen bonds, which persists at low ionic strength (0-100 mM NaCl).

Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.

Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins

We outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships-however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.