Residue-Specific Force Field Improving the Sample of Intrinsically Disordered Proteins and Folded Proteins

Refined Bonded Terms in Coarse-Grained Models for Intrinsically Disordered Proteins Improve Backbone Conformations

Coarse-grained models designed for intrinsically disordered proteins and regions (IDP/Rs) usually omit some bonded potentials (e.g., angular and dihedral potentials) as a conventional strategy to enhance backbone flexibility. However, a notable drawback of this approach is the generation of inaccurate backbone conformations. Here, we addressed this problem by introducing residue-specific angular, refined dihedral, and correction map (CMAP) potentials, derived based on the statistics from a customized coil database. These bonded potentials were integrated into the existing Mpipi model, resulting in a new model, denoted as the "Mpipi+" model. Results show that the Mpipi+ model can improve backbone conformations. More importantly, it can markedly improve the secondary structure propensity (SSP) based on the experimental chemical shift and, consequently, succeed in capturing transient secondary structures. Moreover, the Mpipi+ model preserves the liquid-liquid phase separation (LLPS) propensities of IDPs.

Towards bone regeneration: Understanding the nucleating ability of proline-rich peptides in biomineralisation

Obtaining rapid mineralisation is a challenge in current bone graft materials, which has been attributed to the difficulty of guiding the biological processes towards osteogenesis. Amelogenin, a key protein in enamel formation, inspired the design of two intrinsically disordered peptides (P2 and P6) that enhance in vivo bone formation, but the process is not fully understood. In this study, we have elucidated the mechanism by which these peptides induce improved mineralisation. Our molecular dynamics analysis demonstrated that in an aqueous environment, P2 and P6 fold to interact with the surrounding Ca2+, PO43- and OH- ions, which can lead to apatite nucleation. Although P2 has a less stable backbone, it folds to a stable structure that allows for the nucleation of larger calcium phosphate aggregates than P6. These results were validated experimentally in a concentrated simulated body fluid solution, where the peptide solutions accelerated the mineralisation process compared to the control and yielded mineral structures mimicking the amorphous calcium phosphate crystals that can be found in lamella bone. A pH drop for the peptide groups suggests depletion of calcium and phosphate, a prerequisite for intrinsic osteoinduction, while S/TEM and SEM suggested that the peptide regulated the mineral nucleation into lamella flakes. Evidently, the peptides accelerate and guide mineral formation, elucidating the mechanism for how these peptides can improve the efficacy of P2 or P6 containing devices for bone regeneration. The work also demonstrates how experimental mineralisation study coupled with molecular dynamics is a valid method for understanding and predicting in vivo performance prior to animal trials.

Multiscale Computational Framework for the Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

The reversible assembly of intrinsically disordered proteins (IDPs) to form membraneless organelles (MLOs) is a fundamental process involved in the spatiotemporal regulation in living cells. MLOs formed via liquid-liquid phase separation (LLPS) serve as molecule-enhancing hubs to regulate cell functions. Owing to the complexity and dynamic nature of the protein assembly via a network of weak inter- and intra-molecular interactions, it is challenging to describe and predict the LLPS behavior. We have developed a multiscale computational model for IDPs, using the fused in sarcoma (FUS) protein and its variants as illustrative examples. To simplify the description of protein, FUS is represented as a linear chain of stickers interspaced by spacers, as inspired by the associative polymer theory. Low-complexity aromatic-rich kinked segments (LARKS) available in FUS were identified using LARKSdb and represented as "stickers". The pairwise potential energies of each pair of stickers and their β-sheet-forming propensity were estimated via molecular docking and all atomistic molecular dynamics (AA-MD) simulations. Subsequently, FUS chains were randomly positioned in a cubic lattice as coarse-grained (CG) beads, with the bead assignment based on the Kuhn length estimation of stickers and spacers. Stochastic FUS movements were modeled by Monte Carlo (MC) simulations. In addition to the Metropolis algorithm, discretized pair potential distributions between stickers were considered in the move acceptance criteria. The chosen pair potential represents one of the possible binding energy states, with its probability determined by the frequency of the binding energy distribution histogram. The fluctuations of averaged radial distribution functions (RDFs) in successive MC trial move intervals of equilibrated lattice MC simulations were used to indicate the dynamic nature of assembly/disassembly of the protein chains. This multiscale computational framework provides an economical and efficient way of predicting and describing the LLPS behavior of IDPs.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

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.

