Refined Force Field Research Articles

Absolute Binding Free Energies with OneOPES.

The calculation of absolute binding free energies (ABFEs) for protein-ligand systems has long been a challenge. Recently, refined force fields and algorithms have improved the quality of the ABFE calculations. However, achieving the level of accuracy required to inform drug discovery efforts remains difficult. Here, we present a transferable enhanced sampling strategy to accurately calculate absolute binding free energies using OneOPES with simple geometric collective variables. We tested the strategy on two protein targets, BRD4 and Hsp90, complexed with a total of 17 chemically diverse ligands, including both molecular fragments and drug-like molecules. Our results show that OneOPES accurately predicts protein-ligand binding affinities with a mean unsigned error within 1 kcal mol-1 of experimentally determined free energies, without the need to tailor the collective variables to each system. Furthermore, our strategy effectively samples different ligand binding modes and consistently matches the experimentally determined structures regardless of the initial protein-ligand configuration. Our results suggest that the proposed OneOPES strategy can be used to inform lead optimization campaigns in drug discovery and to study protein-ligand binding and unbinding mechanisms.

Molecular dynamics simulations of nucleosomes are coming of age

AbstractUnderstanding the function of eukaryotic genomes, including the human genome, is undoubtedly one of the major scientific challenges of the 21st century. The cornerstone of eukaryotic genome organization is nucleosomes—elementary building blocks of chromatin about 10 nm in size that wrap DNA around an octamer of histone proteins. Nucleosomes are integral players in all genomic processes, including transcription, DNA replication and repair. They mediate genome regulation at the epigenetic level, bridging the discrete nature of the genetic information encoded in DNA with the analog physical nature of the intermolecular interactions required to access that information. Due to their relatively large size and dynamic nature, nucleosomes are difficult objects for experimental characterization. Molecular dynamics (MD) simulations have emerged over the years as a useful tool to complement experimental studies. Particularly in recent years, advances in computing power, refinement of MD force fields and codes have opened up new frontiers in terms of simulation timescales and quality for nucleosomes and related systems. It has become possible to elucidate in atomistic detail their functional dynamics modes such as DNA unwrapping and sliding, to characterize the effects of epigenetic modifications, DNA and protein sequence variation on nucleosome structure and stability, to describe the mechanisms governing nucleosome interactions with chromatin‐associated proteins and the formation of supranucleosome structures. In this review, we systematically analyzed all‐atom MD simulation studies of nucleosomes and related structures published since 2018 and discussed their relevance in the context of older studies, experimental data, and related coarse‐grained and multiscale studies.This article is categorized under: Software > Molecular Modeling Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Structure and Mechanism > Computational Biochemistry and Biophysics

Evaluating Polarizable Biomembrane Simulations against Experiments.

Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.

Taming Conformational Heterogeneity on Ion Racetrack to Unveil Principles that Drive Membrane Permeation of Cyclosporines.

Our study reveals the underlying principles governing the passive membrane permeability in three large N-methylated macrocyclic peptides (N-MeMPs): cyclosporine A (CycA), Alisporivir (ALI), and cyclosporine H (CycH). We determine a series of conformers required for robust passive membrane diffusion and those relevant to other functions, such as binding to protein targets or intermediates, in the presence of solvent additives. We investigate the conformational interconversions and establish correlations with the membrane permeability. Nuclear magnetic resonance (NMR) and cyclic ion-mobility spectrometry-mass spectrometry (cIMS-MS) are employed to characterize conformational heterogeneity and identify cis-amides relevant for good membrane permeability. In addition, ion mobility selected cIMS-MS and infrared (IR) multiple-photon dissociation (IRMPD) spectroscopy experiments are conducted to evaluate the energy barriers between conformations. We observe that CycA and ALI, both cyclosporines with favorable membrane permeabilities, display multiple stable and well-defined conformers. In contrast, CycH, an epimer of CycA with limited permeability, exhibits fewer and fewer stable conformers. We demonstrate the essential role of the conformational shift from the aqueous cis MeVal11-MeBmt1 state (A1) to the closed conformation featuring cis MeLeu9-MeLeu10 (C1) in facilitating membrane permeation. Additionally, we highlight that the transition from A1 to the all-trans open conformation (O1) is specifically triggered by the presence of CaCl2. We also capture a set of conformers with cis Sar3-MeLeu4, MeLeu9-MeLeu10, denoted as I. Conformationally selected cIMS-MS and IRMPD data of [CycA+Ca]2+ show immediate repopulation of the original population distribution, suggesting that CaCl2 smooths out the energy barriers. Finally, our work presents an improved sampling molecular dynamics approach based on a refined force field that not only consistently and accurately captures established conformers of cyclosporines but also exhibits strong predictive capabilities for novel conformers.

