Visualization of allostery in P-selectin lectin domain using MD simulations

Allostery of P-selectin lectin (Lec) domain followed by an epithelial growth factor (EGF)-like domain is essential for its biological functionality, but the underlying pathways have not been well understood. Here the molecular dynamics simulations were performed on the crystallized structures to visualize the dynamic conformational change for state 1 (S1) or state 2 (S2) Lec domain with respective bent (B) or extended (E) EGF orientation. Simulations illustrated that both S1 and S2 conformations were unable to switch from one to another directly. Instead, a novel S1' conformation was observed from S1 when crystallized B-S1 or reconstructed "E-S1" structure was employed, which was superposed well with that of equilibrated S1 Lec domain alone. It was also indicated that the corresponding allosteric pathway from S1 to S1' conformation started with the separation between residues Q30 and K67 and terminated with the release of residue N87 from residue C109. These results provided an insight into understanding the structural transition and the structure-function relationship of P-selectin allostery.     

P-selectin must extend a sufficient length from the plasma membrane to mediate rolling of neutrophils

Under physiological shear stress, neutrophils roll on P-selectin on activated endothelial cells or platelets through interactions with P- selectin glycoprotein ligand-1 (PSGL-1). Both P-selectin and PSGL-1 are extended molecules. Human P-selectin contains an NH2-terminal lectin domain, an EGF domain, nine consensus repeats (CRs), a transmembrane domain, and a cytoplasmic tail. To determine whether the length of P- selectin affected its interactions with PSGL-1, we examined the adhesion of neutrophils to CHO cells expressing membrane-anchored P- selectin constructs in which various numbers of CRs were deleted. Under static conditions, neutrophils attached equivalently to wild-type P- selectin and to constructs containing from 2-6 CRs. Under shear stress, neutrophils attached equivalently to wild-type and 6 CR P-selectin and nearly as well to 5 CR P-selectin. However, fewer neutrophils attached to the 4 CR construct, and those that did attach rolled faster and were more readily detached by increasing shear stress. Flowing neutrophils failed to attach to the 3 CR and 2 CR constructs. Neutrophils attached and rolled more efficiently on 4 CR P-selectin expressed on glycosylation-defective Lec8 CHO cells, which have less glycocalyx. We conclude that P-selectin must project its lectin domain well above the membrane to mediate optimal attachment of neutrophils under shear forces. The length of P-selectin may: (a) facilitate interactions with PSGL-1 on flowing neutrophils, and (b) increase the intermembrane distance where specific bonds form, minimizing contacts between the glycocalyces that result in cell-cell repulsion.     

Contrasting roles of dynamics in protein allostery: NMR and structural studies of CheY and the third PDZ domain from PSD-95

Allosteric regulation is a ubiquitous phenomenon exploited in biological processes to control cells in a myriad of ways. It is also of emerging interest in the design of functional proteins and therapeutics. Even though allostery was proposed over 50 years ago and has been studied intensively from a structural perspective, many key details of allosteric mechanisms remain mysterious. Over the last decade significant attention has been paid to the “dynamic component” of allostery, as opposed to the analysis of rigid structures. Nuclear magnetic resonance spectroscopy and its ability to detect conformationally dynamic processes at atomic resolution have played an important role in expanding our understanding of allosteric mechanisms and opening up new questions. This article focuses on work that highlights how protein dynamics can factor into allosteric processes in distinct ways. Two cases are contrasted. The first considers the “traditionally allosteric” protein CheY, which undergoes a conformational change as a key element of its allostery. The second considers the more rarely observed “dynamic allostery” in a PDZ domain, in which allosteric behavior arises from changes in internal structural dynamics. Interestingly, the dynamic processes in these two contrasting examples occur on different timescales. In the case of the PDZ domain, subsequent experimental and computational work is reviewed to reveal a more complete picture of this interesting case of allostery.     

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.     

The N-terminal end of the catalytic domain of SRC kinase Hck is a conformational switch implicated in long-range allosteric regulation

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.     

Molecular dynamics simulation of shear- and stretch-induced dissociation of P-selectin/PSGL-1 complex