Balanced Three-Point Water Model OPC3-B for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) play a critical role in many biological processes. Due to the inherent structural flexibility of IDPs, experimental methods present significant challenges for sampling their conformational information at the atomic level. Therefore, molecular dynamics (MD) simulations have emerged as the primary tools for modeling IDPs whose accuracy depend on force field and water model. To enhance the accuracy of physical modeling of IDPs, several force fields have been developed. However, current water models lack precision and underestimate the interaction between water molecules and proteins. Here, we used Monte-Carlo re-weighting method to re-parameterize a three-point water model based on OPC3 for IDPs (named OPC3-B). We benchmarked the performance of OPC3-B compared with nine different water models for 10 IDPs and three ordered proteins. The results indicate that the performance of OPC3-B is better than other water models for both IDPs and ordered proteins. At the same time, OPC3-B possess the power of transferability with other force field to simulate IDPs. This newly developed water model can be used to insight into the research of sequence-disordered-function paradigm for IDPs.

Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

Phosphorylation of proteins plays an important regulatory role at almost all levels of cellular organization. Molecular dynamics (MD) simulation is a promising tool to reveal the mechanism of how phosphorylation regulates many key biological processes at the atomistic level. MD simulation accuracy depends on force field precision, while the current force fields for phospho-amino acids have resulted in notable inconsistency with experimental data. Here, a new force field parameter (named FB18CMAP) is generated by fitting against quantum mechanics (QM) energy in aqueous solution with φ/ψ dihedral potential-energy surfaces optimized using CMAP parameters. MD simulations of phosphorylated dipeptides, intrinsically disordered proteins (IDPs), and ordered (folded) proteins show that FB18CMAP can mimic NMR observables and structural characteristics of phosphorylated dipeptides and proteins more accurately than the FB18 force field. These findings suggest that FB18CMAP performs well in both the simulation of ordered and disordered states of phosphorylated proteins.

Influence of force field choice on the conformational landscape of rat and human islet amyloid polypeptide

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein (IDP) whose abnormal aggregation into toxic soluble oligomers and insoluble amyloid fibrils is a pathological feature in type-2 diabetes. Rat IAPP (rIAPP) differs from hIAPP by only six amino acids yet has a reduced tendency to aggregate or form fibrils. The structures of the monomeric forms of IAPP are difficult to characterize due to their intrinsically disordered nature. Molecular dynamics simulations can provide a detailed characterization of the monomeric forms of rIAPP and hIAPP in near-physiological conditions. In this work, the conformational landscapes of rIAPP and hIAPP as a function of secondary structure content were predicted using well-tempered bias exchange metadynamics simulations. Several combinations of commonly used biomolecular force fields and water models were tested. The predicted conformational preferences of both rIAPP and hIAPP are typical of IDPs, exhibiting dominant random coil structures but showing a low propensity for transient α-helical conformations. Predicted nuclear magnetic resonance Cα chemical shifts reveal different preferences with each force field towards certain conformations, with AMBERff99SBnmr2/TIP4Pd showing the best agreement with the experiment. Comparisons of secondary structure content demonstrate residue-specific differences between hIAPP and rIAPP that may reflect their different aggregation propensities.

Personal Precise Force Field for Intrinsically Disordered and Ordered Proteins Based on Deep Learning

