The molecular sociology of NHERF1 PDZ proteins controlling renal hormone-regulated phosphate transport

Parathyroid hormone (PTH) and fibroblast growth factor-23 (FGF23) control extracellular phosphate levels by regulating renal NPT2A-mediated phosphate transport by a process requiring the PDZ scaffold protein NHERF1. NHERF1 possesses two PDZ domains, PDZ1 and PDZ2, with identical core-binding GYGF motifs explicitly recognizing distinct binding partners that play different and specific roles in hormone-regulated phosphate transport. The interaction of PDZ1 and the carboxy-terminal PDZ-binding motif of NPT2A (C-TRL) is required for basal phosphate transport. PDZ2 is a regulatory domain that scaffolds multiple biological targets, including kinases and phosphatases involved in FGF23 and PTH signaling. FGF23 and PTH trigger disassembly of the NHERF1-NPT2A complex through reversible hormone-stimulated phosphorylation with ensuing NPT2A sequestration, down-regulation, and cessation of phosphate absorption. In the absence of NHERF1-NPT2A interaction, inhibition of FGF23 or PTH signaling results in disordered phosphate homeostasis and phosphate wasting. Additional studies are crucial to elucidate how NHERF1 spatiotemporally coordinates cellular partners to regulate extracellular phosphate levels.

Scribble scrambles parathyroid hormone receptor interactions to regulate phosphate and vitamin D homeostasis

G protein-coupled receptors, including PTHR, are pivotal for controlling metabolic processes ranging from serum phosphate and vitamin D levels to glucose uptake, and cytoplasmic interactors may modulate their signaling, trafficking, and function. We now show that direct interaction with Scribble, a cell polarity-regulating adaptor protein, modulates PTHR activity. Scribble is a crucial regulator for establishing and developing tissue architecture, and its dysregulation is involved in various disease conditions, including tumor expansion and viral infections. Scribble co-localizes with PTHR at basal and lateral surfaces in polarized cells. Using X-ray crystallography, we show that colocalization is mediated by engaging a short sequence motif at the PTHR C-terminus using Scribble PDZ1 and PDZ3 domain, with binding affinities of 31.7 and 13.4 μM, respectively. Since PTHR controls metabolic functions by actions on renal proximal tubules, we engineered mice to selectively knockout Scribble in proximal tubules. The loss of Scribble impacted serum phosphate and vitamin D levels and caused significant plasma phosphate elevation and increased aggregate vitamin D3 levels, whereas blood glucose levels remained unchanged. Collectively these results identify Scribble as a vital regulator of PTHR-mediated signaling and function. Our findings reveal an unexpected link between renal metabolism and cell polarity signaling.

Mutations in an unrecognized internal NPT2A PDZ motif disrupt phosphate transport and cause congenital hypophosphatemia

The Na+-dependent phosphate cotransporter-2A (NPT2A, SLC34A1) is a primary regulator of extracellular phosphate homeostasis. Its most prominent structural element is a carboxy-terminal PDZ ligand that binds Na+/H+ Exchanger Regulatory Factor-1 (NHERF1, SLC9A3R1). NHERF1, a multidomain PDZ protein, establishes NPT2A membrane localization and is required for hormone-inhibitable phosphate transport. NPT2A also possesses an uncharacterized internal PDZ ligand. Two recent clinical reports describe congenital hypophosphatemia in children harboring Arg495His or Arg495Cys variants within the internal PDZ motif. The wild-type internal 494TRL496 PDZ ligand binds NHERF1 PDZ2, which we consider a regulatory domain. Ablating the internal PDZ ligand with a 494AAA496 substitution blocked hormone-inhibitable phosphate transport. Complementary approaches, including CRISPR/Cas9 technology, site-directed mutagenesis, confocal microscopy, and modeling, showed that NPT2A Arg495His or Arg495Cys variants do not support PTH or FGF23 action on phosphate transport. Coimmunoprecipitation experiments indicate that both variants bind NHERF1 similarly to WT NPT2A. However, in contrast with WT NPT2A, NPT2A Arg495His, or Arg495Cys variants remain at the apical membrane and are not internalized in response to PTH. We predict that Cys or His substitution of the charged Arg495 changes the electrostatics, preventing phosphorylation of the upstream Thr494, interfering with phosphate uptake in response to hormone action, and inhibiting NPT2A trafficking. We advance a model wherein the carboxy-terminal PDZ ligand defines apical localization NPT2A, while the internal PDZ ligand is essential for hormone-triggered phosphate transport.

