Prediction of solution structures of the Ca2+-bound gamma-carboxyglutamic acid domains of protein S and homolog growth arrest specific protein 6: use of the particle mesh Ewald method

Role of TAM Receptors in Antimalarial Humoral Immune Response

Immune response against malaria and the clearance of Plasmodium parasite relies on germinal-center-derived B cell responses that are temporally and histologically layered. Despite a well-orchestrated germinal center response, anti-Plasmodium immune response seldom offers sterilizing immunity. Recent studies report that certain pathophysiological features of malaria such as extensive hemolysis, hypoxia as well as the extrafollicular accumulation of short-lived plasmablasts may contribute to this suboptimal immune response. In this review, we summarize some of those studies and attempt to connect certain host intrinsic features in response to the malarial disease and the resultant gaps in the immune response.

Therapeutic targeting of the functionally elusive TAM receptor family

The TAM receptor family of TYRO3, AXL and MERTK regulates tissue and immune homeostasis. Aberrant TAM receptor signalling has been linked to a range of diseases, including cancer, fibrosis and viral infections. Specifically, the dysregulation of TAM receptors can enhance tumour growth and metastasis due to their involvement in multiple oncogenic pathways. For example, TAM receptors have been implicated in the epithelial-mesenchymal transition, maintaining the stem cell phenotype, immune modulation, proliferation, angiogenesis and resistance to conventional and targeted therapies. Therapeutically, multiple TAM receptor inhibitors are in preclinical and clinical development for cancers and other indications, with those targeting AXL being the most clinically advanced. Although there has been notable clinical advancement in recent years, challenges persist. This Review aims to provide both biological and clinical insights into the current therapeutic landscape of TAM receptor inhibitors, and evaluates their potential for the treatment of cancer and non-malignant diseases.

Identification of novel inhibitor of protein tyrosine phosphatases delta: structure-based pharmacophore modeling, virtual screening, flexible docking, molecular dynamics simulation, and post-molecular dynamics analysis

Owing to their unique functions in regulating the synapse activity of protein tyrosine phosphatases delta (PTPδ) that has drawn special attention for developing drugs to autism spectrum disorders (ASDs). In this study, the PTPδ pharmacophore was first established by the structure-based pharmacophore method. Subsequently, 10 compounds contented Lipinski's rule of five was acquired by the virtual screening of the PTPδ pharmacophore against ZINC and PubChem databases. Then, the 10 identified molecules were discovered that had better binding affinity than a known PTPδ inhibitors compound SCHEMBL16375396. Two compounds SCHEMBL16375408 and ZINC19796658 with high binding score, low toxicity were gained. They were observed by docking analysis and molecular dynamics simulations that the novel potential inhibitors not only possessed the same function as SCHEMBL16375396 did in inhibiting PTPδ, but also had more favorable conformation to bind with the catalytic active regions. This study provides a new method for identify PTPδ inhibitor for the treatment of ASDs disease.Communicated by Ramaswamy H. Sarma.

The Emerging Role of TYRO3 as a Therapeutic Target in Cancer

The TAM family (TYRO3, AXL, MERTK) tyrosine kinases play roles in diverse biological processes including immune regulation, clearance of apoptotic cells, platelet aggregation, and cell proliferation, survival, and migration. While AXL and MERTK have been extensively studied, less is known about TYRO3. Recent studies revealed roles for TYRO3 in cancer and suggest TYRO3 as a therapeutic target in this context. TYRO3 is overexpressed in many types of cancer and functions to promote tumor cell survival and/or proliferation, metastasis, and resistance to chemotherapy. In addition, higher levels of TYRO3 expression have been associated with decreased overall survival in patients with colorectal, hepatocellular, and breast cancers. Here we review the physiological roles for TYRO3 and its expression and functions in cancer cells and the tumor microenvironment, with emphasis on the signaling pathways that are regulated downstream of TYRO3 and emerging roles for TYRO3 in the immune system. Translational agents that target TYRO3 are also described.

Gene of the month: Axl

