Three-Dimensional Interaction Homology: Deconstructing Residue-Residue and Residue-Lipid Interactions in Membrane Proteins.

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

ATM phosphorylates PP2A subunit A resulting in nuclear export and spatiotemporal regulation of the DNA damage response.

Ataxia telangiectasia mutated (ATM) is a serine-threonine protein kinase and important regulator of the DNA damage response (DDR). One critical ATM target is the structural subunit A (PR65-S401) of protein phosphatase 2A (PP2A), known to regulate diverse cellular processes such as mitosis and cell growth as well as dephosphorylating many proteins during the recovery from the DDR. We generated mouse embryonic fibroblasts expressing PR65-WT, -S401A (cannot be phosphorylated), and -S401D (phospho-mimetic) transgenes. Significantly, S401 mutants exhibited extensive chromosomal aberrations, impaired DNA double-strand break (DSB) repair and underwent increased mitotic catastrophe after radiation. Both S401A and the S401D cells showed impaired DSB repair (nonhomologous end joining and homologous recombination repair) and exhibited delayed DNA damage recovery, which was reflected in reduced radiation survival. Furthermore, S401D cells displayed increased ERK and AKT signaling resulting in enhanced growth rate further underscoring the multiple roles ATM-PP2A signaling plays in regulating prosurvival responses. Time-lapse video and cellular localization experiments showed that PR65 was exported to the cytoplasm after radiation by CRM1, a nuclear export protein, in line with the very rapid pleiotropic effects observed. A putative nuclear export sequence (NES) close to S401 was identified and when mutated resulted in aberrant PR65 shuttling. Our study demonstrates that the phosphorylation of a single, critical PR65 amino acid (S401) by ATM fundamentally controls the DDR, and balances DSB repair quality, cell survival and growth by spatiotemporal PR65 nuclear-cytoplasmic shuttling mediated by the nuclear export receptor CRM1.

Novel eIF4A1 inhibitors with anti-tumor activity in lymphoma.

BACKGROUND: Deregulated translation initiation is implicated extensively in cancer initiation and progression. It is actively pursued as a viable target that circumvents the dependency on oncogenic signaling, a significant factor in current strategies. Eukaryotic translation initiation factor (eIF) 4A plays an essential role in translation initiation by unwinding the secondary structure of messenger RNA (mRNA) upstream of the start codon, enabling active ribosomal recruitment on the downstream genes. Several natural product molecules with similar scaffolds, such as Rocaglamide A (RocA), targeting eIF4A have been reported in the last decade. However, their clinical utilization is still elusive due to several pharmacological limitations. In this study we identified new eIF4A1 inhibitors and their possible mechanisms. METHODS: In this report, we conducted a pharmacophore-based virtual screen of RocA complexed with eIF4A and a polypurine RNA strand for novel eIF4A inhibitors from commercially available compounds in the MolPort Database. We performed target-based screening and optimization of active pharmacophores. We assessed the effects of novel compounds on biochemical and cell-based assays for efficacy and mechanistic evaluation. RESULTS: We validated three new potent eIF4A inhibitors, RBF197, RBF 203, and RBF 208, which decreased diffuse large B-cell lymphoma (DLBCL) cell viability. Biochemical and cellular studies, molecular docking, and functional assays revealed that thosenovel compounds clamp eIF4A into mRNA in an ATP-independent manner. Moreover, we found that RBF197 and RBF208 significantly depressed eIF4A-dependent oncogene expression as well as the colony formation capacity of DLBCL. Interestingly, exposure of these compounds to non-malignant cells had only minimal impact on their growth and viability. CONCLUSIONS: Identified compounds suggest a new strategy for designing novel eIF4A inhibitors.

From oncoproteins to spike proteins: the evaluation of intramolecular stability using hydropathic force field.

