Balanced Amino-Acid-Specific Molecular Dynamics Force Field for the Realistic Simulation of Both Folded and Disordered Proteins

Atomistic molecular dynamics simulations of intrinsically disordered proteins

Recent years have seen remarkable gains in the accuracy of atomistic molecular dynamics (MD) simulations of intrinsically disordered proteins (IDPs) and expansion in the types of calculated properties that can be directly compared with experimental measurements. These advances occurred due to the use of IDP-tested force fields and the porting of MD simulations to GPUs and other computational technologies. All-atom MD simulations are now explaining the sequence-dependent dynamics of IDPs; elucidating the mechanisms of their binding to other proteins, nucleic acids, and membranes; revealing the modes of drug action on them; and characterizing their phase separation. Artificial intelligence (AI) and machine learning (ML) are further expanding the reach of atomistic MD simulations.

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

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

Coil-Library-Derived Amino-Acid-Specific Side-Chain χ1 Dihedral Angle Potentials for AMBER-Type Protein Force Field

The successful simulation of proteins by molecular dynamics (MD) critically depends on the accuracy of the applied force field. Here, we modify the AMBER-family ff99SBnmr2 force field through improvements to the side-chain χ1 dihedral angle potentials in a residue-specific manner using conformational dihedral angle distributions from an experimental coil library as targets. Based on significant deviations observed for the parent force field with respect to the coil library, the χ1 dihedral angle potentials of seven amino acids were modified, namely, Val, Ser, His, Asn, Trp, Tyr, and Phe. The new force field, named ff99SBnmr2Chi1, was benchmarked against NMR-derived χ1 rotamer populations of denatured proteins, overall resulting in much better agreement and without any noticeable adverse consequences on the quality of the simulation of folded proteins. The new force field should allow more realistic modeling of protein side-chain properties by MD of both folded and unfolded protein systems, such as for the better in-silico characterization of protein-protein and protein-ligand interactions.

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

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

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.

Uncovering the Association Mechanism between Two Intrinsically Flexible Proteins

The understanding of protein-protein interaction mechanisms is key to the atomistic description of cell signaling pathways and for the development of new drugs. In this context, the mechanism of intrinsically disordered proteins folding upon binding has attracted attention. The VirB9 C-terminal domain (VirB9Ct) and the VirB7 N-terminal motif (VirB7Nt) associate with VirB10 to form the outer membrane core complex of the Type IV Secretion System injectisome. Despite forming a stable and rigid complex, VirB7Nt behaves as a random coil, while VirB9Ct is intrinsically dynamic in the free state. Here we combined NMR, stopped-flow fluorescence, and computer simulations using structure-based models to characterize the VirB9Ct-VirB7Nt coupled folding and binding mechanism. Qualitative data analysis suggested that VirB9Ct preferentially binds to VirB7Nt by way of a conformational selection mechanism at lower temperatures. However, at higher temperatures, energy barriers between different VirB9Ct conformations are more easily surpassed. Under these conditions the formation of non-native initial encounter complexes may provide alternative pathways toward the native complex conformation. These observations highlight the intimate relationship between folding and binding, calling attention to the fact that the two molecular partners must search for the most favored intramolecular and intermolecular interactions on a rugged and funnelled conformational energy landscape, along which multiple intermediates may lead to the final native state.

High-fidelity, hyper-accurate, and evolved mutants rewire atomic-level communication in CRISPR-Cas9

The high-fidelity (HF1), hyper-accurate (Hypa), and evolved (Evo) variants of the CRISPR-associated protein 9 (Cas9) endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in recognition subdomain 3 (Rec3) mediate specificity in these variants are poorly understood. Here, solution nuclear magnetic resonance and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic His-Asn-His (HNH) domain. The variants remodel the communication sourcing from the Rec3 α helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Cyclic Peptide Linker Design and Optimization by Molecular Dynamics Simulations

