Structural and Molecular Mechanisms of Cytokine-Mediated Endocrine Resistance in Human Breast Cancer Cells.

Human breast cancers that exhibit high proportions of immune cells and elevated levels of pro-inflammatory cytokines predict poor prognosis. Here, we demonstrate that treatment of human MCF-7 breast cancer cells with pro-inflammatory cytokines results in ERα-dependent activation of gene expression and proliferation, in the absence of ligand or presence of 4OH-tamoxifen (TOT). Cytokine activation of ERα and endocrine resistance is dependent on phosphorylation of ERα at S305 in the hinge domain. Phosphorylation of S305 by IKKß establishes an ERα cistrome that substantially overlaps with the estradiol (E2)-dependent ERα cistrome. Structural analyses suggest that S305-P forms a charge-linked bridge with the C-terminal F domain of ERα that enables inter-domain communication and constitutive activity from the N-terminal coactivator-binding site, revealing the structural basis of endocrine resistance. ERα therefore functions as a transcriptional effector of cytokine-induced IKKß signaling, suggesting a mechanism through which the tumor microenvironment controls tumor progression and endocrine resistance.

A posttranslational modification of the mitotic kinesin Eg5 that enhances its mechanochemical coupling and alters its mitotic function.

Numerous posttranslational modifications have been described in kinesins, but their consequences on motor mechanics are largely unknown. We investigated one of these-acetylation of lysine 146 in Eg5-by creating an acetylation mimetic lysine to glutamine substitution (K146Q). Lysine 146 is located in the α2 helix of the motor domain, where it makes an ionic bond with aspartate 91 on the neighboring α1 helix. Molecular dynamics simulations predict that disrupting this bond enhances catalytic site-neck linker coupling. We tested this using structural kinetics and single-molecule mechanics and found that the K146Q mutation increases motor performance under load and coupling of the neck linker to catalytic site. These changes convert Eg5 from a motor that dissociates from the microtubule at low load into one that is more tightly coupled and dissociation resistant-features shared by kinesin 1. These features combined with the increased propensity to stall predict that the K146Q Eg5 acetylation mimetic should act in the cell as a "brake" that slows spindle pole separation, and we have confirmed this by expressing this modified motor in mitotically active cells. Thus, our results illustrate how a posttranslational modification of a kinesin can be used to fine tune motor behavior to meet specific physiological needs.

Comparative structural dynamic analysis of GTPases.

GTPases regulate a multitude of essential cellular processes ranging from movement and division to differentiation and neuronal activity. These ubiquitous enzymes operate by hydrolyzing GTP to GDP with associated conformational changes that modulate affinity for family-specific binding partners. There are three major GTPase superfamilies: Ras-like GTPases, heterotrimeric G proteins and protein-synthesizing GTPases. Although they contain similar nucleotide-binding sites, the detailed mechanisms by which these structurally and functionally diverse superfamilies operate remain unclear. Here we compare and contrast the structural dynamic mechanisms of each superfamily using extensive molecular dynamics (MD) simulations and subsequent network analysis approaches. In particular, dissection of the cross-correlations of atomic displacements in both the GTP and GDP-bound states of Ras, transducin and elongation factor EF-Tu reveals analogous dynamic features. This includes similar dynamic communities and subdomain structures (termed lobes). For all three proteins the GTP-bound state has stronger couplings between equivalent lobes. Network analysis further identifies common and family-specific residues mediating the state-specific coupling of distal functional sites. Mutational simulations demonstrate how disrupting these couplings leads to distal dynamic effects at the nucleotide-binding site of each family. Collectively our studies extend current understanding of GTPase allosteric mechanisms and highlight previously unappreciated similarities across functionally diverse families.

Navigating the conformational landscape of G protein-coupled receptor kinases during allosteric activation.