Evaluation of the intramolecular stability of proteins plays a key role in the comprehension of their biological behavior and mechanism of action. Small structural alterations such as mutations induced by single nucleotide polymorphism can impact biological activity and pharmacological modulation. Covid-19 mutations, that affect viral replication and the susceptibility to antibody neutralization, and the action of antiviral drugs, are just one example. In this work, the intramolecular stability of mutated proteins, like Spike glycoprotein and its complexes with the human target, is evaluated through hydropathic intramolecular energy scoring originally conceived by Abraham and Kellogg based on the "Extension of the fragment method to calculate amino acid zwitterion and side-chain partition coefficients" by Abraham and Leo in Proteins: Struct. Funct. Genet. 1987, 2:130 - 52. HINT is proposed as a fast and reliable tool for the stability evaluation of any mutated system. This work has been written in honor of Prof. Donald J. Abraham (1936-2021).

FOSL1 promotes metastasis of head and neck squamous cell carcinoma through super-enhancer-driven transcription program.

Previously, we discovered that FOSL1 facilitates the metastasis of head and neck squamous cell carcinoma (HNSCC) cancer stem cells in a spontaneous mouse model. However, the molecular mechanisms remained unclear. Here, we demonstrated that FOSL1 serves as the dominant activating protein 1 (AP1) family member and is significantly upregulated in HNSCC tumor tissues and correlated with metastasis of HNSCC. Mechanistically, FOSL1 exerts its function in promoting tumorigenicity and metastasis predominantly via selective association with Mediators to establish super-enhancers (SEs) at a cohort of cancer stemness and pro-metastatic genes, such as SNAI2 and FOSL1 itself. Depletion of FOSL1 led to disruption of SEs and expression inhibition of these key oncogenes, which resulted in the suppression of tumor initiation and metastasis. We also revealed that the abundance of FOSL1 is positively associated with the abundance of SNAI2 in HNSCC and the high expression levels of FOSL1 and SNAI2 are associated with short overall disease-free survival. Finally, the administration of the FOSL1 inhibitor SR11302 significantly suppressed tumor growth and lymph node metastasis of HNSCC in a patient-derived xenograft model. These findings indicate that FOSL1 is a master regulator that promotes the metastasis of HNSCC through a SE-driven transcription program that may represent an attractive target for therapeutic interventions.

Elucidating the Interaction between Pyridoxine 5'-Phosphate Oxidase and Dopa Decarboxylase: Activation of B6-Dependent Enzyme.

Pyridoxal 5'-phosphate (PLP), the active form of vitamin B6, serves as a cofactor for scores of B6-dependent (PLP-dependent) enzymes involved in many cellular processes. One such B6 enzyme is dopa decarboxylase (DDC), which is required for the biosynthesis of key neurotransmitters, e.g., dopamine and serotonin. PLP-dependent enzymes are biosynthesized as apo-B6 enzymes and then converted to the catalytically active holo-B6 enzymes by Schiff base formation between the aldehyde of PLP and an active site lysine of the protein. In eukaryotes, PLP is made available to the B6 enzymes through the activity of the B6-salvage enzymes, pyridoxine 5'-phosphate oxidase (PNPO) and pyridoxal kinase (PLK). To minimize toxicity, the cell keeps the content of free PLP (unbound) very low through dephosphorylation and PLP feedback inhibition of PNPO and PLK. This has led to a proposed mechanism of complex formation between the B6-salvage enzymes and apo-B6 enzymes prior to the transfer of PLP, although such complexes are yet to be characterized at the atomic level, presumably due to their transient nature. A computational study, for the first time, was used to predict a likely PNPO and DDC complex, which suggested contact between the allosteric PLP tight-binding site on PNPO and the active site of DDC. Using isothermal calorimetry and/or surface plasmon resonance, we also show that PNPO binds both apoDDC and holoDDC with dissociation constants of 0.93 ± 0.07 µM and 2.59 ± 0.11 µM, respectively. Finally, in the presence of apoDDC, the tightly bound PLP on PNPO is transferred to apoDDC, resulting in the formation of about 35% holoDDC.

Structural understanding of 5-(4-hydroxy-phenyl)-N-(2-(5-methoxy-1H-indol-3-yl)-ethyl)-3-oxopentanamide as a neuroprotectant for Alzheimer's disease.

