Ligand-dependent dynamics and intramolecular signaling in a PDZ domain

Allosteric communication is a fundamental process that proteins use to propagate signals from one site to functionally important distal sites. Although allostery is usually associated with multimeric proteins and enzymes, “long-range” communication may be a fundamental property of proteins. In some cases, communication occurs with minimal structural change. PDZ (post-synaptic density-95/discs large/zonula occludens-1) domains are small, protein-protein binding modules that can use multiple surfaces for docking diverse molecules. Furthermore, these domains have long-range energetic couplings that link the ligand-binding site to distal regions of the structure. Here, we show that allosteric behavior in a representative member of the PDZ domain family may be directly detected using side-chain methyl dynamics measurements. The changes in side-chain dynamics parameters in the second PDZ domain from the human tyrosine phosphatase 1E (hPTP1E) were determined upon binding a peptide target. Long-range dynamic effects were detected that correspond to previously observed pair-wise energetic couplings. These results provide one of the first experimental examples for the potential role of ps-ns timescale dynamics in propagating long-range signals within a protein, and reinforce the idea that dynamic fluctuations in proteins contribute to allosteric signal transduction.     

Allosteric Communication Occurs via Networks of Tertiary and Quaternary Motions in Proteins

Allosteric proteins bind an effector molecule at one site resulting in a functional change at a second site. We hypothesize that allosteric communication in proteins relies upon networks of quaternary (collective, rigid-body) and tertiary (residue–residue contact) motions. We argue that cyclic topology of these networks is necessary for allosteric communication. An automated algorithm identifies rigid bodies from the displacement between the inactive and the active structures and constructs “quaternary networks” from these rigid bodies and the substrate and effector ligands. We then integrate quaternary networks with a coarse-grained representation of contact rearrangements to form “global communication networks” (GCNs). The GCN reveals allosteric communication among all substrate and effector sites in 15 of 18 multidomain and multimeric proteins, while tertiary and quaternary networks exhibit such communication in only 4 and 3 of these proteins, respectively. Furthermore, in 7 of the 15 proteins connected by the GCN, 50% or more of the substrate-effector paths via the GCN are “interdependent” paths that do not exist via either the tertiary or the quaternary network. Substrate-effector “pathways” typically are not linear but rather consist of polycyclic networks of rigid bodies and clusters of rearranging residue contacts. These results argue for broad applicability of allosteric communication based on structural changes and demonstrate the utility of the GCN. Global communication networks may inform a variety of experiments on allosteric proteins as well as the design of allostery into non-allosteric proteins.     

Interesting structural and dynamical behaviors exhibited by the AF-6 PDZ domain upon Bcr peptide binding

PDZ (postsynaptic density-95, disks large, zonula occludens-1) domains are small, protein-protein interaction modules that have multiple binding surfaces for the docking of diverse molecules. These domains can propagate signals from ligand-binding site to distal regions of the structure through allosteric communication. Recent works have revealed that picosecond to nanosecond time scale dynamics play a potential role in propagating long-range signals within a protein. Comparison of AF-6 PDZ domain structures in free and complex forms shows a conformation rearrangement of distal surface 2, which is far from the peptide binding groove. The relaxation dispersion experiments detected that the free AF-6 PDZ domain was sampling multiple conformations; millisecond dynamics mapped a network for allostery signal transmission throughout the AF-6 PDZ domain in the weak saturation state, and intramolecular motions were observed in distal surface 1 when the protein was saturated. These results provide evidence that the allosteric process in the AF-6 PDZ domain is not two-state; instead, the millisecond dynamic network provides a mechanism for the transmission of allosteric signals throughout a protein. Interestingly, the two distal surfaces of the AF-6 PDZ domain respond differently to peptide binding; distal surface 1 changes in millisecond dynamics, whereas distal surface 2 undergoes structural rearrangement. The significance of the different response patterns in the signaling pathway and its relevance to the function of the AF-6 PDZ domain should be studied further.     

Crystal structure and molecular dynamics of human POLDIP2, a multifaceted adaptor protein in metabolism and genome stability

