Homology modeling and molecular dynamics simulation of human prothrombin fragment 1

Face-specific binding of prothrombin fragment 1 and human serum albumin to inorganic and urinary calcium oxalate monohydrate crystals

OBJECTIVE

To compare the intracrystalline distributions of prothrombin fragment 1 (PTF1) and human serum albumin (HSA) within inorganic and urinary calcium oxalate (CaOx) monohydrate (COM) crystals and to determine whether binding of PTF1 can be explained by interactions between particular gamma-carboxyglutamic (Gla) residues and atomic arrays on individual faces of the COM crystal. MATERIALS AND METHODS COM: crystals were precipitated from inorganic solutions and ultrafiltered urine containing fluorescent HSA or PTF1 at different relative concentrations and examined by fluorescence microscopy. Accelrys Materials Studio and Discovery Studio were used to model the binding of PTF1 to the top, side and apical faces of the COM crystal.

RESULTS

PTF1 alone always adsorbed predominantly to the COM apical surfaces, while HSA bound principally to the side faces under inorganic conditions, but to the apical faces in urine. In the presence of each other, both proteins competed for adsorption to the apical faces, with attachment of PTF1 dominating over that of HSA. Modelling showed that urinary PTF1 had equal theoretical bonding potential for all three COM surfaces.

CONCLUSIONS

(i) Anisotropic inclusion of HSA and PTF1 into urinary and inorganic COM crystals results from their preferential binding to specific COM faces; (ii) the binding preference of HSA differs under inorganic and urinary conditions; (iii) preferential binding of PTF1 to the apical faces of COM is more complex than can be explained by interactions between Gla groups and surface atomic arrays; (iv) future studies of interactions between urinary proteins and stone mineral crystal surfaces should be performed in urine.

Structure modeling, ligand binding, and binding affinity calculation (LR-MM-PBSA) of human heparanase for inhibition and drug design

Human heparanase is an endo-beta-D-glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three-dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three-dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three-dimensional structure of heparanase using a homology-modeling approach. The homology-built structure is consistent to previous bioinformatics and site-mutation experimental results. The heparanase features a (alpha/beta)(8) TIM-barrel fold with two glutamate residues (Glu225 and Glu343) located in the active-site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM-PBSA approach (LR-MM-PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave-one-out cross-validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein "induced-fit" that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme.

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.

Crystal±macromolecule interactions in urolithiasis: lessons from healthy biomineralization systems

Recent advances in understanding the molecular level mechanisms that control biomineralization and organic-inorganic interactions in urolithiasis are reviewed, and the validity of drawing parallels between them are considered. The multifunctional nature of many biomineralization molecules, particularly proteins, and the different roles they play in and out of solution suggest avenues for resolving the apparent paradox that urinary macromolecules have been observed both to inhibit and to promote stone formation.

Identification of residues in the Gla-domain of human factor IX involved in the binding to conformation specific antibodies

The binding of Ca2+ induces a conformational change in factor IX which can be monitored with conformation specific antibodies. Anti-FIX:Mg(II) antibodies recognize a conformational epitope (FIX') that can be induced by several metal ions such as Ca2+, Mg2+, Mn2+ and Ba2+, while anti-FIX:Ca(II) antibodies recognize a conformational epitope (FIX*) that can be only induced by Ca2+ and Sr2+ ions (Liebman et al., J. Biol. Chem., vol. 262 (1987) pp. 7605-7612). The latter conformation is essential for the function of factor IX. In this study we tried to identify residues in the Gla-domain of factor IX which are involved in binding to anti-factor IX:Mg(II) and anti-factor IX:Ca(II) antibodies. For this we substituted residues in recombinant human factor IX for those of factor X or factor VII. The substitution of residues 1-40 of factor IX by those of factor VII eliminated binding to both types of antibodies. Re-introduction of factor IX specific residues increased the binding to conformation specific anti-factor IX antibodies, but reduced the binding to conformation specific anti-factor VII antibodies, indicating that the structural integrity of the Gla-domain was not seriously affected by the mutations. We provide evidence that residues 33, 39 and 40 of human factor IX are important for binding to anti-factor IX:Mg(II) antibodies, while residues 1-11 are important for binding to anti-factor IX:Ca(II) antibodies.

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

The solution structures of the N-terminal domains of protein S, a plasma vitamin K-dependent glycoprotein, and its homolog growth arrest specific protein 6 (Gas6) were predicted by molecular dynamics computer simulations. The initial structures were based on the x-ray crystallographic structure of the corresponding region of bovine prothrombin fragment 1. The subsequent molecular dynamics trajectories were calculated using the second-generation AMBER force field. The long-range electrostatic forces were evaluated by the particle mesh Ewald method. The structures that stabilized over a 400-ps time interval were compared with the corresponding region of the simulated solution structure of bovine prothrombin fragment 1. Structural properties of the gamma-carboxyglutamic acid (Gla) domains obtained from simulations and calcium binding were found to be conserved for all three proteins. Analysis of the predicted solution structure of the Gla domain of Gas6 suggests that this domain should bind with negatively charged phospholipid surfaces analogous to bovine prothrombin fragment 1 and protein S.

A theoretical model for the Gla-TSR-EGF-1 region of the anticoagulant cofactor protein S: from biostructural pathology to species-specific cofactor activity

Protein S (PS), which functions as a species-specific anticoagulant cofactor to activated protein C (APC), is a mosaic protein that interacts with the phospholipid membrane via its gamma-carboxyglutamate-rich (Gla) module. This module is followed by the thrombin-sensitive region (TSR), sensitive to thrombin cleavage, four epidermal growth factor (EGF)-like modules and a last region referred to as the sex hormone binding globulin (SHBG) domain. Of these, the TSR and the first EGF-like regions have been shown to be important for the species-specific interaction with APC. Difficulties in crystallising PS have so far hindered its study at the atomic level. Here, we report theoretical models for the Gla and EGF-1 modules of human PS constructed using prothrombin and factor X experimental structures. The TSR was built interactively. Analysis of the model linked with the large body of biochemical literature on PS and related proteins leads to suggestions that (i) the TSR stabilises the calcium-loaded Gla module through hydrophobic and ionic interactions and its conformation depends on the presence of the Gla module; (ii) the TSR does not form a calcium binding site but is protected from thrombin cleavage in the calcium-loaded form owing to short secondary structure elements and close contact with the Gla module; (iii) the PS missense mutations in this region are consistent with the structural data, except in one case which needs further investigation; and (iv) the two PS 'faces' involving regions of residues Arg49-Gln52-Lys97 (TSR-EGF-1) and Thr103-Pro106 (EGF-1) may be involved in species-specific interactions with APC as they are richer in nonconservative substitution when comparing human and bovine protein S. This preliminary model helps to plan future experiments and the resulting data will be used to further validate and optimise the present structure.