In our continuing efforts to develop novel neuroprotectants for Alzheimer's disease (AD), a series of analogs based on a lead compound that was recently shown to target the mitochondrial complex I were designed, synthesized and biologically characterized to understand the structure features that are important for neuroprotective activities. The results from a cellular AD model highlighted the important roles of the 4-OH on the phenyl ring and the 5-OCH3 on the indole ring of the lead compound. The results also demonstrated that the ß-keto moiety can be modified to retain or improve the neuroprotective activity. Docking studies of selected analogs to the FMN site of mitochondrial complex I also supported the observed neuroprotective activities. Collectively, the results provide further information to guide optimization and development of analogs based on this chemical scaffold as neuroprotectants with a novel mechanism of action for AD.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues.

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ1, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the ß-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Functional diversification of the chemical landscapes of yeast Sec14-like phosphatidylinositol transfer protein lipid-binding cavities.

Phosphatidylinositol-transfer proteins (PITPs) are key regulators of lipid signaling in eukaryotic cells. These proteins both potentiate the activities of phosphatidylinositol (PtdIns) 4-OH kinases and help channel production of specific pools of phosphatidylinositol 4-phosphate (PtdIns(4)P) dedicated to specific biological outcomes. In this manner, PITPs represent a major contributor to the mechanisms by which the biological outcomes of phosphoinositide are diversified. The two-ligand priming model proposes that the engine by which Sec14-like PITPs potentiate PtdIns kinase activities is a heterotypic lipid-exchange cycle where PtdIns is a common exchange substrate among the Sec14-like PITP family, but the second exchange ligand varies with the PITP. A major prediction of this model is that second-exchangeable ligand identity will vary from PITP to PITP. To address the heterogeneity in the second exchange ligand for Sec14-like PITPs, we used structural, computational, and biochemical approaches to probe the diversities of the lipid-binding cavity microenvironments of the yeast Sec14-like PITPs. The collective data report that yeast Sec14-like PITP lipid-binding pockets indeed define diverse chemical microenvironments that translate into differential ligand-binding specificities across this protein family.

Synthesis of novel 5-substituted-2-aminotetralin analogs: 5-HT1A and 5-HT7 G protein-coupled receptor affinity, 3D-QSAR and molecular modeling.

The serotonin 5-HT7 G protein-coupled receptor (GPCR) is a proposed pharmacotherapeutic target for a variety of central and peripheral indications, albeit, there are no approved drugs selective for binding 5-HT7. We previously reported that a lead analog based on the 5-substituted-N,N-disubstituted-1,2,3,4-tetrahydronaphthalen-2-amine (5-substituted-2-aminotetralin, 5-SAT) scaffold binds with high affinity at the 5-HT7 GPCR, and can treat symptoms of autism in mouse models; subsequently, the lead was found to have high affinity at the 5-HT1A GPCR. Herein, we report the synthesis of novel 5-SAT analogs to develop a 3-dimensional quantitative structure-affinity relationship (3D-QSAR) at the human 5-HT7 receptor for comparison with similar studies at the highly homologous 5-HT1A receptor. We report 35 new 5-SAT ligands, some with very high affinity (Ki ≤ 1 nM) and stereoselectivity at 5-HT7 + or 5-HT1A receptors, several with modest selectivity (up to 12-fold) for binding at 5-HT7, and, several ligands with high selectivity (up to 40-fold) at the 5-HT1A receptor. 3D-QSAR results indicate that steric extensions at the C(5)-position improve selectivity for the 5-HT7 over 5-HT1A receptor, while steric and hydrophobic extensions at the chiral C(2)-amino position impart 5-HT1A selectivity. In silico receptor homology modeling studies, supplemented with molecular dynamics simulations and binding free energy calculations, were used to rationalize experimentally-determined receptor selectivity and stereoselective affinity results. The data from these studies indicate that the 5-SAT chemotype, previously shown to be safe and efficacious in rodent paradigms of neurodevelopmental and neuropsychiatric disorders, is amenable to structural modification to optimize affinity at serotonin 5-HT7 vs. 5-HT1A GPCRs, as may be required for successful clinical translation.