Mutations in an unrecognized internal NPT2A PDZ motif disrupt phosphate transport causing congenital hypophosphatemia

The Na ⁺ -dependent phosphate cotransporter-2A (NPT2A, SLC34A1) is a primary regulator of extracellular phosphate homeostasis. Its most prominent structural element is a carboxy-terminal PDZ ligand that binds Na ⁺ /H ⁺ Exchanger Regulatory Factor-1 (NHERF1, SLC9A3R1). NHERF1, a multidomain PDZ protein,establishes NPT2A membrane localization and is required for hormone-sensitive phosphate transport. NPT2A also possesses an uncharacterized internal PDZ ligand. Two recent clinical reports describe congenital hypophosphatemia in children harboring Arg ⁴⁹⁵ His or Arg ⁴⁹⁵ Cys variants within the internal PDZ motif. The wild-type internal ⁴⁹⁴ TRL ⁴⁹⁶ PDZ ligand binds NHERF1 PDZ2, which we consider a regulatory domain. Ablating the internal PDZ ligand with a ⁴⁹⁴ AAA ⁴⁹⁶ substitution blocked hormone-sensitive phosphate transport. Complementary approaches, including CRISPR/Cas9 technology, site-directed mutagenesis, confocal microscopy, and modeling, showed that NPT2A Arg ⁴⁹⁵ His or Arg ⁴⁹⁵ Cys variants do not support PTH or FGF23 action on phosphate transport. Coimmunoprecipitation experiments indicate that both variants bind NHERF1 similarly to WT NPT2A. However, in contrast to WT NPT2A, NPT2A Arg ⁴⁹⁵ His or Arg ⁴⁹⁵ Cys variants remain at the apical membrane and are not internalized in response to PTH. We predict that Cys or His substitution of the charged Arg ⁴⁹⁵ changes the electrostatics, preventing phosphorylation of the upstream Thr ⁴⁹⁴ , interfering with phosphate uptake in response to hormone action, and inhibiting NPT2A trafficking. We advance a model wherein the carboxyterminal PDZ ligand defines apical localization NPT2A, while the internal PDZ ligand is essential for hormone-triggered phosphate transport.

RGS14 regulates PTH- and FGF23-sensitive NPT2A-mediated renal phosphate uptake via binding to the NHERF1 scaffolding protein

Phosphate homeostasis, mediated by dietary intake, renal absorption, and bone deposition, is incompletely understood because of the uncharacterized roles of numerous implicated protein factors. Here, we identified a novel role for one such element, regulator of G protein signaling 14 (RGS14), suggested by genome-wide association studies to associate with dysregulated Pi levels. We show that human RGS14 possesses a carboxy-terminal PDZ ligand required for sodium phosphate cotransporter 2a (NPT2A) and sodium hydrogen exchanger regulatory factor-1 (NHERF1)-mediated renal Pi transport. In addition, we found using isotope uptake measurements combined with bioluminescence resonance energy transfer assays, siRNA knockdown, pull-down and overlay assays, and molecular modeling that secreted proteins parathyroid hormone (PTH) and fibroblast growth factor 23 inhibited Pi uptake by inducing dissociation of the NPT2A-NHERF1 complex. PTH failed to affect Pi transport in cells expressing RGS14, suggesting that it suppresses hormone-sensitive but not basal Pi uptake. Interestingly, RGS14 did not affect PTH-directed G protein activation or cAMP formation, implying a postreceptor site of action. Further pull-down experiments and direct binding assays indicated that NPT2A and RGS14 bind distinct PDZ domains on NHERF1. We showed that RGS14 expression in human renal proximal tubule epithelial cells blocked the effects of PTH and fibroblast growth factor 23 and stabilized the NPT2A-NHERF1 complex. In contrast, RGS14 genetic variants bearing mutations in the PDZ ligand disrupted RGS14 binding to NHERF1 and subsequent PTH-sensitive Pi transport. In conclusion, these findings identify RGS14 as a novel regulator of hormone-sensitive Pi transport. The results suggest that changes in RGS14 function or abundance may contribute to the hormone resistance and hyperphosphatemia observed in kidney diseases.

