Large Hydrophobic Side Chains Research Articles

An internal docking site stabilizes substrate binding to γ-secretase: Analysis by molecular dynamics simulations

Amyloid precursor protein (APP) is cleaved and processed sequentially by γ-secretase yielding amyloid β (Aβ) peptides of different lengths. Longer Aβ peptides are associated with the formation of neurotoxic plaques related to Alzheimer’s disease. Based on the APP substrate-bound structure of γ-secretase, we investigated the enzyme-substrate interaction using molecular dynamics simulations and generated model structures that represent the sequentially cleaved intermediates during the processing reaction. The simulations indicated an internal docking site providing strong enzyme-substrate packing interaction. In the enzyme-substrate complex, it is located close to the region where the helical conformation of the substrate is interrupted and continues toward the active site in an extended conformation. The internal docking site consists of two non-polar pockets that are preferentially filled by large hydrophobic or aromatic substrate side chains to stabilize binding. Placement of smaller residues such as glycine can trigger a shift in the cleavage pattern during the simulations or results in destabilization of substrate binding. The reduced packing by smaller residues also influences the hydration of the active site and the formation of a catalytically active state. The simulations on processed substrate intermediates and a substrate G33I mutation offer an explanation of the experimentally observed relative increase of short Aβ fragment production for this mutation. In addition, studies on a substrate K28A mutation indicate that the internal docking site opposes the tendency of substrate dissociation due to a hydrophobic mismatch at the membrane boundary caused by K28 during processing and substrate movement toward the enzyme active site. The proposed internal docking site could also be useful for the specific design of new γ-secretase modulators.

Mutation of the eunicellane synthase Bnd4 alters its product profile and expands its prenylation ability.

Bnd4 catalyzes the first committed step in the biosynthesis of the bacterial diterpenoid benditerpenoic acid and was the first eunicellane synthase identified from nature. We investigated the catalytic roles of the aromatic residues in the active site of Bnd4 through a series of mutation studies. These experiments revealed that large hydrophobic or aromatic side chains are required at F162 and Y197 for eunicellane formation and that selected mutations at W316 converted Bnd4 into a cembrane synthase. In addition, the Bnd4Y197A variant expanded the native prenylation ability of Bnd4 from accepting C5 and C10 prenyl donors to C15. This study supports the mechanism of eunicellane formation by Bnd4 and encourages further engineering of terpene synthases into practical and efficient prenyltransferases.

Collagen’s primary structure determines collagen:HSP47 complex stoichiometry

Collagens play important roles in development and homeostasis in most higher organisms. In order to function, collagens require the specific chaperone HSP47 for proper folding and secretion. HSP47 is known to bind to the collagen triple helix, but the exact positions and numbers of binding sites are not clear. Here, we employed a collagen II peptide library to characterize high-affinity binding sites for HSP47. We show that many previously predicted binding sites have very low affinities due to the presence of a negatively charged amino acid in the binding motif. In contrast, large hydrophobic amino acids such as phenylalanine at certain positions in the collagen sequence increase binding strength. For further characterization, we determined two crystal structures of HSP47 bound to peptides containing phenylalanine or leucine. These structures deviate significantly from previously published ones in which different collagen sequences were used. They reveal local conformational rearrangements of HSP47 at the binding site to accommodate the large hydrophobic side chain from the middle strand of the collagen triple helix and, most surprisingly, possess an altered binding stoichiometry in the form of a 1:1 complex. This altered stoichiometry is explained by steric collisions with the second HSP47 molecule present in all structures determined thus far caused by the newly introduced large hydrophobic residue placed on the trailing strand. This exemplifies the importance of considering all three sites of homotrimeric collagen as independent interaction surfaces and may provide insight into the formation of higher oligomeric complexes at promiscuous collagen-binding sites.

Structural basis of different surface-modified fullerene derivatives as novel thrombin inhibitors: insight into the inhibitory mechanism through molecular modelling studies

