Jonathan H. Sheehan

Practically Useful: What the Rosetta Protein Modeling Suite Can Do for You

The objective of this review is to enable researchers to use the software package Rosetta for biochemical and biomedicinal studies. We provide a brief review of the six most frequent research problems tackled with Rosetta. For each of these six tasks, we provide a tutorial that illustrates a basic Rosetta protocol. The Rosetta method was originally developed for de novo protein structure prediction and is regularly one of the best performers in the community-wide biennial Critical Assessment of Structure Prediction. Predictions for protein domains with fewer than 125 amino acids regularly have a backbone root-mean-square deviation of better than 5.0 Å. More impressively, there are several cases in which Rosetta has been used to predict structures with atomic level accuracy better than 2.5 Å. In addition to de novo structure prediction, Rosetta also has methods for molecular docking, homology modeling, determining protein structures from sparse experimental NMR or EPR data, and protein design. Rosetta has been used to accurately design a novel protein structure, predict the structure of protein−protein complexes, design altered specificity protein−protein and protein−DNA interactions, and stabilize proteins and protein complexes. Most recently, Rosetta has been used to solve the X-ray crystallographic phase problem.

Small-molecule ligand docking into comparative models with Rosetta

Structure-based drug design is frequently used to accelerate the development of small-molecule therapeutics. Although substantial progress has been made in X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy, the availability of high-resolution structures is limited owing to the frequent inability to crystallize or obtain sufficient NMR restraints for large or flexible proteins. Computational methods can be used to both predict unknown protein structures and model ligand interactions when experimental data are unavailable. This paper describes a comprehensive and detailed protocol using the Rosetta modeling suite to dock small-molecule ligands into comparative models. In the protocol presented here, we review the comparative modeling process, including sequence alignment, threading and loop building. Next, we cover docking a small-molecule ligand into the protein comparative model. In addition, we discuss criteria that can improve ligand docking into comparative models. Finally, and importantly, we present a strategy for assessing model quality. The entire protocol is presented on a single example selected solely for didactic purposes. The results are therefore not representative and do not replace benchmarks published elsewhere. We also provide an additional tutorial so that the user can gain hands-on experience in using Rosetta. The protocol should take 5–7 h, with additional time allocated for computer generation of models.

The Mode of Action of Centrin BINDING OF Ca2+ AND A PEPTIDE FRAGMENT OF Kar1p TO THE C-TERMINAL DOMAIN

Centrin is an EF-hand calcium-binding protein closely related to the prototypical calcium sensor protein calmodulin. It is found in microtubule-organizing centers of organisms ranging from algae and yeast to man. In vitro, the C-terminal domain of centrin binds to the yeast centrosomal protein Kar1p in a calcium-dependent manner, whereas the N-terminal domain does not show any appreciable affinity for Kar1p. To obtain deeper insights into the structural basis for centrins function, we have characterized the affinities of the C-terminal domain of Chlamydomonas reinhardtii centrin for calcium and for a peptide fragment of Kar1p using CD, fluorescence, and NMR spectroscopy. Calcium binding site IV in C. reinhardtii centrin was found to bind Ca2+ ∼100-fold more strongly than site III. In the absence of Ca2+, the protein occupies a mixture of closed conformations. Binding of a single ion in site IV is sufficient to radically alter the conformational equilibrium, promoting occupancy of an open conformation. However, an exchange between closed and open conformations remains even at saturating levels of Ca2+. The population of the open conformation is substantially stabilized by the presence of the target peptide Kar1p-(239–257) to a point where a single ion bound in site IV is sufficient to completely shift the conformational equilibrium to the open conformation. This is reflected in the enhancement of the Ca2+ affinity in this site by more than an order of magnitude. These data confirm the direct coupling of the Ca2+ binding-induced shift in the equilibrium between the closed and open conformations to the binding of the peptide. Combined with the common localization of the two proteins in the microtubule organizing center, our results suggest that centrin is constitutively bound to Kar1p through its C-terminal domain and that centrins calcium sensor activities are mediated by the N-terminal domain.

Glycosylation of β2 Subunits Regulates GABAA Receptor Biogenesis and Channel Gating