Observation and Thermodynamic Analysis of Reversible Fluid Isobars in Isolated Single Digit Nanopore Carbon Nanotubes

Recent interest in one-dimensionally confined fluids, where confinement approaches molecular dimensions, has demonstrated exceptionally high fluxes from slip flow and large distortions of fluid phase boundaries [1,2]. The Center for Enhanced Nanofluidic Transport (CENT) was recently formed as an intellectual hub for studying extreme fluid confinement in what we label Single Digit Nanopores (SDNs). To this end, in this work, we note that predicting such phenomena for a given conduit dimension has been confounded by a dearth of fundamental thermodynamic measurements and analysis as a function of confinement diameter and wall composition. Herein, we utilize Raman spectroscopy and carbon nanotubes with diameters less than 3 nm to identify and study newly observed, thermally driven, fluid adsorbed phases under conditions at constant ambient pressure, tracing what we label a fluid isobar as a function of temperature. We show that the Raman radial breathing(-like) vibrational modes (RBLMs) can serve as fluid sensors to trace the dependence of the fluid isobar and allow study of basic thermodynamics of fluid phase transition as a function of conduit curvature and confinement volume. We study fluids on the exterior and interior of double- and single- walled carbon nanotubes spanning 0.9 to 3.3 nm in diameter, showing that the fluid location can be disambiguated on double walled tubes by using an adaptation of elastic shell theory that allows deconvolution of inner and outer shell frequencies with fluid location. We find that the RBLM FWHM systematically changes as a unique signature of the location of the fluid phase near the CNT wall. A Buckingham potential is used along with electron diffraction assigned DWNT to provide a high accuracy estimate of the carbon-carbon coupling constant for two elastic shells, an important structural parameter for DWNT. The thermodynamic analysis of isobars measured for more than 25 CNTs under various environmental conditions is conducted by comparison to a new fluid equation of state that provides a continuum description of the fluid isobar under confinement. Our work sets forth a new and reliable approach for studying phase change of different fluids in large numbers of isolated single digit nanopores that may inform the refinement of fluid force fields and motivate theoretical refinements to the calculation of thermodynamic properties under confinement.[1] S. Faucher, N. Aluru, M. Z. Bazant, D. Blankschtein, A. H. Brozena, J. Cumings, J. Pedro De Souza, M. Elimelech, R. Epsztein, J. T. Fourkas, A. G. Rajan, H. J. Kulik, A. Levy, A. Majumdar, C. Martin, M. McEldrew, R. P. Misra, A. Noy, T. A. Pham, M. Reed, E. Schwegler, Z. Siwy, Y. Wang, M. S. Strano. Critical Knowledge Gaps in Mass Transport through Single-Digit Nanopores: A Review and Perspective. J. Phys. Chem. C 123, 21309–21326 (2019).[2] K. V. Agrawal, S. Shimizu, L. W. Drahushuk, D. Kilcoyne, M. S. Strano. Observation of extreme phase transition temperatures of water confined inside isolated carbon nanotubes. Nat. Nanotechnol. 12, 267–273 (2017)

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions.

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Coacervation and adhesion mechanism of PEI/TA hydrogel revealed by molecular dynamics simulations.

Coacervate-derived polyethyleneimine (PEI)/thioctic acid (TA) hydrogels can efficiently repel interfacial water and achieve robust underwater adhesion on diverse substrates, thus representing a promising material for a wide range of biomedical uses such as wound dressings, hemostatic agents, drug carriers, and tissue repair. Understanding the cohesion and adhesion dynamics in this coacervate-to-hydrogel strategy is critical for the development of underwater bioadhesive. However, these mechanisms remain obscure due to the difficulty of obtaining the atomic details. In this work, we first applied quantum calculations and classical mechanics models to obtain accurate force field parameters for PEI and TA. Based on the refined force field, we employed microsecond-long all-atom molecular dynamics (MD) simulations to investigate the coacervation of PEI/TA complex and its adhesion on glass and polypropylene. Our simulations revealed that the phase transition in PEI/TA coacervation was jointly achieved by two types of aggregation, i.e., TA-TA aggregation through hydrophobic associations and PEI-TA aggregation driven by electrostatic interactions. In addition, MD simulations further depicted the distinct adhesion modes of PEI/TA on the two substrates. TA assembled on top of the hydrophobic polypropylene substrate, replacing the interfacial water, whereas PEI amines penetrated the interfacial water to directly bind to the silanols of the hydrophilic glass substrate. The tight adhesion of PEI/TA hydrogel was then quantified by the diffusion coefficient of molecular moieties, which reflected the extremely low mobility of the bound moieties. These findings deepen our understanding of the formation and adhesion of PEI/TA hydrogel at the atomic level and provide guidance for the future design of advanced bioadhesive for potential clinical applications.

