Novel Computational Protein Design Algorithms with Applications to Cystic Fibrosis and HIV

Proteins are essential components of cells and are crucial for catalyzing reactions, signaling, recognition, motility, recycling, and structural stability. This diversity of function suggests that nature is only scratching the surface of protein functional space. Protein function is determined by structure, which in turn is determined predominantly by amino acid sequence. Protein design aims to explore protein sequence and conformational space to design novel proteins with new or improved function. The vast number of possible protein sequences makes exploring the space a challenging problem.

Computational structure-based protein design (CSPD) allows for the rational design of proteins. Because of the large search space, CSPD methods must balance search accuracy and modeling simplifications. We have developed algorithms that allow for the accurate and efficient search of protein conformational space. Specifically, we focus on algorithms that maintain provability, account for protein flexibility, and use ensemble-based rankings. We present several novel algorithms for incorporating improved flexibility into CSPD with continuous rotamers. We applied these algorithms to two biomedically important design problems. We designed peptide inhibitors of the cystic fibrosis agonist CAL that were able to restore function of the vital cystic fibrosis protein CFTR. We also designed improved HIV antibodies and nanobodies to combat HIV infections.

Electrostatic Energetics of Bacillus subtilis Ribonuclease P Protein Determined by Nuclear Magnetic Resonance-Based Histidine pKa Measurements.

The pKa values of ionizable groups in proteins report the free energy of site-specific proton binding and provide a direct means of studying pH-dependent stability. We measured histidine pKa values (H3, H22, and H105) in the unfolded (U), intermediate (I), and sulfate-bound folded (F) states of RNase P protein, using an efficient and accurate nuclear magnetic resonance-monitored titration approach that utilizes internal reference compounds and a parametric fitting method. The three histidines in the sulfate-bound folded protein have pKa values depressed by 0.21 ± 0.01, 0.49 ± 0.01, and 1.00 ± 0.01 units, respectively, relative to that of the model compound N-acetyl-l-histidine methylamide. In the unliganded and unfolded protein, the pKa values are depressed relative to that of the model compound by 0.73 ± 0.02, 0.45 ± 0.02, and 0.68 ± 0.02 units, respectively. Above pH 5.5, H22 displays a separate resonance, which we have assigned to I, whose apparent pKa value is depressed by 1.03 ± 0.25 units, which is ∼0.5 units more than in either U or F. The depressed pKa values we observe are consistent with repulsive interactions between protonated histidine side chains and the net positive charge of the protein. However, the pKa differences between F and U are small for all three histidines, and they have little ionic strength dependence in F. Taken together, these observations suggest that unfavorable electrostatics alone do not account for the fact that RNase P protein is intrinsically unfolded in the absence of ligand. Multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Mass spectrometry-based thermal shift assay for protein-ligand binding analysis.

Described here is a mass spectrometry-based screening assay for the detection of protein-ligand binding interactions in multicomponent protein mixtures. The assay utilizes an oxidation labeling protocol that involves using hydrogen peroxide to selectively oxidize methionine residues in proteins in order to probe the solvent accessibility of these residues as a function of temperature. The extent to which methionine residues in a protein are oxidized after specified reaction times at a range of temperatures is determined in a MALDI analysis of the intact proteins and/or an LC-MS analysis of tryptic peptide fragments generated after the oxidation reaction is quenched. Ultimately, the mass spectral data is used to construct thermal denaturation curves for the detected proteins. In this proof-of-principle work, the protocol is applied to a four-protein model mixture comprised of ubiquitin, ribonuclease A (RNaseA), cyclophilin A (CypA), and bovine carbonic anhydrase II (BCAII). The new protocol's ability to detect protein-ligand binding interactions by comparing thermal denaturation data obtained in the absence and in the presence of ligand is demonstrated using cyclosporin A (CsA) as a test ligand. The known binding interaction between CsA and CypA was detected using both the MALDI- and LC-MS-based readouts described here.

Interaction of Cryptococcus neoformans Rim101 and protein kinase A regulates capsule.