Cyclic peptides are an emerging therapeutic modality that can target protein-protein interaction sites with high affinity and selectivity. A common medicinal chemistry strategy for the optimization of peptide hits is conformational stabilization through macrocyclization. We present a method based on explicit solvent enhanced sampling molecular dynamics simulations for estimating the impact of varying linker lengths and chemistry on the conformational stability of a peptide. The method is demonstrated on three cyclic peptide series that bind to proteins PCSK9, trypsin, and MDM2 adopting loop, β-sheet, and helical secondary structures. In general, the simulations show greater solution stability of the receptor-bound conformation for the higher-affinity peptides, consistent with the idea that preorganizing a ligand for binding can enhance binding affinity. The impact of the force field and sampling is discussed for one series that does not follow this trend. We have successfully applied this method to internal discovery programs to design peptides with increased potency and chemical stability.

High-Fidelity, Hyper-Accurate, and Evolved Mutants Rewire Atomic Level Communication in CRISPR-Cas9

The Cas9-HF1, HypaCas9, and evoCas9 variants of the Cas9 endonuclease are critical tools to mitigate off-target effects in the application of CRISPR-Cas9 technology. The mechanisms by which mutations in the Rec3 domain mediate specificity in these variants are poorly understood. Here, solution NMR and molecular dynamics simulations establish the structural and dynamic effects of high-specificity mutations in Rec3, and how they propagate the allosteric signal of Cas9. We reveal conserved structural changes and peculiar dynamic differences at regions of Rec3 that interface with the RNA:DNA hybrid, transducing chemical signals from Rec3 to the catalytic HNH domain. The variants remodel the communication sourcing from the Rec3 α-helix 37, previously shown to sense target DNA complementarity, either directly or allosterically. This mechanism increases communication between the DNA mismatch recognition helix and the HNH active site, shedding light on the structure and dynamics underlying Cas9 specificity and providing insight for future engineering principles.

Application of molecular dynamics simulation in self-assembled cancer nanomedicine

Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.

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.

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.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

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.

Structural Basis for Reduced Dynamics of Three Engineered HNH Endonuclease Lys-to-Ala Mutants for the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Associated 9 (CRISPR/Cas9) Enzyme

Many bacteria possess type-II immunity against invading phages or plasmids known as the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) system to detect and degrade the foreign DNA sequences. The Cas9 protein has two endonucleases responsible for double-strand breaks (the HNH domain for cleaving the target strand of DNA duplexes and RuvC domain for the nontarget strand, respectively) and a single-guide RNA-binding domain where the RNA and target DNA strands are base-paired. Three engineered single Lys-to-Ala HNH mutants (K810A, K848A, and K855A) exhibit an enhanced substrate specificity for cleavage of the target DNA strand. We report in this study that in the wild-type (wt) enzyme, D835, Y836, and D837 within the Y836-containing loop (comprising E827-D837) adjacent to the catalytic site have uncharacterizable broadened 1H15N nuclear magnetic resonance (NMR) features, whereas remaining residues in the loop have different extents of broadened NMR spectra. We find that this loop in the wt enzyme exhibits three distinct conformations over the duration of the molecular dynamics simulations, whereas the three Lys-to-Ala mutants retain only one conformation. The versatility of multiple alternate conformations of this loop in the wt enzyme could help to recruit noncognate DNA substrates into the HNH active site for cleavage, thereby reducing its substrate specificity relative to the three mutants. Our study provides further experimental and computational evidence that Lys-to-Ala substitutions reduce dynamics of proteins and thus increase their stability.

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.

Structural and dynamic insights into the HNH nuclease of divergent Cas9 species

CRISPR-Cas9 is a widely used biochemical tool with applications in molecular biology and precision medicine. The RNA-guided Cas9 protein uses its HNH endonuclease domain to cleave the DNA strand complementary to its endogenous guide RNA. In this study, novel constructs of HNH from two divergent organisms, G. stearothermophilus (GeoHNH) and S. pyogenes (SpHNH) were engineered from their respective full-length Cas9 proteins. Despite low sequence similarity, the X-ray crystal structures of these constructs reveal that the core of HNH surrounding the active site is conserved. Structure prediction of the full-length GeoCas9 protein using Phyre2 and AlphaFold2 also showed that the crystallographic construct of GeoHNH represents the structure of the domain within the full-length GeoCas9 protein. However, significant differences are observed in the solution dynamics of structurally conserved regions of GeoHNH and SpHNH, the latter of which was shown to use such molecular motions to propagate the DNA cleavage signal. Indeed, molecular simulations show that the intradomain signaling pathways, which drive SpHNH function, are non-specific and poorly formed in GeoHNH. Taken together, these outcomes suggest mechanistic differences between mesophilic and thermophilic Cas9 species.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