Multisite NHERF1 phosphorylation controls GRK6A regulation of hormone-sensitive phosphate transport

The type II sodium-dependent phosphate cotransporter (NPT2A) mediates renal phosphate uptake. The NPT2A is regulated by parathyroid hormone (PTH) and fibroblast growth factor 23, which requires Na⁺/H⁺ exchange regulatory factor-1 (NHERF1), a multidomain PDZ-containing phosphoprotein. Phosphocycling controls the association between NHERF1 and the NPT2A. Here, we characterize the critical involvement of G protein–coupled receptor kinase 6A (GRK6A) in mediating PTH-sensitive phosphate transport by targeted phosphorylation coupled with NHERF1 conformational rearrangement, which in turn allows phosphorylation at a secondary site. GRK6A, through its carboxy-terminal PDZ recognition motif, binds NHERF1 PDZ1 with greater affinity than PDZ2. However, the association between NHERF1 PDZ2 and GRK6A is necessary for PTH action. Ser¹⁶², a PKCα phosphorylation site in PDZ2, regulates the binding affinity between PDZ2 and GRK6A. Substitution of Ser¹⁶² with alanine (S¹⁶²A) blocks the PTH action but does not disrupt the interaction between NHERF1 and the NPT2A. Replacement of Ser¹⁶² with aspartic acid (S¹⁶²D) abrogates the interaction between NHERF1 and the NPT2A and concurrently PTH action. We used amber codon suppression to generate a phosphorylated Ser¹⁶²(pSer¹⁶²)-PDZ2 variant. K_D values determined by fluorescence anisotropy indicate that incorporation of pSer¹⁶² increased the binding affinity to the carboxy terminus of GRK6A 2-fold compared with WT PDZ2. Molecular dynamics simulations predict formation of an electrostatic network between pSer¹⁶² and Asp¹⁸³ of PDZ2 and Arg at position −1 of the GRK6A PDZ-binding motif. Our results suggest that PDZ2 plays a regulatory role in PTH-sensitive NPT2A-mediated phosphate transport and phosphorylation of Ser¹⁶² in PDZ2 modulates the interaction with GRK6A.

Noncanonical Sequences Involving NHERF1 Interaction with NPT2A Govern Hormone-Regulated Phosphate Transport: Binding Outside the Box

Na+/H+ exchange factor-1 (NHERF1), a multidomain PDZ scaffolding phosphoprotein, is required for the type II sodium-dependent phosphate cotransporter (NPT2A)-mediated renal phosphate absorption. Both PDZ1 and PDZ2 domains are involved in NPT2A-dependent phosphate uptake. Though harboring identical core-binding motifs, PDZ1 and PDZ2 play entirely different roles in hormone-regulated phosphate transport. PDZ1 is required for the interaction with the C-terminal PDZ-binding sequence of NPT2A (-TRL). Remarkably, phosphocycling at Ser290 distant from PDZ1, the penultimate step for both parathyroid hormone (PTH) and fibroblast growth factor-23 (FGF23) regulation, controls the association between NHERF1 and NPT2A. PDZ2 interacts with the C-terminal PDZ-recognition motif (-TRL) of G Protein-coupled Receptor Kinase 6A (GRK6A), and that promotes phosphorylation of Ser290. The compelling biological puzzle is how PDZ1 and PDZ2 with identical GYGF core-binding motifs specifically recognize distinct binding partners. Binding determinants distinct from the canonical PDZ-ligand interactions and located "outside the box" explain PDZ domain specificity. Phosphorylation of NHERF1 by diverse kinases and associated conformational changes in NHERF1 add more complexity to PDZ-binding diversity.

Parallel Post-Translational Modification Scanning Enhancing Hydrogen-Deuterium Exchange-Mass Spectrometry Coverage of Key Structural Regions

