Dataset Information

Benchmarking all-atom simulations using hydrogen exchange.

ABSTRACT: Long-time molecular dynamics (MD) simulations are now able to fold small proteins reversibly to their native structures [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520]. These results indicate that modern force fields can reproduce the energy surface near the native structure. To test how well the force fields recapitulate the other regions of the energy surface, MD trajectories for a variant of protein G are compared with data from site-resolved hydrogen exchange (HX) and other biophysical measurements. Because HX monitors the breaking of individual H-bonds, this experimental technique identifies the stability and H-bond content of excited states, thus enabling quantitative comparison with the simulations. Contrary to experimental findings of a cooperative, all-or-none unfolding process, the simulated denatured state ensemble, on average, is highly collapsed with some transient or persistent native 2° structure. The MD trajectories of this protein G variant and other small proteins exhibit excessive intramolecular H-bonding even for the most expanded conformations, suggesting that the force fields require improvements in describing H-bonding and backbone hydration. Moreover, these comparisons provide a general protocol for validating the ability of simulations to accurately capture rare structural fluctuations.

SUBMITTER: Skinner JJ

PROVIDER: S-EPMC4234613 | biostudies-literature | 2014 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Benchmarking all-atom simulations using hydrogen exchange.

Skinner John J JJ Yu Wookyung W Gichana Elizabeth K EK Baxa Michael C MC Hinshaw James R JR Freed Karl F KF Sosnick Tobin R TR

Proceedings of the National Academy of Sciences of the United States of America 20141027 45

Long-time molecular dynamics (MD) simulations are now able to fold small proteins reversibly to their native structures [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520]. These results indicate that modern force fields can reproduce the energy surface near the native structure. To test how well the force fields recapitulate the other regions of the energy surface, MD trajectories for a variant of protein G are compared with data from site-resolved hydrogen exchange ...[more]

PMID: 25349413

Similar Datasets

Project description:Bimolecular rate constants have been measured for reactions that involve hydrogen atom transfer (HAT) from hydroxylamines to nitroxyl radicals, using the stable radicals TEMPO(*) (2,2,6,6-tetramethylpiperidine-1-oxyl radical), 4-oxo-TEMPO(*) (2,2,6,6-tetramethyl-4-oxo-piperidine-1-oxyl radical), di-tert-butylnitroxyl ((t)Bu(2)NO(*)), and the hydroxylamines TEMPO-H, 4-oxo-TEMPO-H, 4-MeO-TEMPO-H (2,2,6,6-tetramethyl-N-hydroxy-4-methoxy-piperidine), and (t)Bu(2)NOH. The reactions have been monitored by UV-vis stopped-flow methods, using the different optical spectra of the nitroxyl radicals. The HAT reactions all have |DeltaG (o)| < or = 1.4 kcal mol(-1) and therefore are close to self-exchange reactions. The reaction of 4-oxo-TEMPO(*) + TEMPO-H --> 4-oxo-TEMPO-H + TEMPO(*) occurs with k(2H,MeCN) = 10 +/- 1 M(-1) s(-1) in MeCN at 298 K (K(2H,MeCN) = 4.5 +/- 1.8). Surprisingly, the rate constant for the analogous deuterium atom transfer reaction is much slower: k(2D,MeCN) = 0.44 +/- 0.05 M(-1) s(-1) with k(2H,MeCN)/k(2D,MeCN) = 23 +/- 3 at 298 K. The same large kinetic isotope effect (KIE) is found in CH(2)Cl(2), 23 +/- 4, suggesting that the large KIE is not caused by solvent dynamics or hydrogen bonding to solvent. The related reaction of 4-oxo-TEMPO(*) with 4-MeO-TEMPO-H(D) also has a large KIE, k(3H)/k(3D) = 21 +/- 3 in MeCN. For these three reactions, the E(aD) - E(aH) values, between 0.3 +/- 0.6 and 1.3 +/- 0.6 kcal mol(-1), and the log(A(H)/A(D)) values, between 0.5 +/- 0.7 and 1.1 +/- 0.6, indicate that hydrogen tunneling plays an important role. The related reaction of (t)Bu(2)NO(*) + TEMPO-H(D) in MeCN has a large KIE, 16 +/- 3 in MeCN, and very unusual isotopic activation parameters, E(aD) - E(aH) = -2.6 +/- 0.4 and log(A(H)/A(D)) = 3.1 +/- 0.6. Computational studies, using POLYRATE, also indicate substantial tunneling in the (CH(3))(2)NO(*) + (CH(3))(2)NOH model reaction for the experimental self-exchange processes. Additional calculations on TEMPO((*)/H), (t)Bu(2)NO((*)/H), and Ph(2)NO((*)/H) self-exchange reactions reveal why the phenyl groups make the last of these reactions several orders of magnitude faster than the first two. By inference, the calculations also suggest why tunneling appears to be more important in the self-exchange reactions of dialkylhydroxylamines than of arylhydroxylamines.

Project description:Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5-7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Project description:Targeted protein degradation (TPD) has emerged as a powerful approach for removing (rather than inhibiting) proteins implicated in diseases. A key step in TPD is the formation of an induced proximity complex where a degrader molecule recruits an E3 ligase to the protein of interest (POI), facilitating the transfer of ubiquitin to the POI and initiating the proteasomal degradation process. Here, we address three critical aspects of the TPD process using atomistic simulations: 1) formation of the ternary complex induced by a degrader molecule, 2) conformational heterogeneity of the ternary complex, and 3) degradation efficiency via the full Cullin Ring Ligase (CRL) macromolecular assembly. We combine experimental biophysical data---in this case hydrogen-deuterium exchange mass spectrometry (HDX-MS, which measures the solvent exposure of protein residues)---with molecular dynamics (MD) simulations aided by enhanced sampling techniques to accurately predict ternary complex structures at atomic resolution. We demonstrate improved efficiency, accuracy, and reliability in the ternary structure predictions of the bromodomain of the cancer target SMARCA2 with the E3 ligase VHL mediated by three different degrader molecules. The simulations accurately reproduce X-ray crystal structures -- including a new structure that we determined in this work (PDB ID: 7S4E) -- with root mean square deviations (RMSD) of 1.1 to 1.6 \r{A}. The simulations also reveal a structural ensemble of low-energy conformations of the ternary complex. Snapshots from these simulations are used as seeds for additional simulations, where we perform 7.1 milliseconds of aggregate simulation time using Folding@home. The detailed free energy surface captures the crystal structure conformation within a low-energy basin and is consistent with solution-phase experimental data (HDX-MS and SAXS). Finally, we graft a structural ensemble of the ternary complexes onto the full CRL and perform enhanced sampling simulations, which suggest that differences in degradation efficiency may be related to the proximity distribution of lysine residues on the POI relative to the E2-loaded ubiquitin. Several of the predicted ubiquitinated lysine residues have been validated through a ubiquitin mapping proteomics experiment. DOI 10.1038/s41467-022-33575-4

Dataset Information

Benchmarking all-atom simulations using hydrogen exchange.

Publications

Benchmarking all-atom simulations using hydrogen exchange.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets