Single, Double and Quadruple Alanine Substitutions at Oligomeric Interfaces Identify Hydrophobicity as the Key Determinant of Human Neutrophil Alpha Defensin HNP1 Function

HNP1 is a human alpha defensin that forms dimers and multimers governed by hydrophobic residues, including Tyr16, Ile20, Leu25, and Phe28. Previously, alanine scanning mutagenesis identified each of these residues and other hydrophobic residues as important for function. Here we report further structural and functional studies of residues shown to interact with one another across oligomeric interfaces: I20A-HNP1 and L25A-HNP1, plus the double alanine mutants I20A/L25A-HNP1 and Y16A/F28A-HNP1, and the quadruple alanine mutant Y16A/I20A/L25A/F28A-HNP1. We tested binding to HIV-1 gp120 and HNP1 by surface plasmon resonance, binding to HIV-1 gp41 and HNP1 by fluorescence polarization, inhibition of anthrax lethal factor, and antibacterial activity using the virtual colony count assay. Similar to the previously described single mutant W26A-HNP1, the quadruple mutant displayed the least activity in all functional assays, followed by the double mutant Y16A/F28A-HNP1. The effects of the L25A and I20A single mutations were milder than the double mutant I20A/L25A-HNP1. Crystallographic studies confirmed the correct folding and disulfide pairing, and depicted an array of dimeric and tetrameric structures. These results indicate that side chain hydrophobicity is the critical factor that determines activity at these positions.


Introduction
Defensins are a family of antimicrobial peptides of innate immunity broadly active against bacteria, viruses and toxins [1,2,3,4,5,6]. They are small (2)(3)(4)(5), stabilized structurally by three disulfide bonds, form oligomers, and are divided into a [6], b [7] and h [8] structural classes. There are six human a-defensins: human neutrophil peptides (HNP) 1-4 [9,10,11,12,13] and human defensins (HD) 5-6 [14,15]. The structural basis for defensin antimicrobial activity against bacteria and viruses is only partially understood, but it is clear that defensins possess two fundamental characteristics: cationicity and hydrophobicity. a-Defensins are positively charged, and their activity against Gram negative bacteria such as Escherichia coli is thought to be mediated by the electrostatic attraction between cationic arginine residues and the anionic phospholipids of bacterial membranes [16]. By contrast, the human a-defensins are more selective for Gram positive strains, and specific, chiral interactions with lipid II molecules have been proposed to be primarily responsible for activity against Gram positive bacterium, Staphylococcus aureus [17], as is the case with several other defensins [18,19,20,21]. Defensins provide two types of protection against bacilli: they are active against Bacillus cereus directly with a high degree of potency [22,23], and they bind to and inhibit bacterial toxins such as anthrax lethal factor (LF) [24,25,26,27,28]. A picture of the myriad functions of HNP1 and HD5 is emerging that identifies hydrophobicity as more critical than cationicity, which can be explained partially by the hydrophobic effect on the a-defensin ''canonical'' dimer. For HNP1, the canonical dimer interface is formed by the antiparallel extensions of the b2 strands stabilized by the reciprocal mainchain hydrogen bonds contributed by Thr 18 and Ile 20 , and the hydrophobic packing of the side chains of Tyr 16 , Tyr 21 , Phe 28 and the Cys 2 -Cys 30 disulfide [27,29]. Underlying these hydrophobic contacts are the aromatic rings of Trp 26 [28], constituting the bulk of the hydrophobic core of the monomer.
Alanine scanning mutagenesis is a powerful tool for exploring the contributions of individual side chains to aspects of protein and peptide function such as binding [30,31,32], stability [33,34,35], and catalysis [36,37]. Alanine scanning mutageneses of both HNP1 [28] and HD5 [26] identified hydrophobic residues as the most important for activity, especially HNP1 Trp 26 and HD5 Leu 29 . In the alanine scan of HNP1, the hydrophobic residues of Tyr 16 , Ile 20 , Leu 25 and Phe 28 were also shown to be significant, since the Y16A, I20A, L25A and F28A mutants all had decreased bactericidal activity against S. aureus [28]. Surface plasmon resonance (SPR) studies supported the antibacterial assays, identifying F28A-HNP1 as the second lowest self-association on the HNP1 surface of all of the alanine scanning mutants, which corresponded to its second worst antibacterial activity, after W26A. However, in spite of the power of alanine scanning, it can generate misleading results when at least two side chains interact, which can cause non-additivity [38,39]. In order to circumvent this potential shortcoming and shed light upon the cumulative effects of side chain interactions, if any, multiple mutants can be constructed. Crystallography identified two pairs of side chains that would be particularly insightful: Tyr 16 /Phe 28 and Ile 20 /Leu 25 . Tyr 16 and Phe 28 mapped to the dimer interface, making side chain contacts with the opposing monomer. The canonical dimer interface of HNP1 is formed by the antiparallel extensions of the b2 strands stabilized by the reciprocal main-chain hydrogen bonds contributed by Thr 18 and Ile 20 , and the hydrophobic packing of the side chains of Tyr 16 , Tyr 21 , Trp 26 , Phe 28 and the Cys 2 -Cys 30 and Cys 4 -Cys 19 disulfides (Fig. 1A) [27,28]. Further analysis of probable quaternary structures by the Protein Interfaces, Surfaces and Assemblies (PISA) software indicates that the HNP1 dimer forms a tetramer or dimer of dimers in a crystal. The dimer-dimer interface consists of the Ile 20 and Leu 25 side chains and the assembly is maintained exclusively through van der Waals interactions (Fig. 1B).
Here we report the x-ray crystal structures of I20A-HNP1, I20A/I25A-HNP1, Y16A/F28A-HNP1 and Y16A/I20A/L25A/ F28A-HNP1. The functional consequences of each mutation are also investigated using SPR, fluorescence polarization, an enzyme kinetic assay to quantify LF inhibition, and the virtual colony count assay to quantify antibacterial activity.