Thrombin is the main target of antithrombotic and anticoagulant therapy. Although recent studies have shown that surface-modified fullerene derivatives can be used as novel nano-anticoagulants to effectively inhibit thrombin, the specific mechanism remains unclear. In this study, we investigated the binding mode and inhibition mechanism of fullerenes with thrombin by using various molecular modelling approaches. The results showed that owing to the hydrophobicity of the fullerene skeleton, the four existing fullerene inhibitors mainly combined with the active site of thrombin through van der Waals interaction and tended to make the surrounding environment hydrophobic. This made the inhibitors with hydrophobic side chains, especially those with large hydrophobic side chains containing benzene rings, have greater binding free energy with protein. However, as more hydrophilic groups can change the surrounding hydrophilic–hydrophobic environment and form more hydrogen bonds, the electrostatic contribution of the polyhydroxy side chain designed in this study was increased further. At the same time, the van der Waals interaction between the fullerene carbon skeleton and protein did not decrease significantly, resulting in a significantly increased total binding free energy and an enhanced inhibition effect. These findings provide a theoretical basis for the improvement of highly effective fullerene anticoagulants.

Structural basis for ion selectivity in TMEM175 K+ channels.

The TMEM175 family constitutes recently discovered K+channels that are important for autophagosome turnover and lysosomal pH regulation and are associated with the early onset of Parkinson Disease. TMEM175 channels lack a P-loop selectivity filter, a hallmark of all known K+ channels, raising the question how selectivity is achieved. Here, we report the X-ray structure of a closed bacterial TMEM175 channel in complex with a nanobody fusion-protein disclosing bound K+ ions. Our analysis revealed that a highly conserved layer of threonine residues in the pore conveys a basal K+ selectivity. An additional layer comprising two serines in human TMEM175 increases selectivity further and renders this channel sensitive to 4-aminopyridine and Zn2+. Our findings suggest that large hydrophobic side chains occlude the pore, forming a physical gate, and that channel opening by iris-like motions simultaneously relocates the gate and exposes the otherwise concealed selectivity filter to the pore lumen.

Separation and Quantification of Superwarfarin Rodenticide Diastereomers-Bromadiolone, Difenacoum, Flocoumafen, Brodifacoum, and Difethialone-in Human Plasma.

Superwarfarins, second-generation long-acting anticoagulant rodenticides, are 4-hydroxycoumarin analogues of warfarin that contain a large hydrophobic side chain. These compounds contain two chiral centers and are synthesized for commercial use as two pairs of diastereomer. To support studies of superwarfarin pharmacokinetics and other efforts to improve clinical care for poisoning victims, a quantitative assay was developed for the measurement of diastereomer of bromadiolone, difenacoum, flocoumafen, brodifacoum, and difethialone in human plasma. Based on ultrahigh-pressure liquid chromatography-tandem mass-spectrometry (UHPLC-MS/MS), this method was validated according to U.S. Food and Drug Administration (FDA) guidelines. Sample preparation involved simple protein precipitation followed by reversed phase UHPLC, which resolved all five pairs of cis/trans diastereomer in less than 10 min. Superwarfarins were measured using negative ion electrospray followed by selected-reaction monitoring on a triple quadrupole mass spectrometer. Calibration curves covered 3-4 orders of magnitude with linear regression coefficients of >0.999. The lower limits of quantitation were from 0.013 to 2.41 ng/mL, and intra-day and inter-day accuracy and precision coefficients of variation were <12%. A 10-min UHPLC-MS/MS assay was developed and validated for the separation and quantitative analysis of the pairs of diastereomer of five superwarfarins in human plasma. This method was used to identify and measure superwarfarins and their cis/trans diastereomers in plasma obtained from patients treated for coagulopathy following consumption of contaminated synthetic cannabinoid products.

Interactions between a subset of substrate side chains and AAA+ motor pore loops determine grip during protein unfolding.