Molecular Basis of Small-Molecule Binding to α-Synuclein

Intrinsically disordered proteins (IDPs) are implicated in many human diseases. They have generally not been amenable to conventional structure-based drug design, however, because their intrinsic conformational variability has precluded an atomic-level understanding of their binding to small molecules. Here we present long-time-scale, atomic-level molecular dynamics (MD) simulations of monomeric α-synuclein (an IDP whose aggregation is associated with Parkinson’s disease) binding the small-molecule drug fasudil in which the observed protein–ligand interactions were found to be in good agreement with previously reported NMR chemical shift data. In our simulations, fasudil, when bound, favored certain charge–charge and π-stacking interactions near the C terminus of α-synuclein but tended not to form these interactions simultaneously, rather breaking one of these interactions and forming another nearby (a mechanism we term dynamic shuttling). Further simulations with small molecules chosen to modify these interactions yielded binding affinities and key structural features of binding consistent with subsequent NMR experiments, suggesting the potential for MD-based strategies to facilitate the rational design of small molecules that bind with disordered proteins.

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

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.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

When Order Meets Disorder: Modeling and Function of the Protein Interface in Fuzzy Complexes

The degree of proteins structural organization ranges from highly structured, compact folding to intrinsic disorder, where each degree of self-organization corresponds to specific functions: well-organized structural motifs in enzymes offer a proper environment for precisely positioned functional groups to participate in catalytic reactions; at the other end of the self-organization spectrum, intrinsically disordered proteins act as binding hubs via the formation of multiple, transient and often non-specific interactions. This review focusses on cases where structurally organized proteins or domains associate with highly disordered protein chains, leading to the formation of interfaces with varying degrees of fuzziness. We present a review of the computational methods developed to provide us with information on such fuzzy interfaces, and how they integrate experimental information. The discussion focusses on two specific cases, microtubules and homologous recombination nucleoprotein filaments, where a network of intrinsically disordered tails exerts regulatory function in recruiting partner macromolecules, proteins or DNA and tuning the atomic level association. Notably, we show how computational approaches such as molecular dynamics simulations can bring new knowledge to help bridging the gap between experimental analysis, that mostly concerns ensemble properties, and the behavior of individual disordered protein chains that contribute to regulation functions.

Membrane-bound KRAS approximates an entropic ensemble of configurations

KRAS4B is a membrane-anchored signaling protein and a primary target in cancer research. Predictions from molecular dynamics simulations that have previously shaped our mechanistic understanding of KRAS signaling disagree with recent experimental results from neutron reflectometry, NMR, and thermodynamic binding studies. To gain insight into these discrepancies, we compare this body of biophysical data to back-calculated experimental results from a series of molecular simulations that implement different subsets of molecular interactions. Our results show that KRAS4B approximates an entropic ensemble of configurations at model membranes containing 30% phosphatidylserine lipids, which is not significantly shaped by interactions between the globular G-domain of KRAS4B and the lipid membrane. These findings revise our understanding of KRAS signaling and promote a model in which the protein samples the accessible conformational space in a near-uniform manner while being available to bind to effector proteins.

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.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-μs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-μs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to μs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S² order parameter. The method is demonstrated for ²H-relaxation of CH₂D methyl moieties and cross-correlated relaxation of CH₃ groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 μs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

Accurate Structure Prediction for Protein Loops Based on Molecular Dynamics Simulations with RSFF2C