Synthesis and Folding of Defensins
All five HNP1 mutants (Table 1) were synthesized on an ABI 433A synthesizer using an optimized HBTU activation/DIEA in situ neutralization protocol developed by Kent and coworkers for Boc chemistry solid phase peptide synthesis (SPPS) [40]. All peptides were purified by C 18 reversed phase high performance liquid chromatography (RP-HPLC), and their molecular masses ascertained by electrospray ionization mass spectrometry (ESI-MS) (data not shown). Correct folding of the Y16A/F28A, I20A/ L25A, I20A, and L25A mutants was achieved at 0.25 mg/mL in 25% N,N-dimethylformamide, 2 M urea, 50 mM Tris/HCl, 3 mM reduced and 0.3 mM oxidized glutathione (pH 8.3) overnight at room temperature [41]. Y16A/I20A/L25A/F28A-HNP1, with the side chains of Cys 9 and Cys 29 acetamidomethyl (Acm) protected, was first oxidized at 0.25 mg/mL in 50 mM Tris/HCl buffer (pH 8.3) by stirring in the open air overnight to form two disulfide bridges, Cys 2 -Cys 30 and Cys 4 -Cys 19 . The desired folding intermediate, confirmed by disulfide mapping through digesting with 0.1 mg/mL chymotrypsin for 1.5 hr at room temperature in 50 mM Tris/HCl, 20 mM CaCl 2 , 0.005% Triton X-100 (pH 8.3), was then treated by 0.5 mM iodine for 45 min at 0.5 mg/mL in an acidic solution containing 0.1 M citric acid, 0.2 M HCl and 20% methanol to deprotect Acm and form the third disulfide bond. The reaction was quenched by 0.2 M ascorbic acid, and the fully folded peptide was purified to homogeneity by RP-HPLC and its molecular mass verified by ESI-MS. Defensin stock solutions prepared with water were quantified spectroscopically at 280 nm using molar extinction coefficients calculated according to the algorithm of Pace et al [42].