The interaction between Axl receptor tyrosine kinase and its main ligand Gas6 has been implicated in the progression of a wide number of malignancies. More recently, overexpression of Axl has emerged as a key molecular determinant underlying the development of acquired resistance to targeted anticancer agents. The activation of Axl is overexpression-dependent and controls a number of hallmarks of cancer progression including proliferation, migration, resistance to apoptosis and survival through a complex network of intracellular second messengers. Axl has been noted to influence clinically meaningful end points including metastatic recurrence and survival in the vast majority of tumour types. With Axl inhibitors having gained momentum as novel anticancer therapies, we provide an overview of the biological and clinical relevance of this molecular pathway, outlining the main directions of research.

The role of Ca(2+) ions in the complex assembling of protein Z and Z-dependent protease inhibitor: A structure and dynamics investigation

We investigated the solution structure and dynamics of the human anti-coagulation protein Z (PZ) in the complex with protein Zdependent protease inhibitor (ZPI) to order to understand key structural changes in the presence and absence of Ca(2+). Structural features of the complete complex of PZ-ZPI are poorly understood due to lack of complete atomic model of the PZ-ZPI complex. We have constructed a model of the complete PZ-ZPI complex and molecular dynamics (MD) simulation of the solvated PZ-ZPI complex with and without Ca(2+) was achieved for 100ns. It is consider that the Ω-loop of GLA domains interacts with negatively charged biological membranes in the presence of Ca(2+) ions. The PZ exerts its role as cofactor in a similar way. However, we used solvent-equilibrated dynamics to show structural features of the PZ-ZPI complex in the presence and the absence of Ca(2+)ions. We observed that the distance between the interacting sites of the ZPI with the PZ and the GLA domain decreases in the presence of Ca(2+) ions. Further, we postulated that the calculated distance between the dominant plane of the Ca(2+) ions and Ser196 of the pseudo-catalytic triad of the PZ is similar to the equivalent distance of FXa. This suggests that the central role of the PZ in the blood coagulation may be to align the inhibitory site of the ZPI with the active site of the FXa, which is depends on the interaction of the calcium bound GLA domain of the PZ with the active membrane.

Molecular dynamics simulations of a reversibly folding beta-heptapeptide in methanol: influence of the treatment of long-range electrostatic interactions

Eight 100-ns molecular dynamics simulations of a beta-heptapeptide in methanol at 340 K (within cubic periodic computational boxes of about 6-nm edge) are reported and compared. These simulations were performed with three different charge-state combinations at the peptide termini, one of them with or without a neutralizing chloride counterion, and using either the lattice-sum (LS) or reaction-field (RF) scheme to handle electrostatic interactions. The choice of the electrostatic scheme has essentially no influence on the folding-unfolding equilibrium when the peptide termini are uncharged and only a small influence when the peptide is positively charged at its N-terminus (with or without inclusion of a neutralizing chloride counterion). However, when the peptide is zwitterionic, the LS scheme leads to preferential sampling of the high-dipole folded helical state, whereas the RF scheme leads to preferential sampling of a low-dipole unfolded salt-bridged state. A continuum electrostatics analysis based on the sampled configurations (zwitterionic case) suggests that the LS scheme stabilizes the helical state through artificial periodicity, but that the magnitude of this perturbation is essentially negligible (compared to the thermal energy) for the large box size and relatively polar solvent considered. The results thus provide clear evidence (continuum electrostatics analysis) for the absence of LS artifacts and some indications (still not definitive because of the limited sampling of the folding-unfolding transition) for the presence of RF artifacts in this specific system.

Gas6 and protein S. Vitamin K-dependent ligands for the Axl receptor tyrosine kinase subfamily

Gas6 and protein S are two homologous secreted proteins that depend on vitamin K for their execution of a range of biological functions. A discrete subset of these functions is mediated through their binding to and activation of the receptor tyrosine kinases Axl, Sky and Mer. Furthermore, a hallmark of the Gas6-Axl system is the unique ability of Gas6 and protein S to tether their non receptor-binding regions to the negatively charged membranes of apoptotic cells. Numerous studies have shown the Gas6-Axl system to regulate cell survival, proliferation, migration, adhesion and phagocytosis. Consequently, altered activity/expression of its components has been detected in a variety of pathologies such as cancer and vascular, autoimmune and kidney disorders. Moreover, Axl overactivation can equally occur without ligand binding, which has implications for tumorigenesis. Further knowledge of this exquisite ligand-receptor system and the circumstances of its activation should provide the basis for development of novel therapies for the above diseases.