Polymerase δ‐interacting protein 2 (POLDIP2, PDIP38) is a multifaceted, “moonlighting” protein, involved in binding protein partners from many different cellular processes, including mitochondrial metabolism and DNA replication and repair. How POLDIP2 interacts with many different proteins is unknown. Towards this goal, we present the crystal structure of POLDIP2 to 2.8 Å, which exhibited a compact two‐domain β‐strand‐rich globular structure, confirmed by circular dichroism and small angle X‐ray scattering approaches. POLDIP2 comprised canonical DUF525 and YccV domains, but with a conserved domain linker packed tightly, resulting in an “extended” YccV module. A central channel was observed, which we hypothesize could influence structural changes potentially mediated by redox conditions, following observation of a modified cysteine residue in the channel. Unstructured regions were rebuilt by ab initio modelling to generate a model of full‐length POLDIP2. Molecular dynamics simulations revealed a highly dynamic N‐terminal region tethered to the YccV‐domain by an extended linker, potentially facilitating interactions with distal binding partners. Models of POLDIP2 complexed with two of its partners, PrimPol and PCNA, indicated that dynamic flexibility of the POLDIP2 N‐terminus and loop regions likely mediate protein interactions.     

Contribution of the CR Domain to P-Selectin Lectin Domain Allostery by Regulating the Orientation of the EGF Domain

The allostery of P-selectin has been studied extensively with a focus on the Lec and EGF domains, whereas the contribution of the CR domain remains unclear. Here, molecular dynamics simulations (MDS) combined with homology modeling were preformed to investigate the impact of the CR domain on P-selectin allostery. The results indicated that the CR domain plays a role in the allosteric dynamics of P-selectin in two ways. First, the CR1 domain tends to stabilize the low affinity of P-selectin during the equilibration processes with the transition inhibition from the S1 to S1’ state by restraining the extension of the bent EGF orientation, or with the relaxation acceleration of the S2 state by promoting the bending of the extended EGF orientation. Second, the existence of CR domain increases intramolecular extension prior to complex separation, increasing the time available for the allosteric shift during forced dissociation with a prolonged bond duration. These findings further our understanding of the structure-function relationship of P-selectin with the enriched micro-structural bases of the CR domain.     

The order of PDZ3 and TrpCage in fusion chimeras determines their properties—a biophysical characterization

Most of the structural proteins known today are composed of domains that carry their own functions while keeping their structural properties. It is supposed that such domains, when taken out of the context of the whole protein, can retain their original structure and function to a certain extent. Information on the specific functional and structural characteristics of individual domains in a new context of artificial fusion proteins may help to reveal the rules of internal and external domain communication. Moreover, this could also help explain the mechanism of such communication and address how the mutual allosteric effect plays a role in a such multi‐domain protein system. The simple model system of the two‐domain fusion protein investigated in this work consisted of a well‐folded PDZ3 domain and an artificially designed small protein domain called Tryptophan Cage (TrpCage). Two fusion proteins with swapped domain order were designed to study their structural and functional features as well as their biophysical properties. The proteins composed of PDZ3 and TrpCage, both identical in amino acid sequence but different in composition (PDZ3‐TrpCage, TrpCage‐PDZ3), were studied using circualr dichroism (CD) spectrometry, analytical ultracentrifugation, and molecular dynamic simulations. The biophysical analysis uncovered different structural and denaturation properties of both studied proteins, revealing their different unfolding pathways and dynamics.     

Characterizing and controlling the inherent dynamics of cyclophilin-A

With the recent advances in NMR relaxation techniques, protein motions on functionally important timescales can be studied at atomic resolution. Here, we have used NMR-based relaxation experiments at several temperatures and both 600 and 900 MHz to characterize the inherent dynamics of the enzyme cyclophilin-A (CypA). We have discovered multiple chemical exchange processes within the enzyme that form a “dynamic continuum” that spans 20–30 Å comprising active site residues and residues proximal to the active site. By combining mutagenesis with these NMR relaxation techniques, a simple method of counting the dynamically sampled conformations has been developed. Surprisingly, a combination of point mutations has allowed for the specific regulation of many of the exchange processes that occur within CypA, suggesting that the dynamics of an enzyme may be engineered.     

