Sparse pseudocontact shift NMR data obtained from a non-canonical amino acid-linked lanthanide tag improves integral membrane protein structure prediction.

A single experimental method alone often fails to provide the resolution, accuracy, and coverage needed to model integral membrane proteins (IMPs). Integrating computation with experimental data is a powerful approach to supplement missing structural information with atomic detail. We combine RosettaNMR with experimentally-derived paramagnetic NMR restraints to guide membrane protein structure prediction. We demonstrate this approach using the disulfide bond formation protein B (DsbB), an α-helical IMP. Here, we attached a cyclen-based paramagnetic lanthanide tag to an engineered non-canonical amino acid (ncAA) using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry reaction. Using this tagging strategy, we collected 203 backbone HN pseudocontact shifts (PCSs) for three different labeling sites and used these as input to guide de novo membrane protein structure prediction protocols in Rosetta. We find that this sparse PCS dataset combined with 44 long-range NOEs as restraints in our calculations improves structure prediction of DsbB by enhancements in model accuracy, sampling, and scoring. The inclusion of this PCS dataset improved the Cα-RMSD transmembrane segment values of the best-scoring and best-RMSD models from 9.57 Å and 3.06 Å (no NMR data) to 5.73 Å and 2.18 Å, respectively.

Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery.

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 µM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.