G protein-coupled receptors (GPCRs) are essential for transferring extracellular signals into carefully choreographed intracellular responses controlling diverse aspects of cell physiology. The duration of GPCR-mediated signaling is primarily regulated via GPCR kinase (GRK)-mediated phosphorylation of activated receptors. Although many GRK structures have been reported, the mechanisms underlying GRK activation are not well-understood, in part because it is unknown how these structures map to the conformational landscape available to this enzyme family. Unlike most other AGC kinases, GRKs rely on their interaction with GPCRs for activation and not phosphorylation. Here, we used principal component analysis of available GRK and protein kinase A crystal structures to identify their dominant domain motions and to provide a framework that helps evaluate how close each GRK structure is to being a catalytically competent state. Our results indicated that disruption of an interface formed between the large lobe of the kinase domain and the regulator of G protein signaling homology domain (RHD) is highly correlated with establishment of the active conformation. By introducing point mutations in the GRK5 RHD-kinase domain interface, we show with both in silico and in vitro experiments that perturbation of this interface leads to higher phosphorylation activity. Navigation of the conformational landscape defined by this bioinformatics-based study is likely common to all GPCR-activated GRKs.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins.

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gαi. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.

Online interactive analysis of protein structure ensembles with Bio3D-web.

Bio3D-web is an online application for analyzing the sequence, structure and conformational heterogeneity of protein families. Major functionality is provided for identifying protein structure sets for analysis, their alignment and refined structure superposition, sequence and structure conservation analysis, mapping and clustering of conformations and the quantitative comparison of their predicted structural dynamics. AVAILABILITY: Bio3D-web is based on the Bio3D and Shiny R packages. All major browsers are supported and full source code is available under a GPL2 license from http://thegrantlab.org/bio3d-web CONTACT: bjgrant@umich.edu or lars.skjarven@uib.no.

Mapping the Processivity Determinants of the Kinesin-3 Motor Domain.

Kinesin superfamily members play important roles in many diverse cellular processes, including cell motility, cell division, intracellular transport, and regulation of the microtubule cytoskeleton. How the properties of the family-defining motor domain of distinct kinesins are tailored to their different cellular roles remains largely unknown. Here, we employed molecular-dynamics simulations coupled with energetic calculations to infer the family-specific interactions of kinesin-1 and kinesin-3 motor domains with microtubules in different nucleotide states. We then used experimental mutagenesis and single-molecule motility assays to further assess the predicted residue-wise determinants of distinct kinesin-microtubule binding properties. Collectively, our results identify residues in the L8, L11, and α6 regions that contribute to family-specific microtubule interactions and whose mutation affects motor-microtubule complex stability and processive motility (the ability of an individual motor to take multiple steps along its microtubule filament). In particular, substitutions of prominent kinesin-3 residues with those found in kinesin-1, namely, R167S/H171D, K266D, and R346M, were found to decrease kinesin-3 processivity 10-fold and thus approach kinesin-1 levels.

Drug hypersensitivity caused by alteration of the MHC-presented self-peptide repertoire.

Idiosyncratic adverse drug reactions are unpredictable, dose-independent and potentially life threatening; this makes them a major factor contributing to the cost and uncertainty of drug development. Clinical data suggest that many such reactions involve immune mechanisms, and genetic association studies have identified strong linkages between drug hypersensitivity reactions to several drugs and specific HLA alleles. One of the strongest such genetic associations found has been for the antiviral drug abacavir, which causes severe adverse reactions exclusively in patients expressing the HLA molecular variant B*57:01. Abacavir adverse reactions were recently shown to be driven by drug-specific activation of cytokine-producing, cytotoxic CD8(+) T cells that required HLA-B*57:01 molecules for their function; however, the mechanism by which abacavir induces this pathologic T-cell response remains unclear. Here we show that abacavir can bind within the F pocket of the peptide-binding groove of HLA-B*57:01, thereby altering its specificity. This provides an explanation for HLA-linked idiosyncratic adverse drug reactions, namely that drugs can alter the repertoire of self-peptides presented to T cells, thus causing the equivalent of an alloreactive T-cell response. Indeed, we identified specific self-peptides that are presented only in the presence of abacavir and that were recognized by T cells of hypersensitive patients. The assays that we have established can be applied to test additional compounds with suspected HLA-linked hypersensitivities in vitro. Where successful, these assays could speed up the discovery and mechanistic understanding of HLA-linked hypersensitivities, and guide the development of safer drugs.