3D interaction homology: The hydropathic interaction environments of even alanine are diverse and provide novel structural insight.

Analyses of the hydropathic environments of protein amino acid residues reveal structural information on multiple levels. The interactions made by each residue are the basis for sidechain (rotamer) conformation and ultimately for secondary, tertiary and even quaternary protein structure. By identifying and characterizing the interactions for each residue type, we are developing a basis set of environmental data that can be used to understand protein structure. This work focuses alanine and its roles. We calculated and analyzed separately backbone-to-environment and sidechain-to-environment 3D maps for over 57,000 alanines that, with respect to hydrophobic and polar interactions, show the environment around each. After binning by backbone ϕ and ψ angles, we clustered each bin with k-means based on calculated map similarities between map-map pairs. Four bins were examined in detail: one in the ß-pleat region, two in the right-hand α-helix (RHα) region and one in the left-hand α-helix region of the Ramachandran plot. All regions indicated a common map motif of hydrophobic-hydrophobic interactions along the CA-CB axis, accounting for 62% in the ß-pleat bin, about one-third in the two RHα bins and 42% in the LHα bin. Another shared motif shows no interactions along the CA-CB axis; this was uncommon (8%) in ß-pleat, but >30% elsewhere. The maps calculated for the two RHα bins are extremely similar (pairwise >0.9787), which suggests that the hydropathic interaction sets or motifs found around each residue are conserved. Altogether, these results are integral to a new paradigm for understanding protein structure and function.

Inhibiting Pneumococcal Surface Antigen A (PsaA) with Small Molecules Discovered through Virtual Screening: Steps toward Validating a Potential Target for Streptococcus pneumoniae.

The pneumococcal surface antigen A (PsaA) metal transporter protein provides manganese to bacterial cells. The X-ray crystal structures of PsaA, in both closed (Mn bound) and open (metal free) conformations, were explored with virtual screening to identify potential inhibitors of manganese transport. We pursued three strategies for inhibition: i) targeting a cavity close to the bound Mn to keep the metal in place; ii) targeting the metal-free Mn site to prevent metal uptake; and iii) targeting a potentially druggable allosteric site involving loops that translate between the conformations. Tiered assays were used to test the resulting 170 acquired hits: i) assay 1 tested the compounds' growth inhibition of the TIGR4 S. pneumoniae strain (ΔPsaA mutant control), yielding 80 compounds (MIC≤250 µm); ii) assay 2 tested if the addition of 20 µm Mn to inhibited cell cultures restored growth, yielding 21 compounds; and iii) assay 3 confirmed that the restored bacterial growth was Mn concentration dependent, as was the restoration of ΔPsaA growth, yielding 12 compounds with MICs of 125 µm or greater. It may be possible for a small molecule to inhibit PsaA, but we have not yet identified a compound with exemplary properties.

Diaminopimelic acid (DAP) analogs bearing isoxazoline moiety as selective inhibitors against meso-diaminopimelate dehydrogenase (m-Ddh) from Porphyromonas gingivalis.

Two diastereomeric analogs (1 and 2) of diaminopimelic acid (DAP) bearing an isoxazoline moiety were synthesized and evaluated for their inhibitory activities against meso-diaminopimelate dehydrogenase (m-Ddh) from the periodontal pathogen, Porphyromonas gingivalis. Compound 2 showed promising inhibitory activity against m-Ddh with an IC50 value of 14.9µM at pH 7.8. The two compounds were further tested for their antibacterial activities against a panel of periodontal pathogens, and compound 2 was shown to be selectively potent to P. gingivalis strains W83 and ATCC 33277 with minimum inhibitory concentration (MIC) values of 773µM and 1.875mM, respectively. Molecular modeling studies revealed that the inversion of chirality at the C-5 position of these compounds was the primary reason for their different biological profiles. Based on these preliminary results, we believe that compound 2 has properties consistent with it being a lead compound for developing novel pathogen selective antibiotics to treat periodontal diseases.

