Factor Va alternative conformation reconstruction using atomic force microscopy.

Protein conformational variability (or dynamics) for large macromolecules and its implication for their biological function attracts more and more attention. Collective motions of domains increase the ability of a protein to bind to partner molecules. Using atomic force microscopy (AFM) topographic images, it is possible to take snapshots of large multi-component macromolecules at the single molecule level and to reconstruct complete molecular conformations. Here, we report the application of a reconstruction protocol, named AFM-assembly, to characterise the conformational variability of the two C domains of human coagulation factor Va (FVa). Using AFM topographic surfaces obtained in liquid environment, it is shown that the angle between C1 and C2 domains of FVa can vary between 40° and 166°. Such dynamical variation in C1 and C2 domain arrangement may have important implications regarding the binding of FVa to phospholipid membranes.

Genotype and phenotype relationships in 10 Pakistani unrelated patients with inherited factor VII deficiency.

Inherited factor VII (FVII) deficiency is one of the commonest rare bleeding disorders. It is characterized by a wide molecular and clinical heterogeneity and an autosomal recessive pattern of inheritance. Factor VII-deficient patients are still scarcely explored in Pakistan although rare bleeding disorders became quite common as a result of traditional consanguineous marriages. The aim of the study was to give a first insight of F7 gene mutations in Pakistani population. Ten unrelated FVII-deficient patients living in Pakistan were investigated (median FVII:C = 2%; range = 2-37%). A clinical questionnaire was filled out for each patient and direct sequencing was performed on the coding regions, intron/exon boundaries and 5' and 3' untranslated regions of the F7 gene. Nine different mutations (eight missense mutations and one located within the F7 promoter) were identified on the F7 gene. Five of them were novel (p.Cys82Tyr, p.Cys322Ser, p.Leu357Phe, p.Thr410Ala, c-57C>T, the last being predicted to alter the binding site of transcription factor HNF-4). Half of the patients had single mutations in Cys residues involved in disulfide bridges. The p.Cys82Arg mutation was the most frequent in our series. Six of seven patients with FVII:C levels below 10% were homozygous in connection with the high percentage of consanguinity in our series. In addition, we graded the 10 patients according to three previously published classifications for rare bleeding disorders. The use of the bleeding score proposed by Tosetto and co-workers in 2006 appears to well qualify the bleeding tendency in our series.

An integrated methodology for data processing in dynamic force spectroscopy of ligand-receptor binding.

Dynamic force spectroscopy (DFS), using atomic force microscopy (AFM), is a powerful tool to study ligand-receptor binding. The interaction mode of two binding partners is investigated by exploring stochastic behaviors of bond rupture events. However, to define a rupture event from force-distance measurements is not conclusive or unique in literature. To reveal the influence of event identification methods, we have developed an efficient protocol to manage tremendous amount of data by implementing different choices of peak selection from the force-distance curve. This data processing software simplifies routinely experimental procedures such as cantilever spring constant and force-distance curve calibrations, statistical treatments of data, and analysis distributions of rupture events. In the present work, we took available experimental data from a complex between a chelate metal compound and a monoclonal antibody as a study system.

Structural insights into protein-uranyl interaction: towards an in silico detection method.

Documenting the modes of interaction of uranyl (UO(2)2+) with large biomolecules, and particularly with proteins, is instrumental for the interpretation of its behavior in vitro and in vivo. The gathering of three-dimensional information concerning uranyl-first shell atoms from two structural databases, the Cambridge Structural Databank and the Protein Data Bank (PDB) allowed a screening of corresponding topologies in proteins of known structure. In the computer-aided procedure, all potentially bound residues from the template structure were granted full flexibility using a rotamer library. The Amber force-field was used to loosen constraints and score each predicted site. Our algorithm was validated as a first stage through the recognition of existing experimental data in the PDB. The coherent localization of missing atoms in the density map of an ambiguous uranium/uranyl-protein complex exemplified the efficiency of our approach, which is currently suggesting the experimental investigation of uranyl-protein binding site.

