Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery.

Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence-structure-function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.

Incorporation of D2O-Induced Fluorine Chemical Shift Perturbations into Ensemble-Structure Characterization of the ERalpha Disordered Region.

Structural characterization of intrinsically disordered proteins (IDPs) requires a concerted effort between experiments and computations by accounting for their conformational heterogeneity. Given the diversity of experimental tools providing local and global structural information, constructing an experimental restraint-satisfying structural ensemble remains challenging. Here, we use the disordered N-terminal domain (NTD) of the estrogen receptor alpha (ERalpha) as a model system to combine existing small-angle X-ray scattering (SAXS) and hydroxyl radical protein footprinting (HRPF) data and newly acquired solvent accessibility data via D2O-induced fluorine chemical shifting (DFCS) measurements. A new set of DFCS data for the solvent exposure of a set of 12 amino acid positions were added to complement previously acquired HRPF measurements for the solvent exposure of the other 16 nonoverlapping amino acids, thereby improving the NTD ensemble characterization considerably. We also found that while choosing an initial ensemble of structures generated from a different atomic-level force field or sampling/modeling method can lead to distinct contact maps even when the same sets of experimental measurements were used for ensemble-fitting, comparative analyses from these initial ensembles reveal commonly recurring structural features in their ensemble-averaged contact map. Specifically, nonlocal or long-range transient interactions were found consistently between the N-terminal segments and the central region, sufficient to mediate the conformational ensemble and regulate how the NTD interacts with its coactivator proteins.

Fusion-associated carcinomas of the breast: Diagnostic, prognostic, and therapeutic significance.

Recurrent gene fusions comprise a class of viable genetic targets in solid tumors that have culminated several recent breakthrough cancer therapies. Their role in breast cancer, however, remains largely underappreciated due to the complexity of genomic rearrangements in breast malignancy. Just recently, we and others have identified several recurrent gene fusions in breast cancer with important clinical and biological implications. Examples of the most significant recurrent gene fusions to date include (1) ESR1::CCDC170 gene fusions in luminal B and endocrine-resistant breast cancer that exert oncogenic function via modulating the HER2/HER3/SRC Proto-Oncogene (SRC) complex, (2) ESR1 exon 6 fusions in metastatic disease that drive estrogen-independent estrogen-receptor transcriptional activity, (3) BCL2L14::ETV6 fusions in a more aggressive form of the triple-negative subtype that prime epithelial-mesenchymal transition and endow paclitaxel resistance, (4) the ETV6::NTRK3 fusion in secretory breast carcinoma that constitutively activates NTRK3 kinase, (5) the oncogenic MYB-NFIB fusion as a genetic driver underpinning adenoid cystic carcinomas of the breast that activates MYB Proto-Oncogene (MYB) pathway, and (6) the NOTCH/microtubule-associated serine-threonine (MAST) kinase gene fusions that activate NOTCH and MAST signaling. Importantly, these fusions are enriched in more aggressive and lethal breast cancer presentations and appear to confer therapeutic resistance. Thus, these gene fusions could be utilized as genetic biomarkers to identify patients who require more intensive treatment and surveillance. In addition, kinase fusions are currently being evaluated in breast cancer clinical trials and ongoing mechanistic investigation is exposing therapeutic vulnerabilities in patients with fusion-positive disease.

Integrative structural modeling of a multidomain polo-like kinase.

Polo-like kinase 1 (PLK1) is a key regulator and coordinator for mitotic signaling that contains two major functional units of a kinase domain (KD) and a polo-box domain (PBD). While individual domain structures of the KD and the PBD are known, how they interact and assemble into a functional complex remains an open question. The structural model from the KD-PBD-Map205PBM heterotrimeric crystal structure of zebrafish PLK1 represents a major step in understanding the KD and the PBD interactions. However, how these two domains interact when connected by a linker in the full length PLK1 needs further investigation. By integrating different sources of structural data from small-angle X-ray scattering, hydroxyl radical protein footprinting, and computational sampling, here we report an overall architecture for PLK1 multidomain assembly between the KD and the PBD. Our model revealed that the KD uses its C-lobe to interact with the PBD via the site near the phosphopeptide binding site in its auto-inhibitory state in solution. Disruption of this auto-inhibition via site-directed mutagenesis at the KD-PBD interface increases its kinase activity, supporting the functional role of KD-PBD interactions predicted for regulating the PLK1 kinase function. Our results indicate that the full length human PLK1 takes dynamic structures with a variety of domain-domain interfaces in solution.