Most AAA+ remodeling motors denature proteins by pulling on the peptide termini of folded substrates, but it is not well-understood how motors produce grip when resisting a folded domain. Here, at single amino-acid resolution, we identify the determinants of grip by measuring how substrate tail sequences alter the unfolding activity of the unfoldase-protease ClpXP. The seven amino acids abutting a stable substrate domain are key, with residues 2-6 forming a core that contributes most significantly to grip. ClpX grips large hydrophobic and aromatic side chains strongly and small, polar, or charged side chains weakly. Multiple side chains interact with pore loops synergistically to strengthen grip. In combination with recent structures, our results support a mechanism in which unfolding grip is primarily mediated by non-specific van der Waal's interactions between core side chains of the substrate tail and a subset of YVG loops at the top of the ClpX axial pore.

Crystallographic and thermodynamic characterization of phenylaminopyridine bisphosphonates binding to human farnesyl pyrophosphate synthase.

Human farnesyl pyrophosphate synthase (hFPPS) catalyzes the production of the 15-carbon isoprenoid farnesyl pyrophosphate. The enzyme is a key regulator of the mevalonate pathway and a well-established drug target. Notably, it was elucidated as the molecular target of nitrogen-containing bisphosphonates, a class of drugs that have been widely successful against bone resorption disorders. More recently, research has focused on the anticancer effects of these inhibitors. In order to achieve increased non-skeletal tissue exposure, we created phenylaminopyridine bisphosphonates (PNP-BPs) that have bulky hydrophobic side chains through a structure-based approach. Some of these compounds have proven to be more potent than the current clinical drugs in a number of antiproliferation assays using multiple myeloma cell lines. In the present work, we characterized the binding of our most potent PNP-BPs to the target enzyme, hFPPS. Co-crystal structures demonstrate that the molecular interactions designed to elicit tighter binding are indeed established. We carried out thermodynamic studies as well; the newly introduced protein-ligand interactions are clearly reflected in the enthalpy of binding measured, which is more favorable for the new PNP-BPs than for the lead compound. These studies also indicate that the affinity of the PNP-BPs to hFPPS is comparable to that of the current drug risedronate. Risedronate forms additional polar interactions via its hydroxyl functional group and thus exhibits more favorable binding enthalpy; however, the entropy of binding is more favorable for the PNP-BPs, owing to the greater desolvation effects resulting from their large hydrophobic side chains. These results therefore confirm the overall validity of our drug design strategy. With a distinctly different molecular scaffold, the PNP-BPs described in this report represent an interesting new group of future drug candidates. Further investigation should follow to characterize the tissue distribution profile and assess the potential clinical benefits of these compounds.

Roles of Leu28 side chain intercalation in the interaction between Cren7 and DNA.

Crenarchaeal chromatin protein Cren7 binds double-stranded DNA in the minor groove, introducing a sharp single-step DNA kink. The side chain of Leu28, a residue conserved among all Cren7 homologs, intercalates into the kinked DNA step. In the present study, we replaced Leu28 with a residue containing a hydrophobic side chain of different sizes (i.e. L28A, L28V, L28I, L28M and L28F). Both the stability of the Cren7-DNA complex and the ability of Cren7 to constrain DNA supercoils correlated well with the size of the intercalated side chain. Structural analysis shows that L28A induces a kink (∼43°), nearly as sharp as that produced by wild-type Cren7 (∼48°), in the bound DNA fragment despite the lack of side chain intercalation. In another duplex DNA fragment, L28F inserts a large hydrophobic side chain deep into the DNA step, but introduces a smaller kink (∼39°) than that formed by the wild-type protein (∼50°). Mutation of Leu28 into methionine yields two protein conformers differing in loop β3-β4 orientation, DNA-binding surface and DNA geometry in the protein-DNA structure. Our results indicate that side chain intercalation is not directly responsible for DNA kinking or bending by Cren7, but plays a critical role in the stabilization of the Cren7-DNA complex. In addition, the flexibility of loop β3-β4 in Cren7, as revealed in the crystal structure of L28M-DNA, may serve a role in the modulation of chromosomal organization and function in the cell.