Vitamin K-dependent Gas6 activates ERK kinase and stimulates growth of cardiac fibroblasts

The protein product of growth arrest specific gene 6 (Gas6), is the biological ligand for the Axl subfamily of receptor tyrosine kinases. We investigated the effects of exogenous Gas6 on growth of cardiac fibroblasts isolated from genetically Gas6-deficient mice. Recombinant Gas6, containing N terminal gamma-carboxyglutamic acid residues formed from a vitamin K-dependent reaction, stimulated both DNA synthesis and proliferation of cardiac fibroblasts under serum-free conditions. Gas6 also markedly enhanced survival of cells during prolonged serum starvation. Gas6 stimulated tyrosine phosphorylation of Axl as well as phosphorylation of ERK kinase. The mitogenic effects of Gas6 were inhibited by neutralising anti-Gas6 antibodies and by a soluble Axl ectodomain fusion protein. In contrast, recombinant Gas6 from cells treated with warfarin, which prevents the gamma-carboxylation reaction, neither stimulated fibroblast proliferation nor activated Axl tyrosine phosphorylation. Gas6-induced cell proliferation was additive to the effects of epidermal growth factor, suggesting activation of discrete signalling pathways. In conclusion, Gas6 appears to be a unique growth factor for fibroblasts and post-translational gamma-carboxylation is necessary for its biological activity. These findings implicate vitamin K-dependent biochemical reactions in growth processes in development and in disease.

Predicted solution structure of zymogen human coagulation FVII

A model solution structure for the complete tissue factor-free calcium ion-bound human zymogen FVII (residues 1-406) (FVII) has been constructed to study possible conformational changes associated with the activation process and tissue factor (TF) binding. The initial structure for the present model was constructed using the X-ray crystallographic structure of human coagulation FVIIa/TF complex bound with calcium ions (Banner et al., Nature 1996, 380, 41-46). This model was subsequently subjected to lengthy molecular dynamics simulations. The Amber force field in conjunction with the PME electrostatic summation method was employed. The estimated TF free solution structure was then compared with the currently available X-ray crystal structures of FVIIa (with or without TF, variable inhibitor bound) to estimate the restructuring of FVII due to TF binding and activation. The solution structure of the zymogen FVII in the absence of TF is predicted to be an extended domain structure similar to that of the TF-bound X-ray crystal structure. An additional extension of the serine protease (SP) domain of the zymogen above a reference lipid surface by approximately 7 A was in agreement with experiment. Significant Gla-EGF1 and EGF1-EGF2 interdomain motions in the zymogen were observed. Carbohydrate dimers attached to Ser-52 and Ser-60 did not cause restructuring in this domain. Minimal restructuring of the SP domain is found upon inference of the zymogen from the activated form. The catalytic triad residues maintain the H-bonded network while Lys-341 occupies the S1 specific site in the zymogen.

Four loops of the catalytic domain of factor viia mediate the effect of the first EGF-like domain substitution on factor viia catalytic activity

The presence of tissue factor is essential for factor VIIa (FVIIa) to reach its full catalytic potential. The previous work in this laboratory demonstrated that substitution of the EGF1 domain of factor VIIa with that of factor IX (FVII((IXegf1))a) results in a substantial decrease in TF-binding affinity and catalytic activity. Supporting simulations of the solution structures of Ca(2+)-bound factor VIIa and FVII((IXegf1))a with tissue factor are provided. Mutants are generated, based on the simulation model, to study the effect of EGF1 substitution on catalytic activity. The simulations show larger Gla-EGF1 and EGF1-EGF2 inter-domain motions for FVII((IXegf1))a than for factor VIIa. The catalytic domain of the chimeric factor VIIa has been disturbed and several surface loops in the catalytic domain of FVII((IXegf1))a (Loop 170s (170-182), Loop 1 (185-188) and Loop 2 (221A-225)) manifest larger position fluctuations than wild-type. The position of Loop 140s (142-152) of FVII((IXegf1))a, near the N terminus insertion site of the catalytic domain, shifts relative to factor VIIa, resulting in a slight alteration of the active site. The results suggest that these four loops mediate the effect of the EGF1 domain substitution on the S1 site and catalytic residues. To test the model, we prepared mutations of these surface loops, including four FVII mutants, D186A, K188A, L144A and R147A, a FVII mutant with multiple mutations (MM3: L144A+R147A+D186A) and a FVII mutant with Loop 170s partially deleted, Loop 170s(del). The catalytic activities towards a small peptidyl substrate decreased 2.4, 4.5 and 9-fold for Loop 170s(del)a (a, activated), L144Aa and D186Aa, respectively, while MM3a lost almost all catalytic activity. The combined results of the simulations and mutants provide insight into the mechanism by which tissue factor enhances factor VIIa catalytic activity.