Intrinsically disordered proteins (IDPs) are proteins without a fixed three-dimensional (3D) structure under physiological conditions and are associated with Parkinson's disease, Alzheimer's disease, cancer, cardiovascular disease, amyloidosis, diabetes, and other diseases. Experimental methods can hardly capture the ensemble of diverse conformations for IDPs. Molecular dynamics (MD) simulations can sample continuous conformations that might provide a valuable complement to experimental data. However, the accuracy of MD simulations depends on the quality of force field. In particular, the evolutionary conservation and coevolution of IDPs introduce that current force fields could not precisely reproduce the conformation of IDPs. In order to improve the performance of force field, deep learning and reweighting methods were used to automatically generate personal force field parameters for intrinsically disordered and ordered proteins. At first, the deep learning method predicted more accuracy φ/ψ dihedral of residue than the previous method. Then, reweighting optimized the personal force field parameters for each residue. Finally, typical representative systems such as IDPs, structure protein, and fast-folding protein were used to evaluate this force field. The results indicate that two personal force field parameters (named PPFF1 and PPFF1_af2) could better reproduce the experimental observables than ff03CMAP force field. In summary, this strategy will provide feasibility for the development of precise personal force fields.

Modeling structural interconversion in Alzheimers' amyloid beta peptide with classical and intrinsically disordered protein force fields

A comprehensive understanding of the aggregation mechanism in amyloid beta 42 (Aβ42) peptide is imperative for developing therapeutic drugs to prevent or treat Alzheimer's disease. Because of the high flexibility and lack of native tertiary structures of Aβ42, molecular dynamics (MD) simulations may help elucidate the peptide's dynamics with atomic details and collectively improve ensembles not seen in experiments. We applied microsecond-timescale MD simulations to investigate the dynamics and conformational changes of Aβ42 by using a newly developed Amber force field (ff14IDPSFF). We compared the ff14IDPSFF and the regular ff14SB force field by examining the conformational changes of two distinct Aβ42 monomers in explicit solvent. Conformational ensembles obtained by simulations depend on the force field and initial structure, Aβ42^α-helix or Aβ42^β-strand. The ff14IDPSFF sampled a high ratio of disordered structures and diverse β-strand secondary structures; in contrast, ff14SB favored helicity during the Aβ42^α-helix simulations. The conformations obtained from Aβ42^β-strand simulations maintained a balanced content in the disordered and helical structures when simulated by ff14SB, but the conformers clearly favored disordered and β-sheet structures simulated by ff14IDPSFF. The results obtained with ff14IDPSFF qualitatively reproduced the NMR chemical shifts well. In-depth peptide and cluster analysis revealed some characteristic features that may be linked to early onset of the fibril-like structure. The C-terminal region (mainly M35-V40) featured in-registered anti-parallel β-strand (β-hairpin) conformations with tested systems. Our work should expand the knowledge of force field and structure dependency in MD simulations and reveals the underlying structural mechanism-function relationship in Aβ42 peptides. Communicated by Ramaswamy H. Sarma.

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

Polarizable Force Field of Intrinsically Disordered Proteins with CMAP and Reweighting Optimization

Intrinsically disordered proteins (IDPs) are highly structurally heterogeneous without a specific tertiary structure under physiology conditions and play key roles in the development of human diseases. Due to the characteristics of diverse conformations, as one of the important methods, molecular dynamics simulation can complement information for experimental methods. Because of the enrichment for charged amino acids for IDPs, polarizable force fields should be a good choice for the simulation of IDPs. However, current polarizable force fields are limited in sampling conformer features of IDPs. Therefore, a polarizable force field was released and named Drude2019IDP based on Drude2019 with reweighting and grid-based potential energy correction map optimization. In order to evaluate the performance of Drude2019IDP, 16 dipeptides, 18 short peptides, 3 representative IDPs, and 5 structural proteins were simulated. The results show that the NMR observables driven by Drude2019IDP are in better agreement with the experiment data than those by Drude2019 on short peptides and IDPs. Drude2019IDP can sample more diverse conformations than Drude2019. Furthermore, the performances of the two force fields are similar to the sample ordered proteins. These results confirm that the developed Drude2019IDP can improve the reproduction of conformers for intrinsically disordered proteins and can be used to gain insight into the paradigm of sequence-disorder for IDPs.

Balanced Force Field ff03CMAP Improving the Dynamics Conformation Sampling of Phosphorylation Site