Cryptococcus neoformans is a prevalent human fungal pathogen that must survive within various tissues in order to establish a human infection. We have identified the C. neoformans Rim101 transcription factor, a highly conserved pH-response regulator in many fungal species. The rim101 multiply sign in circle mutant strain displays growth defects similar to other fungal species in the presence of alkaline pH, increased salt concentrations, and iron limitation. However, the rim101 multiply sign in circle strain is also characterized by a striking defect in capsule, an important virulence-associated phenotype. This capsular defect is likely due to alterations in polysaccharide attachment to the cell surface, not in polysaccharide biosynthesis. In contrast to many other C. neoformans capsule-defective strains, the rim101 multiply sign in circle mutant is hypervirulent in animal models of cryptococcosis. Whereas Rim101 activation in other fungal species occurs through the conserved Rim pathway, we demonstrate that C. neoformans Rim101 is also activated by the cAMP/PKA pathway. We report here that C. neoformans uses PKA and the Rim pathway to regulate the localization, activation, and processing of the Rim101 transcription factor. We also demonstrate specific host-relevant activating conditions for Rim101 cleavage, showing that C. neoformans has co-opted conserved signaling pathways to respond to the specific niche within the infected host. These results establish a novel mechanism for Rim101 activation and the integration of two conserved signaling cascades in response to host environmental conditions.

The G-protein-coupled receptor kinases beta ARK1 and beta ARK2 are widely distributed at synapses in rat brain.

The beta-adrenergic receptor kinase (beta ARK) phosphorylates the agonist-occupied beta-adrenergic receptor to promote rapid receptor uncoupling from Gs, thereby attenuating adenylyl cyclase activity. Beta ARK-mediated receptor desensitization may reflect a general molecular mechanism operative on many G-protein-coupled receptor systems and, particularly, synaptic neurotransmitter receptors. Two distinct cDNAs encoding beta ARK isozymes were isolated from rat brain and sequenced. The regional and cellular distributions of these two gene products, termed beta ARK1 and beta ARK2, were determined in brain by in situ hybridization and by immunohistochemistry at the light and electron microscopic levels. The beta ARK isozymes were found to be expressed primarily in neurons distributed throughout the CNS. Ultrastructurally, beta ARK1 and beta ARK2 immunoreactivities were present both in association with postsynaptic densities and, presynaptically, with axon terminals. The beta ARK isozymes have a regional and subcellular distribution consistent with a general role in the desensitization of synaptic receptors.

Independent beta-arrestin 2 and G protein-mediated pathways for angiotensin II activation of extracellular signal-regulated kinases 1 and 2.

Stimulation of a mutant angiotensin type 1A receptor (DRY/AAY) with angiotensin II (Ang II) or of a wild-type receptor with an Ang II analog ([sarcosine1,Ile4,Ile8]Ang II) fails to activate classical heterotrimeric G protein signaling but does lead to recruitment of beta-arrestin 2-GFP and activation of extracellular signal-regulated kinases 1 and 2 (ERK1/2) (maximum stimulation approximately 50% of wild type). This G protein-independent activation of mitogen-activated protein kinase is abolished by depletion of cellular beta-arrestin 2 but is unaffected by the PKC inhibitor Ro-31-8425. In parallel, stimulation of the wild-type angiotensin type 1A receptor with Ang II robustly stimulates ERK1/2 activation with approximately 60% of the response blocked by the PKC inhibitor (G protein dependent) and the rest of the response blocked by depletion of cellular beta-arrestin 2 by small interfering RNA (beta-arrestin dependent). These findings imply the existence of independent G protein- and beta-arrestin 2-mediated pathways leading to ERK1/2 activation and the existence of distinct "active" conformations of a seven-membrane-spanning receptor coupled to each.

G protein signaling and vein graft intimal hyperplasia: reduction of intimal hyperplasia in vein grafts by a Gbetagamma inhibitor suggests a major role of G protein signaling in lesion development.

Vein grafting results in the development of intimal hyperplasia with accompanying changes in guanine nucleotide-binding (G) protein expression and function. Several serum mitogens that act through G protein-coupled receptors, such as lysophosphatidic acid, stimulate proliferative pathways that are dependent on the G protein betagamma subunit (Gbetagamma)-mediated activation of p21ras. This study examines the role of Gbetagamma signaling in intimal hyperplasia by targeting a gene encoding a specific Gbetagamma inhibitor in an experimental rabbit vein graft model. This inhibitor, the carboxyl terminus of the beta-adrenergic receptor kinase (betaARK(CT)), contains a Gbetagamma-binding domain. Vein graft intimal hyperplasia was significantly reduced by 37% (P<0.01), and physiological studies demonstrated that the normal alterations in G protein coupling phenotypically seen in this model were blocked by betaARK(CT) treatment. Thus, it appears that Gbetagamma-mediated pathways play a major role in intimal hyperplasia and that targeting inhibitors of Gbetagamma signaling offers novel intraoperative therapeutic modalities to inhibit the development of vein graft intimal hyperplasia and subsequent vein graft failure.

beta2-Adrenergic receptor regulation by GIT1, a G protein-coupled receptor kinase-associated ADP ribosylation factor GTPase-activating protein.

