Role of the hexapeptide disulfide loop in the gamma-carboxyglutamic acid domain of protein C in Ca2+-mediated structural and functional properties

The anticoagulant and immunomodulatory effects of protein C (PC) rely on the presence of the N-terminal gamma-carboxyglutamic acid (Gla) domain. This domain is strongly conserved among vitamin K-dependent blood proteins and, in addition to a high relative content of Gla, contains a hexapeptide disulfide loop between Cys residues 17 and 22. In the present study, the contribution of the hexapeptide loop toward Gla domain structure and function was evaluated using wild-type and Cys17/Cys22-alkylated synthetic peptide analogues of the 47-residue Gla domain/helical stack of PC. Circular dichroism and intrinsic fluorescence measurements revealed significant differences in the metal ion-dependent conformations of the two peptides. Disruption of the disulfide loop slightly altered the capacity of the peptide to interact with acidic phospholipid (PL) vesicles. The affinity of the alkylated peptide for soluble endothelial protein C receptor (EPCR), as demonstrated by surface plasmon resonance studies, was increased compared with the wild-type species, although total binding was compromised. These results suggest that the disulfide loop of PC contributes to the overall Ca(2+)-dependent conformation but is not strictly required for PL membrane binding or EPCR recognition.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Multifunctional specificity of the protein C/activated protein C Gla domain

Activated protein C (APC) has potent anticoagulant and anti-inflammatory properties that are mediated in part by its interactions with its cofactor protein S and the endothelial cell protein C receptor (EPCR). The protein C/APC Gla domain is implicated in both interactions. We sought to identify how the protein C Gla domain enables specific protein-protein interactions in addition to its conserved role in phospholipid binding. The human prothrombin Gla domain, which cannot bind EPCR or support protein S cofactor activity, has 22/45 residues that are not shared with the human protein C Gla domain. We hypothesized that the unique protein C/APC Gla domain residues were responsible for mediating the specific interactions. To assess this, we generated 13 recombinant protein C/APC variants incorporating the prothrombin residue substitutions. Despite anticoagulant activity similar to wild-type APC in the absence of protein S, APC variants APC(PT33-39) (N33S/V34S/D35T/D36A/L38D/A39V) and APC(PT36/38/39) (D36A/L38D/A39V) were not stimulated by protein S, whereas APC(PT35/36) (D35T/D36A) exhibited reduced protein S sensitivity. Moreover, PC(PT8/10) (L8V/H10K) displayed negligible EPCR affinity, despite normal binding to anionic phospholipid vesicles and factor Va proteolysis in the presence and absence of protein S. A single residue variant, PC(PT8), also failed to bind EPCR. Factor VIIa, which also possesses Leu-8, bound soluble EPCR with similar affinity to wild-type protein C, collectively confirming Leu-8 as the critical residue for EPCR recognition. These results reveal the specific Gla domain residues responsible for mediating protein C/APC molecular recognition with both its cofactor and receptor and further illustrate the multifunctional potential of Gla domains.     

Protein-protein docking with multiple residue conformations and residue substitutions

The protein docking problem has two major aspects: sampling conformations and orientations, and scoring them for fit. To investigate the extent to which the protein docking problem may be attributed to the sampling of ligand side‐chain conformations, multiple conformations of multiple residues were calculated for the uncomplexed (unbound) structures of protein ligands. These ligand conformations were docked into both the complexed (bound) and unbound conformations of the cognate receptors, and their energies were evaluated using an atomistic potential function. The following questions were considered: (1) does the ensemble of precalculated ligand conformations contain a structure similar to the bound form of the ligand? (2) Can the large number of conformations that are calculated be efficiently docked into the receptors? (3) Can near‐native complexes be distinguished from non‐native complexes? Results from seven test systems suggest that the precalculated ensembles do include side‐chain conformations similar to those adopted in the experimental complexes. By assuming additivity among the side chains, the ensemble can be docked in less than 12 h on a desktop computer. These multiconformer dockings produce near‐native complexes and also non‐native complexes. When docked against the bound conformations of the receptors, the near‐native complexes of the unbound ligand were always distinguishable from the non‐native complexes. When docked against the unbound conformations of the receptors, the near‐native dockings could usually, but not always, be distinguished from the non‐native complexes. In every case, docking the unbound ligands with flexible side chains led to better energies and a better distinction between near‐native and non‐native fits. An extension of this algorithm allowed for docking multiple residue substitutions (mutants) in addition to multiple conformations. The rankings of the docked mutant proteins correlated with experimental binding affinities. These results suggest that sampling multiple residue conformations and residue substitutions of the unbound ligand contributes to, but does not fully provide, a solution to the protein docking problem. Conformational sampling allows a classical atomistic scoring function to be used; such a function may contribute to better selectivity between near‐native and non‐native complexes. Allowing for receptor flexibility may further extend these results.     