Change in allosteric network affects binding affinities of PDZ domains: analysis through perturbation response scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation.     

Molecular and Mechanical Causes of Microtubule Catastrophe and Aging

Tubulin polymers, microtubules, can switch abruptly from the assembly to shortening. These infrequent transitions, termed “catastrophes”, affect numerous cellular processes but the underlying mechanisms are elusive. We approached this complex stochastic system using advanced coarse-grained molecular dynamics modeling of tubulin-tubulin interactions. Unlike in previous simplified models of dynamic microtubules, the catastrophes in this model arise owing to fluctuations in the composition and conformation of a growing microtubule tip, most notably in the number of protofilament curls. In our model, dynamic evolution of the stochastic microtubule tip configurations over a long timescale, known as the system’s “aging”, gives rise to the nonexponential distribution of microtubule lifetimes, consistent with experiment. We show that aging takes place in the absence of visible changes in the microtubule wall or tip, as this complex molecular-mechanical system evolves slowly and asymptotically toward the steady-state level of the catastrophe-promoting configurations. This new, to our knowledge, theoretical basis will assist detailed mechanistic investigations of the mechanisms of action of different microtubule-binding proteins and drugs, thereby enabling accurate control over the microtubule dynamics to treat various pathologies.     

LRRK2 dynamics analysis identifies allosteric control of the crosstalk between its catalytic domains

The 2 major molecular switches in biology, kinases and GTPases, are both contained in the Parkinson disease–related leucine-rich repeat kinase 2 (LRRK2). Using hydrogen–deuterium exchange mass spectrometry (HDX-MS) and molecular dynamics (MD) simulations, we generated a comprehensive dynamic allosteric portrait of the C-terminal domains of LRRK2 (LRRK2_RCKW). We identified 2 helices that shield the kinase domain and regulate LRRK2 conformation and function. One helix in COR-B (COR-B Helix) tethers the COR-B domain to the αC helix of the kinase domain and faces its activation loop, while the C-terminal helix (Ct-Helix) extends from the WD40 domain and interacts with both kinase lobes. The Ct-Helix and the N-terminus of the COR-B Helix create a “cap” that regulates the N-lobe of the kinase domain. Our analyses reveal allosteric sites for pharmacological intervention and confirm the kinase domain as the central hub for conformational control.

The Parkinson’s disease-related protein LRRK2 contains the two major molecular switches in biology; a kinase and a GTPase. This study uses hydrogen-deuterium exchange mass-spectrometry and molecular dynamics simulations to explore the conformational space of the four C-terminal domains of LRRK2, highlighting two essential regulatory helices that control LRRK2 dynamics.     

Interdomain interactions dictate the function of the Candida albicans Hsp110 protein Msi3

Heat shock proteins of 110 kDa (Hsp110s), a unique class of molecular chaperones, are essential for maintaining protein homeostasis. Hsp110s exhibit a strong chaperone activity preventing protein aggregation (the “holdase” activity) and also function as the major nucleotide-exchange factor (NEF) for Hsp70 chaperones. Hsp110s contain two functional domains: a nucleotide-binding domain (NBD) and substrate-binding domain (SBD). ATP binding is essential for Hsp110 function and results in close contacts between the NBD and SBD. However, the molecular mechanism of this ATP-induced allosteric coupling remains poorly defined. In this study, we carried out biochemical analysis on Msi3, the sole Hsp110 in Candida albicans, to dissect the unique allosteric coupling of Hsp110s using three mutations affecting the domain–domain interface. All the mutations abolished both the in vivo and in vitro functions of Msi3. While the ATP-bound state was disrupted in all mutants, only mutation of the NBD-SBDβ interfaces showed significant ATPase activity, suggesting that the full-length Hsp110s have an ATPase that is mainly suppressed by NBD-SBDβ contacts. Moreover, the high-affinity ATP-binding unexpectedly appears to require these NBD-SBD contacts. Remarkably, the “holdase” activity was largely intact for all mutants tested while NEF activity was mostly compromised, although both activities strictly depended on the ATP-bound state, indicating different requirements for these two activities. Stable peptide substrate binding to Msi3 led to dissociation of the NBD-SBD contacts and compromised interactions with Hsp70. Taken together, our data demonstrate that the exceptionally strong NBD-SBD contacts in Hsp110s dictate the unique allosteric coupling and biochemical activities.     

