deletions | additions       diff --git a/subsection_NMR_calculations_and_protein__.tex b/subsection_NMR_calculations_and_protein__.tex  index 0eb89cc..160118a 100644  --- a/subsection_NMR_calculations_and_protein__.tex  +++ b/subsection_NMR_calculations_and_protein__.tex 
...
\subsection{NMR calculations and protein structures used}  In this paper we benchmark the NMR chemical shift predictions on Ubiquitin and GB3. The structures are geometry optimized using PM6-D3H+ \cite{25024918} using the PCM solvation model and the CHARMM22/CMAP force field \cite{Mackerell_2004} using the GB/SA solvation model \cite{Qiu_1997} using with  the 1UBQ \cite{3041007} and 2OED \cite{15369375} structures as starting points. The PM6-D3H+ optimizations are done using the GAMESS program and a convergence criterion of 5 $\times$ 10$^{-4}$ atomic units, while the CHARMM22/CMAP optimizations are done using TINKER with default convergence criterion 0.01 kcal/mole/angstom. In addition the following NMR-derived structural ensembles are used without further refinement: 1D3Z \cite{Cornilescu_1998}, 2K39 \cite{Lange_2008}, 1XQQ \cite{Lindorff_Larsen_2005}, 2LJ5 \cite{Montalvao_2012}, 2KOX \cite{Fenwick_2011}. In all calculation we used charged protonation states for the acidic and basic side chains including histidine. OBPE/6-31G(d,p)//PM6-D3H+ GIAO NMR shielding calculations are performed with Gaussian09 using the CPCM solvation model. ProCS15 calculations where done using a program specifically written for this paper which can be downloaded at xx. CheShift-2 calculations are performed using the web interface at cheshift.com xx?. CamShift, PPM\_One, Sparta+, shAIC, and ShiftX2 calculations are performed using stand-alone codes downloaded and installed on a local cluster. The NMR chemical shielding and shifts are compared to shifts measured for ubiquitin \cite{Cornilescu_1998} (BMRB ID 17769) and GB3 \cite{V_geli_2012} (BMRB ID 18531), respectively. Much of the variation in some of the chemical shifts comes from the nature of the side-chain itself and the side chains before and after in the sequence, which can lead to inflated R values. To separate the contributions of the sequence and the structure we subtract the measured sequence corrected random coil values \cite{Tamiola_2010} from all predicted and experimental values. Note that this does not affect the computed RMSD values.  
...