The crystal structure of the endothelial protein C receptor and a bound phospholipid

The endothelial cell protein C receptor (EPCR) shares approximately 20% sequence identity with the major histocompatibility complex class 1/CD1 family of molecules, accelerates the thrombin-thrombomodulin-dependent generation of activated protein C, a natural anticoagulant, binds to activated neutrophils, and can undergo translocation from the plasma membrane to the nucleus. Blocking protein C/activated protein C binding to the receptor inhibits not only protein C activation but the ability of the host to respond appropriately to bacterial challenge, exacerbating both the coagulant and inflammatory responses. To understand how EPCR accomplishes these multiple tasks, we solved the crystal structure of EPCR alone and in complex with the phospholipid binding domain of protein C. The structures were strikingly similar to CD1d. A tightly bound phospholipid resides in the groove typically involved in antigen presentation. The protein C binding site is outside this conserved groove and is distal from the membrane-spanning domain. Extraction of the lipid resulted in loss of protein C binding, which could be restored by lipid reconstitution. CD1d augments the immune response by presenting glycolipid antigens. The EPCR structure is a model for how CD1d binds lipids and further suggests additional potential functions for EPCR in immune regulation, possibly including the anti-phospholipid syndrome.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Identification of protein functional surfaces by the concept of a split pocket

The function of a protein is often fulfilled via molecular interactions on its surfaces, so identifying the functional surface(s) of a protein is helpful for understanding its function. Here, we introduce the concept of a split pocket, which is a pocket that is split by a cognate ligand. We use a geometric approach that is site‐specific. Specifically, we first compute a set of all pockets in the protein with its ligand(s) and a set of all pockets with the ligand(s) removed and then compare the two sets of pockets to identify the split pocket(s) of the protein. To reduce the search space and expedite the process of surface partitioning, we design probe radii according to the physicochemical textures of molecules. Our method achieves a success rate of 96% on a benchmark test set. We conduct a large‐scale computation to identify ～19,000 split pockets from 11,328 structures (1.16 million potential pockets); for each pocket, we obtain residue composition, solvent‐accessible area, and molecular volume. With this database of split pockets, our method can be used to predict the functional surfaces of unbound structures. Indeed, the functional surface of an unbound protein may often be found from its similarity to remotely related bound forms that belong to distinct folds. Finally, we apply our method to identify glucose‐binding proteins, including unbound structures. Our study demonstrates the power of geometric and evolutionary matching for studying protein functional evolution and provides a framework for classifying protein functions by local spatial patterns of functional surfaces. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Molecular dynamics calculations of wild type vs. mutant protein C: relationship between binding affinity to endothelial cell protein C receptor and hereditary disease

Molecular dynamics simulations of the protein C gamma-carboxyglutamic acid (Gla) domain and endothelial cell protein C receptor (EPCR) complex were performed to determine the effect of a hereditary disease, which results in a mutation (Gla 25 --> Lys) in the protein C Gla domain. Our results suggest that the Gla 25 --> Lys mutation causes a significant reduction in the binding force between protein C Gla domain and EPCR due to destabilization of the helix structure of EPCR and displacement of a Ca2+ ion.     

Molecular Dynamics Calculations of Wild Type vs. Mutant Protein C: Relationship Between Binding Affinity to Endothelial Cell Protein C Receptor and Hereditary Disease

Abstract

Molecular dynamics simulations of the protein C γ-carboxyglutamic acid (Gla) domain and endothelial cell protein C receptor (EPCR) complex were performed to determine the effect of a hereditary disease, which results in a mutation (Gla 25 → Lys) in the protein C Gla domain. Our results suggest that the Gla 25 → Lys mutation causes a significant reduction in the binding force between protein C Gla domain and EPCR due to destabilization of the helix structure of EPCR and displacement of a Ca²⁺ ion.     

Lipid a from lipopolysaccharide recognition: Structure,dynamics and cooperativity by molecular dynamics simulations

Molecular dynamics simulations of Lipid A and its natural precursor Lipid IVA from E.coli have been carried out free in solution, bound to the myeliod differentiation protein 2 (MD2) and in the complex of MD2 with the toll like receptor 4 (TLR4). In addition, simulations of the ligand free MD2 and MD2‐TLR4 complex were performed. A structural and energetic characterization of the bound and unbound states of Lipid A/IVA was generated. As the crystal structures depict, the main driving force for MD2‐Lipid A/IVA are the hydrophobic interactions between the aliphatic tails and the MD2 cavity. The charged phosphate groups do strongly interact with positively charged residues, located at the surface of MD2. However, they are not essential for keeping the lipids in the cavity, indicating a more prominent role in binding recognition and ionic interactions with TLR4 at the MD2/TLR4 interface. Interestingly, in the absence of any ligand MD2 rapidly closes, blocking the binding cavity. The presence of TLR4, though changing the dynamics, was not able to impede the aforementioned closing event. We hypothesize that fluctuations of the H1 region are essential for this phenomenon, and it is plausible that an equilibrium between the open and closed states exists, although the lengths of our simulations are not sufficient to encompass the reversible process. The MD2/Lipid A‐TLR4 complex simulations show that the presence of the ligand energetically stabilizes the complex relative to the ligand‐free structures, indicating cooperativity in the binding process. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Identification of the protein C/activated protein C binding sites on the endothelial cell protein C receptor. Implications for a novel mode of ligand recognition by a major histocompatibility complex class 1-type receptor