Protein loops, connecting the α-helices and β-strands, are involved in many important biological processes. However, due to their conformational flexibility, it is still challenging to accurately determine three-dimensional (3D) structures of long loops experimentally and computationally. Herein, we present a systematic study of the protein loop structure prediction via a total of ∼850 μs molecular dynamics (MD) simulations. For a set of 15 long (10-16 residues) and solvent-exposed loops, we first evaluated the performance of four state-of-the-art loop modeling algorithms, DaReUS-Loop, Sphinx, Rosetta-NGK, and MODELLER, on each loop, and none of them could accurately predict the structures for most loops. Then, temperature replica exchange molecular dynamics (REMD) simulations were conducted with three recent force fields, RSFF2C with TIP3P water model, CHARMM36m with CHARMM-modified TIP3P, and AMBER ff19SB with OPC. We found that our recently developed residue-specific force field RSFF2C performed the best and successfully predicted 12 out of 15 loops with a root-mean-square deviation (RMSD) < 1.5 Å. As an alternative with lower computational cost, normal MD simulations at high temperatures (380, 500, and 620 K) were investigated. Temperature-dependent performance was observed for each force field, and, for RSFF2C+TIP3P, we found that three independent 100-ns MD simulations at 500 K gave comparable results with REMD simulations. These results suggest that MD simulations, especially with enhanced sampling techniques such as replica exchange, with the RSFF2C force field could be useful for accurate loop structure prediction.

Influence of a Ser111-phosphorylation on Rab1b GTPase conformational dynamics studied by advanced sampling simulations

Rab GTPases constitute the largest branch of the Ras protein superfamily that regulate intra-cellular membrane trafficking. Their signaling activity is mediated by the transition between an active GTP-bound state and an inactive GDP-bound state. In the inactive state the switch I and II segments adopt largely disordered flexible conformations, whereas in the active state these regions are in well-defined conformations. The switch I and II states are central for recognition of Rab GTPases by interacting partners. Phosphorylation of the Rab1b-GTPase at residue Ser111 (pS111) results in modulation of the signaling activity due to alterations of the protein interaction interface and also due to modulation of the conformational flexibility. We have studied the flexibility of native and pS111-Rab1b in complex with GTP or GDP using extensive Molecular Dynamics (MD) simulations and an advanced sampling method called DIhedral Angle-biasing potential Replica-Exchange Molecular dynamics (DIA-REMD). The DIA-REMD method promotes backbone and side chain dihedral transitions along a series of replica simulations in selected protein segments and through exchanges also improves sampling in an unbiased reference simulation. Application to the Rab1b system results in significantly enhanced sampling of different switch I/II conformational states in the GDP-bound Rab1b state. The pS111 modification is found to reduce the conformational flexibility even in the presence of GDP, which may influence signaling activities. The stabilizing effect can be attributed to the formation of additional surface salt bridges between Arg-residues and pS111 not present in the native structure. The DIA-REMD method could be a valuable approach for studying also other signaling proteins that contain flexible segments.

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.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

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.

Systematic Differences between Current Molecular Dynamics Force Fields To Represent Local Properties of Intrinsically Disordered Proteins

The prevalence of intrinsically disordered proteins (IDPs) and protein regions in structural biology has prompted the recent development of molecular dynamics (MD) force fields for the more realistic representations of such systems. Using experimental nuclear magnetic resonance backbone scalar ³J-coupling constants of the intrinsically disordered proteins α-synuclein and amyloid-β in their native aqueous environment as a metric, we compare the performance of four recent MD force fields, namely, AMBER ff14SB, CHARMM C36m, AMBER ff99SB-disp, and AMBER ff99SBnmr2, by partitioning the polypeptides into an overlapping series of heptapeptides for which a cumulative total of 276 μs MD simulations were performed. The results show substantial differences between the different force fields at the individual residue level. Except for ff99SBnmr2, the force fields systematically underestimate the scalar ³J(H^N,Hα)-couplings due to an underrepresentation of β-conformations and an overrepresentation of either α- or PP_II conformations. The study demonstrates that the incorporation of coil library information in modern MD force fields, as shown here for ff99SBnmr2, provides substantially improved performance and more realistic sampling of the local backbone dihedral angles of IDPs as reflected by the good accuracy of the computed scalar ³J(H^N,Hα)-couplings with less than 0.5 Hz error. Such force fields will enable a better understanding of how structural dynamics and thermodynamics influence the IDP function. Although the methodology based on heptapeptides used here does not allow the assessment of potential intramolecular long-range interactions, its computational affordability permits well-converged simulations that can be easily parallelized. This should make the quantitative validation of intrinsic disorder observed in MD simulations of polypeptides with experimental scalar J-couplings widely applicable.

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.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.

Targeting Intrinsically Disordered Proteins through Dynamic Interactions

Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule-IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.