Lars Andersen Bratholm edited section_Summary_and_Outlook_In__.tex  over 8 years ago

Commit id: 2791d951af8658d51cfb290d8a9a27d9556f4342

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 and side-chain  dihedral angle changes, dihedral angle-changes angles  of the side-chain a residue  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 of  hydrogen-bonding. ProCS15 reproduces the chemical shielding values computed using PCM/OPBE/6-31G(d,p)//PM6-D3H+ for ubiquitin Ubiquitin  and 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 additivity and the simplified model systems, does  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 Ubiquitin  and GB3 structure. structures.  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 for  C$\alpha$ and C$\beta$ atoms. The RMSD values computed using ProCS for ubiquitin 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. shifts.  ProCS15 is freely available at xx.