Hydrogen-deuterium exchange-mass spectrometry (HDXMS) is a powerful technology to characterize conformations and conformational dynamics of proteins and protein complexes. HDXMS has been widely used in the field of therapeutics for the development of protein drugs. Although sufficient sequence coverage is critical to the success of HDXMS, it is sometimes difficult to achieve. In this study, we developed a HDXMS data analysis strategy that includes parallel post-translational modification (PTM) scanning in HDXMS analysis. Using a membrane-delimited G protein-coupled receptor (vasopressin type 2 receptor; V₂R) and a cytosolic protein (Na⁺/H⁺ exchanger regulatory factor-1; NHERF1) as examples, we demonstrate that this strategy substantially improves protein sequence coverage, especially in key structural regions likely including PTMs themselves that play important roles in protein conformational dynamics and function.

Parathyroid hormone initiates dynamic NHERF1 phosphorylation cycling and conformational changes that regulate NPT2A-dependent phosphate transport

Na⁺-H⁺ exchanger regulatory factor-1 (NHERF1) is a PDZ protein that scaffolds membrane proteins, including sodium-phosphate co-transport protein 2A (NPT2A) at the plasma membrane. NHERF1 is a phosphoprotein with 40 Ser and Thr residues. Here, using tandem MS analysis, we characterized the sites of parathyroid hormone (PTH)-induced NHERF1 phosphorylation and identified 10 high-confidence phosphorylation sites. Ala replacement at Ser⁴⁶, Ser¹⁶², Ser¹⁸¹, Ser²⁶⁹, Ser²⁸⁰, Ser²⁹¹, Thr²⁹³, Ser²⁹⁹, and Ser³⁰² did not affect phosphate uptake, but S290A substitution abolished PTH-dependent phosphate transport. Unexpectedly, Ser²⁹⁰ was rapidly dephosphorylated and rephosphorylated after PTH stimulation, and we found that protein phosphatase 1α (PP1α), which binds NHERF1 through a conserved VxF/W PP1 motif, dephosphorylates Ser²⁹⁰ Mutating ²⁵⁷VPF²⁵⁹ eliminated PP1 binding and blunted dephosphorylation. Tautomycetin blocked PP1 activity and abrogated PTH-sensitive phosphate transport. Using fluorescence lifetime imaging (FLIM), we observed that PTH paradoxically and transiently elevates intracellular phosphate. Added phosphate blocked PP1α-mediated Ser²⁹⁰ dephosphorylation of recombinant NHERF1. Hydrogen-deuterium exchange MS revealed that β-sheets in NHERF1's PDZ2 domain display lower deuterium uptake than those in the structurally similar PDZ1, implying that PDZ1 is more cloistered. Dephosphorylated NHERF1 exhibited faster exchange at C-terminal residues suggesting that NHERF1 dephosphorylation precedes Ser²⁹⁰ rephosphorylation. Our results show that PP1α and NHERF1 form a holoenzyme and that a multiprotein kinase cascade involving G protein-coupled receptor kinase 6A controls the Ser²⁹⁰ phosphorylation status of NHERF1 and regulates PTH-sensitive, NPT2A-mediated phosphate uptake. These findings reveal how reversible phosphorylation modifies protein conformation and function and the biochemical mechanisms underlying PTH control of phosphate transport.

Origins of PDZ Binding Specificity. A Computational and Experimental Study Using NHERF1 and the Parathyroid Hormone Receptor

Na⁺/H⁺ exchanger regulatory factor-1 (NHERF1) is a scaffolding protein containing two PSD95/discs large protein/ZO1 (PDZ) domains that modifies the signaling, trafficking, and function of the parathyroid hormone receptor (PTHR), a family B G-protein-coupled receptor. PTHR and NHERF1 bind through a PDZ-ligand-recognition mechanism. We show that PTH elicits phosphorylation of Thr591 in the canonical -ETVM binding motif of PTHR. Conservative substitution of Thr591 with Cys does not affect PTH(1-34)-induced cAMP production or binding of PTHR to NHERF1. The findings suggested the presence of additional sites upstream of the PDZ-ligand motif through which the two proteins interact. Structural determinants outside the canonical NHERF1 PDZ-PTHR interface that influence binding have not been characterized. We used molecular dynamics (MD) simulation to predict residues involved in these interactions. Simulation data demonstrate that the negatively charged Glu side chains at positions -3, -5, and -6 upstream of the PDZ binding motif are involved in PDZ-PTHR recognition. Engineered mutant peptides representing the PTHR C-terminal region were used to measure the binding affinity with NHERF1 PDZ domains. Comparable micromolar affinities for peptides of different length were confirmed by fluorescence polarization, isothermal titration calorimetry, and surface plasmon resonance. Binding affinities measured for Ala variants validate MD simulations. The linear relation between the change in enthalpy and entropy following Ala substitutions at upstream positions -3, -5, and -6 of the PTHR peptide provides a clear example of the thermodynamic compensation rule. Overall, our data highlight sequences in PTHR that contribute to NHERF1 interaction and can be altered to prevent phosphorylation-mediated inhibition.