By mediating the tethering and rolling of leukocytes on vascular surfaces, the interactions between P-selectin and the P-selectin glycoprotein ligand 1 (PSGL-1) play crucial roles during inflammation cascade. Tensile stretch produced by rolling leukocytes and shear stress exerted by blood flow constitute the two types of mechanical forces that act on the P-selectin/PSGL-1 bond. These forces modulate not only dissociation kinetics of this bond, but also the leukocyte adhesion dynamics. However, the respective contribution of the two forces to bond dissociation and to the corresponding microstructural bases remains unclear. To mimic the mechanical microenvironment, we developed two molecular dynamics approaches; namely, an approach involving the shear flow field with a controlled velocity gradient, and the track dragging approach with a defined trajectory. With each approach or with both combined, we investigate the microstructural evolution and dissociation kinetics of the P-LE/SGP-3 construct, which is the smallest functional unit of the P-selectin/PSGL-1 complex. The results demonstrate that both shear flow and tensile stretch play important roles in the collapse of the construct and that, before bond dissociation, the former causes more destruction of domains within the construct than the latter. Dissociation of the P-LE/SGP-3 construct features intramolecular destruction of the epidermal-growth-factor (EGF) domain and the breaking of hydrogen-bond clusters at the P-selectin-lectin/EGF interface. Thus, to better understand how mechanics impacts the dissociation kinetics of the P-selectin/PSGL-1 complex, we propose herein two approaches to mimic its physiological mechanical environment.     

SH3 Domains Differentially Stimulate Distinct Dynamin I Assembly Modes and G Domain Activity

Dynamin I is a highly regulated GTPase enzyme enriched in nerve terminals which mediates vesicle fission during synaptic vesicle endocytosis. One regulatory mechanism involves its interactions with proteins containing Src homology 3 (SH3) domains. At least 30 SH3 domain-containing proteins bind dynamin at its proline-rich domain (PRD). Those that stimulate dynamin activity act by promoting its oligomerisation. We undertook a systematic parallel screening of 13 glutathione-S-transferase (GST)-tagged endocytosis-related SH3 domains on dynamin binding, GTPase activity and oligomerisation. No correlation was found between dynamin binding and their potency to stimulate GTPase activity. There was limited correlation between the extent of their ability to stimulate dynamin activity and the level of oligomerisation, indicating an as yet uncharacterised allosteric coupling of the PRD and G domain. We examined the two variants, dynamin Iab and Ibb, which differ in the alternately splice middle domain α2 helix. They responded differently to the panel of SH3s, with the extent of stimulation between the splice variants varying greatly between the SH3s. This study reveals that SH3 binding can act as a heterotropic allosteric regulator of the G domain via the middle domain α2 helix, suggesting an involvement of this helix in communicating the PRD-mediated allostery. This indicates that SH3 binding both stabilises multiple conformations of the tetrameric building block of dynamin, and promotes assembly of dynamin-SH3 complexes with distinct rates of GTP hydrolysis.     

The Role of Protein-Ligand Contacts in Allosteric Regulation of the Escherichia coli Catabolite Activator Protein

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. Both experimental and theoretical evidence demonstrate that allostery can be communicated through altered slow relaxation protein dynamics without conformational change. The catabolite activator protein (CAP) of Escherichia coli is an exemplar for the analysis of such entropically driven allostery. Negative allostery in CAP occurs between identical cAMP binding sites. Changes to the cAMP-binding pocket can therefore impact the allosteric properties of CAP. Here we demonstrate, through a combination of coarse-grained modeling, isothermal calorimetry, and structural analysis, that decreasing the affinity of CAP for cAMP enhances negative cooperativity through an entropic penalty for ligand binding. The use of variant cAMP ligands indicates the data are not explained by structural heterogeneity between protein mutants. We observe computationally that altered interaction strength between CAP and cAMP variously modifies the change in allosteric cooperativity due to second site CAP mutations. As the degree of correlated motion between the cAMP-contacting site and a second site on CAP increases, there is a tendency for computed double mutations at these sites to drive CAP toward noncooperativity. Naturally occurring pairs of covarying residues in CAP do not display this tendency, suggesting a selection pressure to fine tune allostery on changes to the CAP ligand-binding pocket without a drive to a noncooperative state. In general, we hypothesize an evolutionary selection pressure to retain slow relaxation dynamics-induced allostery in proteins in which evolution of the ligand-binding site is occurring.     

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BackgroundAccumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.Methods/resultsWe employ two computationally efficient approaches – graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.ConclusionsContact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.General significanceThis study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.     