Crystallization and Data Collection
Lyophilized HNP1 mutant proteins were dissolved in water (20 mg/mL), mixed in a 1:1 ratio (1 ml total) with appropriate precipitant solutions and left to equilibrate in a hanging drop at room temperature. Crystals were flash frozen in liquid nitrogen after briefly soaking in the crystallization condition plus 15-25% glycerol prior to data collection.
For the I20A-HNP1 and I20A/L25A-HNP1 mutants, data were collected using a rotating anode x-ray generator Rigaku-MSC Micromax 7 and a Raxis-4++ image plate detector (at the Xray crystallography core facility, University of Maryland, Baltimore). Diffraction data for HNP1 Y16A/F28A and Y16A/I20A/ L25A/F28A mutants were collected at the Stanford Synchrotron Radiation Light Source (SSRL) BL7-1 beamline on an ADSC Quantum 315 area detector. All data were processed and reduced with HKL2000 [43]. Structures were solved by molecular replacement with Phaser [44] from the CCP4 suite based on the coordinates of the HNP1 monomer (PDB: 3GNY). Refinement was carried out with Refmac [45] and/or Phenix [46] and model building was done with COOT [47]. Data collection and refinement statistics are shown in Table 2. Ramachandran statistics were calculated with Molprobity [48] and illustrations were prepared with Pymol molecular graphics (http://pymol.org) or Molscript [49].

SPR Based wt-HNP1 and gp120 Binding
Experiments were performed at 25uC on a BIAcore T100 System (BIAcore, Inc., Piscataway, NY). The assay running buffer was 10 mM HEPES, 150 mM NaCl, 0.05% surfactant P20, pH 7.4 (63 mM EDTA). 233 response units (RU) of HNP1 or 2770 RU of gp120 were immobilized on CM5 sensor chips using the amine-coupling chemistry. Analytes were introduced into the flow-cells at 30 ml/min in the running buffer. Association and dissociation were assessed for 5 and 10 min, respectively. After each analysis, the sensor chip surfaces were regenerated with 30 mM HCl for HNP1, 10 mM glycine solution (pH 2.0) and 10 mM NaOH for gp120, and equilibrated with the running buffer before the next injection. Binding isotherms were analyzed with BIAevaluation software.

Fluorescence Polarization-based Defensin and N36-gp41 Binding
An N-terminally acetylated N36 (HIV-1 gp160 546-581) peptide derived from the N-terminus of gp41 was synthesized by Boc-chemistry SPPS and purified to homogeneity by preparative C18 RP-HPLC. Succinimidyl ester-activated carboxyfluorescein (FAM-NHS) was covalently conjugated to N-acetylated N36 peptide via its Lys 574 (HIV gp160 numbering) side chain in DMF, and the resultant product N-acetyl-N36-FAM was HPLCpurified and lyophilized. The defensin-N36 binding experiments were performed in 384-well plates on a Tecan Infinite M1000 multimode plate reader. 2-fold serially diluted defensins were prepared in PBS and incubated with 50 nM FAM-labeled N36 in a total volume of 100 ml per well. After a 30 min incubation at room temperature, fluorescence polarization values were measured at exciting and emitted wavelengths of 470 nm and 530 nm, respectively.

LF inhibition kinetics
A 2-fold dilution series of defensin, ranging from 1024 to 1 nM in 20 mM HEPES buffer containing 1 mM CaCl 2 and 0.5% Nonidet P-40 (pH 7.2), was incubated at 37uC for 30 min with 1 mg/mL (,10 nM) of LF. Then, 20 ml of 1 mM LF substrate was added into each well to a final concentration of 100 mM in a total volume of 200 ml. Kinetic measurements of LF enzymatic activity were monitored at 405 nm over 30 min at 37uC on a Tecan Infinite M1000 microplate reader. Data were presented in a plot showing percent inhibition versus defensin concentration, from which IC 50 values (the concentration of defensin that reduced the enzymatic activity of LF by 50%) were derived by a non-linear regression analysis [27,28].

Virtual Colony Count
Antibacterial assays against Escherichia coli strain American Type Culture Collection (ATCC) 25922, Staphylococcus aureus ATCC 29213, and Bacillus cereus ATCC 10876 were conducted using the Virtual Colony Count 96-well kinetic turbidimetric method [22]. Strains were obtained from ATCC (Manassas, VA). A 2-fold dilution series of defensin, ranging from 256 to 1 mg/mL, was incubated with ,5610 5 virtual colony forming units (CFUv)/mL bacteria at 37uC for 2 h in 10 mM sodium phosphate buffer, pH 7.4, 1% tryptic soy broth (TSB), followed by addition of twiceconcentrated Mueller-Hinton broth and kinetic measurements of bacterial growth at 650 nm every 5 minutes over 12 h using a Tecan Infinite M1000 plate reader set to shake 3s orbitally before each read. B. cereus was also assayed using a 2-fold dilution series of defensin ranging from 4-0.016 mg/mL. The 10 mM sodium phosphate incubation buffer included 1% TSB to increase the sensitivity of bacteria to defensin activity [9]. The zero time point was discarded and the optical density at the 5-minute time point was subtracted from subsequent growth kinetic optical density readings, as previously described [23]. Sextuplicate calibration curves were measured at a threshold change in optical density at 650 nm of 0.05. The virtual LD 50 (vLD 50 ), vLD 90 , vLD 99 , and vLD 99.9 were reported as the defensin concentration that resulted in survival rates of 0.5, 0.1, 0.01, and 0.001, respectively. Measurements were done in triplicate on three separate days, except the assay of Y16A/I20A/L25A/F28A-HNP1 against B. cereus from 256 to 1 mg/mL, which was done in triplicate on the same 96-well plate. Statistical p-values were calculated as the output of the paired two-tailed Microsoft Excel T.TEST function.

Total Chemical Synthesis and Oxidative Folding of HNP1 Analogs
Crude HNP1 mutants, after HF cleavage and ether precipitation, gave rise to molecular masses in agreement with the expected values calculated on the basis of the average isotopic compositions of reduced defensins. Using the efficient protocol established for folding of HNP1, the reduced analogs Y16A/F28A-HNP1, I20A/ L25A-HNP1, I20A-HNP1, and L25A-HNP1 productively folded in the presence of 2 M urea and 25% DMF, giving the three disulfide bridges of Cys 2 -Cys 30 , Cys 4 -Cys 19 and Cys 9 -Cys 29 . For Y16A/I20A/L25A/F28A-HNP1, Cys 9 and Cys 29 were strategically selected to be protected by Acm during peptide synthesis. The correctly folded product could then be produced in high yield through a 2-step folding procedure. Two folding species were produced after the first air-oxidation reaction, one of which generated three desired major segments after chymotrypsin digestion: [AC 2 Y][AAC 29

Crystallographic Monomers of HNP1 Analogs are Similar to the Previously Determined Structure of the HNP1 Monomer
The crystal structures of I20A-HNP1, I20A/L25A-HNP1, Y16A/F28A-HNP1 and Y16A/I20A/L25A/F28A-HNP1 were solved by molecular replacement at resolutions between 1.66-2.0 Å . In each case, multiple copies of defensin molecules were present in the asymmetric unit of the crystal ( Table 2). Structural alignment of these crystallographically independent monomers reveals no changes to the core structure, which is the same as in wild type HNP1 (Fig. S1). When superimposed, the root-meansquare (RMS) deviations between 30 equivalent Ca-atoms of mutants and wild-type HNP1 are in the range of 0.51-0.78 Å . There are no changes to the network of three disulfide bridges. Both the disulfide bond connectivity and stereochemistry are identical to wild-type HNP1. In addition to the flexible termini, the major difference is observed for the loop region connecting the second and third b strands. Since the b2-b3 loop of HNP1 is involved in intensive crystal contacts, these differences are attributed to the effect of crystal packing rather than determined by the intrinsic properties of the backbone.

Mutants Exhibit a Variety of Oligomerization Structures
Analysis of HNP1 quaternary structures by PISA indicates that Tyr 16 and Phe 28 are essential for HNP1 dimerization contributing 16 and 11%, respectively, of the solvent-accessible surface buried by each monomer at the dimer interface.
Recently we have shown that single point mutations of Y16A or F28A had no effect on the ability of HNP1 to form stable dimers in crystals [28]. To further test the role of Tyr 16 and Phe 28 in mediating the HNP1 dimer, we prepared the double mutant Y16A/F28A-HNP1. Analysis of intermolecular contacts within the Y16A/F28A-HNP1 crystal indicates that mutant monomers assemble into compact dimers reassembling canonical dimers of wild-type HNP1 (Fig. 2). The replacement of the bulky side chains of Tyr 16 and Phe 28 by linear Ala side chains allows the dimer to pack more tightly and form an additional set of Hbonds on the opposite face of the dimer. The new set of four Hbonds pins the short b1 strand of one monomer to the equivalent strand of the other monomer to form a two-stranded antiparallel b-sheet ( Fig. 2A). These include the reciprocal mainchain H-bonds contributed by the carbonyl oxygen of Cys 2 and the nitrogen of Cys 4 and main chain nitrogen of Ala 1 of monomer a and the hydroxyl of Tyr 3 of monomer b. Formation of this dimer results in the average burial of more than 495 Å 2 of molecular surface per each monomer, which compares to an average of 370 Å 2 for the wild-type HNP1 dimer (PDB:3GNY, [27]). The significant increase in the average value of the molecular surface buried within inter-monomer interactions of Y16A/F28A-HNP1 as compared to the wild-type dimer is attributed to the presence of additional interactions mediated by the b1 strands. Similarly, a more compact dimer was formed by the replacement of the single Phe 28 at the dimer interface [28]. Superimposition analysis (Fig. S2) of the Y16A/F28A-HNP1 and F28A-HNP1 dimers reveals very close similarity as shown by an average RMDS value of 0.6 Å for 60 aligned Ca atoms and an almost identical molecular surface buried within the dimer (495 and 500 Å 2 per monomer for the Y16A/F28A-HNP1 and F28A-HNP1 dimer, respectively). Introduction of mutations at the dimer interface have no effect on the ability of Y16A/F28A-HNP1 to form a stable tetramer, which closely resembles the wild type HNP1 tetramer (Fig. S2).
Surprisingly, although the asymmetric unit of the I20A-HNP1 crystal contained two monomers they were not arranged into a dimer. Analysis of intermolecular contacts within the crystal unambiguously rules out the formation of any quaternary structure for the I20A-HNP1 mutant.
We failed to obtain diffracting crystals of the L25A-HNP1 mutant. I20A/L25A-HNP1 crystallizes as a tetramer with two independent tetramers in the asymmetric unit of the crystal (Fig. 3 and Fig. 4). The RMSD between tetramers is 0.744 Å . Unlike wild type HNP1 and most other HNP1 mutants, the tetramer I20A/ L25A-HNP1 is not formed as a dimer of dimers. Each monomer contributes equally to tetramer stability with a combined buried surface area of approximately 2760 Å 2 (as compared to 2280 Å 2 of the wild type HNP1 tetramer). Tyr 16 and Phe 28 make up the bulk of the hydrophobic core with contributions from Tyr 21 and the adjacent disulfide bonds of Cys 4 -Cys 19 and Cys 2 -Cys 30 . In addition, there are two core sets of hydrogen bonds in the tetramer center that help in its stability, utilizing main chain atoms from Ala 1 , Cys 2 , and Cys 30 and the hydroxyl of Tyr 16 (Fig. 3). The tetramer is further stabilized by a hydrogen bond between Arg 14 and Tyr 21 of one pair of subunits.
Interactions of HNP1 Analogs with N36 Peptide HNP1 inhibits HIV-1 infectivity via multiple mechanisms [50,51]. Recently, HNP1 has been shown to be able to interact with the heptad repeat domains of of HIV-1 gp41 comprising both N36 and C34 peptides, and contribute to the inhibition of fusion by interfering with the formation of the fusogenic gp41 structure [52]. An N-acetyl-N36 peptide labeled with 6-carboxyfluorescein (FAM) at Lys 574 (HIV gp160 numbering) at 50 nM in PBS was incubated at room temperature for 30 min with a 2-fold dilution series of HNP1 mutants (0.024-50 mM) before fluorescence polarization was measured in a triplicate assay on a Tecan Infinite M1000 multimode plate reader. All HNP1 analogs except Y16A/ I20A/L25A/F28A-HNP1 bound N36 in a dose-dependent fashion, reflected by increasing fluorescence polarization values with the increase in defensin concentration. The binding of HNP1 analogs to N36 was also structure-dependent, since the polarization values were in the order HNP1. L25A.I20A.I20A/ L25A.Y16A/F28A.Y16A/I20A/L25A/F28A at almost all of the concentration series. (Fig. 8).

Inhibition of Lethal Factor by HNP1 Analogs
The inhibition of lethal factor by various defensins was quantified at 37uC using an enzyme kinetic assay [27]. As shown as inhibition curves and IC 50 values in Fig. 9, replacement of Ile 20 and/or Leu 25 with Ala slightly weakened inhibition of LF relative to wild type HNP1. However, the inhibitory activity of Y16A/ F28A and Y16A/I20A/L25A/F28A in particular was greatly  suppressed with the latter yielding an IC 50 95-fold higher than that of wild type HNP1.

Antibacterial Activity of HNP1 Analogs
Shown in Fig. 10 are virtual survival curves of E. coli, S. aureus, and B. cereus exposed to each defensin, and in Table 3    remained less than vLD 99.9 against E. coli and its vLD 99.9 against B. cereus was 11-fold higher than HNP1. S. aureus growth curves were not necessarily parallel to non-defensin controls, suggesting that defensins can affect growth kinetics and delay threshold times after the addition of twice-concentrated Mueller-Hinton broth during the 12 h outgrowth phase of the virtual colony count assay.

Discussion
Defensins form a wide array of oligomeric structures. There is considerable structural diversity among the defensins, and the adefensins are not the only group that has been shown to dimerize and multimerize. The authors of a study of the plant defensin NaD1 mutated a lysine to alanine and observed reduced activity, which they attributed to impaired dimer formation [53]. The activity could also have decreased due to a decrease in cationicity. This type of ambiguity does not apply to our present study, since we only mutated hydrophobic residues to alanine. The crystal structure of human b-defensin 2 depicted dimers and octamers, although the authors commented that it is unclear whether the octamer is physiologically relevant given the lack of conservation of the residues at the various subunit interfaces [54]. By contrast, a study of human b-defensins (HBD) 1-3 depicted HBD1 and 2 as monomers, while HBD3 crystallized as a dimer [55]. The enhanced activity of HBD3 compared to HBD1 and HBD2 may be related to its ability to dimerize. The h-defensin retrocyclin-2 trimerized in solution, according to analytical ultracentrifugation and nuclear magnetic resonance results [56].
Starting with the crystal structure of HNP-3 [29], and the nuclear magnetic resonance (NMR) structures of HNP1 [57,58], wild-type human a-defensin structures have typically depicted dimers. However, monomeric mutant forms of HNP1 and HD5 have been generated by N-methylating an amide bond at the dimer interface, thus breaking a hydrogen bond and introducing steric hindrance to the canonical dimerization scheme [26,59]. In both cases, the monomeric mutant forms were as active against E. coli as the respective wild-type defensins. Other mutant forms have emphasized the versatility of the defensin framework. The HD5 canonical dimer with antiparallel b-strands existed in the crystal structure of the L29Abu-HD5 mutant, while the L29Nle-HD5  mutant exhibited a novel dimeric form with b-strands stacked in parallel [26]. Interestingly, the non-canonical L29Nle-HD5 molecule was functional, exhibiting less than wild-type activity against S. aureus but greater than wild-type activity against E. coli [26]. In the present study, we also discovered novel quaternary structures in the crystal, including a form of I20A-HNP1 that appears to be unable to form dimers or tetramers and a novel tetrameric form of I20A/L25A-HNP1. In addition, the canonical dimer formed by the quadruple mutant displayed a larger interface than wild-type. The HNP1 and HD5 results indicate that the presence of a canonical form of dimerization in the crystal does not correlate with a-defensin activity. Structural diversity may be a consequence of the small size of the peptides, leading to agile dynamic movements, combined with the rigidity of the tridisulfide structure that prevents major conformational transitions within the monomer. Defensin structural dynamics has been studied in solution in the context of binding to membranes or micelles [56,60,61]. The structure of HNP1 investigated using solid-state NMR suggested that the loop region between the first and second beta-strands is flexible and may change conformation when exposed to membranes [58]. However, the effect of membranes on HNP1 quaternary structure is poorly understood.
The stability of the dimer apparently varies from one a-defensin to the next. The HNP4 dimer is weakened by smaller hydrophobic side chains than HNP1 [62]. The HD5 dimer is stabilized by an additional short two-stranded beta sheet and hydrogen bond at the interface, yet dimerization is less important for HD5 than for HNP1 as indicated by the results of studies of forced monomeric defensins [26,59]. However, a study of HD5 activity against nonenveloped viruses demonstrated that dimerization and multimerization are important for antiviral activity [63]. There is also diversity in modes of tetramerization among the a-defensins. The Ile 20 /Leu 25 -mediated HNP1 tetramer does not have an equivalent form in the cases of HNP4 or HD5, although there is evidence that HD5 forms tetramers [64]. Even though it has almost no direct antibacterial activity [22], HD6 exhibits the ability to form long chains of monomers [6,62] and nanonets that trap bacteria and protect transgenic mice expressing HD6 from Salmonella enterica serovar Typhimurium infection [65].
Hydrophobicity has repercussions beyond oligomerization. In the case of HNP1, the finding that mutating four hydrophobic residues depicted in the crystal as contacting one another to alanine interferes with HNP1 function indicates that hydrophobicity at these positions is potentially evolutionarily significant. However, these results with HNP1 do not apply to a-defensins unviersally. Hydrophobicity is not conserved at any of the four positions (16, 20, 25, or 28) in an alignment of mammalian adefensins. Radically different examples exist, such as the more basic rabbit neutrophil peptide NP-3a (Swiss-Prot identifier DEF1_RABIT), which has Ser at position 16 and Arg at positions 20, 25 and 28 (HNP1 numbering) [66]. Cationicity could be more significant than hydrophobicity in these cases, whereas for HNP1 and HD5 alanine scanning mutagenesis showed that hydrophobicity is clearly more important than cationicity [26,28]. For HD5, positions 16, 20, 25, and 28 (HNP1 numbering; to translate into HD5 numbering add 1) are occupied by Ser, Glu, Leu and Leu, respectively, and the latter Leu, Leu 29 , has been shown to be the most important residue for antibacterial and LF-inhibiting activity [26]. A defensin like NP-3a would presumably form much different quaternary structures than HNP1, if any, because the electrostatic repulsion of the arginine side chains would make dimer and tetramer structures such as those described in the crystallographic section of this work unlikely.
The cumulative functional effect of the quadruple alanine mutations was similar to the effect of the single mutation, W26A [28]. This result sheds new light upon the role of Trp 26 , since interfering with this residue forming the hydrophobic core probably changes the conformation of several different side chains, including Phe 28 . The same study also showed that F28A had the second least SPR self-association of all HNP1 alanine scanning mutants, and the second lowest activity against S. aureus [28]. Interestingly, the F28A mutation did not prevent HNP1 from inhibiting LF and binding gp120, indicating that self-association and S. aureus activity are more sensitive to perturbations of HNP1 structure at this position. Here, our studies using the double and quadruple mutants that include the F28A mutation extend the previous findings. Unlike the F28A mutation alone, the Y16A/ F28A double mutant and the quadruple mutant LF inhibition IC 50 values increased 6-fold and 95-fold, respectively. In the virtual colony count assay, the double mutant Y16A/F28A-HNP1 had a vLD 99.9 value against S. aureus 7-fold higher than HNP1, and Y16A/I20A/L25A/F28A-HNP1 failed to exhibit activity equivalent to the vLD 90 level against S. aureus. The Y16A/F28A and . Strains were exposed a twofold dilution series of defensins at concentrations varying from 1 to 256 mg/mL, except B. cereus was exposed to all defensins at 0.016 to 4 mg/mL and only Y16A/I20A/L25A/F28A-HNP1 at 1 to 256 mg/mL. Points equivalent to zero survival cannot be plotted on a logarithmic scale, such that against S. aureus virtual survival was zero above 32 mg/mL for HNP1, above 64 mg/mL for L25A-HNP1, and at 256 mg/mL for I20A-HNP1; against B. cereus virtual survival was zero at 4 mg/mL for HNP1. Each point is the mean of triplicate measurements, except HNP1 at 32 mg/mL against S. aureus (mean of two measurements; the third replicate gave a virtual survival value of zero) and L25A-HNP1 at 4 mg/mL against B. cereus (single measurement; the other two replicates gave virtual survival values of zero). Although error cannot easily be read directly from this plot, error values are quantified in Table 3. doi:10.1371/journal.pone.0078937.g010 Y16A/I20A/L25A/F28A mutants were also significantly less potent than HNP1 when assayed against both E. coli and B. cereus.
The functional consequences of the I20A and L25A mutations, alone or in concert, were milder than the Y16A/F28A double mutant or the quadruple mutant. LF inhibition was less than 2fold different than wild-type for the single mutants and about 2fold different for the I20A/L25A double mutant, as measured by IC 50 values. In the virtual colony count assay, whereas the I20A and L25A mutations each made a significant difference against S. aureus, they did not make a significant difference against E. coli or B. cereus even in the I20A/I25A double mutant. Since Ile 20 and Leu 25 function as an isologous binding site for the formation of a dimer of dimers, their effect or lack thereof on E. coli is not surprising, given a previous study of N-methylated Ile 20 that showed that dimerization is dispensible for E. coli activity [59]. The fact that the I20A/L25A double mutations were irrelevant, but the Y16A/F28A double mutations were quite relevant, for E. coli activity is insightful, because it demonstrates that the hydrophobicity of Tyr 16 and Phe 28 has repercussions beyond dimerization. If providing hydrophobic contacts for dimerization were the sole function of Tyr 16 and Phe 28 , we would expect the activity of Y16A/F28A-HNP1 to be the same as wild-type HNP1, just as there was no difference between HNP1 or MeIle 20 -HNP1 when assayed against E. coli [59].
Functional assays correlated with each other and with hydrophobicity. Here we observed the binding affinity of HNP1 and mutants to HNP1 and gp120, as measured by SPR, and the N-terminus of HIV-1 gp41, as measured by fluorescence polarization. In all cases, binding was in agreement with LF inhibition and activity against S. aureus, in the order HNP1. L25A-HNP1. I20A-HNP1. I20A/L25A-HNP1. Y16A/F28A-HNP1. Y16A/I20A/L25A/F28A-HNP1. Correlation between these quantities was expected, given the previously observed correlation for mutants of Trp 26 [28]. Self-association on the HNP1 surface was lessened by the mutations, reinforcing that these four positions are important for dimerization and multimerization. Diminished association with the two HIV-1 proteins also suggests that dimerization and multimerization are important for antiviral function. However, the diversity of monomeric, dimeric and tetrameric crystallographic forms suggest that propensity to form certain quaternary structures in the crystal does not always inform or predict defensin function.
Although these crystallographic results do not support the premise that crystal structures are necessarily uniformly functionally insightful, the convergence of evidence from the wide variety of other methods employed has yielded the ability to rank the contribution of each of the side chains mutated to alanine. According to the Wimley-White experimentally determined hydrophobicity scales for proteins at membrane interfaces [67], the five most hydrophobic residues are Trp.Phe.Tyr.Leu.Ile. Our studies of HNP1 have determined the relative importance of residues to be the same, save for the transposition of Ile and Leu: Trp 26 .Phe 28 . Tyr 16 . Ile 20 . Leu 25 . The Wimley-White scale was determined using large unilamellar vesicles comprised of the zwitterionic phospholipid palmitoyloleoylphosphatidylcholine; the scale might have differed had they used more anionic membranes typical of bacteria. Nevertheless, the results of our alanine scanning mutagenesis studies of HNP1 with single, double, and quadruple alanine mutants emphasize the importance of hydrophobicity as the primary factor that determines activity, and indicate that the importance of a residue is roughly proportional to its hydrophobicity regardless of structural location.
Although these residues are shown to be at oligomeric interfaces in most crystal structures, the consequences of decreasing hydrophobicity at any of these positions on crystallographic results indicates that the canonical dimer is delicate. In the dynamic environment of the liquid phase as opposed to the solid crystal, a loose association between monomers may allow hydrophobic residues to exchange their affinities for one another for interactions with carbon atoms in diverse targets such as bacterial membranes, enzymes and proteins. Weak van der Waals interactions can be interchangeable and do not necessarily result in specificity. Therefore, while the significance of these four residues for function has been clearly established, the degree to which oligomers such as those implied by the canonical dimer contribute to defensin activity remains an open question worthy of further study.  [28]). Both dimers are stabilized by the same network of main chain H-bonds (shown as blue dashes). The H-bond formed between hydroxyl groups of Tyr 3 in the F28A-HNP1 dimer (shown as magenta dashes) is replaced by the H-bond formed between the main chain nitrogen of Ala 1 and the hydroxyl of Tyr 3 in Y16A/F28A-HNP1 dimer (shown as light blue dashes). The molecular surface buried within the Y16A/F28A-HNP1 dimer is 494 Å 2 per monomer, which compares to 500 Å 2 for the F28A-HNP1 dimer. (B) Structural alignment of Y16A/F28A-HNP1 and wild type HNP1 tetramers. Tetramers were aligned based on the AB dimer and residues involved in tetramer formation are shown as balls and sticks. The Y16A/F28A-HNP1 dimers are colored as in (A) and the HNP1 dimers (PDB:3GNY, [27]) are colored red and blue. The molecular surface buried within the Y16A/F28A-HNP1 tetramer is 1420 Å 2 per dimer, which compares to 1140 Å 2 for the HNP1 tetramer. (TIF)