Actin-Sorting Nexin 27 (SNX27)-Retromer Complex Mediates Rapid Parathyroid Hormone Receptor Recycling

The G protein-coupled parathyroid hormone receptor (PTHR) regulates mineral-ion homeostasis and bone remodeling. Upon parathyroid hormone (PTH) stimulation, the PTHR internalizes into early endosomes and subsequently traffics to the retromer complex, a sorting platform on early endosomes that promotes recycling of surface receptors. The C terminus of the PTHR contains a type I PDZ ligand that binds PDZ domain-containing proteins. Mass spectrometry identified sorting nexin 27 (SNX27) in isolated endosomes as a PTHR binding partner. PTH treatment enriched endosomal PTHR. SNX27 contains a PDZ domain and serves as a cargo selector for the retromer complex. VPS26, VPS29, and VPS35 retromer subunits were isolated with PTHR in endosomes from cells stimulated with PTH. Molecular dynamics and protein binding studies establish that PTHR and SNX27 interactions depend on the PDZ recognition motif in PTHR and the PDZ domain of SNX27. Depletion of either SNX27 or VPS35 or actin depolymerization decreased the rate of PTHR recycling following agonist stimulation. Mutating the PDZ ligand of PTHR abolished the interaction with SNX27 but did not affect the overall rate of recycling, suggesting that PTHR may directly engage the retromer complex. Coimmunoprecipitation and overlay experiments show that both intact and mutated PTHR bind retromer through the VPS26 protomer and sequentially assemble a ternary complex with PTHR and SNX27. SNX27-independent recycling may involve N-ethylmaleimide-sensitive factor, which binds both PDZ intact and mutant PTHRs. We conclude that PTHR recycles rapidly through at least two pathways, one involving the ASRT complex of actin, SNX27, and retromer and another possibly involving N-ethylmaleimide-sensitive factor.

Phosphorylation of ezrin-radixin-moesin-binding phosphoprotein 50 (EBP50) by Akt promotes stability and mitogenic function of S-phase kinase-associated protein-2 (Skp2)

The regulation of the cell cycle by the ubiquitin-proteasome system is dependent on the activity of E3 ligases. Skp2 (S-phase kinase associated protein-2) is the substrate recognition subunit of the E3 ligase that ubiquitylates the cell cycle inhibitors p21(cip1) and p27(kip1) thus promoting cell cycle progression. Increased expression of Skp2 is frequently observed in diseases characterized by excessive cell proliferation, such as cancer and neointima hyperplasia. The stability and cellular localization of Skp2 are regulated by Akt, but the molecular mechanisms underlying these effects remain only partly understood. The scaffolding protein Ezrin-Binding Phosphoprotein of 50 kDa (EBP50) contains two PDZ domains and plays a critical role in the development of neointimal hyperplasia. Here we report that EBP50 directly binds Skp2 via its first PDZ domain. Moreover, EBP50 is phosphorylated by Akt on Thr-156 within the second PDZ domain, an event that allosterically promotes binding to Skp2. The interaction with EBP50 causes cytoplasmic localization of Skp2, increases Skp2 stability and promotes proliferation of primary vascular smooth muscle cells. Collectively, these studies define a novel regulatory mechanism contributing to aberrant cell growth and highlight the importance of scaffolding function of EBP50 in Akt-dependent cell proliferation.

Ezrin-anchored protein kinase A coordinates phosphorylation-dependent disassembly of a NHERF1 ternary complex to regulate hormone-sensitive phosphate transport