Missense Mutation R338W in ARHGEF9 in a Family with X-linked Intellectual Disability with Variable Macrocephaly and Macro-Orchidism.

Non-syndromal X-linked intellectual disability (NS-XLID) represents a broad group of clinical disorders in which ID is the only clinically consistent manifestation. Although in many cases either chromosomal linkage data or knowledge of the >100 existing XLID genes has assisted mutation discovery, the underlying cause of disease remains unresolved in many families. We report the resolution of a large family (K8010) with NS-XLID, with variable macrocephaly and macro-orchidism. Although a previous linkage study had mapped the locus to Xq12-q21, this region contained too many candidate genes to be analyzed using conventional approaches. However, X-chromosome exome sequencing, bioinformatics analysis and segregation analysis revealed a novel missense mutation (c.1012C>T; p.R338W) in ARHGEF9. This gene encodes collybistin (CB), a neuronal GDP-GTP exchange factor previously implicated in several cases of XLID, as well as clustering of gephyrin and GABAA receptors at inhibitory synapses. Molecular modeling of the CB R338W substitution revealed that this change results in the substitution of a long electropositive side-chain with a large non-charged hydrophobic side-chain. The R338W change is predicted to result in clashes with adjacent amino acids (K363 and N335) and disruption of electrostatic potential and local folding of the PH domain, which is known to bind phosphatidylinositol-3-phosphate (PI3P/PtdIns-3-P). Consistent with this finding, functional assays revealed that recombinant CB CB2SH3−R338W was deficient in PI3P binding and was not able to translocate EGFP-gephyrin to submembrane microaggregates in an in vitro clustering assay. Taken together, these results suggest that the R338W mutation in ARHGEF9 is the underlying cause of NS-XLID in this family.

The identification of novel Mycobacterium tuberculosis DHFR inhibitors and the investigation of their binding preferences by using molecular modelling.

It is an urgent need to develop new drugs for Mycobacterium tuberculosis (Mtb), and the enzyme, dihydrofolate reductase (DHFR) is a recognised drug target. The crystal structures of methotrexate binding to mt- and h-DHFR separately indicate that the glycerol (GOL) binding site is likely to be critical for the function of mt-DHFR selective inhibitors. We have used in silico methods to screen NCI small molecule database and a group of related compounds were obtained that inhibit mt-DHFR activity and showed bactericidal effects against a test Mtb strain. The binding poses were then analysed and the influence of GOL binding site was studied by using molecular modelling. By comparing the chemical structures, 4 compounds that might be able to occupy the GOL binding site were identified. However, these compounds contain large hydrophobic side chains. As the GOL binding site is more hydrophilic, molecular modelling indicated that these compounds were failed to occupy the GOL site. The most potent inhibitor (compound 6) demonstrated limited selectivity for mt-DHFR, but did contain a novel central core (7H-pyrrolo[3,2-f]quinazoline-1,3-diamine), which may significantly expand the chemical space of novel mt-DHFR inhibitors. Collectively, these observations will inform future medicinal chemistry efforts to improve the selectivity of compounds against mt-DHFR.

Differential roles of tryptophan residues in conformational stability of Porphyromonas gingivalis HmuY hemophore