Intrinsic stability and functional properties of disulfide bond-stabilized coagulation factor VIIIa variants.

BACKGROUND: The utility of purified coagulation factor (F)VIII for treatment of hemophilia A is limited in part by its instability following activation by thrombin, which is caused by spontaneous dissociation of the A2 domain from the activated FVIII (FVIIIa) heterotrimer. To prevent this A2 domain dissociation in FVIIIa, we previously engineered a cysteine pair (C664-C1826) in recombinant FVIII that formed a disulfide bond cross-linking the A2 domain in the heavy chain to the A3 domain in the light chain. This engineered disulfide bond resulted in a more stable FVIIIa. AIMS: Here, we characterize the functional parameters of C664-C1828 FVIII and of a new disulfide bond-stabilized FVIII (C662-C1828 FVIII). METHODS: In order to assess whether these FVIII variants might be good candidates for a new therapeutic agent to treat hemophilia A, we investigated a variety of functional parameters that might affect the in vivo properties of the variants, including half-life of disulfide bond-stabilized FVIII and FVIIIa and the potency of these FVIIIa molecules in the FXase complex. RESULTS: Both disulfide bond-stabilized variants had improved affinity for von Willebrand factor (VWF). In studies of FX activation by purified FIXa and FVIIIa, C662-C1828 FVIIIa had normal activity while C664-C1826 FVIIIa had reduced activity. Both C664-C1826 FVIIIa and C662-C1828 FVIIIa were inactivated by activated protein C (APC) but the rates of inactivation were different. CONCLUSION: Overall, the specific location of the disulfide bridge between the A2 and A3 domains appears to affect functional properties of FVIIIa. In summary, introduction of engineered interdomain disulfides results in FVIIIa variants that resist spontaneous loss of activity while retaining susceptibility to APC proteolytic inactivation and maintaining VWF binding.

An engineered interdomain disulfide bond stabilizes human blood coagulation factor VIIIa.

The blood coagulation disorder, hemophilia A, is caused by deficiency of coagulation factor (F)VIII. Hemophilia A is now treated by infusions of pure FVIII, but the activity of FVIII is limited because it is unstable following activation by thrombin. This instability of activated FVIII is the result of dissociation of the A2 subunit. To obtain increased stability in FVIIIa, a disulfide bond between the A2 domain and the A3 domain, preventing A2 subunit dissociation, has been engineered. Structural analysis of the FVIII A domain homology model allowed us to identify residues 664 and 1826 as a potential disulfide bond pair. A FVIII mutant containing Cys664 and Cys1826 was produced and purified (C664-C1826 FVIII). Immunoblotting showed that a disulfide bond did form to link covalently the A2 and the A3 domains. Following activation of the recombinant C664-C1826 FVIII by thrombin, the mutant FVIIIa had increased stability and retained more than 90% of its clotting activity at a time at which wild-type FVIIIa lost more than 90% of its activity. This remarkably stable C664-C1826 FVIIIa provides a unique approach for studies of the cofactor activity of FVIIIa and also for new, improved therapy for hemophilia A.

F-box protein Grr1 interacts with phosphorylated targets via the cationic surface of its leucine-rich repeat.