γ-Aminobutyric acid type A (GABAA) receptors are heteropentameric glycoproteins. Based on consensus sequences, the GABAA receptor β2 subunit contains three potential N-linked glycosylation sites, Asn-32, Asn-104, and Asn-173. Homology modeling indicates that Asn-32 and Asn-104 are located before the α1 helix and in loop L3, respectively, near the top of the subunit-subunit interface on the minus side, and that Asn-173 is located in the Cys-loop near the bottom of the subunit N-terminal domain. Using site-directed mutagenesis, we demonstrated that all predicted β2 subunit glycosylation sites were glycosylated in transfected HEK293T cells. Glycosylation of each site, however, produced specific changes in α1β2 receptor surface expression and function. Although glycosylation of Asn-173 in the Cys-loop was important for stability of β2 subunits when expressed alone, results obtained with flow cytometry, brefeldin A treatment, and endo-β-N-acetylglucosaminidase H digestion suggested that glycosylation of Asn-104 was required for efficient α1β2 receptor assembly and/or stability in the endoplasmic reticulum. Patch clamp recording revealed that mutation of each site to prevent glycosylation decreased peak α1β2 receptor current amplitudes and altered the gating properties of α1β2 receptor channels by reducing mean open time due to a reduction in the proportion of long open states. In addition to functional heterogeneity, endo-β-N-acetylglucosaminidase H digestion and glycomic profiling revealed that surface β2 subunit N-glycans at Asn-173 were high mannose forms that were different from those of Asn-32 and N104. Using a homology model of the pentameric extracellular domain of α1β2 channel, we propose mechanisms for regulation of GABAA receptors by glycosylation.

NFI transcription factors interact with FOXA1 to regulate prostate-specific gene expression

Androgen receptor (AR) action throughout prostate development and in maintenance of the prostatic epithelium is partly controlled by interactions between AR and forkhead box (FOX) transcription factors, particularly FOXA1. We sought to identity additional FOXA1 binding partners that may mediate prostate-specific gene expression. Here we identify the nuclear factor I (NFI) family of transcription factors as novel FOXA1 binding proteins. All four family members (NFIA, NFIB, NFIC, and NFIX) can interact with FOXA1, and knockdown studies in androgen-dependent LNCaP cells determined that modulating expression of NFI family members results in changes in AR target gene expression. This effect is probably mediated by binding of NFI family members to AR target gene promoters, because chromatin immunoprecipitation (ChIP) studies found that NFIB bound to the prostate-specific antigen enhancer. Förster resonance energy transfer studies revealed that FOXA1 is capable of bringing AR and NFIX into proximity, indicating that FOXA1 facilitates the AR and NFI interaction by bridging the complex. To determine the extent to which NFI family members regulate AR/FOXA1 target genes, motif analysis of publicly available data for ChIP followed by sequencing was undertaken. This analysis revealed that 34.4% of peaks bound by AR and FOXA1 contain NFI binding sites. Validation of 8 of these peaks by ChIP revealed that NFI family members can bind 6 of these predicted genomic elements, and 4 of the 8 associated genes undergo gene expression changes as a result of individual NFI knockdown. These observations suggest that NFI regulation of FOXA1/AR action is a frequent event, with individual family members playing distinct roles in AR target gene expression.

Protocols for Molecular Modeling with Rosetta3 and RosettaScripts

Previously, we published an article providing an overview of the Rosetta suite of biomacromolecular modeling software and a series of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987–2998]. The overwhelming positive response to this publication we received motivates us to here share the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking, small molecule docking, and protein design. This updated and expanded set of tutorials is needed, as since 2010 Rosetta has been fully redesigned into an object-oriented protein modeling program Rosetta3. Notable improvements include a substantially improved energy function, an XML-like language termed “RosettaScripts” for flexibly specifying modeling task, new analysis tools, the addition of the TopologyBroker to control conformational sampling, and support for multiple templates in comparative modeling. Rosetta’s ability to model systems with symmetric proteins, membrane proteins, noncanonical amino acids, and RNA has also been greatly expanded and improved.

Structure of the N-terminal Calcium Sensor Domain of Centrin Reveals the Biochemical Basis for Domain-specific Function

Centrin is an essential component of microtubule-organizing centers in organisms ranging from algae and yeast to humans. It is an EF-hand calcium-binding protein with homology to calmodulin but distinct calcium binding properties. In a previously proposed model, the C-terminal domain of centrin serves as a constitutive anchor to target proteins, and the N-terminal domain serves as the sensor of calcium signals. The three-dimensional structure of the N-terminal domain of Chlamydomonas rheinhardtii centrin has been determined in the presence of calcium by solution NMR spectroscopy. The domain is found to occupy an open conformation typical of EF-hand calcium sensors. Comparison of the N- and C-terminal domains of centrin reveals a structural and biochemical basis for the domain specificity of interactions with its cellular targets and the distinct nature of centrin relative to other EF-hand proteins. An NMR titration of the centrin N-terminal domain with a fragment of the known centrin target Sfi1 reveals binding of the peptide to a discrete site on the protein, which supports the proposal that the N-terminal domain serves as a calcium sensor in centrin.