BackgroundWe have previously shown that the P. gingivalis HmuY hemophore-like protein binds heme and scavenges heme from host hemoproteins to further deliver it to the cognate heme receptor HmuR. The aim of this study was to characterize structural features of HmuY variants in the presence and absence of heme with respect to roles of tryptophan residues in conformational stability.ResultsHmuY possesses tryptophan residues at positions 51 and 73, which are conserved in HmuY homologs present in a variety of bacteria, and a tryptophan residue at position 161, which has been found only in HmuY identified in P. gingivalis strains. We expressed and purified the wildtype HmuY and its protein variants with single tryptophan residues replaced by alanine or tyrosine residues. All HmuY variants were subjected to thermal denaturation and fluorescence spectroscopy analyses. Replacement of the most buried W161 only moderately affects protein stability. The most profound effect of the lack of a large hydrophobic side chain in respect to thermal stability is observed for W73. Also replacement of the W51 exposed on the surface results in the greatest loss of protein stability and even the large aromatic side chain of a tyrosine residue has little potential to substitute this tryptophan residue. Heme binding leads to different exposure of the tryptophan residue at position 51 to the surface of the protein. Differences in structural stability of HmuY variants suggest the change of the tertiary structure of the protein upon heme binding.ConclusionsHere we demonstrate differential roles of tryptophan residues in the protein conformational stability. We also propose different conformations of apo- and holoHmuY caused by tertiary changes which allow heme binding to the protein.

Roles of Hydrophobic Interactions and Hydrogen Bonds in Beta-Sheet Formation

In this presentation, we will discuss interactions of extended conformations of homodimeric peptides made of small (glycine or alanine) and large hydrophobic (valine or leucine) sidechains using all-atom computer simulations to decipher driving forces behind beta-sheet formation. Dimers adopt beta-sheet conformations at short inter-peptide distances (x ∼ 0.5 nm) while at intermediate distances (∼ 0.8 nm), dimers made of valine or leucine assume “cross-beta like” conformations with sidechains interpenetrating each other. These two states are identified as minima in the Potential of Mean Force (PMF). While the number of inter-peptide hydrogen bonds increases with decreasing inter-peptide distance, the total hydrogen bond number in the system does not change significantly, suggesting that formation of beta-sheet structures from extended conformations is not driven by hydrogen bonds. This is confirmed by an increase in electrostatic energy at short inter-peptide distances. A remarkable correlation between the volume of the system and the total electrostatic energy is observed, supporting the view that excluded water regions in proteins have an enthalpic penalty. We will also discuss microscopic mechanisms accounting for beta-sheet formation based on computed enthalpy and entropy and we will show that they are different for peptides with small and large sidechains.

Hydrophobic interactions and hydrogen bonds in β-sheet formation

In this study, we investigate interactions of extended conformations of homodimeric peptides made of small (glycine or alanine) and large hydrophobic (valine or leucine) sidechains using all-atom molecular dynamics simulations to decipher driving forces for β-sheet formation. We make use of a periodic boundary condition setup in which individual peptides are infinitely long and stretched. Dimers adopt β-sheet conformations at short interpeptide distances (ξ ~ 0.5 nm) and at intermediate distances (~0.8 nm), valine and leucine homodimers assume cross-β-like conformations with side chains interpenetrating each other. These two states are identified as minima in the potential of mean force. While the number of interpeptide hydrogen bonds increases with decreasing interpeptide distance, the total hydrogen bond number in the system does not change significantly, suggesting that formation of β-sheet structures from extended conformations is not driven by hydrogen bonds. This is supported by an increase in electrostatic energy at short interpeptide distances. A remarkable correlation between the volume of the system and the total electrostatic energy is observed, further reinforcing the idea that excluding water in proteins comes with an enthalpic penalty. We also discuss microscopic mechanisms accounting for β-sheet formation based on computed enthalpy and entropy and we show that they are different for peptides with small and large side chains.

Insight into structure‐function relationship in phenol‐soluble modulins using an alanine screen of the phenol‐soluble modulin (PSM) α3 peptide