Force field refinement for reproducing experimental infrared spectra of ionic liquids.

We employ polarizable molecular dynamics simulations with the newly developed FFGenOpt parametrization tool to reproduce IR spectra of several ionic liquid cations and anions in the gas phase. Our results show that polarizable force fields in the bulk phase provide a reasonable compromise between computational effort and accuracy for investigating IR spectra when treating the transition from gas to liquid phase carefully. Although collectivity seems to play only a minor role, the liquid phase not only changes the electrostatic environment of the molecules but also introduces friction and intermolecular interactions altering the IR spectrum significantly. In addition to the classical force field approach, we also tested if the additional computational effort of machine learning potentials justifies their application in reproducing IR spectra. However, the main purpose of this work is to improve the quality of polarizable force fields concerning vibrations and not the prediction of IR spectra which can be better done with quantum-mechanical cluster approaches.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes.

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

Refined Force Field for Liquid Sulfolane with Particular Emphasis to Its Transport Characteristics.

An all-atom force field dedicated to capturing the properties of multifunctional sulfolane is necessary. In addition to being an excellent solvent and extractor, sulfolane is also a frequently investigated component for battery electrolytes in recent times. Given this, theoretically capturing its transport properties is essential. However, given the rather high shear viscosity of liquid sulfolane and its polar aprotic nature, formulating an appropriate non-polarizable force field for this compound remains a challenge. Starting from a generic force field, we report a refined force field for sulfolane which quantitatively captures its bulk properties, resulting in significantly improved estimates for self-diffusion constant and shear viscosity of sulfolane in comparison to force fields reported hitherto. Density, self-diffusion constant, and shear viscosity were determined between temperatures (303 and 398) K and at 1 bar pressure. All properties determined from the refined force field are in good agreement with experiments. The refined model employs atomic site charges obtained from the density-derived electrostatic and chemical (DDEC6) method for liquid sulfolane modeled using quantum density functional theory. Lennard-Jones parameters were refined using quantum potential energy scans. Despite possessing a large dipole moment, the large molecular size of sulfolane partially disrupts intermolecular dipolar ordering in liquid sulfolane. Molecular dipoles of near neighbor sulfolane, however, retain a partial preference for antiparallel orientation even at the highest temperatures investigated here.

Force Field Refinement Based on High Resolution Protein Structures

Force Field Refinement Based on High Resolution Protein Structures

Refined Empirical Force Field to Model Protein-Self-Assembled Monolayer Interactions Based on AMBER14 and GAFF.

Understanding protein interaction with material surfaces is important for the development of nanotechnological devices. The structures and dynamics of proteins can be studied via molecular dynamics (MD) if the protein-surface interactions can be accurately modeled. To answer this question, we computed the adsorption free energies of peptides (representing eleven different amino acids) on a hydrophobic self-assembled monolayer (CH3-SAM) and compared them to the benchmark experimental data set. Our result revealed that existing biomolecular force fields, GAFF and AMBER ff14sb, cannot reproduce the experimental peptide adsorption free energies by Wei and Latour (Langmuir, 2009, 25, 5637-5646). To obtain the improved force fields, we systematically tuned the Lennard-Jones parameters of selected amino acid sidechains and the functional group of SAM with repeated metadynamics and umbrella sampling simulations. The final parameter set has yielded a significant improvement in the free energy values with R = 0.83 and MSE = 0.65 kcal/mol. We applied the refined force field to predict the initial adsorption orientation of lysozyme on CH3-SAM. Two major orientations-face-down and face-up-were predicted. Our analysis on the protein structure, solvent accessible surface area, and binding of native ligand NAG3 suggested that lysozyme in the face-up orientation can remain active after initial adsorption. However, because of its weaker affinity (ΔΔG = 7.86 kcal/mol) for the ligand, the bioactivity of the protein is expected to reduce. Our work facilitates the use of MD for the study of protein-SAM systems. The refined force field compatible with GROMACS is available at https://cbbio.cis.um.edu.mo/software/SAMFF .