Human Tau Isoforms Assemble into Ribbon-like Fibrils That Display Polymorphic Structure and Stability

Fibrous aggregates of Tau protein are characteristic features of Alzheimer disease. We applied high resolution atomic force and EM microscopy to study fibrils assembled from different human Tau isoforms and domains. All fibrils reveal structural polymorphism; the “thin twisted” and “thin smooth” fibrils resemble flat ribbons (cross-section ∼10 × 15 nm) with diverse twist periodicities. “Thick fibrils” show periodicities of ∼65–70 nm and thicknesses of ∼9–18 nm such as routinely reported for “paired helical filaments” but structurally resemble heavily twisted ribbons. Therefore, thin and thick fibrils assembled from different human Tau isoforms challenge current structural models of paired helical filaments. Furthermore, all Tau fibrils reveal axial subperiodicities of ∼17–19 nm and, upon exposure to mechanical stress or hydrophobic surfaces, disassemble into uniform fragments that remain connected by thin thread-like structures (∼2 nm). This hydrophobically induced disassembly is inhibited at enhanced electrolyte concentrations, indicating that the fragments resemble structural building blocks and the fibril integrity depends largely on hydrophobic and electrostatic interactions. Because full-length Tau and repeat domain constructs assemble into fibrils of similar thickness, the “fuzzy coat” of Tau protein termini surrounding the fibril axis is nearly invisible for atomic force microscopy and EM, presumably because of its high flexibility.     

Peptide Inhibitors of Dengue-Virus Entry Target a Late-Stage Fusion Intermediate

The mechanism of membrane fusion by “class II” viral fusion proteins follows a pathway that involves large-scale domain rearrangements of the envelope glycoprotein (E) and a transition from dimers to trimers. The rearrangement is believed to proceed by an outward rotation of the E ectodomain after loss of the dimer interface, followed by a reassociation into extended trimers. The ∼55-aa-residue, membrane proximal “stem” can then zip up along domain II, bringing together the transmembrane segments of the C-terminus and the fusion loops at the tip of domain II. We find that peptides derived from the stem of dengue-virus E bind stem-less E trimer, which models a conformational intermediate. In vitro assays demonstrate that these peptides specifically block viral fusion. The peptides inhibit infectivity with potency proportional to their affinity for the conformational intermediate, even when free peptide is removed from a preincubated inoculum before infecting cells. We conclude that peptides bind virions before attachment and are carried with virions into endosomes, the compartment in which acidification initiates fusion. Binding depends on particle dynamics, as there is no inhibition of infectivity if preincubation and separation are at 4°C rather than 37°C. We propose a two-step model for the mechanism of fusion inhibition. Targeting a viral entry pathway can be an effective way to block infection. Our data, which support and extend proposed mechanisms for how the E conformational change promotes membrane fusion, suggest strategies for inhibiting flavivirus entry.     

Alternate pleckstrin homology domain orientations regulate dynamin-catalyzed membrane fission