Phenol-soluble modulins (PSMs) are a family of peptides with multiple functions in staphylococcal pathogenesis. To gain insight into the structural features affecting PSM functions, we analyzed an alanine substitution library of PSMα3, a strongly cytolytic and proinflammatory PSM of Staphylococcus aureus with a significant contribution to S. aureus virulence. Lysine residues were essential for both receptor-dependent proinflammatory and receptor-independent cytolytic activities. Both phenotypes also required additional structural features, with the C terminus being crucial for receptor activation. Biofilm formation was affected mostly by hydrophobic amino acid positions, suggesting that the capacity to disrupt hydrophobic interactions is responsible for the effect of PSMs on biofilm structure. Antimicrobial activity, absent from natural PSMα3, could be created by the exchange of large hydrophobic side chains, indicating that PSMα3 has evolved to exhibit cytolytic rather than antimicrobial activity. In addition to gaining insight into the structure-function relationship in PSMs, our study identifies nontoxic PSMα3 derivatives for active vaccination strategies and lays the foundation for future efforts aimed to understand the biological role of PSM recognition by innate host defense.

Properties and applications of phytepsins from thistle flowers

Aqueous extracts of thistle flowers from the genus Cynara-Cardueae tribe Cass. (Cynareae Less.), Asteraceae Dumortier-are traditionally used in the Mediterranean region for production of artisanal cheeses. This is because of the presence of aspartic proteases (APs) with the ability to coagulate milk. Plant APs, collectively known as phytepsins (EC 3.4.23.40), are bilobed endopeptidases present in an ample variety of plant species with activity mainly at acidic pHs, and have two aspartic residues located on each side of a catalytic cleft that are responsible for catalysis. The cleavage of the scissile peptide-bond occurs primarily between residues with large hydrophobic side-chains. Even when aspartylendopeptidase activity in plants is normally present at relatively low levels overall, the flowers of several species of the Cardueae tribe possess APs with extremely high specific activities in certain tissues. For this reason, in the last two decades, APs present in thistle flowers have been the subject of intensive study. Present here is a compilation of work that summarizes the known chemical and biological properties of these proteases, as well as their biomedical and biotechnological applications.

The α-defensin salt-bridge induces backbone stability to facilitate folding and confer proteolytic resistance