Phosphorylation plays a key role in plant biology, such as the accumulation of plant cells to form the observed proteome. Statistical analysis found that many phosphorylation sites are located in disordered regions. However, current force fields are mainly trained for structural proteins, which might not have the capacity to perfectly capture the dynamic conformation of the phosphorylated proteins. Therefore, we evaluated the performance of ff03CMAP, a balanced force field between structural and disordered proteins, for the sampling of the phosphorylated proteins. The test results of 11 different phosphorylated systems, including dipeptides, disordered proteins, folded proteins, and their complex, indicate that the ff03CMAP force field can better sample the conformations of phosphorylation sites for disordered proteins and disordered regions than ff03. For the solvent model, the results strongly suggest that the ff03CMAP force field with the TIP4PD water model is the best combination for the conformer sampling. Additional tests of CHARMM36m and FB18 force fields on two phosphorylated systems suggest that the overall performance of ff03CMAP is similar to that of FB18 and better than that of CHARMM36m. These results can help other researchers to choose suitable force field and solvent models to investigate the dynamic properties of phosphorylation proteins.

Base-specific RNA force field improving the dynamics conformation of nucleotide

RNA plays a key role in numerous biological processes. Traditional experimental methods have difficulties capturing the structure and dynamic conformation of RNA. Thus, Molecular dynamic simulations (MDs) has become an essential complementary for RNA experiment. However, state-of-the-art RNA force fields have two major limitations of overestimation base stacking propensity and generation of a high ratio of intercalated conformations. Therefore, a two-step strategy was used to optimize the parameters of ff99bsc0χOL3 (named BSFF1) to improve these limitations, which as well adjusted the unbonded parameters of nucleobase heavy atoms and added ζ/α grid-based energy correction map energy term with reweighting. MD simulations of tetranucleotides indicate that BSFF1 can significantly decrease the ratio of intercalated conformations. Tests of single-strand RNA and kink-turn show that BSFF1 force field can reproduce more accurate conformers than ff99bsc0χOL3 force field. BSFF1 can also stabilize the conformers of duplex and riboswitch. The successful ab initio folding of tetraloop further supports the performance of BSFF1. These findings confirm that the newly developed force field BSFF1 can improve the conformer sampling of RNA.

Integration of Experimental Data and Use of Automated Fitting Methods in Developing Protein Force Fields

The development of accurate protein force fields has been the cornerstone of molecular simulations for the past 50 years. During this period, many lessons have been learned regarding the use of experimental target data and parameter fitting procedures. Here, we review recent advances in protein force field development. We discuss the recent emergence of polarizable force fields and the role of electronic polarization and areas in which additive force fields fall short. The use of automated fitting methods and the inclusion of additional experimental solution data during parametrization is discussed as a means to highlight possible routes to improve the accuracy of force fields even further.

Early aggregation mechanism of Aβ16-22 revealed by Markov state models

Aβ_16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aβ_16-22 peptides. The identified dimeric form adopted either globular random-coil or extended β-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel β-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of β secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel β sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aβ_16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel β sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered β-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aβ aggregation and will insight into early oligomerization for large β amyloids.

RNA-Specific Force Field Optimization with CMAP and Reweighting

RNA plays a key role in a variety of cell activities. However, it is difficult to capture its structure dynamics by the traditional experimental methods because of the inherent limitations. Molecular dynamics simulation has become a valuable complement to the experimental methods. Previous studies have indicated that the current force fields cannot accurately reproduce the conformations and structural dynamics of RNA. Therefore, an RNA-specific force field was developed to improve the conformation sampling of RNA. The distribution of ζ/α dihedrals of tetranucleotides was optimized by a reweighting method, and the grid-based energy correction map (CMAP) term was first introduced into the Amber RNA force field of ff99bsc0χOL3, named ff99OL3_CMAP1. Extensive validations of tetranucleotides and tetraloops show that ff99OL3_CMAP1 can significantly decrease the population of an incorrect structure, increase the consistency between the simulation results and experimental values for tetranucleotides, and improve the stability of tetraloops. ff99OL3_CMAP1 can also precisely reproduce the conformation of a duplex and riboswitches. These findings confirm that the newly developed force field ff99OL3_CMAP1 can improve the conformer sampling of RNA.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