The flexibility and specificity of ubiquitin-dependent proteolysis are mediated, in part, by the E3 ubiquitin ligases. One class of E3 enzymes, SKp1/cullin/F-box protein (SCF), derives its specificity from F-box proteins, a heterogeneous family of adapters for target protein recognition. Grr1, the F-box component of SCF(Grr1), mediates the interaction with phosphorylated forms of the G(1) cyclins Cln1 and Cln2. We show that binding of Cln2 by SCF(Grr1) was dependent upon its leucine-rich repeat (LRR) domain and its carboxy terminus. Our structural model for the Grr1 LRR predicted a high density of positive charge on the concave surface of the characteristic horseshoe structure. We hypothesized that specific basic residues on the predicted concave surface are important for recognition of phosphorylated Cln2. We show that point mutations that converted the basic residues on the concave surface but not those on the convex surface to neutral or acidic residues interfered with the capacity of Grr1 to bind to Cln2. The same mutations resulted in the stabilization of Cln2 and Gic2 and also in a spectrum of phenotypes characteristic of inactivation of GRR1, including hyperpolarization and enhancement of pseudohyphal growth. It was surprising that the same residues were not important for the role of Grr1 in nutrient-regulated transcription of HXT1 or AGP1. We concluded that the cationic nature of the concave surface of the Grr1 LRR is critical for the recognition of phosphorylated targets of SCF(Grr1) but that other properties of Grr1 are required for its other functions.

Three-dimensional model of coagulation factor Va bound to activated protein C.

A complete molecular model of blood coagulation factor Va (FVa) bound to anticoagulant activated protein C (APC) and to a phospholipid membrane was constructed. The three homologous A domains and the two homologous C domains of FVA were modeled based on the X-ray crystallographic structures of ceruloplasmin and C2 domain of factor V, respectively. The final arrangement of the five domains in the complete FVa model bound to a membrane incorporated extensive published experimental data. FVa binds the phospholipid membrane through its C2 domain while the A-domain trimer is located from 40 through 100 A above the membrane plane. From our model we infer a probable role for metal ions at the interface between FVa light and heavy chains, provide an explanation for the slower APC cleavage at Arg306 relative to Arg506, and predict specific interactions between positively and negatively charged exosites in APC and FVa, respectively.

Stabilization of bound polycyclic aromatic hydrocarbons by a pi-cation interaction.

Proteins can use aromatic side-chains to stabilize bound cationic ligands through cation-pi interactions. Here, we report the first example of the reciprocal process, termed pi-cation, in which a cationic protein side-chain stabilizes a neutral aromatic ligand. Site-directed mutagenesis revealed that an arginine side-chain located in the deep binding pocket of a monoclonal antibody (4D5) is essential for binding the neutral polynuclear aromatic hydrocarbon benzo[a]pyrene. This Arg was very likely selected for in the primary response, further underscoring the importance of the pi-cation interaction for ligand binding, which should be considered in protein analysis and design when ligands include aromatic groups.

Blood coagulation: The outstanding hydrophobic residues.

Newly determined crystal structures suggest that the membrane-binding C2 domains of blood coagulation cofactors Va and VIIIa bind anionic phospholipids through protruding solvent-exposed hydrophobic residues, aided by a crown of positively charged residues and by specific hydrogen-bonding side chains.

Structural basis for hemophilia A caused by mutations in the C domains of blood coagulation factor VIII.

Three dimensional homology models for the C1 and C2 domains of factor VIII (FVIII) were generated. Each C domain formed a beta-sandwich, and C1 was covalently connected to C2 in a head-to-head orientation. Of the >250 missense mutations that cause FVIII deficiency and hemophilia A, 34 are in the C domains. We used the FVIII C1-C2 model to infer the structural basis for the pathologic effects of these mutations. The mutated residues were divided into four categories: 15 conserved buried residues that affect normal packing of the hydrophobic side chains, 2 non-conserved buried residues that affect structure, 11 conserved exposed residues and 6 non-conserved exposed residues. The effects of all 34 missense mutations can be rationalized by predictable disruptions of FVIII structure while at most four mutations (S2069F, T2154I, R2209Q/G/L and E2181D) may affect residues directly involved in intermolecular interactions of FVIII/VIIIa with other coagulation factors or vWF.

Unraveling the effect of changes in conformation and compactness at the antibody V(L)-V(H) interface upon antigen binding.