EGFR Kinase Domain Duplication (EGFR-KDD) Is a Novel Oncogenic Driver in Lung Cancer That Is Clinically Responsive to Afatinib

UNLABELLED Oncogenic EGFR mutations are found in 10% to 35% of lung adenocarcinomas. Such mutations, which present most commonly as small in-frame deletions in exon 19 or point mutations in exon 21 (L858R), confer sensitivity to EGFR tyrosine kinase inhibitors (TKI). In analyzing the tumor from a 33-year-old male never-smoker, we identified a novel EGFR alteration in lung cancer: EGFR exon 18-25 kinase domain duplication (EGFR-KDD). Through analysis of a larger cohort of tumor samples, we detected additional cases of EGFR-KDD in lung, brain, and other cancers. In vitro, EGFR-KDD is constitutively active, and computational modeling provides potential mechanistic support for its auto-activation. EGFR-KDD-transformed cells are sensitive to EGFR TKIs and, consistent with these in vitro findings, the index patient had a partial response to the EGFR TKI afatinib. The patient eventually progressed, at which time resequencing revealed an EGFR-dependent mechanism of acquired resistance to afatinib, thereby validating EGFR-KDD as a driver alteration and therapeutic target. SIGNIFICANCE We identified oncogenic and drug-sensitive EGFR-KDD that is recurrent in lung, brain, and soft-tissue cancers and documented that a patient with metastatic lung adenocarcinoma harboring the EGFR-KDD derived significant antitumor response from treatment with the EGFR inhibitor afatinib. Findings from these studies will be immediately translatable, as there are already several approved EGFR inhibitors in clinical use.

Second-site suppressors of HIV-1 capsid mutations: restoration of intracellular activities without correction of intrinsic capsid stability defects

BackgroundDisassembly of the viral capsid following penetration into the cytoplasm, or uncoating, is a poorly understood stage of retrovirus infection. Based on previous studies of HIV-1 CA mutants exhibiting altered capsid stability, we concluded that formation of a capsid of optimal intrinsic stability is crucial for HIV-1 infection.ResultsTo further examine the connection between HIV-1 capsid stability and infectivity, we isolated second-site suppressors of HIV-1 mutants exhibiting unstable (P38A) or hyperstable (E45A) capsids. We identified the respective suppressor mutations, T216I and R132T, which restored virus replication in a human T cell line and markedly enhanced the fitness of the original mutants as revealed in single-cycle infection assays. Analysis of the corresponding purified N-terminal domain CA proteins by NMR spectroscopy demonstrated that the E45A and R132T mutations induced structural changes that are localized to the regions of the mutations, while the P38A mutation resulted in changes extending to neighboring regions in space. Unexpectedly, neither suppressor mutation corrected the intrinsic viral capsid stability defect associated with the respective original mutation. Nonetheless, the R132T mutation rescued the selective infectivity impairment exhibited by the E45A mutant in aphidicolin-arrested cells, and the double mutant regained sensitivity to the small molecule inhibitor PF74. The T216I mutation rescued the impaired ability of the P38A mutant virus to abrogate restriction by TRIMCyp and TRIM5α.ConclusionsThe second-site suppressor mutations in CA that we have identified rescue virus infection without correcting the intrinsic capsid stability defects associated with the P38A and E45A mutations. The suppressors also restored wild type virus function in several cell-based assays. We propose that while proper HIV-1 uncoating in target cells is dependent on the intrinsic stability of the viral capsid, the effects of stability-altering mutations can be mitigated by additional mutations that affect interactions with host factors in target cells or the consequences of these interactions. The ability of mutations at other CA surfaces to compensate for effects at the NTD-NTD interface further indicates that uncoating in target cells is controlled by multiple intersubunit interfaces in the viral capsid.

Measurement of Aptamer–Protein Interactions with Back-Scattering Interferometry

We report the quantitative measurement of aptamer-protein interactions using backscattering interferometry (BSI) and show that BSI can determine when distinct binding regions are accessed. As a model system, we utilized two DNA aptamers (Tasset and Bock) that bind to distinct sites of a target protein (human α-thrombin). This is the first time BSI has been used to study a multivalent system in free solution wherein more than one ligand binds to a single target. We measured aptamer equilibrum dissociation constants (K(d)) of 3.84 nM (Tasset-thrombin) and 5.96 nM (Bock-thrombin), in close agreement with the literature. Unexpectedly, we observed allosteric effects such that the binding of the first aptamer resulted in a significant change in the binding affinity of the second aptamer. For example, the K(d) of Bock aptamer binding to preformed Tasset-thrombin complexes was 7-fold lower (indicating higher affinity) compared to binding to thrombin alone. Preliminary modeling efforts suggest evidence for allosteric linkage between the two exosites.