Salt-bridge interactions between acidic and basic amino acids contribute to the structural stability of proteins and to protein-protein interactions. A conserved salt-bridge is a canonical feature of the α-defensin antimicrobial peptide family, but the role of this common structural element has not been fully elucidated. We have investigated mouse Paneth cell α-defensincryptdin-4 (Crp4) and peptide variants with mutations at Arg7 or Glu15 residue positions to disrupt the salt-bridge and assess the consequences on Crp4 structure, function, and stability. NMR analyses showed that both (R7G)-Crp4 and (E15G)-Crp4 adopt native-like structures, evidence of fold plasticity that allows peptides to reshuffle side chains and stabilize the structure in the absence of the salt-bridge. In contrast, introduction of a large hydrophobic side chain at position 15, as in (E15L)-Crp4 cannot be accommodated in the context of the Crp4 primary structure. Regardless of which side of the salt-bridge was mutated, salt-bridge variants retained bactericidal peptide activity with differential microbicidal effects against certain bacterial cell targets, confirming that the salt-bridge does not determine bactericidal activity per se. The increased structural flexibility induced by salt-bridge disruption enhanced peptide sensitivity to proteolysis. Although sensitivity to proteolysis by MMP7 was unaffected by most Arg(7) and Glu(150 substitutions, every salt-bridge variant was degraded extensively by trypsin. Moreover, the salt-bridge facilitates adoption of the characteristic α-defensin fold as shown by the impaired in vitro refolding of (E15D)-proCrp4, the most conservative salt-bridge disrupting replacement. In Crp4, therefore, the canonical α-defensin salt-bridge facilitates adoption of the characteristic α-defensin fold, which decreases structural flexibility and confers resistance todegradation by proteinases.

Cavities in the Hydrophobic Core Govern Pressure Unfolding of Proteins

Since Bridgman's seminal experiments on high-pressure denaturation of albumen in 19141, the origin of pressure unfolding of proteins has remained unresolved. We report here a systematic study of the contribution of cavities to the volume difference between unfolded and folded states (ΔVu), using 10 single point variants of staphylococcus nuclease (SNase). Each mutation is localised in a strategic position and was designed to change a large buried hydrophobic side chain into alanine, thus opening tunable cavities in the SNase structure. For every variant, a crystal structure confirmed the presence of the designed cavity with no detectable presence of water molecules. High-pressure fluorescence experiments show significantly larger ΔVu values for the cavity mutants in comparison to the reference protein. This demonstrates that solvent-excluded cavities make a major contribution to ΔVu. Thus, pressure effects have their origin in a property of the folded states of proteins, unlike temperature and chemical denaturants, whose effects are governed by exposed surface area in the unfolded state. High-pressure NMR experiments on 4 cavity mutants, recording HSQCs peak intensities up to 2500 bar, allowed precise estimations of the apparent ΔVu monitored by more than two-thirds of the residues. An innovative combination of the site-specific NMR data and Go-model simulations revealed significant departures from the apparent two-state folding process for the SNase reference protein and the cavity mutants. This study opens up highly promising perspectives on the use of high pressure for characterization of folding landscapes inaccessible by other methods.1. Bridgman, P. W. J. Biol. Chem. 1914.19, 511–512.

Architecture of a Full-length Retroviral Integrase Monomer and Dimer, Revealed by Small Angle X-ray Scattering and Chemical Cross-linking

We determined the size and shape of full-length avian sarcoma virus (ASV) integrase (IN) monomers and dimers in solution using small angle x-ray scattering. The low resolution data obtained establish constraints for the relative arrangements of the three component domains in both forms. Domain organization within the small angle x-ray envelopes was determined by combining available atomic resolution data for individual domains with results from cross-linking coupled with mass spectrometry. The full-length dimer architecture so revealed is unequivocally different from that proposed from x-ray crystallographic analyses of two-domain fragments, in which interactions between the catalytic core domains play a prominent role. Core-core interactions are detected only in cross-linked IN tetramers and are required for concerted integration. The solution dimer is stabilized by C-terminal domain (CTD-CTD) interactions and by interactions of the N-terminal domain in one subunit with the core and CTD in the second subunit. These results suggest a pathway for formation of functional IN-DNA complexes that has not previously been considered and possible strategies for preventing such assembly.

Identification of Interacting Hot Spots in the β3 Integrin Stalk Using Comprehensive Interface Design

Protein-protein interfaces are usually large and complementary surfaces, but specific side chains, representing energetic "hot spots," often contribute disproportionately to binding free energy. We used a computational method, comprehensive interface design, to identify hot spots in the interface between the stalk regions of the β3 and the complementary αIIb and αv integrin subunits. Using the Rosetta alanine-scanning and design algorithms to predict destabilizing, stabilizing, and neutral mutations in the β3 region extending from residues Lys(532) through Gly(690), we predicted eight alanine mutations that would destabilize the αIIbβ3 interface as well as nine predicted to destabilize the αvβ3 interface, by at least 0.3 kcal/mol. The mutations were widely and unevenly distributed, with four between residues 552 and 563 and five between 590 and 610, but none between 565 and 589, and 611 and 655. Further, mutations destabilizing the αvβ3 and αIIbβ3 interfaces were not identical. The predictions were then tested by introducing selected mutations into the full-length integrins expressed in Chinese hamster ovary cells. Five mutations predicted to destabilize αIIb and β3 caused fibrinogen binding to αIIbβ3, whereas three of four predicted to be neutral or stabilizing did not. Conversely, a mutation predicted to destabilize αvβ3, but not αIIbβ3 (D552A), caused osteopontin binding to αvβ3, but not fibrinogen binding to αIIbβ3. These results indicate that stability of the distal stalk interface is involved in constraining integrins in stable, inactive conformations. Further, they demonstrate the ability of comprehensive interface design to identify functionally significant integrin mutations.