Ortho group activation of a bromopyrrole ester in Suzuki-Miyaura cross-coupling reactions: Application to the synthesis of new microtubule depolymerizing agents with potent cytotoxic activities.

New microtubule depolymerizing agents with potent cytotoxic activities have been prepared with a 5-cyano or 5-oximino group attached to a pyrrole core. The utilization of ortho activation of a bromopyrrole ester to facilitate successful Suzuki-Miyaura cross-coupling reactions was a key aspect of the synthetic methodology. This strategy allows for control of regiochemistry with the attachment of four completely different groups at the 2, 3, 4 and 5 positions of the pyrrole scaffold. Biological evaluations and molecular modeling studies are reported for these examples.

Biological Characterization of an Improved Pyrrole-Based Colchicine Site Agent Identified through Structure-Based Design.

A refined model of the colchicine site on tubulin was used to design an improved analog of the pyrrole parent compound, JG-03-14. The optimized compound, NT-7-16, was evaluated in biological assays that confirm that it has potent activities as a new colchicine site microtubule depolymerizer. NT-7-16 exhibits antiproliferative and cytotoxic activities against multiple cancer cell lines, with IC(50) values of 10-16 nM, and it is able to overcome drug resistance mediated by the expression of P-glycoprotein and the ßIII isotype of tubulin. NT-7-16 initiated the concentration-dependent loss of cellular microtubules and caused the formation of abnormal mitotic spindles, leading to mitotic accumulation. The direct interaction of NT-7-16 with purified tubulin was confirmed, and it was more potent than combretastatin A-4 in these assays. Binding studies verified that NT-7-16 binds to tubulin within the colchicine site. The antitumor effects of NT-7-16 were evaluated in an MDA-MB-435 xenograft model and it had excellent activity at concentrations that were not toxic. A second compound, NT-9-21, which contains dichloro moieties in place of the 3,5-dibromo substituents of NT-7-16, had a poorer fit within the colchicine site as predicted by modeling and the Hydropathic INTeractions score. Biological evaluations showed that NT-9-21 has 10-fold lower potency than NT-7-16, confirming the modeling predictions. These studies highlight the value of the refined colchicine-site model and identify a new pyrrole-based colchicine-site agent with potent in vitro activities and promising in vivo antitumor actions.

Crystal structure of carbonmonoxy sickle hemoglobin in R-state conformation.

The fundamental pathophysiology of sickle cell disease is predicated by the polymerization of deoxygenated (T-state) sickle hemoglobin (Hb S) into fibers that distort red blood cells into the characteristic sickle shape. The crystal structure of deoxygenated Hb S (DeoxyHb S) and other studies suggest that the polymer is initiated by a primary interaction between the mutation ßVal6 from one Hb S molecule, and a hydrophobic acceptor pocket formed by the residues ßAla70, ßPhe85 and ßLeu88 of an adjacent located Hb S molecule. On the contrary, oxygenated or liganded Hb S does not polymerize or incorporate in the polymer. In this paper we present the crystal structure of carbonmonoxy-ligated sickle Hb (COHb S) in the quaternary classical R-state at 1.76Å. The overall structure and the pathological donor and acceptor environments of COHb S are similar to those of the isomorphous CO-ligated R-state normal Hb (COHb A), but differ significantly from DeoxyHb S as expected. More importantly, the packing of COHb S molecules does not show the typical pathological interaction between ßVal6 and the ßAla70, ßPhe85 and ßLeu88 hydrophobic acceptor pocket observed in DeoxyHb S crystal. The structural analysis of COHb S, COHb A and DeoxyHb S provides atomic level insight into why liganded hemoglobin does not form a polymer.

Molecular basis of E. coli L-threonine aldolase catalytic inactivation at low pH.