Glycine substitution in SH3-SH2 connector of Hck tyrosine kinase causes population shift from assembled to disassembled state.

A combination of small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations based on a coarse grained model is used to examine the effect of glycine substitutions in the short connector between the SH3 and SH2 domains of Hck, a member of the Src-family kinases. It has been shown previously that the activity of cSrc kinase is upregulated by substitution of 3 residues by glycine in the SH3-SH2 connector. Here, analysis of SAXS data indicates that the population of Hck in the disassembled state increases from 25% in the wild type kinase to 76% in the glycine mutant. This is consistent with the results of free energy perturbation calculations showing that the mutation in the connector shifts the equilibrium from the assembled to the disassembled state. This study supports the notion that the SH3-SH2 connector helps to regulate the activity of tyrosine kinases by shifting the population of the active state of the multidomain protein independent of C-terminal phosphorylation.

Protein Footprinting: Auxiliary Engine to Power the Structural Biology Revolution.

Structural biology is entering an exciting time where many new high-resolution structures of large complexes and membrane proteins are determined regularly. These advances have been driven by over fifteen years of technology advancements, first in macromolecular crystallography, and recently in Cryo-electron microscopy. These structures are allowing detailed questions about functional mechanisms of the structures, and the biology enabled by these structures, to be addressed for the first time. At the same time, mass spectrometry technologies for protein structure analysis, "footprinting" studies, have improved their sensitivity and resolution dramatically and can provide detailed sub-peptide and residue level information for validating structures and interactions or understanding the dynamics of structures in the context of ligand binding or assembly. In this perspective, we review the use of protein footprinting to extend our understanding of macromolecular systems, particularly for systems challenging for analysis by other techniques, such as intrinsically disordered proteins, amyloidogenic proteins, and other proteins/complexes so far recalcitrant to existing methods. We also illustrate how the availability of high-resolution structural information can be a foundation for a suite of hybrid approaches to divine structure-function relationships beyond what individual techniques can deliver.

A Metastable Contact and Structural Disorder in the Estrogen Receptor Transactivation Domain.

The N-terminal transactivation domain (NTD) of estrogen receptor alpha, a well-known member of the family of intrinsically disordered proteins, mediates the receptor's transactivation function. However, an accurate molecular dissection of NTD's structure-function relationships remains elusive. Here, we show that the NTD adopts a mostly disordered, unexpectedly compact conformation that undergoes structural expansion on chemical denaturation. By combining small-angle X-ray scattering, hydroxyl radical protein footprinting, and computational modeling, we derive the ensemble-structures of the NTD and determine its ensemble-contact map revealing metastable long-range contacts, e.g., between residues I33 and S118. We show that mutation at S118, a known phosphorylation site, promotes conformational changes and increases coactivator binding. We further demonstrate via fluorine-19 (19F) nuclear magnetic resonance that mutations near I33 alter 19F chemical shifts at S118, confirming the proposed I33-S118 contact in the ensemble of structural disorder. These findings extend our understanding of how specific contact metastability mediates critical functions of disordered proteins.

Alleviated Inhibition of Single Enzyme in Confined and Crowded Environment.

Most proteins perform functions in intracellular milieu. The crowding, compartmentalized cytosol environment affects the protein structure, folding, conformational stability, substrate diffusion, and substrate-enzyme binding. Moreover, enzymes are available at single or very low copy numbers in a cell, and thus the conformation fluctuations of a single enzyme in a crowding environment could also greatly influence its kinetics. However, the crowding effect is poorly understood in the kinetical aspect of enzymatic reactions. In the present study, individual horseradish peroxidase (HRP) is encapsulated in a liposome containing crowding reagents as mimics of viscous cytosol. The confined crowding environment possesses a profound influence on both the catalytic activity and the product inhibition of enzymes. By analyzing the correlation between product generation and product inhibition, we find that the allosteric noncompetitive inhibition of HRP is alleviated in the crowded and confined milieu. Small-angle X-ray scattering experiments provide straightforward proofs of structural changes of enzymes in crowding environments, which are responsible for the reduced enzyme activity and increased enzyme-substrate affinity. We expect that this work may deepen the understanding of correlations between enzymatic conformations and activity performance in real cellular environments.

Reactive cysteine residues in the oxidative dimerization and Cu2+ induced aggregation of human γD-crystallin: Implications for age-related cataract.