The Role of Strong Electrostatic Interactions at the Dimer Interface of Human Glutathione Synthetase

The obligate homodimer human glutathione synthetase (hGS) provides an ideal system for exploring the role of protein–protein interactions in the structural stability, activity and allostery of enzymes. The two active sites of hGS, which are 40 Å apart, display allosteric modulation by the substrate γ-glutamylcysteine (γ-GC) during the synthesis of glutathione, a key cellular antioxidant. The two subunits interact at a relatively small dimer interface dominated by electrostatic interactions between S42, R221, and D24. Alanine scans of these sites result in enzymes with decreased activity, altered γ-GC affinity, and decreased thermal stability. Molecular dynamics simulations indicate these mutations disrupt interchain bonding and impact the tertiary structure of hGS. While the ionic hydrogen bonds and salt bridges between S42, R221, and D24 do not mediate allosteric communication in hGS, these interactions have a dramatic impact on the activity and structural stability of the enzyme.     

CR1/CR2 interactions modulate the functions of the cell surface epidermal growth factor receptor

Recent crystallographic data on the isolated extracellular domain of the epidermal growth factor receptor (EGFR) have suggested a model for its activation by ligand. We have tested this model in the context of the full-length EGFR displayed at the cell surface, by introducing mutations in two regions (CR1 and CR2) of the extracellular domain thought to be critical for regulation of receptor activation. Mutations in the CR1 and CR2 domains have opposing effects on ligand binding affinity, receptor dimerization, tyrosine kinase activation, and signaling competence. Tyr(246) is a critical residue in the CR1 loop, which is implicated in the positioning and stabilization of the receptor dimer interface after ligand binding; mutations of Tyr(246) impair or abolish receptor function. Mutations in CR2, which weaken the interaction that restricts the receptor to the tethered (inactive) state, enhance responsiveness to EGF by increasing affinity for the ligand. However, weakening of the CR1/CR2 interaction does not result in spontaneous activation of the receptors' kinase. We have used an antibody (mAb 806), which recognizes a transition state of the EGF receptor between the negatively constrained, tethered state and the fully active back-to-back dimer conformation, to follow conformational changes in the wild-type and mutant EGF receptors after ligand binding. Our results suggest that EGFR on the cell surface can be untethered, but this form is inactive; thus, untethering of the receptor is not sufficient for activation, and ligand binding is essential for the correct positioning of the two receptor subunits to achieve kinase activation.     

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.     

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct “enhanced network models”, describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. “Ab-initio network models” combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.     

Positive allosteric modulators of lecithin: Cholesterol acyltransferase adjust the orientation of the membrane-binding domain and alter its spatial free energy profile

Lecithin:cholesterol acyltransferase protein (LCAT) promotes the esterification reaction between cholesterol and phospholipid-derived acyl chains. Positive allosteric modulators have been developed to treat LCAT deficiencies and, plausibly, also cardiovascular diseases in the future. The mechanism of action of these compounds is poorly understood. Here computational docking and atomistic molecular dynamics simulations were utilized to study the interactions between LCAT and the activating compounds. Results indicate that all drugs bind to the allosteric binding pocket in the membrane-binding domain in a similar fashion. The presence of the compounds in the allosteric site results in a distinct spatial orientation and sampling of the membrane-binding domain (MBD). The MBD’s different spatial arrangement plausibly affects the lid’s movement from closed to open state and vice versa, as suggested by steered molecular dynamics simulations.     

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.     

Hinge Twists and Population Shifts Deliver Regulated Catalysis for ATP-PRT in Histidine Biosynthesis