Balanced Solvent Model for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) have no fixed three-dimensional (3D) structures under physiological conditions, with the content being about 51% in human proteomics. IDPs are associated with many human diseases, such as cancer, diabetes, and neurodegenerative diseases. Because IDPs do not crystallize and have diverse conformers, traditional experimental methods such as crystallization and NMR can hardly capture their conformation ensemble and just provide average structural characters of IDPs. Therefore, molecular dynamics (MD) simulations become a valuable complement to the experimental data. However, the accuracy of molecular dynamics simulation for IDPs depends on the combination of force fields and solvent models. Recently, we released an environment-specific force field (ESFF1) for IDPs, which can well reproduce the local structural properties (such as J-coupling and secondary chemical shifts). However, there is still a large deviation for the radius of gyration (R_g). Therefore, a solvent model combined with ESFF1 is necessary to capture the local and global characters for IDPs and ordered proteins. Here, we investigated the underestimation or overestimation of the solvent interaction for four solvent models (TIP3P, TIP4P-Ew, TIP4P-D, OPC) under ESFF1 and found the important ε parameter of the solvent model to play a key role in scaling R_g. A near-linear relationship between the simulation R_g and the ε parameter was used to develop the new solvent model, named TIP4P-B. The results indicate that the simulated R_g with TIP4P-B is in better agreement with the experimental observations than the other four solvent models. Simultaneously, TIP4P-B can also maintain the advantages of the ESFF1 force field for the local structural properties. Additionally, TIP4P-B can successfully sample the conformation of ordered proteins. These findings confirm that TIP4P-B is a balanced solvent model and can improve sampling R_g performance for folded proteins and IDPs.

Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions

Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.

The structural heterogeneity of α-synuclein is governed by several distinct subpopulations with interconversion times slower than milliseconds

α-Synuclein plays an important role in synaptic functions by interacting with synaptic vesicle membrane, while its oligomers and fibrils are associated with several neurodegenerative diseases. The specific monomer structures that promote its membrane binding and self-association remain elusive due to its transient nature as an intrinsically disordered protein. Here, we use inter-dye distance distributions from bulk time-resolved Förster resonance energy transfer as restraints in discrete molecular dynamics simulations to map the conformational space of the α-synuclein monomer. We further confirm the generated conformational ensemble in orthogonal experiments utilizing far-UV circular dichroism and cross-linking mass spectrometry. Single-molecule protein-induced fluorescence enhancement measurements show that within this conformational ensemble, some of the conformations of α-synuclein are surprisingly stable, exhibiting conformational transitions slower than milliseconds. Our comprehensive analysis of the conformational ensemble reveals essential structural properties and potential conformations that promote its various functions in membrane interaction or oligomer and fibril formation.

Catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase

Methane is among the most potent of the greenhouse gases, which plays a key role in global climate change. As an excellent carbon and energy source, methane can be utilized by anaerobic methane oxidizing archaea and aerobic methane oxidizing bacteria. The previous work shows that an anaerobic thermophilic enrichment culture composed of dense consortia of archaea and bacteria apparently uses partly similar pathways to oxidize the C4 hydrocarbon butane. However, the catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase is still unknown. Therefore, molecular dynamics (MD) simulation was used to investigate the dynamics differences of catalytic mechanism between methane coenzyme M reductase (MCR) and alkyl-coenzyme M reductase (ACR). At first, the binding pocket of ACR is larger than that of MCR. Then, the complex of butane and ACR is more stable than that of methane and ACR. Protein conformation cloud suggests that the position of methane is dynamics and methane escapes from the binding pocket of ACR during most of the simulation time, while butane tightly binds in the pocket of ACR. The hydrophobic interactions between butane and ACR are more and stronger than those between methane and ACR. At the same time, the binding free energy between butane and ACR is significantly lower than that between methane and ACR. The dynamics correlation network indicates that the transformation of information flow for ACR-butane is smoother than that for ACR-methane. The shortest pathway for ACR-butane is from Gln144, Ala141, Hie135, Ile133, Ala160, Arg206, Asp97, Met94, Tyr347 to Phe345 with synergistic effect for two butane molecules. This study can insight into the catalytic mechanism for butane/ACR complex.