The endothelial cell protein C receptor (EPCR) is an endothelial cell-specific transmembrane protein that binds both protein C and activated protein C (APC). EPCR regulates the protein C anticoagulant pathway by binding protein C and augmenting protein C activation by the thrombin-thrombomodulin complex. EPCR is homologous to the MHC class 1/CD1 family, members of which contain two alpha-helices that sit upon an 8-stranded beta-sheet platform. In this study, we identified 10 residues that, when mutated to alanine, result in the loss of protein C/APC binding (Arg-81, Leu-82, Val-83, Glu-86, Arg-87, Phe-146, Tyr-154, Thr-157, Arg-158, and Glu-160). Glutamine substitutions at the four N-linked carbohydrate attachment sites of EPCR have little affect on APC binding, suggesting that the carbohydrate moieties of EPCR are not critical for ligand recognition. We then mapped the epitopes for four anti-human EPCR monoclonal antibodies (mAbs), two of which block EPCR/Fl-APC (APC labeled at the active site with fluorescein) interactions, whereas two do not. These epitopes were localized by generating human-mouse EPCR chimeric proteins, since the mAbs under investigation do not recognize mouse EPCR. We found that 5 of the 10 candidate residues for protein C/APC binding (Arg-81, Leu-82, Val-83, Glu-86, Arg-87) colocalize with the epitope for one of the blocking mAbs. Three-dimensional molecular modeling of EPCR indicates that the 10 protein C/APC binding candidate residues are clustered at the distal end of the two alpha-helical segments. Protein C activation studies on 293 cells that coexpress EPCR variants and thrombomodulin demonstrate that protein C binding to EPCR is necessary for the EPCR-dependent enhancement in protein activation by the thrombin-thrombomodulin complex. These studies indicate that EPCR has exploited the MHC class 1 fold for an alternative and possibly novel mode of ligand recognition. These studies are also the first to identify the protein C/APC binding region of EPCR and may provide useful information about molecular defects in EPCR that could contribute to cardiovascular disease susceptibility.     

Druggable Protein Interaction Sites Are More Predisposed to Surface Pocket Formation than the Rest of the Protein Surface

Despite intense interest and considerable effort via high-throughput screening, there are few examples of small molecules that directly inhibit protein-protein interactions. This suggests that many protein interaction surfaces may not be intrinsically “druggable” by small molecules, and elevates in importance the few successful examples as model systems for improving our fundamental understanding of druggability. Here we describe an approach for exploring protein fluctuations enriched in conformations containing surface pockets suitable for small molecule binding. Starting from a set of seven unbound protein structures, we find that the presence of low-energy pocket-containing conformations is indeed a signature of druggable protein interaction sites and that analogous surface pockets are not formed elsewhere on the protein. We further find that ensembles of conformations generated with this biased approach structurally resemble known inhibitor-bound structures more closely than equivalent ensembles of unbiased conformations. Collectively these results suggest that “druggability” is a property encoded on a protein surface through its propensity to form pockets, and inspire a model in which the crude features of the predisposed pocket(s) restrict the range of complementary ligands; additional smaller conformational changes then respond to details of a particular ligand. We anticipate that the insights described here will prove useful in selecting protein targets for therapeutic intervention.     

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.     

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

Monte Carlo replica‐exchange based ensemble docking of protein conformations

A replica‐exchange Monte Carlo (REMC) ensemble docking approach has been developed that allows efficient exploration of protein–protein docking geometries. In addition to Monte Carlo steps in translation and orientation of binding partners, possible conformational changes upon binding are included based on Monte Carlo selection of protein conformations stored as ordered pregenerated conformational ensembles. The conformational ensembles of each binding partner protein were generated by three different approaches starting from the unbound partner protein structure with a range spanning a root mean square deviation of 1–2.5 Å with respect to the unbound structure. Because MC sampling is performed to select appropriate partner conformations on the fly the approach is not limited by the number of conformations in the ensemble compared to ensemble docking of each conformer pair in ensemble cross docking. Although only a fraction of generated conformers was in closer agreement with the bound structure the REMC ensemble docking approach achieved improved docking results compared to REMC docking with only the unbound partner structures or using docking energy minimization methods. The approach has significant potential for further improvement in combination with more realistic structural ensembles and better docking scoring functions. Proteins 2017; 85:924–937. © 2016 Wiley Periodicals, Inc.