Allosteric regulation plays an important role in the control of metabolic flux in biosynthetic pathways. In microorganisms, many enzymes in these pathways adopt different strategies of allostery to allow the tuning of their activities in response to metabolic demand. Thus, it is important to uncover the mechanism of allosteric signal transmission to fully comprehend the complex control of enzyme function and its evolution. ATP-phosphoribosyltransferase (ATP-PRT), as the first enzyme in the histidine biosynthetic pathway, is allosterically regulated by histidine and offers a good platform for the study of allostery. Two forms of ATP-PRT, namely long and short forms, were discovered that show different arrangements of their regulatory machinery. Crystal structures of the long-form ATP-PRT have revealed overall conformational changes in the inhibited state, but the observed changes in the active state are quite subtle, making the elucidation of its allosteric mechanism difficult. Here, we combine computational methods (ligand docking, quantum mechanics/molecular mechanics optimization, and molecular dynamic simulations) with experimental studies to probe the signal transmission between remote allosteric and active sites. Our results reveal that distinct conformational ensembles of the catalytic domain with different dynamic properties exist in the ligand-free and histidine-bound enzymes. These ensembles display different capabilities in supporting the catalytic and allosteric function of ATP-PRT. The findings give insight into the underlying mechanism of allostery and allow us to propose that the hinge twisting within the catalytic domain is the key for both enhancement of catalysis and provision of regulation in ATP-PRT enzymes.     

Allosteric communication in dihydrofolate reductase: signaling network and pathways for closed to occluded transition and back

Escherichia coli dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate. During the catalytic cycle, DHFR undergoes conformational transitions between the closed (CS) and occluded (OS) states that, respectively, describe whether the active site is closed or occluded by the Met20 loop. The CS→OS and the reverse transition may be viewed as allosteric transitions. Using a sequence-based approach, we identify a network of residues that represents the allostery wiring diagram. Many of the residues in the allostery wiring diagram, which are dispersed throughout the adenosine-binding domain as well as the loop domain, are not conserved. Several of the residues in the network have been previously shown by NMR experiments, mutational studies, and molecular dynamics simulations to be linked to equilibration conformational fluctuations of DHFR. To further probe the nature of events that occur during conformational fluctuations, we use a self-organized polymer model to monitor the kinetics of the CS→OS and the reverse transitions. During the CS→OS transition, coordinated changes in a number of residues in the loop domain enable the Met20 loop to slide along the α-helix in the adenosine-binding domain. Sliding is triggered by pulling of the Met20 loop by the βG-βH loop and the pushing action of the βG-βH loop. The residues that facilitate the Met20 loop motion are part of the network of residues that transmit allosteric signals during the CS→OS transition. Replacement of M16 and G121, whose C^α atoms are about 4.3 Å in the CS, by a disulfide cross-link impedes that CS→OS transition. The order of events in the OS→CS transition is not the reverse of the forward transition. The contact Glu18-Ser49 in the OS persists until the sliding of the Met20 loop is nearly complete. The ensemble of structures in the transition state in both the allosteric transitions is heterogeneous. The most probable transition-state structure resembles the OS (CS) in the CS→OS (OS→CS) transition, which is in accord with the Hammond postulate. Structures resembling the OS (CS) are present as minor (∼ 1-3%) components in equilibrated CS (OS) structures.     

Integrating molecular dynamics and co-evolutionary analysis for reliable target prediction and deregulation of the allosteric inhibition of aspartokinase for amino acid production

Deregulation of allosteric inhibition of enzymes is a challenge for strain engineering and has been achieved so far primarily by random mutation and trial-and-error. In this work, we used aspartokinase, an important allosteric enzyme for industrial amino acids production, to demonstrate a predictive approach that combines protein dynamics and evolution for a rational reengineering of enzyme allostery. Molecular dynamic simulation of aspartokinase III (AK3) from Escherichia coli and statistical coupling analysis of protein sequences of the aspartokinase family allowed to identify a cluster of residues which are correlated during protein motion and coupled during the evolution. This cluster of residues forms an interconnected network mediating the allosteric regulation, including most of the previously reported positions mutated in feedback insensitive AK3 mutants. Beyond these mutation positions, we have successfully constructed another twelve targeted mutations of AK3 desensitized toward lysine inhibition. Six threonine-insensitive mutants of aspartokinase I-homoserine dehydrogenase I (AK1-HD1) were also created based on the predictions. The proposed approach can be widely applied for the deregulation of other allosteric enzymes.