Cysteine (Cys) residues are major causes of crystallin disulfide formation and aggregation in aging and cataractous human lenses. We recently found that disulfide linkages are highly and partly conserved in ß- and γ-crystallins, respectively, in human age-related nuclear cataract and glutathione depleted LEGSKO mouse lenses, and could be mimicked by in vitro oxidation. Here we determined which Cys residues are involved in disulfide-mediated crosslinking of recombinant human γD-crystallin (hγD). In vitro diamide oxidation revealed dimer formation by SDS-PAGE and LC-MS analysis with Cys 111-111 and C111-C19 as intermolecular disulfides and Cys 111-109 as intramolecular sites. Mutation of Cys111 to alanine completely abolished dimerization. Addition of αB-crystallin was unable to protect Cys 111 from dimerization. However, Cu2+-induced hγD-crystallin aggregation was suppressed up to 50% and 80% by mutants C109A and C111A, respectively, as well as by total glutathionylation. In contrast to our recently published results using ICAT-labeling method, manual mining of the same database confirmed the specific involvement of Cys111 in disulfides with no free Cys111 detectable in γD-crystallin from old and cataractous human lenses. Surface accessibility studies show that Cys111 in hγD is the most exposed Cys residue (29%), explaining thereby its high propensity toward oxidation and polymerization in the aging lens.

Multidomain architecture of estrogen receptor reveals interfacial cross-talk between its DNA-binding and ligand-binding domains.

Human estrogen receptor alpha (hERα) is a hormone-responsive nuclear receptor (NR) involved in cell growth and survival that contains both a DNA-binding domain (DBD) and a ligand-binding domain (LBD). Functionally relevant inter-domain interactions between the DBD and LBD have been observed in several other NRs, but for hERα, the detailed structural architecture of the complex is unknown. By utilizing integrated complementary techniques of small-angle X-ray scattering, hydroxyl radical protein footprinting and computational modeling, here we report an asymmetric L-shaped "boot" structure of the multidomain hERα and identify the specific sites on each domain at the domain interface involved in DBD-LBD interactions. We demonstrate the functional role of the proposed DBD-LBD domain interface through site-specific mutagenesis altering the hERα interfacial structure and allosteric signaling. The L-shaped structure of hERα is a distinctive DBD-LBD organization of NR complexes and more importantly, reveals a signaling mechanism mediated by inter-domain crosstalk that regulates this receptor's allosteric function.

Cyclin-dependent kinase 7 (CDK7)-mediated phosphorylation of the CDK9 activation loop promotes P-TEFb assembly with Tat and proviral HIV reactivation.

The HIV trans-activator Tat recruits the host transcription elongation factor P-TEFb to stimulate proviral transcription. Phosphorylation of Thr-186 on the activation loop (T-loop) of cyclin-dependent kinase 9 (CDK9) is essential for its kinase activity and assembly of CDK9 and cyclin T1 (CycT1) to form functional P-TEFb. Phosphorylation of a second highly conserved T-loop site, Ser-175, alters the competitive binding of Tat and the host recruitment factor bromodomain containing 4 (BRD4) to P-TEFb. Here, we investigated the intracellular mechanisms that regulate these key phosphorylation events required for HIV transcription. Molecular dynamics simulations revealed that the CDK9/CycT1 interface is stabilized by intramolecular hydrogen bonding of pThr-186 by an arginine triad and Glu-96 of CycT1. Arginine triad substitutions that disrupted CDK9/CycT1 assembly accumulated Thr-186-dephosphorylated CDK9 associated with the cytoplasmic Hsp90/Cdc37 chaperone. The Hsp90/Cdc37/CDK9 complex was also present in resting T cells, which lack CycT1. Hsp90 inhibition in primary T cells blocked P-TEFb assembly, disrupted Thr-186 phosphorylation, and suppressed proviral reactivation. The selective CDK7 inhibitor THZ1 blocked CDK9 phosphorylation at Ser-175, and in vitro kinase assays confirmed that CDK7 activity is principally responsible for Ser-175 phosphorylation. Mutation of Ser-175 to Lys had no effect on CDK9 kinase activity or P-TEFb assembly but strongly suppressed both HIV expression and BRD4 binding. We conclude that the transfer of CDK9 from the Hsp90/Cdc37 complex induced by Thr-186 phosphorylation is a key step in P-TEFb biogenesis. Furthermore, we demonstrate that CDK7-mediated Ser-175 phosphorylation is a downstream nuclear event essential for facilitating CDK9 T-loop interactions with Tat.

Charge Neutralization Drives the Shape Reconfiguration of DNA Nanotubes.

