Modeling and molecular dynamics simulations of the V33 variant of the integrin subunit β3: Structural comparison with the L33 (HPA-1a) and P33 (HPA-1b) variants

Abstract

The human plateletalloantigen (HPA)-1 system, the first cause of alloimmune thrombocytopenia in Caucasians, results from leucine-to-proline substitution (alleles 1a and 1b) of residue 33 in β3 subunit of the integrin αIIbβ3. A third variant with a valine (V33) has been described. Although leucine and valine share similar physicochemical properties, sera containing alloantibodies to the HPA-1a antigen variably reacted with V33-β3, suggesting structural alterations of β3. To analyze the effect of the L33V transition, molecular dynamics simulations were performed on a 3D structural model of the V33 form of the whole β3 extracellular domain (690 residues). Dynamics of the PSI (carrying residue 33), I-EGF-1, and I-EGF-2 domains of β3 were compared to previously obtained dynamics of HPA-1a structure and HPA-1b structural model using classical and innovative developments (a structural alphabet). Clustering approach and local structure analysis showed that L33-β3 and V33-β3 mostly share common structures co-existing in different dynamic equilibria. The L33V substitution mainly displaces the equilibrium between common structures. These observations can explain the variable reactivity of anti-HPA-1a alloantibodies suggesting that molecular dynamic plays a key role in the binding of these alloantibodies. Unlike the L33P substitution, the L33V transition would not affect the structure flexibility of the β3 knee, and consequently the functions of αIIbβ3.