Jan Jensen edited section_Summary_and_Outlook_In__.tex  over 8 years ago

Commit id: 4ccf5553ab9b669c4a08347a1df674a689b99529

deletions | additions      

       

\section{Summary and Outlook}  In this paper we present ProCS15: a program that computes the isotropic chemical shielding values of backbone atoms and C$\beta$ given a protein structure in less than a second. ProCS accounts for the effect of backbone dihedral angle changes, dihedral angle-changes of the side-chain and the two closest neighbors, hydrogen bonding to the backbone amide group and H$\alpha$ as well as ring-current effects \cite{Christensen_2011} on the hydrogen atoms and assumes that these effects are additive. The backbone, side-chain and hydrogen bonding terms are based on ca 2.35 million OPBE/6-31G(d,p)//PM6 calculations on tripeptides and small structural models hydrogen-bonding.  ProCS15 reproduces the chemical shielding values computed using PCM/OPBE/6-31G(d,p)//PM6-D3H+ for ubiquitin and Protein G GB3  with RMSD values (after linear regression) of up to 2.5 ppm for carbon atoms, 0.9 ppm for hydrogen atoms, and 0.4 ppm for nitrogen. These deviations, which presumably result from the assumption of additivity, doe not appear to preclude equal or better accuracy when comparison to experiment because the accuracies of the chemical shifts computed using ProCS15 (based on linear regression of the chemical shifts, cf. Eq \ref{eqn:scaling}) are very similar to the corresponding DFT calculations using single ubiquitin and Protein G GB3  structure. The largest RMSD values observed for carbon, hydrogen, and nitrogen are, respectively, 2.2 (2.8) ppm, 0.7 (0.6) ppm, and 4.7 (4.6) ppm for ProCS15 (PCM/OPBE/6-31G(d,p)). These accuracies are very similar to DFT-based predictions made by other researchers (e.g. \cite{Zhu_2012},\cite{Zhu_2013}, \cite{Exner_2012}) as well as CheShift-2 \cite{24082119}, which is another DFT-based chemical shift predictor C$\alpha$ and C$\beta$ atoms. The RMSD values computed using ProCS for ubiquitin can be reduced by as much as 0.7, 0.1, and 0.7 ppm for carbon, hydrogen, and nitrogen by using NMR-derived structural ensembles. Similar increase in accuracy is also observed for CheShift-2 (for C$\alpha$ and C$\beta$) while for empirical chemical shift predictors the increase in accuracy is at most 0.3 ppm. The latter observation is another indication that empirical chemical shift predictors are less sensitive to small structural changes, which may make them less suitable for chemical shift-guided refinement of protein structure compared to DFT-based predictors. Christensen and co-workers \cite{24391900} have already demonstrated that this is the case for amide hydrogen bonding geometries using a previous incarnation of ProCS limited to amide proton chemical shift predictions and we are now planning similar refinement studies using all backbone atoms and C$\beta$ chemical sifts. ProCS15 is freely available at xx.