We have analyzed conformational changes that occur at the interface between the light (V(L)) and heavy (V(H)) chains in antibody variable fragments upon binding to antigens. We wrote and applied the Tiny Probe program that computes the buried atomic contact surface area of three-dimensional structures to evaluate changes in compactness of the V(L)-V(H) interface between bound and unbound antibodies. We found three categories of these changes, which correlated with the size of the antigen. Upon binding, medium-sized nonprotein antigens cause an opening of the V(L)-V(H) interface (less compact), small antigens or haptens cause a closure of the interface (more compact), whereas large protein antigens have little effect on the compactness of the V(L)-V(H) interface. The largest changes in the atomic buried contact surface area at the V(L)-V(H) interface occur in residue pairs providing two 'shock absorbers' between the edge beta-strands of the V(L) and V(H) beta-sheets forming the antibody binding site. Importantly, the correlation between the size of antigens and conformational changes indicates that the V(L)-V(H) interface in antibodies plays a significant role in the antigen binding process. Furthermore, as the energy involved in such a motion is significant (up to 3 kcal/mol), these results provide a general mechanism for how residues distant from the combining site can significantly alter the affinity of an antibody for its antigen. Thus, mutations introduced at the V(L)-V(H) interface can be used to change antibody binding affinity with antigens. Due to the tightly packed V(L)-V(H) interface, the introduction of random mutations is not advisable. Rather our analysis suggests that concerted mutations of residues preceding CDRL2 and following CDRH3 or residues preceding CDRH2 and at the end of CDRL3 are most likely to alter or improve antigen binding affinity.

Biological sensors: More than one way to sense oxygen.

Recently determined structures of the oxygen-sensing heme domain of the bacterial protein FixL have revealed a new binding environment and signal transduction mechanism for heme; they have also provided new insights into the diverse 'PAS' domain superfamily.

Substrate specificity of prostate-specific antigen (PSA).

BACKGROUND: The serine protease prostate-specific antigen (PSA) is a useful clinical marker for prostatic malignancy. PSA is a member of the kallikrein subgroup of the (chymo)trypsin serine protease family, but differs from the prototypical member of this subgroup, tissue kallikrein, in possessing a specificity more similar to that of chymotrypsin than trypsin. We report the use of two strategies, substrate phage display and iterative optimization of natural cleavage sites, to identify labile sequences for PSA cleavage. RESULTS: Iterative optimization and substrate phage display converged on the amino-acid sequence SS(Y/F)Y decreases S(G/S) as preferred subsite occupancy for PSA. These sequences were cleaved by PSA with catalytic efficiencies as high as 2200-3100 M-1 s-1, compared with values of 2-46 M-1 s-1 for peptides containing likely physiological target sequences of PSA from the protein semenogelin. Substrate residues that bind to secondary (non-S1) subsites have a critical role in defining labile substrates and can even cause otherwise disfavored amino acids to bind in the primary specificity (S1) pocket. CONCLUSION: The importance of secondary subsites in defining both the specificity and efficiency of cleavage suggests that substrate recognition by PSA is mediated by an extended binding site. Elucidation of preferred subsite occupancy allowed refinement of the structural model of PSA and should facilitate the development of more sensitive activity-based assays and the design of potent inhibitors.

Photoactive yellow protein: a structural prototype for the three-dimensional fold of the PAS domain superfamily.

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the approximately 50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical approximately 150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of approximately 125-150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.

Homology models of the C domains of blood coagulation factors V and VIII: a proposed membrane binding mode for FV and FVIII C2 domains.

We present homology models of the C domains of coagulation factors V (FV) and VIII (FVIII). Using a threading approach, we identified the binding domain of galactose oxidase as an appropriate template for each C domain. The C1 and C2 domains of FV associate to form an elongated cylinder of 80A long and 30A diameter. The folding unit is a beta-sandwich with a long axis of 40A and a diameter of 30A. The current model allows us to propose a membrane binding mode for the C2 domains of FV and FVIII with three major characteristics: 1) solvent-exposed hydrophobic side chains from three loops at one end of the beta-sandwich are buried in the hydrophobic layer of the outer phospholipid leaflet; 2) a crown of positively charged residues is located in the polar zone of the phospholipid head groups; and 3) the long axis of the beta-sandwich of the C2 domain is perpendicular to the plane of the membrane. This proposal satisfies experimentally observed characteristics of membrane binding for the C2 domain and the light chain of FVa.