G protein-coupled receptor activation leads to the membrane recruitment and activation of G protein-coupled receptor kinases, which phosphorylate receptors and lead to their inactivation. We have identified a novel G protein-coupled receptor kinase-interacting protein, GIT1, that is a GTPase-activating protein (GAP) for the ADP ribosylation factor (ARF) family of small GTP-binding proteins. Overexpression of GIT1 leads to reduced beta2-adrenergic receptor signaling and increased receptor phosphorylation, which result from reduced receptor internalization and resensitization. These cellular effects of GIT1 require its intact ARF GAP activity and do not reflect regulation of GRK kinase activity. These results suggest an essential role for ARF proteins in regulating beta2-adrenergic receptor endocytosis. Moreover, they provide a mechanism for integration of receptor activation and endocytosis through regulation of ARF protein activation by GRK-mediated recruitment of the GIT1 ARF GAP to the plasma membrane.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.

cDNA for the human beta 2-adrenergic receptor: a protein with multiple membrane-spanning domains and encoded by a gene whose chromosomal location is shared with that of the receptor for platelet-derived growth factor.

We have isolated and sequenced a cDNA encoding the human beta 2-adrenergic receptor. The deduced amino acid sequence (413 residues) is that of a protein containing seven clusters of hydrophobic amino acids suggestive of membrane-spanning domains. While the protein is 87% identical overall with the previously cloned hamster beta 2-adrenergic receptor, the most highly conserved regions are the putative transmembrane helices (95% identical) and cytoplasmic loops (93% identical), suggesting that these regions of the molecule harbor important functional domains. Several of the transmembrane helices also share lesser degrees of identity with comparable regions of select members of the opsin family of visual pigments. We have localized the gene for the beta 2-adrenergic receptor to q31-q32 on chromosome 5. This is the same position recently determined for the gene encoding the receptor for platelet-derived growth factor and is adjacent to that for the FMS protooncogene, which encodes the receptor for the macrophage colony-stimulating factor.

Wee1-regulated apoptosis mediated by the crk adaptor protein in Xenopus egg extracts.

Many of the biochemical reactions of apoptotic cell death, including mitochondrial cytochrome c release and caspase activation, can be reconstituted in cell-free extracts derived from Xenopus eggs. In addition, because caspase activation does not occur until the egg extract has been incubated for several hours on the bench, upstream signaling processes occurring before full apoptosis are rendered accessible to biochemical manipulation. We reported previously that the adaptor protein Crk is required for apoptotic signaling in egg extracts (Evans, E.K., W. Lu, S.L. Strum, B.J. Mayer, and S. Kornbluth. 1997. EMBO (Eur. Mol. Biol. Organ.) J. 16:230-241). Moreover, we demonstrated that removal of Crk Src homology (SH)2 or SH3 interactors from the extracts prevented apoptosis. We now report the finding that the relevant Crk SH2-interacting protein, important for apoptotic signaling in the extract, is the well-known cell cycle regulator, Wee1. We have demonstrated a specific interaction between tyrosine-phosphorylated Wee1 and the Crk SH2 domain and have shown that recombinant Wee1 can restore apoptosis to an extract depleted of SH2 interactors. Moreover, exogenous Wee1 accelerated apoptosis in egg extracts, and this acceleration was largely dependent on the presence of endogenous Crk protein. As other Cdk inhibitors, such as roscovitine and Myt1, did not act like Wee1 to accelerate apoptosis, we propose that Wee1-Crk complexes signal in a novel apoptotic pathway, which may be unrelated to Wee1's role as a cell cycle regulator.

Suppression of conformational heterogeneity at a protein-protein interface.

Staphylococcal protein A (SpA) is an important virulence factor from Staphylococcus aureus responsible for the bacterium's evasion of the host immune system. SpA includes five small three-helix-bundle domains that can each bind with high affinity to many host proteins such as antibodies. The interaction between a SpA domain and the Fc fragment of IgG was partially elucidated previously in the crystal structure 1FC2. Although informative, the previous structure was not properly folded and left many substantial questions unanswered, such as a detailed description of the tertiary structure of SpA domains in complex with Fc and the structural changes that take place upon binding. Here we report the 2.3-Å structure of a fully folded SpA domain in complex with Fc. Our structure indicates that there are extensive structural rearrangements necessary for binding Fc, including a general reduction in SpA conformational heterogeneity, freezing out of polyrotameric interfacial residues, and displacement of a SpA side chain by an Fc side chain in a molecular-recognition pocket. Such a loss of conformational heterogeneity upon formation of the protein-protein interface may occur when SpA binds its multiple binding partners. Suppression of conformational heterogeneity may be an important structural paradigm in functionally plastic proteins.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.