The self-assembling GTPase dynamin catalyzes endocytic vesicle scission via membrane insertion of its pleckstrin homology (PH) domain. However, the molecular mechanisms underlying PH domain–dependent membrane fission remain obscure. Membrane-curvature–sensing and membrane-curvature–generating properties have been attributed, but it remains to be seen whether the PH domain is involved in either process independent of dynamin self-assembly. Here, using multiple fluorescence spectroscopic and microscopic techniques, we demonstrate that the isolated PH domain does not act to bend membranes but instead senses high membrane curvature through hydrophobic insertion into the membrane bilayer. Furthermore, we use a complementary set of short- and long-distance Förster resonance energy transfer approaches to distinguish PH-domain orientation from proximity at the membrane surface in full-length dynamin. We reveal, in addition to the GTP-sensitive “hydrophobic mode,” the presence of an alternate, GTP-insensitive “electrostatic mode” of PH domain–membrane interactions that retains dynamin on the membrane surface during the GTP hydrolysis cycle. Stabilization of this alternate orientation produces dramatic variations in the morphology of membrane-bound dynamin spirals, indicating that the PH domain regulates membrane fission through the control of dynamin polymer dynamics.     

Emerging computational approaches for the study of protein allostery

Allosteric regulation of protein function is key in controlling cellular processes so its underlying mechanisms are of primary concern to research in areas spanning protein engineering and drug design. However, due to the complex nature of allosteric mechanisms, a clear and predictive understanding of the relationship between protein structure and allosteric function remains elusive. Well established experimental approaches are available to offer a limited degree of characterization of mechanical properties within proteins, but the analytical capabilities of computational methods are evolving rapidly in their ability to accurately define the subtle and concerted structural dynamics that comprise allostery. This review includes a brief overview of allostery in proteins and an exploration of relevant experimental methods. An explanation of the transition from experimental toward computational methods for allostery is discussed, followed by a review of existing and emerging methods.     

A generalized allosteric mechanism for cis-regulated cyclic nucleotide binding domains

Cyclic nucleotides (cAMP and cGMP) regulate multiple intracellular processes and are thus of a great general interest for molecular and structural biologists. To study the allosteric mechanism of different cyclic nucleotide binding (CNB) domains, we compared cAMP-bound and cAMP-free structures (PKA, Epac, and two ionic channels) using a new bioinformatics method: local spatial pattern alignment. Our analysis highlights four major conserved structural motifs: 1) the phosphate binding cassette (PBC), which binds the cAMP ribose-phosphate, 2) the “hinge,” a flexible helix, which contacts the PBC, 3) the β_2,3 loop, which provides precise positioning of an invariant arginine from the PBC, and 4) a conserved structural element consisting of an N-terminal helix, an eight residue loop and the A-helix (N3A-motif). The PBC and the hinge were included in the previously reported allosteric model, whereas the definition of the β_2,3 loop and the N3A-motif as conserved elements is novel. The N3A-motif is found in all cis-regulated CNB domains, and we present a model for an allosteric mechanism in these domains. Catabolite gene activator protein (CAP) represents a trans-regulated CNB domain family: it does not contain the N3A-motif, and its long range allosteric interactions are substantially different from the cis-regulated CNB domains.     

Evaluation of energetic and dynamic coupling networks in a PDZ domain protein

A number of computational and experimental studies have identified intramolecular communication "pathways" or "networks" important for transmitting allostery. Here, we have used mutagenesis and NMR relaxation methods to investigate the scope and nature of the communication networks found in the second post-synaptic density-95/discs large/zonula occludens-1 (PDZ) domain of the human protein tyrosine phosphatase 1E protein (hPTP1E) (PDZ2). It was found that most mutations do not have a significant energetic contribution to peptide ligand binding. Three mutants that showed significant changes in binding also displayed context-dependent dynamic effects. Both a mutation at a partially exposed site (H71Y) and a buried core position (I35V) had a limited response in side-chain (2)H-based dynamics when compared to wild-type PDZ2. In contrast, a change at a second core position (I20F) that had previously been shown to be part of an energetic and dynamic network, resulted in extensive changes in side-chain dynamics. This response is reminiscent to that seen previously upon peptide ligand binding. These results shed light on the nature of the PDZ2 dynamic network and suggest that position 20 in PDZ2 acts as a "hub" that is energetically and dynamically critical for transmitting changes in dynamics throughout the PDZ domain.     

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or “cross-talk.” An example is human Pin1, an essential mitotic regulator consisting of a Trp–Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-μs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1''s numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW–PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.