Does conformational free energy distinguish loop conformations in proteins?

Limitations in protein homology modeling often arise from the inability to adequately model loops. In this paper we focus on the selection of loop conformations. We present a complete computational treatment that allows the screening of loop conformations to identify those that best fit a molecular model. The stability of a loop in a protein is evaluated via computations of conformational free energies in solution, i.e., the free energy difference between the reference structure and the modeled one. A thermodynamic cycle is used for calculation of the conformational free energy, in which the total free energy of the reference state (i.e., gas phase) is the CHARMm potential energy. The electrostatic contribution of the solvation free energy is obtained from solving the finite-difference Poisson-Boltzmann equation. The nonpolar contribution is based on a surface area-based expression. We applied this computational scheme to a simple but well-characterized system, the antibody hypervariable loop (complementarity-determining region, CDR). Instead of creating loop conformations, we generated a database of loops extracted from high-resolution crystal structures of proteins, which display geometrical similarities with antibody CDRs. We inserted loops from our database into a framework of an antibody; then we calculated the conformational free energies of each loop. Results show that we successfully identified loops with a "reference-like" CDR geometry, with the lowest conformational free energy in gas phase only. Surprisingly, the solvation energy term plays a confusing role, sometimes discriminating "reference-like" CDR geometry and many times allowing "non-reference-like" conformations to have the lowest conformational free energies (for short loops). Most "reference-like" loop conformations are separated from others by a gap in the gas phase conformational free energy scale. Naturally, loops from antibody molecules are found to be the best models for long CDRs (> or = 6 residues), mainly because of a better packing of backbone atoms into the framework of the antibody model.

Predicting antigenic determinants in proteins: looking for unidimensional solutions to a three-dimensional problem?

In a recent review, Hopp (Peptide Research 6:183-190, 1993) claimed that the Hopp and Woods hydrophilicity method for locating antigenic determinants is superior to all other existing methods for predicting the B cell epitopes of proteins but that it is not useful to aid the investigator in producing peptide-protein cross-reactive antisera. In this article, we challenge both these assertions. Most investigators utilize antigenicity prediction algorithms because they wish to produce anti-peptide antibodies capable of cross-reacting with the intact protein. All prediction methods are based on propensity scales for the 20 amino acids, which describe the tendency of each residue to be associated with properties such as hydrophilicity, surface accessibility or segmental mobility. When we compared the prediction efficacy of 22 different scales, taking into account both correct and incorrect predictions, we found that none of the scales gave a level of correct prediction higher than about 50%-60%. If no antigenicity was found in a particular region of the protein, we took the view that hydrophilicity peaks located in that region amounted to wrong predictions. The much higher success rate reported by Hopp for this method stems from the way he assesses prediction efficacy, i.e., by counting the number of known epitopes located inside and outside hydrophilicity peaks. Reasons for the low success rate of antigenicity prediction are discussed. In most cases, it is unrealistic to try to reduce the complexity of discontinuous, conformational epitopes to simple, linear peptide models.

The Wilson disease gene is a copper transporting ATPase with homology to the Menkes disease gene.

Wilson disease (WD) is an autosomal recessive disorder characterized by the toxic accumulation of copper in a number of organs, particularly the liver and brain. As shown in the accompanying paper, linkage disequilibrium & haplotype analysis confirmed the disease locus to a single marker interval at 13q14.3. Here we describe a partial cDNA clone (pWD) which maps to this region and shows a particular 76% amino acid homology to the Menkes disease gene, Mc1. The predicted functional properties of the pWD gene together with its strong homology to Mc1, genetic mapping data and identification of four independent disease-specific mutations, provide convincing evidence that pWD is the Wilson disease gene.