Kinesin-5 allosteric inhibitors uncouple the dynamics of nucleotide, microtubule, and neck-linker binding sites.

Kinesin motor domains couple cycles of ATP hydrolysis to cycles of microtubule binding and conformational changes that result in directional force and movement on microtubules. The general principles of this mechanochemical coupling have been established; however, fundamental atomistic details of the underlying allosteric mechanisms remain unknown. This lack of knowledge hampers the development of new inhibitors and limits our understanding of how disease-associated mutations in distal sites can interfere with the fidelity of motor domain function. Here, we combine unbiased molecular-dynamics simulations, bioinformatics analysis, and mutational studies to elucidate the structural dynamic effects of nucleotide turnover and allosteric inhibition of the kinesin-5 motor. Multiple replica simulations of ATP-, ADP-, and inhibitor-bound states together with network analysis of correlated motions were used to create a dynamic protein structure network depicting the internal dynamic coordination of functional regions in each state. This analysis revealed the intervening residues involved in the dynamic coupling of nucleotide, microtubule, neck-linker, and inhibitor binding sites. The regions identified include the nucleotide binding switch regions, loop 5, loop 7, ?4-?5-loop 13, ?1, and ?4-?6-?7. Also evident were nucleotide- and inhibitor-dependent shifts in the dynamic coupling paths linking functional sites. In particular, inhibitor binding to the loop 5 region affected ?-sheet residues and ?1, leading to a dynamic decoupling of nucleotide, microtubule, and neck-linker binding sites. Additional analyses of point mutations, including P131 (loop 5), Q78/I79 (?1), E166 (loop 7), and K272/I273 (?7) G325/G326 (loop 13), support their predicted role in mediating the dynamic coupling of distal functional surfaces. Collectively, our results and approach, which we make freely available to the community, provide a framework for explaining how binding events and point mutations can alter dynamic couplings that are critical for kinesin motor domain function.

Integrating protein structural dynamics and evolutionary analysis with Bio3D.

BACKGROUND: Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution. RESULTS: Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case. CONCLUSIONS: The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/ .

Electrostatically biased binding of kinesin to microtubules.

The minimum motor domain of kinesin-1 is a single head. Recent evidence suggests that such minimal motor domains generate force by a biased binding mechanism, in which they preferentially select binding sites on the microtubule that lie ahead in the progress direction of the motor. A specific molecular mechanism for biased binding has, however, so far been lacking. Here we use atomistic Brownian dynamics simulations combined with experimental mutagenesis to show that incoming kinesin heads undergo electrostatically guided diffusion-to-capture by microtubules, and that this produces directionally biased binding. Kinesin-1 heads are initially rotated by the electrostatic field so that their tubulin-binding sites face inwards, and then steered towards a plus-endwards binding site. In tethered kinesin dimers, this bias is amplified. A 3-residue sequence (RAK) in kinesin helix alpha-6 is predicted to be important for electrostatic guidance. Real-world mutagenesis of this sequence powerfully influences kinesin-driven microtubule sliding, with one mutant producing a 5-fold acceleration over wild type. We conclude that electrostatic interactions play an important role in the kinesin stepping mechanism, by biasing the diffusional association of kinesin with microtubules.

Mapping the structural and dynamical features of kinesin motor domains.