Congenital defects in the Na/H exchanger regulatory factor-1 (NHERF1) are linked to disordered phosphate homeostasis and skeletal abnormalities in humans. In the kidney, these mutations interrupt parathyroid hormone (PTH)-responsive sequestration of the renal phosphate transporter, Npt2a, with ensuing urinary phosphate wasting. We now report that NHERF1, a modular PDZ domain scaffolding protein, coordinates the assembly of an obligate ternary complex with Npt2a and the PKA-anchoring protein ezrin to facilitate PTH-responsive cAMP signaling events. Activation of ezrin-anchored PKA initiates NHERF1 phosphorylation to disassemble the ternary complex, release Npt2a, and thereby inhibit phosphate transport. Loss-of-function mutations stabilize an inactive NHERF1 conformation that we show is refractory to PKA phosphorylation and impairs assembly of the ternary complex. Compensatory mutations introduced in mutant NHERF1 re-establish the integrity of the ternary complex to permit phosphorylation of NHERF1 and rescue PTH action. These findings offer new insights into a novel macromolecular mechanism for the physiological action of a critical ternary complex, where anchored PKA coordinates the assembly and turnover of the Npt2a-NHERF1-ezrin complex.

Structural basis for NHERF1 PDZ domain binding

The Na(+)/H(+) exchange regulatory factor-1 (NHERF1) is a scaffolding protein that possesses two tandem PDZ domains and a carboxy-terminal ezrin-binding domain (EBD). The parathyroid hormone receptor (PTHR), type II sodium-dependent phosphate cotransporter (Npt2a), and β2-adrenergic receptor (β2-AR), through their respective carboxy-terminal PDZ-recognition motifs, individually interact with NHERF1 forming a complex with one of the PDZ domains. In the basal state, NHERF1 adopts a self-inhibited conformation, in which its carboxy-terminal PDZ ligand interacts with PDZ2. We applied molecular dynamics (MD) simulations to uncover the structural and biochemical basis for the binding selectivity of NHERF1 PDZ domains. PDZ1 uniquely forms several contacts not present in PDZ2 that further stabilize PDZ1 interactions with target ligands. The binding free energy (ΔG) of PDZ1 and PDZ2 with the carboxy-terminal, five-amino acid residues that form the PDZ-recognition motif of PTHR, Npt2a, and β2-AR was calculated and compared with the calculated ΔG for the self-association of NHERF1. The results suggest that the interaction of the PTHR, β2-adrenergic, and Npt2a involves competition between NHERF1 PDZ domains and the target proteins. The binding of PDZ2 with PTHR may also compete with the self-inhibited conformation of NHERF1, thereby contributing to the stabilization of an active NHERF1 conformation.

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.

Structure and energetics of channel-forming protein-polysaccharide complexes inferred via computational statistical thermodynamics

The ion channel protein alpha-hemolysin (alphaHL) forms supramolecular complexes with the polysaccharide beta-cyclodextrin (betaCD). This system has potential uses in nanoscale device engineering. It has been found recently that betaCD formed longer- or shorter-lived complexes with some engineered alphaHL mutants then with a wild type protein (Gu et al. J. Gen. Physiol. 2001, 118, 481-493). However, how changes in the protein sequence affect complex lifetime was not completely understood in part due to the lack of knowledge of structures of these metastable complexes. In this paper, we present an extensive molecular modeling study of the betaCD-alphaHL and selected mutant complexes to gain insights into the betaCD-alphaHL interaction mechanisms and to predict possible structures and energetics of the complexes. Thermodynamic integration (TI) and umbrella sampling (US) techniques (with the weighted histogram analysis method (WHAM)) were used to calculate the relative binding affinities of the complexes formed with the wild type alphaHL and the M113N, M113E, M113A, and M113V mutants. Our results are in excellent agreement with experiment. While betaCD-M113N and betaCD-M113A complexes were stable in the configuration of the wild type complex, the equilibrium configuration of the betaCD-M113V and betaCD-M113E complexes was significantly different. In these cases, TI alone was insufficient to accurately calculate the corresponding free energy differences. By utilizing a TI/US combination in a novel manner, we were able to accurately calculate free energy changes in these flexible systems. The betaCD-M113A and betaCD-M113E complexes, which exhibited shorter lifetimes than other complexes in an experiment, in simulations exhibited greater flexibility and higher water solvation of the betaCD adapter. MD simulations of the betaCD-M113N complex with betaCD in a downward orientation were also performed.

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.