An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

Summary X-ray crystallography is the primary tool for biomolecular structural determination. However, contacts formed through the crystal lattice are known to affect structures, especially for small and flexible molecules such as DNA oligomers, by introducing significant structural changes in comparison to solution. Furthermore, why molecules crystallize in certain symmetry groups, which role crystallization additives play, and whether they are just innocuous and unspecific crystallization catalysts remain unclear. By using one of the currently best-performing DNA force fields and applying significant computational effort, we described the nature of intermolecular forces that stabilize B-DNA crystals in various symmetry groups and solvent environments with an unprecedented level of detail. We showed a tight coupling between the lattice stability and the type of crystallization additives and that certain symmetry groups are stable only in the presence of a specific additive. Additives and crystal contacts induce small but non-negligible changes in the physical properties of DNA.

Systematic evaluation and refinement of existing all-atom force fields for the simulation of liquid acetonitrile.

The reliability of a molecular dynamics (MD) simulation study mainly depends on the accuracy of the applied force field. Unlike the ability of some potential models for reasonably predicting the thermodynamic properties of acetonitrile (ACN), simulated dynamical properties such as self-diffusion are generally underestimated compared to experimental values. The present work focuses on the evaluation and refinement of several available all-atom force fields for ACN and proposes a refined flexible six-site potential model. The main modification is related to the reduction of intermolecular LJ parameters (σ and ɛ) for hydrogen atoms, especially ɛ, significantly affecting the dynamical behavior. Besides, the adjustment of σ and ɛ for nitrile carbon and nitrogen atoms helps reach optimum results. Our refined model shows an excellent agreement with the experiment for self-diffusion coefficient and thermodynamic quantities as well as providing a qualitative description of the microscopic structure of liquid ACN. © 2018 Wiley Periodicals, Inc.

Noncanonical α/γ Backbone Conformations in RNA and the Accuracy of Their Description by the AMBER Force Field.

The sugar-phosphate backbone of RNA can exist in diverse rotameric substates, giving RNA molecules enormous conformational variability. The most frequent noncanonical backbone conformation in RNA is α/γ = t/t, which is derived from the canonical backbone by a crankshaft motion and largely preserves the standard geometry of the RNA duplex. A similar conformation also exists in DNA, where it has been extensively studied and shown to be involved in DNA-protein interactions. However, the function of the α/γ = t/t conformation in RNA is poorly understood. Here, we present molecular dynamics simulations of several prototypical RNA structures obtained from X-ray and NMR experiments, including canonical and mismatched RNA duplexes, UUCG and GAGA tetraloops, Loop E, the sarcin-ricin loop, a parallel guanine quadruplex, and a viral pseudoknot. The stability of various noncanonical α/γ backbone conformations was analyzed with two AMBER force fields, ff99bsc0χOL3 and ff99bsc0χOL3 with the recent εζOL1 and βOL1 corrections for DNA. Although some α/γ substates were stable with seemingly well-described equilibria, many were unstable in our simulations. Notably, the most frequent noncanonical conformer α/γ = t/t was unstable in both tested force fields. Possible reasons for this instability are discussed. Our work reveals a potentially important artifact in RNA force fields and highlights a need for further force field refinement.

Combining Simulations and Solution Experiments as a Paradigm for RNA Force Field Refinement.

Recent computational efforts have shown that the current potential energy models used in molecular dynamics are not accurate enough to describe the conformational ensemble of RNA oligomers and suggest that molecular dynamics should be complemented with experimental data. We here propose a scheme based on the maximum entropy principle to combine simulations with bulk experiments. In the proposed scheme, the noise arising from both the measurements and the forward models used to back-calculate the experimental observables is explicitly taken into account. The method is tested on RNA nucleosides and is then used to construct chemically consistent corrections to the Amber RNA force field that allow a large set of experimental data on nucleosides and dinucleosides to be correctly reproduced. The transferability of these corrections is assessed against independent data on tetranucleotides and displays a previously unreported agreement with experiments. This procedure can be applied to enforce multiple experimental data on multiple systems in a self-consistent framework, thus suggesting a new paradigm for force field refinement.

Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides.

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.

Empirical Corrections to the Amber RNA Force Field with Target Metadynamics.