Kinesin motor proteins drive intracellular transport by coupling ATP hydrolysis to conformational changes that mediate directed movement along microtubules. Characterizing these distinct conformations and their interconversion mechanism is essential to determining an atomic-level model of kinesin action. Here we report a comprehensive principal component analysis of 114 experimental structures along with the results of conventional and accelerated molecular dynamics simulations that together map the structural dynamics of the kinesin motor domain. All experimental structures were found to reside in one of three distinct conformational clusters (ATP-like, ADP-like and Eg5 inhibitor-bound). These groups differ in the orientation of key functional elements, most notably the microtubule binding α4-α5, loop8 subdomain and α2b-ß4-ß6-ß7 motor domain tip. Group membership was found not to correlate with the nature of the bound nucleotide in a given structure. However, groupings were coincident with distinct neck-linker orientations. Accelerated molecular dynamics simulations of ATP, ADP and nucleotide free Eg5 indicate that all three nucleotide states could sample the major crystallographically observed conformations. Differences in the dynamic coupling of distal sites were also evident. In multiple ATP bound simulations, the neck-linker, loop8 and the α4-α5 subdomain display correlated motions that are absent in ADP bound simulations. Further dissection of these couplings provides evidence for a network of dynamic communication between the active site, microtubule-binding interface and neck-linker via loop7 and loop13. Additional simulations indicate that the mutations G325A and G326A in loop13 reduce the flexibility of these regions and disrupt their couplings. Our combined results indicate that the reported ATP and ADP-like conformations of kinesin are intrinsically accessible regardless of nucleotide state and support a model where neck-linker docking leads to a tighter coupling of the microtubule and nucleotide binding regions. Furthermore, simulations highlight sites critical for large-scale conformational changes and the allosteric coupling between distal functional sites.

Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins.

Heterotrimeric G proteins are conformational switches that turn on intracellular signaling cascades in response to the activation of G-protein-coupled receptors. Receptor activation by extracellular stimuli promotes a cycle of GTP binding and hydrolysis on the G protein α-subunit (Gα). Important conformational transitions occurring during this cycle have been characterized from extensive crystallographic studies of Gα. However, the link between the observed conformations and the mechanisms involved in G-protein activation and effector interaction remain unclear. Here we describe a comprehensive principal component analysis of available Gα crystallographic structures supplemented with extensive unbiased conventional and accelerated molecular dynamics simulations that together characterize the response of Gα to GTP binding and hydrolysis. Our studies reveal details of activating conformational changes as well as the intrinsic flexibility of the α-helical domain that includes a large-scale 60° domain opening under nucleotide-free conditions. This result is consistent with the recently reported open crystal structure of Gs, the stimulatory G protein for adenylyl cyclase, in complex with the α2 adrenergic receptor. Sets of unique interactions potentially important for the conformational transition are also identified. Moreover simulations reveal nucleotide-dependent dynamical couplings of distal regions and residues potentially important for the allosteric link between functional sites.

LeuT conformational sampling utilizing accelerated molecular dynamics and principal component analysis.

Monoamine transporters (MATs) function by coupling ion gradients to the transport of dopamine, norepinephrine, or serotonin. Despite their importance in regulating neurotransmission, the exact conformational mechanism by which MATs function remains elusive. To this end, we have performed seven 250 ns accelerated molecular dynamics simulations of the leucine transporter, a model for neurotransmitter MATs. By varying the presence of binding-pocket leucine substrate and sodium ions, we have sampled plausible conformational states representative of the substrate transport cycle. The resulting trajectories were analyzed using principal component analysis of transmembrane helices 1b and 6a. This analysis revealed seven unique structures: two of the obtained conformations are similar to the currently published crystallographic structures, one conformation is similar to a proposed open inward structure, and four conformations represent novel structures of potential importance to the transport cycle. Further analysis reveals that the presence of binding-pocket sodium ions is necessary to stabilize the locked-occluded and open-inward conformations.

Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics.

Periplasmic binding proteins (PBPs) are a large family of molecular transporters that play a key role in nutrient uptake and chemotaxis in Gram-negative bacteria. All PBPs have characteristic two-domain architecture with a central interdomain ligand-binding cleft. Upon binding to their respective ligands, PBPs undergo a large conformational change that effectively closes the binding cleft. This conformational change is traditionally viewed as a ligand induced-fit process; however, the intrinsic dynamics of the protein may also be crucial for ligand recognition. Recent NMR paramagnetic relaxation enhancement (PRE) experiments have shown that the maltose binding protein (MBP) - a prototypical member of the PBP superfamily - exists in a rapidly exchanging (ns to µs regime) mixture comprising an open state (approx 95%), and a minor partially closed state (approx 5%). Here we describe accelerated MD simulations that provide a detailed picture of the transition between the open and partially closed states, and confirm the existence of a dynamical equilibrium between these two states in apo MBP. We find that a flexible part of the protein called the balancing interface motif (residues 175-184) is displaced during the transformation. Continuum electrostatic calculations indicate that the repacking of non-polar residues near the hinge region plays an important role in driving the conformational change. Oscillations between open and partially closed states create variations in the shape and size of the binding site. The study provides a detailed description of the conformational space available to ligand-free MBP, and has implications for understanding ligand recognition and allostery in related proteins.