Modeling zymogen protein C

A solution structure for the complete zymogen form of human coagulation protein C is modeled. The initial core structure is based on the x-ray crystallographic structure of the gamma-carboxyglutamic acid (Gla)-domainless activated form. The Gla domain (residues 1-48) is modeled from the x-ray crystal coordinates of the factor VII(a)/tissue factor complex and oriented with the epidermal growth factor-1 domain to yield an initial orientation consistent with the x-ray crystal structure of porcine factor IX(a). The missing C-terminal residues in the light chain (residues 147-157) and the activation peptide residues 158-169 were introduced using homology modeling so that the activation peptide residues directly interact with the residues in the calcium binding loop. Molecular dynamics simulations (Amber-particle-mesh-Ewald) are used to obtain the complete calcium-complexed solution structure. The individual domain structures of protein C in solution are largely unaffected by solvation, whereas the Gla-epidermal growth factor-1 orientation evolves to a form different from both factors VII(a) and IX(a). The solution structure of the zymogen protein C is compared with the crystal structures of the existing zymogen serine proteases: chymotrypsinogen, proproteinase, and prethrombin-2. Calculated electrostatic potential surfaces support the involvement of the serine protease calcium ion binding loop in providing a suitable electrostatic environment around the scissile bond for II(a)/thrombomodulin interaction.

Probing the structural changes in the light chain of human coagulation factor VIIa due to tissue factor association

The crystallographic structure of human coagulation factor VIIa/tissue factor complex bound with calcium ions was used to model the solution structure of the light chain of factor VIIa (residues 1-142) in the absence of tissue factor. The Amber force field in conjunction with the particle mesh Ewald summation method to accommodate long-range electrostatic interactions was used in the trajectory calculations. The estimated TF-free solution structure was then compared with the crystal structure of factor VIIa/tissue factor complex to estimate the restructuring of factor VIIa due to tissue factor binding. The solution structure of the light chain of factor VIIa in the absence of tissue factor is predicted to be an extended domain structure similar to that of the tissue factor-bound crystal. Removal of the EGF1-bound calcium ion is shown by simulation to lead to minor structural changes within the EGF1 domain, but also leads to substantial relative reorientation of the Gla and EGF1 domains.

Trans-cis isomerization of proline 22 in bovine prothrombin fragment 1: a surprising result of structural characterization

The calcium ion-mediated interaction of bovine prothrombin (BF1) with negatively charged phospholipid membranes is assumed to be largely via the Gla domain of BF1 with the fold of the Gla domain essential for binding. It has been reported that Pro22 undergoes classical trans to cis isomerization in the presence of calcium ions with the cis conformation of Pro22 of BF1 responsible for membrane binding [Evans, T. C., Jr., and Nelsestuen, G. L. (1996) Biochemistry 35, 8210-8215]. However, Pro22 was found to be in the trans conformation in the crystal structure of BF1. In the present work, we have used molecular dynamics simulations to investigate the relative importance of the two conformations of Pro22 to the structural and dynamical properties of BF1. The initial trans conformation of Pro22 in BF1 was slowly converted to cis-Pro22 using constrained dynamics. The second-generation AMBER force field in conjunction with the particle mesh Ewald method to accommodate long-range interaction was employed in the trajectory calculations. Comparison of the BF1(trans-Pro22) and BF1(cis-Pro22) equilibrated structures reveals surprisingly that the overall structural changes associated with the trans-cis isomerization is minimal and only minor modifications to the hydrogen bond network and the network of N-terminus Ala1 take place. The calculated electrostatic potential energy surfaces of the two protein structures also appear to be very similar, indicating the near equality of the local interaction site environments in the protein prior to lipid binding.