The computational study of conformational transitions in nucleic acids still faces many challenges. For example, in the case of single stranded RNA tetranucleotides, agreement between simulations and experiments is not satisfactory due to inaccuracies in the force fields commonly used in molecular dynamics simulations. We here use experimental data collected from high-resolution X-ray structures to attempt an improvement of the latest version of the AMBER force field. A modified metadynamics algorithm is used to calculate correcting potentials designed to enforce experimental distributions of backbone torsion angles. Replica-exchange simulations of tetranucleotides including these correcting potentials show significantly better agreement with independent solution experiments for the oligonucleotides containing pyrimidine bases. Although the proposed corrections do not seem to be portable to generic RNA systems, the simulations revealed the importance of the α and ζ backbone angles for the modulation of the RNA conformational ensemble. The correction protocol presented here suggests a systematic procedure for force-field refinement.

Role of Electrostatic Interactions in Binding of Peptides and Intrinsically Disordered Proteins to Their Folded Targets: 2. The Model of Encounter Complex Involving the Double Mutant of the c-Crk N-SH3 Domain and Peptide Sos.

In the first part of this work (paper 1, Xue, Y. et al. Biochemistry 2014 , 53 , 6473 ), we have studied the complex between the 10-residue peptide Sos and N-terminal SH3 domain from adaptor protein c-Crk. In the second part (this paper), we designed the double mutant of the c-Crk N-SH3 domain, W169F/Y186L, with the intention to eliminate the interactions responsible for tight peptide-protein binding, while retaining the interactions that create the initial electrostatic encounter complex. The resulting system was characterized experimentally by measuring the backbone and side-chain (15)N relaxation rates, as well as binding shifts and (1)H(N) temperature coefficients. In addition, it was also modeled via a series of ∼5 μs molecular dynamics (MD) simulations recorded in a large water box under an Amber ff99SB*-ILDN force field. Similar to paper 1, we have found that the strength of arginine-aspartate and arginine-glutamate salt bridges is overestimated in the original force field. To address this problem we have applied the empirical force-field correction described in paper 1. Specifically, the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across Arg-to-Asp/Glu salt bridges has been increased by 3%. This modification led to MD models in good agreement with the experimental data. The emerging picture is that of a fuzzy complex, where the peptide "dances" over the surface of the protein, making transient contacts via salt-bridge interactions. Every once in a while the peptide assumes a certain more stable binding pose, assisted by a number of adventitious polar and nonpolar contacts. On the other hand, occasionally Sos flies off the protein surface; it is then guided by electrostatic steering to quickly reconnect with the protein. The dynamic interaction between Sos and the double mutant of c-Crk N-SH3 gives rise to only small binding shifts. The peptide retains a high degree of conformational mobility, although it is appreciably slowed down due to its (loose) association with the protein. Note that spin relaxation data are indispensable in determining the dynamic status of the peptide. Such data can be properly modeled only on a basis of bona fide MD simulations, as shown in our study. We anticipate that in future the field will move away from the ensemble view of protein disorder and toward more sophisticated MD models. This will require further optimization of force fields, aimed specifically at disordered systems. Efforts in this direction have been recently initiated by several research groups; the empirical salt-bridge correction proposed in our work falls in the same category. MD models obtained with the help of suitably refined force fields and guided by experimental NMR data will provide a powerful insight into an intricate world of disordered biomolecules.

RNA Conformational Ensembles: Narrowing the GAP between Experiments and Simulations with Metadynamics

The computational study of conformational transitions in nucleic acids still faces many challenges. For example, in the case of single stranded RNA tetranucleotides, agreement between simulations and experiments is not satisfactory due to inaccuracies in the force fields commonly used in molecular dynamics. Improvement of force fields is however hindered by the difficulties of decoupling those errors from the statistical errors caused by insufficient sampling. We here tackle both problems simultaneously by introducing a novel enhancing sampling method and using experimental data to improve RNA force fields.In this novel method, concurrent well-tempered metadynamics are integrated in a Hamiltonian replica-exchange scheme. The ladder of replicas is built with different strength of the bias potential exploiting the tunability of well-tempered metadynamics so as to scale barriers on individual collective variables [1]. At the same time, the metadynamics algorithm is modified so as to allow enforcing a target distribution of backbone and sugar-base torsion angles taken from experimental structures, using a procedure related to two recently introduced techniques. Replica-exchange simulations of several RNA tetranucleotides with experimental corrections show significantly better agreement with NMR experimental data and suggest a systematic procedure for force field refinement.[1] Gil-Ley, A.; Bussi, G. J. Chem. Theory Comput. 2015, 11, 1077-1085.