Large-scale conformational changes of Trypanosoma cruzi proline racemase predicted by accelerated molecular dynamics simulation.

Chagas' disease, caused by the protozoan parasite Trypanosoma cruzi (T. cruzi), is a life-threatening illness affecting 11-18 million people. Currently available treatments are limited, with unacceptable efficacy and safety profiles. Recent studies have revealed an essential T. cruzi proline racemase enzyme (TcPR) as an attractive candidate for improved chemotherapeutic intervention. Conformational changes associated with substrate binding to TcPR are believed to expose critical residues that elicit a host mitogenic B-cell response, a process contributing to parasite persistence and immune system evasion. Characterization of the conformational states of TcPR requires access to long-time-scale motions that are currently inaccessible by standard molecular dynamics simulations. Here we describe advanced accelerated molecular dynamics that extend the effective simulation time and capture large-scale motions of functional relevance. Conservation and fragment mapping analyses identified potential conformational epitopes located in the vicinity of newly identified transient binding pockets. The newly identified open TcPR conformations revealed by this study along with knowledge of the closed to open interconversion mechanism advances our understanding of TcPR function. The results and the strategy adopted in this work constitute an important step toward the rationalization of the molecular basis behind the mitogenic B-cell response of TcPR and provide new insights for future structure-based drug discovery.

Induced fit or conformational selection? The role of the semi-closed state in the maltose binding protein.

A full characterization of the thermodynamic forces underlying ligand-associated conformational changes in proteins is essential for understanding and manipulating diverse biological processes, including transport, signaling, and enzymatic activity. Recent experiments on the maltose binding protein (MBP) have provided valuable data about the different conformational states implicated in the ligand recognition process; however, a complete picture of the accessible pathways and the associated changes in free energy remains elusive. Here we describe results from advanced accelerated molecular dynamics (aMD) simulations, coupled with adaptively biased force (ABF) and thermodynamic integration (TI) free energy methods. The combination of approaches allows us to track the ligand recognition process on the microsecond time scale and provides a detailed characterization of the protein's dynamic and the relative energy of stable states. We find that an induced-fit (IF) mechanism is most likely and that a mechanism involving both a conformational selection (CS) step and an IF step is also possible. The complete recognition process is best viewed as a "Pac Man" type action where the ligand is initially localized to one domain and naturally occurring hinge-bending vibrations in the protein are able to assist the recognition process by increasing the chances of a favorable encounter with side chains on the other domain, leading to a population shift. This interpretation is consistent with experiments and provides new insight into the complex recognition mechanism. The methods employed here are able to describe IF and CS effects and provide formally rigorous means of computing free energy changes. As such, they are superior to conventional MD and flexible docking alone and hold great promise for future development and applications to drug discovery.

The distinct conformational dynamics of K-Ras and H-Ras A59G.

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species.

The Bio3D packages for structural bioinformatics.

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.

Conformational selection in G-proteins: lessons from Ras and Rho.

The induced fit model has traditionally been invoked to describe the activating conformational change of the monomeric G-proteins, such as Ras and Rho. With this scheme, the presence or absence of the γ-phosphate of GTP leads to an instantaneous switch in conformation. Here we describe atomistic molecular simulations that demonstrate that both Ras and Rho superfamily members harbor an intrinsic susceptibility to sample multiple conformational states in the absence of nucleotide ligand. By comparing the distribution of conformers in the presence and absence of nucleotide, we show that conformational selection is the dominant mechanism by which Ras and Rho undergo nucleotide-dependent conformational changes. Furthermore, the pattern of correlated motions revealed by these simulations predicts a preserved allosteric coupling of the nucleotide-binding site with the membrane interacting C-terminus in both Rho and Ras.