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: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