Extensive evaluation of environment-specific force field for ordered and disordered proteins

Intrinsically disordered proteins (IDPs) have no fixed tertiary structure under physiological conditions and are associated with many human diseases. Because IDPs have the characteristic of possessing diverse conformations, current experimental methods cannot capture all the conformations of IDPs. However, molecular dynamics simulation can sample these atomistically diverse conformations as a valuable complement to experimental data. To accurately describe the properties of IDPs, the environment-specific precise force field (ESFF1) was successfully released to reproduce the conformer character of ordered and disordered proteins. Here, three typical IDPs and thirteen folded proteins were used to further evaluate the performance of this force field. The results indicate that the NMR observables of ESFF1 better approach experimental data than do those of ff14SB for IDPs. The sampling conformations by ESFF1 are more diverse than those of ff14SB. For folded proteins, these force fields have comparable performances for reproducing conformers. Therefore, ESFF1 can be used to reveal the model of sequence-disorder-function for IDPs.

Dissecting the Conformational Free Energy of a Small Peptide in Solution

The free energy surface of a small peptide was analyzed based on an unbiased microsecond molecular dynamics simulation. The peptide sampled disordered conformational ensembles of distinct compactness, and its free energy was decomposed into separate contributions from the intramolecular potential energy, conformational entropy, and solvation free energy. The latter was further broken down into enthalpic and entropic contributions due to peptide-water and water-water interactions. This decomposition was enabled by a generalized linear response relation between the peptide-water interaction energy and the solvation free energy, which was empirically parametrized by explicit solvation free energy calculations for representative peptide conformations. This full dissection of the peptide free energy identifies individual contributions that stabilize and destabilize compact and extended peptide conformational ensembles and reveals the origin of a free energy barrier associated with transitions between them.

Assessment of transferable forcefields for protein simulations attests improved description of disordered states and secondary structure propensities, and hints at multi-protein systems as the next challenge for optimization

Continuous assessment of transferable forcefields for molecular simulations is essential to identify their weaknesses and direct improvement efforts. The latest efforts focused on better describing disordered proteins while retaining proper description of folded domains, important because forcefields of the previous generations produce overly compact disordered states. Such improvements should additionally alleviate the related problem of over-stabilized protein–protein interactions, which has been largely overlooked. Here we evaluated three state-of-the-art forcefields, current flagships of their respective developers, optimized for ordered and disordered proteins: CHARMM36m with its recommended corrected TIP3P* water, ff19SB with the recommended OPC water, and the 2019 a99SBdisp forcefield by D. E. Shaw Research with its modified TIP4P water; plus ff14SB with TIP3P as an example of the former generation of forcefields. Our evaluation entailed simulations of (i) multiple copies of a protein that is highly soluble yet undergoes weak dimerization, (ii) a disordered peptide with low, well-characterized alpha helical propensity, and (iii) a peptide known to form insoluble β-aggregates. Our results recapitulate ff14SB-TIP3P over-stabilizing aggregates and secondary structures and place a99SBdisp-TIP4PD at the other end i.e. predicting overly weak intermolecular interactions despite reasonably predicting secondary structure propensities. In-between, CHARMM36m-TIP3P* still over-stabilizes aggregates but predicts residue-wise alpha helical propensities in solution slightly better than ff19SB-OPC, while ff19SB-OPC poses the best prediction of weak dimerization of the soluble protein still predicting aggregation of the β-peptides. This independent assessment shows that the claimed forcefield improvements are real, but also that a right balance between noncovalent attraction and repulsion has not yet been reached. We thus propose developers to consider systems like those tested here in their forcefield tuning protocols. Last, the good performance of CHARMM36m-TIP3P* further shows that tuning 3-point water models might still be an alternative to the more costly 4-point models like OPC and TIP4PD.

Folding and self-assembly of short intrinsically disordered peptides and protein regions

Proteins and peptide fragments are highly relevant building blocks in self-assembly for nanostructures with plenty of applications. Intrinsically disordered proteins (IDPs) and protein regions (IDRs) are defined by the absence of a well-defined secondary structure, yet IDPs/IDRs show a significant biological activity. Experimental techniques and computational modelling procedures for the characterization of IDPs/IDRs are discussed. Directed self-assembly of IDPs/IDRs allows reaching a large variety of nanostructures. Hybrid materials based on the derivatives of IDPs/IDRs show a promising performance as alternative biocides and nanodrugs. Cell mimicking, in vivo compartmentalization, and bone regeneration are demonstrated for IDPs/IDRs in biotechnological applications. The exciting possibilities of IDPs/IDRs in nanotechnology with relevant biological applications are shown.

Recent Force Field Strategies for Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are widely distributed across eukaryotic cells, playing important roles in molecular recognition, molecular assembly, post-translational modification, and other biological processes. IDPs are also associated with many diseases such as cancers, cardiovascular diseases, and neurodegenerative diseases. Due to their structural flexibility, conventional experimental methods cannot reliably capture their heterogeneous structures. Molecular dynamics simulation becomes an important complementary tool to quantify IDP structures. This review covers recent force field strategies proposed for more accurate molecular dynamics simulations of IDPs. The strategies include adjusting dihedral parameters, adding grid-based energy correction map (CMAP) parameters, refining protein-water interactions, and others. Different force fields were found to perform well on specific observables of specific IDPs but also are limited in reproducing all available experimental observables consistently for all tested IDPs. We conclude the review with perspective areas for improvements for future force fields for IDPs.

Disaggregation mechanism of prion amyloid for tweezer inhibitor

The aggregation of amyloid has been an important event in the pathology of amyloidogenicity. A number of small molecules have been designed for Amyloidosis treatment. Molecular tweezer CLR01, a potential drug for misfolded β-amyloids inhibition, was reportedly bind directly to Lysine residues and interrupt oligomerization. However, the disaggregation mechanism of amyloid for this inhibitor is unclear. Here we used long timescale of molecular dynamic simulation to reveal the mechanism of disaggregation for pentamer prion amyloid. Molecular docking and molecular dynamics simulation demonstrate that CLR01 is attached with Lysine₂₂₂ nitrogen by π-cation interaction of its nine aromatic rings and formation of salt bridge/hydrogen bond of one of the two rotatable peripheral anionic phosphate groups. Upon CLR01 binding, we found a major shifting occurs in initial conformation of the oligomer and stretch out the N-terminal chain A from the rest of the amyloid which seems to be the first stage of disaggregated the fibrils slowly yet efficiently. Moreover, the CLR01 remodelled the pentamer Prion_220-272 into a compact structure which might be the resistant conformation for further oligomerization. Our work will contribute to better understand the interaction and deterioration mechanism of molecular tweezer for prions and similar amyloids, and offer significant insights into therapeutic development for Amyloidosis treatment.

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.

Conformational Preferences of an Intrinsically Disordered Protein Domain: A Case Study for Modern Force Fields

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/E_b, ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

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.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins

Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 μs (five replicates of 200 ns, 300 ns, 1 μs, or 2 μs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, β amyloid (Aβ)-40, Aβ-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of β sheets for β-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of β-amyloid 42 starting from two different initial structures (Aβ42^1Z0Q and Aβ42^model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 μs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.

Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins

The need for accurate and efficient force fields for modeling 3D structures of macrobiomolecules and in particular intrinsically disordered proteins (IDPs) has increased with recent findings to associate IDPs and human diseases. However, most conventional protein force fields and recent IDP-specific force fields are limited in reproducing accurate structural features of IDPs. Here, we present an environmental specific precise force field (ESFF1) based on CMAP corrections of 71 different sequence environments to improve the accuracy and efficiency of MD simulation for both IDPs and folded proteins. MD simulations of 84 different short peptides, IDPs, and structured proteins show that ESFF1 can accurately reproduce spectroscopic properties for different peptides and proteins whether they are disordered or ordered. The successful ab initio folding of five fast-folding proteins further supports the reliability of ESFF1. The extensive analysis documented here shows that ESFF1 is able to achieve a reasonable balance between ordered and disordered states in protein simulations.