Reconfiguration of membrane protein channels for gated transport is highly regulated under physiological conditions. However, a mechanistic understanding of such channels remains challenging owing to the difficulty in probing subtle gating-associated structural changes. Herein, we show that charge neutralization can drive the shape reconfiguration of a biomimetic 6-helix bundle DNA nanotube (6HB). Specifically, 6HB adopts a compact state when its charge is neutralized by Mg2+ ; whereas Na+ switches it to the expanded state, as revealed by MD simulations, small-angle X-ray scattering (SAXS), and FRET characterization. Furthermore, partial neutralization of the DNA backbone charges by chemical modification renders 6HB compact and insensitive to ions, suggesting an interplay between electrostatic and hydrophobic forces in the channels. This system provides a platform for understanding the structure-function relationship of biological channels and designing rules for the shape control of DNA nanostructures in biomedical applications.

A Practical Guide to iSPOT Modeling: An Integrative Structural Biology Platform.

Integrative structure modeling is an emerging method for structural determination of protein-protein complexes that are challenging for conventional structural techniques. Here, we provide a practical protocol for implementing our integrated iSPOT platform by integrating three different biophysical techniques: small-angle X-ray scattering (SAXS), hydroxyl radical footprinting, and computational docking simulations. Specifically, individual techniques are described from experimental and/or computational perspectives, and complementary structural information from these different techniques are integrated for accurate characterization of the structures of large protein-protein complexes.

Structural basis of ligand interaction with atypical chemokine receptor 3.

Chemokines drive cell migration through their interactions with seven-transmembrane (7TM) chemokine receptors on cell surfaces. The atypical chemokine receptor 3 (ACKR3) binds chemokines CXCL11 and CXCL12 and signals exclusively through ß-arrestin-mediated pathways, without activating canonical G-protein signalling. This receptor is upregulated in numerous cancers making it a potential drug target. Here we collected over 100 distinct structural probes from radiolytic footprinting, disulfide trapping, and mutagenesis to map the structures of ACKR3:CXCL12 and ACKR3:small-molecule complexes, including dynamic regions that proved unresolvable by X-ray crystallography in homologous receptors. The data are integrated with molecular modelling to produce complete and cohesive experimentally driven models that confirm and expand on the existing knowledge of the architecture of receptor:chemokine and receptor:small-molecule complexes. Additionally, we detected and characterized ligand-induced conformational changes in the transmembrane and intracellular regions of ACKR3 that elucidate fundamental structural elements of agonism in this atypical receptor.

Thermodynamics of Hydrophobic Amino Acids in Solution: A Combined Experimental-Computational Study.

We present a joint experimental-computational study to quantitatively describe the thermodynamics of hydrophobic leucine amino acids in aqueous solution. X-ray scattering data were acquired at a series of solute and salt concentrations to effectively measure interleucine interactions, indicating that a major scattering peak is observed consistently at q = 0.83 Å-1. Atomistic molecular dynamics simulations were then performed and compared with the scattering data, achieving high consistency at both small and wider scattering angles (q = 0-1.5 Å-1). This experimental-computational consistence enables a first glimpse of the leucine-leucine interacting landscape, where two leucine molecules are aligned mostly in a parallel fashion, as opposed to antiparallel, but also allows us to derive effective leucine-leucine interactions in solution. Collectively, this combined approach of employing experimental scattering and molecular simulation enables quantitative characterization of effective intermolecular interactions of hydrophobic amino acids, critical for protein function and dynamics such as protein folding.

Theoretical modeling of multiprotein complexes by iSPOT: Integration of small-angle X-ray scattering, hydroxyl radical footprinting, and computational docking.

Structural determination of protein-protein complexes such as multidomain nuclear receptors has been challenging for high-resolution structural techniques. Here, we present a combined use of multiple biophysical methods, termed iSPOT, an integration of shape information from small-angle X-ray scattering (SAXS), protection factors probed by hydroxyl radical footprinting, and a large series of computationally docked conformations from rigid-body or molecular dynamics (MD) simulations. Specifically tested on two model systems, the power of iSPOT is demonstrated to accurately predict the structures of a large protein-protein complex (TGFß-FKBP12) and a multidomain nuclear receptor homodimer (HNF-4α), based on the structures of individual components of the complexes. Although neither SAXS nor footprinting alone can yield an unambiguous picture for each complex, the combination of both, seamlessly integrated in iSPOT, narrows down the best-fit structures that are about 3.2Å and 4.2Å in RMSD from their corresponding crystal structures, respectively. Furthermore, this proof-of-principle study based on the data synthetically derived from available crystal structures shows that the iSPOT-using either rigid-body or MD-based flexible docking-is capable of overcoming the shortcomings of standalone computational methods, especially for HNF-4α. By taking advantage of the integration of SAXS-based shape information and footprinting-based protection/accessibility as well as computational docking, this iSPOT platform is set to be a powerful approach towards accurate integrated modeling of many challenging multiprotein complexes.