L-Threonine aldolases (TAs), a family of enzymes belonging to the fold-type I pyridoxal 5'-phosphate (PLP) dependent enzymes, play a role in catalyzing the reversible cleavage of l-3-hydroxy-α-amino acids to glycine and the corresponding aldehydes. Threonine aldolases have great biotechnological potential for the syntheses of pharmaceutically relevant drug molecules because of their stereospecificity. The pH-dependency of their catalytic activity, affecting reaction intermediates, led us to study the effect of low-pH on Escherichia coli TA (eTA) structure. We report here a low-pH crystal structure of eTA at 2.1 Å resolution, with a non-covalently bound uncleaved l-serine substrate, and a PLP cofactor bound as an internal aldimine. This structure contrasts with other eTA structures obtained at physiological pH that show products or substrates bound as PLP-external aldimines. The non-productive binding at low-pH is due to an unusual substrate serine binding orientation in which the α-amino group and carboxylate group are in the wrong positions (relative to the active site residues) as a result of protonation of the α-amino group of the serine, as well as the active site histidines, His83 and His126. Protonation of these residues prevents the characteristic nucleophilic attack of the α-amino group of substrate serine on C4' of PLP to form the external aldimine. Our study shows that at low pH the change in charge distribution at the active site can result in substrates binding in a non-productive orientation.

Design, synthesis, and biological evaluation of substrate-competitive inhibitors of C-terminal Binding Protein (CtBP).

C-terminal Binding Protein (CtBP) is a transcriptional co-regulator that downregulates the expression of many tumor-suppressor genes. Utilizing a crystal structure of CtBP with its substrate 4-methylthio-2-oxobutyric acid (MTOB) and NAD(+) as a guide, we have designed, synthesized, and tested a series of small molecule inhibitors of CtBP. From our first round of compounds, we identified 2-(hydroxyimino)-3-phenylpropanoic acid as a potent CtBP inhibitor (IC50=0.24µM). A structure-activity relationship study of this compound further identified the 4-chloro- (IC50=0.18µM) and 3-chloro- (IC50=0.17µM) analogues as additional potent CtBP inhibitors. Evaluation of the hydroxyimine analogues in a short-term cell growth/viability assay showed that the 4-chloro- and 3-chloro-analogues are 2-fold and 4-fold more potent, respectively, than the MTOB control. A functional cellular assay using a CtBP-specific transcriptional readout revealed that the 4-chloro- and 3-chloro-hydroxyimine analogues were able to block CtBP transcriptional repression activity. This data suggests that substrate-competitive inhibition of CtBP dehydrogenase activity is a potential mechanism to reactivate tumor-suppressor gene expression as a therapeutic strategy for cancer.

3d interaction homology: The structurally known rotamers of tyrosine derive from a surprisingly limited set of information-rich hydropathic interaction environments described by maps.

Sidechain rotamer libraries are obtained through exhaustive statistical analysis of existing crystallographic structures of proteins and have been applied in multiple aspects of structural biology, for example, crystallography of relatively low-resolution structures, in homology model building and in biomolecular NMR. Little is known, however, about the driving forces that lead to the preference or suitability of one rotamer over another. Construction of 3D hydropathic interaction maps for nearly 30,000 tyrosines reveals the environment around each, in terms of hydrophobic (π-π stacking, etc.) and polar (hydrogen bonding, etc.) interactions. After partitioning the tyrosines into backbone-dependent (ϕ, ψ) bins, a map similarity metric based on the correlation coefficient was applied to each map-map pair to build matrices suitable for clustering with k-means. The first bin (-200° ≤ ϕ < -155°; -205° ≤ ψ < -160°), representing 631 tyrosines, reduced to 14 unique hydropathic environments, with most diversity arising from favorable hydrophobic interactions with many different residue partner types. Polar interactions for tyrosine include surprisingly ubiquitous hydrogen bonding with the phenolic OH and a handful of unique environments surrounding the tyrosine backbone. The memberships of all but one of the 14 environments are dominated (>50%) by a single χ(1)/χ(2) rotamer. The last environment has weak or no interactions with the tyrosine ring and its χ(1)/χ(2) rotamer is indeterminate, which is consistent with it being composed of mostly surface residues. Each tyrosine residue attempts to fulfill its hydropathic valence and thus, structural water molecules are seen in a variety of roles throughout protein structure.