Accurate optimization of amino acid form factors for computing small-angle X-ray scattering intensity of atomistic protein structures.

Structure modelling via small-angle X-ray scattering (SAXS) data generally requires intensive computations of scattering intensity from any given biomolecular structure, where the accurate evaluation of SAXS profiles using coarse-grained (CG) methods is vital to improve computational efficiency. To date, most CG SAXS computing methods have been based on a single-bead-per-residue approximation but have neglected structural correlations between amino acids. To improve the accuracy of scattering calculations, accurate CG form factors of amino acids are now derived using a rigorous optimization strategy, termed electron-density matching (EDM), to best fit electron-density distributions of protein structures. This EDM method is compared with and tested against other CG SAXS computing methods, and the resulting CG SAXS profiles from EDM agree better with all-atom theoretical SAXS data. By including the protein hydration shell represented by explicit CG water molecules and the correction of protein excluded volume, the developed CG form factors also reproduce the selected experimental SAXS profiles with very small deviations. Taken together, these EDM-derived CG form factors present an accurate and efficient computational approach for SAXS computing, especially when higher molecular details (represented by the q range of the SAXS data) become necessary for effective structure modelling.

Quantitative protein topography analysis and high-resolution structure prediction using hydroxyl radical labeling and tandem-ion mass spectrometry (MS).

Hydroxyl radical footprinting based MS for protein structure assessment has the goal of understanding ligand induced conformational changes and macromolecular interactions, for example, protein tertiary and quaternary structure, but the structural resolution provided by typical peptide-level quantification is limiting. In this work, we present experimental strategies using tandem-MS fragmentation to increase the spatial resolution of the technique to the single residue level to provide a high precision tool for molecular biophysics research. Overall, in this study we demonstrated an eightfold increase in structural resolution compared with peptide level assessments. In addition, to provide a quantitative analysis of residue based solvent accessibility and protein topography as a basis for high-resolution structure prediction; we illustrate strategies of data transformation using the relative reactivity of side chains as a normalization strategy and predict side-chain surface area from the footprinting data. We tested the methods by examination of Ca(+2)-calmodulin showing highly significant correlations between surface area and side-chain contact predictions for individual side chains and the crystal structure. Tandem ion based hydroxyl radical footprinting-MS provides quantitative high-resolution protein topology information in solution that can fill existing gaps in structure determination for large proteins and macromolecular complexes.

Quantitative mapping of protein structure by hydroxyl radical footprinting-mediated structural mass spectrometry: a protection factor analysis.

Measurements from hydroxyl radical footprinting (HRF) provide rich information about the solvent accessibility of amino acid side chains of a protein. Traditional HRF data analyses focus on comparing the difference in the modification/footprinting rate of a specific site to infer structural changes across two protein states, e.g., between a free and ligand-bound state. However, the rate information itself is not fully used for the purpose of comparing different protein sites within a protein on an absolute scale. To provide such a cross-site comparison, we present a new, to our knowledge, data analysis algorithm to convert the measured footprinting rate constant to a protection factor (PF) by taking into account the known intrinsic reactivity of amino acid side chain. To examine the extent to which PFs can be used for structural interpretation, this PF analysis is applied to three model systems where radiolytic footprinting data are reported in the literature. By visualizing structures colored with the PF values for individual peptides, a rational view of the structural features of various protein sites regarding their solvent accessibility is revealed, where high-PF regions are buried and low-PF regions are more exposed to the solvent. Furthermore, a detailed analysis correlating solvent accessibility and local structural contacts for gelsolin shows a statistically significant agreement between PF values and various structure measures, demonstrating that the PFs derived from this PF analysis readily explain fundamental HRF rate measurements. We also tested this PF analysis on alternative, chemical-based HRF data, showing improved correlations of structural properties of a model protein barstar compared to examining HRF rate data alone. Together, this PF analysis not only permits a novel, to our knowledge, approach of mapping protein structures by using footprinting data, but also elevates the use of HRF measurements from a qualitative, cross-state comparison to a quantitative, cross-site assessment of protein structures in the context of individual conformational states of interest.